Ecología

 

Factors affecting woody plant species diversity of fragmented seasonally dry oak forests in the Mixteca Alta, Oaxaca, Mexico

 

Factores que afectan la diversidad de especies leñosas en fragmentos de bosque de encino estacionalmente seco en la Mixteca alta oaxaqueña, México

 

Remedios Aguilar-Santelises* and Rafael F. del Castillo

 

Centro Interdisciplinario de Investigación para el Desarrollo Integral Regional, Unidad Oaxaca (CIIDIR-Oaxaca). Instituto Politécnico Nacional. Calle Hornos 1003, Col. Nochebuena, 71230 Santa Cruz Xoxocotlán, Oaxaca, México. *mraguilars@ipn.mx

 

Recibido: 28 febrero 2012
Aceptado: 03 diciembre 2012

 

Abstract

We explored the relationship between fragment area, topographic heterogeneity, and disturbance intensity with tree and shrub species diversity in seasonally dry oak forest remnants in the Mixteca Alta, Oaxaca, Mexico. The fragments are distributed in a matrix of eroded lands and crop fields, have a complex topography, and are disturbed by plant extraction and trail opening. Sampling was conducted in 12 fragments from 12-3 211 ha. Topographic heterogeneity was estimated by the fragment's standard deviation in slope-aspect, slope, and altitude. The density of stumps and roads were used as estimators of disturbance intensity. Fisher's αdiversity ranked from 0.95 to 4.55 for the tree layer; and 2.99 to 8.51, for the shrub layer. A structural equation model showed that the diversity of woody plants increases with topographic heterogeneity and disturbance in the remnants. When these 2 variables were considered, diversity tended to decrease with fragment size probably because smaller fragments have a greater perimeter-to-area ratio and therefore proportionally offer more opportunities for pioneer species colonization. Indeed, the tree-to shrub-layer diversity ratio increased with fragment size. Conservation strategies in fragmented forests must consider the fragment's environmental heterogeneity, the disturbance type and intensity, and the species to be preserved.

Key words: fragmentation, seasonally dry oak forest, human disturbance, species-area relationship, topographic heterogeneity, structural equation modeling.

 

Resumen

Exploramos la relación entre el área, la heterogeneidad topográfica y el disturbio en remanentes de bosque de encino estacionales en la Mixteca Alta, Oaxaca, México. Una matriz de suelo erosionado y cultivos rodea los fragmentos, que están afectados por extracción vegetal y caminos y presentan topografía compleja. Muestreamos la vegetación en doce fragmentos de 12 a 3 211 ha. Estimamos la heterogeneidad ambiental con las desviaciones estándar en pendiente, orientación y altitud del fragmento, y la intensidad de disturbio, por la densidad de tocones y el área afectada por caminos. La diversidad αde Fisher varió entre 0.95 y 4.55 para el estrato arbóreo y 2.99 y 8.51 para el arbustivo. Un modelo de ecuaciones estructurales lineales mostró que la diversidad aumenta con la heterogeneidad topográfica y la perturbación humana. Al considerar estas dos variables, la diversidad disminuyó con el tamaño del fragmento probablemente porque proporcionalmente los fragmentos pequeños tienen mayor perímetro que los grandes y favorecen a las pioneras. La razón entre la diversidad del estrato arbóreo y el arbustivo aumentó con el tamaño del fragmento. Las estrategias de conservación en bosques fragmentados deben considerar la heterogeneidad ambiental, el disturbio y las especies que deben ser conservadas.

Palabras clave: fragmentación, bosques de encino estacionalmente secos, disturbio, relación especies-área, heterogeneidad ambiental, modelaje de ecuaciones estructurales.

 

Introduction

Fragmentation processes involve habitat losses and the splitting of the remaining habitat into pieces of various sizes and degrees of isolation (Laurance, 2008). Currently, a large part of the land surface is being affected by human activities, causing ecosystem fragmentation and jeopardizing biodiversity through habitat reduction, increased isolation, and alterations in biotic and abiotic factors in the remaining fragments (Saunders et al., 1991; Fahrig, 2003; Wade et al., 2003; Otálora et al., 2011). Several factors have been associated with biodiversity in fragmented landscapes. These include fragment size per se, based on the species-area theory (Arrhenius,1921; Preston, 1962; MacArthur and Harris, 1984; Tjorve, 2003), environmental heterogeneity (Gaston, 2000; Tews et al., 2004; Clarke and Gaston, 2006), and disturbances, both of natural and anthropogenic origin (Bustamante and Grez, 1995; Williams-Linera et al., 2002; Wade et al., 2003; Davis, 2004).

A positive relationship between biodiversity and fragment or habitat area has been identified for nearly a century by the widely-known species-area relationship (Arrhenius, 1921; Preston, 1962; Bustamante and Grez, 1995). This relationship can be described in a probabilistic model following a geometric (Arrhenius, 1921) or logarithmic series (Preston, 1962), enabling the researcher to estimate the biodiversity of an ecosystem from a known area. Several theoretical and empirical studies in fragmented landscapes have found a close relationship between patch biodiversity and patch size (Hill and Curran, 2001, 2003; Echeverría et al., 2007; Pincheira-Ulbrich et al., 2008).

Although the species-area relationship is one of the main subjects in biogeography (Hill and Curran, 2001, 2003; Echeverría et al., 2007; Pincheira-Ulbricht et al., 2008; Blakely and Didham, 2010), it is not clear which mechanisms are at work. Numerous studies have suggested that environmental heterogeneity, which is usually positively correlated with area, is the main factor explaining biodiversity (Boecklen, 1986; Freemark and Merriam, 1986; Baz and García-Boyero, 1995; Brose 2001; Aström et al., 2007; Blakely and Didham, 2010). Indeed, heterogeneous environments offer greater diversity of niches for the establishment of different species (Boecklen, 1986; Baz and García-Boyero, 1995; Peterson et al., 1997; Tews et al., 2004; Hannus and von Numers, 2008). Complex topography is one of the most distinctive features of environmental heterogeneity in mountain ecosystems by altering soil depth, moisture content, stoniness, compaction, and permeability, among other environmental properties, thereby creating more niches per area than those occurring on a flat surface (Bunting, 1964; Balvanera and Aguirre, 2006; Aström et al., 2007).

Disturbance is another important factor affecting diversity. Disturbance has been defined as a more or less discrete event in time and space, altering the structure of populations, communities or ecosystems, causing drastic changes in resource availability or in the physical environment (Saunders et al., 1991; Bustamante and Grez, 1995; Laurance, 2004; di Bella et al., 2008), facilitating the spreading of short-lived early successional species (Saunders et al., 1991), and the invasion of exotic species that compete with native species for resources (Santos and Tellería, 2006; Stevenson and Rodriguez, 2008). Although disturbance is an important component of many ecosystems, there is no consensus on how it impacts biodiversity (Miller et al., 2011). Disturbance can be of natural origin, such as storms, telluric events, and tree falls, or anthropogenic, as is the case of human settlements, roads, deforestation, and fire. In this study we will focus on human disturbances, which often are difficult to measure directly, but can be estimated by their effects on the fragments of natural ecosystems. The occurrence of roads and stumps are signals of human disturbance (López, 2001; Williams-Linera et al., 2002; Herrera et al., 2004; Rudas et al., 2007), and they may modify the ecosystem structure and composition by: a) affecting microclimatic conditions (Gucinski et al., 2001); b) promoting the invasion of exotic species (Brown et al. 2004, 2006); c) allowing the uncontrolled extraction of natural products (Young, 1994; Verburg et al., 2004); d) setting up barriers between populations that may decrease gene flow and dispersal (Forman and Alexander, 1998), and e) reducing seed production (SEMARNAT, 2005; Alelign et al., 2007), all of which may jeopardize species persistence.

Disturbance is also related to fragment area and habitat heterogeneity. For instance, native species richness per unit area may decrease significantly in small-sized and highly disturbed fragments (Ross et al., 2002; Echeverría et al., 2007). Furthermore, road construction and deforestation induce habitat fragmentation, promote changes in the physical environment, and alter the biota balance (Saunders et al., 1991; Fahrig, 2003; Wade et al., 2003; Otálora et al., 2011). Environmental changes in fragmented communities are more dramatic at the edges than at the center of the fragments (Murcia, 1995; Laurance et al., 2000; Forero-Molina and Finegan, 2004). Some studies have shown that basal area significantly declines with decreasing patch size (Lezcano et al., 2004; Echeverria et al., 2007). The fragment species composition is also affected. The shrub layer diversity to tree layer diversity ratio could be an indicator of disturbance since these life forms usually have different environmental requirements. Shrubs, for instance, tend to have greater survivorship and biomass in open microsites (Asbjornsen et al., 2004b).

In summary, biodiversity may be directly or indirectly associated with fragment area, topographic heterogeneity, and disturbance intensity. For both conservation purposes and ecological studies, it is important to identify the major factors that influence biodiversity in remnant fragments. Few studies have explored simultaneously the role that each of the above factors plays on fragment diversity. Most of these studies have explored only fragment size and environmental heterogeneity (Freemark and Merriam, 1986; Baz and García-Boyero, 1995; Boecklen, 1986; Brose, 2001; Graham and Blake, 2001), but very few have included the effects of area and disturbance on biodiversity (Aström et al., 2007). We are only aware of one exploratory study that considers simultaneously the relationship between habitat size, environmental heterogeneity, and disturbance with diversity (Blakely and Didham, 2010). This study, carried out with insects, surprisingly found a negative relationship between biodiversity and habitat size under experimentally controlled conditions, due explicitly to experimental manipulation in which smaller habitats were modified to be more heterogeneous than larger habitats. This was accomplished by experimentally reversing resource concentration and enhancing drought disturbance, while holding constant colonization-extinction dynamics and habitat heterogeneity. The Mixteca Alta in southern Mexico provides a suitable landscape for exploring these relationships but in a natural habitat and with woody plants. This region is highly fragmented (Asbjornsen et al., 2004b; Martínez and Noriega, 2006), has a complex topography (González-Leyva, 2007), and has been affected by disturbances associated with the presence of nearby human settlements (Asbjornsen et al., 2004a). In this region, we aimed to explore the possible relationships between fragment area, topographic heterogeneity, and intensity of anthropogenic disturbance on fragment biodiversity of trees and shrubs. Based on the empirical evidence and theoretical studies described above, we expected (1) a positive relationship between fragment biodiversity and both area and topographic heterogeneity; and (2), a negative relationship between fragment biodiversity and disturbance intensity.

 

Materials and methods

Study site and sampling design. The study site is located in the Nochixtlán District, Oaxaca, Mexico, at 17°0'-17°50' N, 97°0'97°25' W, between 1 800 and 2 800 m. The study area is mountainous, with a complex geology and topography. The climate is temperate and semi-humid. Annual rainfall varies between 500 and 800 mm. A seasonally dry oak forest comprises most of the vegetation above 1 500 m (Asbjornsen et al., 2004a). The main species are: Quercus liebmannii, Q. acutifolia, and Q. laurina. Endemic species are also relatively abundant (García-Mendoza et al., 1994). Paleontological evidence shows that the Mixteca Alta has been populated since the late Holocene by people who based their use of resources on a wise water management (Guerrero-Arenas et al., 2010). The Spanish conquest was accompanied by the introduction of sheep, goats, and diverse crops, causing an intense process of deforestation. After 500 years, deforestation has resulted in a highly fragmented landscape: 80% of its soils are affected by water erosion (González-Leyva, 2007; Guerrero-Arenas et al., 2010).

Data collection. We selected and characterized our study fragments using Landsat ETM+ satellite imagery (2005, path 24, Row 48, Band 4/7/1, pixel= 30 m). The geographic projection was UTM WGS84 14n zone. Geographic corrections were conducted with control points from digitalized 1: 50 000 road maps, and a second degree polynomial model (Cayuela et al., 2006b). Accuracy ranged from 0.25-0.45 pixels, corresponding to 7.5-13.5 m. Atmospheric corrections were performed using the Chávez reflectivity model (Chuvieco, 2006; Cayuela et al., 2006b), which transforms the original digital numbers into reflectivity values in the corrected images. Elevation digital models were generated using ENVI 4.3 software. The topographic correction was performed using the Teillet et al. (1982) and Riaño et al. (2000) semi-empirical method and the PCI Geomatics software version 7.0. Classification was supervised with PCI Geomatics software version 7.0 using the maximum likelihood criterion. Six class signatures were obtained: 1) bare land, 2) water body, 3) grasslands-shrublands, 4) croplands, 5) urban areas, and 6) native forest. The obtained classification was checked with 300 independent control points located in the field with Google Earth software (image dates from 2004 to 2007). We haphazardly selected 12 seasonally dry oak forest fragments with contrasting areas ranging from 12 to 3 211 ha (Table 1). The area of the selected fragments was estimated using Fragstats (McGarigal and Marks, 1995). Sampling plots were randomly selected within each fragment using the extension Random Sites (Arc View 3.X, public domain), with the restriction that each sampling point should be located at least 70 m from other sampling points or from the fragment edge to avoid overlapping and to decrease edge effects and the probability of autocorrelation between nearby sampling points (Fig. 1).

Vegetation sampling was conducted on 216 plots distributed among the 12 selected fragments, using 4 to 50 sampling points per fragment, depending on fragment size. All plots were geo-referenced, using a GPS (GARMIN 60csx) with a 5 m resolution. In order to analyze the diversity and structure of the vegetation, we sampled the individuals of both the tree and the shrub layer. All woody plants ≥ 2.5 cm DBH and ≥ 2.5 m height found in the sampling plots were included in the tree layer; whereas all woody plants < 2.5 cm DBH or < 2.5 m height were included in the shrub layer. Sampling plots of 102.06 m² and 12.56 m² were used for the tree and the shrub layer, respectively. Specimens of all species found in the sampling plots were deposited at the Herbarium OAX. The tree and shrub diversity per fragment was assessed using Fisher's αbecause it is relatively insensitive to sample size (Fisher et al., 1943; Magurran, 2004).

Topographic heterogeneity was evaluated by the standard deviation (sd) of altitude (m), slope (%), and slope-aspect (-cos ϕ, where ϕis the slope angle in radians) of each fragment. Slope-aspect is defined as main compass direction that a slope faces (Physical Geography Dictionary, 2012). We used the Hawthtools extension of ArcGis to select randomly 50 points (pixel= 30 m), at least 70 m apart. The selected points were overlapped on the slope, altitude, and slope-aspect layers from the digital elevation model to obtain the respective values for each point. Anthropogenic disturbance was estimated by assessing the intensity of logging and the proportion of the fragment area expected to be affected by roads or zone of influence of roads within each fragment. Logging intensity was evaluated by the density of stumps (ha^-1) in each plot. Previous works have found that road effects on biodiversity depend on species, topography, and road type, but usually range between 100 and 200 m on each side of the road (Forman et al., 1997). Based on these studies, we defined a buffer area of 150 m width on both sides of the roads in the study area, to estimate the proportion of the fragment expected to be affected by roads. For this purpose, we used EPS data from Inegi (2011, scale 1:50 000). Road-effect zone was estimated as the ratio of the road buffer area to the total area of the fragment. In our study sites, all roads were of similar width, suggesting the same intensity of use.

Data analysis. In order to disentangle the relationships between fragment area, topographic heterogeneity, anthropogenic disturbance, and woody plant species diversity, we developed a model based on structural equation modeling, using the CALIS procedure from the SAS 9.1 software package (SAS Institute, 1989). Since large fragments are probably both more heterogeneous and less disturbed than small fragments, it is important to explore to which extent biodiversity is directly affected by fragment area, topographic heterogeneity, and disturbance, or indirectly through the associations among these explanatory variables. Structural equation modeling (SEM) allows the testing of complex relationships among variables, partitioning direct and indirect effects, and making quantitative predictions about the relative contribution of each variable in the model (Grace and Pugesek, 1997). This method, based on covariance analysis, can be used to model multivariate relations and to test multivariate hypotheses (Bollen, 1989). An important attribute of structural equation modeling is that it allows the estimation of conceptual unmeasured (latent) variables based on a set of measurable (manifest) variables (Grace and Pugesek, 1997). The accepted models obtained from such analysis can indicate the role of different factors in a system and the strength of their relationships (Spitale et al., 2009). Further details of the SEM methodology can be found in Grace and Pugesek (1997), Spitale et al. (2009), Hayduk (1987) and Reed et al. (2009). In our model, biodiversity is the endogenous (response) latent variable, estimated by the tree and shrub diversity (manifest variables). Topographic heterogeneity, anthropogenic disturbance, and fragment area are our exogenous (explanatory) variables. Topographic heterogeneity is a latent variable estimated by the standard deviation of altitude, slope, and slope-aspect. Anthropogenic disturbance is a latent variable estimated by the exogenous manifest variables: road-effect zone, as defined above, and stump density (Fig. 2).

All the response and predictor variables were standardized (mean= 0, sd= 1) and did not show evidence of deviations from normality. Direct relations among variables (single-headed arrows in figure 2 and figure 4) were estimated as standardized coefficients from the covariance matrix. Non-directional standardized correlation coefficients were also calculated among explanatory manifest variables (double-headed arrows in figure 2). The initial structural model was reduced to the most parsimonious model by means of a stepwise specification search, eliminating in each step the path with the lowest coefficient (in absolute value) until all the remaining coefficient paths were significant (Hayduk, 1987; Grace and Pugesek, 1997; Reed et al., 2009; Blakely and Didham, 2010). The resulting model in each reduction step was checked by its goodness-of-fit index (GFI), its chi-square probability value (p), and its Akaike's information criterion (AIC) value (Hayduk, 1987). In each step, the fitted indices were compared against the previous model. The best model was the one with the GFI nearest to 0.9; the greatest p value, which should be > 0.1, and the lowest AIC value (Mulaik et al., 1989; Stoelting, 2002).

 

Results

The tree canopy sample was composed of 3 301 specimens from 46 species (Appendix 1), and the shrub layer was composed of 7 453 specimens from 116 species (Appendix 2). Fisher's αdiversity ranged from 0.95 to 4.55 for the tree layer, and 2.99 to 8.51 for the shrub layer (Table 1). The shrub layer-to-tree layer diversity ratio decreased significantly with the size of the fragment (r² = 0.315, p< 0.01, Fig. 3). The final and most parsimonious model for species diversity in the seasonally dry oak forest remnants of the Mixteca Alta region had acceptable goodness-of-fit indices (GFI= 0.899, p= 0.171, AIC= -1.595), following Hayduck (1987 [Fig. 4]).

Our SEM analysis revealed that disturbance, habitat heterogeneity, and fragment area have significant relationships with species diversity. Of these drivers of species diversity, habitat heterogeneity, here estimated in terms of topographic variables, is the most important and has a positive effect on species diversity. When habitat heterogeneity is considered, the size of the fragment has a significant but negative effect on diversity. Two out of the 3 explanatory variables used to estimate topographic heterogeneity were significant: the slope-aspect heterogeneity and the altitude heterogeneity. According to our SEM analysis, the Fisher's αdiversity in the fragment tends to increase in fragments with high variation in elevation and low variation in slope-aspect. Human disturbance also affected significantly and positively the diversity of the fragments. Of the manifest variables used to estimate disturbance, only stump density, an indicator of the intensity of plant extraction, was significant.

 

Discussion

Structural equation modeling revealed the effect of anthropogenic disturbance, fragment area, and topographic heterogeneity in woody plant species diversity in remnants of seasonally dry oak forests in the Mixteca Alta, southern Mexico. In accordance with our hypotheses, woody plant species diversity per fragment can be explained directly by the topographic heterogeneity, the intensity of anthropogenic disturbance, and the fragment area, as has been shown in other studies (Ross et al., 2002; Cayuela et al., 2006a; Echeverría et al., 2007). Our results indicate that, when topographic heterogeneity and human disturbance are taken into account, the effect of fragment size on diversity is negative in the seasonally dry oak forest of the Mixteca region. More specifically, smaller fragments with similar altitude and slope-aspect and human disturbance level tend to be more diverse than large fragments.

Edge effects provide a possible explanation for this result since the perimeter to area ratio is greater in small fragments. Small fragments provide more opportunities for light tolerant species to become established, favoring a greater diversity of species. Some species may find the habitat of the edges of the fragment more suitable for survival and reproduction than the center of the fragment (e.g., Bernabe et al., 1999; Asbjornsen et al., 2004b). Indeed, in our study shrubs were proportionally more diverse than trees in smaller than in larger fragments. Shrubs are usually more light-demanding than trees, and their abundance has been found to decrease significantly with the abundance of trees in oak forest in the adjacent Sierra Norte (Zacarías-Eslava and del Castillo, 2010). More attention should be paid to edge effects in future studies to explore the possible role of edge effects on this result.

The negative effect of area on diversity obtained in our analysis suggests that we have not omitted any important diversity driver that is positively associated with fragment area, in which case, a positive, not a negative effect of area on diversity would be obtained. Our results, therefore, suggest that environmental factors associated with topography are among the most important diversity drivers for woody plants in the seasonally dry oak forest remnants of the Mixteca Region.

Both heterogeneity in altitude and heterogeneity in slope-aspect within the fragments showed a significant relationship with habitat heterogeneity. However, their combined effect on diversity effects are opposite, according to our SEM analysis. This result suggests that a greater effect of topographic heterogeneity on species diversity can be achieved with a combination of high heterogeneity in altitude with low heterogeneity in slope-aspect. Thus, these 2 variables should be considered together when analyzing the impact of habitat heterogeneity on diversity. The involvement of climatic effects affecting species diversity and associated with topography may help to interpret this result. The mean annual temperature, for instance, is well known to decrease linearly with altitude (e.g., Zacarías-Eslava and del Castillo, 2010). Throughout the same slope-aspect, a given mean temperature is expected to be found only at a unique elevation point, ignoring microclimatic differences caused by variations in shading by vegetation or micro topography. However, the same mean temperature can be found at different elevations on a mountain if the orientation of the slope changes. North-facing slopes, for instance, are usually colder than south-facing slopes, at the same elevation in the Northern Hemisphere. The same temperature that is found on the north face of a mountain at a given elevation is expected to be found at a higher elevation at other slope orientations. In this way, the combination of high heterogeneity in both slope-aspect and elevation may reduce the total environmental variation of the fragment because different combinations of altitude and orientation can render the same climate.

Disturbance is a factor that undoubtedly alters ecosystem biodiversity, even if there is no consensus on how it works (Mackey and Currie, 2000, 2001). Most studies have developed indices to assess the total effect of disturbance without distinguishing the partial effect of each source of disturbance on biodiversity (Ross et al., 2002). Stump density reveals logging activities in the fragments, whereas road-effect zone is an indicator of the accessibility of the fragment to anthropogenic activities (Forman et al., 1997), as well as potential invasion by exotic species. Only logging, as a disturbance indicator, was significant in our study. In the study site, logging is manual and selective. Logged trees are scattered over the landscape, creating small gaps in the fragments. In seasonally dry oak forests, each small opening may decrease soil moisture, creating inappropriate conditions for native sapling development (Asbjornsen et al., 2004a, 2004b; Brown et al., 2004). In turn, the new conditions in these open spaces may allow the establishment of resilient and short-lived species, such as pioneer species. The above leads to species turnover and an increased biodiversity, by allowing, to a certain extent the coexistence of pioneer and shade-tolerant species in the same fragment. The effect of roads on diversity in the study area was not significant, probably because the roads in the area are used primarily for communication between indigenous villages, which are characterized by very low population densities. As a result, the effect of human disturbance in this area is probably due mainly to logging, either for fuel or for the small-scale production of wood products for local construction or tool manufacturing.

Conservation implications. Our results provide evidence of the importance of selecting fragments with high variation in topographic heterogeneity, which may favor a great diversity of species. The kinds and intensities of disturbance are also crucial since they may not have a common effect on species. The consideration of the species to be preserved is also crucial in developing conservation strategies, since different species may have different requirements, and some strategies may benefit only a limited number of species and harm others.

Using data of 12 fragments and structural equation modeling techniques, we were able to test and confirm a 3 factor model that characterized the diversity of woody plant species of a seasonally dry oak forest in the Mixteca Alta, Oaxaca, Mexico. Topographic heterogeneity, human disturbance, and fragment size, in that order of importance, play a significant role on woody plant diversity of these fragments. Topographic heterogeneity had a positive relationship with diversity. Fisher's α diversity increased significantly with fragment heterogeneity in slope-aspect or altitude. Disturbance, here estimated as a fragment's stump density, also showed a positive relationship with diversity. Small fragments with similar levels of topographic heterogeneity and disturbance tend to be more diverse than large fragments, probably because small fragments have a greater perimeter-to-area ratio and therefore convey more opportunities for the successful establishment of species that are benefited by open, less dense habitats such as edge, resilient, and pioneer species. Structural equation modeling was shown to be an appropriate technique for disentangling the contribution of several factors related with biodiversity. Conservation strategies of fragmented landscapes must consider not only fragment size, but the type and intensity of disturbance affecting the fragments and the species that need to be preserved.

 

Acknowledgments

We thank Raúl Rivera for field and cartographic work, and Adriadna Ferrer, Soledad Iglesias, Remedios Martínez, Hugo Nieto, Adriana Pacheco, and Victoria Pérez for field work. Suzanna Elkin helped us with the English version. The project was funded by grants from the Instituto Politécnico Nacional (SIP and COFAA), and the European Commission INCO V Program (REFORLAN FP6-2004-INCO-DEV-3 032132). We acknowledge the valuable comments of María Luisa Martínez and two anonymous reviewers. We are especially grateful to the people of Santiago Huauclilla, Santiago Tilantongo, and San Antonio Nduayaco (Santiago Apoala municipal agency), Nochixtlán District villages, who allowed us to carry out this work in their municipalities.

 

Literature cited

Alelign, A. D., D. Teketay, Y. Yemshaw and S. Edwards. 2007. Diversity and status of regeneration of woody plants on the peninsula of Zegie, northwestern Ethiopia. Journal of Tropical Ecology 48:37-49.         [ Links ]

Arrhenius, O. 1921. Species and area. Journal of Ecology 9:95-99.         [ Links ]

Asbjornsen, H., M. A. Ashton, D. J. Vogt and S. Palacios. 2004a. Effects of habitat fragmentation on the buffering capacity of edge environments in a seasonally dry tropical oak forest ecosystem in Oaxaca, Mexico. Agriculture, Ecosystems et Environment 103:481-495.         [ Links ]

Asbjornsen, H., K. A. Vogt, and M. S. Ashton. 2004b. Synergistic responses of oak, pine and shrub seedlings to edge environments and drought in a fragmented tropical highland oak forest, Oaxaca, Mexico. Forest Ecology and Management 192:313-334.         [ Links ]

Aström, M., J. Dynesius, K. Hylander and C. Nilsson. 2007. Slope aspect modifies community responses to clear-cutting in boreal forests. Ecology 88:749-758.         [ Links ]

Balvanera, P. and E. Aguirre. 2006. Tree diversity, environmental heterogeneity, and productivity in Mexican tropical dry forest. Biotropica 38:479-491.         [ Links ]

Baz, A. and A. García-Boyero. 1995. The effects of forest fragmentation on butterfly communities in central Spain. Journal of Biogeography 22:129-140.         [ Links ]

Bernabe, N., G. Williams-Linera and R. M. Palacios. 1999. Tree ferns in the interior and at the edge of a Mexican cloud forest remnant: spore germination and sporophyte survival and establishment. Biotropica 31:83-88.         [ Links ]

Blakely, T. J. and R. K. Didham. 2010. Disentangling the mechanistic drivers of ecosystem-size effects on species diversity. Journal of Animal Ecology 79:1204-1214.         [ Links ]

Boecklen, W. J. 1986. Effects of habitat heterogeneity on the species-area relationships of forest birds. Journal of Biogeography 13:59-68.         [ Links ]

Bollen, K. A. 1989. Structural equations with latent variables. John Wiley and Sons. New York-Chichester-Brisbane-Toronto-Singapore. 514 p.         [ Links ]

Brose, U. 2001. Relative importance of isolation, area and habitat heterogeneity for vascular plant species richness of temporary wetlands in east-German farmland. Ecography 24:722–730.         [ Links ]

Brown, K. A., J. Gurevitch and R. Sokal. 2004. Long-term impacts of logging on forest diversity in Madagascar. Proceedings of the National Academy of Sciences of the United States of America 101:6045-6049.         [ Links ]

Brown, K. A., F. N. Scatena and J. Gurevitch. 2006. Effects of an invasive tree on community structure and diversity in a tropical forest in Puerto Rico. Forest Ecology and Management 226:145-152.         [ Links ]

Bunting, B. T. 1964. Slope development and soil formation on some British sandstones. Geographical Journal 130:73-79.         [ Links ]

Bustamante, R. and A. A. Grez. 1995. Consecuencias ecológicas de la fragmentación de los bosques nativos. Ambiente y Desarrollo 11:58-63.         [ Links ]

Cayuela, L., D. J. Golicher, J. M. Rey-Benayas, M. González-Espinosa and N. Ramírez-Marcial N. 2006a. Fragmentation, disturbance and tree diversity conservation in tropical montane forests. Journal of Applied Ecology 43:1172-1181.         [ Links ]

Cayuela, L., D. J. Golicher, J. Salas Rey and J. M. Rey-Benayas. 2006b. Classification of a complex landscape using Dempster-Shafer theory of evidence. International Journal of Remote Sensing 27:1951-1971.         [ Links ]

Chuvieco, E. 2006. Teledetección ambiental. 2da edn. Ariel S. A. Barcelona. 594 p.         [ Links ]

Clarke, A. and K. J. Gaston. 2006. Climate, energy and diversity. Proceedings of the Royal Society B: Biological Sciences 273:2257-2266.         [ Links ]

Davis, S. K. 2004. Area sensitivity in grassland passerines: effects of patch size, patch shape, and vegetation structure on bird abundance and occurrence in southern Saskatchewan. The Auk 121:1130-1145.         [ Links ]

di Bella, C. M., G. Posse, M. E. Beget, M. A. Fischer, N. Mari and S. Veron. 2008. La teledetección como herramienta para la prevención, seguimiento y evaluación de incendios e inundaciones. Ecosistemas 17:39-52.         [ Links ]

Echeverría, C., A. C. Newton, A. Lara, J. M. Rey-Benayas and D. A. Coomes. 2007. Impacts of forest fragmentation on species composition and forest structure in the temperate landscape of southern Chile. Global Ecology and Biogeography 16:426-439.         [ Links ]

Fahrig, L. 2003. Effects of habitat fragmentation on biodiversity. Annual Review of Ecology, Evolution and Systematics 34:487-515.         [ Links ]

Fisher, R. A., A. S. Corbet and C. B. Williams. 1943. The relation between the number of species and the number of individuals in a random sample of an animal population. Journal of Animal Ecology 12:42-58.         [ Links ]

Forero-Molina, A. and B. Finega. 2004. Efectos de borde en la vegetación de remanentes de bosque muy húmedo tropical en el norte de Costa Rica y sus implicaciones para el manejo y la conservación. Revista Forestal Centroamericana 38:39-43.         [ Links ]

Forman, R. T. T., D. S. Friedman, D. Fitzhenry, J. S. Martin, A. S. Chen and L. E. Alexander. 1997. Ecological effects of roads: toward three summary indices and an overview for North America. In Habitat fragmentation and infrastructure: Proceedings, K. Canters (ed.). Ministry of Transport, Public Works and Water Management. Delft, The Netherlands. 40-54.         [ Links ]

Forman, R. T. T. and L. E. Alexander. 1998. Roads and their major ecological effects. Annual Review of Ecology, Evolution and Systematics 29: 207-231.         [ Links ]

Freemark, K. E. and H. G. Merriam. 1986. Importance of area and habitat heterogeneity to bird assemblages in temperate forest fragments. Biological Conservation 36:115-141.         [ Links ]

García-Mendoza, A., P. Tenorio and J. Reyes. 1994. El endemismo en la flora fanerogámica de la Mixteca Alta, Oaxaca-Puebla, México. Acta Botanica Mexicana 27:53:74.         [ Links ]

Gaston, D. J. 2000. Global patterns in biodiversity. Nature 405:220-227.         [ Links ]

González-Leyva, A. 2007. Geografía, lingüística, arqueología e historia de la Mixteca alta antes de la conquista española. Universidad Nacional Autónoma de México. Facultad de Filosofía y Letras, Anuario de Historia 1:45-66.         [ Links ]

Grace, J. B. and B. H. Pugesek. 1997. A structural model of plant species richness and its application to a coastal wetland. The American Naturalist 149:436-460.         [ Links ]

Graham, C. H. and J. G. Blake. 2001. Influence of patch- and landscape- level factors on bird assemblages in a fragmented tropical landscape. Ecological Applications 11:1709-1721.         [ Links ]

Gucinski, H., M. H. Brookes, M. J. Furniss and R. R. Ziemer (eds.). 2001. Forest roads: a synthesis of scientific information. United States Department of Agriculture. Forest Service. USA. http://www.fs.fed.us/pnw/pubs/gtr509.pdf; last access: 28.IX.2011.         [ Links ]

Guerrero-Arenas, R., E. Jiménez-Hidalgo and H. Santiago-Romero. 2010. La transformación de los ecosistemas de la Mixteca Alta oaxaqueña desde el Pleistoceno tardío hasta el Holoceno. Ciencia y Mar 14:61-68.         [ Links ]

Hannus, J. J. and M. von Numers. 2008. Vascular plant species richness in relation to habitat diversity and island area in the Finish Archipelago. Journal of Biogeography 35:1077-1086.         [ Links ]

Hayduk, L. A. 1987. Structural equation modeling with LISREL. The Johns Hopkins University Press. Baltimore and London. No. de págs??        [ Links ]

Herrera, M. A., L. A. Olaya and F. Castro. 2004. Incidencia de la perturbación antrópica en la diversidad, la riqueza y la distribución de Eleutherodactylus (Anura: Leptdactylidae) en un bosque nublado del suroccidente colombiano. Caldasia 26:265-274.         [ Links ]

Hill, J. L. and P. J. Curran. 2001. Species composition in fragmented forest: conservation implications of changing forest area. Journal of Applied Geography 21:157-174.         [ Links ]

Hill, J. L. and P. J. Curran. 2003. Area, shape and isolation of tropical forest fragments: effects on tree species diversity and implications for conservation. Journal of Biogeography 30:1391-1403.         [ Links ]

Inegi (Instituto nacional de estadística y geografía). 2011. Conjunto de datos vectoriales. http://mapserver.inegi.gob.mx/geografia/espanol/prodyserv/bdgweb/wizzr.cfm; last access: 5.IX.2011.         [ Links ]

Laurance, W. F., H. L. Vasconcelos and T. E. Lovejoy. 2000. Forest loss and fragmentation in the Amazon: implications for wildlife conservation. Oryx 34:39-45.         [ Links ]

Laurance, W. F. 2004. Forest-climate interactions in fragmented tropical landscapes. Philosophical Transactions of the Royal Society B: Biological Sciences 359:345-352.         [ Links ]

Laurance, W. F. 2008. Theory meets reality: how habitat fragmentation research has transcended island biogeographic theory. Biological Conservation 141:1731-1744.         [ Links ]

Lezcano, H., B. Finegan, R. Condit and D. Delgado. 2004. Variación de las características de la comunidad vegetal en relación al efecto de borde en fragmentos de bosque Las Pavas, Cuenca del Canal de Panamá. Revista Forestal Centroamericana 38:33-38.         [ Links ]

López, T. M. 2001. Urban expansion and the loss of prime agricultural lands in Puerto Rico. Ambio 30:49-54.         [ Links ]

MacArthur, R. H. and E. O. Wilson. 1967. The Theory of Island Biogeography. Monographs in Population Biology. Princeton University Press, New Jersey. 203 p.         [ Links ]

Mackey, R. L. and D. J. Currie. 2000. A re-examination of the expected effects of disturbance on diversity. Oikos 88:483-493.         [ Links ]

Mackey, R. L. and D. J. Currie. 2001. The diversity-disturbance relationship: is it generally strong and peaked? Ecology 82:3479-3492.         [ Links ]

Magurran, A. E. 2004. Measuring biological diversity. Blackwell Publishing, Oxford. 256 p.         [ Links ]

Martínez, J. and M. Noriega (eds.). 2006. Manejo del agua y restauración productiva en la región indígena Mixteca de Puebla y Oaxaca. The International Bank for Reconstruction and Development/The World Bank. Washington, D. C.-México, D. F. 104 p.         [ Links ]

McGarigal, K. and B. Marks. 1995. FRAGSTATS: a spatial pattern analysis program for quantifying landscape structure v2.0 (unpublished computer program user manual and guide). Oregon State University. Corvallis. USA.         [ Links ]

Miller, A. D., S. H. Roxburgh and K. Shea. 2011. How frequency and intensity shape diversity-disturbance relationships. Proceedings of the National Academy of Sciences of the United States of America 108:5643-5648.         [ Links ]

Mulaik, S. A., L. R. James, J. Van-Alstine, N. Bennett, S. Lind and C. D. Stilwell. 1989. Evaluation of goodness of fit indices for structural equation models. Psychological Bulletin 105:430-445.         [ Links ]

Murcia, C. 1995. Edge effects in fragmented forests: implications for conservation. Trends in Ecology and Evolution 10:52-62.         [ Links ]

Otálora, M., I. Martínez, R. Belinchón, I. Widmer, G. Aragón, A. Escudero and C. Scheidegger. 2011. Remnants fragments preserve genetic diversity of the old forest lichen Lobaria pulmonaria in a fragmented Mediterranean mountain forest. Biodiversity and Conservation 20:1239-1254.         [ Links ]

Peterson, D. l., E. G. Scheriner and N. M. Buckingham. 1997. Gradients, vegetation and climate: spatial and temporal dynamics in the Olympic Mountains, U.S.A. Global Ecology and Biogeography 6:7-17.         [ Links ]

Pincheira-Ulbrich, J., J. R. Rau and E. Hauenstein. 2008. Diversidad de árboles y arbustos en fragmentos de bosque nativo en el sur de Chile. Phyton-International Journal of Experimental Botany 77:321-326.         [ Links ]

Physical Geography Dictionary. 2012. http://www.physicalgeography.net/physgeoglos/s.html; lats access: 10.XI.2012.         [ Links ]

Preston, F. W. 1962. The canonical distribution of commonness and rarity: parts I and II. Ecology 43:185-215 and 43:410-432.         [ Links ]

Reed, D. N., T. M. Anderson, J. Dempewolf, K. Metzger and S. Serneels. 2009. The spatial distribution of vegetation types in the Serengeti ecosystem: the influence of rainfall and topographic relief on vegetation patch characteristics. Journal of Biogeography 36:770-782.         [ Links ]

Ross, K. A., B. J. Fox and M. D. Fox. 2002. Changes to plant species richness in forest fragments: fragment age, disturbance and fire history may be as important as area. Journal of Biogeography 29:749-765.         [ Links ]

Rudas, G., D. Marcelo, D. Armenteras, N. Rodríguez, M. Morales, L. C. Delgado and A. Sarmiento. 2007. Biodiversidad y actividad humana: relaciones en ecosistemas de bosque subandino en Colombia. Instituto de Investigación de Recursos Biológicos Alexander von Humboldt, Bogotá D.C., Colombia. 128 p.         [ Links ]

Santos, T. and J. L. Tellería. 2006. Pérdida y fragmentación del hábitat: efecto sobre la conservación de las especies. Ecosistemas 15:3-12.         [ Links ]

SAS Institute. 1989. SAS/STAT User's Guide. Version 6, Fourth Edition, Volume 1, Cary, NC: SAS Institute Inc.         [ Links ]

Saunders, D. A., R. J. Hobbs and C. R. Margules. 1991. Biological consequences of ecosystem fragmentation: a review. Conservation Biology 5:18-32.         [ Links ]

SEMARNAT (Secretaría de Medio Ambiente y Recursos Naturales). 2005. Indicadores Básicos del Desempeño Ambiental de México 2005. http://148.223.105.188:2222/gif/snif_portal/index.php?option=com_contentettask=viewetid=49etItemid=63; last access: 18.I.2011.         [ Links ]

Spitale, D., A. Petraglia and M. Tomaselli. 2009. Structural equation modelling detects unexpected differences between bryophyte, and vascular plant richness along multiple environmental gradients. Journal of Biogeography 36:745-755.         [ Links ]

Stevenson, P. R. and M. E. Rodríguez. 2008. Determinantes de la composición florística y efecto de borde en un fragmento de bosque en el Garviare, Amazonia Colombiana. Revista Colombiana Forestal 11:5-17.         [ Links ]

Stoelting, R. 2002. Structural Equation Modeling/Path Analysis. http://userwww.sfsu.edu/~efc/classes/biol710/path/SEMwebpage.htm; last access: 20.X.2010.         [ Links ]

Teillet, P. M., B. Guindon and D. G. Goodeonugh. 1982. On the slope-aspect correction of multispectral scanner data. Canadian Journal of Remote Sensing 8:84-106.         [ Links ]

Tews, J., U. Brose, V. Grimm, K. Tielbörger, M. C. Wichmann, M. Schwager and F. Jeltsch. 2004. Animal species diversity driven by habitat heterogeneity/diversity: the importance of keystone structures. Journal of Biogeography 31:79-92.         [ Links ]

Tjorve, E. 2003. Shapes and functions of species-area curves: a review of possible models. Journal of Biogeography 30:837-835.         [ Links ]

Verburg, P. H., K. P. Overmass and N. Witte. 2004. Accessibility and land-use patterns at the forest fringe in the Northeastern part of the Philippines. The Geographical Journal 170:238-255.         [ Links ]

Wade, T. G., K. H. Riiters, J. D. Wickham and K. B. Jones. 2003. Distribution and causes of global forest fragmentation. Conservation Ecology 7:7. [online] http://www.consecol.org/vol7/iss2/art7/; last access: 20.X.2010.         [ Links ]

Williams-Linera, G., R. H. Manson and E. Isunza. 2002. La fragmentación del bosque mesófilo de montaña y patrones de uso del suelo en la región oeste de Xalapa, Veracruz, México. Madera y bosques 8:73-89.         [ Links ]

Young, K. R. 1994. Roads and the environmental degradation of tropical montane forests. Conservation Biology 8:972-976.         [ Links ]

Zacarías-Eslava, Y. and R. F. del Castillo. 2010. Comunidades vegetales templadas de la Sierra Juárez Oaxaca: pisos altitudinales y sus posibles implicaciones ante el cambio climático. Boletín de la Sociedad Botánica de México 87:13-28.         [ Links ]