More than three-quarters of the world's terrestrial biodiversity occurs in the forests and constitutes an invaluable socio-economical resource for millions of people worldwide, particularly in rural areas (FAO 2018). On a global level, forests are under pressure from biological processes and anthropic exploitation, even when they are located in Protected Natural Areas, Mexico is no exception (FAO 2018, Regil-García et al. 2020). Without proper management, anthropogenic pressures can produce adverse effects on forest structure and functionality (Gómez Díaz &Villalobos 2020).
In Mexico, coniferous forests have a wide floristic and ecological diversity, that commonly occur in temperate and cold climates (Rzedowski 2006). In particular, pine forests cover 38 % (6.15 million ha) of the country's forested land surface (Quintero-Gradilla et al. 2019). They are widely distributed throughout the mountain ranges of the Sierra Madre Occidental, the Sierra Madre Oriental, the Trans-Mexican Volcanic Belt, the Sierra Madre del Sur, and the Sierra de Oaxaca, covering a latitudinal gradient from Baja California to the Yucatan peninsula (Sánchez-González 2008). The pines are distributed from 100 to 4,200 m asl and P. hartwegii reaches up to 4,200 m asl (Galicia et al. 2015). The Trans-Mexican volcanic belt functions as a zone of diversification, endemism, and biogeographic transition between the Nearctic and Neotropical regions. With very particular geomorphic, geological and biotic characteristics, the complexity of its ecosystems favors the evolution of a great variety of floral endemism. It possesses a large proportion of pine species and 21 out of 43 the endemic species of Mexico stands there (Rzedowski 1978, Fa & Morales 1991, Farjón et al. 1997, Villaseñor et al. 2007, Gernandt & Pérez-de la Rosa 2014, Halffter & Morrone 2017, Morrone 2019).
La Malinche National Park (hereinafter LMNP), a mountain located in the Trans-Mexican Volcanic Belt, is home to a wide diversity of pine species, including Pinus hartwegii Lindl., Pinus leiophylla Schiede ex Schltdl. & Cham., P. montezumae Lamb., P. pseudostrobus Lindl., and P. teocote Schiede ex Schltdl. Furthermore, reforestation programs have introduced P. strobiformis Engelm and P. patula Schiede ex Schltdl. & Cham. to the National Park. These species are distributed over the mountain, generally in monospecific fragmented stands or in clusters at elevations of between 2,000 and 4,000 m. P. leiophylla is characteristic of pine forests at elevations below 2,500 m, mostly as mixed forest. P. montezumae is the dominant species between 2,500 and 3,100 m and it may be found with P. pseudostrobus, P. teocote and various species of Abies and Quercus. At elevations of 3,500 to 4,200 m, there is a community of P. hartwegii (López-Domínguez & Acosta Pérez 2005, Rojas-García & Villers-Ruíz 2005, Suárez-Mota et al. 2013, Arriola Padilla et al. 2014). The most exploited species in the LMNP are pines, because the habitants of rural settlements use these species to obtain different resources, such as firewood, wood, bark, resin, and ocote (resinous wood chips), among others (Espejel 1996, Espejel Rodríguez et al. 2001).
Over the last 60 years, the combined effects of illegal and authorized commercial logging of conifers in the LMNP has led to a reduction of over 50 % of the 30,000 ha of forested land surface that existed in 1936. Currently, there are only 15,000 a of forest and 77 % of the vegetation exhibits certain degree of damage due to the impact of human activities, mainly logging, intentional and unintentional fires, and wood-stripping, a traditional practice known as "ocoteo” (Werner 1990, Espejel 1996, Espejel Rodríguez et al. 2001, Espejel Rodríguez et al. 2004). Ocoteo consists of manually cutting off strips of wood and bark impregnated with resin (Rzedowski 2006). A series of cuts are made in the trunk until a large wound with a depth of 5 to 10 cm is formed, more than two meters above ground level. Additionally, this practice often involves setting a fire at the base of the tree to increase resin flow from the wound (Veblen 1982). This traditional practice has been used in Mexico for decades, and it is an essential economic activity for many rural communities (Montero López et al. 2016, Molina-Luna & Arellanes Cancino 2016). The demand for "ocote" is due to its effectiveness as a fire starter in traditional stoves, ovens (Martínez 1992, Espejel 1996) and, more recently, tourist campfires. Records of this practice can be found since pre-Hispanic times in many lowland Mayan sites such as Tikal and Chinikiha in Chiapas and Belize (Montero López et al. 2016).
Although ocoteo is common in Mexico's pine forests, no study has evaluated the damage it causes, and it is uncertain if the persons which practices this activity have a preference for trees with specific dendrometric characteristics (DBH, height, and age). Some authors report that the purposes of ocoteo are to (1) obtain ocote for domestic use, (2) sell ocote in local markets, and (3) cause mechanical damage to trees to gradually weaken them until they die, in an attempt to find a "legal” way of simulating a natural death and extracting the natural resource (Werner 1990, Espejel Rodríguez et al. 2001). It is unknown to what extent the practice of ocoteo affects forests. However, this practice is deemed harmful and it represents an important factor in the loss of forest cover (Martínez 1992, Werner 1990, Rzedowski 2006).
Furthermore, studies have shown that when pine trees are mechanically damaged, traumatic ducts form in the wood as a reaction. There are fewer of these ducts in trees that have not been damaged (Nagy et al. 2004, Franceschi et al. 2005, Luchi et al. 2005, Krokene et al. 2008, Rodríguez-García et al. 2014).
The purpose of this study is to assess the damage done by ocoteo to different pine species found in the LMNP by making: (a) an estimation of ocoteo intensity at different altitudes according to the distribution of pine species and their associations on the mountain; (b) a comparison of dendrometric variables between wood-stripped (WS) and non-wood-stripped (NWS) trees; (c) a histological comparison of the presence and number of traumatic ducts in WS and NWS trees.
The research questions were: Which pine species are subjected to wood-stripping in the LMNP? What are the dendrometric characteristics of WS trees compared to NWS trees? Do WS trees present more traumatic ducts than NWS? Does the number of trees subjected to wood-stripping in the LMNP increase or decrease with altitude?
Materials and methods
Study area. The LMNP is located in the States of Tlaxcala and Puebla in the eastern central region of Mexico at 19º 06´ 30”- 19º 20´ 1” N and 97º 55´ 32”- 98º 09´ 55” W, at 2,300-4,461 m asl. It is one of the central volcanic landforms of the Trans-Mexican Volcanic Belt. It has an area of 46,093 ha, 33,161 ha of which are in the state of Tlaxcala and 12,932 ha in the State of Puebla. The Park contains about 37 % of the forested land in the State of Tlaxcala. Abies, Pinus, and Quercus are the dominant genera in forest stands in the lower parts of the LMNP. The climate is mild sub-humid, with summer rains, an average annual rainfall between 1,000 and 1,200 mm, and soils with variable composition derived from volcanic eruptions where lithosol, andosol, regosol, fluvisol, and luvisol predominate (Espejel-Rodríguez et al. 1999, López-Domínguez & Acosta Pérez 2005, Villers Ruiz et al. 2006).
Fieldwork. The LMNP was divided into three altitudes: low (2,000 to 2,500 m), mid (2,500 to 3,000 m), and high (3,000 to 3,500 m; Figure 1) with 33 randomly selected circular plots within the three altitudes. Each plot had an area of 400 m2 and a radius of 11.28 m (Schlegel et al. 2001, CONAFOR 2010). Sampling was done in 2017 and 2018. Each plot's center point was georeferenced, and all trees with a diameter ≥ 10 cm at breast-height (1.3 m) were counted (Schlegel et al. 2001, Vargas Larreta 2013, Kometter 2005). Tree height was measured using a Suunto clinometer (PM-5/360 degree/percentage) (Kometter 2005, Romahn de la Vega & Ramírez 2010).
To count the growth rings, cores were obtained from 193 trees located in the 33 plots of the LMNP; the extractions were made at a breast height of 1.3 m using a 5.15 mm Pressler increment borer (Haglöf Inc., Madison, MS, USA). After obtaining the cores, the drilled holes were sealed with Campeche wax (Rojas-García & Villers-Ruíz 2005, Gutiérrez & Ricker 2014). Then each core was placed in a plastic straw, marked, and taken to the laboratory for analysis.
From the 193 samples, four NWS individuals of each species across the different altitudes (a total of 28 samples) were sampled to study the traumatic ducts in the wood. Sample was taken at a breast height of 1.3 m, bark was removed and using a saw a piece of wood removed. In the case of WS trees, samples from the middle of the ocoteo wound were taken from 1-4 individuals of each species (20 samples) across the different altitudes (Rodríguez-García et al. 2014). All samples were fragmented to a size of 3-5 cm3 and were fixed in a formalin 10 ml (37 % formaldehyde)-glacial acetic acid 5 ml-ethanol (96 %) 50 ml- water 35 ml, mixture solution (Krommenhoek et al. 1986). The list of species is shown in Table 1.
Altitude | Number of plots | Number of pine trees | Total WS trees | WS species | Number of WS trees by species | Average of WS individuals by plot | Specific percentage of WS trees | Absolute percentage of WS trees |
---|---|---|---|---|---|---|---|---|
Low (2,000-2,500 m) |
7 | 79 | 8 | P. leiophylla | 8 | 1.3 ± 0.1 | 100 % | 10.1 % |
Medium (2,500-3,000 m) |
P. montezumae | 22 | 3.1 ± 0.9 | 81.4 % | 15.7 % | |||
14 | 171 | 27 | P. leiophylla | 2 | 1 ± 0 | 7.4 % | ||
P. pseudostrobus | 3 | 1 ± 0.4 | 11.1 % | |||||
High (3,000-3,500 m) | P. montezumae | 9 | 1.7 ± 0.2 | 75.0 % | 8.0 % | |||
12 | 149 | 12 | P. pseudostrobus | 9 | 1.7 ± 0.2 | 75.0 % | ||
P. teocote | 4 | 1 ± 0 | 8.3 % |
(WS) wood-stripped.
Laboratory work. Taxonomic keys available in the literature were used for taxonomical determination botanical species. Two samples of leaves and strobiles were taken from each WS tree (Martínez 1992, Villers Ruíz et al. 2006, Ramos-Dorantes et al. 2017). Samples were housed in the National Herbarium (MEXU).
Tree age was estimated by counting growth rings (Fritts 1976). For this process, the collected cores were dried in an oven (Ecoshel-HV-50, VelaQuin, México) at 23 °C for 48 hours following the criteria of Speer (2010) and Gutiérrez & Ricker (2014). Once dried, the cores were mounted on wooden bases and polished to discriminate between growth ring boundaries; afterwards, the identified rings were counted (Speer 2010).
Transverse sections of 25 µm were obtained using a sliding microtome (pfm AG-50996-Germany). Sections were stained with safranin and mounted in synthetic resin (Krommenhoek et al. 1986). The ducts were counted in 25 fields of 1 mm2 for each individual, without discriminating between growth rings
Analyses. For the evaluation of ocoteo, the dependent variables were intensity, percentage, and dendrometric characteristics (DBH, height, age, and the number of traumatic ducts) of the WS and NWS trees in different elevations.
The intensity of wood-stripping (ocoteo) on pine species at each altitude was estimated based on the relative percentage of WS trees per altitude, which was calculated with the following equation: Absolute percentage of WS trees = (Number of WS trees/Total number of registered trees) × 100. The following equation was used for wood-stripping intensity for different species in each elevation: Specific percentage of WS trees = (Number of WS trees of each species/Total number of WS trees) × 100.
The relative frequency of ocoteo in each species per altitude was calculated by: hi = fi/N, where hi = Relative frequency of the ith observation, fi= Absolute frequency of the ith observation, N = Total number of observations.
To determine whether the dendrometric characteristics of WS and NWS trees (DBH, height, and age) had any influence on ocoteo in the different elevations on the mountain, different statistical analyses were performed according to the type of the response variable. The DBH data was normalized with a logarithmic transformation and was compared using one-way ANOVA and Fisher's least-significant-difference method (statistically significant differences were considered at P < 0.05) followed by post-hoc Tukey-Kramer test. Height data were analyzed using a non-parametric Kruskal-Wallis test with the Dunnett test's respective multiple comparison tests (Zar 1999).
Tree age analysis was performed using a Generalized Linear Model (GLM) with Poisson error distribution and a logarithmic link function. The predictor variables were the altitude (low, medium, high), the tree species (P. montezumae, P. leiophylla, P. pseudostrobus, and P. teocote), and the presence or absence of wood-stripping (ocoteo) marks. Wood-stripping in WS pine species was evaluated using the same type of GLM. In these models, the resin duct predictive variables were: Species, Wood-stripping (ocoteo), and Individuals. The Species variable includes P. montezumae, P. leiophylla, P. pseudostrobus, and P. teocote. P. hartwegii was removed from the model because no WS specimens were found. The wood-stripping (ocoteo) variable includes both WS and NWS trees of the same species. The Individual variable (variation between individuals). Of the different models generated, the simplest or most parsimonious was selected so that non-significant explanatory variables were eliminated during the process. The statistical significance for each factor was determined by means of a deviance analysis and evaluation of how well the data fit the model.
Statistical analyses of the dendrometric variables (DBH and height) of WS and NWS trees were performed done using NCSS Statistical Program v10.0, NCSS 2015 (Hintze 2008). Analysis of the number of traumatic ducts and the species and ages of trees was done using R language v3.01, (R Core Team 2015, Crawley 2013).
Results
Intensity of wood-stripping (ocoteo). The intensity of the practice of ocoteo in the LMNP was higher in the mid altitude, with an absolute percentage of 15.7 %; values of 8.0 % and 10.1 % were recorded at high and low altitudes, respectively (Table 1), P. montezumae was the species most affected ocoteo at high and mid altitudes, no ocoteo was observed in P. hartwegii, only found at the highest altitudes. In the lowest altitudes, ocoteo was observed only by P. leiophylla trees (Table 1).
Comparison of dendrometric characteristics in WS and NWS trees. The DBH diameter of WS and NWS trees increased with altitude. The analysis of variance showed significant differences (F5,5 = 5.07, P < 0.001) in the DBH diameter of WS and NWS trees in high and mid altitudes (Figure 2A). In the LMNP, the tallest trees were those subjected to wood-stripping (Table 2). Significant differences were found between the heights of WS and NWS trees χ( = 51.6 and P < 0.001) in high and mid altitudes (Figure 2B). The average age of WS trees in mid and low altitudes was greater than NWS trees. At high altitudes WS trees were youngest (Table 2). Deviance analysis revealed no significant differences (Deviance = 0.14, df = 2, P > 0.05) between tree age among altitudes.
Altitude | DBH (cm) | Height (m) | Age (years) | ||||
---|---|---|---|---|---|---|---|
NWS | WS | NWS | WS | NWS | WS | ||
n | 277 | 47 | 277 | 47 | 147 | 47 | |
Low (2,000-2,500 m) |
26 ± 2 | 36 ± 6 | 12 ± 1 | 18 ± 3 | 48 ± 3 | 65 ± 7 | |
Medium (2,500-3,000 m) |
32 ± 1 | 46 ± 3 | 19 ± 0.6 | 26 ± 1 | 40 ± 2 | 54 ± 4 | |
High (3,000-3,500 m) |
27 ± 1 | 57 ± 5 | 14 ± 0.7 | 24 ± 2 | 48 ± 2 | 39 ± 6 |
(DBH) diameter at breast height; average values (N = 324) and standard error, (WS) wood-stripped, (NWS) non-wood-stripped.
Comparison of dendrometric characteristics between tree species subjected to wood-stripping (ocoteo). DBH. The WS species with the greatest diameter were P. pseudostrobus found at mid altitude, followed by P. montezumae at high altitude. The WS species with the smallest diameter was P. leiophylla (Table 3). No significant differences were found between species' DBH among altitudes (Figure 3A).
Altitude | Species | DBH (cm) | Height (m) | Age (years) | n |
---|---|---|---|---|---|
Low (2,000-2,500 m) |
Pinus leiophylla | 36 ± 6 | 19 ± 2 | 62 ± 6 | 8 |
Medium (2,500-3,000 m) |
Pinus pseudostrobus | 83 ± 10 | 25 ± 4 | 53 ± 10 | 3 |
Pinus montezumae | 42 ± 4 | 25 ± 1 | 49 ± 4 | 22 | |
Pinus leiophylla | 22 ± 13 | 28 ± 5 | 104 ± 12 | 2 | |
High (3,000-3,500 m) |
Pinus montezumae | 62 ± 6 | 23 ± 2 | 38 ± 6 | 9 |
Pinus teocote | 38 ± 18 | 28 ± 7 | 62 ± 18 | 2 | |
Pinus pseudostrobus | 45 ± 12 | 26 ± 5 | 35 ± 12 | 2 |
(DBH) diameter at breast height; average values (N = 48 trees); sample size (n) and standard error.
Height.- The tallest WS trees were P. leiophylla trees found at mid altitude. The smallest WS trees were P. leiophylla, but located at low altitude (Table 3), and Kruskal-Wallis showed no significant differences (Figure 3B).
Age.- The oldest trees subjected to wood-stripping (ocoteo) were P. leiophylla, found at mid and low altitudes. The youngest WS trees were located at high altitude and belong to P. pseudostrobus and P. montezumae (Table 3). Deviance analysis showed significant differences in the age of WS species (Deviance = 0.003, df = 3, P < 0.01), in the interaction of Elevation × Species (Deviance = 0.031, df = 3, P < 0.01) and in the interaction Species × Wood-stripping (ocoteo) (Deviance = 0.003, df = 3, P < 0.01).
Traumatic resin ducts.- Many more traumatic ducts were recorded in WS trees than NWS ones (Table 4). The highest number of traumatic ducts per mm2 were present in P. teocote followed by P. leiophylla, P. montezumae, and P. pseudostrobus (Table 4 and Figure 4A-I). Deviance analysis showed significant differences in the number of traumatic ducts between WS and NWS trees (Deviance = 8.96 e-13, df = 3, P < 0.001) and between species subjected to WS (Deviance < 2.2e -16, df = 3, P < 0.001).
Species | Number of traumatic resin ducts in Pinus species | |
---|---|---|
NWS | WS | |
Pinus montezumae | 19 ± 7 | 39 ± 7 |
Pinus leiophylla | 23 ± 5 | 51 ± 8 |
Pinus pseudostrobus | 12 ± 5 | 13 ± 2 |
Pinus teocote | 11 ± 0.9 | 98 ± 18 |
Pinus hartwegii | 16 ± 16 | *** |
*** The species was identified in the LMNP. However, there were no found WS trees.
Discussion
This study showed that the practice of wood-stripping in the LMNP mainly affects to P. leiophylla, P. montezumae, P. pseudostrobus, and P. teocote. No trees of P. hartwegii were found to have been subjected to wood-stripping. P. montezumae and P. leiophylla were the species with the highest number of trees showing evidence of wood-stripping, although in the latter species the practice was less intense. The results showed that the practice of wood-stripping in the LMNP is associated with the natural distribution of the Pinus species from 2,000 to 4,000 m; also with their abundance at different altitudes.
Rojas-García & Villers-Ruíz (2008) have reported that the most abundant species in the LMNP are P. montezumae and P. leiophylla, in association with P. pseudostrobus and P. teocote. Our study supported that report; the presence of P. montezumae at the mid altitude explains why the largest number of WS trees belong to this species. P. leiophylla is found below 2,900 m asl and predominantes at mid and low altitudes. According to López-Domínguez & Acosta Pérez (2005) P. leiophylla predominates below 2,500 m asl, which explains the results obtained showing that P. leiophylla was the only species subjected to wood-stripping at the low altitude.
According to Zamora-Campos et al. (2007), P. montezumae has wide natural distribution throughout the Mexican Republic. It is also the most chosen source of timber in the country, mainly in the Trans-Mexican Volcanic Belt. Wood-stripping is very important to the people who live in the vicinity of the volcano in the LMNP, which was evident in this study (Figure 5A-D).
The results of this study show that the practice of ocoteo in the LMNP is performed predominantly on trees with greater DBH (low 36 ± 6, mid 46 ± 3, high 57 ± 5), with greater height (low 18 ± 3, mid 26 ± 1, high 24 ± 2), and with respect to age only at low and mid altitudes (low 65 ± 7, mid 54 ± 4). These results suggest that the people who use this resource empirically select the most robust and tallest trees because they offer ocote strips of higher quality. According to people living in the vicinity, the amount of resin is higher in these trees (personal communication). Particularly, López Rucuch (1996) found a positive correlation between trunk diameter and the amount of resin produced by for P. oocarpa and P. pseudostrobus.
Other studies have positively related wound size to resin production in pine trees (Schweingruber 1996, Mumm & Hilker 2006, Granados-Sánchez et al. 2008, Bohlmann 2008, Hadiyane et al. 2015, Reyes-Ramos et al. 2016). The results obtained in this study showed that wood-stripping causes an increase in the number of traumatic ducts. The species with the highest number of traumatic ducts was P. teocote followed by P. leiophylla. The increase in the number of traumatic ducts is a natural response from Pinus trees. The number of traumatic ducts increases after the trees have gone through a period of stress caused by mechanical damage (Luchi et al. 2005, Franceschi et al. 2005, Moreira Tomé et al. 2008, Sampedro et al. 2011, Ferrenberg et al. 2014, Zas et al. 2014). According to Sampedro et al. (2011) the increase in the formation of traumatic ducts induces the use of a large amount of energy to drive defense systems and other functions associated with growth and reproduction.
The results obtained reveal that the trees selected for wood-stripping in the LMNP are more robust and taller individuals, while tree age varies according to altitude. Younger trees at a higher altitude and older trees at mid and low altitudes are selected for wood-stripping. The fact that age is not a tangible characteristic for the selection of individual trees for ocoteo may explain the lack of a clear relationship between tree age and the presence of wood-stripping marks.