Ecuación 1
Recibido: 20 de enero de 2021; Aceptado: 17 de septiembre de 2021
Resumen
La degradación de los ecosistemas en los Andes, una de las regiones con mayor biodiversidad del mundo, es atribuida principalmente a la deforestación y la fragmentación. La implementación de áreas protegidas (AP) es una estrategia de conservación efectiva para mitigar los efectos de la transformación del paisaje sobre la biodiversidad. En la Cordillera Occidental de Colombia, en Antioquia, se estableció recientemente el Corredor del Oso de Anteojos (COA, ~4169 km2), como una estrategia para la conservación de la biodiversidad. En este estudio, a partir de información derivada de sensores remotos y bases de datos, localizamos las áreas con mayor biodiversidad en este corredor, y analizamos si estás áreas están siendo protegidas. Los resultados indican que la deforestación (162.37 km2) y la fragmentación (> 90 %) en el COA durante los últimos 19 años se concentran en las áreas más biodiversas, al interior y por fuera de las AP. Si bien las AP alcanzan 30% del área total del COA, las áreas de mayor biodiversidad están representadas en menos de un 17%. Estos resultados indican el riesgo de pérdida de biodiversidad y la creciente necesidad de fortalecer las AP y delimitar nuevas áreas biodiversas.
Palabras clave:
Fragmentación., Heterogeneidad climática, Índices Dinámicos de hábitat, Pérdida de bosque, Vulnerabilidad..Abstract
The degradation of ecosystems in the Andes, one of the most biodiverse regions in the world, is caused by deforestation and forest fragmentation. The implementation of protected areas (PA) is an effective conservation strategy to mitigate the effects of landscape transformation on biodiversity. In the western mountain range of Colombia, in Antioquia, the Spectacled Bear Corridor (COA, ~ 4169 km2) was recently established. We describe the state of forest conservation in the COA and analyze whether areas of high biodiversity are being protected or not, using information obtained from remote sensors and large databases. The results indicate that deforestation (162.37 km2) and fragmentation (> 90%) in the COA during the last 19 years are concentrated in the areas of greatest biodiversity, and occur both inside and outsidePA. Although protected areas reach 30% of the total area of the COA, the areas with the greatest biodiversity represent less than 17%. These results indicate the risk of biodiversity loss and the growing need to strengthen PA and to delimit new biodiverse areas.
Keywords:
Climatic heterogeneity, Dynamic Habitat Indices, Forest loss, Fragmentation, Vulnerability.Introducción
La pérdida de biodiversidad se ha incrementado globalmente durante los últimos 20 años, debido, principalmente, a las presiones antropogénicas, como la deforestación, la fragmentación, los incendios, la sobreexplotación de recursos y el cambio climático, entre otras (Laurance et al., 2012). Este fenómeno ha llevado, a su vez, a la disminución de servicios ecosistémicos indispensables para el sostenimiento y el bienestar de la humanidad (Cardinale et al., 2012). En los Andes tropicales, la región con mayor biodiversidad del mundo (Myers et al., 2000) y que provee servicios ecosistémicos fundamentales para la supervivencia de cerca del 40% de la población de Suramérica (Banco Mundial, 2013), la degradación de sus ecosistemas naturales es atribuida a dos causas principales: (1) la deforestación (~2.6 millones ha. de bosque perdidas anualmente; FAO, 2020) y (2) la fragmentación de los remanentes de bosque (Armenteras & Villareal, 2002; Zhai et al. 2015; Barlow et al., 2016).
Una de las estrategias globales de conservación de la diversidad biológica más efectivas para mitigar los efectos de la transformación del paisaje es la implementación de áreas protegidas (AP) (Brooks et al., 2006; Venter et al., 2018). Actualmente, más de 700 millones de hectáreas de bosque a nivel global se encuentran en AP legalmente establecidas (FAO, 2020). En los últimos 15 años, los mecanismos globales de conservación exigen la identificación, el reconocimiento y la protección de nuevos sitios prioritarios o de importancia para la conservación de la biodiversidad (CBD, 2004; Jackson & Kennedy, 2009). Actualmente, la delimitación de AP requiere incorporar las medidas de importancia biológica (e.g., número de especies endémicas, riqueza de especies, diversidad beta (Myers et al., 2000; Radeloff et al ., 2019) y medidas de vulnerabilidad del hábitat (e.g., pérdida de hábitat, deforestación y fragmentación; Hansen et al., 2013).
Adicionalmente a las medidas de importancia biológica y de vulnerabilidad, en los análisis de planificación de las AP en regiones con amplia heterogeneidad ambiental como los Andes tropicales (Kattan et al., 2004; Ramirez-Villegas et al., 2014), es importante incluir medidas que permitan conocer el cambio en la composición de especies a lo largo de gradientes ambientales de la región (i.e., diversidad beta o complementariedad). No obstante, en regiones tan poco muestreadas como los trópicos, medidas de heterogeneidad climática son útiles para reconocer el recambio de especies, al asumir que dos unidades del paisaje climáticamente diferentes tendrán un conjunto de especies disímiles, particularmente en montañas tropicales, donde los cambios en condiciones climáticas generan un marcado gradiente de diversidad (Ferrier, 2002; Soria-Auza & Kessler, 2008; Slik et al., 2015; Gardner et al., 2020).
Si bien las medidas antes mencionadas (importancia biológica y vulnerabilidad) han sido utilizadas para la delimitación de áreas prioritarias para la conservación de la biodiversidad (Myers et al., 2000; Bruner et al., 2001; Ferrier, 2002; Orme et al., 2005; Brooks et al., 2006), aún hay AP que se encuentran ubicadas en lugares poco representativos de la biodiversidad local (Armenteras & Villareal, 2002; Andam et al., 2008; Forero-Medina & Joppa, 2010; Venter et al., 2018) o tienen un fuerte impacto por la transformación del paisaje (Armenteras & Villareal, 2002; Leisher et al., 2013). En ese sentido, evaluar el estado de conservación y la representación de las áreas más biodiversas en las AP es indispensable para la gestión y la toma de decisiones que lleven a reducir la pérdida de la biodiversidad, bajo el escenario actual de cambio climático y la fuerte transformación del paisaje. En este estudio describimos el estado de conservación, en términos de vulnerabilidad (i.e. deforestación y fragmentación de los bosques), en el Corredor del Oso de Anteojos (Tremarctos ornatus) en el suroeste y occidente de Antioquia, Colombia, y analizamos la representatividad de las áreas de alta importancia biológica en las AP. El Corredor del Oso de Anteojos (COA), ubicado en una de las regiones mejor preservadas de Antioquia, nació como una estrategia social, que ha sido promovida a nivel gubernamental por medio de una ordenanza departamental para la conservación de la fauna silvestre y la flora de la región desde el año 2017.
Materiales y métodos
Área de estudio
El COA delimita ~4169.7 km2 en el suroeste y occidente del departamento de Antioquia, en la Cordillera Occidental de Colombia (Figura S1a). Está ubicado en una de las regiones con mayor extensión de bosques y, aunque ha sido poco estudiado, se considera una de las regiones con mayor diversidad biológica en Antioquia (González-Caro & Vásquez, 2017; Triviño et al., 2018). El COA presenta una precipitación media anual de ~2400 mm.año-1 y una temperatura media anual de ~15 °C (Karger et al., 2017). El corredor abarca principalmente ecosistemas de bosques andinos y altoandinos, en un rango de elevación entre 2600 y 3900 msnm. En el 40% del área del COA (~1676.57 km2) se encuentran seis figuras de protección, que corresponden a trece AP: tres Distritos Regionales de Manejo Integrado (DRMI), un Parque Nacional Natural (PNN), un Parque Natural Regional (PNR), tres Reservas Forestales Protectoras Nacionales (RFPN) y cinco Reservas Naturales de la Sociedad Civil (RNSC) (Figura S1a, Tabla S1). Las figuras de protección de estas AP difieren principalmente en la permisibilidad para la ejecución de actividades económicas dentro de su jurisdicción (Decreto 2372 del 2010, Ministerio de Medio Ambiente y Desarrollo Territorial de Colombia). Los límites de las AP en formato vector fueron obtenidos de los datos abiertos de Parques Nacionales Naturales de Colombia.
Importancia biológica
Para estimar la importancia biológica, combinamos la diversidad local, derivada a partir de los Índices Dinámicos de Hábitat (IDH), y la heterogeneidad climática, como indicador de la variación en la composición biológica. Estos índices han sido usados para describir patrones globales de biodiversidad; sin embargo, a escala de paisaje pueden presentar sesgos, ya que pueden no reflejar apropiadamente los cambios en diversidad. No obstante, la escasez de muestreos en el área de estudio no permite realizar validaciones de los IDH a esta escala. Por otro lado, los IDH son un proxy para conocer los cambios en diversidad a lo largo del tiempo a diferentes escalas espaciales, y son útiles para evaluaciones rápidas de la representatividad y efectividad de las AP (Hobi et al., 2017; Hobi et al., 2021; Politi et al., 2021).
Primero, utilizamos los IDH como una aproximación de la biodiversidad en el COA (Radeloff et al., 2019), los cuales fueron obtenidos del repositorio SILVIS LAB para el promedio de los años 2003 a 2014, con base en el Índice de Vegetación de Diferencia Normalizada (NDVI, por sus siglas en inglés) en formato ráster con una resolución espacial de 1 km x 1 km. Los IDH permiten conocer la riqueza de especies, basados en medidas de productividad bajo los siguientes supuestos (Radeloff et al., 2019): 1) la diversidad de especies aumenta con la productividad (hipótesis de energía; Wright, 1983); IDH1 estimado como la sumatoria de la productividad anual según el NDVI; 2) disminuye hacia los ambientes extremos (hipótesis de filtro ambiental; Currie et al., 2004); IDH2 estimado como el valor mínimo anual del NDVI; y 3) aumenta con la estabilidad ambiental (hipótesis de estacionalidad; Janzen, 1967); IDH3 estimado como la desviación estándar anual del NDVI. Estos tres índices son complementarios entre sí y proporcionan una mirada integral de la riqueza de especies de una región. Para realizar posteriores análisis, se hizo una división del COA en cuadrantes de 1 km x 1 km y se extrajo el valor de cada uno de los IDH para cada cuadrante.
Segundo, estimamos la heterogeneidad climática a partir de las variables bioclimáticas (i.e., temperatura media anual, precipitación anual, estacionalidad de la temperatura, precipitación del cuarto más seco y variabilidad de la precipitación) de WorldClim (Fick & Hijmans, 2017) a una resolución espacial de 1 km x 1 km. A partir de esta información, estimamos tres grupos climáticos para la región de estudio, a través del método de aglomeración conocido como k-medias. Este método establece grupos con varianzas homogéneas distinguibles de otros grupos (Altman & Krzywinski, 2017; Gardner et al., 2020). Finalmente, para validar los grupos climáticos, realizamos comparaciones de la composición de especies de aves, mamíferos y angiospermas entre grupos, usando listados de especies derivados de los registros depositados en el Repositorio Mundial de Información sobre Biodiversidad (GIBIF, por sus siglas en inglés) (Beck et al., 2014) y cuantificamos el porcentaje de similitud entre ellos (i.e., Diversidad Beta, número de especies compartidas/la suma de especies entre los dos sitios; Krebs, 1999).
Por último, calculamos el índice de importancia biológica ubicando los sitios con mayores valores del IDH1 dentro de cada grupo climático, para localizar las áreas de mayor interés. Con el fin de establecer comparaciones entre los grupos climáticos, se escalaron los IDH1 de cada grupo por su mayor valor; de esta forma todos varían entre 0 y 1, siendo 1 los sitios con mayor importancia biológica (i.e., píxeles con mayor riqueza de especies).
Vulnerabilidad
Se estimó la vulnerabilidad a la pérdida de hábitat y conectividad por medio de la cuantificación de la deforestación y la fragmentación. La cuantificación de la deforestación y la fragmentación de los bosques dentro y fuera de las AP se realizó usando información satelital derivada de Landsat (Hansen et al., 2013). Los datos de deforestación fueron obtenidos de los resultados del análisis de series de tiempo de imágenes Landsat realizado por Hansen et al. (2013), depositados en el geoportal Global Forest Change. La capa ráster utilizada fue Year of gross forest cover loss event (i.e., lossyear), la cual tiene una resolución espacial aproximada de 30 m x 30 m y proporciona información sobre eventos de deforestación anual. El ráster de deforestación fue transformado en unos y ceros, donde uno representa un evento de deforestación para cualquier año, y cero, ningún evento de deforestación. Se utilizaron los cuadrantes de 1 km x 1 km para extraer la suma de los pixeles deforestados para cada uno de estos cuadrantes, lo cual nos permitió conocer el valor de los IDH en cada sitio que ha sido deforestado y en los que no ha sido deforestado. Finalmente, con el uso del vector de AP para el COA y el ráster de deforestación, extrajimos la deforestación dentro de las AP y por fuera de ellas. La transformación del ráster de deforestación fue realizada utilizando las funciones de los paquetes Raster y rgdal del programa R.
Para estimar la fragmentación, en primer lugar, obtuvimos el ráster de cobertura de bosque para el año 2019, calculado a partir de un ráster de cobertura del año 2000 (i.e., ráster de referencia) obtenido de la base de datos de Global Forest Change (treecover2000; Hansen et al., 2013). El ráster de referencia tiene una resolución espacial de aproximadamente 30 m x 30 m por píxel y contiene valores de 0 a 100 % de cobertura. Los valores del ráster de referencia fueron transformados a cero y uno, donde los píxeles con cobertura de bosque mayores o iguales a 80 % tomaron valores de uno. Posteriormente, sustrajimos los píxeles que fueron deforestados (i.e., ráster anteriormente obtenido) desde el año 2001 hasta el año 2019. El análisis fue realizado utilizando las funciones del paquete Raster del programa R.
El índice de fragmentación f en el COA se calculó a partir del ráster de cobertura, contando el número de píxeles que contiene bosque (i.e., píxeles iguales a 1,n ); el número de grupos (i.e., uno o más píxeles continuos con bosque, g); el número de píxeles conectados (i.e., aquellos que tienen más de dos píxeles aledaños con cobertura boscosa, c); y el número de píxeles totales (N), dentro de una ventana de análisis de 3 km x 3 km (Ecuación 1). Este índice tiene valores de cero cuando la fragmentación es alta y valores cercanos a uno cuando la fragmentación es nula (González-Caro & Vásquez, 2017). El índice fue obtenido utilizando las funciones de los paquetes sp, raster y maptools del programa R.
Al igual que con el ráster de deforestación, para la fragmentación se utilizaron los cuadrantes de 1 km x 1 km para extraer el valor de la fragmentación para cada uno de estos cuadrantes. Esto nos permite conocer el valor de los IDH y su grado de fragmentación en un área específica. Finalmente, con el uso del vector de AP para el COA y el ráster de fragmentación extrajimos la fragmentación dentro de las AP y por fuera de ellas.
Análisis estadísticos
Para analizar si las áreas de mayor importancia biológica en el COA están representadas en las AP, consideramos las metas de conservación reconocidas globalmente bajo el Convenio de Diversidad Biológica (CDB). Se esperaba que para el año 2020 al menos 17 % de la tierra fuera declarada como AP (meta Aichi 12). Sin embargo, diferentes estudios han demostrado que solo 17 % del territorio puede ser insuficiente para mantener poblaciones viables de muchas especies, por lo cual se ha sugerido incrementar el área hasta 30 % aproximadamente (Venter et al., 2018). De esa manera, calculamos si las áreas con mayor importancia biológica están representadas en al menos 17 % y 30 % en las AP del COA. Para esto ordenamos los valores del mayor al menor y seleccionamos los píxeles correspondientes al 17 % y 30 % del área de estudio, desde el valor más alto para cada uno de los tres grupos climáticos.
Para evaluar el efecto de la deforestación y la fragmentación dentro del COA, comparamos los IDH obtenidos entre 2003 y 2014 entre áreas deforestadas (durante 2014 y 2019) y áreas no deforestadas. De igual forma se hizo para áreas fragmentadas y no fragmentadas. Para esto utilizamos la prueba de Kolmogórov-Smirnov, la cual nos permite comparar distribuciones no paramétricas. Para realizar la prueba, utilizamos la función ks.test del paquete stats del programa R (Lopes et al., 2007; Young, 1977). Este análisis se realizó separando AP y áreas no protegidas, donde los valores perdidos de bosque (en hectáreas) dentro de las AP significan deficiencia en la protección y manejo de dichas áreas y fuera de ellas indica pérdida de conectividad regional. Todos los análisis fueron ejecutados usando R 3.5.1.
Resultados
El área de cobertura boscosa actual en el COA fue de 2720.6 km2, es decir, ~65.24 % de su extensión total son bosques (Figura S1c). De estos, aproximadamente 1265.25 km2 (30,34 % de los bosques actuales y el 40.21 % del área total) se encuentran dentro de alguna figura de protección (Tabla S1).
Importancia biológica
En promedio, la productividad anual acumulativa (i.e., IDH1, derivado del NDVI banda 1) fue de 15.04 ± 3.30, y el 65.42 % (km2) del COA presentó valores mayores o iguales a 15, es decir, áreas de alta productividad y diversidad de especies. Por otro lado, la productividad anual mínima (i.e., IDH2, derivado del NDVI banda 2) fue 0.42 ± 0.23, con solo el 29.84 % del corredor presentando valores de alta diversidad de especies. Finalmente, la variación estacional de la productividad, (i.e., IDH3, derivado del NDVI banda 3) fue 0.19 ± 0.17, donde ~ 86 % del corredor tuvo valores menores o iguales a 0.2, lo cual significa una elevada diversidad de especies en casi toda su extensión.
La heterogeneidad climática del COA nos permitió diferenciar tres zonas climáticas (grupos 1, 2 y 3; Figura S1b). El grupo 2 fue el de mayor extensión con ~1828.66 km2, seguido del grupo 3 (1593.68 km2) y, por último, el grupo 1 (745.70 km2). Aproximadamente el 45.32 % de la extensión del grupo 3 estuvo dentro de un área protegida, a diferencia de los otros dos grupos que tuvieron una menor representación de AP dentro de sus límites, 37.03 % y 36.09 %, grupo 1 y grupo 2, respectivamente.
En cuanto a las comparaciones de la composición de especies entre los tres grupos climáticos, encontramos que estos se diferencian principalmente en su composición de especies de plantas (angiospermas). Como era de esperarse, los grupos taxonómicos con mayor movilidad como aves o mamíferos son muy similares, pues el área del COA es relativamente pequeña. Sin embargo, los grupos climáticos 2 y 3 se diferencian en la composición de mamíferos. Es importante resaltar que el grupo climático 2 es el que tiene un mayor número de registros biológicos, mientras que los otros están pobremente estudiados, lo que reduce la confiabilidad de la comparación de los listados.
Por último, los sitios con mayores valores de importancia biológica para cada grupo climático tuvieron valores de 0.889 ± 0.05, 0.883 ± 0.04 y 0.885 ± 0,05, respectivamente (Figura S1c). Aproximadamente el 49.80 %, 64.06 % y 44.29 % de cada grupo climático, en orden, representan los sitios más relevantes (i.e., valores cercanos a 1) para conservar.
Vulnerabilidad
En los últimos 19 años, el COA perdió en cobertura boscosa alrededor de 162.37 km2 (8.55 ± 5.94 km2.año-1), siendo el 2017 el año que presentó una deforestación cuatro veces más grande que el promedio anual. Específicamente, la deforestación por dentro y fuera de las AP del COA fue de 37.93 km2, equivalente al 2.2% del área total protegida (0.22 ± 0.36 km2.año-1; Figura 1a, c y e; Tabla S1) y 124.43 km2, 4.9% del área (6.55 ± 4.20 km2.año-1; Figura 1g, i y k ), respectivamente. El área protegida con mayor tasa de deforestación (~0.64 km2.año-1) fue la RFPN Farallones del Citará ubicada en el suroccidente del COA. Es importante resaltar que la deforestación estuvo fuertemente correlacionada con la extensión de las AP (r = 0.63, P < 0.01). Es decir, a mayor área se evidenció un incremento proporcional. Adicionalmente, la deforestación, tanto dentro de las AP como por fuera de ellas se dio en los sitios de mayor riqueza de especies, como lo confirma la prueba de Kolmogórov-Smirnov (P < 0.001) (Figura 1a, e, g, i y k).
Figura 1: Distribución de los Índices Dinámicos de Hábitat en sitios deforestados y fragmentados dentro y fuera de las AP del Corredor del Oso de Anteojos (COA), en Antioquia, Colombia. Las gráficas enseñan los histogramas de distribución de los Índices Dinámicos de Hábitat (IDH1, IDH2, IDH3; barras de color verde), donde a, b, c, d, g, h, i y j corresponden a los IDH1 y IDH2, en los cuales, a medida que aumenta el valor del índice, hay mayor diversidad de especies; y e, f, k y l corresponden al IDH3, donde a mayor valor, menor diversidad de especies. Los valores de cada índice en los sitios deforestados entre 2001 y 2019 se representan en las barras grises, es decir, el valor de los IDH para cada sitio que fue deforestado. Las barras de color blanco representan los sitios con alta fragmentación por fuera y dentro de las AP, es decir, el valor de los IDH en los sitios con alta fragmentación (< 0.7). D es el estadístico de la prueba de Kolmogórov-Smirnov.
En cuanto a la fragmentación, solo el 4.27 % de la cobertura boscosa por fuera de las AP presentó bajo grado de fragmentación (i.e., valores entre 0.8 y 1.0; Figura S1d), principalmente en el centro y noreste del COA y en lugares con baja riqueza de especies. De igual forma, dentro de las AP, solo el 16 % de la cobertura boscosa presentó valores bajos de fragmentación. La fragmentación en el área restante (i.e., 84 %, valores entre 0.0 y 0.7) se dio en los sitios con mayor riqueza de especies (prueba de Kolmogórov-Smirnov P < 0.001; Figura 1b, d, f, h, j, l) lo que favorece al aislamiento de las AP. Por otra parte, encontramos 26 remanentes de bosque sin ningún grado de fragmentación en el COA, de los cuales, el 72.50 % se encuentran en las AP (Figura S1d). Es de resaltar que más del 50 % de los remanentes de bosque tienen pendientes mayores a 20° (de fuertemente inclinadas a escarpadas). Los remanentes poseen áreas entre 99.36 km2 y 0.33 km2 y están distribuidos en su mayoría en el noroeste del COA (57.70 % de ellos) dentro de las AP PNN Las Orquídeas, RFPN de Urrao y RFPN Carauta.
Representación dentro de las AP
Las AP representan el 40.21 % del área total del COA. No obstante, al evaluar la representatividad de las áreas de mayor importancia biológica para cada grupo climático dentro de las AP del COA, según los objetivos de conservación mundial (17 % meta Aichi 12 y 30 %), identificamos que los grupos climáticos no están apropiadamente representados por las AP dentro del COA (i.e., ningún grupo alcanza el 100 % para cumplir con los objetivos de conservación mínimos 17 % y 30 %; Tabla 1). El grupo menos representado es el uno (menos del 15 % en ambos casos), el cual está distribuido hacia la vertiente occidental (Chocó biogeográfico). Los otros dos grupos se encuentran representados en aproximadamente la mitad de sus áreas prioritarias. Aunque el PNN Las Orquídeas no está completamente incluido en el COA, el aporte de este no es significativo para el área de estudio, comparado con otras figuras de conservación. Las reservas forestales de Carauta, Urrao y Farallones de Citará tienen un aporte importante, contribuyendo a la conservación de los grupos climáticos simultáneamente, con valores significativos (Tabla 1). Los DMI tienen mayor impacto en la conservación del 30 % que del 17 %, lo que implica que están localizados hacia áreas de menor importancia que los otros esquemas de conservación. La representatividad de las reservas de la sociedad civil es baja, debido a su poca área, comparada con otras figuras de conservación.
* AP con una extensión considerablemente pequeña que no alcanzan a proteger áreas de importancia biológica de manera representativa.
Tabla 1: Porcentaje de representatividad de las áreas de importancia biológica para cada grupo climático, dentro de las de AP del Corredor del Oso de Anteojos (COA), en Antioquia, Colombia, según los objetivos de conservación mundial (17 % meta Aichi 11 y 30 %). En negrilla están resaltados los valores más altos por grupo y objetivo de conservación. DRMI: Distrito Regional de Manejo Integrado; PNN: Parque Nacional Natural; PNR: Parques Naturales Regionales; RFPN: Reserva Forestal Protectora Nacional, RNSC: Reserva Natural de la Sociedad Civil. El total representa el porcentaje de cumplimiento de la meta propuesta.
Discusión
A partir de información derivada de sensores remotos y grandes bases de datos, logramos realizar una identificación rápida de la distribución de la perdida de cobertura vegetal, fragmentación de los bosques y su relación sobre la diversidad en el COA en Antioquia. El COA, ubicado en una de las últimas regiones relativamente mejor conservadas en el departamento de Antioquia, con un gran auge en la consolidación de diversas organizaciones ambientales locales y regionales (Corporación Gaia, VerdeAgua, Nodo de reservas de la sociedad Civil, Corantioquia, entre otras), presenta una de las coberturas de bosque más altas del departamento (65% COA vs el 58% Antioquia, Morales et al., 2020). Por tanto, se están conservando los bosques al compararlo con otras zonas en el departamento, como en las subregiones del Bajo Cauca, Nordeste y Magdalena Medio, donde se encuentran los municipios con mayor deforestación en Antioquia y donde se ha deforestado entre los años 2000 y 2019 aproximadamente 3000 km2 (Morales et al., 2020), 18 veces más que la deforestación acumulada en el COA (162.37 km2). No obstante, el COA no es ajeno a los motores de transformación del paisaje que enfrenta la región. Si bien presenta una tasa de deforestación relativamente baja para el departamento de Antioquia con solo 4% de pérdida acumulada en los últimos 19 años, los procesos de fragmentación son altos (i.e. ~ 90% de la cobertura boscosa presenta valores entre 0.0 y 0.7) y las coberturas boscosas continuas grandes son escasas. Sumado a que los procesos de fragmentación se concentran en los sitios con mayor importancia biológica (Figura 1 b, f, h, j y l ).
Adicional a la consecuente pérdida de hábitat, la deforestación en sí misma genera una mayor sensibilidad a las perturbaciones al interior de los remanentes de bosque, tales como incendios forestales y tala selectiva, lo que a su vez induce a cambios en la estructura del bosque y la composición de especies (Barlow et al., 2016). Igualmente, aunque no existe un consenso sobre los efectos de la fragmentación sobre la biodiversidad (Haddad et al., 2015; Piquer-Rodríguez et al., 2015), varios estudios han encontrado que ésta afecta negativamente la riqueza de especies, la retención de nutrientes y las dinámicas tróficas de los ecosistemas (Zhai et al., 2014; Haddad et al., 2015; Piquer-Rodríguez et al., 2015). Por tanto, de no mitigarse la fragmentación de los bosques, podríamos enfrentar procesos graves de defaunación en una de las últimas regiones conservadas del departamento.
¿Se están conservando los bosques y las áreas con mayor biodiversidad en el Corredor del Oso de Anteojos en la Cordillera Occidental de Antioquia?
Sorpresivamente, el sistema de AP dentro de la región (~ 40 % del área total del COA) no representa de manera apropiada su importancia biológica. Es decir, el 17 % del territorio con mayor concentración potencial de especies no se encuentra cubierto por el sistema de AP. Al analizar cada uno de los grupos climáticos obtenidos, ninguno de ellos se encuentra representado en las AP, al menos en la mitad de lo esperado (Tabla 1). En particular, uno de los grupos climáticos que contienen la vertiente occidental de la cordillera Occidental en el COA (Chocó biogeográfico; Figura S1b) está representado en las AP en menos del 15 % (2.55 % del territorio). Otros ecosistemas, como los bosques altoandinos y páramos, los cuales se encuentran en los otros dos grupos climáticos, están mejor representados en el sistema de AP.
La biodiversidad del Chocó biogeográfico en Antioquia, vertiente occidental de la Cordillera Occidental, es altamente desconocida. Para esta área, el número de registros reportado en GBIF corresponde en su gran mayoría a aves (probablemente el grupo taxonómico más estudiado) con solo 45 registros, lo cual demuestra el vasto desconocimiento que se tiene de esta región. Sin embargo, es conocido que la Cordillera Occidental, y probablemente su vertiente occidental, está entre las regiones con mayor endemismo del mundo (Kattan et al., 2004; Sandel et al., 2011). Si a esto sumamos la escasa representatividad de esta región en las AP del COA y en general del departamento de Antioquia, se puede concluir que esta zona debería ser una de las principales prioridades de conservación en Antioquia, y merecería ser incluida en la planificación territorial a nivel del COA y departamental. También, algunas áreas climáticamente similares a los páramos en la cuchilla de Jardín-Támesis podrían incluirse como ecosistemas estratégicos, para incrementar la eficacia en su conservación y reforzar el sistema de AP dentro del COA.
Por otro lado, detectamos procesos de desforestación y fragmentación en las AP dentro del COA (Figura 1a, c y e ; Figura S1b). Es importante resaltar que la pérdida y degradación del hábitat converge con las áreas donde los IDH muestran una mayor riqueza de especies potencial (Figura 1a-f), lo cual puede estar relacionado al uso de las áreas con condiciones óptimas para la producción agrícola o ganadera. Estos cambios en el uso del suelo son permitidos en algunas de las estrategias de conservación que involucran desarrollo sostenible, como los DMI. Sin embargo, en AP estrictas, como los PNN, también encontramos registros de deforestación (Tabla S1). Lo anterior, ha sido mostrado a escala nacional, y Clerici et al. (2020) reportaron pérdida de bosques en el sistema de parques nacionales en los últimos años, atribuido a los cambios en los territorios relacionados, entre otros a procesos de acaparamiento de tierra que se han generado posterior a la firma del Acuerdo de Paz de Colombia en los años 2016 y 2017. De igual manera, a nivel departamental y en el COA, observamos un incremento significativo en la deforestación para este periodo (Morales et al., 2020).
Tanto la pérdida de bosque como la fragmentación son mayores en las zonas por fuera de las AP del COA (Figura 1 g-l ; Figura S1b y c). Sin embargo, la fragmentación tiene un mayor impacto en las regiones por fuera de las AP. La fragmentación se midió a través del grado de cohesión de los píxeles de bosque dentro de una ventana de 3 km x 3 km, donde se observa que la mayoría del área del COA no contiene bosques continuos de este tamaño, es decir, gran parte de los bosques fuera del sistema de AP son fragmentos de menos de 1 km2 (100 ha; Figura S1d). En términos de la viabilidad de las especies, se puede pensar que estos pequeños parches de bosque no son importantes, pero se ha demostrado que, en los Andes, estos son claves para mantener relictos de poblaciones de especies endémicas o raras, y además facilitan la conectividad entre grandes áreas de bosque (Wilson & Rhemtulla, 2018). Es por ello que mantenerlos, y preferiblemente incrementarlos, es un asunto prioritario. Adicionalmente, estos fragmentos pueden tener un valor agregado en la agricultura local (Ricketts et al., 2004), manteniendo poblaciones de polinizadores (González-Chaves et al., 2020). Por lo cual, intensificar el estudio de las poblaciones en estos fragmentos y sus efectos en la conectividad de las poblaciones naturales es una línea de investigación clave para el desarrollo territorial sostenible.
Perspectivas y limitaciones
La importancia biológica que implementamos, basada en la combinación de grupos climáticos y los IDH, nos permitió, de manera relativamente fácil, hacer una evaluación preliminar del estado de los bosques y su representatividad el área de estudio, y su dinámica entre 2003 y 2019. Así mismo, nos posibilitó identificar que las áreas de pérdida en este periodo han sido áreas importantes para la conservación, probablemente con un mayor valor que áreas conservadas actualmente. Adicionalmente, la delimitación de los grupos climáticos nos permitió maximizar la cobertura de diferentes ecosistemas dentro del área de estudio de manera sistemática y enfocada al contexto biogeográfico de la misma. Es decir, que al diferenciar en el COA áreas similares climáticamente, nos pudimos enfocar en la cobertura de cada área y realizar un análisis por separado para cada grupo. Aunque esta metodología puede ser poco precisa en generar regiones bióticamente homogéneas, como suponemos, puede ser muy útil en áreas con escasa información de su biodiversidad, como es el caso del COA, o de muchas regiones del trópico en general; no obstante, es importante calibrar estos resultados con estudios de línea base de biodiversidad (Arbeláez-Cortés et al., 2017). Estos grupos climáticos fueron validados mediante la comparación de listas de especies derivadas de GBIF, soportando que los grupos climáticos faciliten la identificación de grupos bióticamente homogéneos. Sin embargo, la validación de nuestros resultados con muestreos en campo y sitios de monitoreo constante nos pueden permitir calibrar de manera apropiada los IDH para futuros estudios de manera más precisa.
Aunque actualmente podemos conocer los procesos de deforestación, idealmente se deberían incluir predicciones futuras del comportamiento de la pérdida de hábitat para tener un escenario más amplio de los cambios que puede tener una región (Morelli et al., 2020). No obstante, los modelos para realizar predicciones de deforestación en escenarios futuros pueden tener diferentes limitaciones (Morelli et al., 2020). Por otra parte, además de medir el efecto de la deforestación y la fragmentación, es importante incluir el efecto que podría tener el cambio climático y sus sinergias sobre los ecosistemas, puesto que solamente el efecto del cambio climático podría reducir drásticamente el hábitat de grandes especies sombrilla, como el Oso de Anteojos.
Hemos incorporado en este análisis dos variables principales para que sean incluidos al sistema de conservación actual del departamento: representatividad de las áreas de mayor importancia biológica en las AP, incluyendo la biodiversidad, y la evaluación de la vulnerabilidad basado en fragmentación y deforestación, puesto que el COA busca proteger una amplia variedad de ecosistemas en los cuales habita el oso de anteojos. Las acciones de protección buscan maximizar la biodiversidad protegida, su integridad (áreas amplias con poca fragmentación) y que estas acciones incluyan las variaciones en diversidad (representatividad). Probablemente hemos excluido análisis minuciosos que requieren corroboración en campo, pero ante la realidad de nuestras regiones, en las cuales los recursos para realizar inventarios y monitores de la biodiversidad son escasos, consideramos que el uso de sistemas de información geográfico puede guiar los esfuerzos para priorizar exploraciones en áreas concretas y así avanzar en nuevas delimitaciones. Por otro lado, si bien nuestros resultados nos señalan la debilidad de las AP en el COA en proteger la cobertura vegetal, no obstante, es la estrategia más efectiva a nivel mundial. Sin embargo, en aquellas AP donde los niveles de vulnerabilidad son más altos, se podría considerar incrementar las restricciones. Ahora bien, como lo señalan Corantioquia y la Corporación Gaia, la conservación de este territorio debe incorporar múltiples opciones de desarrollo sostenible, que incluyan procesos de educación y desarrollo económico integrados, como, por ejemplo, el programa Abrazando Montañas, sumado al acompañamiento de las múltiples iniciativas locales, para garantizar la conservación del territorio de la mano de sus habitantes (Molina & Aristizábal, 2018).
Conclusión
Aunque el área del Corredor del Oso de Anteojos se encuentra en un 40% dentro de AP, éstas no necesariamente fueron diseñadas para maximizar la protección de las especies en la región. Por esto, es una prioridad incluir nuevas áreas, principalmente en la vertiente occidental de esta cordillera (Chocó biogeográfico). También se encontraron pérdidas de bosque dentro de las AP restrictivas, como los parques nacionales, lo cual hace un llamado a los entes administrativos a incrementar las estrategias para reducirla, ya que, si bien la deforestación es relativamente menor en las AP, los procesos de fragmentación son significativos, es decir, que los motores de trasformación permean las AP. Finalmente, se identificó que la fragmentación y pérdida de pequeños parches a lo largo del corredor puede disminuir sustancialmente la conectividad entre los parches y, por ende, contravenir el objetivo principal de esta iniciativa de conservación de mantener las poblaciones de osos de anteojos. En ese sentido, es relevante ejecutar acciones que aumenten la efectividad de las AP y sus alrededores, tales como programas de compensación a comunidades locales, los cuales han demostrado ser efectivos contra la explotación forestal y la deforestación en los bosques tropicales.
