Introducción

Los incendios forestales son uno de los factores más importantes que dan forma y modifican el ecosistema y la biodiversidad del Planeta desde hace millones de años. En múltiples regiones del mundo, estos eventos son causados por una combinación de fenómenos naturales como los rayos y la disponibilidad de combustible vegetal debido a las sequías pronunciadas en ciertos meses del año. Por ejemplo, en México 40 % de los tipos de vegetación son mantenidos por incendios (^{Rodríguez, 2014}), y en algunos bosques tropicales ocurren en cada temporada de secas, por lo que las especies arbóreas exhiben rasgos de adaptación a los incendios, y estos desempeñan un papel ecológico (^{Nasi et al., 2002}). En gran medida, los seres humanos son responsables, directa o indirectamente, de manera deliberada o por descuido, de entre 75 % y 96 % de ellos (^{Hirschberger, 2016}).

La expresión del fuego como factor ecológico en los ecosistemas forestales ocurre a partir de regímenes de fuego, es decir, los patrones de incendios repetidos y expresados como frecuencia, estación, tipo, gravedad y extensión en un paisaje (^{Scott et al., 2014}). Se estima que 61 % de las ecorregiones del Planeta tienen regímenes de fuego degradados o muy degradados, de modo que el fuego desempeña un papel fundamental en el mantenimiento de muchos ecosistemas; debido a las actividades antrópicas, el fuego se comporta, actualmente, de manera diferente que en cualquier otra época de la historia de la humanidad (^{Shlisky et al., 2007}).

A una escala global, la alteración de los regímenes de fuego es una importante fuente de emisiones de gases de efecto invernadero. En el ámbito regional, los incendios impactan las reservas de biomasa, el ciclo hidrológico y la salud de las personas, además afectan significativamente la biodiversidad en los bosques. ^{Nasi et al. (2002)} señalan que a finales del siglo XX, los cambios en la dinámica entre el hombre, el fuego y el aumento en la frecuencia de El Niño-Oscilación del Sur, han originado una situación en la que los incendios constituyen una importante amenaza para muchos bosques y su biodiversidad. En particular, El Niño-Oscilación del Sur ha generado condiciones inusualmente cálidas y secas en muchas regiones del mundo susceptibles a los incendios, lo que puede resultar en un incremento de la superficie quemada (^{Burton et al., 2020}).

Desde el punto de vista de la salud de las personas, el humo de los incendios forestales tiene el potencial de afectar a millones de personas, por lo que constituye un problema de salud pública considerable (^{Chen et al., 2006}; ^{Dennekamp y Abramson, 2011}).

Con los avances tecnológicos actuales es posible monitorear grandes superficies mediante el uso extensivo de los datos satelitales. Una aplicación es la detección de los incendios activos (focos de calor) que, generalmente, se basa en la región espectral infrarroja media (MIR) (3-5 μm), la cual implica una radiancia espectral de la temperatura a la que se quema la vegetación (de 500 a 1 000 K). Esta es superior a las temperaturas promedio de la Tierra (de 300 K), por lo que es posible distinguir los incendios activos (^{Chuvieco, 2008}). El sensor más importante es el Espectrorradiómetro de Imágenes de Resolución Moderada (MODIS) de los satélites Aqua y Terra, el cual identifica los “píxeles de fuego” del tamaño de 1 km² (en promedio). ^{Müller et al. (2013)} refieren que el tamaño de los incendios detectables depende de la temperatura del fuego, del área del incendio, de la cobertura vegetal y del ángulo de visión del sensor.

La observación y el registro de los focos de calor también permite estudiar, analizar y responder preguntas científicas a partir de datos espacio-temporales. La información disponible a la fecha no respalda un aumento general de la superficie quemada, ni de la gravedad global de los incendios, pero hay evidencias, basadas en algunas escalas regionales, de que se ha incrementado el número de eventos y el área total afectada (^{Doerr y Santín, 2016}). Por lo tanto, cabe suponer que los incendios forestales se pueden analizar y modelar a nivel regional.

En México hay algunos esfuerzos de investigación encaminados a comprender diversos aspectos de los incendios forestales, desde los puntos de vista de la causalidad (Ávila-Flores et al., 2010; ^{Antonio y Ellis, 2015}; ^{Pompa-García et al., 2018}; ^{Zúñiga-Vásquez y Pompa-García, 2019}) hasta de la predicción (^{Perez-Verdin et al., 2014}; ^{Ibarra-Montoya y Huerta-Martínez, 2016}; ^{Galván y Magaña, 2020}; ^{Monjarás-Vega et al., 2020}; ^{Ruíz-García et al., 2022}).

Oaxaca es el estado con mayor biodiversidad en todo México; alberga casi la mitad de las especies de plantas y de tipos de vegetación del país, así como 40 % de los taxa de mamíferos, 63 % de los de aves, 26 % de los reptiles y 23 % de los peces de río (^{Oviedo, 2002}). El estado se ubica en el suroeste de México y colinda con Puebla, Veracruz, Chiapas, Guerrero, y con el océano Pacífico al sur (Figura 1). El bosque seco tropical, en el que viven como endémicas hasta 70 % de las especies de plantas, cruza el estado.

Figura 1 Ubicación geográfica del estado de Oaxaca en México. 

Con base en los resultados obtenidos a partir de datos de los años 2000 a 2012, Oaxaca es una de las primeras cinco entidades con más emisiones de carbono negro (4 557-6 309 t⋅año^-1) y de carbono orgánico (48 441-70 663 t⋅ año^-1) derivados de incendios forestales (^{Cruz et al., 2014}). En el último informe técnico de la Comisión Nacional Forestal (Conafor), en 2019 Oaxaca figuraba entre los 10 estados con mayor número de incendios forestales y se hallaba entre los primeros cinco con más superficie dañada (^{Conafor, 2020}).

^{Doerr y Santín (2016)} señalan la necesidad de realizar estudios a escalas local y regional, y dada la gran biodiversidad existente en Oaxaca, en el presente estudio se analizan retrospectivamente los eventos de focos de calor ocurridos en los últimos 19 años (2001-2019) en el estado. El objetivo principal fue analizar y caracterizar el patrón espacio-temporal de los casos de focos de calor para conocer si sus patrones varían en el espacio y en el tiempo, o si los eventos se dispersan al azar en ambas escalas. Se consideró que en el área de estudio los focos de calor no están dispersos al azar en el espacio, sino que se forman agrupamientos, los cuales se explican, parcialmente, como una función de ciertas variables ambientales como la temperatura, la precipitación, la pendiente del terreno, la velocidad de viento y el tipo de vegetación. También, se planteó la hipótesis de que los focos de calor tienden a incrementarse significativamente en el tiempo y el espacio.

En México, se han realizado muchos estudios sobre incendios forestales, pero sin un enfoque unificado que modele en el espacio y el tiempo e integre y cuantifique los efectos de ciertas variables ambientales en el riesgo de incendios. Es posible incluir en investigaciones futuras variables relacionadas con las actividades humanas que igual influyen en la distribución de los focos de calor. Cabe señalar que en el presente estudio se analizan, por primera vez, los datos de los focos de calor de dos décadas de registros satelitales en Oaxaca.

Materiales y Métodos

Resultados y Discusión

En el Cuadro 1 se presentan la media a posteriori, la desviación estándar a posteriori y el intervalo de credibilidad a posteriori a 95 % de cada elemento de (β) y para α. Solo cuatro de siete variables ambientales fueron estadísticamente significativas (aquellas cuyos intervalos no contienen el cero), para explicar la tasa logarítmica de incendios forestales. Estas fueron TempCS e IVM con efecto positivo, y la velocidad del viento y el fenómeno de El Niño, con efectos negativos. Resulta fácil interpretar los parámetros en la escala original; es decir, exp⁡(β^). Por ejemplo, un aumento de una unidad estandarizada de TempCS o una unidad de IVM implica un incremento de 78.1 % ((1.781-1) ×100)) y de 34.9 % en el riesgo de desarrollarse incendios forestales, respectivamente.

En otros estudios también se documenta que la temperatura aumenta el riesgo de incendios forestales (^{Perez-Verdin et al., 2014}; ^{Antonio y Ellis, 2015}). Por otra parte, los valores promedio del IVM proporcionan información acerca de las reservas de biomasa que potencialmente podrían estar disponibles como combustible forestal en la temporada seca del año; por lo tanto, conforme se incrementan los valores de IVM es mayor el efecto sobre el peligro de incendios forestales.

Contrariamente a lo que se esperaría, el aumento de una unidad estandarizada en la velocidad del viento disminuye el riesgo de un incendio forestal en 29.4 % ((1-0.706) ×100)). Una posible explicación es que el viento constituye un factor de peligro de incendio, es decir, a medida que aumenta su velocidad, también se incrementa la tasa de propagación del fuego, pero el viento solamente causa incendios directamente cuando los vientos fuertes generan incendios secundarios durante los incendios forestales. Además, la velocidad del viento puede interactuar con otras variables no consideradas en este análisis.

El efecto asociado a El Niño requiere de especial atención. Cuando se observa el fenómeno de El Niño, al año siguiente se reduce en 40.7 % el riesgo de que se presenten focos de calor. Lo anterior se debe a que El Niño ocasiona, en general, más precipitación en México durante el invierno (^{Bravo-Cabrera et al., 2017}), y en consecuencia, hay menos combustible vegetal seco durante la temporada de secas (en la primavera del año siguiente).

En la Figura 2 se representa la moda a posteriori del efecto espacial principal ζi=expui+νi. Se observan dos agrupamientos con los valores más altos de dicho efecto espacial. El agrupamiento corresponde a la banda que se extiende desde el noroeste hacía el sureste del estado de Oaxaca, que inicia en la sierra Mazateca, cruza la Sierra Juárez cerca del centro del estado, y se extiende hasta la sierra Mixe y el Istmo de Tehuantepec. Es importante señalar que esta zona es el final de la Sierra Madre Oriental que cruza todo el territorio nacional; en estas latitudes, el bosque es del tipo trópico seco y hasta 70 % de las especies vegetales que la habitan son endémicas. Es probable que los valores más altos del efecto espacial principal (ζ=max⁡{ζi}) respondan a los siguientes efectos antropogénicos: prácticas de roza, tumba y quema, cría de ganado, explotación forestal ilegal, en general, a cambios de uso del suelo (^{Galván y Magaña, 2020}), y no a variables ambientales.

Figura 2 Moda a posteriori del efecto espacial principal ζi=expui+νi de eventos de focos de calor en Oaxaca, México. 

Algunas celdas con valores altos de ζ cercanas al estado de Veracruz (zona inferior de la cuenca del río Papaloapan) se explican por la presencia de agricultura extensiva y cría de ganado. En el estado vecino de Chiapas, la composición arbórea alterada de muchas selvas caducifolias es el resultado de incendios forestales y de la cría de ganado bovino y caprino (^{Miranda, 2015}).

En México, muchas selvas caducifolias originales son una mezcla de especies arbóreas adaptadas al fuego y sensibles a este, de tal manera que el fuego tiende a seleccionar a las primeras y las segundas se eliminan (^{Rodríguez et al., 2019}). Es importante señalar que la selva de Chimalapas, región con gran biodiversidad, tiene valores del efecto espacial pequeños a pesar de que las zonas aledañas presentan valores altos. Probablemente, los esfuerzos de conservación de su población indígena nativa expliquen los valores bajos. En 1998, año de incendios extremos en México, la región estuvo severamente afectada por incendios de gran magnitud.

El segundo agrupamiento va del suroeste (limítrofe con Guerrero) al sureste del estado, una zona que comprende la Sierra Madre del Sur, la cual está constituida por bosques cálidos secos y cordilleras templadas; es probable que los valores altos del efecto espacial principal se deban también a efectos antropogénicos. Ambos agrupamientos coinciden con los identificados por ^{Zúñiga-Vásquez et al. (2017)}, quienes utilizaron datos de todo el territorio mexicano.

En el Cuadro 2 se presenta un resumen de los hiperparámetros. El cálculo de fspatial=σu2/σν2=0.9998 indicó que 99.98 % de la variación observada puede explicarse por el término de estructura espacial u.

Las Figuras 3a y 3b muestran la moda a posteriori de los efectos temporales estructurados γt y no estructurados ϕt, respectivamente, registrados en la escala original; esto es, exponenciando. Para ambas cantidades, se presenta el intervalo de credibilidad a posteriori (ICA) de 95 %. En la Figura 3a es evidente que los efectos estructurados tuvieron una tendencia a incrementarse de 2001 a 2005. Después, durante un periodo de 10 años, esta se estabilizó, incluso disminuyó marginalmente; de 2015 a 2019 se observó una nueva tendencia ascendente. El incremento de 2001 a 2019 de los efectos estructurados fue de 42.2 %. Sería importante confirmar esta señal de alarma en el futuro, ya que el cambio climático y el calentamiento global están afectando los patrones de incendios en algunas regiones, quizás con aumentos a nivel local, tanto espacial como temporalmente.

a) Efecto temporal estructurado expγt; y b) Efecto no estructurado expϕt.

Figura 3 Moda a posteriori e intervalo de credibilidad de la tendencia temporal para los eventos de focos de calor en Oaxaca, México. 

Las proyecciones de ^{Senande-Rivera et al. (2022)} apuntan a que el área global susceptible de experimentar incendios frecuentes podría incrementarse en 29 % (25 % en las zonas templadas) a finales del siglo XXI. Un efecto del aumento en la incidencia de incendios forestales en el futuro, sería una mayor pérdida de cobertura en zonas boscosas que alteraría los procesos hidrológicos (^{Ruíz-García et al., 2022}), además de tener impactos sociales y ambientales. Por otra parte, el término del efecto no estructurado (Figura 3b) presentó cierta fluctuación en torno a uno como se esperaba, debido a que representa la variación temporal (residual) no explicada por el modelo.

Las probabilidades de excedencia se representan en las Figuras 4a-4d. Las probabilidades de excedencia corresponden a la probabilidad de que la tasa relativa de que en cierta celda i sea más alta que el valor c, matemáticament Pζi>c|y. Se observa que casi la mitad del estado tiene una probabilidad de excedencia superior a 0.75 (75 % de probabilidades) de presentar al menos dos eventos de focos de calor en un año (Figura 4a). En los casos más extremos, pocas celdas en el Istmo de Tehuantepec registraron por lo menos una probabilidad de 0.75 de exceder 12 eventos por año (Figura 4d).

a) Para c=1; en b) para c=4; en c) para c=8; y en d) para c=12.

Figura 4 Probabilidad de excedencia Pζi>c|y, para la tasa relativa de eventos de focos de calor en Oaxaca, México. 

En la Figura 5 se representa la moda a posteriori para la interacción espacio-temporal en la escala original; es decir, expδit. Cuando expδit>1, implica que los registros de eventos de focos de calor de la celda i en el año t fueron superiores al promedio (α^=1.165). Por otra parte, cuando expδit<1, los registros de focos de calor fueron menores que el promedio.

Figura 5 Moda a posteriori de la interacción espacio-temporal expδit. 

A través de los años, las celdas naranja y rojas (con valores altos de expδit) cubren año con año la mayoría de la superficie del estado. El incremento en el número de celdas naranja y rojas es pequeño, pero perceptible y concuerda con la tendencia temporal principal que se observa en la Figura 3a; es decir, el número de eventos está aumentando en espacio y en tiempo, lo cual es una señal de alarma que exige adoptar políticas de manejo de incendios más eficientes, tanto a los gobiernos locales como al federal; de modo que en Oaxaca, se reduzca y mitigue el impacto de los regímenes alterados de incendios.

Conclusión

Se realiza un análisis espacio-temporal de los eventos de focos de calor registrados en el estado de Oaxaca en los últimos 19 años mediante datos del sensor MODIS, con el propósito de investigar si los incendios forestales ocurren en forma dispersa, o si se forman agrupamientos. Los resultados evidencian que las variables ambientales como la temperatura media del cuatrimestre más seco, los valores promedio del Índice de Vegetación Mejorada, velocidad media del viento, y la incidencia de la Oscilación del Sur El Niño, explican algunas de las variaciones espaciales observadas en los eventos de focos de calor. De acuerdo al componente espacial del modelo estadístico empleado, los focos de calor forman agrupamientos, sobre todo a lo largo de la Sierra Juárez y hasta el Istmo de Tehuantepec, pero también a lo largo de toda la Sierra del Sur (hacia la costa del estado). A partir del análisis, se identifica una señal de alarma: el término temporal muestra una tendencia no lineal que se ha incrementado en 42.2 % en el periodo analizado, lo que altera aún más los regímenes de incendios. El parámetro de las interacciones espacio-temporales también confirma esa tendencia, y demuestra que el número de eventos está incrementándose en el espacio y el tiempo. De continuar el aumento, en las próximas décadas cientos de hectáreas de bosque y su biodiversidad estarán amenazadas; ello incluirá las actividades económicas y la salud de las personas que viven tanto en áreas rurales, como urbanas.

Este trabajo puede ser un punto de partida para investigaciones posteriores encaminadas a entender mejor la distribución espacial y temporal de los incendios forestales y los regímenes alterados en Oaxaca. La identificación de los principales agrupamientos permitirá a la Conafor diseñar una mejor logística para actuar de manera oportuna, y así reducir los impactos ambientales y socioeconómicos de los incendios forestales, además de optimizar los recursos económicos y técnicos que se requieren para combatir los incendios.

La probabilidad de excedencia revela algunas regiones muy locales en las que se espera un número alto de eventos por año (>4 en un año), por lo que sería interesante responder desde un punto de vista causal por qué ocurre esto, así como diseñar acciones de mitigación eficaces. La tendencia temporal creciente puede revertirse de muchas maneras, como la concientización mediante campañas completas y bien planeadas que se centren en las poblaciones rurales que más interactúan con el bosque. También se requiere proporcionar a los agricultores capacitación técnica para cambiar a prácticas agrícolas que conlleven un menor riesgo de provocar incendios forestales. La explotación ilegal de los bosques y el cambio de uso del suelo se neutralizan eficazmente a partir de la generación del desarrollo rural a través de actividades rentables como el ecoturismo, la agroforestería, la incorporación de las comunidades rurales en los programas de pago por servicios ambientales, entre otros.

Desde un punto de vista estadístico, en el futuro se propone ampliar el análisis mediante la incorporación de otras variables ambientales, además de las de tipo socioeconómico para cuantificar indirectamente los efectos antropogénicos de la pobreza y de la desigualdad social en el riesgo de incendios forestales y de la alteración de sus regímenes.