Abstract

This research examines the diversity-productivity relationship in a semiarid scrubland, initially under late successional conditions and subsequently under early successional conditions created by experimental clearing, to explore the roles that productivity and stochastic mortality play in species exclusion in this environment. A total of fifteen plots were studied by measuring environmental conditions and biomass components of shrubs and seedlings. These stands were distributed along a productivity gradient across five different landforms. A hypothesis about the stochastic self-thinning mortality model along the gradient was evaluated with the diversity-productivity-environment data. The diversity-productivity relationship was linear and reversed between the early and late succession stages. The hypothesis of stochastic mortality of species exclusion was rejected in the early stages of succession and partially accepted in the mature stage of succession. Species exclusion was negatively related to productivity gradients, suggesting that strong interspecific competition occurs in high productivity plots and that a larger number of species can survive in higher abiotic stress landscapes. Further research is needed to understand the temporal and spatial variations of the ecological interactions that shape this plant community.

1. Introduction

Environmental gradients offer ideal conditions for exploring the interactions between evolutionary adaptations and ecological processes over short spatial distances. Productivity [1], temperature [2], evapotranspiration [3], and the energy and water balance [2] are potential drivers of species diversity that can vary significantly over short distances, within the dispersal range of a species pool. However, local regional peculiarities such as fires, land use, or drought may obscure interpretation of the response of species diversity to environmental gradients.

Productivity is probably the single most important and most studied driver of species diversity because it is a function of several other critical factors such as the energy and water balance, evapotranspiration, and temperature [13]. Productivity is defined as the amount of plant mass produced over a certain time in a given area. Productivity is influenced by genetic adaptations and site environmental factors. Plant productivity is technically defined as net primary productivity or NPP, which is the difference between photosynthesis (gross primary productivity) and the total respiration of plants in an ecosystem. Both components are difficult to measure in the field at the spatial scales of interest [4]. Therefore, the net plant growth has been considered an estimator of NPP at short temporal scales where recruitment and mortality are absent.

Net primary productivity typically changes over time as plants increase in size, and species composition changes over the course of succession [5]. A linear relationship between productivity and species diversity has been reported for grasslands in several regions. However, many studies across a wide range of productivity conditions find a unimodal or hump-backed pattern in which diversity increases with productivity, reaches a maximum at intermediate levels of productivity, and decreases in highly productive sites [1, 57]. Abrams [8] pointed out that there are good theoretical reasons to predict a variety of differential relationships between diversity and productivity. Mittelbach et al. [1] summarized 171 observations on diversity-productivity relationships and found the unimodal pattern was the most common, followed by the monotonic and the -form models. Grime [7], Huston [5], Rosenzweig [9] and Waide et al. [10] indicated that the hump-back pattern is frequently observed in multiple taxas and at several scales of measurement.

The relationship between diversity and productivity has been studied primarily in grasslands or ruderal vegetal communities [6, 11, 12] because these short-lived species provide a quick response to environmental changes. Mittelbach et al. [1] recorded only five studies using shrubs and trees exclusively. Mittelbach et al. [1] suggested the need for observing the diversity-productivity relationship under field conditions in different ecosystems in order to understand and model the relationship more accurately. This would help to provide conceptual bases for the sustainable management and conservation of forest ecosystems.

The mechanisms that create the diversity-productivity patterns are not fully resolved. It seems that diversity increases as productivity increases, and this trend grows monotonically or decays in more productive sites, in agreement with the unimodal model. Along many productivity gradients, there is a clear tendency for diversity to diminish under highly productive conditions. That is, sites with intermediate productivity have greater species diversity than more productive sites [6, 7, 11, 13].

Two hypotheses have been put forth to explain the species diversity reduction where productivity is high: competitive exclusion by inter-specific competition (CEI), and the −3/2 law of the self-thinning mechanism at the level of the species assembly (ST) [11, 14]. In the first mechanism, inter-specific competition for limited resources promotes species exclusion [6, 7, 13]. In the second mechanism, the species disappearance is explained by random processes associated with the stand density reduction as individual size increases with time. In this stochastic mortality process, rare species have a probability of disappearance equal to individuals of more abundant species, and so are more likely to disappear because of their small population size. These effects of productivity on species diversity can reverse completely depending on the relative levels of disturbance and growth rates [13].

In spite of the extensive management of the Tamaulipan thornscrub plant community for ranching and forestry practices since early 1700’s [15]; contemporary land use changes [16] are: (a) modifying the local hydrology of this ecosystem [17, 18], (b) increasing the risk of extinction of several endemic and already endangered flora and fauna species [19, 20], and (c) contributing to potential global climate change by releasing carbon stocks contained in soils and biomass [16]. Due to high contemporary deforestation rates by gas mining in the region, restoration efforts of the plant community are underway [20] but information is required on the spatial and temporal changes of tree diversity and what may be the potential drivers of these shifts. Therefore, this research had the following objectives: (i) to determine the relationship between diversity and productivity in the early and mature stages of succession of plant communities in the Tamaulipan thornscrub plant community, (ii) to identify the possible mechanisms that control these relationships, including testing the hypothesis of stochastic mortality, and (iii) to understand the spatiotemporal patterns of these relationships and the mechanisms that may cause them to change over time.

2. Materials and Methods

2.1. Study Area

This research was conducted in the Tamaulipan thornscrub forest reserve of the Faculty of Forest Sciences, UANL, located in the city of Linares, N. L. in the Coastal Plain of the northern Gulf of Mexico. This physiographic region spans an area of 200,000 km2 in northeast Mexico and the southern United States. In northeastern Mexico, this ecosystem is in a state of transition due to rapid land use changes [16, 21]. The study area is located at 360 masl and is characterized by a warm, semi-arid, subtropical climate with an annual average precipitation of 745 mm (standard deviation, 250 mm) over an average of 80 rain events, of which 5% exceed 20 mm of rainfall. The annual average temperature is 22°C [22, 23].

The soils are derived from the Sierra Madre Oriental mountain range, which is composed of limestones and other sedimentary rocks uplifted by the orogeny that began near the end of the Upper Cretaceous and continued into the Paleocene period. These movements occurred simultaneously with the extension of the coastal plain of the present Gulf of Mexico. In general, the soil formation is influenced by its topographic position. The soil texture varies from loamy-sandy to clayey soils, with a wide range of colors, from gray to greenish and yellowish brown [24].

The Tamaulipan thornscrub plant community is diverse and is found on a wide variety of landforms. The floodplains in the lowest parts of the landscape are characterized by their proximity to intermittent or ephemeral streams, deep soils, and high availability of soil moisture during and shortly after the rainy season. Here, plants are characterized by larger size than anywhere else in the Tamaulipan thornscrub community. Following the gradient upslope, most of the landscape consists of valleys or open plains with deep to moderately deep clay soils, good fertility, and vegetation of medium height. Hills are conspicuous in the landscape and have high slopes and are usually crowned by mesas or plateaus that have shallow soils and low fertility levels. Finally, the elevated mesas or plateaus appear as isolated components of the landscape; they display shallow, thin, and stony soils, with scarce, short plant cover.

In general, the plant cover is characterized by species such as Ebenopsis ebano, Zanthoxylum fagara, Acacia farnesiana, Cordia boissieri, Caesalpinia mexicana, and Celtis pallida. These taxa have high abundance and cover [25]. The vegetative cover is dominated by shrubs (72.8%) and herbaceous plants (27.2%) [19]. Historically, the vegetation of the study area has been managed with conventional pasturing, as the land was extensively used for cattle ranching. However, the study area has been protected from this practice for approximately 30 years.

2.2. Sampling

Sampling consisted of establishing fifteen plots of , with three located on each of five landforms: (1) intermittent floodplain, (2) ephemeral floodplain, (3) valleys, (4) high-sloped terrain, and (5) mesas or plateaus. Replicates on each landform were separated by a distance of at least 20 meters. Therefore, the experimental design was stratified by landform with replicate plots positioned randomly within each landform. The five landforms created a productivity gradient within a linear distance of 1.5 km and an altitudinal range of 45 m.

In each plot, data collection for estimating diversity consisted of measuring basal diameter, top height, and cover for each shrub present in the stand (by measuring the vertical projection of the canopy). Total aboveground biomass was estimated by harvesting all existing vegetation, freshly weighed for each individual shrub; samples of biomass components (foliage, branches, and boles) were taken, freshly weighed and oven-dried, and the ratio of dry/fresh samples was used for estimation of the dry matter by biomass component by shrub, by species and finally by plot. Aboveground productivity was estimated by giving the harvested plots a year to grow seedlings, either by seed or by sprouting, and then repeating the same procedure for estimating the diversity and seedling biomass. With this experimental design, information was first collected on tree diversity and aboveground biomass for the later stages of succession and, later following a year of regrowth after complete harvesting, on the biomass, diversity and productivity for the early stages of succession. Total aboveground biomass was correlated with productivity; henceforth, an estimate of productivity for the later stages of succession was derived. In this manner, potential changes in the relationship between diversity and productivity over time were determined.

Other measured productivity correlates included soil nitrogen content, soil organic matter, and soil organic carbon in the upper 30 cm of soil. Organic matter content was measured by wet combustion according to a modified Walkley/Black methodology [24]. Soil nitrogen content was estimated by means of the digestion of total nitrogen according to the procedure of Kjeldahl. The soil texture was also measured by the Bouyoucos methodology.

2.3. Statistical Analysis

A correlation analysis was conducted between diversity ( = species richness), density (Den = density), and environmental (Text = texture, N = nitrogen content, CC = carbon content, and CN = carbon/nitrogen ratio) variables for both early and late successional stages. The diversity-productivity relationship was evaluated by fitting quadratic, linear, and second-degree polynomial regressions, where the aboveground productivity, , is the independent variable that predicts species diversity, , for example, []. The probabilities of the statistical parameters , , and define which is the best model. For example, the statistical significance of , , and can indicate a hump-back or -form relationship depending on the signs of the coefficients and . A lack of significance for indicates that the diversity-productivity relationship fits the monotonic model.

If the unimodal or form model defines the relationship between productivity and diversity, the species loss in high productive gradients suggests two possibilities: (a) mortality by inter-specific competition or (b) random, stochastic mortality. In order to test the first hypothesis, studies on the physiology and ecology of the present species are required, but this is not the focus of this research. The second hypothesis can be tested by fitting empirical models. To do so, it is necessary to establish the relationships between density, productivity, and diversity based on density. The first relationship is described by the following function (1): where is total density (number of stems per ha), is the parameter of productivity, measured as total aboveground NPP (Mg ha−1 year−1), and and are empirical constants.

The relationship between diversity and density is described by (2) where is diversity with stand density , is stand density (number of shrubs per ha), and and are empirical constants.

Equations (1) and (2) are combined in (3), which defines the increase or decrease in species diversity as a function of productivity as follows: where , , , and are empirical constants.

The size (basal diameter, top height) distributions of the species, the number of individuals, and the productivity of individual species or individuals of different sizes plotted in the altitudinal gradient allow us to draw some inferences about the role of competition, as well.

3. Results

The species richness of seedlings and sprouts after a year of growth following complete harvest was 37 species, distributed in 10 families and 31 orders (Table 1). The growth forms of these species are characterized as 70% shrubs, 18.9% small shrub, 5.4% herbaceous, and 2.7% trees and lianas (Table 1). In this early stage of succession, the legume family was the most abundant and most dominant. The mean number of species per 25 m2 across landforms was 8, with a standard deviation of 2 species. There was an observed tendency for increased species diversity from the plateaus to the floodplains, with means and standard deviations of in the plateau, in the slopes, ( in the valley, in the ephemeral floodplain, and in the intermittent floodplain. The total number of species in all three sites per landform also increased from the plateau (7) to the intermittent floodplain (18).

The total number of species increased nonlinearly as a function of productivity. However this relationship fitted an opposite nonlinear trend for the mature stages of succession. Total aboveground biomass and stand density increased from dry to riparian plots and so did mean basal diameter and mean top height but only for the mature succession stages (Table 2).

The correlation matrices for early (Table 3(a)) and late (Table 3(b)) successional stages between productivity, diversity, and environmental variables show S is better and statistically related to N in the early as well as to NPP, Text, and CN ratio for late successional stages. NPP is statistically related to Text and CN ratio for both early and late successional stages.

Species diversity was high in the early succession stage, in contrast to the species diversity present at the time of harvesting vegetation along the topographic gradient (Figure 1; see nomenclature of Table 1). After one year of postclearing growth, the relative portion of perennials (sprouts) and ephemerals (seedlings) can be attributed to the form of reproduction and to the different survival strategies of these life forms.

The species Amyris madrensis, Amyris texana, Havardia pallens, Bernardia myricaefolia, Condalia hookeri, Randia laetevirens, Sideroxylum langinosum, Forestiera angustifolia, Celtis pallida, Fraxinus greggi, Ehretia anacua, and Celtis laevigata are characteristic of the most productive sites (intermittent and ephemeral floodplains), whereas Lippia graveolens, Acacia rigidula, Acacia berlandieri, Calliandra conferta, Turnera diffusa, Krameria lanceolata, Helietta parvifolia, and Leucophyllum frutescens are typical shrub species found in low productivity and semi-arid sites.

In general, the most productive mature sites have a mean of only nine established shrub species. In the early stages of succession, mean species richness per plot increased to 25 species. This represents a species loss of 74% from the early stage to the late succession stage (9/25). In contrast, in the least productive landforms, the number of species present before harvesting consisted of a mean of 13 species, whereas in the early phases of succession, 20 species appeared during the course of one year. The rate of species loss in these sites was only 35% (13/20). These results suggest that the rate of local extinction is higher in the most productive sites, twice as high as that in the least productive sites. The differential growth rhythms and forms of germination provide a competitive advantage to tall species and species with asexual reproduction since these species provide strong competition for light and soil water in the early stages.

Seedling density in early succession is ten times greater than shrub density in the late stages of succession. Average seedling density one year after harvesting shrubs was 52,425 seedlings ha−1, whereas before harvesting, shrub density averaged 4,900 ha−1. These dramatic changes are due to the fact that shrubs grow over time, attaining larger dimensions and displacing or shading other shrubs by overcrowding. The 3/2 thinning law may operate at both intra- and inter-specific levels allowing competition to modulate both stand density and diversity. Intra- and inter-specific competition can cause mortality in individuals of the same species and individuals of different species with similar adaptive strategies. Mortality or species exclusion by factors associated with inter-specific competition must also be present, but discussing the mechanisms that cause this requires more frequent measurements between the different stages of succession. Additional data is also required in order to develop or to fit available models of succession [5].

The stocking reduction is less variable than species exclusion across the productivity gradient. The average survival rate (confidence interval) defined as (where = seedling density; = shrub density) is 11% (4%), and this is weakly related to stand productivity. That is, the mortality of individuals in the transition is high and somehow related, but with a small mean survival range across productivity levels (from 13% in high to 8% in low productivity sites, resp.). Shrub survival rate is higher in the least dense stands and lower in the densest stands. High productivity sites have smaller initial and final stand densities with higher survival rates. The inverse is true for plateaus and dry sites.

The soil variables display high dispersion for the different measures of productivity. At the plateau and slope sites, the texture is sandy-silty-clayey, with contrasting organic matter and nitrogen contents, while the mean organic matter content remains quite stable between productivity sites; the average soil nitrogen content is high (14%) in high and small (5%) in low productivity sites. The amount of organic matter and nitrogen is similar to that of the valley and intermittent floodplain and generally greater than those of the plateau; nevertheless, the soil textures are different, being clay-silty-sandy and clayey-silty, respectively. The ephemeral floodplain presents a clayey texture, with medium organic matter and organic carbon contents and the highest values of soil nitrogen (Table 4).

Aboveground NPP during the first year of regrowth correlates with total aboveground biomass in late succession across all landforms, with a tendency for increased aboveground biomass from the plateau toward the floodplains. This suggests that annual NPP can be estimated on the basis of late-successional biomass, using the equation: productivity = 0.85 + 0.034 total aboveground biomass; ; .

Productivity across plots varies from 0.62 to 5.35 Mg ha−1 year−1, with an average of 2.39 Mg ha−1 year−1 and a standard deviation of 1.32 Mg ha−1 year−1 (Figure 2). It tends to increase from the plateau toward the floodplains almost linearly with aboveground biomass stocks.

The soil nitrogen-productivity relationship fits better a power-law (; ) than a linear (; ) equation, although an -test provided insufficient evidence (; ) that the power-law model is statistically different from the linear model. Therefore, the productivity rate increases monotonically with soil nitrogen content (Figure 2). The classic triangular distribution has an upper boundary given by 3 Mg ha−1 y−1 and 0.15% of soil nitrogen content. Other limiting elements, such as soil carbon content, soil water content, and sunlight must control the lower aboveground productivity.

Aboveground biomass and soil organic carbon content fitted a negative linear relationship (TAB = CC; ). Aboveground biomass productivity is a log-normal function of the soil carbon/nitrogen (CN) ratio () (Figure 3). Low CN ratios are associated with high availability of mineral nitrogen, which is the form available to plants. In sites with reduced aboveground biomass productivity, soil nitrogen content seems to be the limiting factor. The CN ratio has a large effect on the microbial processes that regulate N availability in the soil. This effect is typical in sites of a low accumulation of organic matter as the carbon source, such as in the plateaus of the Coastal Plain of the northern Gulf of Mexico. A lack of other minerals appears to be the limiting factor at sites of high organic matter accumulation in plots established in the high productivity intermittent and ephemeral floodplains where decomposition rates are high in warm moist environments with low CN ratios of the plant material.

Shrub diversity shows a positive linear relationship to soil nitrogen content in the pioneer stages of succession (Figure 4). Unsurprisingly, the proportion of leguminous species is smaller at sites with high soil nitrogen content, like floodplains (5%) than at sites with low soil nitrogen content, such as the slopes and mesas (11%).

The relationship of soil nitrogen to species richness shifts from positive to negative between the early and late successional stages (Figure 4). That is, diversity diminishes monotonically with increasing soil nitrogen content. Soils of the harvested sites remained undisturbed (by, e.g., plowing, trampling, or grazing) during the one-year period of plant regrowth; therefore, soil organic content remained quite constant throughout the study period. This partially explains the increase in seedling diversity with increasing soil nitrogen in the pioneering stages, in contrast to the shrub diversity in the later stages of succession, during which shrub diversity decreased with increasing soil nitrogen content.

The diversity of shrubby species is also negatively related in a power fashion to the carbon/nitrogen ratio, although this relationship is not statistically significant (). Shrub diversity and calculated aboveground productivity are significantly related through monotonic models. That is, no humpback patterns fit the data of the late stages of succession in the Tamaulipan thornscrub plant community for measured as well as modeled data (Figure 5). Shrub diversity diminishes incrementally with increasing aboveground biomass productivity. This finding is consistent with the reduced shrub diversity in highly productive gradients demonstrated above.

The hypothesis of diversity reduction by the stochastic mortality of shrub species due to self-thinning was tested using the statistical relationships between density and productivity and diversity and density for the primary and late stages of succession (Models 2 and 3). These relationships were not statistically significant in the early stages; only the second relationship displayed statistical significance in mature scrublands. Therefore, the exclusion of shrub species by stochastic mortality explained less than 11% of the total S variation and therefore the hypothesis of self-thinning for the plant community was not supported by data from the Tamaulipan thornscrub community.

4. Discussion

Our average estimate of aboveground NPP in the early stages of succession 2.39 (±1.32) Mg ha−1 a−1 is higher than the mean current annual productivity for Tamaulipan thornscrub forests of northeastern Mexico reported by Návar et al. [26] of 1.61 Mg ha−1 a−1. Higher productivity values during the early phases of succession are consistent with the theory of highest growth rates during the initial stages of succession in most plant populations or communities [27, 28].

Neither the monotonic () nor the humpback () models fit the relationship between diversity and productivity in the early or mature stages of succession. The monotonic model shows significant changes in diversity with increasing productivity. The relationship between shrub diversity and aboveground biomass productivity differs in the early and late stages of succession. In the former, the relationship is positive, whereas in the latter stages, the relation is negative. Note that this shift is predicted by Smith and Huston [29] and discussed in detail by Huston [5]. Both theoretical models of plant competition and applied models of plant growth and competition predict that the pattern of diversity in relation to productivity reverses between early and late succession. This is a general result and should be found in nearly all vegetation types. This finding requires more information to test the hypothesis of unimodality; however, it does show that there is no single relationship that describes the diversity-productivity of the Tamaulipan thornscrub plant community at all successional stages. It is more likely that this relationship differs between the different succession stages of plant communities. These observations have been predicted earlier by several authors [5, 13, 29].

No species was found in all sites at both stages of succession. Only four and two species distributed in four sites and most species (17 and 6) were found in one single site for early and mature stages of succession, respectively. The number of species occupying different sites coalesces in 5 species distributing at three specific sites. The distribution of widespread species along this gradient fit a linear (; ) and a second-degree polynomial (; ) adjusted the relationship between the number of species () distributing along sites () for early and mature successional stages. The potential or physiological niche (e.g., [13]) and ecological niche converge at 2 species (D. texana and C. boissieri, in the early and H. pallens and A. rigidula in the mature stages of succession) distributed across four sites.

The model defining diversity versus productivity in the late stages of succession with stochastic species exclusion explained only 11% of the total variation in shrub species diversity, but showed a tendency toward reduced shrub diversity as aboveground biomass productivity increased. This information is consistent with the models by Grime [6] and Stevens and Carson [11, 14] that emphasize deterministic mortality caused by similitude in competitive ability of species to establish and last through the mature stages of succession. Therefore, it can be concluded that in the Tamaulipan thornscrub community, the rate of change of shrub species seems to only partially follow the mortality by stochastic processes that happen in monocultures fitting the 3/2 law of thinning. It is more likely that the mechanisms of intra- and inter-specific density dependence are simultaneously operating with a differential magnitude in different life stages and modified by disturbances (e.g., growth and droughts, among others). This seems to be due to the fact that inter-specific competition, although it begins in the early stages of succession, is still not as important as the processes of adaptation to environmental conditions that prevails on-site. Mortality by nonrandom processes also plays an important role and can explain the rest of the variation in the model of mortality by demography in the late stages of succession. As the species develop in dimensions and density is extremely high, both intra- and inter-specific competition reduce the number of individuals and seem to reduce the number of species as the succession advances as well. The inter-specific competition gradient moves from the plateau to the riparian community, with quasi null competition in dry and strong inter-specific competition in riparian places.

The range from a low species survival rate of 36% in highly productive environments to a high of 65% in low productive environments suggests that inter-specific competition is stronger in more productive landscapes and less pronounced in less productive sites in the Tamaulipan thornscrub forests. In fact, the relationship between species survival and aboveground biomass productivity was negatively and highly significant (; ; ; ) where SR is rate of species survival (number of species in late/number of species in early stages of succession) and P is productivity (Mg ha−1 y−1). The rate of species exclusion is . Although a similar relationship was observed between rate of species survival and soil nitrogen content, it was not statistically significant. That is, all abiotic factors that control productivity, but in particular soil nitrogen content, should be considered stressors, given by the deficiency of resources. These findings are partially consistent with theoretical models that have predicted the relative importance of facilitation, and competition will vary inversely across gradients of abiotic stress, with facilitation being the dominant interaction under high abiotic stress conditions [30, 31].

Stand density is reduced more homogeneously than species diversity along the productivity gradient, with values ranging between 6 and 16% for the pioneer to the late succession stages. However, stand density is not related to stand aboveground biomass productivity, but only to aboveground biomass, which is the long term sum of the productivity. Since species exclusion is less pronounced in low productivity sites, intra-specific competition is more acute in high abiotic stress gradients since stand density survival reduces in a negative power fashion with productivity. In other words, intra-specific density dependence plays an important role in higher abiotic stress gradients or it could be simply abiotic stress with no intra-specific interactions. It appears that a combination of both intra- and inter-specific density dependence is more important in more productive areas. The former may be explained by the temporal and spatial overlap in the use of limited resources of individuals of the same species.

In late succession stages, shrub diversity is higher in less productive sites and so must be the variety of survival mechanisms they use during development. This can be explained by the variety of strategies they use in the temporal and spatial exploitation of limited resources (soil water and fertility). Stable coexistence in harsh and variable environments requires mechanisms that distinguish species ecologically, which may be a contrast to Huston’s [13] findings (such differences are more important in stable productive environments and less important in harsh and variable environments). This is a prerequisite for intra-specific density dependence to be stronger than inter-specific density dependence. In high productivity landscapes or environments with lower abiotic stress this is an issue primarily in equilibrium competition analysis such as Lotka-Volterra models or Tilmans models [13]. Under nonequilibrium conditions, disturbances and environmental variability interrupt the process of competitive exclusion, such as disturbances on the intermittent and ephemeral floodplains. Intra- and inter-specific competition are better balanced under average site conditions since species exclusion reaches 36% and self-thinning mechanisms achieve above 16% of the original stand density in early successional stages.

Finally, both negative interactions (e.g., competition, [6, 7, 13, 32]) and positive interactions (e.g., facilitation, [13, 30, 33, 34]) are likely to be present in this plant community. These occur simultaneously, but their dominance depends on the successional stage of development of the plant community. The relative dominance of these interactions must control the relationship between species diversity and productivity. This is a matter for further research.

5. Conclusions

Tree diversity changes in both space and time. In pioneering succession stages it declines from the high to the low productivity gradients but it reverses in the late phases of succession. The rate of species survival increases from 36% in the transition from the pioneering to the late stages in highly productive gradients to 65% in low productivity environments. Stem density is also greater in the pioneering rather than the late stages of succession. The mortality of individuals and species by a variety and magnitude of processes is important in this plant community. After this short period of study, the colonization and permanence of species are attributable to several factors, within which intra- and inter-specific competition and facilitation are factors that modulate the exclusion and coexistence of individuals and species. This quantitative evidence suggests that inter-specific competition is more pronounced in more productive sites than in less productive sites, although fewer shrub species colonize the latter due to low productivity levels in the early stages of succession. Stochastic mortality explained less than 11% of the species exclusion. It seems that the deterministic mechanisms related to competitive exclusion where there are similarities in the use of resources, physiology, and exploitation of soil and light, among others, can overlap between several species in productive sites and that this similarity vanishes partially over time in less productive sites. It is also likely that abiotic stress may be related to species mortality. Further research is required in order to further understand the mechanisms of species exclusion and coexistence.

Conflict of Interests

The authors have no conflict of interests to declare.

Acknowledgments

The authors of this report thank Raul Bretado and Maricela Peña for their help with data collection. We also thank the numerous reviewers who improved the scientific understanding and readability of this manuscript. The Instituto Politécnico Nacional authors are COFAA-IPN scholars.