Servicios Personalizados
Revista
Articulo
Indicadores
-
Citado por SciELO
-
Accesos
Links relacionados
-
Similares en SciELO
Compartir
Revista mexicana de ciencias forestales
versión impresa ISSN 2007-1132
Rev. mex. de cienc. forestales vol.8 no.43 México sep./oct. 2017
Article
Sampling unit for determining the spatial variability of the surface area burnt by forest fires
1Instituto Nacional de Investigaciones Agrícolas y Pecuarias
The strategies against forest fires are directed to priority areas, cartographically defined according to such criteria as fire hazard (inception, propagation, difficulty of control and impact). This is reflected, among other things, on the size of the burnt surface. Accordingly, the burnt surface can be considered as a support criterion for defining priority areas for protection against forest fires. However, the processes to determine their spatial variation must be standardized; otherwise, there may be a mismatch of their estimates between different zones. For this reason, among others, a common sampling unit must be defined in order to obtain comparable and compatible information. Thus, the objective of the present work was to determine a sample unit size that enables the capture of the spatial variability of the surface burnt by forest fires. Reference sites (RSs) were analyzed, considering: a) four sampling intensities (100, 300, 500 and 1 000 sites); and b) twelve RS sizes (1, 2, 4, 8, 10, 15, 30, 50, 70, 100, 150 and 200 Km2). The average burnt area in each RS was determined based on fire information for the state of Jalisco during the 2005-2015 period. The variability of the average burnt area diminished as the size of the RS increased, until an asymptotic behavior occurred. Considering that the data did not define a normal distribution, the Kruskal-Wallis test determined that there is a significant difference between RS sizes. An analysis between paired independent samples (Mann-Whitney test) defined that there is a difference between the sizes of neighboring RSs in the area where the asymptote begins. Finally, an RS of 100 Km2 was defined as the size of the common sampling unit to be used to generate a cartography of the spatial variability of the area burnt by forest fires.
Keywords: Priority areas; fuel loads; forest fire management; forest fire hazard; forest fire risk; reference places
Las estrategias contra incendios forestales se dirigen a áreas prioritarias, determinadas cartográficamente con base en criterios como el peligro de incendio (inicio, propagación, dificultad de control e impacto). Lo anterior se refleja, entre otras cosas, en la superficie que llega a quemarse. De acuerdo a esto, puede considerarse la superficie quemada como un criterio de soporte para definir áreas prioritarias de protección contra incendios forestales. Sin embargo, los procesos para conocer su variación espacial deben estar estandarizados, de lo contario puede implicar una falta de correspondencia de sus estimaciones entre diferentes zonas. Para ello, entre otros aspectos, se tiene que establecer una unidad de muestreo común, para obtener información comparable y compatible. En ese contexto, el objetivo del presente trabajo fue determinar un tamaño de unidad de muestreo que permita captar la variabilidad espacial de la superficie quemada por incendios forestales. Se analizaron sitios de referencia (SIR), y consideraron: a) cuatro intensidades de muestreo (100, 300, 500 y 1 000 sitios); y b) doce tamaños de SIR (1, 2, 4, 8, 10, 15, 30, 50, 70, 100, 150 y 200 Km2). Se calculó la superficie quemada promedio en cada SIR, con base a partir de la información de incendios del período 2005-2015, en el estado de Jalisco. La variabilidad de la superficie promedio incendiada disminuyó a medida que aumentaba el tamaño del SIR, hasta alcanzar un comportamiento asintótico. Dado que los datos no presentaron una distribución normal, la prueba Kruskal-Wallis demostró que existe diferencia significativa entre los tamaños de SIR. Un análisis entre pares de muestras independientes (prueba Mann-Whitney), evidenció diferencias entre los tamaños de SIR vecinos en el área de inicio de la asíntota. Finalmente, del SIR de 100 km2 resultó el tamaño de la unidad de muestreo común, para la generación de cartografía de la variabilidad espacial de la superficie quemada por incendios forestales.
Palabras clave: Áreas prioritarias; carga de combustible; manejo del fuego; peligro de incendio forestal; riesgo de incendio; sitios de referencia
Introduction
Every year, the surface area affected by forest fires across the world and the number of these events are very variable, depending on the weather and the social and economic conditions existing in the damaged areas. In this regard, the average number of forest fires per year in Mexico is 8 000; these fires manifest in various magnitudes. However, their average size is 30 ha (Conafor, 2013). Consequently, the impact that they cause due to their frequency and intensity is an important ecological disturbance factor, as they bring about transformations at various levels in several forest ecosystems (Jardel et al., 2010).
Likewise, they also have social and economic effects. Therefore, fire prevention and firefighting must be contemplated within a fire management plan that includes, among other aspects, the zoning of priority areas against forest fires (Conafor, 2010). This would contribute to an efficient use of human, material and financial resources, all of which are limited (Carrillo et al., 2012).
In general, three criteria are analyzed in order to determine priority areas of protection against forest fires: fire risks, fire hazard and potential damage or value (Flores et al., 2016a); these are assessed and weighted in terms of a series of variables (Rojo et al., 2001). Thus, the location and sizing of surfaces based on certain levels of hazard support the definition of strategies for both fire prevention and firefighting.
Few studies have been carried out to estimate forest fire hazards (Carrillo et al., 2012). For example, the spatial variation of fire hazard is modeled based on interpolations (kriging) of fuel loads (Rodríguez et al., 2011). Furthermore, variables that exert an influence on the hazard are assessed through multicriterion analysis in order to determine their corresponding model (Muñoz et al., 2005). Another assessment method is based on teledetection (images taken by the MODIS sensor) and meteorological indices, correlated with the temporal variation of the fuel’s moisture content (Yebra et al., 2005).
Nevertheless, we must take into account the fact that, in order to determine the forest fire hazard, variables affecting their inception, propagation, difficulty of control and impact need to be assessed (Rodríguez et al., 2011).
This is reflected, among other things, on the eventually burnt surface, since it implies that: a) there were conditions for their inception; b) the rapidity of their propagation determines the surface area burnt during a given time period; c) in general, the larger the surface, the greater the difficulty to control the fires, and d) the effect of a high impact can be exponential in direct relationship to the affected surface. Therefore, the burnt surface area is a support criterion for determining priority areas of protection against forest fires. However, a sampling strategy must first be established in which, for instance, an adequate site size is determined for the assessment of the burnt surface.
From the point of view of forestry, this site size is determined for various purposes (biomass estimation, tree density, diversity, etc.); a minimum size is sought that involves savings in terms of both time and costs, as well as the guarantee of statistically accurate estimates. This requires the site to display a tendency to reduce the variance of estimations, taking into account that small sampling sites favor a high variation in the estimations, and that as the size increases the variation decreases (Aguirre et al., 1997). Nevertheless, this tendency can eventually determine an asymptotic behavior in which the increase in the size of the site does not involve a relevant reduction of the variance. Therefore, the adequate site size is where the said asymptotic behavior begins (Gallegos et al., 2006).
Since the definition of a common site size may imply a lack of correspondence between the estimations of the burnt surface areas resulting from fires in different areas, the purpose of this work was to define one that will allow capturing the spatial variation of this surface. Thus, the following hypothesis is proposal: the variability of the burnt surface diminishes as the size of the site increases to a point in which this variability determines an asymptotic behavior.
Materials and Methods
Description of the study area
The present research uses information from the state of Jalisco, located in western central Mexico (Figure 1), with a surface area of 78 588 Km2, of which 68 % has a warm subhumid climate along the coast and its central area, while 18 % is subhumid temperate, in the higher parts of the mountain chains, and 14 % is dry and semi-dry, in the north and northeast of the state. The mean annual precipitation is approximately 850 mm, although in the coastal zones it is over 1 000 mm per year. According to Rzedowski’s classification (1986), Jalisco has 13 types of vegetation, most of which consists of pine oak forests and deciduous tropical forests (Ramos et al., 2007).
![](/img/revistas/remcf/v8n43//2007-1132-remcf-8-43-00117-gf5.jpg)
Figure 1 Study area of the sampling unit for determining the spatial variability of the surface area burnt by forest fires in the state of Jalisco, Mexico.
In Jalisco, the average annual burnt surface ranges between 15 000 and 20 761.58 hectares (Conafor, 2015; Semadet, 2016). The temperate forest is the most affected, with an annual average of nearly 12 008 hectares; it is followed by the lowland forest and undetermined vegetation, with an average of 1 903.9 hectares; grasslands, with 729.57 hectares; tropical forests, with 463.32 hectares; shrubs, with 190.55 hectares, and montane cloud forests, with 24.27 hectares (Conafor, 2016).
Burnt surface
In general, fire hazard is defined as the probability that an event may propagate, prosper and cause damages in the vegetation, as a result of various factors, such as temperature, relative humidity, characteristics of the fuels, land conditions, and speed and direction of the wind (Chandler et al., 1983). Based on the above, the hazard is reflected, among other aspects, on the surface area that burns during a forest fire (Vilchis et al., 2015). In general, a fire may be considered to be more dangerous in direct proportion to the size of the burnt surface. Therefore, in this research we considered the analysis of the average damaged surface within a given area. However, the burnt surface is not synonymous with the forest fire hazard.
Surface area of reference sites
The most convenient minimum sampling surface area for capturing the variability of the average burnt surface was defined in the actual project. For this purpose, a series of circular polygons of twelve different surfaces known as reference sites (RSs) ‒‒1, 2, 4, 8, 10, 15, 30, 50, 70, 100, 150 and 200 Km2‒‒ were analyzed. It was decided that the polygons should be circular because this shape has been used in various forest samplings. Furthermore, circular polygons require only one control point, which facilitates the distribution of the sites (McRoberts et al., 1992). Any geometric figure can be used, as what matters is the size of the site, not its shape.
In order to capture the variability of the average burnt surface, which may occur due to the density of sampling points, four sampling intensities were established: 100, 300, 500 and 1 000. Their location was selected at random, and therefore any point had the same probability of being selected; furthermore, this probability is independent of other points.
Although the spatial independence of the data is outside the scope of this study, it must be analyzed, as through it we seek a low heterogeneity of the variable burnt surface that will ensure the non-existence of spatial autocorrelation. Otherwise, more complex statistic methods integrating the modeling of the spatial structure of the data would have to be applied (Zas et al., 2008). The twelve circular polygons were subsequently located in each reference site (RS).
The next step was to determine how many forest fires have taken place in each RS (Figure 2). For this purpose, statistic information on the forest fires of the 2005-2015 period (Conafor, 2015) was utilized to determine the average burnt surface corresponding to each RS.
Variability analysis
As stated above, for each of the four sampling intensities, the average burnt surface area per RS was recorded in order to estimate the descriptive statistics by RS and by sampling intensity. The analysis of the variability of the burnt surface in relation to the size of the RS was carried out by charting the corresponding variation coefficients by sampling intensity. The most adequate surface of the sampling polygon (RS) was subsequently determined, using as a criterion the threshold in which the variability (variation coefficient) began to display an asymptotic behavior.
In order to ensure that the selected RS size differed from the adjoining RSs, the kurtosis values and the asymmetry coefficient were analyzed. This made it possible to determine that the data do not tend to define a normal distribution. The Kruskal-Wallis (non-parametric) test was applied (Fowler et al., 1998) in order to determine whether or not there was a significant difference between RS sizes. Finally, an analysis between independent paired samples was carried out using the Mann-Whitney test (Fowler et al., 1998) to determine whether or not there was a significant difference between the sizes of adjacent RSs.
Thematic cartography
Once the most adequate RS size was obtained, the corresponding thematic map was generated, illustrating the spatial variation between the surface burnt by forest fires across the state of Jalisco. The interpolation technique known as “weighted inverse distance” --in which the value of an unsampled point is the mean weighted inverse distance between the values of the sampled points-- was used for this purpose (Burrough and McDonnell, 1998). Its interpolation based on the weighted inverse distance is represented by means of the following linear function (Flores, 2001):
Where:
* (x
0
) = Estimated value at an unsampled site
x 0 = Location referred to a coordinate system
* (x
i
) = Observed value at a sampled site
x i, d i = Distances between each of the sample sites and the unsampled point
p = Exponent of the distance (weighting)
n = Number of sampled sites. In this work, a weighting value (p) of 2 was used, as it best represented the spatial variation of the average burnt surface.
Results
General statistics
Based on the four sampling intensities (100, 300, 500 and 1 000 sites), the statistics corresponding to the various site sizes were assessed (Table 1). The minimum value, regardless of the size of the site, was observed to be zero burnt hectares, while the maximum burnt surface ranged between 104.9 has (in 1 km2) and 866.66 ha (in 200 km2). In terms of the means and the modes, it was inferred that at most sites there were no burnt surfaces.
Sampling unit
The measure of description was the variation coefficient as a dispersion measure, as it describes the amount of variability (in relation to the mean) without being based on the sample size. Thus, unlike the standard deviation, it was possible to compare the dispersion of the utilized sampling intensities, regardless of the difference between their means. The variability of the average burnt surface decreased in direct proportion to the RS size (Figure 3). This tendency was adjusted to the correlation models corresponding to the sampling intensities (Table 2). This occurred in every case, until an asymptote was reached in which the values of the variation coefficient tended to become stable, which occurs within an interval approximately between 70 and 100 km2. Thus, the RS size finally considered as the breaking point of the variability, was 100 km2 (radius = 5.650 m), defined as the sampling unit (SU) for the estimation of the spatial variation of the average surface burnt by forest fires.
Table 2 Statistic values to which the tendency of the variation coefficient is adjusted in relation to the site size, for different sampling intensities.
![](/img/revistas/remcf/v8n43//2007-1132-remcf-8-43-00117-gt6.jpg)
The Kruskal-Wallis test, a non-parametric equivalent of ANOVA, was utilized to determine the existence of significant differences between the twelve RS sizes for each sampling intensity. Unlike ANOVA, instead of means, the intervals of the RSs are analyzed. First, the data for all RS sizes were combined; subsequently, they were ordered by size, from smaller to larger. The average interval for the data for each RS size was then calculated. The decision criterion applied to accept or discard the null hypothesis that there is no significant difference is the probability value (P). When this is below 0.05, the null hypothesis is discarded.
The resulting probability values are shown in Table 3; at all the intensities, the value of P is below 0.05, which indicates that there is a statistically significant difference between the various RS sizes.
The comparative analysis of paired samples using the Mann-Whitney test determined that there are no differences in the burnt surfaces between RS sizes above and below 100 km2 (Table 4), which supports their selection as the breaking point for the asymptotic behavior of the variability (variation coefficient).
![](/img/revistas/remcf/v8n43//2007-1132-remcf-8-43-00117-gf7.jpg)
Coeficiente de variación = Variation coefficient; Tamaño de sitio = Size site. The vertical line marks the beginning of the asymptotic behavior of the variation coefficient (breaking point [100 km2]).
Figure 3 Tendency of the variation coefficient in relation to site size for different sampling intensities: (A) 100, (B) 300, (C) 500, (D) 1 000.
Table 3 Results of the Kruskal-Wallis test.
RS | MR-100 | MR-300 | MR-500 | MR-1000 |
---|---|---|---|---|
1 | 373.92 | 1 155.5 | 2023.58 | 4 325.38 |
2 | 380.83 | 1 244.12 | 2 108.3 | 4 556.05 |
4 | 408.88 | 1 304.37 | 2287.86 | 4910.26 |
8 | 478 215.00 | 1 457.30 | 2 461. 56 | 5 279.23 |
10 | 505 035.00 | 1 496.85 | 2553.21 | 5 433.18 |
15 | 551 005.00 | 1 661.44 | 2 801.95 | 5 762.49 |
30 | 634 455.00 | 1 860.19 | 3186.03 | 6457.73 |
50 | 697 695.00 | 2 048.03 | 3 435.15 | 6 844.41 |
70 | 739915.00 | 2189.47 | 3605.55 | 7 012.20 |
100 | 745 075.00 | 2 296.35 | 3 776.34 | 7 130.64 |
150 | 821 955.00 | 2 415.04 | 3 862.75 | 7204.18 |
200 | 869.02 | 2 477.35 | 3 903.74 | 7 090.27 |
H | 341.91 | 903.539 | 1 235.66 | 1 528.94 |
P | 0.000 | 0.000 | 0.000 | 0.000 |
MR = Interval mean (the number indicates the simple size); H= Kruskal-Wallis statistic; P= Probability value.
Spatial distribution of the burnt surface
Figure 4 represents the results of this process, which shows a larger burnt surface area (in 100 km2) in the areas where the density of forest fires is higher. Conversely, there were regions with a low or null burnt surface. At the same time, we should note that the density of the number of fires and the spatial variation of the burnt surface area was observed. The highest concentration of fires occurred in central Jalisco, where the highest concentrations of population are located.
Discussion
Although several researches define the hazard of fire forests (Muñoz et al., 2005; Yebra et al., 2005; Rodríguez et al., 2011), the information (basically in thematic cartography) is generated using various procedures. Notably, the absence of definition of common strategies limits the compatibility and comparability of the results. In order to avoid this, a statistically robust process must be established to determine the burnt surface. For this purposes, it should begin by defining a sampling unit; the present work does so successfully.
Although there are no records of similar researches, it is important to point out that some of the studies specify certain sampling units (Carrillo et al., 2012). Nevertheless, these are not defined statistically but by circumstantial situations. For example, in the work by Carrillo et al., (2012) addressing the fire hazard criterion, a logistic regression analysis is carried out to establish the most important variables through the random sampling of 10 km2 units. However, this unit size is determined only by the spatial resolution of the source information utilized and therefore cannot be used as a comparative sampling unit.
Although there are several studies in which a sampling unit is determined for purposes of forest inventories (O'Regan and Arvanitis, 1966; Zeide, 1980; McRoberts et al., 1992; Aguirre et al., 1997), there are no bases to use this unit in a sampling design to estimate the spatial variability of the surface burnt by forest fires so as to ensure a better location for and sizing of fire management activities. For, taking into account that the task of fighting forest fires is very costly and dangerous, such a unit is very helpful for the management and allocation of resources (Díaz et al., 1998; Rojo et al., 2001).
A sampling unit allowing the estimation of the spatial variation of the burnt surface in a standardized manner was set. This makes it possible to carry out comparisons between several regions in order to determine priority attention areas. Likewise, the methodology may be used for other types of comparative analyses, e.g. of the temporary variations in the frequency and location of forest fires (Rodríguez, 2012).
The size of the site is defined based on the tendency to the decrease in the variation of the variability (coefficient of variation) of the average burnt surface as the size of RSs increases until it reaches an asymptotic behavior (Gallegos et al., 2006), which occurs in a similar way in the various tested sampling intensities. Thus, an area of 100 km2 (radius = 5.650 m) was established as sampling unit (SU) for the estimation of the surface burnt by forest fires. Yet, according to the results, what was determined was not a specific breaking point but an interval (between 70 and 100 km2). In order to standardize a methodological process and ensure that the SU is above the lower limit, an area of 100 km2 was selected.
In the process of defining the sampling unit, official statistics for a 10-year period are used (Conafor, 2013). Thus, the georeferentiation of forest fires makes it possible to locate those areas with the highest incidence. These evidences show that, every year, forest fires do not occur in a totally random manner but follow a self-organized critical behavior (Torres et al., 2007; Ávila et al., 2010; Pérez et al., 2013). However, a more formal analysis is required to determine an occurrence pattern (Ávila et al., 2010). It may be nevertheless assumed that the average burnt surface does not vary greatly from one year to another. Still, predominant factors such as the variation in the climate conditions and the accumulation of forest fuels must be taken into account. The information on the burnt surface defined using the heat points’ technology can also be incorporated, though with certain reservations, as there is not always a direct relationship between these points and the burnt surfaces (Tansey et al., 2008).
Since the highest frequency of the estimated average burnt surfaces corresponded to low values (<30 ha), and a lower frequency, to 400 ha, the distribution is not normal for all the tested sampling intensities. Notably, a comparative analysis of paired samples established the existence of a significant difference between site sizes within the interval which determines the beginning of the asymptotic behavior of the variability (variation coefficient). Thus, we corroborate that the variability is better captured with the 100 km2 site size than with 70 km2. This defines a breaking point that is closer to the asymptotic behavior of the variability.
As for the spatial variation of the burnt surface, in locations with a higher density of forest fires, the defined SU was observed to capture a larger number of fires affecting the average burnt surface. The resulting cartography indicates that the largest surfaces are located mainly in central Jalisco. Two important regions to the south and southwest of that state have also been identified. In all of these areas there are large concentrations of population and an important surface is devoted to agricultural and livestock breeding activities; this is relevant, if we consider that most forest fires are caused by anthropic actions. This, in turn, determines geographic patterns of occurrence of fires, based on which priority areas can be delimited in relation to the surface burnt by forest fires (Flores et al., 2016a).
The information thus generated can be utilized both to reinforce and to validate the current determinations of priority hazard areas, which are used for medium- and long-term planning purposes, based on the integral analysis of various parameters, including the temperature, exposure, slope, behavior of the fire, effect of the fire, etc. (Flores et al., 2016b).
Nevertheless, it should be noted that the notion of forest fire hazard must not be mistaken for that of hazard index, which is used for short-term (daily) estimations for operational purposes.
Certain authors suggest defining priority forest fire areas in terms of the spatial variation of their density, using specific processes such as the kernel density function (de la Riva et al., 2004; Kuter et al., 2011), which is the basis for the estimation of the probability of a high or low density.
Conclusions
The reduction of the variability (variation coefficient) of the average burnt surface observed when the RS size is increased until it displays an asymptotic behavior suggests the acceptance of the initially posed hypothesis. Hence, the sampling unit for the estimation of the variability of the surface burnt by forest fires is 100 km2.
The determination of a common sampling unit in order to estimate the spatial variability of the surface burnt by forest fires prevents the lack of correspondence between estimations for different areas.
The results for several areas can be not only compatible but also comparable. I.e., based on the sampling unit (SU), the average burnt surface of a specific area to be sampled can be located, and comparisons can be carried out between several areas. Likewise, the SU makes it possible to share information between different areas, allowing the integration of information in order to determine fire control and prevention strategies in the state of Jalisco or at a regional level.
It is possible to use the determined sampling unit to support the determination of standardized strategies for the validation of forest fire hazard areas, which may be defined using other parameters, such as the temperature, the exposure, the slope, the behavior of the fire, the effect of the fire, etc.
The non-parametric analysis of the data evidences the existence of statistically significant differences between the various site sizes, allowing the determination of the breaking point for the beginning of the asymptotic behavior of the variability, as well as to the establishment of an adequate site size.
As for the spatial variation of the burnt surface, we conclude that the interpolation technique, based on the weighted inverse distance, defines an adequate distribution, as there are spatial coincidences between the highest values for the burnt surface and those areas where the largest concentrations of fire occur.
Although the sample unit was defined using an example at a state level (in Jalisco, Mexico), it can be used in areas with similar conditions. Otherwise, specific sampling units will have to be generated, at a regional, state or municipal level, using the methodology described in this paper. It is suggested to work with statistical data on forest fires for a minimum period of five years.
Referencias
Aguirre C., O. A., J. Jiménez P., E. J. Treviño G. y B. Meraz A. 1997. Evaluación de diversos tamaños de sitio de muestreo en inventarios forestales. Madera y Bosques 3 (1):71-79. [ Links ]
Ávila F., D. Y., M. Pompa G. y E. Vargas P. 2010. Análisis espacial de la ocurrencia de incendios forestales en el estado de Durango. Revista Chapingo. Serie Ciencias Forestales y del Ambiente 16 (2): 253-260. [ Links ]
Burrough, P. A. and R. McDonnell. 1998. Principles of Geographical Information Systems. Oxford University Press. New York, NY USA. 333 p. [ Links ]
Carrillo G., R. L. et al. 2012. Análisis espacial de peligro de incendios forestales en Puebla, México. Interciencia 37 (9): 678-683. [ Links ]
Chandler, C., P. Cheney, P. Thomas, L. Trabaud and D. Williams. 1983. Fire in Forestry. Vol. II. Wiley. New York, NY USA. 298 p. [ Links ]
Comisión Nacional Forestal (Conafor). 2010. Procedimiento para la elaboración de un mapa de áreas de atención prioritaria contra incendios forestales. Gerencia de Protección Contra Incendios Forestales. Comisión Nacional Forestal. Semarnat. México, D.F., México. 49 p. [ Links ]
Comisión Nacional Forestal (Conafor). 2013. Incendios Forestales. Campaña 2013. http://www.conafor.gob.mx:8080/documentos/docs/7/4339Campa%C3%B1a%20de%20contra%20incendios%202013.pdf (30 de mayo de 2016). [ Links ]
Comisión Nacional Forestal (Conafor). 2015. Incendios. http://www.conafor.gob.mx/web/temas-forestales/incendios/ (4 de agosto de 2015). [ Links ]
Comisión Nacional Forestal (Conafor). 2016. Reporte semanal de resultados de incendios forestales 2016. Programa Nacional de Prevención de Incendios Forestales. Gerencia de Protección Contra Incendios Forestales. Comisión Nacional Forestal. Semarnat. México, D.F., México. 22 p. [ Links ]
de la Riva J., F. Pérez C., N Lana R. and N Koutsias. 2004. Mapping wildfire occurrence at regional scale. Remote Sensing of Environment 92 (3): 363-369. [ Links ]
Díaz D., R., R. Salvador, J. Valeriano y X. Pons. 1998. Detección de superficies forestales en Cataluña mediante imágenes de satélite durante el período 1975-1995. Aplicación para la caracterización del régimen de incendios y los procesos de regeneración de la vegetación. Serie Geográfica 7: 129-138. [ Links ]
Flores, J. G. 2001. Modeling the spatial variability of forest fuels arrays. Doctoral dissertation. Department of Forest Sciences. Colorado State University. Fort Collins, CO, USA. 184 p. [ Links ]
Flores G., J. G. et al. 2016a. Descripción de variables para definición de riesgo de incendios forestales en México. INIFAP-CIRPAC. Campo Experimental Centro Altos de Jalisco. Tepatitlán de Morelos, Jal., México. Folleto técnico Núm. 161 p. [ Links ]
Flores G., J. G. et al. 2016b. Descripción de variables para la definición de peligro de incendios forestales en México. Folleto técnico Núm. 3. INIFAP-CIRPAC. Campo Experimental Centro Altos de Jalisco. Tepatitlán de Morelos, Jal., México. 53 p. [ Links ]
Fowler, J., L. Cohen and P. Jarvis. 1998. Practical statistics for field biology. 2nd edition. Wiley. England, United Kingdom. 259 p. [ Links ]
Gallegos, A., O. A. Aguirre y G. A. González. 2006. Optimización de inventarios forestales para manejo forestal de un bosque tropical en Jalisco, México. In: IV Simposio Internacional sobre Manejo Sostenible de los Recursos Forestales. 19 al 22 de abril 2006. Pinar del Río, Cuba. pp. 171-184. [ Links ]
Jardel P., E. J. et al. 2010. Prioridades de investigación en manejo del fuego en México. CIECO, Fondo Mexicano para la Conservación de la Naturaleza, UAS, UNAM, Universidad de Guadalajara, USAID. México, D.F., México. 41 p. [ Links ]
Kuter N., F. Yenilmez and S. Kuter. 2011. Forest fire risk mapping by kernel density estimation. Croatian Journal of Forest Engineering 32 (2): 599-609. [ Links ]
McRoberts, R. E., E. O. Tomppo y R. L. Czaplewski. 1992. Diseños de muestreo de las evaluaciones forestales nacionales. Antología de conocimiento para la evaluación de los recursos forestales nacionales. FAO. Roma, Italia. 21 p. [ Links ]
Muñoz R., C. A. et al. 2005. Desarrollo de un modelo espacial para la evaluación del peligro de incendios forestales en la Sierra Madre Oriental de México. Investigaciones geográficas 56: 101-117. [ Links ]
O'Regan, W. G. and L. G. Arvanitis. 1966. Cost efectiveness in forest sampling. Forest Science 12 (4): 406-414. [ Links ]
Pérez V., G., M. A. Márquez L., A. Cortés O. y M. Salmerón M. 2013. Análisis espacio-temporal de la ocurrencia de incendios forestales en Durango, México. Madera y Bosques 19 (2): 37-58. [ Links ]
Ramos V., I., S. Guerrero V. y F. M. Huerta M. 2007. Patrones de distribución geográfica de los mamíferos de Jalisco, México. Revista mexicana de biodiversidad 78 (1): 175-189. [ Links ]
Rodríguez M., A. 2012. Cartografía multitemporal de quemas e incendios forestales en Bolivia: Detección y validación post-incendio. Ecología en Bolivia 47 (1): 53-71. [ Links ]
Rodríguez T., D. A. et al. 2011. Modelaje del peligro de incendio forestal en las zonas afectadas por el huracán Dean. Agrociencia 45 (5): 593-608. [ Links ]
Rojo M., G. E., J. Santillán P., H. Ramírez M. y B. Arteaga M. 2001. Propuesta para determinar índices de peligro de incendios forestal en bosque de clima templado en México. Revista Chapingo Serie Ciencias Forestales y del Ambiente 7 (1): 39-48. [ Links ]
Rzedowski, J. 1986. Vegetación de México. Limusa. México, D.F., México. 432 p. [ Links ]
Secretaria de Medio Ambiente y Desarrollo Territorial (Semadet). 2016. Incendios Forestales. Secretaria del Medio Ambiente y Desarrollo Territorial Secretaria del Medio Ambiente y Desarrollo Territorial http://incendios.semadet.jalisco.gob.mx/estadisticas (6 de junio de 2016). [ Links ]
Tansey, K. et al. 2008. Relationship between MODIS fire hot spot count and burned area in a degraded tropical peat swamp forest in Central Kalimantan, Indonesia. Journal Geophysical Research 113: 1-8. [ Links ]
Torres R., J. M., O. S. Magaña T. y G. A. Ramírez F. 2007. Índice de peligro de incendios forestales de largo plazo. Agrociencia 41 (6): 663-674. [ Links ]
Vilchis F., A. Y. et al. 2015. Modelado espacial para peligro de incendios forestales con predicción diaria en la cuenca del río Balsas. Agrociencia 49 (7): 803-820. [ Links ]
Yebra A., M., A. de Santis y E. Chuvieco. 2005. Estimación del peligro de incendios a partir de teledetección y variables meteorológicas: variación temporal del contenido de humedad del combustible. Recursos Rurais 1 (1): 9-19. [ Links ]
Zas A., R., P. Martíns G. y R. de la Mata P. 2008. Autocorrelación espacial: un problema común…mente olvidado. Cuadernos de la Sociedad Española de Ciencias Forestales 24: 139-145. [ Links ]
Zeide, B. 1980. Plot size optimization. Forest Science 26 (2): 251-257. [ Links ]
Received: May 04, 2017; Accepted: August 31, 2017