Datos de contacto:

/

logo

Introducción

Monterrey hace parte de los 19 municipios que tiene el departamento de Casanare,limita al norte con Paéz (Boyacá) y Tauramena (Casanare), al sur con Tauramena y Villanueva (Casanare), al oriente con Tauramena y al occidente con Sabanalarga (Casanare). Su área es de aproximadamente 75.900 hectáreas.

Para introducir los usos y cambios de suelo a nivel regional, se efectuó un análisis preliminar del producto ofrecido por The Europian Space Agency (ESA) World cover a una resolución de 10m de las imagenes satélitales de Sentinel 1 y Sentinel 2. Se efectuó una mascara para el departamento de Casanare y posteriormente para llevar el análisis al contexto local, se hizo el mismo ejercicio para el municipio de Monterrey. Los siguientes mapas visualizan el estado de intervención del territorio. Si bien, Monterrey representa tan solo el 1.7% del área total del departamento, por su posición geográfica, en términos de producción hidríca, las alteraciones ecológicas, especialmente la perdida de bosque en las zonas montañosas, repercuten en un desequilibrio del balance hidrológico regional.

Se deriva entonces, datos representativos. De acuerdo con ESA, en el área total del departamento de Casanare, aproximadamente el 22% es cobertura Boscosa. El aporte municipal de Monterrey en términos de bosque esta representado por un 3% del total departamental. Este 3% corresponde a la zona de piedemonte productora de agua para Casanare. Esta es la ecorregión de Bosques subandinos los cuales son estratégicos porque es donde se hace la mayor recarga para la red hídrica del afluente del río Meta, el cual aporta más del 70% del caudal que este río entrega al Orinoco (WFF Colombia, 2016).

Cobertura % area respecto al departamento
Agua 1.42
Bosque 3.31
Pasto 1.49
Vegetación inundable 0.07
Cultivos 0.48
Matorrales 1.43
Zona urbana 6.63
Matorrales 1.64
Nubes 0.81

El municipio presenta dos climas diferenciadores, por un lado en estribaciones de la cordillera oriental presenta un clima medio y muy húmedo comprendido entre las altitudes de 1000 a 2000 msnm. En mayor proporción del gradiente geográfico del municipio se encuentra el clima cálido y húmedo propio de altitudes menores a 1000 msnm.

De acuerdo con la evolución geológica y física el municipio cuenta con tres unidades de paisaje en un gradiente comprendido entre los 188 a 2409 m.s.n.m., siendo estas Montaña, piedemonte y sabana. El municipio a la fecha no reporta actos administrativos para áreas protegidas, distritos de manejo integrado, parques municipales, parques nacionales naturales, reserva forestal natural protectora, reserva municipal, reservas naturales de la sociedad civil, reservas naturales y/o patrimonio ecológico ni reserva forestal hidrogeográfica.

Bosques

En la presente sección se evalúa el estado de la cobertura boscosa de Monterrey - Casanare, esto a partir de información satelital proveniente de Landsat y procesado por la Universidad de Maryland (Estados Unidos), el análisis conlleva el conteo del área boscosa inicial para el año 2000 y temporalmente se visualiza y calcula cómo se ha venido perdiendo / ganando esta cobertura boscosa hasta el año 2020. Esta metodología podría ser extrapolada a cualquier zona del mundo, entre tanto, se tiene disponibilidad para todos los países.

El insumo empleado de acceso gratuito (Global Forest Change) con una resolución aproximada de 0.00027 segundos de arco (~30 mt) en formato espacial (tiff) se validó usando las capas de bosque / no bosque de la metodología diseñada por el IDEAM. Este análisis fue diseñado en su totalidad con el lenguaje de programación de R (las librerías utilizadas fueron raster, rgdal, rgeos, stringr, sf, tidyverse, terra, fs, glue), y usando la IDE RStudio.

Los pasos grosso modo fueron los siguientes:

La capa de cobertura boscosa del IDEAM contiene tres distintas categorías dentro de un mismo archivo espacial raster:

Por su parte para los datos de cobertura boscosa de Global Forest Change se tienen dos distintas capas

  1. Treecover, contiene valores de 0 a 100, indicando el porcentaje del área boscosa de cada píxel.
  2. Loss, contiene dentro de cada valor de píxel el año en el que hubo posible pérdida del área boscosa.

Con el fin de evitar las zonas “Sin información” del IDEAM (estas dadas por el problema del Bandeado de Landsat 7), se utilizó la información de Global Forest Change, y para poder identificar la línea base de bosque

Mientras que los datos de Global Forest Change contiene distintas categorías

  1. Extracción de los valores que toma el raster de treecover para las zonas de bosque del IDEAM

  1. Binarización del raster de treecover para el año 2000 (línea base)

Producto de este paso se tiene el siguiente mapa que ilustra las zonas de bosque (color verde) y las zonas de no bosque (en color gris) para el año 2000.

  1. Conteo de los píxeles en los cuales hubo pérdida de cobertura boscosa, haciendo uso de la capa loss

  1. Cálculo de la pérdida de cobertura boscosa por unidad político administrativa a nivel de vereda 2000 - 2019, en el siguiente gráfico se presentan estos datos.

  1. Generación de salidas gráficas / mapas

De acuerdo a los análisis de cobertura boscosa en el municipio de Monterrey, se encontró que entre el período 2000-2019 se estima una pérdida de bosque de aprox 500 ha en 20 años. De acuerdo a los mapas presentados, es posible derivar dos procesos de alta importancia a nivel espacial. La primera es la visualización de zonas de importancia de protección ambiental por el estado de conservación y, la segunda corresponde a evidenciar aquellas veredas que requieren construir corredores de conectividad debido a la fragilidad de los remanentes de bosque.

La fragmentación del bosque es diferenciada geoespacialmente. En terminos veredales, de las 21 veredas que componen el municipio de Monterrey, tan solo dos presentan cobertura boscosa mayor al 50% (Guayabal y Piñalera), seguida de Guadualito con porcentaje de bosque del 25% y las 18 restantes presentan coberturas de interés ecológico inferiores al 20%, llegando a valores críticos de hasta presencias menores al 5% en la vereda La Vigía.

En general, durante el período evaluado la degradación y perdida de bosque en el municipio esta asociado a picos de deforestación para la primera decena entre 2000 a 2019. Para éste, la vereda Palmira reporta el pico más alto de perdida de bosque, para inicios del período con 80 ha/año de deforestación. En menor proporción se encuentran las veredas Caño Rico y Tacuya con áreas cercanas a las 40 ha/año de pérdida de bosque y el restante con cifras iguales o inferiores a 25 ha año.

Para la segunda decena del período evaluado, la deforestación decrece y permanece constante con valores cercanos a 1 ha/año de deforestación, incluso en algunas veredas como Marenao, La Vigia y Guaneyes se registran valores de no deforestación.

Puntos de calor

Mediante el uso de dos productos disponibles proporcionados por MODIS Collection 6.1 temporal coverage y VIIRS S-NPP (acceso libre en https://firms.modaps.eosdis.nasa.gov/download/) , producidas para el período de 19 años registrados durante de 1 de Noviembre de 2000 al 31 de Dic de 2019. Se agrego la información de las dos capas disponibles para el municipio de Monterrey Casanare y se produjo el conteo asociado a puntos de calor en el terriotio, derivando análisis veredales. Aquí se hizo uso de la librería raster, sf, y tidyverse.

En la siguiente figura se ilustra la ubicación de los incendios ocurridos en el territorio de Monterrey. A la izquierda los registrados por el sensor MODIS y a la derecha los puntos de calor registrados por el sensor VIIRS.

El análisis de conteos de puntos de calor muestra una diferenciación geográfica de los posibles incendios registrados en los 19 años estudiados para las 21 veredas, siendo la vereda Palmira quien reporta un conteo de 44 focos de calor en el terriotrio. seguida de Iguaro con 37 focos de calor, Brisas del Llano con 28 puntos de calor, Bella Vista con 24 puntos, Caño Rico con 23 y Tacuya con 20 puntos de calor.

En menor proporción, pero aún con registros de posibles incendios se encuentran las veredas Palo Negro, Barbasco, Garrabal, la Piñalera y Guayabal con incidencias menores a 16 puntos de calor a nivel municipal.

Las dinámicas de frecuencia de deforestación se pueden explicar en el gráfico participación de puntos de calor que estiman deforestación. La metodología consiste en recrear zonas buffer de 365 metros alrededor de cada foco de calor y que posiblemente fueron afectadas por los incendios ocasionados. Una vez estimadas estas áreas se produce un conteo de los píxeles que fueron afectados en el período de análisis. Entre mayor cercanía geográfica están los puntos de calor, mayor frecuencia de conteo.

Siendo así, se obtuvo que las diez veredas que registran mayores frecuencias de deforestación son Palmira > 25% de estimación de deforestación, Iguaro y Brisas del Llano con porcentajes cercanos al 20% e inferiores a este valor Palo Negro, Barbasco, Tacuya, Caño Rico, Garrabal, Piñalera y Guayabal. Nótese que esta frecuencia no esta asociada directamente a los puntos de calor, sino a los píxeles cercanos que son asociados a la deforestación producida por los focos de incendio.

También fue posible caracterizar la incidencia de los puntos de calor a nivel temporal, tal como se muestra en el gráfico Cantidad de puntos de calor-mes. Para el período evaluado se encuentran un total de 264 focos calientes que se presentaron en su gran mayoría durante la estación de verano o tiempo de baja intensidad de lluvias, siendo el mes de Marzo el que presenta 100 puntos de calor seguidos de los meses de febrero y enero.

Este fenómeno está estrechamente relacionado con el proceso cultural que se lleva a cabo en la región denominado roza, tumba y quema para expansión de la frontera agrícola. Los meses que menos registran puntos de calor son Junio y Julio explicados por los picos altos de pluviosidad que se registran en estos meses. No se evidencia mes que no tenga al menos un registro de punto de calor en el territorio.

Ecosistemas con posible afectación por puntos de calor

En un posterior análisis nos permitió identificar aquellos ecosistemas que han sido modificados por la presencia de zonas de actividad de incendios. De acuerdo a la capa mapa de ecosistemas continentales y costeros de Colombia v. 2.1 escala 1:100.000 del año 2017 generada por el IDEAM y el cruce con los puntos de calor del período en estudio, se obtuvo que cerca del 50% de los incendios han ocurrido en Agroecosistemas, seguido por intervenciones menores al 20% por ecosistemas transicionales o transformados y Bosques. Aunque en menor proporción (<10%), se registran incendios en las rondas hidráulicas de ríos, vegetación secundaria, sabanas, bosque fragmentado y herbazal.

Las dinámicas son heterogéneas en el gradiente municipal. Por ejemplo, Se evidencia que para la vereda con mayor registro de incendios, La Palmira hay focos donde existen las categorías de uso adecuado o sin conflicto y subutilización del suelo. Sin embargo, si nos movemos a la vereda Tacuya y Bella Vista, los puntos de calor están asociados en mayor proporción a Sobreutilización del suelo. Es importante destacar, que la zona Norte del municipio corresponde a la zona de recarga hídrica, la cual tiene uno de los servicios ecológicos principales dentro del municipio. Esta zona presenta un mosaico de puntos que varían desde la sobreutilización hasta el uso adecuado del suelo.

Ecología del paisaje - Métricas del paisaje

Para el cálculo de ecología del paisaje y las métricas del paisaje se hizo uso de las librerias: raster, rgdal, rgeos, tidyverse y quizá la más importante es landscapemetrics, esta última librería se utilizó para el cálculo de las métricas del paisaje, tanto para las zonas de bosque como para las coberturas del paisaje que componen el municipio de Monterrey. La fuente de las capas de cobertura de la tierra son del IDEAM, descargadas desde su portal web para los años 2000, 2005, 2010 y 2018.

Se hicieron uso de 18 métricas del paisaje para analizar los patrones de distribución del paisaje en Monterrey con características de zona rural con asentamientos urbanos dispersos. Es decir, el análisis del paisaje se efectuó para determinar el grado de alteración antrópica dentro del territorio. Las métricas del paisaje utilizadas se pueden caracterizar en dos grandes grupos, por clase y por parche. La clase refiere a los tipos de ecosistemas presentes en el área administrativa municipal y los de parche corresponde a las características que presentan los grupos de parche dentro de cada clase presente.

En la siguiente tabla se ilustra el listado de las métricas realizadas.
Clase Parche
Índice de agregación Área del parche
Área media del parche Índice de continuidad
Radio de giro (área y borde métrico) Coeficiente de variación de la distancia euclidiana al vecino más cercano (métrica de agregación)
Índice de aglomeración (métrica de agregación) Índice de dimensión fractal del coeficiente de variación (métrica de forma)
Porcentaje de área central del paisaje (CPLAND) Relación perímetro-área (forma métrica)
Coeficiente de variación del área del parche Radio de giro (área y borde métrico)
Índice de parche más grande (métrica de área y borde) Índice de forma (métrica de forma)
Densidad del área central disjunta (métrica del área central)
Porcentaje de adyacencias similares (métrica de agregación)
Dimensión fractal de área perimetral (métrica de forma)
Borde total (clase) (métrica de área y borde)

En la siguiente figura se ilustran las distintas coberturas para los años mencionados anteriormente.

Índice de agregación

En unidades porcentuales muestra para los cuatro años estudiados y las 12 clases de cobertura de suelo analizadas valores cercanos a 100%, en consecuencia, en el mosaico geográfico se puede decir que las clases están altamente agregadas.

Área media del parche

En unidades de hectáreas, esta métrica indica cuáles son las clases de mayor presencia dentro del municipio. Para los años 2000, 2005 y 2010 la mayor representación estaba dada por el mosaico de espacios naturales y cultivos, sin embargo, el año 2015 muestra una fuerte disminución de esta cobertura y alcanza valores mayores la cobertura boscosa, bosque ripario y herbazal.

Cabe destacar, que para el año 2018, la entidad cultivos fue diferenciada a mayor especificidad en la capa del IDEAM y por esto, también sugiere que la clase de mosaicos naturales y cultivos esta desagregada. Es importante destacar una clase importante en términos agrícolas, los pastos también presentan una disminución en el 2018 respecto a los años predecesores, y aparecen nuevas clases en mayor proporción: cultivos, herbazales y piscícolas.

Índice de aglomeración (métrica de agregación)

Este índice sin unidades, permite estimar que tan agregados están los parches dentro del paisaje. La gráfica muestra que en general todas las clases tienen un comportamiento cercano a 1 lo que indica que las clases dentro del paisaje tienden a estar agregadas. Si fuesen valores cercanos a -1, el paisaje sería altamente desagregado y si estuviesen cercanos a 0 esto indicaría una conformación de clases aleatoria dentro del paisaje.

Porcentaje de área central del paisaje

Esta medida expresada en porcentaje, indica que las clases de vegetación secundaria, bosque denso y ripario, tejido urbano, plantaciones forestales, suelo desnudo, ríos, herbazales, cultivos y piscícolas están distribuidos de manera simple dentro del paisaje. Parches pequeños representados por valores porcentuales inferiores a 10%, mientras que pastos y mosaicos naturales con pastos y cultivos con valores superiores al 35% y 15% respectivamente y a lo largo de la serie de tiempo del análisis, son más grandes en términos de forma dentro del paisaje.

Coeficiente de variación del área del parche

Unidad que representa métrica tanto de área como borde del parche, expresado en hectáreas, indica que existe una alta variación de tamaño de todos los parches que conforman las 12 clases del mosaico del paisaje a excepción del suelo desnudo y las plantaciones forestales. Las clases que presentan picos más altos de variación en comparación con los años evaluados entre una misma clase, son heterogéneas. Esto indica, por ejemplo, para el caso de los ríos que las zonas de aluvión han variado en tamaño y forma para el año 2018 respectos a los años predecesores. Otro hallazgo similar es la clase herbazal, que por las dinámicas de intervención varia respecto a la serie de tiempo trabajada. Los bosques presentan una variación positiva, ya que incrementan las hectáreas positivamente para el final de la serie, es decir, para el año 2018. ## Cambio climático

Coeficiente de variación de la distancia euclidiana al vecino más cercano (métrica de agregación)


Este coeficiente expresado en metros pero para facilidad de gráfico lo simplificamos a unidades en kilómetros, muestra la distancia que puede existir de un parche a su vecino de la misma clase en el paisaje analizado. En general para la serie de tiempo, las distancias entre parches cercanos son inferiores a 1 km, pero para el año 2018 las distancias entre parches de la misma clase se distancian más, llegando a valores de 4 km. La clase de cobertura denominada tejido urbano, muestra que para inicios del período de estudio los parches se encontraban más distantes entre sí; en 2018 esta distancia euclidiana se disminuye, lo cual puede ser explicado por una dinámica de urbanización más acelerada.

Índice de parche más grande (métrica de área y borde)


Sin lugar a dudas esta medida porcentual es una de las métricas más importantes para describir el paisaje. Nos permite saber que comportamiento porcentual ha sufrido cada cobertura del suelo a lo largo del tiempo. Los pastos muestran menor fragmentación y por lo tanto, parches más grandes dentro de la unidad de paisaje, sin embargo, para el año final de estudio 2018 presenta una disminución mayor a tres puntos porcentuales.

Las dinámicas agrícolas asociadas con cultivos son heterogéneas, disminuyen al inicio del período, luego hay un momento constante y posteriormente incrementan al finalizar 2018. Esto indica que hay parches pequeños y discontinuos dentro del paisaje. En términos negativos, no se evidencia recuperación de bosque denso al final de la serie de tiempo, pero si un ligero incremento para las zonas de ronda hidráulica.

En general no hay una marcada tendencia de dominancia de las clases dentro de la unidad de paisaje, pero los parches que tienden a ser cada vez más pequeños son las plantaciones y los suelos desnudos. No hay evidencia de contrastante para la clase piscicola ya que solo aparece para la capa del año 2018.

Densidad del área central disjunta (métrica del área central)


Métrica expresada en Número de parches por cada 100 hectáreas, derivadas del área central de cada parche. Para la serie de tiempo en general, los parches se comportan de una manera más compleja o fragmentada y dispersa dentro de la unidad de paisaje. Casos particulares que pueden ser descritos son los pastos, bosque ripario, ríos y herbazales que para el año 2018 tienden a comportarse más grandes respecto a los otros años comparativos.

Porcentaje de adyacencias similares (métrica de agregación)


En unidades porcentuales es importante destacar que las 12 clases muestran que los píxeles de cada parche están acompañados de parches similares, es decir, están agregados ya que los valores tienden a ser cercanos a 100%. Sin embargo, desde el 2010 la mayoría de las clases disminuyen la forma de agregación, indicando mayor intervención dentro de las coberturas del suelo.

Dimensión fractal de área perimetral (métrica de forma)


Métrica adimensional que varía entre 1 y 2 que indica la forma de los parches en términos de la unidad de paisaje. Para la serie de tiempo, lo que se puede interpretar es que el mosaico de parches fluctúa entre formas de parches simples, es decir, valores cercanos a 1 y, en menor proporción a formas complejas que es el comportamiento natural de los ríos.

Borde total (clase) (métrica de área y borde)

Es una medida expresada en metros que muestra para toda la serie de tiempo y las 12 coberturas valores mayores a 1 y sin límite, indicando que el paisaje en general esta fragmentado.

Área del parche


En hectáreas la métrica muestra las áreas de parche por clases presentes en la serie de estudio. Puede incrementar sin limite de acuerdo al tamaño de los parches. Para el paisaje del municipio de Monterrey, durante los primeros tres años evaluados el mosaico de espacios naturales, pastos y cultivos predomina con parches que alcanzan hasta las 350 ha, sin embargo, al finalizar el período decrece considerablemente hasta alcanzar tamaños de parche de hasta 150 ha. Lo que destaca este cambio es la variación de la clase en mención por la clase herbazales, que suponen zonas de descanso de pastos e inician los procesos de sucesión secundaria. Una ventaja a lo largo del tiempo, es que se denota una ligera recuperación de área de bosque denso y ripario en este mosaico que se encuentra altamente intervenido y fragmentado. Por ejemplo, para el bosque denso la media de área de parche esta sobre 50 ha y la de bosque ripario cada se fragmenta más ya que inicialmente se observaba medias que alcanzaban las 70 ha y al finalizar el período evaluado se registran parches menores a 50 ha. En términos agrícolas, la clase pastos tiende a mostrar parches de menor tamaño, incluso con medias menores a los 5 ha, cuando para el año 2000 la media de área podría ser de hasta 30 ha. Estas dinámicas han ido cambiando por tres fenómenos diferentes, la construcción de un tejido urbano más denso, la aparición de cultivos que pueden presentar picos extremos mayores a 150 ha y para el final de período de estudio, la aparición de piscícolas con áreas de intervención superiores a la 50 ha.


Índice de continuidad


Este índice adimensional expresa valores cercanos a cero cuando los parches en el paisaje están conectados a parches de una categoría de clase específica. La gráfica indica que hay valores extremos cercanos a cero para el primer cuartil de la serie, especialmente para aquellas clases en el paisaje que son más predominantes, como pastos, bosques densos, bosques de galería, ríos y herbazales. Esto quiere decir, que hay un porcentaje relativamente alto en el paisaje (al menos 25%) para estas clases que no están conectadas. Para clases con parches de menor área, la conectividad del vecino es más elevada, es decir, es posible que para estas clases como mosaicos y tejidos urbanos haya una tendencia a expandirse o estar muy cercanos entre ellos.

Índice de dimensión fractal del coeficiente de variación (métrica de forma)


Siendo adimensional, esta métrica muestra la complejidad que puede presentar cada parche dentro de una misma clase, siendo el valor 0 la forma menos compleja dentro del paisaje. En general, la variación que presenta el grupo de parches seccionado por clases es superior a 1, e incrementa sin limite a valores extremos dentro del grupo de datos ligeramente. En general, todos los valores expresan una ligera complejidad dentro del paisaje de Monterrey, esto puede entenderse debido a la variación de tamaño del parche dentro de cada clase, pero variaciones provocadas por la intervención cultural de las coberturas del suelo.

Relación perímetro-área (forma métrica)


Esta métrica describe la forma del parche, buscando identificar la complejidad de las formas presentes en el paisaje. Lo que se obtiene en el gráfico es el cambio del radio directamente relacionado con el tamaño del parche e independiente de su forma. En general los valores incluso extremos de los años evaluados sugieren formas poco complejas que varían entre 0.001 y 0.15.


Radio de giro (área y borde métrico)


En general, para el histórico de años y clases analizadas, el radio de giro muestra que a pesar de tener dentro del mosaico áreas de parche grandes, no necesariamente indica que están compactados ya que las unidades de la métrica superan grandemente los valores a 0. Si los parches dentro del paisaje lo estuvieran, la métrica sería cercana a 0 metros.


Índice de forma (métrica de forma)


La descripción del índice adimensional indica la descripción del radio entre el actual perímetro del parche y el perímetro mínimo hipotético. Así, para valores equivalentes a 1, presentes en siete de las doce clases presentes en toda la serie de tiempo, las formas de parche tienden a ser cuadradas. Sin embargo, cuando se evalúa el diagrama de caja, desde el segundo cuartil se puede identificar para todas las clases una tendencia dentro del conjunto de datos a ser cada vez más complejas las formas dentro del mosaico del paisaje.


Porcentaje de área central del paisaje


Esta medida expresada en porcentaje, indica que las clases de vegetación secundaria, bosque denso y ripario, tejido urbano, plantaciones forestales, suelo desnudo, ríos, herbazales, cultivos y piscícolas están distribuidos de manera simple dentro del paisaje. Parches pequeños representados por valores porcentuales inferiores a 10%, mientras que pastos y mosaicos naturales con pastos y cultivos con valores superiores al 35% y 15% respectivamente y a lo largo de la serie de tiempo del análisis, son más grandes en términos de forma dentro del paisaje.