Determinación del riesgo a la erosión potencial hídrica en la zona cafetera del Quindío, Colombia
Determination of the risk to the potential erosion by water in the coffee zone of the Quindio, Colombia
A avaliação de risco para o potencial de erosão hídrica na região do café de Quindio, Colombia
Andrés Felipe Castro Quintero1, Luz Adriana Lince Salazar2 & Orlando Riaño Melo3
1Ingeniero Agrónomo. 2Geóloga, Magister en Desarrollo Sostenible y Medio Ambiente. 3Matemático, Ingeniero de Sistemas y Computación, Especialista en Sistemas de Información Geográfica, Magister en Geomática
1, 2Centro Nacional de Investigaciones de Café- CENICAFE. Manizales. Colombia.
3Maestría en Teledetección. Universidad Católica de Manizales. Colombia
1andres.castro@cafedecolombia.com, 2luz.lince@cafedecolombia.com, 3orianom@unal.edu.co
Resumen
En la presente investigación se estimó la erosión potencial hídrica –EP-, la cual causa degradación y afectación a la fertilidad del suelo reduciendo el rendimiento del cultivo e incrementando los costos en producción. Se utilizó la Ecuación Universal de Pérdida de Suelo (USLE/RUSLE), con la formula A= (R)*(K)*(LS), se realizó muestreo aleatorio simple, en once Unidades de suelo como repeticiones. Se trabajó con el factor de erosividad de la lluvia R reportados por Lince (2015), con precipitación mensual de 19 estaciones meteorológicas para el periodo 1999-2011. El cálculo del Factor (K) se cuantificó utilizando un simulador de lluvia, los datos obtenidos fueron interpolados con kriging, y por último se halló el factor topográfico –LS-, mediante el tratamiento a la imagen Aster GDEM, propiedad de la Nasa. Luego se procedió a utilizar el álgebra de mapas. Se encontró que la EP en los municipios de Buenavista, Calarcá, Córdoba, Pijao y Génova (Quindío, Colombia) presentaron más de 27 % del área de estudio con valores de calificación de pérdida de suelo Moderada, Alta y muy Alta que van desde las 26 a más de 100 t/ha-1.año-1, y los municipios de Armenia, Circasia, Filandia, Montenegro, Quimbaya y La Tebaida, presentan EP con calificaciones consideradas como Ligera a Baja, representando hasta el 73% del área de estudio con valores que van desde 10 a 25 t/ha-1.año-1, y correspondiendo a pendientes poco pronunciadas, lo que representa pérdidas potenciales considerables en la capa orgánica del suelo.
Palabras clave: erosión potencial, interpolación Kriging, ponderación por distancia y simulador de lluvias.
Abstract
This research estimated erosion potential - EP-, which causes degradation and affecting the fertility of the soil by reducing the crop yield and increasing costs in production. Used the Universal equation of soil loss (USLE/RUSLE), with the formula A = (R) * (K) * (LS), simple random sampling, in eleven ground units such as replications was conducted. He worked with R rain rainfall erosivity factor reported by Lynx (2015), with monthly precipitation of 19 weather stations for the period 1999-2011. The calculation of Factor (K) was quantified using a rainfall Simulator, data were interpolated with kriging, and finally found the topographic factor - LS-, through treatment in the Aster GDEM image, owned by Nasa. They then proceeded to use the algebra of maps. Found that the EP in the municipalities of Buenavista, Pijao, Calarca, Córdoba, Genoa (Quindío, Colombia) presented more than 27% of the study area with moderate soil loss rating values, high and very high ranging from 26 to more than 100 t/ha-1.ano-1, and the municipalities of Armenia, Circassia, Filandia, Montenegro, Quimbaya, and La Tebaida they present EP with qualifications considered light to low, representing up to 73% area of study with values that range from 10 to 25 t/ha-1.ano-1, and corresponding to little steep slopes which represents significant potential losses in the organic layer of the soil.
Key-words: potential erosion, interpolation Kriging, weighting by distance, rain simulators.
Introducción
La erosión hídrica de suelos en Colombia es un problema que afecta directamente la fertilidad del suelo, ocasionando baja producción de los cultivos e incrementando los costos en producción Rivera & Gómez (1991). En la zona cafetera hay procesos naturales que favorecen erosión como lo son las fuertes lluvias, material parental susceptible y las altas pendientes. En este sentido Castillo (2004), menciona que la erosión hídrica puede convertir el terreno en áreas degradadas. Mientras Sagredo (2005), afirma que los índices máximos se presentan cuando existe un uso y manejo inadecuado de estos a favor de la pendiente. Por lo tanto, se hace necesario estimar la pérdida del suelo debido a la erosión hídrica en el departamento del Quindío, y de esta manera identificar zonas donde ocurre este proceso. Para calcular la pérdida de suelo en este caso para erosión potencial se hace por medio de la implementación de un modelo matemático empírico para el cálculo de la pérdida de suelo USLE/ RUSLE, desarrollado por Wischmeier & Smith (1978) y revisado por Renard, Meyer & Foster (1989).
Se aplicó la metodología de Hincapié & Rivera (2003) obteniendo los datos numéricos para determinar la erosividad de la lluvia (R) en MJ.mm.ha-1, adquiriendo datos de precipitación de 19 estaciones climática, aplicando la fórmula matemática del índice de Fournier IFM, reportado por Arnoldus, Boodt & Gabriels (1980) y el EI30 propuesto por Rivera & Gómez (1991) y Rivera et al. (2010) para la zona centro cafetera. La determinación del índice de Erodabilidad (Ki) en t/ha/MJ-1.mm-1.ha-1, el cual es cuantificado utilizando el simulador de lluvia, determinando la susceptibilidad del suelo a erosionarse por acción de la gota de la lluvia, por último se tuvieron en cuenta los factores topográficos LS, utilizando la metodología citada por Cocuyame & Salazar (2015). Donde, para ello se utilizó un modelo digital de elevación (ASTER GDEM), calculando la longitud y ángulo de la pendiente. En este estudio se asume que no existe cobertura alguna (C=1) y que el suelo se encuentra sin prácticas de manejo de conservación (P=1), y están a favor de la pendiente (Montes-León, Uribe-Alcantara, García-Celis, 2011), la cual permite estimar cuantitativamente la erosión potencial, representada por la pérdida de suelo en t/ha/año.
Los datos obtenidos fueron proyectados, utilizando tecnologías de información Geográfica TIG, que se convierte en herramienta de gestión ambiental y control de riesgos para analizar y prevenir la pérdida de suelos productivos (Rodríguez et al., 2004). Además en la actualidad para la región cafetera de Colombia, se desconoce cuánta es la pérdida potencial de suelo por erosión, por tal motivo la información debe conocerse a escala regional o local, con el fin de diseñar las estrategias acordes a dichas condiciones, que en muchos casos la baja adopción de las prácticas de conservación, puede obedecer a recomendaciones generales que no corresponden a las condiciones ni a las necesidades de los productores. Por lo anterior, en esta investigación se plantea llevar a cabo una predicción de erosión potencial en la zona cafetera del departamento del Quindío. Siendo necesario demostrar el deterioro real para crear conciencia, y proponer soluciones desde el interior de las comunidades.
Materiales y métodos
Área de estudio.
El departamento del Quindío está ubicado en la región centro-occidental de Colombia; está dividido en 12 municipios: Armenia, Buenavista, Calarcá, Circasia, Córdoba, Filandia, Génova, La Tebaida, Montenegro, Pijao, Quimbaya y Salento, cuenta con una superficie de 1.845 km2. En la mayoría de la zona cafetera del departamento, desde Filandia en el norte, hasta Génova en el sur. Cuenta con una temperatura media entre 18 y 24°C, y un promedio anual de lluvias de 2.000 a 4.000 mm.
Procesamiento
En el procesamiento se seleccionaron suelos desarrollados bajo distintas condiciones climáticas a lo largo de gradientes, para cada sitio se identificó el material parental de los suelos, teniendo en cuenta los sitios por rangos altitudinales (1200 a 1400 msnm – 1400 a 1600 msnm y 1600 a 1800 msnm, sembrados con cultivos de café. Se elaboró un mapa de muestreo aleatorio simple identificando once unidades de suelos representativas de la zona de estudio (Figura 1), distribuida de la siguiente forma: unidad Montenegro (11 muestras), Quindío (6 muestras), Génova (5 muestras), Chinchiná (3 muestras), Buenavista (3 muestras), Malabar (2 muestras), Filandia (2 muestras), La Concha (1 muestras), Patiobonito (1 muestras), Aluvial (1 muestras) y Sedimentario (1 muestras), el número de repeticiones está determinado por el área de la unidad cartográfica.
Figura 1. Ubicación de la zona de estudio.
Fase de Campo. Se recolectaron 36 muestras de suelo a una profundidad de 0 a 10 cm, para un total aproximado de 30 kg por muestra, estas fueron empacadas en bolsas herméticas, llevadas a laboratorio y procesadas para realizar análisis físico y químico a los suelos.
Fase de Laboratorio. Las muestras fueron secadas al aire libre bajo sombra, se separó 1kg de suelo para análisis de propiedades físicas y 1 kg para análisis químico, el suelo sobrante fue pasado por un tamiz con diámetro , para ser utilizado en construcción de microparcelas del simulador de lluvias.
Recopilación de datos primarios. Como mapa base, se empleó una imagen ASTER (Global Digital Elevation Model –GDEM), en formato raster con un tamaño de píxel 30m por 30 m. Dicha información geoespacial se adquirió de la NASA, (https://asterweb.jpl.nasa.gov/gdem.asp) Para la precipitación se obtuvo información histórica contenida en la red climática de Federación Nacional de Cafeteros Cenicafe y datos de erodabilidad.
Características del simulador de lluvias. El equipo para determinar el factor de erodabilidad (k) es del tipo Purdue Programable, que cuenta con Boquilla Vee–Jet 80 – 100; el cual se calibro con una precipitación de 80 mm/hora, y una presión de 8 PSI. Las microparcelas se graduaron con inclinación del 70% que hace referencia a la pendiente promedio de la zona cafetera colombiana.
Descripción del procedimiento. Para calcular la pérdida de suelos por erosión hídrica se aplicó la ecuación universal de pérdida de suelo Wischmeier & Smith (1978), usando el proceso en el álgebra de mapas, con el software ArcGIS (ESRI, 2008; Figueroa-Jáuregui et al., 2011; Flores López et al., 2003), La cual consiste en crear un mapa raster para cada uno de los factores de la ecuación, la importancia radica en que puede ser usada en cualquier región geográfica modificando sus factores de acuerdo a las características propias de la región de estudio y así cuantificar la erosión potencial, usando una relación funcional de varios factores expresados en la ecuación USLE/RUSLE (ecuación 1):
Siendo: A la pérdida de suelo por unidad de superficie (t/ha.año), R es el factor lluvia o índice de erosión pluvial (MJ. mm /ha h), K es el factor erodabilidad del suelo (t.ha.h / ha.MJ.mm).determinado en simulador de lluvia, L es el factor longitud de pendiente (adimensional) y S es la pendiente (adimensional). No se tienen en cuenta los factores C y P, para el cálculo de la erosión potencial EP.
Determinación de los factores de la ecuación
El Factor de erosividad (R) lo definen varios investigadores como la agresividad climática, siendo la fuerza de las gotas de agua lluvia, las que causan la erosión, y una forma de estimar el factor R es emplear la expresión utilizada por Ramírez et al. (2007) y Lince & Castro (2015), los cuales trabajan con índice de Fournier modificado (IMF), el cual corresponde al promedio mensual propuesto por Arnoldus et al. (1980) y el índice de erosividad o capacidad de las lluvias de erosionar el suelo (EI30), recomendado por Rivera (1991) (ecuación 2). La precipitación de la zona de estudio se caracterizó mediante información histórica del periodo de 15 años, para ello se seleccionaron 19 estaciones meteorológicas. Para estimar el factor R se emplea la extensión Spatial Analyst, el método de interpolación distancia inversa ponderada IDW. Se obtuvo una superficie estimada a partir de este conjunto de datos los puntos dispersos en el área geográfica. El factor R utilizado fue de datos aportados en investigaciones anteriores para el departamento de Quindío por Lince & Castro (2015).
Luego de ello se obtienen los valores de EI30 para la zona de estudio (ecuación 3)
Factor (K) o Índice de Erodabilidad. Según (Figueroa et al., 1991) se define como la susceptibilidad de un suelo a la erosión, este factor se realizó en 11 unidades de suelo, mediante el simulador de lluvia, el cual, cuantificó la pérdida de suelo arrastrado por escorrentía, salpique y depósito en canal y agua escorrentía e infiltrada, el suelo recolectado en cada prueba se dejó decantar durante 48 horas, luego fue secado en estufa a 105°C durante 48 h, y pesado en una balanza de precisión de 0.0001 g; los datos obtenidos fueron modelizados con métodos geoestadísticos con interpolación (KRIGING).
Factor topográfico (LS). Es la distancia horizontal desde el punto más alto, donde se origina el flujo superficial al punto más bajo donde comienza la deposición que fluye por escorrentía a un canal, mediante la ecuación usada por USLE/RUSLE (Renard et al., 1997), implementando el modelo de dirección de flujo y acumulación de flujo con la extensión de arcHydro para Arcgis. (Foster et al., 1977, citado por Barrios & Quiñonez, 2000). Se implementa el mapa de pendientes y se reclasifica, luego se procede a obtener la distribución espacial de LS. En el cálculo del factor L, fue necesario aplicar un relleno al DEM para evitar espacios vacíos en el raster, esto se elaboró mediante la herramienta “Fill” del software empleado, con la imagen resultante se creó un raster de pendientes –Slope- el mismo que fue utilizado para crear el factor empleando las ecuaciones 4 y 5.
Obtenido el Factor β se procedió a obtener el factor M que también es un dato a considerar en la obtención del Factor L, para ello se emplea la ecuación 6 en “Raster Calculator”
Para obtener el cálculo del factor L con el área de drenaje aportadora se procedió a realizar la dirección de flujo y acumulación del flujo con las herramientas Flow Direction y Flow Accumulation respectivamente. Obtenidas estas dos imágenes se halló el Factor L en “Raster Calculator” empleando la Ecuación 6 y formula Desmet & Govers (1996) citado por (Shi et al., 2004).
Para el cálculo del Factor S se empleara la ecuación 7.
De igual modo se obtiene el ráster subfactor de (S) la cual es la pendiente del terreno, donde el ángulo β se toma como el ángulo medio a todos los subgrids en la dirección de mayor pendiente según (Mccool et al., 1987), citado por Barrios & Quiñonez, 2000). Menciona Velásquez (2008) que, cuando se aplica esta fórmula en el Raster Calculator de ArcGIS se debe tomar en cuenta que el ángulo que deberá ser convertido a radianes (1 grado sexagesimal=0,01745 radianes), para poder ser multiplicado por los demás componentes de las ecuaciones. Una vez obtenidos todos los factores anteriores procedió a calcular el factor LS que es objeto de la presente metodología. Para ello se utilizó la ecuación 8.
Erosión Hídrica Potencial
En el caso de la erosión hídrica potencial, se multiplican los factores de la ecuación, hallados en las ecuaciones anteriores.
La clasificación de pérdida de suelo potencial se comparó con la propuesta por Pérez (2001).
Programación de la ecuación.
Para fines prácticos se estructuro en la herramienta ModelBuilder, el modelo de geoprocesamiento automatizando mediante análisis espacial y de administración de datos, además se representa con un diagrama que concadena secuencias de procesos y herramientas de geoprocesamiento utilizando la salida y entrada de las variables de la ecuación (Figura 2) en esta metodología se adaptó para reemplazar símbolo por F, para facilidad el modelo (Quintana et al., 2009).
Figura. 2. Diagrama de la estructura con Modelbuilder para EP.
Resultados y discusión
Factor de erodabilidad (K). En la Figura 3 se ilustran los valores del factor K de erodabilidad en cada municipio, los cuales van desde 0.0002 hasta 0.0077 t/ha/MJ-1.mm-1.ha-1. Siendo este el rango de susceptibilidad natural como lo indican Rivera & Gómez (1991). En los municipios Buenavista, Calarcá, Córdoba, Génova, Montenegro, Pijao y Quimbaya; se presentan valores de susceptibilidad altos que van desde 0.0041 a 0.0077 t/ha/MJ-1.mm-1.ha-1; siendo los municipios de Buenavista y Pijao, los que presentan mayor porcentaje del área susceptible con valores de 90 y 83 % en área total de superficie. Esto corresponde a suelos con pendientes superiores al 30 % calificada como fuerte.
En los municipios de Génova, Montenegro y Quimbaya muestran valores de susceptibilidad Moderados que van desde 0.0041 a 0.0057 t/ha/MJ-1.mm-1.ha-1, con porcentajes susceptible con valores 5, 25 y 60% en su área de estudio. Mientras que los municipios de Armenia, Circasia, Filandia y La Tebaida corresponden a los valores susceptibilidad Baja con valores que oscilan entre 0.0002 a 0.0041 t/ha/MJ-1.mm-1.ha-1. Presentando el menor porcentaje del área susceptible comprometida con valores del 10, 15, 8 y 30 % de su área total de la superficie; esto se puede deber a que las pendientes de estas zonas están clasificadas como ligera a bajas 0 a 15 %.
Figura 3. Mapa de erodabilidad del departamento de Quindío
Pérdida de suelo potencial (EP). En la Figura 4, como se muestra en la parte baja del mapa, se observa la pérdida de suelo con valores que van de 9.431 a 202.17 t/ ha-1. año-1. Los cuales se clasificaron con las pérdidas de suelos propuesta por Pérez (2001). Es importante mencionar que se realizó la proyección con los promedios de pérdidas potenciales de suelos en la zona cafetera del departamento del Quindío que es el resultado de multiplicar, el factor R agresividad de la lluvia, Factor K de erodabilidad o susceptibilidad de los suelos a la erosión y por último longitud y gradiente de la pendiente, llevándolo a condiciones extremas donde no existirían cobertura ni prácticas de conservación.
Los municipios de Buenavista, Calarcá, Córdoba, Génova, y Pijao, presentaron el 27 % de área total de superficie con riesgo a la erosión potencial. Esto municipios se clasificaron de Moderada, Alta y Muy Alta con pérdidas que van desde 26 o más de 100 t/ ha-1. año-1. En los municipios de Armenia, Circasia, Filandia, Montenegro, Quimbaya y La Tebaida presentan calificación de Ligera a Baja en un porcentaje del 73 % de su área total de superficie con valores que van desde 10 a 25 t/ ha-1. año-1. En los municipios de Armenia, Montenegro y La Tebaida fueron los que menores pérdidas presentan. Aunque las pérdidas fueron menores puede llegar a causar daños irreparables en la capa orgánica del suelo, donde se encuentra la mayor cantidad de raíces del café, encargadas de la nutrición de la planta. El espesor del suelo puede perder de 2 a 3 mm por cada 30 t.ha.año como lo reportan estudios previos de CENICAFE en 1975, lo que significa que 27 % del área de estudio puede llegar perder entos 2 mm. ha . año y en un periodo de 10 años por hectárea puede llegar a perder 2 cm de espesor de la capa fértil del suelo, lo cual genera infertilidad. Por lo tanto, para cada uno de los municipios se clasificó la pérdida de suelo según la Tabla 1.
Tabla 1
Clasificación de la pérdida de suelo de acuerdo con la calificación
Figura 4. Mapa de Erosión Potencial. Pérdida de Suelo potencial en el departamento del Quindío
Conclusiones
El factor de erodabilidad más alta y que cubre mayor área de estudios se presenta en los municipios de Buenavista y Pijao, con valores de 0.0041 a 0.0077 t/ha/MJ-1.mm-1.ha-1, lo cual coincide con los municipios que se clasifican con riesgo potencial Moderada, Alta y Muy Alta y compromete más del 90 y 83 % de su área en estudio, lo cual es un indicador para la toma de decisiones a la hora de implementar prácticas de conservación.
La ecuación de pérdida de suelo USLE/RUSLE es un herramienta que permite generar los mapas de riesgos de la erosión potencial y así identificar las posibles zonas que pueden llegar a presentar baja fertilidad y desde el tema ambiental la que genera mayores aporte de sedimentos.
El material parental, la erosividad y las pendientes jugaron un papel fundamental en las pérdidas de suelo, como se pudo observar en los municipios de Montenegro, Circasia, La Tebaida y Armenia, donde las pérdidas fueron clasificadas de Ligera a Baja con valores menores 25 t.ha.año, que correspondían a pendientes del 15 y 30 %.
Agradecimientos
Los autores agradecen al Comité Departamental de Cafeteros de Quindío y a las disciplinas de suelos y agroclimatología de Cenicafé.
Literatura citada