nueva página del texto (beta)
Inglés (pdf)
Artículo en XML
Referencias del artículo
Traducción automática
Enviar artículo por email
Citado por SciELO 
Similares en
SciELO 
Permalink
Although patterns in ecological communities change over spatio-temporal scales, ecological theory has mostly focused on a single scale, assuming that communities are closed and stable entities over time (Leibold et al. 2004). However, ecologists are increasingly focusing on meta-communities. In this approach, communities are considered as open entities that maintain different levels of functional interconnection in a regional environment (Holyoak et al. 2005). Thus, metacommunities are defined as "... a set of local communities that are connected by dispersal of multiple species that can potentially interact" (Wilson 1992). The meta-community level allows us to incorporate large space-time scales to community ecology. For example, the meta-community approach can range from microarthropod communities in lichen patches (small scale), to studies across biogeographic provinces or continents (large scale) (Mouquet & Loreau 2003, Chase & Bengtsson 2010). On the other hand, metacommunity dynamics is apparently closely related to environmental filters (abiotic interactions), biotic interactions and to the dispersal of species between local communities (Holyoak et al. 2005, Puth & Post 2005, Hutchinson 2008, Logue et al. 2011). These factors change in space and time (Devlaeminck et al. 2005, Leibold & McPeek 2006). Therefore, it is important to consider the environmental heterogeneity between patches and the relatively unfavorable matrix in which to advance knowledge, conservation and management of ecosystems (Leibold et al. 2004, Datry et al. 2015).
Based on metacommunities theory, attempt have been made to analyze the elements of metacommunity structure (EMS) (Leibold & Mikkelson 2002) as an exploratory examination of possible changes in metacommunity structure (Keith et al. 2011). EMS are used to characterize the emergent properties of metacommunities using an incidence matrix with three measures: (1) coherence, or the degree of continuity of species across a major gradient; (2) turnover, or species replacement from site to site; and (3) boundary clumping, the extent to which the limits of the ranges of different species are clustered together (Leibold & Mikkelson 2002, Presley et al. 2010). Based on this approach, ecologists have described theoretical patterns for the distribution of species between patches in a metacommunity (Leibold & Mikkelson 2002) (Figure 1). These are: a) Nested subsets (positive coherence and negative turnover). There are fewer embedded absences and less species turnover than those expected by chance. Communities with few species form nested subsets of increasingly speciose communities, due to the tolerance of each species to abiotic conditions (Patterson & Atmar 1986). b) Checkerboards (negative coherence). There are significantly more embedded absences than expected by chance. Strong competition may result in this pattern due to species pairs with mutually exclusive distributions (Diamond 1975). c) Clementsian (positive coherence, positive turnover, and positive clumping). Discrete communities replace each other along gradients and share few species with other types of communities; coexistence depends on interaction with other species (Clements 1916). d) Gleasonian (positive coherence, positive turnover, and non-significant clumping). This pattern differs from the Clementsian, in that the grouping of boundaries is not different from chance. Species are distributed along environmental gradients independently of other species and coexistence results from casual similarities in requirements or tolerances (Gleason 1926). e) Evenly spaced gradients (positive coherence, positive turnover, and negative clumping). There is significantly a lower grouping of species’ range limits than expected by chance. This model was proposed by Tilman in 1982, the assumption is that there is strong competition between species, so the species ranges are evenly arranged along gradients with no discrete communities. f) Random (all EMS are indistinguishable from a random expectation). In this model there is no pattern in species distributions (Simberloff 1983).
Figure 1 Graphical representation of the possible combinations of EMS to differentiate between the six idealized structures (Leibold & Mikkelson 2002) and the associated quasi-structures (Presley et al. 2010). Modified from Presley et al. (2010).
Additional six quasi-structures have been proposed (Presley et al. 2010) (Figure 1): g) quasi-nested clumped species loss, h) quasi-nested stochastic species loss, i) quasi-nested hyperdispersed species loss, j) quasi-Clementsian, k) quasi-Gleasonian, and l) quasi-evenly spaced. These quasi-structures have the same characteristics as their associated basic structures, except that turnover does not differ significantly from a random expectation. For example, a nested metacommunity may exhibit loss of hyperdispersed species, with stochastic distribution or grouped between sites (Presley et al. 2010).
The EMS is based on an analysis of the distribution patterns of species along a major compositional-variation axis assumed to correspond to an important environmental or geographical gradient (Leibold & Mikkelson 2002, Leibold et al. 2004). An explicit consideration of environmental gradients and of the dispersal ability of different species between local communities would allow unraveling the role of stochastic and/or deterministic factors in structuring the metacommunity (Ozinga et al. 2004, Meynard et al. 2013, Denelle et al. 2019). A few previous studies have evaluated changes in the structure of metacommunities over time. Ellis et al. (2006) analyzed the structural changes in space and time (from 1978 to 2003) of the mosquito communities that inhabit water-filled holes in trees in southeastern Florida, USA, finding that temporal variation was greater than the spatial one. Keith et al. (2011) found in Dorset, England, that the structure of the plant metacommunity was Clementsian and that, despite the recorded loss of species, the metacommunity remained unchanged over 70 years.
There are few long-term studies that quantify changes in plant metacommunity structure in intertropical drylands (Miller et al. 2012, Stewart et al. 2014, Abella et al. 2019). These environments are important due to their high diversity and endemism of plant species. Approximately 60 % of the plants that grow in the semiarid areas of Mexico are endemic (Rzedowski 1978). In addition, arid and semi-arid plant communities are characterized by patchy spatial arrangements allowing for a straightforward metacommunity analysis (Wilson 1992, Aguiar & Sala 1994, Holyoak et al. 2005). In this arrangement, the spatial colonization between patches can be affected by geographical distance, interspecific dispersal differences, and environmental filters. These factors may be determinants of the structure, dynamics (MacArthur & Wilson 1967, Nathan & Muller-Landau 2000, Pueyo et al. 2008) and persistence of local communities (Holyoak et al. 2005, Hubbell 2005, Venable et al. 2008, Jones et al. 2015).
Zavala-Hurtado (1982) and Sandoval-Palacios (2010) analyzed the changes in composition and abundance of perennial plant species in 26 patches from three vegetation types (Table 1) over a 23-year period (1980-2003) in the semi-arid Zapotitlán valley, Mexico. Sites classified in 1980 as cardonal and shrubland did not show significant changes in their composition of species in 2003. However, 62.5 % (10 out of 16) sites occupied by tetechera changed to shrubland, suggesting that tetechera is unstable through time (Sandoval-Palacios 2010). These changes appeared to involve an “invasion” of one community by another and a trend towards a homogenization of the metacommunity (Sandoval-Palacios 2010). This led us to ask whether these changes have an effect on the structure of the spatial distribution of the species in the metacommunity, and how much does the environment and geographical distance influence the structure of the metacommunity.
Table 1 Three vegetation types recognized by Zavala-Hurtado (1982) and Sandoval-Palacios (2010) on the 1980 and 2003 historical surveys, respectively. The number of sampling units covered by each vegetation type on each survey and characteristic (dominant) species are indicated.
| Number of sampling units | |||
|---|---|---|---|
| Vegetation type | 1980 | 2003 | Characteristic species |
| cardonal | 7 | 7 |
|
| tetechera | 16* | 4 |
|
| shrubland | 4 | 15 |
|
* One of these patches was transformed to crop field. So it was removed from the RA.
Using the species’ presence-absence data from Zavala-Hurtado (1982) and Sandoval-Palacios (2010), we analyzed the structure of the metacommunity in the Zapotitlán valley. We expected the turnover and agglomeration of range limits would have decreased due to homogenization, resulting in a change from a Clementsian or quasi-Clementsian in 1980 toward a Gleasonian or quasi-Gleasonian structure in 2003. In addition, we expected to find changes in the effect of environmental (edaphic) and spatial (latitude, slope and geographical distance) gradients on the metacommunity structure. Thus, we expected that environmental and spatial variables at the local community level would explain, at least partially, the structure of the metacommunity. If the structure of the metacommunity in 2003 implied greater independence between the distribution ranges of the species, we would expect environmental filters to be of lesser importance that year.
Materials and methods
Study site. The study area is located in the semi-arid Zapotitlán valley in central Mexico (Figure 2), an intertropical local basin that covers an area of 413.89 km2 (located in 18° 11′ -18° 25′ N and 97° 39′ -97° 22′ W) in the Tehuacan-Cuicatlan Biosphere Reserve (Zavala-Hurtado 1982). The altitude ranges from 1,280 to 2,720 m asl. (Zavala-Hurtado 1982). The average temperature is 28 °C and the average annual rainfall is 380 mm (Zavala-Hurtado 1982). Soils are rocky and shallow, derived mainly from sedimentary and metamorphic rocks. The valley is mainly covered by tropical xerophytic scrub with a physiognomic predominance of columnar cacti (Rzedowski 1978, Valiente-Banuet et al. 1991).
Figure 2 Distribution of 26 sampling sites in the semiarid Zapotitlán valley (color squares), where the studied metacommunity was defined (Zavala-Hurtado 1982, Sandoval-Palacios 2010). The color on the left side of each square represents the type of community registered in the 1980 survey and the color on the right side to that of 2003: red = cardonal; green = tetechera; blue = scrubland.
Historical vegetation surveys. Through 1980, Zavala-Hurtado (1982) carried out a characterization of the vegetation of the Zapotitlán valley from 30 sampling sites. The selection of the sites was preferential and tried to cover different physiognomic modalities of vegetation accessed from extensive routes throughout the valley. The location of each site was marked on a topographic chart and referenced in field notebooks, and the locations of roads, towns and outstanding geomorphological features in the vicinity of each site were identified. A 50 × 5 m transect was chosen at each site and all the individuals of perennial plants present in each transect were registered and identified. The consideration of only perennial species allowed the sampling to be carried out throughout the year. In this survey, 85 perennial plant species were recorded at the sampling sites (Zavala-Hurtado 1982). Multivariate classification analyses (cluster analysis using Ward’s method and multiple discriminant analysis) allowed the characterization of three main types of vegetation in the Zapotitlán valley (Table 1). Cardonal, characterized by the presence of the giant unbranched columnar cactus, Cephalocereus columna-trajani. Tetechera, dominated by the giant branched columnar cactus, Cephalocereus tetetzo. Shrubland, dominated by thorny trees and shrubs, legumes, agaves, and small cacti (Zavala-Hurtado 1982).
In 2003, 23 years later, Sandoval-Palacios (2010) returned to the original location of 28 of the original 30 sampling sites because two sites were inaccessible at that time. Furthermore, two of these 28 sites were eliminated from the study because they had been completely transformed and lacked vegetation cover. So, the second survey included only 26 of the original sites. In these stands, the same sampling procedures of the environmental and biological variables were repeated as in 1980 (Sandoval-Palacios 2010). A total of 104 perennial plant species was recorded in 2003 and the same three types of vegetation as described in 1980 were recognized by using the same multivariate techniques.
In this work, to verify the consistence of the three recognized types of local communities, we performed another multivariate classification analysis (cluster analysis, using Ward’s method) on the joint incidence matrix of the two surveys (1980 and 2003). Additionally, we performed a principal component analysis (PCA) of the correlation matrix (Kovach 2013) on the same joint incidence matrix to visualize the magnitude of the compositional changes in the local communities after 23 years. Multivariate analyzes for this work were performed using the MVSP software (version 3.22, Kovach Computing Services, Pentraeth, Wales, UK).
Spatial variables. Zavala-Hurtado (1982) in the 1980 survey, at each sampling site and at the middle point of each transect, recorded the altitude (m asl.) using an analog Thommen Classic altimeter, and the slope inclination (º) using a Brunton Pocket transit compass. These measurements were repeated in 2003 (Sandoval-Palacios 2010) to verify the location of the sampling sites. In 2003, the altitude (m asl.) and the location of the sampling sites (UTM WGS84) were recorded using a Garmin 12CX GPS and the slope of terrain (º) using the same Brunton Pocket transit as that used in 1980. As a proxy for stand’s isolation, related to dispersal limitation (Nekola & White 1999, Moritz et al. 2013), we estimated the distance (km) between stand from the longitudinal and latitudinal UTM coordinates for each site using BaseCamp 4.7.0 (Garmin 2018).
Environmental variables. During the two historical surveys, at 5 m from the longitudinal edges and in the center of each transect, soil depth was measured, and three soil samples were taken at a depth of 15 cm. These samples were mixed to obtain a composite sample for each transect, in order to evaluate differences between sites. The variables texture, pH, CII, organic matter and the nutrients Na, K, P, Ca and Mg were quantified for each sample in the Edaphology Laboratory of the National Institute of Research on Biotic Resources (currently Institute of Ecology, AC) for the 1980 survey and in the Soils Fertility Laboratory of the University of Chapingo for the 2003 survey. Standard extraction and quantification techniques were used as described in Sandoval-Palacios (2010).
Metacommunity structure analysis. There are several methods for ordering an incidence matrix so that the arrangement obtained preserves the similarities and differences in the ranges of species distribution and in the composition of species between sites (Leibold & Mikkelson 2002). We used a method that evaluates the ordination of species and sites in gradients extracted by using a multivariate analysis of reciprocal averaging (RA), because it leads to a stable solution that does not depend on the initial arrangement of sites. Additionally, the ordination obtained by this method is not distorted by involution, at least in the first axis of variation, compared to other techniques such as PCA (Gauch et al. 1977). Furthermore, it allows recognition of joint variations in the structure of the metacommunity and environmental factors since ecological ordination seeks an understanding of gradient relationships between communities through the ordering of data along implicit environmental gradients (Whittaker 1967, Austin 1976).
The structure of the metacommunity was evaluated using presence-absence data of all the species recorded in the 26 patches sampled in 1980 and 2003. We used the metacom package V.1.5.1, from Dallas (2014), in R code V.3.6.1 (R core team 2019). This software provides the coherence, turnover and boundary clumping functions for the analysis of the elements of the metacommunity structure within a framework proposed by Leibold & Mikkelson (2002) and Presley et al. (2010). These functions are evaluated from the ordination of the species and sites in gradients extracted by means of a reciprocal averaging (RA) multivariate analysis (Gauch et al. 1977) on incidence matrices. This maximizes the proximity of species with similar distributions and the proximity of sites with compositions of similar species. The null model used to assess the significance of coherence and species’ turnover was the “r1”, which is the recommended one when species are distributed along environmental or spatial gradients, and because it is less sensitive to type I and type II errors (Presley et al. 2010). This model preserves the species richness of a site (row totals) and fills the species ranges (columns) based on their marginal probabilities. For each analysis, 1000 null matrices were simulated obtaining the ordination of species and samples on the first RA axis. For the boundary clumping analysis, metacom calculates the Morisita Index (Morisita 1962) on the presence-absence matrix, using a chi-squared test for significance. A significant value > 1 indicates that the limits are clumped, and a significant value < 1 indicates that the limits are not grouped. Thus, the idealized distribution patterns (nested, Clementsian, Gleasonian, evenly spaced, random, and checkerboard patterns, as well as the quasi-patterns) within the metacommunity of the study were determined for both 1980 and 2003. We compared the metacommunity structures between the two years to assess possible temporal changes.
Relationship between the structure of the metacommunity and environmental and spatial variables. The EMS analysis allows one to describe the structure of the metacommunity in a domain delimited by an indirect gradient defined by the first RA axis. However, this analysis does not permit an identification of the factors that could be responsible for generating the observed patterns. To analyze the effect of each factor (environmental, spatial, environmental + spatial, and those not explained) in the metacommunity structure, we used a variance partition approach by applying partial canonical correspondence analysis (PCCA) (Borcard et al. 1992). This is a multivariate equivalent to partial linear regression (Borcard et al. 2011) and applies the varpart function of the vegan package V.6.1 for R V.3.6.1 (R core team 2019). To assess the significance of particular variables within each explanatory matrix (environmental and spatial), we used canonical correspondence analysis (CCA) (Borcard et al. 2011). Before performing the PCCA and CCA analyzes, a multicollinearity test was applied to the soil variables to avoid increases in variance due to a possible correlation between the variables. For this test we used the pairs function in the car package v.3.0-3 for R v.3.6.1 (R core team 2019).
Results
Temporal changes in local community types. The results of the classification (Figure 3A) and ordination analyses (Figure 3B) on the joint matrix of the two temporal surveys are consistent with the previous analyses made with the matrices for each separate survey date (Zavala-Hurtado 1982, Sandoval-Palacios 2010). There is a clear distinction between the cardonal communities and those of the tetechera and scrubland in that although they form different groups, some of the sites appear intermingled in the boundary between these community types in the ordination space. The characterization of these three types of community is based on the recurring presence of representative species characterized from their scores on the first two axes of the PCA (Figure 3B): cardonal (Cephalocereus columna-trajani, Mascagnia seleriana, Jatropha dioica, Agave kerchovei, Euphorbia cymbifera, Bursera arida, Hechtia liebmannii); tetechera (Cephalocereus tetetzo, Lantana camara, Varronia curassavica, Plumeria rubra, Opuntia tomentosa); and scrubland (Prosopis laevigata, Myrtillocactus geometrizans, Parkinsonia praecox, Gymnosperma glutinosum, Opuntia pilifera, Celtis pallida, Lemaireocereus hollianus, Dalea carthagenensis, Agave marmorata). It also can be seen that the cardonal sites, though were displaced within the territory of this type of community, maintaining their identity in the two temporal surveys. In contrast, most of the sites originally identified as tetechera were reclassified as scrubland, and those that remained identified as tetechera in 2003 were now placed on the border between these two types of community.
Figure 3 Dendrogram from a Ward’s method Cluster analysis (A) and Principal Components Analysis (B) on a 1980-2003 joint incidence matrix (104 species in 26 × 2 stands) of a plant metacommunity in the semiarid Zapotitlán valley, Mexico. Symbols in the dendrogram are: diamonds = 1980 survey; squares = 2003 survey. Different colors indicate the previously assigned type of community (Zavala-Hurtado 1982, Sandoval-Palacios 2010): red = cardonal; green = tetechera; blue = scrubland.
Metacommunity structure. When considering the two survey times in the light of an analysis of the three metacommunity structural elements, there were significantly fewer embedded absences than those anticipated by the null model; that is, there was positive coherence. The species’ turnover was significantly greater than that expected from the null model; that is, there was no nested structure. Finally, the metacommunity showed significant clumping of species’ ranges, the values of the Morisita Index being > 1 and different from those expected from chance (Table 2). So, the metacommunity showed a Clementsian structure that had not changed after 23 years (Figure 4). Nevertheless, the individual elements did show some changes. While coherence and range boundary clumping showed decrement, turnover had increased 23 years after the first survey (Table 2).
Figure 4 Reciprocal Averaging ordination of incidence matrices: (A) 84 species in 26 stands (1980), and (B) 104 species in 26 stands (2003) of a plant metacommunity in the semiarid Zapotitlán valley, Mexico. Colored squares represent incidence of species in each stand. Different colors indicate previously assigned type of community (Zavala-Hurtado 1982, Sandoval-Palacios 2010): red = cardonal; green = tetechera; blue = scrubland.
Table 2 Results of the analyses of coherence, species turnover and boundary clumping for the plant desert metacommunity of the Zapotitlán valley in two historical surveys (1980 and 2003). Metacommunity elements were extracted from the primary ordination axis (reciprocal averaging, RA). Mean and (SD) values were estimated from 1,000 generated null matrices, based on the r1 null model. EA = embedded absences, SR = species replacements, MI = Morisita Index. Significant results (p < 0.05) are in bold.
Structuring mechanisms: environmental and spatial factors. Of the 12 soil variables initially considered as possible structuring factors of the metacommunity, those that were severely multicollinear (variance inflation factor > 10) and those that showed a non-significant and low correlation with the first CCA axis were eliminated. This was done for both survey times, and the number of soil variables entered into the analysis was six in 1980 and four in 2003. None of the three spatial variables showed multicollinearity. The variation partition analysis used to assess the independent effect of spatial and environmental gradients on the spatial structure of the metacommunity for the 1980 survey indicated that the environmental effects, after removing the spatial effects, explained 15.1 % of the total variation. The effects associated with space, without the environmental effects, accounted for 8.2 % of the total variation. The variation shared by the effects of environmental and spatial variables was 0.07 %. The unexplained variation (the residuals) was 76.5 %, and this represented the un-assessed effects (of chance and noise) that do not depend on the environment, on the space, or on their joint effect (Figure 5A; Table 3). For the 2003 survey, the variation explained by spatial factors was 9.2 %, whereas the contribution of the environmental filters was 9.3 %. The variation explained by environmental and spatial effects was negative (-1.8 %), and such cases are regarded as zero values (Borcard et al. 2004). The variance not explained by any of the gradients evaluated (environmental and spatial) was 83.3 % (Figure 5B; Table 3).
Figure 5 Venn diagrams showing variance partitioning of species distribution in the metacommunity in 1980 (A) and 2003 (B) surveys with environmental (soil) and spatial explanatory variables using partial CCA analyses. X1 = environmental (soil) matrix, X2 = spatial matrix. The overlapped area between the two circles is the variance explained jointly by the two matrices.
Table 3 Correlation of the species distribution in the metacommunity with the environment (soil) and spatial variables using partial CCA analyses. Zero values are regarded as R < 0. NT = no testable. Significant results (p < 0.05) are in bold.
| PCCA 1980 | PCCA 2003 | |||
|---|---|---|---|---|
| Presence-absence matrix vs. | Adj. R (p value) | Presence-absence matrix vs. | Adj. R (p value) | |
| Environment | 0.15 (0.001) | Environment | 0.093 (0.008) | |
| Spatial | 0.08 (0.001) | Spatial | 0.092 (0.002) | |
| Environment-Spatial | 0.0006 (NT) | Environmental-Spatial | 0 (NT) | |
To identify the individual variables that influenced the results described above, we assessed the correspondence between the species incidence matrix and the environmental and spatial matrices in separate CCAs for each survey. In the 1980 survey, five environmental variables showed a significant correlation with the first CCA axis: Na, K, Mg, P and organic matter (Table 4). In contrast, at the second survey (2003), only two environmental variables were significantly correlated with the first canonical component, Ca and P (Table 4).
Table 4 Canonical coefficients and the intraset correlation of environmental and spatial variables with the first two axes of canonical correspondence analysis (CCA). Significant results (p < 0.05) are in bold.
| CORRELATION COEFFICIENTS | |||
|---|---|---|---|
| 1980 | Axis 1 | Axis 2 | p |
| Environmental variables | |||
|
|
|
|
|
| Spatial variables | |||
|
|
|
|
|
| 2003 | |||
| Environmental variables | |||
|
|
|
|
|
| Spatial variables | |||
|
|
|
|
|
In 1980, three spatial variables were significantly correlated with the first of the CCA axes: slope, latitude, and longitude. In 2003, only the spatial variables latitude and longitude showed significant correlation with the first canonical axis (Table 4).
Discussion
Temporal changes in local community types. After combining the samples from 1980 and 2003, we found a pattern consistent with previous analyses performed separately that clearly separated the three recognized types of local communities (shrubland, cardonal and tetechera) (Zavala-Hurtado 1982, Sandoval-Palacios 2010). Although the sites identified as cardonal remained as such over time (Sandoval-Palacios 2010), there were dramatic changes in the tetechera community, suggesting a scrubland encroachment into lands previously occupied by tetechera. This change could be driven by anthropogenic disturbances that involve the loss of species characteristic of the tetechera, particularly the giant columnar cactus Cephalocereus tetetzo, and thus would constitute a stage in a long-term successional process (Sandoval-Palacios 2010).
Metacommunity structure. The plant metacommunity of the Zapotitlán valley showed a Clementsian structure in the 1980 survey, which was maintained 23 years later, despite the documented changes in the composition and abundance of the species present in local communities. For this inter-tropical semi-arid region it is the first work on the dynamics of the structure of plant metacommunities and, as far as we know, also for Mexico, where studies on the dynamics of the structure of metacommunities have only been developed for bats (López-González et al. 2012) and amphibians (Ochoa-Ochoa & Whittaker 2014). This Clementsian arrangement is consistent with the discrete segregation of the three previously defined types (Table 1) for the 26 local communities in their ordination along the RA gradient with the sequential ordering of cardonal-tetechera-shrubland sites (Figure 4). Although there was no temporary change in the Clementsian structure, the coherence and grouping of the species distribution limits showed a decrease, although not statistically significant (Table 2). However, this allows us to suggest that the decrease in these two elements of the metacommunity structure is possibly an ongoing trend towards homogenization. This can be seen in the RA ordination diagram in 2003 (Figure 4B), where two local communities determined as shrubland were interspersed among the tetechera communities and where the number of shrubland increased. However, caution should be exercised with these assumptions because we only have the data of two discrete points in time. In order to get a more complete idea of the dynamics of these patches, it is necessary to consider more points over time.
According to the theory of metacommunities (Leibold & Mikkelson 2002, Presley 2010), the Clementsian structure of the metacommunity of the Zapotitlán valley implies that it has a coherent spatial distribution, with replacement of species between local communities and with grouping of the range limits of distribution of the species. That is, it is organized in different discrete types of communities and the constituent species respond similarly to the biotic and abiotic environmental factors of each local community (Clements 1916).
Other studies that have evaluated the dynamics of structure over time have also found temporal stability (Ellis et al. 2006, Azeria & Kolasa 2008, Keith et al. 2011, Newton et al. 2012, Erös et al. 2014). However, only the work of Keith et al. (2011) and Newton et al. (2012) were related to plant communities (Southern England forest and Dorset calcareous grasslands, respectively). They also found a Clementsian structure that remained stable after 70 years, despite a loss in diversity recorded in that time interval. Although these are different ecosystems from the one we studied, these stability results are consistent with ours. Apparently, although there are temporary changes in the composition and abundance of species in local communities, both in plants (Keith et al. 2011) and in other groups of organisms (Ellis et al. 2006, Azeria & Kolasa 2008, Erös et al. 2014), the different types of communities within the metacommunity are maintained. In the case of the Zapotitlán valley, the structure of the plant metacommunity did not show changes after 23 years, although there were dramatic changes in the composition and relative abundance in some of the local communities. Thus, the three types of communities (shrubland, tetechera and cardonal) that make up the meta-community persisted during the time covered by the study, although the number of tetechera communities decreased and the number of shrubland communities increased (Sandoval-Palacios 2010).
Structuring mechanisms: environmental and spatial factors. The variance partition of the effect of the environmental gradient (edaphic variables) and spatial (distance between patches), suggest that in 1980 the environmental filter had a greater effect than the spatial one as structuring mechanisms of the metacommunity (Figure 5A). This is in accordance with the Clementsian pattern for this year. For the 2003 survey, the environmental and spatial gradients explained the same percentage of variability (Figure 5B), while the Clementsian structure was maintained. This suggests that, in the plant metacommunity of the Zapotitlán valley, both factors would be important in the dynamics of species distribution. Our results are similar to those of Qian & Shimono (2012), who found, for a metacommunity of Tibetan alpine meadows in China, that the distribution of species correlated with the joint effect of geographical and environmental distance, which explained almost the same percentages of variability (4.8 and 3.5 %, respectively). It has also been reported that in the Amazon rainforest both geographical distance and environmental factors contribute in almost equal proportion to the assembly of local communities (Tuomisto et al. 2003). In contrast to what we found in the studied semi-arid environment, Hassler et al. (2010), in a study on the savannah in northwestern Namibia, did not find a relationship with either geographical distance or with environmental gradients. This suggests that deterministic and stochastic factors are not exclusive or static, and that at different points in time and space one or the other may have greater influence (Chase & Myers 2011). Additionally, they can act together to determine patterns of distribution of species at the level of metacommunities, as seen in the Zapotitlán valley.
Because it is very difficult to measure dispersal directly, one way to assess the effect of dispersal limitation on the metacommunity structure is through geographical distance between communities, which can be used as a proxy of connectivity between patches and can give an idea of the limits of regional dispersal (Nekola & White 1999, Hassler et al. 2010, Moritz et al. 2013). That is, closer patches have a greater chance of exchanging species than farther patches (Wang et al. 2011, Siefert et al. 2012, Qian & Shimono 2012, Moritz et al. 2013). Dispersal limitation has been recognized as a fundamental process in the dynamics of metacommunities (MacArthur & Wilson 1967, Hubbell et al. 1999, Nathan & Muller-Landau 2000, Jones et al. 2015). When considering the dynamics of semiarid metacommunities, Shreve (1929) and Shreve & Hinckley (1937) suggested that observed changes in the structure of local Arizona communities were due only to species movement between patches. This would indicate a trend toward a more uniform representation of species throughout a study area. Thus, the spatial and temporal changes observed in the Zapotitlán valley, given by increased species’ frequency/incidence in certain local communities, might be due to the arrival of plant species from other nearby patches or from the seed bank. However, environmental filters must be overcome because the patches within the Zapotitlán valley are immersed in a relatively unfavorable matrix. So the species must be able to reach new patches and fall into suitable sites for their establishment (Silvertown & Wilson 1994). In Zapotitlán, as well as in other arid environments, it has been found that these favorable sites often tend to occur under the canopies of shrubs that generate a microenvironment, where the radiation, temperature and humidity conditions are tempered, compared to open spaces (Valiente-Banuet et al. 1991). It has also been reported that under these shrubs there is greater accumulation of nutrients in the soil (Franco & Nobel 1989, Valiente-Banuet & Ezcurra 1991, Perroni-Ventura et al. 2010). Although it was not the objective of this work, we suggest that the conceptual paradigm of metacommunities (Holyoak et al. 2005) that could be adjusted to the semi-arid metacommunity of the Zapotitlán valley would be the ‘species sorting’ one. This, because the metacommunity has an apparently stable Clementsian structure in the analyzed time window. In addition, environmental and spatial filtering showed a small but significant effect on the structure of the metacommunity.
Finally, we found no evidence to support our hypothesis that the influence of changes in the composition and abundance of species at the level of local communities could affect the meta-community structure, tending towards homogenization. The Clementsian structure with discrete local communities was maintained over time. Our results also did not support the hypothesis about the importance of structuring mechanisms. From the recognized Clementsian distribution, we expected that the soil characteristics would have more influence than the distance between the patches; however, both gradients explained almost the same proportion of variability. This shows us that the Zapotitlán valley has a Clementsian structure that did not change despite changes in species diversity in local communities after 23 years. This suggests that the Clementsian structure is robust to local changes. Besides, it seems that there are multiple mechanisms and processes that influence the assembly and structure of metacommunities, since the gradients we consider in this work only explain 23 % (1980) and 18 % (2003) of the total variability. Many other biotic and abiotic factors must be acting in the dynamics of the communities. Because of this, multiple environmental gradients need to be considered, including temporary variations in climatic and spatial factors, as well as longer-term monitoring at the study sites.
