ISSN e: 2007-4018 / ISSN print: 2007-4018

English | Español

     

 
 
 
 
 
 
 
 

Vol. XXII, issue 2 May - August 2016

ISSN: ppub: 2007-3828 epub: 2007-4018

Scientific article

Probability of mortality by fire damage of young Pinus hartwegii Lindl. trees in the Izta-Popo National Park

http://dx.doi.org/10.5154/r.rchscfa.2015.08.034

Robles-Gutiérrez1, César A. 1 ; Velázquez-Martínez1*, Alejandro 1 * ; Rodríguez-Trejo2, Dante A. 2 ; Reyes-Hernández1, Valentín J. 1 ; Etchevers-Barra1, Jorge D. 1

  • 1Colegio de Postgraduados, Campus Montecillo. Carretera México-Texcoco km 36.5. C. P. 56230. Montecillo, Texcoco, Estado de México Mexico.
  • 2Universidad Autónoma Chapingo, División de Ciencias Forestales. Carretera México-Texcoco km 38.5. C. P. 56230. Chapingo, Texcoco, Estado de México.Mexico

*Corresponding author. Email: alejvela@colpos.mx. Tel.: (595) 9520200 ext. 1470.

Received: August 04, 2015; Accepted: February 15, 2016

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

Abstract:

Great efforts are carried out at a global level in order to determine the effect of fire on the mortality of tree species, like the Pinus genus. This paper evaluated the influence of fire on the probability of mortality of young Pinus hartwegii trees in the Izta-Popo National Park, a protected natural area in central Mexico. The effects of a medium to high intensity fire that occurred in March 2013, a low intensity prescribed burning applied in April of the same year, and a control area without recent fire presence were evaluated and compared. The results showed that after 18 months of the applied treatments, the tree mortality in the area with prescribed burning was 13 %, whereas tree mortality in the area with fire and in the control area was 28 and 4 %, respectively. Six logistic models were obtained in order to predict mortality; the percentage of damaged canopy, scar height, and normal bark diameter were the significant predictive variables (P < 0.05). The results suggest that a decrease in the photosynthetic rate and nutrients,transport associated with damage to the canopy and cambium, are closely related to the probability of mortality.

Keywords:Prescribed burning; conflagration; percentage of damaged canopy; logistic model.

Introduction

In Mexico between 1998 and 2014, there were 8,500 fires on average with an estimated damage of 300,000 ha·year-1 (Comisión Nacional Forestal [CONAFOR], 2014). The greatest percentage of fires (98 %) was attributed to human activities, and the rest drived from natural causes such as electrical discharges and volcanic eruptions (Ressl & Cruz, 2012). Although it is apparent that human beings have favored the occurrence and frequency of fires due to their domestication, it is also evident that there have been modifications in fire regimes, to which some ecological systems have adapted (ecosystems with natural-anthropic regimes) (Pyne, 1995; Rodríguez, 2014). Therefore, and in accordance to fire response, ecosystems are classified as follows: 1) dependent, 2) sensible, 3) independent, and 4) influenced (Hardesty, Myers, & Fulks, 2005).

The fire dependent ecosystems require this ecological factor in order to survive in the landscape; in Mexico, these ecosystems are dominated by Pinus hartwegii Lindl. The ecosystems of this species present a natural-anthropic fire regime. The regime is characterized by high frequency (every five years on average), low intensity (Rodríguez, 2001; Rodríguez & Fulé, 2003) defined by speed, flame height, and longitude (Willson & Sorenson, 1979), and low to moderate severity (Jardel, Alvarado, Morfín, Castillo, & Flores, 2009) characterized by tree mortality, biomass consumption, and basal area and canopy reduction with the subsequent modification of stand composition and structure (Fulé & Covington, 1997; Glitzenstein, Streng & Wade, 2003; Regelbrugge & Conrad, 1993).

Pinus hartwegii is a species with adaptations to fire such as: thick bark, grass stage, basal sprouts, recovery of burned canopy, natural prune, and capacity to be regenerated in burned areas (Rodríguez & Fulé, 2003). However, poor fire management strategies as fire suppression, as well as the large human contribution to these disturbances, have resulted in the high frequency or lack of fire management. Such causes have modified the fire regime in P. hartwegii ecosystems, which are associated with high mortality indices of young trees (Rodríguez, Martínez, & Ortega, 2004). It has been indicated that the resistance of a tree to fire-caused mortality depends primarily on the useful morphological characteristics to protect its vital tissues (which increase with age), the ability to recover from injuries, and the season in which the events occur (Regelbrugge & Conard, 1993). In the case of annual species, fires are less severe before the opening of the buds, mainly due to the existing carbohydrate reserve in the trees (Van Wagner, 1973).

The importance of predicting tree mortality resides in that it is a planning tool for carrying out operations such as tree regeneration or prescribed burnings (Regelbrugge & Conard, 1993). When modeling the mortality or survival of trees, the dependent variable is a binomial character (i. e., a tree is alive or dead after a period of time); therefore, the logistic regression is the most utilized mathematical expression (Hosmer & Lemeshow, 2000). The tree mortality models must include variables that reflect the severity of the fire damages, in this case, proportion of crown damage is one of the most used characteristics (Hood, McHugh, Ryan, Reinhardt, & Smith, 2007; Van Wagner, 1973). Likewise, the inclusion of the cambium damage reflects the increase in the probability of mortality for a given crown damage level (Hood et al., 2007; Peterson, 1983).

This paper seeks to contribute to the understanding of fire as a natural process that frequently operates as an integral part of the ecosystem where it occurs and therefore, fire can be utilized as a tool for the integral management of the dependent ecological systems. The objective was to evaluate and model the P. hartwegii trees probability of mortality in seedling and sapling stage development, under a low intensity prescribed burning and to a medium to high fire intensity, in a forest under conservation conditions. The main hypothesis considers that the individual probability of mortality is greater when there is a larger proportion of crown damage, greater fire scar height, and lower bark diameter.

Materials and methods

Area of study

The area of study is located in the Iztaccíhuatl-Popocatépetl National Park, found in the eastern central part of the Trans-Mexican Volcanic Belt in central Mexico. The site is known as Telapón hill and is found between the extreme coordinates 19º 22’ 8.4” -19° 22’ 44” LN and 98º 42’ 32.40” - 98º 43’ 1.2” LO, at an average altitude of 3,800 m (Figure 1). The climate is temperate, sub-humid, with rain in the summer and an annual average temperature of 15 °C (García, 1988). The soil is Andosol with a sandy-loam texture, a dark brown to black color, and with 2 to 8 % organic matter (Secretaría de Medio Ambiente y Recursos Naturales [SEMARNAT], 2013). The dominant vegetation is the P. hartwegii forest and among the most common species in the undergrowth are Lupinus montanus Kunth, Muhlenbergia quadridentata (Kunth) Trin., Festuca spp., Calamagrostis spp. and Acaena elongate L. (Rzedowski, 1978).

Figure 1. Location of the area of study for the probability of mortality evaluation of Pinus hartwegii in the Izta- Popo National Park.

Treatments

On March 3, 2013, a fire occurred affecting 59 ha of the study area. Subsequently, on March 15 of the same year, preliminary observations were made of the percentage of trees cown damage, which served as the basis to choice the area of evaluation. The main assumption was that with a greater percentage of crown damaged in the trees of this area, a better contrast would be obtained in comparison with prescribed burning, and that by utilizing the trees of the two zones in the survival modeling, a higher response interval would be obtained for different crown damage ratios. On the same date, the appropriate sites were determined to carry out the prescribed burning and for the establishment of the control area; the location of the forest roads was a priority. A clinometer (Suunto® PM5/1520, Finland) and a compass were utilized for the homogenization of the areas with regard to gradient and exposure. The preselected locations were georeferenced using a GPS (Garmin GPSMAP® 78, USA) with a maximum precision of 3 m for their subsequent location in the map. The points located in the field were verified with the program ArcGIS 9.3® (Environmental Systems Research Institute [ESRI], 2009) and then proceeded to the selection of the previously standardized areas. A digital elevation model was utilized for the area of study at a 1:50,000 scale (Instituto Nacional de Estadística y Geografía [INEGI], 2013), with which gradient (30 to 40 %), exposure (E-SE), and altitude (3,780 to 3,900 m) maps were obtained.

Finally, according to previous obtained values, the optimal surfaces for the delimitation of the burned area, the establishment of the prescribed burning and the control area were determined using the algebraic map technique (De Mers, 2002). Figure 2 shows the three selected areas.

Figure 2. Delimitation of the areas of study for the probability of mortality evaluation of Pinus hartwegii in the Izta-Popo National Park. The numbered points represent the sampling sites.

Four ha were delimited of the total burned area, The preliminary evaluation showed that in this area wind and slope favors the fire, with an average flame height of 2.5 m. The flame height was calculated through random sampling from the scar height of 50 trees in the evaluated area. The prescribed burning was carried out on April 13, 2013, with drip torches as a mean of ignition and backing strips as the burning technique. The wind speed at medium flame was 10 km·h-1, average relative humidity of 34 %, average temperature of 13 °C, and fire speed was 5 m·min-1. The activity finished when wind gusts reached 15 km·h-1, the relative humidity decreased to 29 %, the temperature increased to 16 °C, and the flame height reached 1 m, resulting in a treated surface area of 1 ha. The burning was done based on the standard established in the NOM- 015-SEMARNAT/SAGARPA-2007 (SEMARNAT, 2007) and with the collaboration of specialized fire workers (Universidad Autónoma Chapingo, Comisión Nacional Forestal (CONAFOR) and the Comisión Nacional de Áreas Naturales Protegidas (CONANP) crews), as well as two community crews (Figure 3). The flame height for the prescribed burning and wild fire were the fire intensity indicators (Willson & Sorenson, 1979). The control area remained delimited to 3.5 ha.

Figure 3. Prescribed burning for the probability of mortality evaluation of Pinus hartwegii: a) start of the prescribed burning (ignition of belt one), b) fire behavior, c) ignition of belt two, d) experimental site moments after the prescribed burning.

Evaluation of the treatments

In the control and fire areas, a systematic sampling was carried out. Due to its dimensions, the prescribed burning was evaluated by simple random sampling. In both cases, the site was the sampling unit. The systematic arrangement was done with a 50 m separation between lines and sites. 15 sites were established in the fire area, 6 in the prescribed burning area, and 10 in the control (Figure 2). The sampling sites were 100 m2 square plots, where the trees with a height between 0.30 to 3.0 m were measured.

After the treatments were applied, the following dendrometric data was obtained from each sampling site: basal diameter (Db, cm), normal bark diameter at a height of 1.30 m (DNCC, cm), total height (At, m), scar height (ACtz, m), canopy longitude (Lco, m), longitude of the crown damaged (Lcd, m), and proportion of crown damage (PCD, %). This last variable was calculated with the equation proposed by González-Rodríguez and Rodríguez-Trejo (2004):

Finally, the trees of the quadrants were evaluated 18 months after the application of prescribed burning or fire occurrence for the mortality analysis after two growth periods of primary meristems. All the trees with all the aerial part or part of it alive, or with the presence of basal sprouts were classified as “alive” and were assigned a value of “0” (zero), whereas the trees that did not show these characteristics were classified as “dead” and were assigned a value of “1”.

The registered information in the field was refined through At dispersion graphs from the DNCC and Db, thus detecting potential erroneous data. The slimness index (TI) was calculated through its average and confidence interval, as well as with the dispersion graph, trees data was verified and those that had an allometric relation outside the norm were eliminated (i.e., a tree that visually appeared as an erroneous data and that did not fit in the confidence interval of the TI); thus, a proportional relation was not presented between the DNCC or Db and its height. The three areas (fire, prescribed burning, and control) were refined by classifying the trees by its height. The categories were seedlings (trees between 0.30 to 1.29 m in height) and saplings (trees between 1.3 to 6.0 m in height). The averages of the dendrometric variables per height category are shown in Table 1.

Table 1. Comparative table of the dendrometric variables obtained after three treatments applied in order to determine the probability of mortality of Pinus hartwegii.

Treatments / Tratamientos Monte bravo (1.3-6.0 m) / Seedlings (1.3-6.0 m) Brinzales (0.30-1.29 m) / Saplings (0.30-1.29 m)
Density (trees·ha-1) / Densidad (árboles·ha-1) DNCC (cm) At (m) ACtz (m) PCD (%) Densidad (árboles·ha-1) / Density (trees·ha-1) Db (cm) At (m) ACtz (m) PCD
Fire / Incendio 809 (156.40) 5.09 (0.28) 3.00 (0.12) 1.84 (0.07) 81.71 (2.16) 1673 (525.70) 2.26 (0.06) 0.69 (0.02) 0.68 (0.02) 99.87 (0.12)
Prescribed burning / Quema prescrita 900 (287.50) 6.09 (0.47) 3.12 (0.18) 0.80 (0.03) 48.93 (5.04) 780 (215.40) 3.40 (0.17) 0.93 (0.04) 0.46 (0.03) 97.33 (0.94)
Testigo / Control 920 (193.10) 5.63 (0.37) 3.16 (0.14) 0 0 910 (250.50) 2.87 (0.14) 0.74 (0.03) 0 0
The numbers in parenthesis correspond to the standard error. DNCC = Normal bark diameter; Db = Basal diameter; At = Total height; ACtz = Scar height; PCD = Percentage of damaged canopy.

Statistical analysis

The mortality proportion was determined for the three treatments in order to have an impact estimator of the trees. Subsequently, the individual probability of tree mortality less than 6.0 m in height with some type of fire damage was determined through the adjustment of a logistic model:

The effects of the independent variables were determined through adjusted models for each treatment and as a group, using the PROC LOGISTIC of the SAS® statistical packet version 9.4 (Statistical Analysis System [SAS], 2013). The independent variables analyzed were Db, DNCC (only for stem exclusion), At, PCD, and ACtz. With the significance of parameters, parsimonious models were obtained (i.e., the simplest models possible with the greatest explicative or predictive ability) that included only significant variables (Ӽ2 P < 0.05) by treatment and by group. Furthermore, it was determined that the selected mathematical expressions did not include the unit inside the Wald confidence interval. Finally, a percentage of accordance for the predicted probabilities was obtained with regard to those observed.

The models were classified according to the quantity of treatments they contained and with the height category (seedling and sapling). Thus, the “total” models included in their structure all of the trees of their respective height category, both from the prescribed burning as well as the fire, whereas the models for burning and fire were separately considered as “partial”. Furthermore, a general model was obtained that included the trees from 0.30 m up to 6.0 m in height for both treatments.

Results and discussion

In general, low mortality indices were shown for P. hartwegii in the three treatments. Mean mortality was 28 % in the fire area, 13 % in the prescribed burning, and 3 % in the control.

The low mortality in the fire area is attributed to a short residence time at high or lethal temperatures fires a characteristic related to fire adaptation P. hartwegii (Rodríguez et al., 2004), and to the temporality in which the fire occurred. Vera and Rodríguez (2007) evaluated the same species in the Ajusco zone in the Distrito Federal and reported lower mortality indices for high intensity fires that occurred in the March (13.9 %) compared to those that occurred on May (52 %) two years after treatments application.

The high mortality presented in the prescribed burning probably was due to the application timing (April). Vera and Rodríguez (2007) reported a mortality of 67 % for a high intensity prescribed treatment (that emulated wild fire), 4.4 % for a low intensity burning, and 4.2 % for the control area in treatments established in March. The differences between the study by Vera and Rodríguez and this paper could be attributed to the variability in the fire behavior and to the temporality in which this factor presented itself (Ryan & Reinhardt, 1988; Vera & Rodríguez, 2007). Vera and Rodríguez (2007) also report that the low intensity prescribed burnings in the month of May can reach up to 14.5 % mortality, due to the accumulated dryness during the fire season (evaluation carried out two years after applying the treatments). Other influencing factors are the tree vigor before the prescribed burning (Ryan & Reinhardt, 1988) and the seasonality of bud growth (Wagener, 1961). Both of them, can be important and defining factors, given that in the study area, dwarf mistletoe (Arceuthobium sp.) is present; furthermore, it was visually confirmed that the beginning of bud elongation coincided with the application of prescribed burning, exposing the sapling buds to the fire effects.

The probability of mortality models were adjusted from a final basis of 549 trees, obtaining six models in accordance with the parsimony and significance of parameters (Ӽ2 P < 0.05). None of the generated models include the unit inside the Wald confidence interval and, in general, all of them showed good concordance between the percentage of the predicted probabilities and the observed probabilities. It can be seen in Table 2 that the minimum concordance values are close to 60 % and that the models structured from PCD showed greater concordance to the prediction with regard to the appearance of the mortality event. The most consistent explicative variables in the models of the sapling category (1.3 to 6.0 m) were PCD and DNCC, although this latter variable was only consistent in the case of the total modeling. In the case of the seedlings, the probability of mortality was explained in a greater percentage due to the scar height and, finally, for the general model, the most explicative variables were both the PCD and the scar height (Table 2).

Table 2. Significance of parameters for the logistic mortality models of Pinus hartwegii.

Model Verisimilitude radius (P-value) Concordance of predicted and observed probabilities (%) Parameters Estimator Standard error Ӽ2 P-value Wald confidence intervals (95 %)
Inferior Superior
Less than 1.30m 0.0192 59.5 Intercept -1.926 0.409 <0.0001
ACtz 1.285 0.550 0.0195 1.229 10.631
Greater than 1.29 m <0.0001 77.7 Intercept -9.915 2.305 <0.0001
PCD 0.100 0.024 <0.0001 1.054 1.159
Greater than 1.29 m Total 0.0376 59.6 Intercept -0.496 0.405 0.2212
DNCC -0.146 0.074 0.0473 0.747 0.998
Greater than 1.29 m <0.0001 75.2 Intercept -12.205 3.218 0.0001
Conflagration PCD 0.124 0.033 0.0002 1.061 1.210
Greater than 1.29 m 0.0001 81.8 Intercept -6.787 2.610 0.0093
Burning PCD 0.065 0.028 0.0208 1.010 1.128
General <0.0001 72.2 Intercept -11.548 2.589 <0.0001
ACtz 0.871 0.230 0.0002 1.523 3.753
PCD 0.100 0.025 <0.0001 1.052 1.161
ACtz = Scar height; PCD = Percentage of damaged canopy; DNCC = Normal bark diameter.

Figure 4 shows the mortality behavior of P. hartwegii in each of the generated logistic models. Based on the models, the maximum probability of mortality of seedlings with a scar height of 1.29 m is 44 % (Figure 4a). In the sapling category, the DNCC and PCD were the most important explicative variables. According to this, the maximum individual probability of mortality estimated with the total model on the individuals with DNCC of 1 cm was 13 % (Figure 4b). Regarding the PCD, the young trees showed a maximum probability of mortality of 52 % in the total model, whereas the maximum values in the prescribed burning and conflagration were 43 and 57 %, respectively (Figure 4c). Finally, the general model was determined by the PCD and ACtz, variables that favored the probability of mortality as the values increased (Figure 4d).

Figure 4. Individual probability of mortality of the Pinus hartwegii trees in function of: a) scar height in saplings, b) normal bark diameter in stem exclusion, c) damaged canopy by treatment and total in stem exclusion, and d) damaged canopy and scar height for all trees.

The probabilities of mortality of P. hartwegii obtained with the generated models correspond to those reported by other authors. In general, the probability of mortality increases when the measure of the trees is smaller, in this case trees less than 1.3 m in height, due to the damages caused to the cambium (Peterson & Ryan, 1986), independent of the quantity of the damaged canopy. This is basically due to the younger trees having a tendency to have thinner bark (Costa, Oliveira, Viexas, & Neto, 1991; Regelbrugge & Conard, 1993). Wyant, Omi, and Laven (1986) determined that the damage to the cambium is correlated to the ACtz. In this study, the ACtz was the characteristic that showed the best results in predicting mortality in saplings; as this variable increases, the individual probability of mortality is greater (Figure 4a).

In the case of trees greater than 1.29 m in height, mortality is influenced mainly by the PCD and DNCC (Table 2). The DNCC-mortality relation is notably inverse (Figure 4b), which coincides with the results of Vera and Rodríguez (2007). The PCD is one of the indicators most utilized in determining mortality (Hood et al., 2007; Peterson, 1985; Ryan, Peterson, & Reinhart, 1988). The models generated in our study show that at a greater PCD there is a greater probability of mortality, which coincides to that established for Pinus ponderosa Dougl. ex Laws (Wyant et al., 1986) and P. lambertiana Dougl. (Mutch & Pearson, 1998). The mortality of these species is influenced mainly by the substantial crown reduction (Wyant et al., 1986) and, therefore, a decrease in the production of photosynthates. Figure 4c shows that in the sapling category, the model predicts the least probability of mortality in the area with prescribed burning when there is 100 % damaged canopy, in comparison with fire or the total model. Vera and Rodríguez (2007) also found that P. hartwegii mortality in prescribed burnings in comparison to forest fires is much less during fire season. The inflexion points (around 70 % damaged canopy) of the curves defined by the PCD models coincide with those generated by Stephens and Finney (2002) for P. ponderosa. This suggests that the deficit in the production of photosynthates behaves in an exponential way after the loss of three fourths of the canopy. The general model shows the combined effect of PCD and damage to the cambium, represented by the ACtz (Figure 4d); it is observed that tree mortality increases when damage to the canopy is accompanied by damage to the conductive tissue.

Conclusions

The probability of mortality increases with a greater percentage of crown damage and thinner bark diameter in trees classified as saplings; in the case of the seedlings, the probability increases with a greater scar height. Therefore, the probability of mortality in trees of 0.3 to 6.0 m increases by the interaction of a greater scar height and percentage of damaged canopy. The mortality was less in prescribed burning than in the fire area, which is attributed to a greater fire intensity in the later. The prescribed burning presented 54 % less mortality than the fire area, but four times high mortality than the control area, due mainly to the tree vigor before the burning and to the timing treatment. Prescribed burnings are a viable tool in the management of fire dependent ecosystems (i.e. P. hartwegii); however, it is recommended that future studies evaluate the characteristics of the behavior and temporality of fire, as well as the physiology of the trees, with the purpose of providing more elements to the forest managers such as when and how to implement this tool and the possible responses of the species to the stress of fire.

References

Comisión Nacional Forestal (CONAFOR). (2014). Reporte semanal de resultados de incendios forestales 2014: Del 01 de enero al 31 de diciembre del 2014. Jalisco, México. Comisión Nacional Forestal.

Costa, J. J., Oliveira, L. A., Viexas, D. X., & Neto, L. P. (1991). On the temperature distribution inside a tree under fire conditions. International Journal Wildland Fire, 1(2), 87-96. doi: 10.1071/WF9910087

De Mers, M. N. (2002). GIS modeling in raster. New York, USA: John Wiley & Sons.

Enviromental Systems Research Institute (ESRI). (2009). ArcGIS (versión 9.3). Software de procesamiento digital de imágenes satelitales. EUA: Autor.

Fulé, P. Z., & Convington, W. W. (1997). Fire regimes and forest structure in the Sierra Madre Occidental, Durango, México. Conservation Biology, 13(3), 640-652. Obtenido de http://www.redalyc.org/articulo.oa?id=57404107

García, E. (1988). Modificaciones al sistema de clasificación climática de Köppen. México: UNAM.

Glitzenstein, J. F., Streng, D. R., & Wade, D. D. (2003). Fire frequency effects on longleaf pine (Pinus palustris P. Miller) vegetation in South Carolina and Northeast Florida, USA. Natural Areas Journal, 23, 22-37. Obtenido de http://www.srs.fs.usda.gov/pubs/ja/ja_glitzenstein001.pdf?

González-Rodríguez, A., & Rodríguez-Trejo, D. A. (2004). Efecto del chamuscado de copa en el crecimiento en diámetro de Pinus hartwegii Lindl. en el Distrito Federal, México. Agrociencia, 38, 537-544. Obtenido de http://www.redalyc.org/articulo.oa?id=30238508

Hardesty, J., Myers, R. L., & Fulks, W. (2005). Fire, ecosystems, and people: A preliminary assessment of fire as a global conservation issue. The George Wright Forum, 22, 78-87. Obtenido de http://www.georgewright.org/224hardesty.pdf

Hood, S. M., McHugh, C. W., Ryan, k. C., Reinhardt, E., & Smith, S. L. (2007). Evaluation of a post-fire tree mortality model for western USA conifers. International Journal of Wildland Fire, 16, 679-689. doi: 10.1071/WF06122

Hosmer, D. W., & Lemeshow, S. (2000). Applied logistic regression. New York, USA: Wiley.

Instituto Nacional de Estadística y Geografía (INEGI). (2013). Modelo digital de elevación, E14B31. 1:50000. Aguascalientes, México: Autor. Obtenido de http://www3.inegi.org.mx/sistemas/biblioteca/ficha.aspx?upc=702825733124

Jardel, P. E. J., Alvarado, J. E., Morfín, R. E., Castillo, N. F., & Flores, G. J. G. (2009). Regímenes de fuego en ecosistemas forestales de México. In G. J. G. Flores (Ed.), Impacto ambiental de incendios forestales (pp. 73-100). México: Instituto Nacional de Investigaciones Forestales Agrícolas y Pecuarias-Colegio de Postgraduados.

Mutch, L. S., & Parsons, D. J. (1998). Mixed conifer forest mortality and establishment before and after prescribed fire in Sequoia National Park, California. Forest Science, 44(3), 341-355. Obtenido de http://www.nps.gov/seki/learn/nature/upload/fsdoc.pdf

Peterson, D. L. (1983). Predicting fire-caused mortality in four northern Rocky Mountain conifers. New forests for a changing world: Proceedings of the 1983 convention of the Society of American Foresters (pp. 276-280). Maryland, USA: Society of American Foresters.

Peterson, D. L. (1985). Crown scorch volume and scorch height: Estimates of post-fire tree condition. Canadian Journal of Forest Research, 15, 596-598. doi: 10.1139/X06-315

Peterson, D. L., & Ryan, K. C. (1986). Modeling postfire conifer mortality for long- range planning. Environmental Management, 10(6), 797-808. doi: 10.1007/BF01867732

Pyne, S. J. (1995). World fire: The culture of fire on earth. USA: University of Washington Press.

Regelbrugge, J. C., & Conard, S. G.(1993). Modeling tree mortality following wildfire in Pinus ponderosa forests in the Central Sierra Nevada of California. International Journal Wildland Fire, 3(3), 139-148. doi: 10.1071/WF9930139

Ressl, R., & Cruz, I. (2012). Detección y monitoreo de incendios forestales mediante imágenes de satélite. CONABIO. Biodiversitas, 100, 12-13. Obtenido de http://www.biodiversidad.gob.mx/Biodiversitas/Articulos/biodiv100art7.pdf

Rodríguez, T. D. A. (2001). La ecología del fuego en el ecosistema de Pinus hartwegii. Revista Chapingo Serie Ciencias Forestales y del Ambiente, 7(2), 145-151. Obtenido en http://www.chapingo.mx/revistas/forestales/contenido.php?id_revista_numero=24

Rodríguez, T. D. A. (2014). Incendios de vegetación. Su ecología manejo e historia. México: Biblioteca básica de agricultura.

Rodríguez, T. D. A., & Fulé, P. Z. (2003). Fire ecology of mexican pines and a fire management proposal. International Journal Wildland Fire, 12, 23-37. doi: 10.1071/WF02040

Rodríguez, T. D. A., Martínez, H. H. C., & Ortega, B. V. (2004). Ecología del fuego en bosques de Pinus hartwegii. In R. L. Villers, & B. J. López (Eds.), Incendios forestales en México: Métodos de evaluación (pp. 103-120). México: UNAM .

Ryan, K. C., & Reinhardt, E. D. (1988). Predicting postfire mortality of seven western conifers. Canadian Journal of Forest Research, 18, 1291-1297. doi: 10.1139/x88-199

Ryan, K. C., Peterson D. L., & Reinhart, E. D. (1988). Modeling longterm fire caused mortality of douglas-fir. Forest Science, 34, 190-199. Obtenido en http://www.fs.fed.us/rm/pubs_other/rmrs_1988_ryan_k001.pdf

Rzedowski, J. (1978). Vegetación de México. México: Limusa.

Secretaría de Medio Ambiente y Recursos Naturales (SEMARNAT). (2009). Norma Oficial Mexicana NOM- 015-SEMARNAT/SAGARPA-2007 que establece las especificaciones técnicas de métodos de uso del fuego en los terrenos forestales y en los terrenos de uso agropecuario. México: Diario Oficial de la Federación. Obtenido de http://www.profepa.gob.mx/innovaportal/file/3331/1/nom-015-semarnat-sagarpa-2007.pdf

Secretaría de Medio Ambiente y Recursos Naturales (SEMARNAT). (2013). Programa de manejo del Parque Nacional Iztaccíhuatl-Popocatépetl. México: Autor. Obtenido de http://iztapopo.conanp.gob.mx/documentos/programa_de-manejo_izta_popo.pdf

Statistical Analysis System (SAS Institute Inc.). (2013). SAS/ STAT 9.4 user’s guide. Cary, NC, USA: Author

Stephens, S. L., & Finney, M. A. (2002). Prescribed fire mortality of Sierra Nevada mixed conifer tree species: Effects of crown damage and forest floor combustion. Forest Ecology and Management, 162(2002), 261-271. doi: 10.1016/S0378-1127(01)00521-7

Van Wagner, C. E. (1973). Height of crown scorch in forest fires. Canadian Journal of Forest Research, 3, 373-378. doi: 10.1139/X73-055

Vera, V. V., & Rodríguez, T. D. A. (2007). Supervivencia e incremento en altura de Pinus hartwegii a dos años de quemas prescritas e incendios experimentales. Agrociencia, 41, 219-230. Obtenido de http://www.colpos.mx/agrocien/Bimestral/2007/feb-mar/art-9.pdf

Wagener, W. W. (1961). Guidelines for estimating the survival of fire-damaged trees in California. Berkeley, CA, USA: Pacific Southwest Forest and Range Experiment Station- U. S. Forest Service.

Willson, C., & Sorenson, J. (1979). Algunos factores comunes acerca del comportamiento del fuego, en caso de incendios forestales trágicos y casi trágicos. EUA: USDA Forest Service.

Wyant, J. G., Omi, P. N., & Laven, R. D. (1986). Fire induced tree mortality in a Colorado ponderosa pine/douglas-fir stand. Forest Science, 32(1), 49-59.

Figures:

Figure 1. Location of the area of study for the probability of mortality evaluation of Pinus hartwegii in the Izta- Popo National Park.
Figure 2. Delimitation of the areas of study for the probability of mortality evaluation of Pinus hartwegii in the Izta-Popo National Park. The numbered points represent the sampling sites.
Figure 3. Prescribed burning for the probability of mortality evaluation of Pinus hartwegii: a) start of the prescribed burning (ignition of belt one), b) fire behavior, c) ignition of belt two, d) experimental site moments after the prescribed burning.
Figure 4. Individual probability of mortality of the Pinus hartwegii trees in function of: a) scar height in saplings, b) normal bark diameter in stem exclusion, c) damaged canopy by treatment and total in stem exclusion, and d) damaged canopy and scar height for all trees.

Tables:

Table 1. Comparative table of the dendrometric variables obtained after three treatments applied in order to determine the probability of mortality of Pinus hartwegii.
Treatments / Tratamientos Monte bravo (1.3-6.0 m) / Seedlings (1.3-6.0 m) Brinzales (0.30-1.29 m) / Saplings (0.30-1.29 m)
Density (trees·ha-1) / Densidad (árboles·ha-1) DNCC (cm) At (m) ACtz (m) PCD (%) Densidad (árboles·ha-1) / Density (trees·ha-1) Db (cm) At (m) ACtz (m) PCD
Fire / Incendio 809 (156.40) 5.09 (0.28) 3.00 (0.12) 1.84 (0.07) 81.71 (2.16) 1673 (525.70) 2.26 (0.06) 0.69 (0.02) 0.68 (0.02) 99.87 (0.12)
Prescribed burning / Quema prescrita 900 (287.50) 6.09 (0.47) 3.12 (0.18) 0.80 (0.03) 48.93 (5.04) 780 (215.40) 3.40 (0.17) 0.93 (0.04) 0.46 (0.03) 97.33 (0.94)
Testigo / Control 920 (193.10) 5.63 (0.37) 3.16 (0.14) 0 0 910 (250.50) 2.87 (0.14) 0.74 (0.03) 0 0
The numbers in parenthesis correspond to the standard error. DNCC = Normal bark diameter; Db = Basal diameter; At = Total height; ACtz = Scar height; PCD = Percentage of damaged canopy.
Table 2. Significance of parameters for the logistic mortality models of Pinus hartwegii.
Model Verisimilitude radius (P-value) Concordance of predicted and observed probabilities (%) Parameters Estimator Standard error Ӽ2 P-value Wald confidence intervals (95 %)
Inferior Superior
Less than 1.30m 0.0192 59.5 Intercept -1.926 0.409 <0.0001
ACtz 1.285 0.550 0.0195 1.229 10.631
Greater than 1.29 m <0.0001 77.7 Intercept -9.915 2.305 <0.0001
PCD 0.100 0.024 <0.0001 1.054 1.159
Greater than 1.29 m Total 0.0376 59.6 Intercept -0.496 0.405 0.2212
DNCC -0.146 0.074 0.0473 0.747 0.998
Greater than 1.29 m <0.0001 75.2 Intercept -12.205 3.218 0.0001
Conflagration PCD 0.124 0.033 0.0002 1.061 1.210
Greater than 1.29 m 0.0001 81.8 Intercept -6.787 2.610 0.0093
Burning PCD 0.065 0.028 0.0208 1.010 1.128
General <0.0001 72.2 Intercept -11.548 2.589 <0.0001
ACtz 0.871 0.230 0.0002 1.523 3.753
PCD 0.100 0.025 <0.0001 1.052 1.161
ACtz = Scar height; PCD = Percentage of damaged canopy; DNCC = Normal bark diameter.