Skip to main content

Structure, composition and carbon stocks of woody plant community in assisted and unassisted ecological succession in a Tamaulipan thornscrub, Mexico



Ecological restoration is a process that helps the recovery of an ecosystem that has been degraded, damaged or destroyed. The success of restoration actions depend on the understanding of the processes, mechanisms and factors that guide vegetation dynamics. The restoration of plant communities can be made by unassisted (passive restoration, where the community recover by itself) and assisted (active restoration) ecological succession. It is imperative to know the scope of both types of activities to guide future restoration actions by evaluating the properties and functioning of the intervened communities.


Variance analysis of the carbon (C) stocks, basal area, canopy area, Shannon–Weaver index values, specific richness and abundance of three Tamaulipan thornscrub communities (assisted natural succession area, unassisted natural succession area and control area) was performed. Furthermore, a similarity analysis between the sampling areas using the floristic composition (abundance) was performed.


In total, 11 families, 17 genera and 20 species of vascular plants were registered. The richness of species and abundance ranked as follows per area: assisted ecological natural succession > control > unassisted ecological natural succession. The species composition between sampling areas showed a low number of common species between plant communities.


The values of species richness, diversity, abundance, basal area, and canopy area of the assisted natural succession vegetal community was statistical similar to the control plan community. The values of C stocks showed that assisted ecological succession could recover not only structure and composition attributes but also this key ecosystem property.


Ecological restoration is a process that helps the recovery of an ecosystem that has been degraded, damaged or destroyed [61]. Aims to recover the biodiversity and ecosystem services, considering the ecological, socio-economic and cultural processes [56]. Recently, the recovery of disturbed landscapes and biological communities is becoming a priority in many regions of the world [5]. However, the results of restoration are different and depends largely on the conditions and type of restoration action used to achieve the desired success.

In general terms, ecosystem ecological restoration can be classified into two types according to the strategy used: passive restoration where the ecosystem recovery by itself, like unassisted natural succession, and active restoration, where restoration activities are undertaken to accelerate the recovery of the ecosystem, like planting species to accelerate recovery in an assisted natural succession [12, 42, 61]. The unassisted ecological succession (passive restoration) can achieve good results in certain cases, but usually takes decades or centuries and generally does not guarantee a composition of species like the original state [1]. On the contrary, the assisted ecological succession (active restoration), usually shows better results, allowing the establishment of a higher canopy that improves the ecological conditions of the ecosystems in less than ten years [56, 66]. Regardless of the restoration chosen, it is imperative to carry out studies and actions that help ecological recovery in regions whose biological and ecosystem biodiversity is strongly threatened [32, 57, 58].

As mentioned on Linding-cisneros research on the challenges of a multicultural and megadiverse country like Mexico [30, 64], the government has promoted the restoration in natural protected areas and degraded communities. In this aspect, Tobon et al. developed a conceptual framework on the restoration of all terrestrial ecosystems, including xeric shrublands [64]. These plant communities occupy large areas in northern Mexico and in many cases the original vegetation has been eliminated to establish agricultural or livestock areas [6]. For this reason, restoration actions on shrublands such as the Tamaulipan thornscrub are key activities to recover and conserve vegetation structure and diversity [4, 50, 65].

The Tamaulipan thornscrub is a plant community usually found on deep clay soils, plateaus, hillocks and plains between 200 and 500 m altitude, under the influence of warm and semi-warm climates [60]. It is restricted to northeastern Mexico and southern Texas, and its original distribution in these areas has been reduced by 70% and 90%, respectively, due to changes in land use for agricultural and livestock purposes [7]. This scrub has a structure formed mainly by two strata, one shrubby (2–5 m high) usually composed of thorny elements and another herbaceous that does not exceed the height of a meter [60]. It has a high species richness (S = 21; [33, 37], S = 29 [24]) including species such as: Celtis pallida Torr. (Cannabaceae), Opuntia leptocaulis DC. (Cactaceae), Prosopis glandulosa Torr. (Fabaceae) Acacia rigidula Benth. (Fabaceae) and Ruellia nudiflora (Engelm. & A. Gray) Urb. (Acanthaceae), being some of these, important forestry and silvopastoral species [19].

Although it has been observed that the ecological restoration increases biodiversity and ecosystem services in general, results can be different for the different types of restoration and ecosystems [10, 26, 56]. Thus, studies that systematically evaluate the results of active and passive ecological restoration in the same location after the same disturbance, related to the recovery of structure, diversity and provision of ecosystem services such as carbon sequestration, are important to be undertaked in differents ecosystem [26, 39, 53, 56]. Furthermore, although the Tamaulipan thornscrub has been intensively studied in recent years [24, 33, 53, 59], only a few investigations have evaluated the structure and composition of woody plan community after agricultural activity in this ecosystem. Until now, research efforts on the Tamaulipan thornscrub have focused on timber volume production [20, 47], species survival [4, 65] or leaf and litter biomass [48].

The objective of this study is to evaluate the structure and composition such as species richness, diversity, abundance, basal area, canopy area and carbon (C) stocks of woody plant communities of the Tamaulipan thornscrub and to compare the differences in these ecosystem attributes between assisted/unassisted ecological succession and control areas. Our hypothesis is that the assisted ecological succession and control areas will show similar vegetation structure and composition, as well as in terms of the ecosystem service carbon sequestration.


Study area

The study area is located in the municipality of Pesquería, Nuevo Leon (Northeast Mexico, Fig. 1), 25°45´17.78´´ north latitude and 99° 58´01.40´´ west longitude. The climate is very dry to semi-warm (BWhw) according to the classification of Köppen (modified by Enriqueta Garcia) [23], the average annual temperature is between 20 and 21 °C. Soil types on the site are: xerosol, castañozem, feozem, regosol and minimally, fluvisol, vertisol, and rendzina [4], The average annual rainfall on the study area is 550 mm and an altitude of 330 m.a.s.l. The representative plant community consists mainly of Tamulipan thornscrub with shrubs as secondary vegetation mainly dominated by Acacia farnesiana [3].

Fig. 1
figure 1

Study area location

Vegetation analysis

In November of 2017, the structure and composition of woody plant communities were investigated through a floristic composition and diversity evaluation on three areas: a control area, an assisted ecological succession area and an unassisted ecological succession area. The species registered were identified at the species level by qualified personnel of the Faculty of Forest Science of the Universidad Autonoma de Nuevo León (UANL). Their scientific names were validated against names in Tropicos ( and IPNI ( databases.

The control area is composed of a Tamaulipan thornscrub climax community where there are not records of productive activity involving total or partial removal of vegetation (such as: livestock, slaughter, agriculture etc.). The assisted ecological succession areas have a history of agricultural use and received an assisted ecological succession treatment after abandonment in September of 2009. The unassisted ecological succession area has the same history of use and abandonment as the previous mentioned site but assisted ecological succession actions were not performed. The agricultural activity mentioned on both assisted and unassisted areas consisted of plantations of oats, wheat and corn for approximately two decades (1989–2009). The crops were sowed in a rotating manner and irrigated by a gravity-driven water channel. Planting and harvesting activities were made using a farm tractor. All areas had similar characteristics: vertisol type of soil, 310 m. in elevation, the slope of < 3%.

The assisted ecological succession consisted of reforestation with native species in April of 2010. The species used for reforestation were rescued from a land-use changed area, authorized by SEMARNAT (Secretariat of Environment and Natural Resources) due the absence of nurseries producing the high diversity of native species of Tamaulipan thornscrub. The rescue started with the selection of abundant individuals of smaller size: arboreal and shrub species with heights from 0.5 to 2.0 m and basal diameters from 0.5–1.5 cm were considered, and each individual was labeled with: the common name, serial number and coordinates in where it was found. The specimen extraction process started with the pot preparation, using plastic paper to cover the roots to prevent from exposure and soil crumbling, then, the individuals were moved to a different site by trailer. When the plantation actions started, the specimens were placed in the strains and the plastic paper removed, fungicide and root promoter were applied, and plant tutor sticks were placed followed by the first irrigation. The plantation area was periodically cleaned of weeds. The plantation was made using a triangular plantation system with a distance of 1.5 m between individuals and lines having a plantation density of 992 n ha−1. The process scheme can be observed in Additional file 1: Appendix.

An analysis of preliminary samples obtained for the study area was carried out seven years before (September 2017) to obtain a value of sampling units needed to register 95% of the richness of each community. For doing so, we used the equation of Mostacedo and Fredericksen [43]:


where n is the adequate number of sampling areas, E is the error (20%), t is the extracted value of t student tables (p < 0.05), N is the total sampling unites and CV the variation coefficient. According to this mathematical model, 12 sampling plots of 10 × 10 m (100 m2) were randomly established in each area. A census of all plant species higher than 10 cm was carried out. Dendrometric measurements of total height (h), basal diameter (d0.10), which was estimated from the basal areas of all the ramifications, (cross sectional area) and canopy diameter (dcanopy) were carried out on each individual. Diametric classes (5, 10, 15, 20, 25, 30, 35 and > 35) were established.

Data analysis

We analyzed the relative abundance, dominance, and frequency of plant species to estimate their importance value index [43, 44]. The absolute abundance was calculated using the number of individuals found on the plots and was extrapolated to hectare, using the following formula: Ai = Ni/S, where Ai is the absolute abundance of species i, Ni is the number of individuals of species i, and S the sample surface (in hectares). The relative abundance of species was calculated using the equation: Ari = (AiAi)*100, where Ari is the relative abundance and ΣAi is the total abundance of species i. To estimate dominance the basal area for each individual was calculated. The basal area was calculated with the formula A = π/4*d2, where A is the basal área, π is a constant value of 3.14, and d is the basal diameter (d0.10). The absolute dominance is Di = Abi/S, where Di is the absolute dominance of species i, Abi is the basal area (in square meters) of species i, and S the sample surface (in hectares). The relative dominance was calculated using the equations Dri = (DiDi)*100, where Dri is the relative dominance of species i total dominance over and ΣDi is the total dominance of species. The relative frequency is Fri = (FiFi)100, where Fri is the relative frequency of species i over the total frequency, Fi is the number of plots that present the species i and ΣFi is the total frequency of all species. The importance value index (IVI) was calculated as follows: IVI = (Σ Ari + Dri + Fri)/3, where Ari is relative abundance, Dri is relative dominance and Fri is relative frequency [43, 44].

To determine the alpha diversity of the studied communities, two indices were used: Margalef (DMg), which is based on the quantification of the number of species present (specific wealth), and Shannon (H´), which refers to the numerical structure of the community or proportional distribution of the abundance of each species [40]. The Margalef index was calculated with the following equation DMg = (S-1) / ln (N), where S is the number of species present and N the total number of individuals.

The similarity assessment between the sampling areas was performed using the floristic composition (abundance). For this purpose, a hierarchic cluster analysis was carried out using the block-distance measure and Ward’s clustering method, graphing it on a similarity dendrogram. Statistical analyzes were performed using the software R ver. 2.6.0. with cluster package [54]. In addition, the communities (control, assisted ecological succession and unassisted ecological succession area) were explored by Bray–Curtis ordination analysis [9]. The analysis is based on an algorithm allowing analyzing the similarity of the communities sampled by calculating the percentage of similarity (0–100%) in a distance measure [8]. The results were presented in a dendrogram and the statistical analyzes were performed using the software package BioDiversity Professional Ver. 2 [38].

Carbon stocks

To determine the carbon stocks of the arboreal and shrub species of Tamaulipan thornscrub, first, the aerial biomass was determined using the local allometric equation (r = 0.80). [45].


where, BT = total aerial biomass, d = basal diameter and h = total height.

The arboreal and shrub species in the Tamaulipas thornscrub do not have allometric equations that assessment of carbon stock, therefore, the carbon stock was calculated using the 0.454 factor according to [68], where it is interpreted that 45.4% of the biomass aerial corresponded to carbon stocks.

Statistical tests using IBM SPSS Statistic ver. 19 © [69] were carried out to verify that the assumptions of the residuals were fulfilled, normality using the Shapiro–Wilk test and homogeneity of variances by the Levene test, both with a significance value (p < 0.05). Since they fulfilled the assumptions, a generalized linear model (glm) was used with a Gaussian distribution function for Shannon–Weaver index values, specific richness, abundance and an analysis of variance (ANOVA) of one factor (area) was performed for carbon stock, basal area and canopy area. The Duncan test was used to determine if there were significant differences (p < 0.05) of the variables analyzed per area.


In total, 11 families, 17 genera and 20 species of vascular plants were registered. Fabaceae was the family with more species with 9, followed by Boraginaceae with two. The Additional file 2: Appendix II shows the floristic inventory of the study areas in detail.

The control area had 13 species, where, Cordia Boissieri A. DC., Parkinsonia texana (A.Gray) S. Watsonhad and Acacia rigidula were dominant. The assisted ecological succession area had 12 species with species like Acacia farnesiana, Leucophyllum frutescens (Berland.) I.M. Johnst, Cordia boissieri. The unassisted ecological succession area registered six species; Acacia farnesiana, Leucophyllum frutescens, and Prosopis glandulosa were the most representative species in this area.

The richness of species using the Margalef index showed significant differences between the studied areas. The assisted ecological succession area showed greater richness with 5.58 ± 0.43 (average ± standard deviation), followed by the control area (4.41 ± 0.37) and unassisted ecological succession area (0.91 ± 0.25). The values of the Shannon index in the assisted ecological succession areas (0.64 ± 0.04) and control area (0.54 ± 0.04) showed no significant differences, and both were higher than the unassisted ecological succession (0.30 ± 0.01) (Fig. 2).

Fig. 2
figure 2

A) Margalef index (DMg), Shannon index (H ') Abundance (m2 ha−1), Basal area (m2 ha−1), Plant cover (m2 ha−1) and carbon stock (Mg ha−1) of the evaluated communities. Means (± standard error, n = 12) followed by different letters (a, b, c) show significant differences at P = 0.05 (Duncan test)

The Table 1 showed the results of the GLM, χ2-test and the freedom degrees for each sampling area. The Shannon Index (81.429), Species richness (Margalef index) (68.411) and Abundance (155.394) the values of χ2 were less than χ2-table (27.59) at alpha (α = 0.05), therefore, the null hypothesis (equality of means) was rejected. This means that the observed data fit the theoretical distribution. The results of ANOVA for basal area and canopy area were showed in the Table 1 and the F-statistic value were 15.839 and 33.407, respectively; therefore, rejected H0 (H0: μ1 = μ2 = μ3) because these values were > 3.32 (F-table). Thus, there were difference statistically significant at α = 0.05 between sampling areas (Control area, assisted ecological succession and unassisted ecological succession area). In the post-hoc test and their p-values, the mean values of the Control area and assisted ecological succession do not have significantly different, but the mean values of the assisted ecological succession and unassisted ecological succession area have significantly different for Shannon Index, Species richness (Margalef index), Abundance and basal area (Fig. 2A, B, C and D). Only canopy area had significantly different between three sampling areas (Fig. 2F). Regarding the analysis of diametric classes of each area, the control area was the only one showing individuals equal to or greater than 35 cm and registered the biggest number of individuals in the 10 cm category (more than 50% of the individuals), while the assisted ecological succession showed more than 50% in the 5 cm category. In the three areas, the greatest number of individuals was concentrated in the categories of 5, 10 and 15 cm respectively (Table 2).

Table 1 Generalized linear model (glm) for Shannon–Weaver index values, specific richness, abundance, and analysis of variance (ANOVA) of one factor (area) was performed for carbon stock, basal area and canopy area
Table 2 Number of individuals per hectare in each diametric category

Absolute density, dominance and importance value index of all the studied areas are shown in Table 3. The abundance values of the control area (1175 ± 121 N ha−1) and assisted ecological succession (1416 ± 178 N ha−1) areas were higher than those of the unassisted ecological succession area (100 ± 27). The basal area of the vegetation communities of the assisted ecological succession and control area showed statistical similarity (7.86 ± 1.30 and 8.31 ± 1.03 m2 ha−1 respectively) and were significantly higher than the unassisted ecological succession community (0.85 ± 0.74 m2 ha−1).

Table 3 Absolute density (number per hectare), absolute dominance (square meters per hectare), and importance value index (IVI, percent) of the three areas evaluated

The canopy area of the individuals in the assisted ecological succession and control area showed statistical similarity. The control area showed overlap in the canopy area with 13,365 ± 1470 m2 ha−1 and the assisted ecological succession with 9298 ± 880 m2 ha−1, which is equivalent to 92% coverage. The unassisted ecological succession community presented 1039 ± 777 m2 ha−1, equivalent to the 10% coverage.

The composition of species per sampling area is showed in Fig. 3. After clustering control area (plots 1–12) five sites were grouped in unassisted ecological succession area (plots 1, 6–8 and 11) and one in the assisted ecological succession area (4). While the unassisted ecological succession area (plots 25–36) showed similarity with the other areas, adding plots 1, 6–8, 11, 14–18, 22–23. The assisted ecological succession area (13 to 24) showed a lower similarity in species composition, only having four of the original sampling plots (13, 19, 20,21), and adding another one (4). This means that the sampling areas contain de majority of the sites with their original category (Control area, assisted ecological succession and unassisted ecological succession area), for this reason the analysis of similarity between communities was carried out. The species composition analysis showed few common species between plant communities. From the numerical classification, a group formed by control and assisted ecological succession areas was distinguished with 23% of similarity, while the similarity with unassisted ecological succession area is 16% (Fig. 4).

Fig. 3
figure 3

Dendrogram of similarity by hierarchic cluster analysis using the block-distance measure, and Ward’s clustering method between the different evaluated areas based on abundance, on the analysis. Red box: Control area, green box: Unassisted ecological succession area and blue box: Assisted ecological succession area

Fig. 4
figure 4

Dendrogram of similarity between the different evaluated communities based on the Bray–Curtis analysis

Carbon stocks

The results of ANOVA for carbon stocks are shown in the Table 1. The F-statistic value was 14.573; therefore, we rejected H0 (H0: μ1 = μ2 = μ3) because 14.573 > 3.32 (F-table). Thus, there was a statistically difference significant at α = 0.05 between sampling areas. In the post-hoc test, the mean values of the control area (10.35 ± 1.75 Mg ha−1) and assisted ecological succession (7.45 ± 0.92 Mg ha−1) do not showed significant statistical differences, but the mean values of the assisted ecological succession and unassisted ecological succession area (0.99 ± 0.88 Mg ha-1) showed significant statistical differences for Carbon stock (Fig. 2G). More than 50% of the concentration of carbon stock registered in both control and assisted ecological succession area was centered in diameter categories from 5 to 20 cm, in the unassisted ecological succession area the higher concentration of carbon was found in the 30 cm category.


Fabaceae showed the largest number of species in the study area. The representativeness of the Fabaceae family was also reported in previous studies carried out in North Mexico, where the family was identified as a big diversity group [15] in conserved areas, in secondary vegetation communities and assisted ecological succession communities [2, 25, 35, 39].

The diversity of Fabaceae is attributed to the ecological plasticity of the group, which contains species with different seed sizes, types of dormancy, germination rates [17, 27,28,29], longevities [11, 14], growth rates and patterns [14, 21] and tolerance to water scarcity and to nitrogen [3, 13, 24, 50]. All above mentioned factors allow members of the Fabaceae family to be present in different environmental conditions [16, 46, 55] influencing the existence of different compositions within the same Tamaulipan thornscrub [15, 55]. Factors such as propagules availability and ecological performance, along with the interactions between other taxonomic groups presented in the Tamaulipan thornscrub community are responsible of determining vegetation performance [52].

The comparison of richness and diversity of species between sites were similar to the values reported on previous studies [2, 25, 35]. The unassisted ecological succession area low richness and diversity of species can be largely attributable to the dominance of one or two taxa over others, which can be due the action of environmental filters limiting the presence of other species (e.g., low nitrogen content, insolation and drought) as suggested by [31].

The importance value index of the control area is similar to the values reported by [39] showing the characteristics of a mature community [18]. The presence of Cordia boissieri and Parkinsonia texana in a mature secondary vegetation community like in the control area has been reported by [15, 45]. Those species are considered as pioneer species of early successional stages [11, 14], showing a life spectrum of 30–50 years [11, 14] which could explain why they exist in the control area. Besides, the rapid growth of P. texana, along with the high germination rate in share light conditions with C. boissieri [29], and the establishment capacity under shadow conditions suggested by its seed size. [30] could explain the Importance Value registered. Acacia rigidula has a higher germination rate than C. boissieri [27, 30], but the size of its seed suggests low tolerance to shade and rapid growth [22, 34], which could explain why it occupies the third important value index of the control area.

Regarding the secondary natural succession trajectory, the presence of Acacia farnesiana and Prosopis glandulosa in the vegetation community of the unassisted ecological succession area match the results of [15, 25] registering the species in areas that have suffered some type of disturbance due to anthropogenic activity. A. farnesiana is a pioneer species that presents a high germination rate in temperature ranges from 12 to 28 0C [17, 29, 30] and like P. glandulosa, it is characterized by their rapid growth [16, 21] and its tolerance to edaphic drought which is a common characteristic of the study area. This species also has higher values in hydric potential of the xylem on the pre-dawn and in midday, under conditions of lack of water, compared to other species [13, 24]. Its presence is also associated with low availability of nitrogen in the soil, which results in a limiting factor for other groups and an advantage for the species of this family. These aspects suggest that unfavorable habitat conditions advantage the establishment and dominance [3, 50], making it an important group in the initial stages of succession and suitable for restoration activities.

The plantation abundance of the assisted ecological succession was 992 N ha−1, and after seven years it was 1416 ± 178 N ha−1, which indicates that in addition to the established individuals there was a regeneration of native species. The composition of this assisted community was partly a product of the species selected for revegetation, the restoration activities, as well as of the ecological performance of the species observed at the end of the study.

The highest importance value found in A. farnesiana can be explained through some factors such as: the advantage of its ecological performance in harsh environments, the low mortality of seedlings and saplings in plantation conditions [21] and the high capacity for canopy development in plantation conditions (which also explains the higher basal area). The rapid development of the canopy can strongly limit the germination, establishment, and development of light-demanding species [22], which can, on the other hand, benefit the species tolerance to shading conditions like C. boissieri, that is also included in the actions of restoration, and that occupies the third value of importance in this area.

Leucophyllum frutescens, a pioneer species with the second most important value in the assisted ecological succession area is commonly associated with low productivity areas, which could explain the ability to withstand limiting condition areas [47]. P. glandulosa was the species found in all studied areas, alongside the tolerance to adverse environmental and edaphic conditions, show low mortality in plantations with seedlings surviving rate above 85% [21]. However, as growing time passes, the plantations are susceptible to pest attacks, which can increase the mortality, this could explain the low importance value in the assisted ecological succession area (IVI = 7%).

Ebenopsis ebano (Berland.) Barneby & J.W Grimes and Acacia rigidula species of moderate growth [14, 28] presented on both control and the assisted ecological succession area. These species can withstand the stress conditions [21], which also makes them an adequate option for restoration actions.

Similar to what was reported in this study, other authors [36, 38] showed that the assisted ecological succession treatment was effective in terms of recovering the richness and diversity of tree and shrub species compared to passive restoration actions in northeastern Mexico, which explains to a large extent the similarity between the control and the assisted ecological succession site, unlike unassisted ecological succession area.

The dissimilarity in the composition of plant communities observed in our study is similar to the results recorded by Alanis et al. [3] on vegetation communities of the Tamaulipan thornscrub with a different productive use and obtained inferior values of similarity. The differences in species composition could be explained as differences in growing of the communities in control and unassisted ecological succession area, also the influence of the origin of the structure and composition of the control area compared to the assisted ecological succession area, which is strongly limited to the species which it was intervened. Additionally, the dissimilarity of the communities in our study and the study of Alanis et al. [3] can be explained through the heterogeneity of Tamaulipan thornscrub vegetation that usually shows different species composition according to the land use, topographic and edaphic conditions y differences in altitude [55].

The differences between assisted ecological succession and control area plant communities are determined by the dissimilarity of common species and high-density values. In the control area, the species are Cordia boissieri (243 N ha−1) and Parkinsonia texana (201 N ha−1) while in the assisted ecological succession area Acacia farnesiana (374 N ha−1) and Leucophyllum frutescens (263 N ha−1).

Carbon stocks

It is relevant to analyze the carbon stock in the three areas, since CO2 is the main greenhouse gas, and that plant species play a fundamental role in the absorption and fixation of carbon [67]. In this regard, the highest value of captured carbon registered for the control area and assisted ecological succession communities compared to the capture of unassisted ecological succession area is mainly attributable to the higher number of individuals and the diametric categories observed in these areas (Table 2). The unassisted ecological succession area concentrates practically all the individuals of the smaller diameter category, which is logical because to reach similar values as the control area must take several years for the development of the vegetation and the recovery of the ecological and dasometric normal parameters.

Compared to the dasometric values of the control area (mature thornscrub community), the ecological, dasometrical and C Stock variables of the assisted ecological succession shrubland area were recovered in seven years after the restoration the activities. Similar distribution pattern of the individuals in the diametric classes, and thus of C stocks, were observed in other successional studies of the Tamaulipan thornscrub communities, where early successional stages concentrated the individuals of the smaller diameter category and had a low C stock [36, 67]. This pattern suggests that vegetation from assisted ecological succession shrubland follows an optimal growth trend [67].


Our results provide valuable information about the restoration and regeneration condition and structure and composition of woody plant community after seven years of assisted ecological succession treatment in a Tamaulipan thornscrub located in northeastern Mexico.

The values of species richness, diversity, abundance, basal area and canopy area of the assisted ecological succession vegetal community presented statistical equality with the control plant community.

The values of captured carbon showed that assisted ecological succession could recover not only structure and composition attributes but also this ecological service. These results indicate that the assisted ecological succession treatment was effective on recovering the richness and diversity of tree and shrub species compared to passive restoration actions in the Tamaulipan thornscrub studied.

Availability of data and materials


  1. Aide TM, Ruiz-Jaen MC, Grau HR. What is the state of tropical montane cloud forest restoration? In: Bruijnzeel LA, Scatena FN, Hamilton LS, editors. Tropical montane cloud forests science for conservation and management. Cambridge: Cambridge University Press; 2010. p. 101–10.

    Google Scholar 

  2. Alanís E, Jiménez J, Aguirre O, Treviño E, Jurado E, González M. Efecto del uso del suelo en la fitodiversidad del matorral espinoso tamaulipeco. Ciencia UANL. 2008;11(1):56–62.

    Google Scholar 

  3. Alanís E, Jiménez J, González MA, Yerana JI, Cuellar LG, Mora-Olivo A. Análisis de la vegetación secundaria del matorral espinoso tamaulipeco, México. Phyton Int J Exper Bot. 2013;82:185–91.

    Google Scholar 

  4. Alanís E, Molina VM, Rechy L, Alcalá AG, Marín JD, Pequeño MA. Composición, diversidad y sobrevivencia de un área restaurada en el Complejo Siderúrgico de Ternium, Pesquería, México. In: Experiencias mexicanas en la restauración de los ecosistemas.  Primera edición. Cuernavaca, UNAM. 2016;255–274. ISBN: 9786078434763.

  5. Aronson J, Alexander S. Ecosystem restoration is now a global priority: time to roll up our sleeves. Restor Ecol. 2013;21(3):293–6.

    Article  Google Scholar 

  6. Arriaga L. Implicaciones del cambio de uso de suelo en la biodiversidad de los matorrales xerófilos: un enfoque multiescalar. Investigación Ambiental. 2009;1(1):6–16.

    Google Scholar 

  7. Arriaga-Cabrera L, Espinoza-Rodríguez JM, Aguilar-Zúñiga C, Martínez-Romero E, Gómez-Mendoza L, Loa-Loza E. Regiones terrestres prioritarias de México. Comisión Nacional para el conocimiento y uso de la biodiversidad. Mexico. 2000. ISBN: 970–9000–16–0

  8. Beals EW. Bray-Curtis ordination: an effective strategy for analysis of multivariate ecological data. Adv Ecol Res. 1984;14:1–56.

    Article  Google Scholar 

  9. Bray JR, Curtis JT. An ordination of the upland forest communities of Southern Wisconsin. Ecol Monogr. 1957;27:325–49.

    Article  Google Scholar 

  10. Bullock JM, Aronson J, Newton AC, Pywell RF, Rey-Benayas JM. Restoration of ecosystem services and biodiversity: conflicts and opportunities. Trends Ecol Evol. 2011;26(10):541–9.

    Article  PubMed  Google Scholar 

  11. Burch DG, Will AA. Trees that know their place. Proc Florida State Horticult Soc. 1990;193:366–8.

    Google Scholar 

  12. Chazdon RL, Guariguata MR. Natural regeneration as a tool for large-scale forest restoration in the tropics: prospects and challenges. Biotropica. 2016;48:716–30.

    Article  Google Scholar 

  13. Domínguez GT, González H, Ramírez R, Estrada E, Cantú I, Gómez M, Villarreal J, Alvarado S, lvarado, Alanís G. Diversidad estructural del matorral espinoso tamaulipeco durante las épocas seca y húmeda. Revista Chapingo , Serie Ciencias Forestales y del Ambiente. 2013;4(17):106–22.

    Google Scholar 

  14. Elizondo R. Guía de árboles y otras plantas nativas en la zona metropolitana de Monterrey. Monterrey: Fondo Editorial de Nuevo León; 2009. p. 320.

    Google Scholar 

  15. Estrada E, Yen AD, Villarreal J. Leguminosas del centro del estado de Nuevo León, México. Anales del Instituto de Biología, Universidad Nacional Autónoma de México. Serie Botanica. 2004;75:73–85.

    Google Scholar 

  16. Ewing K, Best C. South Texas tamaulipan thornscrub restoration experiment measures growth of planted woody vegetation. Ecol Restor. 2004;22(1):11–7.

    Article  Google Scholar 

  17. Flores J, Jurado E. Germination and early growth traits of 14 plant species native to Northern Mexico. Southwesthern Natur . 1998;43(1):40–6.

    Google Scholar 

  18. Foroughbakhch R, Maginot NH. Structure and Floristic Composisions of Ecosystems of Northeast of Mexico. Chron Bioresour Manage. 2017;1(2):060–4.

    Google Scholar 

  19. Foroughbakhch R, Alvarado-Vázquez MA, Carrillo Parra A, Hernández-Piñero JL, Guzmán Lucio MA. Floristic diversity of a shrubland in northeastern Mexico. Phyton. 2013;82:175–84.

    Google Scholar 

  20. Foroughbakhch R, Carrillo Parra A, Hernández Piñero JL, Alvarado Vázquez MA, Rocha Estrada A, and Cardenas ML. Wood volume production and use of 10 woody species in semiarid zones of Northeastern Mexico. Int J Forestry Res 2012;1–7. doi:

  21. Foroughbakhch R, Hernández-Piñero JL, Carrillo-Parra A. Adaptability, growth and firewood volume yield of multipurpose tree species in semiarid regions of Northeastern Mexico. Int J Agric Policy Res. 2014;2(12):444–53.

    Article  Google Scholar 

  22. Foroughbakhch R. Establishment and growth potential of fuelwood species in northeastern Mexico. Agrofor Syst. 1992;19:95–108.

    Article  Google Scholar 

  23. García E. Modificaciones al sistema de clasificación climática de Köppen. México: Instituto de Geografía, UNAM; 1964.

    Google Scholar 

  24. González-Rodríguez H, Cantú I, Ramírez RG, Gómez MV, Pando M, López JM. Potencial hídrico xilemático en cuatro especies arbustivas nativas del noreste de México. Revista Chapingo Serie Ciencias Forestales y del Ambiente. 2011;17:97–109.

    Article  Google Scholar 

  25. Jiménez J, Alanís E, González MA, Aguirre OA, Treviño EJ. Characterizing regeneration of woody species in areas with different land-history tenure in the Tamaulipan thornscrub Mexico. Southwestern Natur. 2013;58(3):299–304.

    Article  Google Scholar 

  26. Jones HP, Jones PC, Barbier EB, Blackburn RC, Rey Benayas JM, Holl KD, Mateos DM. Restoration and repair of Earth’s damaged ecosystems. Proc Royal Soc B: Biol Sci. 2018;285(1873):20172577.

    Article  Google Scholar 

  27. Jurado E, Aguirre O, Flores J, Návar J, Villalón H, Wester D. Germination in tamaulipan thornscrub of north-eastern Mexico. J Arid Environ. 2000;46:413–24.

    Article  Google Scholar 

  28. Jurado E, Estrada E, Moles A. Characterizing plant attributes with particular emphasis on seeds in Tamaulipan thornscrub in semi-arid Mexico. J Arid Environ. 2001;48:309–21.

    Article  Google Scholar 

  29. Jurado E, Návar JH, Villalón H, Pando M. Germination associated with season and sunlight for Tamaulipan thornscrub plants in north-eastern Mexico. J Arid Environ. 2001;49:833–41.

    Article  Google Scholar 

  30. Lindig-Cisneros R. Ecological restoration in Mexico: the challenges of a multicultural megadiverse country. Ecol Restor. 2010;28(3):232–3.

    Article  Google Scholar 

  31. Lohbeck M, Poorter L, Martínez-Ramos M, Rodríguez-Velazquez J, van Breugel M, Bongers F. Functional ecology changing drivers of species dominance during tropical forest succession. Funct Ecol. 2014;28(4):1052–8.

    Article  Google Scholar 

  32. López-Barrera F, Martínez-Garza C, Ceccon E. Ecología de la restauración en México: estado actual y perspectivas. Revista Mexicana de Biodiversidad. 2017;88:97–112.

    Article  Google Scholar 

  33. Maiti R, González R, Kumari CA. Wood carbon and nitrogen of 37 woody shrubs and trees in Tamaulipan thorn scrub, Northeastern Mexico. Pakistan J Bot. 2019;51(3):979–84.

    Article  CAS  Google Scholar 

  34. Marañón T, Grubb PJ. Physiological basis and ecological significance of the seed size and relative growth rate relationship in Mediterranean annuals. Funct Ecol. 1993;7:591–9.

    Article  Google Scholar 

  35. Marroquín-Castillo JJ, Alanís-Rodríguez E, Jiménez-Pérez J, Aguirre-Calderón OA, Mata-Balderas JM, Rubio-Camacho EA, Chávez-Costa AC. Composición florística y diversidad de un área restaurada post-minería en el matorral espinoso tamaulipeco. Acta Botánica Mexicana. 2017;120:7–20.

    Article  Google Scholar 

  36. Martínez DD, Alanís E, Jiménez J, Yerena JI, Gómez MM. Contenido de carbono en un sistema agroforestal pastorilsilvícola del matorral espinoso tamaulipeco. Ciencia UANL. 2015;18(74):58–62.

    Google Scholar 

  37. Mata M, Treviño EJ, Jiménez J, Aguirre OA, Alanís E, Foroughbakhch R. Prácticas de rehabilitación en un ecosistema semiárido afectado por el establecimiento de un banco de material, en el noreste de México. Ciencia UAT. 2014;8(2):32–43.

    Article  Google Scholar 

  38. McAleece N, Lambshead J, Paterson GLJ, and Gage JD. BioDiversity pro, Version 2. The Natural History Museum, London, and The Scottish Association of Marine Science, Oban, Scotland. 1997.

  39. Meli P, Holl KD, Rey Benayas JM, Jones HP, Jones PC, Montoya D. and Moreno Mateos D. A global review of past land use, climate, and active vs. passive restoration effects on forest recovery. Plos one, 2017;12(2).

  40. Mora CA, Alanís E, Jiménez J, González MA, Yerena JI, Cuellar LG. Estructura, composición florística y diversidad del matorral espinoso tamaulipeco, México. Ecología  Aplicada. 2013;12(1):29–34.

    Article  Google Scholar 

  41. Moreno CE. Métodos para medir la biodiversidad. M&T Manuales y Tesis SEA España. 2001;1:84.

    Google Scholar 

  42. Morrison EB, Lindell CA. Active or passive forest restoration? Assessing restoration alternatives with avian foraging behavior. Restor Ecol. 2010;19(201):170–7.

    Article  Google Scholar 

  43. Mostacedo B and Fredericksen TS. Manual de Métodos Básicos de Muestreo y Análisis en Ecología Vegetal. Proyecto de Manejo Forestal Sostenible. Santa Cruz de la Sierra. 2000.

  44. Mueller-Dombois D, Ellemberg H. Aims and methods of vegetation ecology. New York: Wiley; 1974.

    Google Scholar 

  45. Návar J, Mendez E, Graciano J, Dale V, Parresol B. Biomass equations for shrub species of Tamaulipan thornscrub of northeastern Mexico. J Arid Environ. 2004;59:657–74.

    Article  Google Scholar 

  46. Návar JF, Rodríguez-Flores J, Dominguez-Calleros PA, Pérez-Verdin G. Diversity-productivity relationship in the northeastern tamaulipan thornscrub forest of Mexico. Int J Ecol. 2010;1:11.

    Article  Google Scholar 

  47. Ngangyo M, Foroughbakhch R, Carrillo-Parra A, Lidia-Rosaura SC. Estimation of timber production of five species of the tamaulipas thorny shrubs growing in native stands and plantations. Open J Forestry. 2014;4(3):239.

    Article  Google Scholar 

  48. Ngangyo-Heya M, Foroughbakhch-Pournavab R, Carrillo-Parra A, Salas-Cruz LR, Alvarado-Vázquez MA, Bustamante-García V. Foliar biomass production and litterfall pattern of five timber species in forest plantations of semi-arid lands of the northeastern Mexico. Bot Sci. 2017;95(2):295–305.

    Article  Google Scholar 

  49. Panshin AJ, Zeeuw CD. Textbook of Wood Technology. 4th ed. New York: McGraw-Hill; 1980.

    Google Scholar 

  50. Pequeño LM, Alanís E, Jiménez J, González MA, Yerena I, Cuellar LG, Mora A. Análisis de la restauración pasiva post-pecuaria en el matorral espinoso tamaulipeco del noreste de México. Ciencia UAT. 2012;24(2):48–53.

    Article  Google Scholar 

  51. Pequeño-Ledezma MA, Alanís-Rodríguez E, Jiménez-Pérez J, González-Tagle MA, Yerena-Yamallel JI, Cuellar-Rodríguez LG, Mora-Olivo A. Análisis de la restauración pasiva post-pecuaria en el matorral espinoso tamaulipeco del noreste de México. Ciencia UAT. 2012;4(2):48–53.

    Article  Google Scholar 

  52. Pickett STA, McDonnell MJ. Changing perspectives in community dynamics: a theory of successional forces. Trends Ecol Evol. 1989;8(4):241–5.

    Article  Google Scholar 

  53. Pontifes PA, García-Meneses PM, Gómez-Aíza L, Monterroso-Rivas AI, Caso-Chávez M. Land use/land cover change and extreme climatic events in the arid and semi-arid ecoregions of Mexico. Atmósfera. 2018;31(4):355–72.

    Article  Google Scholar 

  54. R core Team . RStudio: Integrated Development for R. Boston, MA: RStudio, Inc. 2015. Available from:

  55. Reid N, Stafford Smith DM, Beyer-Münzel P, Marroquín J. Floristic and structural variation in the Tamaulipan thornscrub, northeastern Mexico. J Veg Sci. 1990;1:529–38.

    Article  Google Scholar 

  56. Rey Benayas JM, Barral MP, Meli P. Lecciones de cuatro meta-análisis globales sobre la restauración de la biodiversidad y los servicios ecosistémicos. Ecol Austral. 2017;27:193–8.

    Article  Google Scholar 

  57. Ruwanza S, Gaertnera M, Eslera KJ, Richardson DM. The effectiveness of active and passive restoration on recovery of indigenous vegetation in riparian zones in the Western Cape, South Africa: A preliminary assessment. S Afr J Bot. 2013;88:132–41.

    Article  Google Scholar 

  58. Sanchez-Castillo L, Kosugi K, Masaoka N, Kubota T. Eco-morphological characteristics of fern species for slope conservation. J Mountain Sci. 2019;16:504–15.

    Article  Google Scholar 

  59. Sarmiento-Muñoz TI, Alanís-Rodríguez E, Mata-Balderas JM, Mora-Olivo A. Estructura y diversidad de la vegetación leñosa en un área de matorral espinoso tamaulipeco con actividad pecuaria en Nuevo León. México CienciaUAT. 2019;14(1):31–44.

    Article  Google Scholar 

  60. Secretaría de Programación y Presupuesto (SPP). 1983. Síntesis Geográfica del Estado de Tamaulipas. Instituto Nacional de Estadística, Geografía e Informática. México, D.F. 158 pp. Available from:

  61. SERI. The SER International Primer on Ecological Restoration. URL 2006.

  62. Shannon CE. A mathematical theory of communication. Bell Syst Techn J. 1948;27(379–423):623–56.

    Article  Google Scholar 

  63. Society of Ecological Restoration. Principios de SER Internacional sobre restauración ecológica. Available from:

  64. Tobón WT, Urquiza-Haas P, Koleff M, Schröter R, Ortega-Álvarez J, Campo R, Lindig-Cisneros J, Sarukhán AB. Restoration planning to guide Aichi targets in a megadiverse country. Conserv Biol. 2017;31(5):1086–97.

    Article  PubMed  Google Scholar 

  65. Vega-López JA, Alanís-Rodríguez E, Molina-Guerra VM, Buendía-Rodríguez E. Selección de especies arbóreas y arbustivas para la restauración del matorral espinoso tamaulipeco. Aridociencia. 2017;2(1):3–10.

    Google Scholar 

  66. Wilson SJ, Rhemtulla JM. Acceleration and novelty: community restoration speeds recovery and transforms species composition in Andean cloud forest. Ecol Appl. 2016;26:203–18.

    Article  PubMed  Google Scholar 

  67. Yerena JI, Jiménez J, Alanís E, Aguirre OA, González MA, Treviño EJ. Dinámica de la captura de carbono en pastizales abandonados del noreste de México. Trop Subtrop Agroecosyst. 2014;17:113–21.

    Google Scholar 

  68. Yerena-Yamallel JI, Jiménez-Pérez J, Aguirre-Calderón OA, Treviño-Garza EJ. Concentración de carbono en la biomasa aérea del matorral espinoso tamaulipeco. Revista Chapingo Serie Ciencias Forestales y del Ambiente. 2011;17:283–91.

    Article  Google Scholar 

  69. Zar JH. Biostatistical analysis. 5th ed. New Jersey: Prentice Hall; 2010. p. 663.

    Google Scholar 

Download references


Not applicable.

Author information

Authors and Affiliations



EAR, VMMG, ACCC, EBR, AMO, LSC, AGAR: all authors conceived the study, designed the experiments, contributed to the data analysis, interpretation, discussion of results and writing. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Laura Sánchez-Castillo.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Appendix I

. Scheme of the plantation process of the assisted ecological succession.

Additional file 2: Appendix II

. Scheme of the plantation process of the assisted natural succession.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Alanís-Rodríguez, E., Molina-Guerra, V.M., Collantes-Chávez-Costa, A. et al. Structure, composition and carbon stocks of woody plant community in assisted and unassisted ecological succession in a Tamaulipan thornscrub, Mexico. Rev. Chil. de Hist. Nat. 94, 6 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: