SciELO - Scientific Electronic Library Online

 
vol.25 número3Influencia de la temperatura e irradiación en la concentración de almidón en Carya illinoinensis K. Koch variedades Wichita y WesternParticipación comunitaria para la medición de carbono en bosques de la Reserva de la Biosfera Mariposa Monarca, México índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista Chapingo serie ciencias forestales y del ambiente

versión On-line ISSN 2007-4018versión impresa ISSN 2007-3828

Rev. Chapingo ser. cienc. for. ambient vol.25 no.3 Chapingo sep./dic. 2019  Epub 19-Feb-2021

https://doi.org/10.5154/r.rchscfa.2018.12.093 

Scientific article

Tree diversity and stored carbon in cocoa (Theobroma cacao L.) agroforestry systems in Soconusco, Chiapas, Mexico

Gicli M. Suárez-Venero1 

Carlos H. Avendaño-Arrazate2  * 

Pablo A. Ruíz-Cruz2 

Paulina Estrada-de los Santos3 

1Universidad de Guantánamo (UG). Av. Che Guevara, km 1.5 carretera Jamaica. C. P. 95100. Guantánamo, Cuba.

2Instituto Nacional de Investigaciones Forestales Agrícolas y Pecuarias (INIFAP), Campo experimental Rosario Izapa. km 18 carretera Tapachula-Cacahoatán. C. P. 30870. Tuxtla Chico, Chiapas, México.

3Instituto Politécnico Nacional (IPN), Escuela Nacional de Ciencias Biológicas. Prolongación de Carpio y Plan de Ayala s/n, col. Santo Tomás. C. P. 11340. Miguel Hidalgo, Ciudad de México, México.


Abstract

Introduction:

Cocoa (Theobroma cacao L.) yields in Chiapas, Mexico remain low, but can be improved with agroforestry designs that consider tree composition to ensure diversity and carbon sequestration.

Objective:

To determine tree diversity and stored carbon at three elevational levels of cocoa agroforestry systems in Soconusco, Chiapas.

Materials and methods:

Representative 50 x 20 m plots were established at three elevational levels (0 to 50 m, 51 to 100 m and ≥101 m). Species were counted and classified taxonomically. Diameter at breast height, height of each species, diversity, similarity between pairs of heights, and carbon storage were estimated.

Results and discussion:

Richness was 35 tree species selected by the producers to shade the cocoa crop. These were grouped into 32 genera and 22 families. Abundance was 199 trees, with abundance per plot and diversity being the highest at ≥101 m. Lower diversity was associated with heights with less equity between abundance and richness. Carbon stored in aerial biomass ranged from 224.9 to 362.1 Mg·ha-1; the amount was highest at lower elevations (0 to 50 m).

Conclusions:

Agroforestry systems had medium to high diversity. Species distribution, by elevational level, responds to the interests of the producers and to the composition of the natural forest. The amount of stored carbon was high, due to the abundance, richness and dominance of the species.

Keywords: Elevational level; species richness; tree composition; carbon sequestration; shade trees

Resumen

Introducción:

Los rendimientos de cacao (Theobroma cacao L.) en Chiapas, México, aún son bajos, pero pueden ser mejorados con diseños agroforestales que consideren la composición arbórea para garantizar la diversidad y captura de carbono.

Objetivo:

Determinar la diversidad arbórea y el carbono almacenado en tres niveles altitudinales de sistemas agroforestales con cacao en el Soconusco, Chiapas.

Materiales y métodos:

Se establecieron parcelas representativas de 50 x 20 m en tres niveles altitudinales (0 a 50 m, 51 a 100 m y ≥101 m). Las especies se contabilizaron y clasificaron taxonómicamente. Se estimaron el diámetro a la altura del pecho, la altura de cada especie, la diversidad y similitud entre pares de alturas, así como el almacenamiento de carbono.

Resultados y discusión:

La riqueza fue de 35 especies arbóreas seleccionadas por los productores para sombra del cultivo de cacao. Estas se agruparon en 32 géneros y 22 familias. La abundancia fue de 199 árboles; la abundancia por parcela y diversidad fueron mayores en la altura ≥101 m. La menor diversidad estuvo asociada a las alturas con menor equidad entre abundancia y riqueza. El carbono almacenado en la biomasa aérea varió de 224.9 a 362.1 Mg·ha-1; la cantidad fue superior a más baja altura (0 a 50 m).

Conclusiones:

Los sistemas agroforestales presentaron diversidad media a alta. La distribución de las especies, por nivel altitudinal, responde a los intereses de los productores y a la composición del bosque natural. La cantidad de carbono almacenado fue alta, debido a la abundancia, riqueza y dominancia de las especies.

Palabras clave: Nivel altitudinal; riqueza de especies; composición arbórea; captura de carbono; árboles de sombra

Introduction

The states of Tabasco and Chiapas are the most important producers of cocoa (Theobroma cacao L.) in Mexico. In Chiapas, the main production regions are the North, Central, Soconusco and Selva-Norte areas, with average yields ranging from 0.20 to 0.54 t·ha-1 (Servicio de Información Agroalimentaria y Pesquera [SIAP], 2018); some of the factors involved in yields are plant health and climate change.

The global climate is changing due to a progressive increase in concentrations of greenhouse gases, such as carbon dioxide (CO2), whose emissions increase every year and contribute to a global increase in temperature (Intergovernmental Panel on Climate Change [IPCC], 2007). Land-use change, deforestation, and massive use of fossil fuels for industrial purposes and transport are the main factors that induce these emissions (Useros, 2013) and affect human health, food security and natural ecosystems (Comisión Económica para América Latina y el Caribe [CEPAL], 2009; Eguren, 2004; IPCC, 2007). Climate change could be mitigated by establishing agroforestry systems (AFS) capable of capturing CO2 from the atmosphere and storing it in aerial and underground biomass for long periods (Gayoso & Guerra, 2005). In this regard, Shibu (2009), Schoeneberger (2009), and Casanova, Caamal, Petit, Solorio, and Castillo (2010) support the fact that AFSs, even if not designed primarily for carbon sequestration, offer an opportunity to increase carbon stocks in the Earth's biosphere.

Most cocoa plantations in the world are established with shade tree species, in some cases under an agroforestry management design. In this context, some authors recognize the need to grow cocoa under shade (Almeida & Valle, 2007; Silva, Orozco, Rayment, & Somarriba, 2013; Somarriba & Quesada, 2005). Therefore, it is necessary to achieve a tree composition for shade within cocoa AFSs, where forest species are arranged in such a way that they allow the entry of light according to the requirements of the crop and positively influence carbon storage.

In Mexico, some studies have been carried out on tree composition in the AFSs of Soconusco, Chiapas (Roa-Romero, Salgado-Mora, & Álvarez-Herrera, 2009; Salgado-Mora, Ibarra, Macías-Sámano, & López-Báez, 2007) and in Tabasco (Ramírez, García, Obrador, Ruiz, & Camacho, 2013); however, such research does not allude to the influence of this tree composition on carbon storage. Therefore, the objective of this research was to determine the tree diversity and stored carbon at three elevational levels in cocoa AFSs in Soconusco, Chiapas, Mexico. The information generated will allow producers to adopt shade management alternatives at the plantations in order to improve the sustainability of these agroecosystems.

Materials and methods

Description of the study area

The research was carried out during the first half of 2018 in three elevational zones that grouped 13 of the 15 municipalities of the Soconusco region, located in the south of the state of Chiapas, Mexico, between coordinates 15° 19' NL and 92° 44' WL covering 4 605.4 km2, which represents 6.28 % of the State area.

Climates in the region are warm and semi-warm; warm sub-humid weather predominates with summer rains, followed by warm humid weather with abundant summer rains. Therefore, according to the Köppen classification modified by García (1973), the predominant climate is type Aw2(w)Ig with average relative humidity of 79.4 % and average annual temperature of 26.8 °C (Comisión Nacional del Agua [CONAGUA], 2015). The Soconusco region and its coastal plain has eight soil types: Litosol, Acrisol, Regosol, Solonchak, Andosol, Luvisol, Nitosol and Cambisol (Comisión Nacional para el Conocimiento y Uso de la Biodiversidad [CONABIO], 2015). However, in cocoa production areas, the predominant soils are generally Andosol, Cambisol and Luvisol.

Sampling

An exploratory tour was carried out and interviews were conducted with the producers in 90 % of the cocoa plantations in each municipality of the Soconusco region, to determine the sampling areas. From this tour, the elevational ranges to be studied, the criteria for the selection of the municipalities to be examined, and the productive entities (farms) to be evaluated were defined.

The selection criteria of the municipalities were chosen based on the representativeness at the State level: 1) area planted with cocoa; 2) number of producers; 3) contribution to production (%); and 4) geographical representativeness by elevational levels (from 0 to 500 m).

The tour identified that the largest area of cocoa plantations is at elevations below 150 m, which ensured representativeness by selecting the largest number of farms at low elevations; therefore, the following elevational levels were defined: 0 to 50 m, 51 to 100 m and ≥101 m, consisting of seven, three and three farms, respectively (Table 1).

Table 1 Locations and heights evaluated in the cocoa agroforestry system in the Soconusco region, Chiapas. 

Municipalities Locations Coordinates X Coordinates Y Elevation (m)
Height 1 (0-50 m)
Huehuetán Rancho Esquipulas 559941 1657255 17
Suchiate Manuel Ávila Camacho 582086 1618980 19
Huixtla El Arenal 554524 1666228 19
Mapastepec La Fronterita 514271 1700423 29
Tuzantán Tercer Cantón 558966 1669483 33
Acapetahua Rancho San Antonio 530773 1692492 34
Escuintla El Triunfo 544935 1697062 36
Height 2 (51-100 m)
Tapachula Raymundo Enriquez 571924 1643444 60
Frontera Hidalgo La Primavera 587577 1635447 83
Metapa Los Cacahuatales 586045 1637418 91
Height 3 (≥101 m)
Acacoyagua Los Cacaos 536672 1701221 398
Tuxtla Chico C. E. Rosario Izapa 590658 1655363 433
Cacahoatán Santa Martha 586167 1658130 500

Response variables

Tree composition

Sampling was carried out in 13 plots of 20 m x 50 m (1 000 m2) corresponding to one plot per municipality, representing a total sampling area of 1.3 ha. In each plot, an inventory was made of the shade species present and the common name of each tree was recorded for taxonomic identification at the family, genus and species level, based on the manual of Central American trees (Barrance et al., 2003) and tropical trees of Mexico (Pennington & Sarukhán, 2005). The basal area of each tree was obtained from the diameter measured at breast height (DBH: 1.3 m above ground level) with a 5 m diameter tape and from the height obtained with a clinometer (SSUNTO brand).

The tree composition and structure of the plots were evaluated at all three heights, by means of the importance value index (IVI) according to Curtis and McIntosh (1951), from the relative density per unit area, the relative basal area and the relative frequency of species, using the equation IVI = relative density + relative dominance + relative frequency. To ensure homogeneity of information, abundance and richness calculations were made per plot at each elevational level.

Diversity and similarity analysis

The diversity of species, genera and families of plants associated with the cocoa AFS at the three selected elevations was analyzed by calculating the Shannon and Simpson indices. For the former, the following equation was used:

H'= -∑ pi(lnpi)

pi = ni/n

where,

H'

Shannon index

ni

number of trees of the ith species

n

total number of trees of all species

∑pi

proportional abundance of the ith species

D = ∑ pi 2

where,

pi

proportion of trees of each species relative to the total.

To convert this probability to a measure of diversity, the complement to Simpson's original measure was used:

1 - D = 1 - ∑pi 2

The similarity between pairs of heights was analyzed with the Jaccard coefficient and the Czekanowski-Sørensen method, qualitative methods explained by Polo (2008). The Jaccard coefficient (Cj) is based on the presence-absence ratio between the number of species in each system and the total number of species (Stiling, 1999), given by the following equation:

Cj=CA+B-C

where,

A

number of species found in system A

B

number of species found in system B

C

number of species common to both systems.

On the other hand, the Czekanowski-Sørensen method is also based on the presence-absence ratio between the number of species shared or not in each system and the total number of species of the two sites under comparison. The following equation was used:

2CS1+S2

where,

C

number of species shared between the two sites

S1

number of species at site 1

S2

number of species at site 2.

Similarity was evaluated with the Morisita-Horn quantitative method, which considers the values ​​of the species shared or not between the two sites under comparison and reflects the similarity of two sites in structure, both in composition and relative abundance. For this, the following equation was used:

2(nia * nib)da+db*(Na*Nb)

where,

nia

number of trees of species i at site a

nib

number of trees of species i at site b

Na

number of trees at site a

Nb

number of trees at site b

For site a:

da=nia2N2a

For site b:

db=Σnib2N2b

Stored carbon

Stored carbon (Mg·ha-1) was determined using the rapid estimation methodology proposed by Segura and Andrade (2008). Only species recorded with DBH1.3 m ≥ 10 cm at each of the studied heights were considered. The species were classified by type of DBH every 5 cm and compared with the tabulated storage values proposed in this methodology for the conditions of cocoa AFSs. The stored carbon values by DBH type were added to a land-use constant for cocoa (17.2). The final carbon storage values were classified by low, medium and high levels, according to the classification of Somarriba, Andrade, Segura, and Villalobos (2008), expressed in the methodology itself.

Statistical analysis

The recorded species data were systematized in a database created with the Excel version 9.1 package, where the frequency and average of species were calculated. With the information originating from the number of trees per species and the abundance to determine the diversity indices, an analysis was made with Fisher’s least significant differences (LSD) test at P ≤ 0.05. The means were compared by multiple range contrast with the statistical package STATGRAPHICS plus version 5.1 (Batanero & Díaz, 2008).

Results

Tree composition

Table 2 shows the 35 tree species selected by producers for shade in cocoa agroforestry plantations in Soconusco, Chiapas. The species were grouped into 32 genera and 22 families. The most represented genera were Citrus (8.5 %) and Tabebuia (5.71 %), which together accounted for 25.12 % of the trees (199 in total). The most represented families were Fabaceae, Moraceae and Rutaceae with four, three and three species, respectively, followed by Apocynaceae, Bignoniaceae, Euphorbiaceae, Malvaceae, Meliaceae and Sapotaceae with two species each. These families grouped 54.1 % of the inventoried species.

Table 2 Species used as shade in cocoa cultivation in 13 sampled municipalities of Soconusco, Chiapas. 

Family Species Common name Uses
Anacardiaceae Mangifera indica L. Mango Fruit
Apocynaceae Aspisperma megalocarpon Müll. Arg Chiche Timber
Apocynaceae Stemmadenia donnell-smithii (Rose) Woodson Chapón Timber
Bignoniaceae Tabebuia donnell-smithii Rose Spring Timber
Bignoniaceae Tabebuia rosea (Bertol.) DC. Pink poui Timber
Boraginaceae Cordia alliodora (Ruiz & Pav.) Oken Laurel Timber
Bombacaceae Ceiba pentandra (L.) Gaertn. Kapok Timber
Burseraceae Bursera simaruba (L.) Sarg. Gumbo-limbo Timber
Cecropiaceae Cecropia obtusifolia Bertol. Trumpet tree Other uses
Chrysobalanaceae Couepia polyandra (Kunth) Rose Ram Timber
Combretaceae Terminalia amazonia (J. F. Gmel.) Exell Guavo Timber
Clusiaceae Garcinia humilis (Vahl) C. D. Adams Achachayru Fruit
Euphorbiaceae Sapium macrocarpum Müll. Arg Chonte Timber
Euphorbiaceae Hippomane mancinella L. Manchineel Timber
Fabaceae Inga micheliana Harms Chalum Timber
Fabaceae Schizolobium parahyba (Vell.) S. F. Blake Zope Timber
Fabaceae Poeppigia procera C. Presl Tepemistle Timber
Fabaceae Erythrina fusca Lour Whistle Other uses
Lauraceae Persea americana Mill. Avocado Fruit
Malvaceae Chiranthodendron pentadactylon Larreat Canaque Timber
Malvaceae Theobroma bicolor Humb. & Bonpl Pate Fruit
Meliaceae Cedrela odorata L. Spanish cedar Timber
Meliaceae Guarea glabra Vahl. Cedrillo Timber
Moraceae Castilla elastica Sesse Panama rubber tree Timber
Moraceae Ficus sp. L. Kill stick Other uses
Moraceae Maclura tinctoria (L.) Steud Old fustic Timber
Muntingiaceae Muntingia calabura L. Capulin Other uses
Rutaceae Citrus nobilis Lour Tangor Fruit
Rutaceae Citrus sinensis L. Orange Fruit
Rutaceae Citrus aurantifolia Swingle Key lime Fruit
Sapindaceae Nephelium lappaceum L. Rambutan Fruit
Sapotaceae Manilkara zapota (L.) van Royen Sapodilla Fruit
Sapotaceae Pouteria sapota (Jacq.) H. E. Moore & Stearn Mamey sapote Fruit
Sterculiaceae Sterculia apelata (Jacq.) Karst Panama tree Timber
Tiliaceae Apeiba tibourbou Aubl. Mico comb Timber

A richness of 35 species was determined at all three elevations. The greatest richness was found at the elevation of ≥101 m with a total record of 21 species which meant, on average, seven species per plot, represented by the families Boraginaceae and Apocynaceae, followed by the elevation between 0 and 50 m with 26 species (on average 3.72 species per plot) represented mainly by the families Boraginaceae, Sapotaceae and Meliaceae; the lowest richness was recorded at elevations between 51 and 100 m with six species in total, averaging two species per plot represented by four families, among them Boraginaceae.

Regarding species abundance, a total of 199 trees were determined in the three elevations combined; 101 trees were recorded within the elevational range from 0 to 50 m, representing 14.42 trees per plot; 47 trees were recorded at the elevational level between 51 and 100 m, representing 15.66 trees per plot, and 51 trees were recorded at heights greater than 101 m, averaging 17 trees per plot. According to Figure 1, the laurel (Cordia alliodora [Ruiz & Pav.] Oken), pink poui (Tabebuia rosea [Bertol.] DC.) and mamey sapote (Pouteria sapota [Jacq.] H. E. Moore & Stearn.) were the most abundant species with 29, 26 and 23 trees, respectively. Of the recorded species, 13 had low abundance with only one tree; of these, eight were found mostly located at the lowest elevational level (0 to 50 m). However, the most abundant species were recorded at all three elevations.

Figure 1 Distribution pattern of species abundance in cocoa agroforestry systems in Soconusco, Chiapas.  

Regarding the IVI, the families Boraginaceae, Fabaceae, Sapotaceae and Meliaceae were among the five most important, which were represented by species located in at least two of the heights studied. The family Boraginaceae was represented by high IVI species at the three elevational levels, which suggests that they are successful taxa with species capable of developing in different environments and heights (Figure 2). In contrast, the family Fabaceae was not represented by any higher IVI species, but it did have species with a lower IVI, such as the zope (Schizolobium parahyba [Vell.] S. F. Blake) and chalum (Inga micheliana Harms), at two elevational levels.

Figure 2 Most important families and species at three elevational levels in Soconusco, Chiapas. IVI: importance value index. 

Unlike the analysis carried out in the families, no important species were recorded developing at the same time at all three height levels, only in two of them (Figure 2). However, some species were found at all three heights, but with lesser or greater importance in one or the other; such is the case of laurel, primavera, mamey sapote, Spanish cedar (Cedrela odorata L.) and avocado (Persea americana Mill.).

Diversity

Table 3 shows that the diversity of species and families, estimated by means of the Shannon index, was significantly (P ≤ 0.05) different among elevational levels; the highest diversity values were found at heights of 51 to 100 m and greater than 101 m. On the other hand, in the diversity estimated through the Simpson index, no significant (P > 0.05) differences were found among the three elevations.

Table 3 Diversity indexes by species and families of the cocoa agroforestry system in the studied heights of Soconusco, Chiapas. 

Height (m) Statistics Species Families
Shannon Index Simpson Index Shannon Index Simpson Index
0-50 Average 1.340 ab 0.670 a 1.330 a 0.660 a
Standard deviation 0.472 0.210 0.510 0.198
Standard error 0.197 0.082 0.190 0.083
Lower limit 0.032 0.542 0.031 0.536
Upper limit 0.647 0.801 1.631 0.797
51-100 Average 0.650 b 0.690 a 0.410 b 0.630 a
Standard deviation 0.701 0.285 0.426 0.324
Standard error 0.323 0.125 0.291 0.127
Lower limit 0.147 0.514 0.041 0.429
Upper limit 1.153 0.909 0.875 0.830
≥101 Average 1.790 a 0.800 a 1.660 a 0.680 a
Standard deviation 0.669 0.148 0.552 0.134
Standard error 0.323 0.125 0.291 0.127
Lower limit 0.290 0.598 2.118 0.486
Upper limit 2.296 0.993 2.156 0.887

Mean values of the Shannon index and Simpson index with the same letter do not present significant statistical differences at elevational level, according to the multiple range test or Fisher’s least significant difference (LSD) test (P > 0.05).

Similarity

According to Table 4, the Jaccard and Sørensen qualitative indices indicate that the greatest similarity was found between tree compositions from 0 to 50 m and ≥101 m; of the total number of species recorded, 12 were common in the three elevational ranges. However, according to Morisita-Horn's quantitative method, the similarity between heights was greater in AFSs located at elevational levels greater than 51 m.

Table 4 Similarity indices of cocoa agroforestry system species by pairs of heights studied in Soconusco, Chiapas. 

Pairs of heights (m) Qualitative methods Quantitative method
Jacard Sørensen Morisita-Horn
0-50 and 51-100 0.2307 0.3750 0.4074
0-50 and ≥101 0.3428 0.5100 0.4325
51-100 and ≥101 0.2272 0.3700 0.4965

Stored carbon

The land used by the cocoa AFSs, in the three elevational ranges studied in Soconusco, stores high levels of carbon with species that accumulate between 3.11 Mg·ha-1 and 205 Mg·ha-1. Based on the information in Table 5, species growing at low and high elevations stored 39.5 % and 36 % of the total, respectively, while at medium elevation (51-100 m) only 24.5 % is stored.

Table 5 Estimation of the stored carbon in the biomass of the cocoa agroforestry system species at the studied heights in Soconusco, Chiapas.  

Height (m) Number of trees (DBH ≥ 10 cm) Stored carbon (Mg·ha-1) Classification
0-50 89 362.1 High
51-100 39 224.9 High
≥ 101 49 329.4 High

DBH: diameter at breast height

Discussion

Tree composition

The greatest abundance per plot was found in the elevational range ≥101 m, mainly due to the number of species recorded in the municipalities of Tuxtla Chico and Acapetagua, with dominance of trumpet tree (Cecropia obtusifolia Bertol.), chapón (Stemmadenia donnell-smithii [Rose] Woodson) and laurel. The tree composition of the three elevational ranges presented typical species of the agroecosystems where cocoa is cultivated, the record of which was similar to that reported in this same region by Salgado et al. (2007) and Roa-Romero et al. (2009).

The family Fabaceae, as one of the most represented, grouped species of the cocoa AFSs, as in the case of chalum (I. micheliana) in the municipalities of Tuxtla Chico, Acacoyagua and Suchiate. This species is commonly found as a shade tree in perennial crops, due to the architecture of its umbrella-shaped crown that allows the homogeneous entrance of light in the plantations. As a service tree it provides nitrogen due to its fixing capacity, and fertility through the pruning residues that are used in the form of dead cover (Barrance et al., 2013).

When examining the regularity of species abundance values when placed in decreasing order, it was observed that the distribution pattern corresponded to the mathematical models explained by Magurran (1989), Margalef (1995) and Krebs (1999), whose curve allowed interpreting the abundance information. This pattern began with a progression similar to a geometric series where few species (13) are dominant and practically rare; that is, a constant proportionality between abundances and species was assumed, so that the series was observed as a straight line in a logarithmic scale. Subsequently, an adjustment to a series with normal logarithmic distribution was observed, where there was a small number of abundant species and a large proportion of species with low abundance, which determined that the curve that it describes has the shape of an inverted jota. This distribution expresses with greater clarity the relationship of trees by species and where species with intermediate abundance were recorded, with respect to the total, and which became more frequent. Finally, the distribution pattern curve showed the conditions of the broken stick model in which species are equally abundant and can be organized into abundance classes.

Six of the 10 species with the highest IVI, at one or more of the elevational levels studied in Soconusco, were found at the elevational level of 51 to 100 m; this is due to the abundance, frequency and dominance of their basal diameter. Mamey sapote was among the species with the highest importance level at the three elevational levels, coinciding with the records of Roa-Romero et al. (2009). This is due to the fact that this species is a fruit tree of great acceptance for its contributions to the economic income of the producer. Commonly, the tree is part of the tree structure of coffee and cocoa plantations, established or naturally regenerated with or without agroforestry design. The location of the mamey sapote is random, as recorded in the plots of Tapachula, Huehuetán and Acapetagua. On the other hand, the mango (Mangifera indica L.), with greater presence in the height from 0 to 50 m, was one of the species with a lower importance value among the 10 most important species; however, this fruit tree is also of great importance to the producers because of the income it provides in light of the main crop’s low prices. Mango was found in cocoa plantations in four municipalities, a situation similar to that recorded by Ramírez et al. (2013) in a 30-year-old cocoa AFS in Tabasco.

Among the timber species that were commonly recorded, at different heights and with a high importance index, is the laurel. Its popularity lies in the value of its wood for local use; it can be combined with annual and perennial crops in an AFS, and it is also used as a medicinal plant (Barrance et al., 2003). The presence of laurel in cocoa AFSs provides an alternative source of income for the producer. Sánchez, Pérez-Flores, Obrador, Sol, and Ruiz-Rosado (2016) also identified this and a large number of other timber species that can potentially be sustainable in cocoa AFSs.

Diversity and similarity

It was assumed that there is low diversity of species and families in the Soconusco region, considering that, in general, Shannon index values range from 1.5 to 3.5 (Zak & Willig, 2004) and that values greater than 3 mean high diversity. The lowest values were calculated at levels between 51 and 100 m, being able to influence the Tapachula municipality plot, whose diversity was null, since it had only one species (mamey sapote) with eight trees. However, according to the calculated Simpson index and considering that the range of values goes from 0 (low diversity) to 1 (high diversity), the cocoa AFSs presented medium diversity with a tendency to high diversity, since the values are higher than 0.50; the highest indices were found at heights greater than 100 m, which is related to their greater richness of species. In this context, it can be deduced that diversity is influenced by a greater height range (>100 m), since the sampled area was the same as in the elevational range from 51 to 100 m, where there was less species diversity and richness. These results may be due to the fact that at heights greater than 100 m there are still many remains of the original forest, whose species diversity has been maintained over time despite the application of shade management techniques for the development of the crop. It also influenced the existence of plantations in municipalities that have greater equity (relationship between richness and relative abundance of each species) and diversity, such as Tuxtla Chico. In this municipality, exploratory tours showed that producers have maintained an adequate diversity of species as shade, with the aim of obtaining additional income, and that they are part of traditional management; these results are similar to those obtained by Roa-Romero et al. (2009). The Simpson index is a measure of the probability that two trees taken randomly belong to the same species and is strongly related to the most abundant species in the sample, which are less sensitive to species richness. Therefore, the closer the index is to unity, the greater the probability of dominance of a species (Polo, 2008).

In the Soconusco region, the low diversity in cocoa AFSs at heights between 51 and 100 m is due to a lower equity between abundance and richness, where an individual-species ratio of 5.7 was calculated; this is not the case at heights greater than 100 m, where a tree composition equivalent to 1.52 trees per species was recorded, followed by the composition recorded at the lowest height studied (0 to 50 m), which was 2.51. Therefore, there is more equity in the cocoa AFSs that are located at heights greater than 100 m, because the ratio between richness and relative abundance of each species is greater.

With respect to the importance of species diversity as shade, the results indicated high variability of illumination in cocoa plantations at all elevational levels, given primarily by the number and type of species per plot, the canopy size of each species and poor management of tree diversity. According to Silva et al. (2013), if this diversity is properly managed, important benefits for the crop can be achieved.

On the other hand, from the qualitative point of view, the similar diversity in tree composition between the lowest and highest height studied is due to the presence of a large number of shared species and to greater diversity. However, from the quantitative point of view, the results of greater similarity between the highest heights may be related to a smaller difference between the number of recorded trees, in addition to greater similarity of the edaphoclimatic conditions with rainfall between 1 600 and 4 000 mm·year-1, an average monthly temperature between 22 and 24 °C and Andosol and Acrisol soil types.

The diversity and similarity results showed that the agroecosystems where cocoa AFSs are developed show signs of low levels of disturbance, despite anthropogenic intervention in these natural environments; the abundance, richness and diversity of species at the three height levels have directly influenced the benefits they bring to producers, regardless of the income earned from the main crop (cocoa). In this sense, producers in the Soconusco region have a diversity of species for timber, energy, medicinal and food uses that make possible the existence of cocoa plantations with economic, ecological and social potential. These similarity measures have already been reported in other cocoa agroforestry environments by Hervé and Vidal (2008) and Zapfack et al. (2002) in Cameroon, who indicate that the type of management is a factor that determines a lower percentage of similarity, a situation inherent to each cocoa region.

According to the results, the establishment of fruit and timber species, based on the interests of the producer for the economic benefits they provide, influenced the tree composition, diversity and similarity of cocoa agroecosystems by elevational levels. This was related to the region’s traditional culture and the intentional selection of producers in search of greater services in conjunction with the main crop. On the other hand, it also influenced the composition of the natural forest since many species originating from the primary forest have developed with a certain dominance and acceptance on the part of the producers. Also, as a result of the climate and soil variation with respect to height, the diversity of the vegetation has changed. In ascending order, we observed medium-height tropical forest, high-height tropical forest, and oak, mountain cloud and pine forests, which are natural vegetation of the Sierra Madre.

Carbon storage

A direct relationship was determined between aerial biomass carbon storage and species abundance. The greatest amount of carbon stored and estimated in the biomass of the species was found at the lowest elevational level (0 to 50 m), due to the greater number of recorded trees per total area, including very abundant species such as pink poui, mamey sapote and Spanish cedar. At the same elevational level (0 to 50 m), dominant species with total DBH of 473 m2·ha-1 were recorded, while at an elevation of 51 to 100 m the lowest abundance and species with lower DBH were recorded, which may have influenced the results. Cerda, Espin, and Cifuentes (2013) determined that carbon storage capacity correlates with the basal area of timber and fruit trees. This means that tree thickness and size is as important as tree abundance and species richness.

In the present research, carbon storage was greater than that determined by Segura and Andrade (2008) in a cocoa AFS in Talamanca, Costa Rica, where they recorded 122 ± 24 Mg C·ha-1; the difference may be due to the tree composition of each territory. In this context, Somarriba et al. (2008) stated that an average total carbon level between 80 and 120 Mg·ha-1 would be the most appropriate in order not to harm cocoa production; quantities above this value are considered high levels. Ortiz, Riascos, and Somarriba (2008) determined that cocoa-laurel systems fixed between 43 and 62 Mg C·ha-1 in 25 years and that in one year they fixed between 1.7 and 2.5 Mg C·ha-1, which corresponds to low storage levels. These same authors refer to research carried out by CATIE (Tropical Agricultural Research and Teaching Center) in Talamanca in 2005, where they report carbon accumulation rates in cocoa plantations, over a 25-year period, of between 2.1 and 2.8 Mg C·ha-1·year-1, whose densities ranged between 100 and 150 trees·ha-1. In the present research, the high volumes of carbon accumulated in the biomass of the trees correspond to densities between 127 and 163 trees·ha-1.

The capacity of agroforestry ecosystems to store carbon in the form of aerial biomass varies according to the age of the trees, and the diameter and height of the shade components. Andrade and Ibrahim (2003) state that AFSs can fix and store between 12 and 228 Mg·ha-1, including soil organic carbon, which represents between 20 and 46 % of the carbon sequestered in primary forests.

From the knowledge of tree composition and species diversity by elevational levels in the cocoa AFSs of Soconusco, producers and decision-makers will be able to evaluate shade management alternatives in plantations, with the aim of designing efficient systems in each edaphoclimatic and physiographic condition. These alternatives may consider species of importance to producers according to their traditions and socio-economic needs, based on the uses of the diverse forest potential in the region. Likewise, strategies with an ecological and agronomic approach can be adopted, with the aim of preserving the diversity of existing species and applying integrated management technology that benefits both cocoa cultivation and forest species, in order to improve the sustainability of these agroecosystems.

Conclusions

The cocoa agroforestry systems in Soconusco, Chiapas had medium to high tree composition and diversity, whose distribution by elevational levels responds to the interests of producers and to the composition of the natural forest in each edaphoclimatic condition. The amount of carbon stored in the aerial biomass in the cocoa agroecosystems was high, mainly due to the abundance, richness and dominance of the species selected by the producer.

References

Almeida, A., & Valle, R. (2007). Ecophysiology of the cocoa tree. Brazilian Journal of Plant Physiology, 19(4), 425-448. doi: 10.1590/S1677-04202007000400011 [ Links ]

Andrade, H., & Ibrahim, M. (2003). ¿Cómo monitorear carbono en sistemas silvopastoriles? Agroforestería en las Américas, 10, 109-116. Retrieved from http://www.fao.org/tempref/docrep/nonfao/lead/x6378s/x6378s00.pdfLinks ]

Barrance, A., Beer, J., Boshier, O. H., Chamberlain, J., Cordero, J., Detlefsen, G., ...Stewart, J. (2003). Árboles de Centroamérica. Un manual para extensionistas. Costa Rica: OFI/CATIE. Retrieved from https://www.arbolesdecentroamerica.infoLinks ]

Batanero, C. & Díaz, M. C. (2008). Análisis de datos con Statgraphics. Granada: Universidad de Granada. Retrieved from https://www.ugr.es/~batanero/pages/ARTICULOS/anadatos.pdfLinks ]

Casanova, L. F., Caamal, M. J. A., Petit, A. J. C., Solorio, S. F. J., & Castillo, C. J. B. (2010). Acumulación de carbono en la biomasa de Leucaena leucocephala y Guazuma ulmifolia asociadas y en monocultivo. Revista Forestal Venezolana, 54(1), 45-50. Retrieved from http://www.saber.ula.ve/bitstream/handle/123456789/31644/articulo4.pdf;jsessionid=6F015D28B656479526937EC3CF352132?sequence=1Links ]

Cerda, R., Espin, T., & Cifuentes, M. (2013). Carbono en sistemas agroforestales de cacao de la Reserva Indígena Bribri de Talamanca, Costa Rica. Agroforestería de las Américas, 49, 33-41 Retrieved from http://infocafes.com/portal/wp-content/uploads/2017/01/345.pdfLinks ]

Comisión Económica para América Latina y el Caribe (CEPAL). (2009). Cambio climático y desarrollo en América Latina y el Caribe: una reseña. Retrieved from https://www.cepal.org/es/publicaciones/3640-cambio-climatico-desarrollo-america-latina-caribe-resenaLinks ]

Comisión Nacional del Agua (CONAGUA). (2015). Actualización de la disponibilidad media anual de agua en el acuífero Soconusco (0710) estado de Chiapas. Retrieved from http://www.gob.mx/cms/uploads/attachment/file/102911/DR_0710.pdfLinks ]

Comisión Nacional para el Conocimiento y Uso de la Biodiversidad (CONABIO). (2015). Región del Soconusco. Retrieved November 22, 2018 from http://www.conabio.gob.mx/conocimiento/regionalizacion/doctos/rhp_032.htmlLinks ]

Curtis, J. T., & McIntosh, R. P. (1951). An upland forest continuum in the prairie-forest border region of Wisconsin. Ecology, 32(3), 476-496. doi: 10.2307/1931725 [ Links ]

Eguren, L. (2004). El mercado de carbono en América Latina y el Caribe: balance y perspectivas. Santiago, Chile: CEPAL. Retrieved from https://repositorio.cepal.org/bitstream/handle/11362/5620/1/S043136_es.pdfLinks ]

García, E. (1973). Modificaciones del sistema de clasificación climática de Köppen para adaptarlo a las condiciones de la república mexicana (2.a ed.). México: Instituto de Geografía, Universidad Nacional Autónoma de México (UNAM). Retrieved from http://www.igeograf.unam.mx/sigg/utilidades/docs/pdfs/publicaciones/geo_siglo21/serie_lib/modific_al_sis.pdfLinks ]

Gayoso, J., & Guerra, J. (2005). Contenido de carbono en la biomasa aérea de bosques nativos en Chile. Bosque, 26(2), 33-38. doi: 10.4067/S0717-92002005000200005 [ Links ]

Hervé, B. D., & Vidal, S. (2008). Plant biodiversity and vegetation structure in traditional cocoa forest gardens in southern Cameroon under different management. Biodiversity and Conservation, 17(8), 1821-1835. doi: 10.1007/s10531-007-9276-1 [ Links ]

Intergovernmental Panel on Climate Change (IPCC). (2007). Synthesis report: Climate change 2007. Switzerland: Author. [ Links ]

Krebs, J. (1999). Ecological methodology. USA: Edic. Addison Wesley Longman. [ Links ]

Magurran, A. E. (1989). Diversidad ecológica y su medición. Barcelona, España: Ed. Vedra. Retrieved from https://dialnet.unirioja.es/servlet/libro?codigo=150880Links ]

Margalef, R. (1995). Ecología. Barcelona, España: Ed. Omega. Retrieved from http://www.ediciones-omega.es/ecologia/47-ecologia-978-84-282-0405-7.htmlLinks ]

Ortiz, A., Riascos, L., & Somarriba, E. (2008). Almacenamiento y tasas de fijación de biomasa y carbono en sistemas agroforestales de cacao (Theobroma cacao L.) y laurel (Cordia alliodora). Agroforestería de las Américas, 46, 26-29. Retrieved from http://repositorio.bibliotecaorton.catie.ac.cr/bitstream/handle/11554/5764/Almacenamiento_y_tasas_de_fijacion_de_biomasa.pdf?sequence=1&isAllowed=yLinks ]

Pennington, T. D., & Sarukhán, K. (2005). Árboles tropicales de México: manual para la identificación de las principales especies. México: UNAM. [ Links ]

Polo, C. S. (2008). Índices más comunes en biología. Segunda parte, similaridad y riqueza beta y gamma. Revista Facultad de Ciencias Básicas, 4(1), 135-142. doi: 10.18359/rfcb.2239 [ Links ]

Ramírez, A., García, E., Obrador, J. J., Ruiz, O., & Camacho, W. (2013). Diversidad florística en plantaciones agroforestales de cacao en Cárdenas, Tabasco, México. Universidad y Ciencia, 29(3), 215-230. Retrieved from http://www.scielo.org.mx/scielo.php?script=sci_arttext&pid=S0186-29792013000300001Links ]

Roa-Romero, H. A., Salgado-Mora, M. G., & Álvarez-Herrera, J. (2009). Análisis de la estructura arbórea del sistema agroforestal de cacao (Theobroma cacao L.) en el Soconusco, Chiapas - México. Acta Biológica Colombiana, 14(3), 97-110. doi: 10.15446/abc [ Links ]

Salgado-Mora, M. G., Ibarra, G., Macías-Sámano, J. E., & López-Báez, O. (2007). Diversidad arbórea en cacaotales del Soconusco, Chiapas, México. Interciencia, 3(11), 763-768. Retrieved from http://www.redalyc.org/pdf/339/33901107.pdfLinks ]

Sánchez, F., Pérez-Flores, J., Obrador, J. J., Sol, A., & Ruiz-Rosado, O. (2016). Árboles maderables en el sistema agroforestal de cacao en Cárdenas, Tabasco, México. Revista Mexicana de Ciencias Agrícolas, 14, 2711-2723. doi: 10.29312/remexca.v0i14.440 [ Links ]

Schoeneberger, M. (2009). Agroforestry: working trees for sequestering carbon on agricultural lands. Agroforestry Systems, 75, 27-37. doi: 10.1007/s10457-008-9123-8 [ Links ]

Shibu, J. (2009). Agroforestry for ecosystem services and environmental benefits: an overview. Agroforestry Systems, 76, 1-10. doi: 10.1007/s10457-009-9229-7 [ Links ]

Segura, M., & Andrade, H. (2008). ¿Cómo estimar rápidamente el carbono almacenado en la biomasa aérea de los sistemas agroforestales indígenas de Talamanca, Costa Rica? Revista Agroforestería de las Américas, 46, 97-103. Retrieved from http://www.worldcocoafoundation.org/wp-content/uploads/files_mf/segura2008b.pdfLinks ]

Servicio de Información Agroalimentaria y Pesquera (SIAP). (2018). Intención de cosecha por estados de México en ciclos perennes. Retrieved from http://infosiap.siap.gob.mx/opt/agricultura/intension/Intencion_cosechaPerenne_estado2018.pdfLinks ]

Silva, C., Orozco, L., Rayment, M., & Somarriba, E. (2013). Conocimiento local sobre los atributos deseables de los árboles y el manejo del dosel de sombra en los cacaotales de Waslala, Nicaragua. Revista Agroforestería de las Américas, 49, 51-60. Retrieved from http://repositorio.bibliotecaorton.catie.ac.cr/bitstream/handle/11554/7050/RAFA_49_Completa.pdf;sequence=3Links ]

Somarriba, E., Andrade, H., Segura, M., & Villalobos, M. (2008) ¿Cómo fijar carbono atmosférico, certificarlo y venderlo para complementarlos ingresos de productores indígenas de Costa Rica? Revista Agroforestería en las Américas, 46, 81-88. Retrieved from http://orton.catie.ac.cr/repdoc/A3126e/A3126e.pdfLinks ]

Somarriba, E., & Quesada, F. (2005). El diseño y manejo de la sombra en el cacaotal. Costa Rica: CATIE. Retrieved from http://www.aprocacaho.com/wp-content/uploads/2013/07/SOMBRA-DEL-CACAO-CATIE.pdfLinks ]

Stiling, P. D. (1999). Ecology: Theories and applications. USA: Ed. Prentice Hall. [ Links ]

Useros, J. L. (2013). El cambio climático: Sus causas y efectos medioambientales. Anales de la Real Academia de Medicina y Cirugía de Valladolid, 50, 71-98. Retrieved from https://dialnet.unirioja.es/descarga/articulo/4817473.pdfLinks ]

Zak, J. C., & Willig, M. R. (2004). Fungal biodiversity patterns. In G. M. Mueller, G. Bills, & M. S. Foster (Eds.), Biodiversity of fungi: Inventory and monitoring methods (pp. 59-75). Amsterdam: Elsevier, Academic Press. [ Links ]

Zapfack, L., Engwald, S. S., Achoundong, G., & Madong, B. (2002). The impact of land conversion on plant biodiversity in the forest zone of Cameroon. Biodiversity and Conservation, 11(11), 2047-2061. doi: 10.1023/A:1020861925294 [ Links ]

Received: December 15, 2018; Accepted: May 30, 2019

*Corresponding author: avendano.carlos@inifap.gob.mx, tel.: +52 (962) 121 0383.

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License