Composition, structure and diversity of a mesquite in Pesquería (Northeastern Mexico)

Abstract

Background

Although the mesquite (mesquital or mezquital in Spanish) is one of the representative ecosystems of the landscapes in the north of Mexico, it is also one of the least studied. This study evaluated the structure (horizontal and vertical) and diversity of a plant community of mesquite in Northeastern Mexico. Three plots of 1,600 m² each were established. All trees and shrubs with a basal diameter (d_0.10) ≥ 0.5 cm were recorded, and total height (h) and crown diameter (d_crown) were measured.

Results

There were 8 families, 12 genera and 14 species. The genus presenting the most species was Acacia (three species). The most representative family was Fabaceae with seven species. The evaluated community presents a density of 375 N/ha and a crown area of 6,600 m²/ha. The species with the highest values on the Importance Value Index (IVI) were Prosopis glandulosa (15.95%), Acacia amentacea (14.50%), Havardia pallens (14.27%) and Acacia farnesiana (11.22%). These four species account for 55.94% of IVI. The value obtained from the Vertical Species Profile Index (A) was 3.03, with an A_max of 3.74 and an A_rel of 81.15%, indicating high structural diversity in the high strata. The evaluated plant community had a Margalef Diversity Index value of D _Mg  = 2.50 and a Shannon Index value of H′ = 2.28, values which are intermediate and considered to be common in the scrublands of Northeastern Mexico.

Conclusions

1) The studied community presents intermediate values that are considered as common in comparison to other arid and semi-arid vegetation associations of Northeastern Mexico. 2) The abundance curve of the species was well adjusted to the geometric model, and the distribution is associated with adverse environments such as semi-arid. 3) The family with greater importance for its contribution to the community is Fabaceae, while the genus with more species was Acacia. The research generated quantitative information of the plant community of a mesquite which is in a phase of mature ecological succession.

Background

Mesquites are spiny trees or shrubs of the pea family that grow preferentially on flat deep soils in arid and semi-arid zones of Asia, Africa and mainly America [1,2,3,4]. These vegetal communities have been called different names such as mesquite forest, mesquite woodland, thorny forest, thorny deciduous forest, high thorny scrub or low thorny evergreen forest [4,5,6,7]. However, what predominantly characterizes mesquites is the dominance of a species of the genus Prosopis in their upper strata [8].

The distribution of mesquites within the New World and particularly Mexico is irregular due to the environmental and geographical conditions wherever they are present. In some cases, they occupy reduced and isolated areas when associated with other types of vegetation such as saltbush, thorny scrub, tropical forest, or submontane scrub [9]. Furthermore, it is also located in extensive and continuous areas such as the Coastal Plains of Mexico [4].

In the Northeastern Mexico, mesquites have been described in Tamaulipas [10], Coahuila [8], and Nuevo León [11,12,13]. Such communities with Prosopis glandulosa Torr., P. laevigata (Humb. & Bonpl. ex Willd.) M.C. Johnst. or P. tamaulipana Burkart have been studied because of their importance from the ecological point of view and the potential use for their flora. Unfortunately, it is common for such work to note the advanced state of deterioration of mesquites that is attributed to changes in and use for agriculture and livestock, to the exploitation of Prosopis spp. as raw material for timber, fuel, fodder and other uses [2].

For the case of Nuevo Leon, the descriptions of the mesquites are more related to floristic attributes or forestry point of view, which are the reasons data on community structure are scarce. Although Briones and Villarreal [11] describe a thorny thicket of Prosopis-Acacia in the north of the state, this does not include ecological data that describe the structure of this plant community in detail. In this context, and considering the need of studies that describe in detail the current state and the structure of the mesquites that still persist in the Northeastern Mexico, this work was carried out. The objective of this study was to know the structure (horizontal and vertical) and floristic diversity of a mesquite community located in Northeastern Mexico, and particularly in the center of the state of Nuevo Leon.

Methods

Area of study

The present research was carried out in the conservation area of the Ternium Industrial Center of México, which has an area of 46 hectares and is located within the municipality of Pesquería, Nuevo León, in Northeastern Mexico (Fig. 1). It is located between latitude 25° 45’ N and longitude 99° 58’ W, at an average altitude of 306 masl, which belongs to the physiographic region of the North Gulf Coastal Plain [6]. The predominant climate is very dry and semi-warm (Bwhw), with an average annual temperature in the range from 20 to 21 °C. The types of soils present in the majority of cases are xerosol, castañozem, feozem, regosol and in the minority of cases are fluvisol, vertisol and rendzine. The average annual rainfall is 550 mm. The vegetation of the area corresponds to a mesquital with a history of cattle and hunting use and is currently disturbed [14].

Fig. 1

Location of the study area. On the left side, the Mexican Republic is indicated, on the upper right side the superior right the State of Nuevo Leon is indicated, and on the lower right side the Municipality of Pesquería is indicated

Analysis of the vegetation

In order to fulfill the objective, the “c” sub basin of the Pesquería River with in a low part of the San Juan River Basin (Rio Bravo, Rio Bravo) was selected, where the conservation area is located and where intermittent surface runoff occurs. Vegetation evaluation was carried out using three randomly distributed sampling sites. The dimensions of each site was 40 m × 40 m (1,600 m²). A census was taken of all shrub and tree species with a basal diameter (d_0.10) ≥ 5 cm. The species of each individual was identified and recorded, taking the measurements of total height (h) and crown diameter (d_crown) in north-south and east-west directions. For the nomenclature of families, orders and species we followed the APG III [15], and the scientific names and families were corroborated in the database of Tropicos [16].

Data analysis

Structure

In order to evaluate the horizontal structure of the species in the study community, we used the following structural variables: abundance, dominance, frequency, with which we calculated the Importance Value Index (IVI) which was calculated from the following mathematical equations [17, 18]:

$$ {A}_i= N i/ S $$

$$ A{R}_i = \left(\raisebox{1ex}{${A}_i$}\!\left/ \!\raisebox{-1ex}{${\displaystyle \sum }{A}_i$}\right.\right)*100 $$

$$ i= 1\dots n $$

Where AR _i is the relative abundance of species i, with respect to total abundance (A _i ); Ni is the number of individuals of species i, and S is the surface (ha).

To estimate the relative dominance we used:

$$ {D}_i = \raisebox{1ex}{$ A{b}_i$}\!\left/ \!\raisebox{-1ex}{$ S(ha)$}\right. $$

$$ D{R}_i = \left(\raisebox{1ex}{${D}_i$}\!\left/ \!\raisebox{-1ex}{${\displaystyle \sum }{D}_i$}\right.\right)*100 $$

$$ i= 1\dots n $$

Where DR _i is the relative dominance of species i, with respect to total dominance (D _i ); Ab _i is the crown area of species i, and S is the surface (ha).

$$ {F}_i = \raisebox{1ex}{${P}_i$}\!\left/ \!\raisebox{-1ex}{$ NS$}\right. $$

$$ F{R}_i = \left(\raisebox{1ex}{${F}_i$}\!\left/ \!\raisebox{-1ex}{${\displaystyle \sum }{F}_i$}\right.\right)*100 $$

$$ i= 1\dots n $$

Where FR _i is the relative frequency of species i with respect to the total frequency (F _i ); P _i is the frequency of species i at sampling sites, and NS is the total number of sampling sites. The Importance Value Index (IVI) is defined as:

$$ IVI = \frac{{\displaystyle {\sum}_n^{i=1}}\left( A{R}_i, D{R}_i, F{R}_i\right)}{3} $$

where AR _i is the relative abundance; DR _i is the relative dominance, and FR _i is the relative frequency.

A graph of height classes was generated in order to evaluate the vertical structure of the community. The Vertical Distribution Index of Pretzsch [19] was calculated for three zones of height: zone I: 80–100% of the maximum height of the population, zone II: 50–80% of the maximum height of the population, zone III: 0–50% of the maximum height of the population [19]. In this study was the high strata (7.20 – 9.00 m), medium strata (4.50 – 7.19 m) and low strata (<4.50 m). The Vertical Distribution Index was calculated according to the following mathematical formula:

$$ A= - {\displaystyle \sum_{i=1}^S}{\displaystyle \sum_{j=1}^Z}{p}_{i j}*\ \ln \left({p}_{i j}\right) $$

Where S is the number of species present; Z is the number of height zones and pij is the proportion of species in each height zone:

$$ p i j= n i j/ N $$

where nij is the number of individuals of the same species (i) in the zone (j) and N is total number of individuals.

In order to compare the Pretzsch Index it is necessary to standardize it and this is undertaken by the value of A _max , which is calculated in the following manner:

$$ Amax= In\left( S* Z\right) $$

Then the value of A can be standardized according to:

$$ Arel=\frac{A}{In\left( S* Z\right)}*100 $$

Diversity

To estimate species diversity, the Shannon Index [20] and the Margalef index [21], respectively, were estimated. The Shannon Index was estimated by using the following equation:

$$ {H}^{\prime }=-{\displaystyle \sum_{i=1}^S}{p}_i* \ln \left({p}_i\right) $$

where S is the number of species present, ln is natural logarithm and pi is the proportion of species. P _i   = n _i /N, where n _i is the number of individuals of species i and N is the total number of individuals. With the same meaning of the variables being common, the Margalef Diversity Index (D _a ) was estimated using the following equation:

$$ Da = \frac{\left( S-1\right)}{ \log N} $$

Species abundance curves

Species density was analyzed using species abundance curves. With these curves it is possible to make inferences about the state of the ecosystems, besides using for descriptions in the form of mathematical models.

In this study, the species abundance curves were fitted to known mathematical models. Currently, there are many models that are used to describe species diversity in a community. However, in this work only three of the best fit models are analyzed: the geometric model [22], the Poisson model of the logarithmic normal series [23], and the Neutral Model of Alonso and Mckane [24, 25].

For the selection of the best model, the Akaike Information Criterion (AIC) was used to compare the selected models, taking into account their fit and complexity. When comparing models using this method, the selection of the best model is based on the lowest value in the AIC. In addition, we used the delta AIC criterion (dAIC) which, when it has a value that is less than 2, indicates that the comparative models similarly explain the trend of the data, (that is, there are no differences between one and the other). To determine the goodness of fit of the models χ ² was used, as recommended by Magurran [26].

Adjusting the models

The models were adjusted using the maximum likelihood method with software R version 3.1.2 [17], with the support of RStudio version 0.99 [27] and also running routines by Prado et al. [25].

Results

Composition

The presence of 8 families, 12 genera and 14 species was recorded. The genus with more species present in the study area was Acacia with three species. The most representative family was Fabaceae with seven species, and the rest presented one species each (Table 1). According to the biological form, nine species are shrubs and five species are trees.

Table 1 Scientific name, common name, family and biological form of the registered species

Community structure

The abundance of the evaluated plant species was 375 N/ha. The most abundant species within the community were Acacia amentacea with 85 N/ha, Havardia pallens with 73 N/ha, and Prosopis glandulosa with 42 N/ha. When added these three species equal 200 N/ha, representing 53.33% of the total abundance of the evaluated plant community.

The evaluated community crown area was 6,600 m²/ha, representing 66.00% coverage. The species that presented the greatest dominance were Prosopis glandulosa with 26.30%, Acacia farnesiana with 17.42%, and Havardia pallens with 13.00%, amounting to 56.72% of the total coverage. The most frequent species had some presence in the three evaluated sites, and were Acacia amentacea, Cordia boissieri, Havardia pallens, Prosopis glandulosa, and Sideroxylon celastrinum.

Importance value

The species with the highest Importance Index were Prosopis glandulosa (15.95%), Acacia amentacea (14.50%), Havardia pallens (14.27%), and Acacia farnesiana (11.22%; Table 2). These four species account for 55.94% of IVI.

Table 2 Abundance, dominance, frequency, and Importance Value Index (IVI) of the species recorded (ranked according to their IVI value)

Vertical structure

The total height of the evaluated individuals fluctuated between 2.3 m and 9.00 m. The highest abundance of individuals is represented in the height class 5–7 m with 173 N/ha which is equivalent to 46% of the evaluated population, followed by class 3–5 m with 118 N/ha representing 31% of the evaluated population (Fig. 2).

Fig. 2

Absolute abundance (N/ha) of individuals according to height class (m)

According to the classification proposed by Pretzsch, where three high strata are established, the high stratum is composed of 60 N/ha (16.11% of the population) belonging to five species, the mean stratum by 221 N/ha (58.89% of the population) of 12 species, and the low stratum by 94 N/ha (25% of the population) with 13 species.

The species found in the three strata of height are Acacia amentacea, Acacia farnesiana, Havardia pallens, Parkinsonia texana, and Prosopis glandulosa. Those species that are well represented in the vertical structure of the community are also those that presented the highest IVI values in the horizontal structure.

The value obtained from the Vertical Species Profile Index (A) was 3.03, with an A _max of 3.74, and an A _rel of 81.15%, indicating high vertical structural diversity. The values of A _rel are close to 100% indicating that species and abundances are evenly distributed in the three height strata (Table 3).

Table 3 Absolute Abundance (N/ha) and proportional abundance (of the total, and in the zone) with respect to height strata of the registered species. High strata (7.20 – 9.00 m), medium strata (4.50 – 7.19 m) and low strata (< 4.50 m)

Species diversity

The evaluated plant community had a Margalef Diversity Index value of D _Mg  = 2.50 and a Shannon Index value of H′ = 2.28.

Species rank abundance curves

The known models were fitted and the best was a geometric model (Fig. 3), followed by the Poisson LogNormal model (AIC = 123.86, dAIC = 2.34, x ² = 7.5581, df = 13, p = 0.87), which is followed by the Neutral model (AIC = 131.25, dAIC = 9.72, x ² = 184.89; df = 13, p < 0.001).

Fig. 3

Range of species abundance in the study area, where: N/ha = Number of trees per hectare, AIC = Akaike’s Information Criterion, χ ² = Chi-square test, df = degree of freedom test χ ², p = probability value of the test χ ²

Discussion

Composition

The composition observed in our study resembles that reported by Briones and Villarreal [11], in the thorny bush of Prosopis-Acacia. According to these authors, the composition and physiognomy correspond to the extra-desert mesquite described by Rezedowski [4] for the state of San Luis Potosí, and the spiny and evergreen forest [7]. However, the composition differs from other studies in which Prosopis sp. (mesquite) dominates the arboreal strata in a mono-specific manner [8]. There are associations similar to those reported in our study, but with lower heights, which are described by Briones and Villarreal [11], such as the medium thorny scrub of Prosopis-Acacia.

The most representative family of the mesquite was Fabaceae, and the others families only presented one species. This record is similar to that obtained by Montaño et al. [28], who report the same number of families in spite of having a greater number of species, also with the most representative being the Fabaceae, a characteristic situation for the xerophilous scrub of this region. The most representative species (Prosopis glandulosa, Acacia amentacea, Havardia pallens, and Acacia farnesiana) in this study are reported as being characteristic of this type of community [29]. The mesquite woodland includes forests dominated by Prosopis glandulosa var. torreyana and var. glandulosa, which are recorded in areas that have been disturbed [30], and these are frequently associated with species such as Forestiera angustifolia, Leucophyllum frutescens, Ziziphus obtusifolia, Acacia rigidula (= A. amentacea), A. berlandieri, Celtis pallida, and Karwinskia humboldtiana [13], the same species as recorded in the present investigation.

Structure

The three species with high abundance were Acacia amentacea, Havardia pallens, and Prosopis glandulosa, and these results are in agreeance with those of Estrada et al. [29], who notes that these species are the most abundant. Montaño et al. [28] note that in Mexico’s xerophilous scrub, mesquite is an abundant species and possibly a key species in disturbed scrubland. Rojas-Mendoza [13], INEGI [6] and Estrada et al. [29] note that it is a species that resists the effects of disturbance and is thus associated with disturbed zones.

The evaluated community has a coverage of 66% (6,600 m²/ha), which means that 34% of the area is devoid of vegetation. This may be due to the fact that only individuals greater than 5 cm in diameter were considered, since for similar communities in Northeastern Mexico that considered smaller individuals, coverage greater than 100% (13,973 m²/ha) was recorded, which denotes a considerable overlap of crowns [31].

The most frequent species are the genus Prosopis and Acacia, similar to those reported by Montaño et al. [28], being of greater presence in the mesquites. Similarly, Mora-Donjuan et al. [32], note that the genus Prosopis has a greater frequency in the microphyllous desert scrublands.

Vertical structure

The total height of the evaluated individuals fluctuated between 2.3 m and 9.00 m. This height is considered high compared to other plant communities of the xerophilous scrub [32]. The value obtained from the Vertical Species Profile Index (A) was 3.03, with an A _max of 3.74 and an A _rel of 81.15%, indicating high vertical structural diversity. Pretzsch [33] notes out that high vertical and horizontal heterogeneity are strongly related to high diversity and ecological stability. This suggests that the species and their abundances are evenly distributed in the three height strata, and that the system is in a good state of conservation. These results are better than those presented by Alanís-Rodríguez et al. [34], Jiménez-Pérez et al. [35], Villavicencio et al. [36], and Mora-Donjuan et al. [32] who determined that the evaluated xerophilous scrub communities, many of which were previously subjected to some type and intensity of use, present only one or two highly dominant strata.

The presence of Acacia amentacea and Parkinsonia texana in the three strata shows the level of adaptation required for these environments. These species are well adapted to water stress, and they use mechanisms to avoid dehydration of their tissues, and to adjust their morpho-physiological characteristics to cope with the long season of dryness, as well as the fixation of active nitrogen [30], which also explains the possibility of their behaving as colonizing species.

Species diversity

The specific richness that was recorded differs from that documented by Montaño et al. [28], in a larger sampling area which registered 36 species, 14 of which were documented in the present study. This situation is presented by the environmental factors in each research area, Montaño et al. [28] studied a semi-arid spiny scrubland of the Mezquital Valley, where precipitation of 520 mm and temperature range of 16–24 °C were similar to the present study precipitation 550 mm and temperature range of 20–21 °C. The altitude factor positively influences the development of the species, whereby Montaño et al. [28] recorded an altitude of 2013 m.a.s.l compared to the altitude of 306 m.a.s.l in our study. This suggests that environmental factors play a positive role in the species diversity of the study area, although most of the recorded species are distributed exclusively or preferably in the arid or semi-arid zones of Mexico [37, 4].

The low number of species is associated with that the mesquites harboring the lowest number of taxa for the northern zone of the state [29], being that this community has the smaller registered presence of Caesalpiniaceae and Fabaceae. This is largely limited to edaphic features, and mesquites are distributed in the northwestern end of the state, where sandy soils dominate [29].

Regarding the diversity of the community described in the present study, the data presented in other studies do not account for the properties of the studied communities, and therefore, the comparison with other mesquite communities is not feasible. However, the evaluated area showed values of D _Mg  = 2.50 and H′ = 2.28, and these values are similar to those recorded for other xerophilic scrubs. Mora-Donjuan et al. [31] evaluated a reference area of Tamaulipan thorn scrub (TTS) and recorded D _Mg  = 2.26 and Diversity H = 1.94. They are also similar to those recorded by Jiménez-Pérez et al. [38], who evaluated the diversity of a regenerated TTS community with agricultural history in Northeastern Mexico and obtained values of D _Mg  = 2.17 and alpha diversity H′ = > 2.27. The values recorded in this study are greater than those of Pequeño-Ledezma et al. [39], who registered values of D _Mg  = 1.40 and alpha diversity of H = 1.27 in a post-livestock regeneration area in the Tamaulipan thorn scrub.

As reported in the existing mesquite studies, these communities are mostly reported as monospecific to Prosopis sp., or floristically not very diverse but mainly dominated by Prosopis and Acacia. This qualitative observation, as reported in other studies (e.g. [11, 8]), is consistent with the quantitative results obtained in our study, where high values for Prosopis glandulosa, Acacia amentacea, Havardia pallens, and Acacia farnesiana - the species that distinguish the association. The dominance of these taxa over others is also observed in the species abundance curve of the studied community which was adjusted to the geometric model.

Species abundance curves

In the present study, the fit to the geometric model suggests a typical structure of strongly stressful environments. Magurran [26] describes the species abundance model adjusted to the geometric model as the model associated with adverse environments or initial successive stages, where the limiting conditions prevail. In the case of mesquite, the environment corresponds to that of arid zone [4]. The dominant species (Prosopis glandulosa, Acacia amentacea, Havardia pallens, and Acacia farnesiana) present different adaptations to conditions of water stress and light saturation [40,41,42].

Conclusions

According to the results obtained in the present study, the following conclusions are highlighted. 1) The studied community presents intermediate values that are considered common in comparison to other arid and semi-arid vegetation associations of Northeastern Mexico. 2) The abundance curve of the species was well adjusted to the geometric model, and the distribution is associated with adverse environments such as semi-arid. 3) The family with greater importance for its contribution to the community is Fabaceae, while the genus with more species was Acacia. The research generated quantitative information of the plant community of a mesquite which is in a phase of mature ecological succession.