Spanish (pdf)
Article in xml format
Article references
Automatic translation
Send this article by e-mail
Cited by SciELO 
Similars in
SciELO 
Permalink
1 Introducción
El aumento en la frecuencia e intensidad de los incendios forestales en el mundo ha generado una creciente preocupación debido a sus impactos ecológicos, económicos y sociales (Westerling et al., 2006; Jones et al., 2022). Para abordar esta problemática, la identificación y monitoreo de áreas propensas a incendios se ha convertido en una herramienta crucial. Los avances en tecnología satelital, como los satélites Terra y Aqua de la NASA, han permitido recopilar un acervo histórico de 24 años de datos valiosos sobre incendios forestales (Justice et al., 2002; Wooster et al., 2021). El potencial de estos datos, junto con la información de otros satélites como Suomi NPP y NOAA-20, puede alimentar modelos predictivos que generan mapas de probabilidad de ocurrencia de incendios forestales (Schroeder et al., 2014; Csiszar et al., 2014; Wooster et al., 2021).
Entre estos modelos predictivos se encuentra la generación de mapas de probabilidad de ocurrencia de incendios forestales antropogénicos, los cuales pueden mejorar la prevención de incendios y la planificación de su combate (Chuvieco et al., 2010; González et al., 2023). Al mejorar la planificación, estos mapas permiten una asignación de recursos más adecuada (Finney, 2005), así como la identificación de áreas vulnerables, como las habitadas, la infraestructura crítica y los ecosistemas que requieren protección (Preisler et al., 2004). Adicionalmente, los mapas de probabilidad de ocurrencia de incendios forestales son herramientas para la educación y generación de conciencia pública, pues sensibilizan a las comunidades sobre la importancia de las prácticas seguras del manejo de fuego y las medidas preventivas (Radeloff et al., 2005). Las políticas basadas en datos precisos y actualizados sobre la probabilidad de incendios son más efectivas para reducir el riesgo y gestionar los recursos de manera sostenible (Calkin et al., 2011).
Sin embargo, muchos países en vías de desarrollo, como Costa Rica, carecen de cartografía detallada sobre la probabilidad de incendios, a pesar de enfrentar una creciente problemática en este ámbito. En el país, desde la década de los cincuenta, las políticas socioambientales, impulsadas por los subsidios internacionales, fomentaron la expansión agrícola y ganadera, lo cual fragmentó el paisaje (Morera et al., 2018). Se sabe que esta fragmentación, caracterizada por la división de ecosistemas forestales en parches aislados, incrementa la probabilidad de ignición la propagación del fuego debido al aumento de bordes, la baja humedad y la proximidad con áreas de manejo antrópico (Laurance y Williamson, 2001; Cochrane, 2002; Silva-Junior et al., 2018). Los incendios, a su vez, agravan la fragmentación, pues crean un ciclo de retroalimentación positiva que conduce a degradar continuamente el bosque (Bowman et al., 2009). Este ciclo se ha intensificado en las últimas décadas en zonas tropicales, incluyendo Costa Rica (Nepstad et al., 1999; Cochrane y Laurance, 2008; Archibald et al., 2013; Duran et al., 2023).
Comprender la relación entre las probabilidades de ocurrencia de los incendios forestales y la configuración del paisaje es fundamental en la lucha por la prevención y su combate. Explícitamente, la cartografía de la probabilidad de ocurrencia de incendios forestales se erige como una herramienta fundamental en la gestión integral de estos eventos, lo cual impacta de forma positiva en diversas áreas clave: prevención, vigilancia, educación y combate. La capacidad para identificar zonas de mayor riesgo permite una asignación estratégica de recursos, acciones preventivas y de vigilancia en dichas áreas. Esto permitiría, en un mediano plazo, la implementación de programas educativos y de concientización dirigidos a las comunidades locales, así como en el fortalecimiento de la vigilancia y el control de actividades potencialmente causantes de los incendios forestales (Finney, 2005; Cochrane y Laurance, 2008).
Por tanto, el presente estudio tuvo como objetivo estimar la probabilidad de ocurrencia de incendios forestales antropogénicos en la Península de Nicoya, actualmente Área de Conservación Tempisque considerando el Corredor Biológico de Chorotega (CBCh), utilizando el método de los pesos de evidencia. Así como también evaluar su relación con la configuración del paisaje en términos de número de fragmentos y densidad de borde, representado a través de la cobertura forestal nativa y el uso agropecuario. La cartografía resultante proporcionará información crucial para diseñar estrategias de prevención de incendios a corto y mediano plazo, así como impulsar la resiliencia comunitaria en un contexto de calentamiento global.
2. Materiales y métodos
2.1 Área de estudio
La península de Nicoya, ubicada en las coordenadas 10°05 ‘46.9 “N, 85°27’ 57.1”W (Figura 1), es una zona de conservación definida como el Área de Conservación Tempisque (ACT), dentro de la ecorregión del mismo nombre en el Pacífico noroeste. Dicha área es parte del Sistema Nacional de Áreas de Conservación (SINAC), con una extensión de 521,048 hectáreas. En la cual existen 29 áreas silvestres protegidas que se conectan a través del CBCh. Su territorio comprende las provincias de Guanacaste y Puntarenas. Con los cantones de Santa Cruz, Nicoya, Hojancha, Nandayure y Puntarenas. En particular, en la provincia de Guanacaste, al norte de la Península de Nicoya, existe un foco de atención debido a la problemática de los incendios forestales. Los datos revelan que esta zona ha sido particularmente afectada, por lo que concentra cerca del 75 % de las tierras quemadas por incendios forestales en Costa Rica durante el período 19982018 (Picado y Cruz, 2021).
Fuente: Elaboración propia.
Figura 1 Ubicación del Área de Conservación Tempisque, en la Península de Nicoya, Costa Rica.
El ACT se encarga de preservar y proteger la cuenca media y baja del río Tempisque. Presenta dos estaciones a lo largo del año: la seca y la lluviosa. La temporada seca abarca desde diciembre hasta mediados de mayo, con una duración de seis meses en el noroeste de Guanacaste (Santa Rosa) y cinco meses en el sureste, a lo largo del río Tempisque (Gillespie et al., 2000; Kappelle, 2016). Las precipitaciones anuales varían entre 1.200 y 2.400 mm (Fallas et al., 2015). En el sur de la península, las lluvias son más copiosas y los períodos de sequía son menos severos que en el resto de la Región Norte. Según Mateo-Vega (2001), la temperatura media anual en la cuenca del río Tempisque es de 27.4 °C. Por otro lado, en la parte norte de la península se registran las temperaturas más altas, superando los 30 °C, con una media mensual de 36 °C en lugares como Santa Elena o Santa Rosa, a lo largo de la costa del Pacífico (Asch et al., 2000).
Gómez y Herrera (1986) describen un total de nueve tipos principales de vegetación, que incluyen manglares, áreas pantanosas, bosques estacionales caducifolios y bosques perennifolios, estos últimos se ubican principalmente en el extremo sur de la Península de Nicoya. En las áreas bajas del Pacífico Norte del ACT, se pueden observar bosques secos y otros tipos de vegetación estacional, que se extienden desde el nivel del mar hasta altitudes medias de 400 m, e incluso hasta los 700 m en algunas áreas elevadas locales. Dicha vegetación está inmersa en una matriz de pastizales, principalmente de la especie exótica Hypharhernia rufa, llamada localmente jaragua.
Finalmente, el único ejemplo previo de vegetación adaptada al fuego se halla en la localidad de Santa Elena. Allí se encuentra un tipo de vegetación herbácea a arbustiva, semiárida y muy distintiva, la cual se caracteriza por matorrales espinosos y árboles pequeños y dispersos (Elizondo y Jiménez, 1988). Esta vegetación ha sido, en gran parte, generada por quemas antropogénicas que han ocurrido durante al menos cuatro siglos (Kappelle, 2016).
2.2 Puntos de calor del sensor MODIS
En el presente estudio se emplearon los productos derivados del espectrorradiómetro de imágenes de resolución moderada MOD04A1 y MYD04A1, diseñados para detectar anomalías térmicas, derivados del satélite MODIS para el periodo 2000-2024 de la NASA Earth data Cloud (https://earthdata.nasa.gov/earth-observation-data/near-real-time/frms/active-fredata). Estos productos han sido previamente utilizados en diversos estudios para identificar la ubicación de los incendios forestales con una resolución de 1 km (Farfán et al., 2021; Eskandari et al., 2020). Fueron seleccionados solo los puntos de calor que coincidían con los sitios de cubierta forestal y se excluyeron aquellos que estuvieran en zonas agropecuarias. Para ello, se obtuvo el mapa de vegetación y uso del suelo del año 2021 del ACT, elaborado por el SINAC, escala 1:20000.
2.3 Estimación de la probabilidad de ocurrencia de los incendios forestales
2.3.1 Variables antropogénicas como fuentes de ignición
Se integró información cartográfica relacionada con seis variables socioambientales en el ACT (Tabla 1). Estas han sido empleadas en otros estudios para estimar la probabilidad de incendios forestales en México y Brasil (Silvestrini et al., 2011; Farfán et al., 2018; Sahagún et al., 2022). Se utilizaron los mapas de vegetación y uso del suelo del 2012 y 2021 a escala 1:20000. Ambos fueron elaborados por el SINAC, el primero en el 2015 y el segundo en el 2021. Asimismo, el modelo digital del territorio, descargado del portal web del Sistema Nacional de Información Territorial, tiene una resolución de tamaño de píxel de 30 metros (SNIT).
Tabla 1 Variables cartográficas socioambientales empleadas en la calibración del modelo de probabilidad para el ACT, con una resolución de 30 metros.
| Insumo cartográfico | Descripción | Fuente |
| Modelo digital de elevación (MDE) | Recorte del modelo original obtenido del sensor ASTER. | https://asterweb.jpl.nasa.gov/gdem.asp (2009) |
| Pendiente | Inclinación del terreno calculada en grados. | Derivado del Modelo Digital de Elevación. |
| Distancia a vías de comunicación | Distancia euclidiana a las principales vías de comunicación. | Elaboración propia con base en el mapa de carreteras del SNIT (2016). https://www.snitcr.go.cr/ |
| Distancia a los asentamientos humanos | Distancia euclidiana a los asentamientos humanos. | Elaboración propia con base en el mapa del SNIT (2015). https://www.snitcr.go.cr/ |
| Distancia al uso de suelo agropecuario | Distancia euclidiana al uso de suelo agropecuario. | Elaboración propia con base en el mapa del SINAC (2021). https://www.sinac.go.cr/ES/bimapas/Paginas/mabosq2021. aspx |
| Distancia al borde del bosque (2021) | Distancia euclidiana al borde del bosque. | Elaboración propia con base en el mapa del SINAC (2021). https://www.sinac.go.cr/ES/bimapas/Paginas/mabosq2021. aspx |
Fuente: Elaboración propia.
Todas las variables socioambientales se trabajaron en el sistema de información geográfica QGIS versión 3.28.13 Firenze, en formato ráster con una resolución de 30 metros en proyección UTM, zona 16 norte.
2.3.2 Estimación de la probabilidad de ocurrencia de los incendios forestales
Se empleó el método de los pesos de evidencia, incorporado en la versión 7.8.0 del programa DINAMICA EGO (Soares-Filho et al., 2006).1 Para ello, se tomaron en cuenta las variables socioambientales enumeradas en la tabla 1, las cuales se integraron y ajustaron durante el proceso de calibración para la construcción del modelo espacialmente explícito. También se emplearon los mapas de cubierta forestal de los años 2012 y 2021, siendo el tiempo inicial y final de calibración del modelo. Se utilizaron también los puntos de calor del año 2012 para la estimación de los pesos de evidencia. Dicho método está basado en las probabilidades condicionales y es un algoritmo bayesiano (Bonham-Carter, 1994; Soares-Filho et al., 2004). Se le llama probabilidad condicional a la probabilidad de que un suceso se cumpla habiéndose cumplido otro previamente. El peso positivo de evidencia, w+, asociado a la presencia de la condición considerada, se calcula según la ecuación (1):
donde 𝑃(𝐶|𝐷) es la probabilidad de ocurrir un evento 𝐶, dado la condición de la presencia del evento 𝐷; y 𝑃(𝐶|¬𝐷) es la probabilidad de que ocurra un evento 𝐶, dado la condición de la ausencia del evento 𝐷. Cuando la presencia de la condición tiende a aumentar la probabilidad del evento el valor del peso w+ es positivo. Por otro lado, cuando la presencia de la condición disminuye, la probabilidad del evento tiene un valor negativo. El cálculo de la probabilidad condicional, al tomar en cuenta varias condiciones, suma los pesos de evidencia. Para esto, el supuesto de independencia condicional entre las variables es requerido, por lo que se calculó el coeficiente de correlación de Cramer entre las variables explicativas.
2.3.3 Evaluación de la probabilidad de ocurrencia de incendios forestales
La evaluación del modelo se centró en determinar la concordancia espacial entre las predicciones de probabilidad a incendios forestales estimadas y los puntos de calor correspondientes a los incendios forestales ocurridos en el año 2022. Esta información se consideró independiente a la fase de calibración del modelo. Este proceso de evaluación se llevó a cabo mediante el análisis estadístico de la curva característica operativa del receptor (ROC), una técnica comúnmente empleada en estudios de predicción de incendios forestales (Sahagún et al., 2022; Farfán et al., 2021; Mas et al., 2013). Los resultados de la curva ROC, que varían entre 0.5 y 1, fueron interpretados de la siguiente manera: un valor de 0.5 corresponde al desempeño de un modelo nulo, mientras que los valores superiores a 0.5 indican una mejora en el rendimiento del modelo. Valores de la curva ROC entre 0.5 y 0.6 sugieren un modelo poco efectivo; entre 0.6 y 0.7, una capacidad de predicción baja; de 0.7 a 0.8, una capacidad de predicción aceptable; de 0.8 a 0.9, una capacidad de predicción muy buena; y de 0.9 a 1, una capacidad de predicción excelente (Mas et al., 2013).
2.4 Evaluación de la fragmentación del paisaje
La pérdida y fragmentación de la cobertura forestal, así como la exposición a bordes, incrementan la ocurrencia e intensidad de incendios forestales en los bosques tropicales (Armenteras et al., 2017; Silva-Junior et al., 2018; Guedes et al., 2020). Por esta razón, se estimaron dos métricas de composición (los porcentajes de cobertura del bosque tropical y de los pastizales ganaderos) y dos métricas de configuración (densidad de borde y de parches) del paisaje a lo largo del ACT (Tabla 2). Las métricas fueron calculadas a partir de una gradilla de 7 ha, la cual representa el tamaño promedio de las fincas de manejo agropecuario en el ACT. Para expresar los resultados de dicho análisis se consideró el diseño del Corredor Biológico Peninsular dentro del ACT.
Tabla 2 Métricas de paisaje estimadas para la evaluación de la configuración del paisaje en el ACT.
| Métrica | Tipo | Fórmula | Relación |
|---|---|---|---|
| Cobertura de bosque | Composición |
 a
ij
: Área (m2) del parche j de la cobertura i en el paisaje. A: Área total del paisaje (m2). Unidades: porcentaje. |
Hay una relación inversa entre la cobertura de bosque y la probabilidad de ignición. |
| Uso agropecuario | Composición | El manejo de tierras agropecuarias en países de América Latina se relaciona con el uso del fuego para su creación o mantenimiento. | |
| Densidad de borde | Configuración |
 E: Perímetro de borde en el paisaje. A: Área total del paisaje (m2). Unidades: metros/hectárea. |
Existe una relación directa entre la exposición a bordes y la desecación de la vegetación arbórea, con lo cual se promueve la propagación del fuego al interior del bosque. |
| Densidad de parches | Configuración |
 n
i
: Número de parches de la cobertura i en el paisaje. A: Área total del paisaje (m2). Unidades: parches/100 hectáreas. |
Paisajes tropicales fragmentados favorecen la propagación del fuego a lo largo del paisaje. |
Fuente: Elaboración propia.
Se analizaron las diferencias en los patrones de estructura del paisaje entre aquellos con alta y baja probabilidad de incendios. Se consideraron paisajes con alta probabilidad de incendios aquellos que tuvieran una probabilidad promedio ≥ 0.47 (valor que representa el 75 % de los valores más altos). Las diferencias se evaluaron con una prueba t de Welch para varianzas heterogéneas (Anderson y Burnham, 2002). Posteriormente, se analizó la relación lineal entre las probabilidades antropogénicas de incendios forestales y las métricas de paisaje a lo largo del ACT a través de correlaciones de Spearman.
3. Resultados
El análisis de los puntos de calor del sensor MODIS entre el 2000 y el 2024 en zonas forestales del ACT reveló un aumento del 62.96 % en el número total de igniciones durante la última década. Esta tendencia ascendente se refleja también en la regresión lineal presentada en la figura 2B. En la figura 3, se muestran las probabilidades estimadas de incendios forestales de origen antropogénico en el Área de Conservación Tempisque (ACT).
Fuente: Elaboración propia.
Figura 2 A) Igniciones forestales anuales para el periodo 2000-2024 registradas por sensor MODIS para el ACT. B) Regresión lineal estimada con las igniciones forestales registradas por el sensor MODIS (20002024) en el ACT.
Destaca la Zona Norte, Guanacaste, con un riesgo significativamente mayor al alcanzar el 100 %, en comparación con la zona sur. Además, se presentan cuatro clases de riesgo a incendios forestales en el Corredor Biológico Chorotega (CBCh) y áreas silvestres protegidas asociadas. Siete áreas silvestres protegidas presentan un riesgo muy alto: Refugio Nacional de Vida Silvestre Iguanita, Refugio Nacional de Vida Silvestre el Conchal, Parque Nacional Barra Honda, Refugio Nacional de Vida Silvestre Hacienda El Viejo, Humedal Palustrino Corral de Piedra, Refugio de Vida Silvestre Laguna Mata Redonda, Parque Nacional Marino las Baulas de Guanacaste y el Parque Nacional el Diriá. En contraste, las áreas de menor riesgo son la Zona Protectora Península de Nicoya, Zona Protectora Monte Alto y Reserva Natural Absoluta Cabo Blanco.
Fuente: Elaboración propia.
Figura 3 A) Probabilidad antropogénica de ocurrencia a incendios forestales en el ACT. B) Clases de riesgo a incendios forestales antropogénicos en el Corredor Biológico Chorotega.
Asimismo, se identificó una relación entre la probabilidad de ocurrencia de incendios y la estructura del paisaje. Paisajes con una alta probabilidad de incendios presentaron un mayor porcentaje de uso de suelo agropecuario, así como un menor porcentaje de bosques (Figura 4). Se estimó que el tamaño del parche más pequeño con alta probabilidad de ignición antropogénica es de 6.26 hectáreas. Además, la configuración del paisaje también estuvo asociada con la probabilidad de incendios, pues los paisajes con una alta probabilidad de ignición antropogénica tuvieron mayores densidades de borde y número de fragmentos (Figura 4C y D).
Fuente: Elaboración propia.
Figura 4 Comparaciones entre probabilidades altas y bajas y estructura del paisaje en términos de a) porcentaje de cobertura de bosque b) porcentaje de uso agropecuario, c) densidad de borde y d) densidad de parche (*p<0.001).
La Figura 5 muestra el grado de fragmentación y densidad de borde a lo largo del corredor biológico peninsular, siendo mayor en el norte que en el sur. Se encontró solamente una porción del CBP que conecta con las áreas silvestres protegidas Monte Alto y el Diriá, con baja fragmentación y densidad de borde en todo el ACT.
Fuente: Elaboración propia.
Figura 5 A) Grado de fragmentación en el CBCh en el ACT. B) Densidad de borde en los parches que integran el CBCh en el ACT.
Los análisis de correlación mostraron una relación inversa entre la probabilidad de incendio forestal antrópico y el porcentaje de cubierta forestal, contrastando con la relación directa observada con el uso de suelo agropecuario (Figura 6 A y B). Además, se identificó una asociación positiva entre la probabilidad de incendio forestal antrópico con el grado de fragmentación del bosque, así como con la densidad de borde. Esto sugiere que un aumento en la probabilidad de incendio forestal se relaciona con un mayor grado de fragmentación y densidad de borde en el bosque (Figura 6 C y D).
Fuente: Elaboración propia.
Figura 6 A. Correlación negativa entre la probabilidad antropogénica de incendios forestales y el porcentaje de cubierta forestal. B. Correlación positiva entre la probabilidad antropogénica de incendios forestales y el porcentaje de uso del suelo agropecuario. C y D. Correlaciones positivas entre la probabilidad antrópica de incendios forestales con la densidad de borde y la densidad de fragmentos
3.1. Evaluación del modelo de probabilidad antropogénica de incendios forestales
Con la prueba estadística de la curva ROC se obtuvo un valor de área bajo la curva (AUC) de 0.72, lo cual indica una concordancia del 72 % entre las probabilidades más altas predichas por el modelo y los incendios ocurridos en el año 2022 (Figura 7).
4. Discusión
El territorio del Área de Conservación Tempisque (ACT) en Costa Rica, ubicado en la región Chorotega y Puntarenas, destaca por su creciente ocurrencia de incendios forestales. En la última década, se ha registrado un incremento de más del 50 % en el número de incendios forestales. Está tendencia coincide con la observada en otros ecosistemas tropicales como el Amazonas, así como en regiones el oeste de Estados Unidos y el sureste de Australia (Jones et al., 2022).
Esta problemática no solo representa una amenaza directa para la biodiversidad y los servicios ecosistémicos de la región, sino también la regulación del ciclo hidrológico, la protección del suelo y la captura de carbono (Turner, 2010). Además, la creciente incidencia de incendios forestales plantea desafíos importantes para gestionar la prevención y el control del fuego, tanto en el CBCh como en el interior de las 29 áreas silvestres protegidas.
Las probabilidades antropogénicas de ocurrencia de los incendios forestales más altas se distribuyeron en la zona norte del ACT, en donde el riesgo a incendios para el CBCh es alto. Janzsen (1986) documentó, para esta región de Guanacaste, cómo la conversión de bosques secos caducifolios en pastizales fue un proceso que incrementó la ocurrencia de incendios forestales para dicha zona. Esta transformación del paisaje, con la introducción de especies altamente inflamables como el pasto jaragua (Hyparrhenia rufa) para el ganado y la presencia de Typha dominguensis en el humedal de Palo Verde, han creado las condiciones para que el fuego se propague durante la época seca del año (Rojas-Chavez et al., 2015). Si bien las probabilidades antropogénicas de ocurrencia de incendios forestales son menores en la zona sur del ACT, región de Puntarenas, esta presenta condiciones antropogénicas de uso de fuego como la quema de los pastizales para el ganado y un nivel medio de fragmentación del bosque que requieren impulsar una gestión integrada del manejo del fuego. Esta estrategia debe priorizar acciones de paisaje que integren fragmentos de bosque de al menos seis hectáreas con alta densidad de borde, al enfocarse en la regeneración natural y la eliminación de pastizales, especialmente en sitios del Corredor Biológico de Chorotega (CBCh). Además, es crucial promover prácticas agrícolas y ganaderas sostenibles que no dependan de la quema del bosque. Otra propuesta importante se relaciona con los planes de manejo del fuego para las áreas silvestres protegidas con estrategias para el manejo de combustibles, la construcción de brechas corta fuego y la sensibilización de la población frente a esta problemática (Pérez-Salicrub et al., 2020).
Los resultados evidencian una clara relación entre la fragmentación del paisaje y el incremento en la probabilidad de incendios forestales dentro del ACT. Estos hallazgos corroboran la hipótesis de que la alteración de la estructura del paisaje, derivada de las prácticas de uso del suelo, facilita la propagación del fuego hacia los fragmentos de bosque remanentes. La fragmentación y la consecuente exposición a bordes aumentan la desecación del bosque, lo cual provoca condiciones propicias para la propagación del fuego en su interior (Laurance, 2004). Estudios previos han estimado que el fuego puede propagarse hasta 2.5 km en bosques con bordes expuestos (Cochrane, 2002), lo que subraya la importancia de considerar la estructura del paisaje en las estrategias de prevención de los incendios forestales.
De acuerdo con Morera-Beita et al. (2018, 2021), entre el 2000 y 2015, el 80 % de los corredores biológicos de Costa Rica aumentaron en la fragmentación del paisaje. Se encontraron coincidencias en el patrón observado para la zona norte del CBCh entre la fragmentación evaluada para el año 2015 y el grado de fragmentación estimada en el 2021.
A menos de que se modifiquen las prácticas actuales de uso de la tierra y el manejo del fuego en el ACT, los incendios forestales podrían transformar vastas extensiones de bosque tropical tanto en el CBCh como al interior de las áreas silvestres protegidas en ecosistemas degradados. Este fenómeno, conocido como “sabanización”, ya está ocurriendo en la Amazonía y representa una amenaza significativa para la biodiversidad, el almacenamiento de carbono y la regulación del clima (Driscoll et al., 2021). Finalmente, además de las actividades agropecuarias ya desarrolladas en el ACT, el desarrollo de nuevos sitios turísticos como marinas y otros proyectos inmobiliarios podrían eventualmente propiciar escenarios de riesgo más complejos a incendios forestales como son lo de interfaz urbano-forestal.
5. Conclusión
En la última década, los incendios forestales en el Área de Conservación Tempisque (ACT) han aumentado en la última década, más del 50 %, según datos del sensor MODIS. El modelo de probabilidad de ocurrencia de incendios antrópicos revela un riesgo del 100 % en la Zona Norte, lo cual pone en riesgo la efectividad para la conservación del CBCh y las áreas silvestres protegidas de la región. Aunque la Zona Sur presenta una menor probabilidad de ignición antropogénica, y un menor riesgo a incendios forestales, la quema de pastizales para ganado y la fragmentación del bosque requieren una gestión integrada del fuego para evitar un escenario similar al de la Zona Norte.
El análisis integrado de la cartografía de probabilidad de incendios antrópicos y la configuración del paisaje en el ACT revela una estrecha relación entre ambos factores en la ocurrencia de incendios forestales, siendo los fragmentos de bosque de 6 hectáreas con un 70 % de probabilidad de ignición antropogénica los más vulnerables a los incendios forestales. Se sabe que el fuego interactúa con la fragmentación a través del efecto de borde del bosque, desecándolo y favoreciendo su propagación. En particular, en los bosques tropicales cuando el fuego y la fragmentación se influyen mutuamente, es posible que se produzcan ciclos de retroalimentación que pueden llevar a la conversión del bosque en pastizales (Driscoll et al., 2021).
Esta información constituye una herramienta fundamental para la gestión del manejo de fuego en el ACT, pues permite diseñar planes específicos para cada corredor biológico y área silvestre protegida, con lo cual se priorizan acciones en zonas de alto riesgo para el CBCh. Estos planes deben incluir estrategias como el manejo de combustibles, la creación de brechas cortafuego y la sensibilización de la población sobre la prevención de incendios. Finalmente, la cartografía aquí presentada es un pilar fundamental para la construcción de una gestión efectiva del fuego para el ACT y el CBCh que contribuya a mitigar el cambio climático al prevenir la liberación de grandes cantidades de carbono almacenado en los bosques tropicales. Futuras investigaciones explorarán la sinergia entre el aumento de la temperatura, las sequías y las igniciones antropogénicas. Los resultados aquí presentados son cruciales para la planificación territorial de las estrategias que fomenten la prevención y el combate de incendios forestales, así como la resiliencia frente al cambio climático en el ACT y el CBCh.
