Spanish (pdf)
Article in xml format
Article references
Automatic translation
Send this article by e-mail
Cited by SciELO 
Cited by Google 
Similars in
SciELO 
Similars in Google 
Permalink
1. Introducción
El paisaje, según Bertrand (1968), es una porción del espacio caracterizado por un tipo de combinación dinámica de elementos geográficos diferenciados, abióticos, bióticos y antrópicos; es también altamente dinámico y mayormente cambiante a medida que se aleja de un punto de equilibrio en términos económicos, productivos o ambientales. Frente a lo anterior, resulta claro que, para el entendimiento de la dinámica del paisaje, los estudios deben enmarcarse dentro de un contexto temporal, que permita evaluar y entender sus niveles de transformación. Las transformaciones del paisaje, según Aldana (2011) se refieren a la variación que el paisaje puede sufrir en su estructura, composición y funcionamiento durante un periodo de tiempo determinado, el cual está dado por el cambio de estado de una o más variables; estos cambios pueden darse por múltiples factores. Al respecto, Gómez-Zotano y Riesco Chueca (2010), manifiestan que dichos cambios se agudizan durante un número reducido de periodos cortos, separados por largos periodos de estabilización.
Uno de los paisajes que a nivel mundial resulta importante explorar y analizar por sus funciones ecológicas y por su importancia social son los páramos, los cuales a pesar de estar inmersos en modelos de protección como son los parques nacionales, las gestiones políticas y administrativas quizá no han sido lo suficientemente eficaces para detener la intervención antrópica; lo cual implicado unas transformaciones constantes, y con ello repercusiones negativas ambientales y sociales. Un ejemplo de estos páramos es el de Paja Blanca, ubicado en el departamento de Nariño, Colombia, decretado como área protegida desde el año 1989 y como Parque Natural Regional en el año 2015.
Corponariño (2010) afirma que en el páramo de Paja Blanca son 113 hectáreas que se deforestan anualmente y que a pesar de que esta entidad ha realizado tres planes de manejo diferentes, y ha venido implementando desde 1997 diversos mecanismos encaminados a detener el avance del sector agropecuario, como son proyectos productivos pecuarios, reforestaciones con especies nativas, fomento de estufas a gas, letrinización, entre otros, los resultados no han sido los mejores ni contundentes; los procesos de transformación del paisaje, reflejados en los cambios de coberturas, aún continúan y con ello las alteraciones ambientales, como es la reducción de los caudales que alimentan los 45 acueductos y los cinco distritos de riego de los siete municipios que lo rodean.
En la actualidad, se han desarrollado varios métodos de análisis y modelamiento espacial para explicar y predecir los cambios que se producen o se producirán en un determinado paisaje dependiendo del comportamiento de las coberturas en el tiempo y espacio. Al respecto, Sandoval y Oyarzun (2004) manifiestan que la modelación espacial del cambio en el uso del suelo se deriva de la combinación y uso de técnicas cartográficas, Sistemas de Información Geográfica y modelos estadísticos multivariantes. En este sentido, en el presente estudio se usaron metodologías de percepción remota y modelamiento espacial para explicar las transformaciones del paisaje del páramo de Paja Blanca en los últimos 30 años, así como su prospectiva para los 30 próximos, de continuar con las tendencias actuales; información que será muy útil para los procesos de toma de decisiones de las entidades encargadas de su administración y manejo.
2. Área de estudio
La zona de estudio corresponde al páramo de Paja Blanca, que se ubica en el departamento de Nariño, en la parte central del nudo de los Pastos entre las coordenadas 0° 50'-1° 06' de latitud norte y 77° 45'-77° 25' de longitud oeste (Figura 1). Según el documento elaborado por Corponariño en el año 2010, denominado "Declaratoria del Parque Natural Regional Páramo de Paja Blanca, Territorio Sagrado del pueblo de los Pastos", tiene un área de 4640 hectáreas, superficie tomada en cuenta para este estudio. Presenta una distribución altitudinal entre los 3000 y 3600 m s. n. m. La cumbre del páramo de Paja Blanca corresponde a un cráter antiguo con productos volcánicos extensos. Según Corponariño (2010), corresponde a un sistema aislado del complejo de páramos suroccidentales del departamento de Nariño; tal situación y la altura de los páramos circunvecinos (Chiles, Cumbal y Azufral) han condicionado una estructura ecosistémica de características atípica-mente más xéricas, con menos influencia húmeda o perhúmeda de las corrientes pacíficas y amazónicas respecto a los otros páramos suroccidentales. Las temperaturas medias anuales oscilan entre los 10,2 °C y 11,6 °C, y las precipitaciones entre 1094 y 1242 mm. (Corponariño, 2007a). Este páramo cubre los sectores altos de los municipios de Pupiales, Sapuyes, Ospina Iles, Guachucal, Contadero y Gualmatán. En relación con los eventos volcánicos, ha recibido, por su relativa cercanía con el sistema del volcán Azufral y en diferentes períodos geológicos, grandes cantidades de material de tipo volcánico, que en la actualidad se puede hacer evidente en los cortes y capas deposicionales presentes en la geología del páramo Paja Blanca, extendiéndose sobre todo hacia el sector nororiente, en Sapuyes (Corponariño, 2007b).
3. Metodología
3.1. Procesamiento digital de imágenes
Para establecer las coberturas del año 1984, se utilizaron imágenes satelitales Landsat 7 TM (Path 010, Row 059). Para el año 2013, se adquirieron imágenes del sensor Aster mediante el aplicativo de la NASA disponible en internet (https://wist.echo.nasa.gov/). Para la calibración geométrica de las imágenes, se utilizó el programa Erdas 2011, en el cual mediante la recolección de puntos GPS en campo con el actual Datum de la cartografía colombiana (MAGNA-Sirgas), se realizó la georreferenciación (RMS < 1). La calibración radiométrica y atmosférica, proceso que tiene como objetivo corregir problemas mecánicos en el sensor que generan valores erróneos en píxeles concretos y eliminar las distorsiones que la atmósfera introduce en los valores de radiancias que llegan al sensor desde la superficie terrestre, implicó la conversión de niveles digitales (ND) a radiancias (Chuvieco, 1996); conversión de niveles de radiancia a reflectividad aparente o reflectancia (Wooster, 1996; Thome, Palluconi, Takashima y Masuda, 1998; Chander y Markham, 2003). La calibración atmosférica se hizo a través del modelo propuesto por Chávez (1996), el cual consiste en detectar en una banda el valor mínimo (m), con el fin de sustraer a cada ND de la imagen dicho valor.
3.2. Clasificación de las coberturas
Después de los anteriores procedimientos, en el programa Erdas 2011 se procedió a realizar una clasificación de las coberturas presentes en cada año; para ello se aplicó la metodología de clasificación supervisada, la cual exige una misma leyenda para las dos fechas, con el propósito de que sean realmente comparables. La clasificación de las categorías de uso se efectuó teniendo en cuenta la metodología Corine Land Cover; y todo el proceso de vectorización así como la producción final de los mapas de coberturas se realizaron en el programa ArcGIS 10. Las categorías utilizadas fueron: bosque natural denso (BND), bosque natural fragmentado (BNF), mosaico de pastos y cultivos (MPC) y coberturas de páramo (PA).
3.3. Dinámica de los cambios de coberturas del suelo ocurridos entre 1984 y 2013
Una vez obtenidos los mapas de usos y coberturas, se realizó un análisis estadístico descriptivo espacio-temporal de los cambios ocurridos entre 1984 y 2013, utilizando para ello el módulo Land Change Modeler (LCM) del programa IDRISI ANDES; este módulo permite realizar de manera rápida un análisis de las ganancias y pérdidas de superficies experimentadas por los distintos usos y coberturas de suelo, sus cifras netas y transiciones; datos que se usaron para construir una matriz de tabulación cruzada, que se describe como una tabla con arreglos simétricos que contiene en uno de sus ejes las categorías de uso de suelo y vegetación del año base y en el otro eje las mismas categorías pero del segundo tiempo.
3.4. Identificación y análisis de los principales factores que han incidido en los procesos de transformación del paisaje
En el presente estudio se optó por tomar las causas más relevantes que propone la literatura basadas en la caracterización de los diferentes flujos de cambio en la cobertura y el uso del suelo. De esta manera, se analizaron los factores demográficos como principales causantes de las transformaciones paisajísticas (Henríquez-Dole, 2012; Bonilla, 1985; Harthshorn, 1983; Pérez y Protti, 1978). Se tomaron datos estadísticos del Departamento Administrativo Nacional de Estadística (DANE), para con ello analizar las tendencias de crecimiento en los siete municipios que conforman el ámbito de estudio; de igual forma, los datos de crecimiento poblacional se utilizaron para generar mapas de densidad poblacional. Con el objetivo de averiguar si dichos valores de densidad poblacional predisponen una mayor ocurrencia de unas categorías temáticas de ocupación del suelo, antes que otras, se utilizó el estadístico Cramer's V., cuyos valores varían entre 0 y 1; un valor de 0 implica independencia de los atributos y un valor de 1 denota asociación perfecta.
Considerando las recomendaciones de los anteriores autores, se realizaron también distancias euclidianas y regresiones logísticas con el fin de descubrir la relación entre la estructura vial y las viviendas, con las áreas que se transformaron durante 1984 y 2013. Los mapas de distancias euclidianas, que corresponde a la distancia desde las vías y viviendas a cualquier punto de las zonas de estudio, se construyeron con el comando "Distance" del programa Idrisis Selva. La regresión logística es un método similar a la regresión lineal múltiple, pero adaptado a modelos en donde es condición necesaria que la variable a explicar sea binaria, pudiendo tomar únicamente valores de 0 y 1 (Pineda-Jaimes, Bosque-Sendra, Gómez-Delgado y Plata-Rocha, 2009). El estadístico utilizado para evaluar el modelo de regresión fue el ROC (Relative Operating Charasteric), sus valores oscilan entre 0 y 1, de modo que valores de 0,5 o menores no indican ningún grado de correlación y valores cercanos a 1 indican un buen ajuste entre variables; todo el proceso se realizó con el módulo Logisticreg del programa Idrisi Selva.
Al igual que las variables socioeconómicas, también se hicieron regresiones logísticas para encontrar relaciones entre algunas variables físico-bióticas y las transformaciones presentadas en los últimos 30 años. Las variables que se escogieron fueron: el relieve, las pendientes y la distancia a los ríos. Para el relieve, se utilizó un modelo digital de terreno de 12,5 metros de resolución, proveniente del sensor ALOS PALSAR, las pendientes se obtuvieron del modelo digital de terreno calculadas en porcentaje con el comando "SLOPE" del programa Idrisi Selva, y para la distancia a los ríos se usaron distancias euclidianas.
3.5. Simulación
La simulación de los paisajes del páramo de Paja Blanca para los próximos 30 años se realizó con base en el enfoque de Redes Neuronales Artificiales (RNA) (Isasi y Galván, 2004), las cuales son estructuras de mapeo no lineal basadas en el funcionamiento del cerebro humano. Se utilizó para ello el módulo Land Change Modeler (LCM) que viene incorporado en el programa Idrisi Selva. Dicho enfoque requirió efectuar los siguientes pasos:
Selección y análisis de variables de transición: con base en los resultados del acápite anterior, se seleccionaron las variables cuyo patrón espacial manifieste algún tipo de comportamiento que se crea que puede tener relación con el patrón de localización de cada uno de los cambios de coberturas que se identificaron. Las variables seleccionadas fueron sometidas a la prueba Cramer's V para calcular la fuerza de la relación entre estas y la distribución de las categorías de usos del suelo de la fecha más reciente (2013). Si bien un alto índice Cramer's V no garantiza una buena modelización, un bajo índice sí es un indicativo de la conveniencia de no incluirlas en la modelación (Clarks Labs, 2006). LCM permite realizar de manera muy fácil la prueba de Cramer's V, para tomar la decisión de incluir o no una variable en el modelo.
Identificación de las transiciones potenciales: en este paso, se definieron en primera instancia cuáles son las transiciones o cambios más importantes, para lo cual se escogieron aquellas que afectarían a más de dos hectáreas de superficie. Con las transiciones definidas, en el módulo LCM de Idrisi Selva fueron agrupadas en un submodelo, al cual se agregaron las variables de transición, para posteriormente generar mapas potenciales de transición, que son la plasmación espacio-temporal del potencial de cambio que oscilan de 0 a 1, indicando el potencial que cada píxel presenta para cambiar de una categoría de cobertura a otra. Este proceso se basa en el entrenamiento de las redes neuronales a través de la minimización del error cuadrático medio, variando la tasa de aprendizaje, el número de nodos en la capa intermedia, el número total de iteraciones y la revisión gráfica de la evolución del error cuadrático medio.
Estimación de la cantidad de cambio y localización de los cambios: después de obtener las matrices de transición de cada categoría de cobertura, el módulo LCM de Idrisi predice los cambios para los próximos 30 años (2042). Este modelo, en primer lugar, genera una matriz de probabilidad de cambio para el periodo asignado, luego ofrece dos clases de predicción: una imagen con la predicción de un escenario específico y una imagen con índices de la vulnerabilidad al cambio. La primera se presenta con la clasificación final basado en un modelo de localización competitiva del suelo, y la segunda representa las áreas que presentan las mejores condiciones para cambiar en una escala de 0 a 1 (Henríquez-Dole, 2012).
3.6. Validación
Para validar el modelo obtenido, se hizo una simulación para el año 2016, obteniéndose un mapa raster de coberturas para dicho año; posteriormente se realizó en campo un muestreo aleatorio simple, en el que se seleccionaron al azar 50 puntos donde se registró su posición GPS y el tipo de cobertura presente. Estos puntos fueron ingresados al programa ArcGIS 10 con sus respectivas coordenadas geográficas y sus atributos de coberturas observadas; posteriormente se ingresó el mapa raster de las coberturas simuladas para el año 2016, y se usó la herramienta "Extract Values to Points" de ArcGIS para capturar las coberturas simuladas. Como resultado se obtiene una tabla donde cada punto contiene el tipo de cobertura observada y simulada. Con estos puntos se construyó una matriz de confusión, que, según Congalton, Lunetta, Fenstermaker, Jensen, McWire y Tinney (1991), consiste en la colocación de filas y columnas relacionadas con las categorías de uso o cobertura relacionada con el mapa. En las columnas se ubican las categorías deducidas por la simulación, mientras que en las filas se colocan las clases de referencia (verdad del terreno). El error asociado a cada columna se denomina "error de inclusión o comisión" y representa aquellas áreas que fueron asignadas a dicha categorías sin pertenecer a ella. En las filas se presentan las categorías verdaderas. El error asociado a cada fila se denomina "error de exclusión u omisión" y representa aquellas áreas que no fueron asignadas a dicha categoría aun cuando pertenecían a ella. También se evaluó la clasificación con el índice de Kappa, que es otra forma altamente difundida de medir la exactitud de un mapa. Conocida también como "k", puede medir la exactitud de manera más precisa que la matriz de confusión, porque incluye dentro del cálculo todos los valores de la matriz y no solamente sus extremos; su fórmula según Cohen (1960) es la siguiente:
Donde:
i = dimensión de la matriz (número de clases);
aii = número de observaciones en la línea i, columna i;
ai. y a.i = total marginal de línea i y de columna i,
n = número total de observaciones
4. Resultados y discusión
4.1 Dinámica de los cambios de coberturas del suelo ocurridos entre 1984 y 2013
La Figura 2 muestra las coberturas del suelo que se obtuvieron mediante el procesamiento de las imágenes Landsat y Aster para los años 1984 y 2013.
Fuente: Elaboración propia.
Figura 2 Coberturas del suelo en el páramo de Paja Blanca: años 1984 y 2013
La Tabla 1 muestra la cantidad de hectáreas que se ganaron y se perdieron en cada categoría de cobertura del suelo durante el periodo 1984-2013.
Tabla 1 Tabulación cruzada con ganancia, pérdidas y persistencia de coberturas del suelo en el Páramo Paja Blanca para el periodo 1984-2013
| 1984/2013 | BND | BNF | MPC | PA | Total | Pérdidas |
|---|---|---|---|---|---|---|
| Bosque natural denso (BND) | 18,5 | 0,0 | 3,8 | 0 | 22,3 | 3,8 |
| Bosque natural fragmentado (BNF) | 12,3 | 1641,2 | 348,8 | 33,2 | 2035,5 | 394,3 |
| Mosaicos de pastos y cultivos (MPC) | 0,2 | 292,3 | 1350,0 | 0 | 1642,5 | 292,5 |
| Páramo (PA) | 0 | 77,8 | 0 | 862,1 | 939,9 | 77,8 |
| Total | 31,0 | 2011,3 | 1702,6 | 895,3 | 4640,2 | |
| Ganancias | 12,5 | 370,1 | 352,6 | 33,2 |
Fuente: Elaboración propia.
Los resultados muestran que en el periodo 1984-2013 fueron 768 ha que se transformaron entre las diferentes categorías. Las transformaciones más críticas corresponden al bosque natural fragmentado (BNF), puesto que en dicho periodo se perdieron 394 hectáreas que pasaron a ser dedicadas en su mayoría (348 has) a mosaico de pastos y cultivos (MPC). Cabrera y Ramírez (2014) sostienen al respecto que el incremento de las actividades antrópicas en los páramos ha sido reportado desde mucho tiempo atrás, y que en la actualidad los MPC es el tipo de cobertura de uso antrópico que más predominan en los páramos colombianos; esto se debe, según Morales-Betancourt y Estévez (2006), a que, a diferencia de las selvas húmedas, los páramos pueden ser adaptados fácilmente para el cultivo y la ganadería con solo la quema de predios en los que el fuego se esparce fácilmente por acción de los fuertes vientos.
Hubo también 292 hectáreas de MPC que mediante procesos de protección y regeneración natural pasaron a ser en el año 2013 bosque natural fragmentado (BNF).
En bosque natural denso (BND), se presentaron más ganancias que pérdidas, que desde el punto de vista ambiental es bastante positivo, reflejando quizá que las gestiones de manejo del páramo por parte de Corponariño han logrado que las comunidades preserven los bosques fragmentados hasta el punto de que lleguen a un estado de BND.
Las coberturas de páramo (PA), ya sean herbáceas o arbustivas, sufrieron pocas transformaciones, presentándose casi 78 hectáreas de páramo que pasaron a ser BNF.
Cabrera y Ramírez (2014), como producto del resultado de la actualización del Atlas de páramos de Colombia, establecen unas categorías según el grado de intervención antrópica que han sufrido; de acuerdo con esto, Paja Blanca presenta un grado de transformación baja, ya que las coberturas transformadas antrópicamente representan el 7,5 % del total del páramo; siendo similar a otros páramos colombianos clasificados dentro de esta categoría y que presentan porcentajes casi similares como son: Guanacas-Puracé-Coconucos (7,34 %), Chilí-Barragán (6,89 %) y Chingaza (6,72 %). Según estos autores, en estos páramos los esfuerzos de restauración deberían enfocarse a la contención y reconversión de la transformación del paisaje que contribuya a disminuir los procesos de fragmentación.
4.2 Principales factores que han incidido en los procesos de transformación del paisaje en el páramo de Paja Blanca
En el páramo de Paja Blanca, las comunidades interactúan directamente con los recursos naturales que ofrece este ecosistema, para solventar necesidades alimenticias y económicas; por ello, Corponariño (2007a) afirma que, en la zona de estudio, el crecimiento poblacional ha sido un factor determinante en las transformaciones del paisaje. La Figura 3 indica cómo ha sido la tendencia del crecimiento poblacional de los siete municipios que conforman el páramo de Paja Blanca desde el año 1985.
Fuente: DANE (2013).
Figura 3 Crecimiento poblacional de los siete municipios que integran el páramo de Paja Blanca 1985-2013
Según el censo poblacional de 1985, el total de población en los municipios del páramo de Paja Blanca fue de 58 788 habitantes, representando una densidad poblacional promedio de 127 habitantes/ km2, que es bastante superior a la densidad media del departamento de Nariño (45 hab/km2), y mucho mayor a la media nacional (36.3 hab/km2). Los municipios de Guachucal y Pupiales presentan la mayor población, doblando su cantidad al resto de municipios, sin embargo, la mayor densidad poblacional se presenta en los municipios de Iles y Sapuyes, con 155 y 190 hab/km2, respectivamente. Para el posterior censo, en el año 2005, el total de población pasó a 71 111 personas, es decir, hubo un aumento del 20 % (12 323 habitantes), y para el 2013, el DANE proyecta una población total de 71 519 habitantes.
El área de influencia directa del páramo Paja Blanca se compone de 27 veredas, con una población total de 6950 personas, que corresponde al 13,08 % del total de la población rural. Según Corponariño (2007a), sobre los 3200 m. s. n. m. habitan 305 familias, con un total de 1400 personas aproximadamente, que representa el 2,4 % del total de la población rural. De acuerdo con Corponariño (2010), el 56 % de la población que habita en el interior del área del páramo tiene 30 años o menos, y su distribución por género es de 51,06 % hombres y 48,94 % mujeres.
Posiblemente la presencia de gente en las partes altas del páramo obedezca a una situación socioeconómica relacionada con la tenencia de la tierra; las partes bajas ya han sido lo suficientemente explotadas, lo que ha obligado a algunos campesinos a buscar tierras más productivas hacia zonas de más altura, donde las condiciones edafológicas son más favorables. Al respecto, Corponariño (2007a) concluye que "la totalidad del páramo de Paja Blanca está ocupado, no existen baldíos, sus propietarios en más del 70 % habitan en el páramo y por ende de acuerdo a las condiciones económicas, sociales y culturales deben su sustento de lo que producen en el mismo" (p. 154).
Los resultados con respecto a si los valores de densidad poblacional de cada uno de los siete municipios que hacen parte de Paja Blanca predisponen una mayor ocurrencia de unas categorías temáticas de ocupación del suelo, antes que otras, no mostraron relación alguna; se presentaron coeficientes de Cramer´s V bastante bajos, lo que indica independencia entre estas dos variables. Las regresiones logísticas entre las vías y viviendas (Figura 4) y su relación con las áreas transformadas (Figura 5) mostraron un ROC bastante alto (Tabla 2).
Fuente: Elaboración propia.
Figura 4 Distancias euclidianas de viviendas (a) y vías (b) en el páramo de Paja Blanca
Fuente: Elaboración propia.
Figura 5 Áreas transformadas o de cambio entre 1984 y 2013 en el páramo de Paja Blanca
Tabla 2 Valores de regresión logística entre áreas transformadas y distancia a vías y viviendas en el Páramo de Paja Blanca
Fuente: Elaboración propia.
Los resultados indican que a pesar de mostrar un ROC alto, el comportamiento es diferente en estas variables; así, la distancia a las viviendas presenta un coeficiente positivo, es decir, un aumento en una unidad de estas distancias aumenta las áreas transformadas, lo que revela que los cambios en el paisaje no se dan necesariamente en las cercanías a las viviendas. Por el contrario, en la distancia a las vías, el coeficiente es negativo, indicando que hay más transformaciones en áreas cercanas a estas. Este comportamiento es de esperarse, pues tal como se observa en la Figura 4(b), la distribución de vías, como son caminos y carreteras, atraviesan casi toda el área del páramo, facilitando de esta manera la accesibilidad y también el transporte de madera, leña y productos agrícolas. Resultados similares han sido encontrados en otros estudios como son Abadía (2011), Henríquez-Dole (2012) y Díaz-Pacheco y Hewitt (2013), demostrando que las vías son causantes de las transformaciones de los paisajes. La Tabla 3 muestra los resultados de regresión logística que se obtuvieron entre las áreas transformadas y las variables biofísicas relieve, distancia a ríos y pendientes.
Tabla 3 Valores de regresión logística entre áreas transformadas y las variables: relieve, distancia a río y pendientes en el Páramo de Paja Blanca
Fuente: Elaboración propia.
La Tabla 3 indica ajustes significativos en los modelos, donde el relieve y las pendientes al presentar un coeficiente positivo están mostrando que las transformaciones se acentúan en zonas altas y de bastante pendiente; esto se debe a que las condiciones de suelo son más apropiadas en estos lugares, lo que favorece para que se den procesos de transformación, ya sea para la agricultura y los pastos o también para el aprovechamiento de madera o leña.
En lo que respecta a la distancia a los ríos y quebradas, el modelo muestra buenos ajustes con coeficiente negativo, es decir, que las transformaciones se presentan en áreas cercanas a este componente biofísico. Este se explica en el sentido que el agua de las quebradas favorece para que se desarrollen actividades agropecuarias, como son los pastos y los cultivos. Resultados similares se han encontrado en otros ámbitos, como, por ejemplo, el estudio realizado por Gómez y Botero (2014), donde encontraron que las transformaciones del paisaje en las regiones del Bajo Cauca y Nordeste del departamento de Antioquia en Colombia fueron más acentuadas en zonas cercanas a los ríos.
4.3 Simulación del paisaje del páramo de Paja Blanca para los próximos 30 años
Una vez identificados y analizados los principales cambios de cobertura entre los periodos de 1984 y 2013, se procedió en primera instancia a la selección de variables, cuyo patrón espacial manifieste algún tipo de comportamiento que se crea que puede tener relación con el patrón de localización de cada uno de los cambios de coberturas que se identificaron. Las variables seleccionadas se indican en la tabla 4 con sus correspondientes valores del test de Cramer's V.
La Tabla 4 muestra que las pendientes, el relieve, la distancia a los ríos y las distancias a las viviendas aledañas ejercen cierto grado de influencia en la ubicación espacial de los mosaicos de pastos y cultivos. En las coberturas de páramo, las variables que ejercen mayor influencia fueron la pendiente y el relieve. El bosque natural denso no presentó relación con ninguna variable; y el bosque natural fragmentado manifestó tener influencia del relieve y la distancia a las viviendas.
Tabla 4 Índice de Cramer’s V para las variables biofísicas y las diferentes categorías de coberturas del suelo en el Páramo de Paja Blanca
Fuente: Elaboración propia.
4.3.1. Transiciones potenciales
Las transiciones más importantes de cambios de cobertura entre 1984 y 2013 fueron las siguientes:
Bosque natural fragmentado (BNF) a mosaico de pastos y cultivos (MPC)
Bosque natural fragmentado (BNF) a páramo (PA)
Mosaico de pastos y cultivos (MPC) a bosque natural fragmentado (BNF)
Páramo (PA) a bosque natural fragmentado (BNF).
Con estas transiciones, y considerando las variables anteriores, en el módulo LCM de Idrisi se realizaron los mapas potenciales de transición con interacciones que variaron entre 5.000 y 10 000, y niveles de precisión mayores al 80 %; los resultados se pueden apreciar en la Figura 6.
Fuente: Elaboración propia.
Figura 6 Mapas potenciales de transición de BNFMPC (a), BNF-PA (b), MPC-BNF (c) y PA-BNF (d), en el páramo de Paja Blanca, obtenidos en el módulo LCM de Idrisi
En estos mapas podemos observar que los valores del potencial en la escala de 0 a 1 son más amplios en la transición de BNF a MPC, y se acentúa hacia el sector occidental, en las veredas San José de Chillanquer y Villanueva, del municipio de Guachucal. Los valores de menos potencial se presentan para la transición PA-BNF.
Con los anteriores mapas potenciales de transición, en el programa Idrisi Selva, con el módulo LCM se simularon las coberturas para los próximos 30 años (1942); como resultado se ha obtenido una matriz de probabilidades de cambio (Tabla 5), una matriz de áreas de cambio (Tabla 6) y un mapa con la localización de las coberturas simuladas para el año 2042 (Figura 7).
Tabla 5 Matriz de probabilidades de cambio de coberturas del suelo en Páramo de Paja Blanca para el periodo 2013-2042
Fuente: Elaboración propia.
Tabla 6 Matriz de áreas de cambio de coberturas del suelo en el Páramo de Paja Blanca para el periodo 2013-2042
Fuente: Elaboración propia.
Fuente: Elaboración propia.
Figura 7 Mapa de coberturas del suelo del páramo de Paja Blanca, simuladas para el año 2042.
Observamos en la Tabla 5 que, para el año simulado, las probabilidades de permanencia de cada categoría de cobertura son casi más del 70 %; también se evidencia que existe una probabilidad del 28 % y 30 % de que las coberturas de BNF y BND respectivamente cambien a MPC, es decir, de continuar con las tendencias actuales, los bosques que aún predominan pueden llegar a alterarse por factores antrópicos. Sin embargo, se presenta una probabilidad del 30 % de que algunas áreas de MPC cambien a BNF, lo cual demuestra que la tendencia a futuro es a que, mediante procesos de conservación, restauración ecológica y regeneración natural, se reconviertan áreas de uso antrópico a bosques naturales, lo cual es importante para la estabilidad ecológica del páramo.
La Tabla 6 indica las áreas (hectáreas) que posiblemente pueden cambiar de una categoría de cobertura a otra en el periodo simulado; en ella se evidencia que sumando las áreas de BND, BNF y PA que cambian a otras categorías mediante procesos de deforestación, obtenemos 486,8 ha afectadas antrópicamente, que corresponde a una tasa de deforestación anual de 16 has/año; es decir, la tendencia a futuro con el comportamiento actual es de mantenerse la misma tasa de deforestación anual que se presenta para el periodo 1984-2013. De esta manera, la simulación advierte que el páramo de Paja Blanca se mantendrá como un páramo clasificado según Cabrera y Ramírez (2014) como de baja transformación.
Se muestra también en la Tabla 6 casi 400 hectáreas de BNF transformadas a MPC, indicando que ahora que Paja Blanca ha sido declarado como Parque Natural Regional, sería de mucha importancia así como lo plantea Morales-Betancourt y Estévez (2006), que las instituciones encargadas el manejo del páramo adelanten programas y proyectos en los cuales sus habitantes no sean catalogados como "problema" o como "generadores de conflictos", sino que se los involucre en procesos de recuperación del ecosistema y en la educación ambiental de los visitantes. Estos autores manifiestan que "[...] debe reconocerse que los habitantes del páramo son parte integral de éste y deben ser entendidos como parte de la comunidad biótica paramuna, no como un accesorio que debe ser removido o del que puede prescindirse" (p. 44).
La Figura 7 presenta la ubicación espacial de las coberturas del suelo que se simularon para los próximos 30 años. Según este mapa, las áreas de BNF que se afectan por procesos antrópicos está más acentuada en los municipios de Guachucal y Sapuyes, especialmente en las veredas San José de Chllanquer, Villanueva y Maramba, localizadas al extremo noroccidental del páramo de Paja Blanca. También indica que en el municipio de Sapuyes, en las veredas La Campana y Los Monos, se presentan importantes áreas de MPC que se convierten a BNF; situación similar ocurre en la vereda Villa Sur del municipio de Ospina, y en la vereda Chorrera Negra del municipio del Cotandero.
Cabe resaltar que en el año 2015, dos años después de que se formuló el presente proyecto, el límite oficial de Paja Blanca, en su declaración oficial como Parque Natural Regional, fue diferente al tomado en este estudio, ya que el actual límite considera una superficie de 3.107 ha, es decir, se excluyeron todas las áreas de mosaico de pastos y cultivos, manteniéndose las áreas de bosque natural denso, bosque natural fragmentado y páramo. Esto quiere decir que los datos hallados en este estudio no dejan de evidenciar una situación preocupante al estimar 400 ha de bosque natural fragmentado que se encuentran en el actual límite que pasaran a ser mosaico de pastos y cultivos. De igual forma, muestra la importancia de que los límites oficiales sean claramente delimitados en el terreno y, sobre todo, de que se adelante un programa de compra de predios para evitar su alteración limítrofe.
Desde el año 1997 hasta la fecha se han realizado tres (3) propuestas de delimitación, cada una de ellas condicionada por algunos factores fisicobióticos, especialmente las coberturas del suelo. Si no se establecen políticas claras y efectivas de conservación y manejo, es posible que dicho límite cada vez sea menor, como producto de los cambios en las coberturas que las simulaciones indican.
4.3.2. Validación del modelo
En vista de que no fue posible adquirir una imagen satelital actual (2016) para realizar un mapa y corroborar la exactitud del modelo, se decidió hacer un levantamiento de 50 puntos GPS en campo con el fin de registrar las coberturas observadas y con ello compararlas con un mapa simulado para el año 2016, a través de una matriz de confusión. La Tabla 7 indica los valores de la matriz de confusión, donde se puede apreciar una exactitud global del 76 %, y un índice de Kappa del 64 %, que según López de Ullibarri-Galparsoro y Pita-Fernández (1999) es considerado como bueno.
5. Conclusiones
En el periodo de 1984-2013, las transformaciones más sobresalientes corresponden al cambio de bosque natural fragmentado a mosaico de pastos y cultivos, y viceversa; indicando que hubo procesos de deforestación, pero también se presentó la reconversión de praderas y áreas agrícolas a bosques, que desde el punto de vista ambiental es bastante importante para la estabilidad del páramo.
La presencia de coberturas de páramo propiamente dicho, como son herbáceas y arbustivas, aproximadamente ocupa el 20 % de todo el territorio de Paja Blanca, cuyas transformaciones o cambio hacia otras categorías no fueron significativas o alarmantes, pues solo se perdieron 78 hectáreas. Esto quizá sea un indicador de efectividad de las gestiones adelantadas por Corponariño en pro de la protección y conservación del páramo.
Con los resultados obtenidos en este estudio, se presenta una discrepancia con respecto a la tasa de deforestación anual que reporta Corponariño (2007b), de 11 has/año, ya que el análisis multitemporal realizado indica una pérdida de cobertura de bosque de 16 hectáreas al año, donde el factor antrópico es el causante de dicho comportamiento.
La densidad poblacional no mostró ser un factor causante de las transformaciones en el páramo de Paja Blanca, ni siquiera la distancia a las viviendas; pero se encontró una estrecha relación entre la distancia a las vías y áreas transformadas. Esto demuestra que las actividades antrópicas que causan el cambio no necesariamente provienen de quienes habitan directamente en el páramo. La tendencia futura del páramo de Paja Blanca es a mantenerse la misma tasa de deforestación (16 has/año), situación que lleva a pensar que se necesitan medidas de gestión más contundentes que logren bajar esta tasa.
El modelo de simulación obtenido presentó un buen ajuste y precisión, convirtiéndose en una herramienta muy importante para los procesos de toma de decisiones de quienes se encargan de la administración y gestión del páramo; sin embargo, se necesitan explorar otras variables socioeconómicas y políticas que puedan tener influencia en los cambios de cobertura, y que lleven a obtener un modelo con mucha mayor precisión.
