Resumen:
Contextualización: las plantas invasoras, como el retamo espinoso (Ulex europaeus), impactan negativamente la biodiversidad y el desarrollo normal de los ecosistemas en los territorios donde está presente. Así mismo, estas plantas disminuyen el área disponible para siembra, lo cual conlleva a la disminución de las cosechas.
Vacío de conocimiento: se hace pertinente desarrollar metodologías de bajo costo y efectivas para identificar las zonas afectadas por este tipo de malezas, con el fin de que exista información para disminuir el impacto de estas al medio ambiente y a las comunidades campesinas de la zona.
Propósito: implementar una propuesta para la identificación de las zonas afectadas por retamo espinoso en el municipio de Sibaté Cundinamarca, debido a su presencia en grandes extensiones.
Metodología: primero, se realizó un trabajo con la comunidad campesina del municipio, identificando y georreferenciando algunas zonas donde la planta hace presencia; seguidamente, usando una imagen satelital y multiespectral [PlanetScope] de los lugares identificados y georreferenciados, se calcularon los valores de intensidad espectral de esta. A continuación, se analizaron estadísticamente los percentiles de la imagen y se cuantificaron los percentiles 25 y 75; después, se calcularon los píxeles de la imagen con algún grado de probabilidad de presencia de la planta invasora, utilizando el algoritmo de “Hot Spot” o “zonas calientes. Para finalizar, se validaron las zonas con presencia probable de la planta, utilizando vehículos aéreos no tripulados [UAV].
Resultados y conclusiones: se identificaron puntos y zonas donde posiblemente hay presencia de retamo espinoso en Sibaté [Cundinamarca]. Se hallaron dos áreas en las que se aplicó la metodología propuesta en el centro del municipio, con un 99 % de probabilidad de presencia de la planta invasora; a la vez, se evidenciaron ocho puntos con una probabilidad de 95 % de presencia de retamo espinoso y, por último, se encontraron ocho puntos con probabilidad del 90 % de la presencia de la planta invasora. Finalmente, se corroboraron parte de las áreas y puntos con posible presencia de la planta, implementando vuelos no tripulados. Se observó una alta correlación espacial entre las áreas georreferenciadas y las áreas corroboradas en campo, lo cual es relevante para procesos de restauración ecológica.
Palabras clave: retamo espinoso,teledetección espacial,imágenes satelitales,espectro radiometría satelital.
Abstract:
Contextualization: Invasive plants, such as thorny broom (Ulex europaeus), have a negative impact on biodiversity and on the normal development of the ecosystems in territories where it is present. Also, those plants reduce the area available for cultivation, which leads to decrease of the harvests.
Knowledge gap: it is pertinent to develop low-cost and effective methodologies to identify the areas affected by this type of weeds; that way, the obtained information could be used to reduce the impact of those plants on the environment and in the peasant communities of the area.
Purpose: To implement a proposal for the identification of the areas affected by the invasive thorny broom plant in the municipality of Sibaté, Cundinamarca.
Methodology: First, with the peasant community of the municipality, were identified and georeferenced some areas where the plant is present; next, using a satellite and multispectral image [from PlanetScope] of the identified and georeferenced areas, the spectral intensity values of the image were calculated. Then, statistical analyzes of the percentiles were realized, and the 25th and 75t. percentiles were quantified; after that, the pixels with some degree of probability of the presence of the invasive plant were calculated, using the "Hot Spot" algorithm or “hot zones”. Finally, the areas with probable presence of thorny brooms were validated, using unmanned aerial vehicles (UAV).
Results and conclusions: It was possible to identify points and areas where there is, possibly, presence of thorny brooms in the municipality of Sibaté, Cundinamarca. There are two areas where the proposed methodology was applied in the center of the municipality, with a 99% probability of having the invasive plant; at the same time, eight points with a probability of 95% of presence of the plant were evidenced; finally, were found eight points with a probability of 90% of the presence of the invasive plant. At the end, part of the areas and points with possible presence of the plant were corroborated, implementing unmanned flights. A high spatial correlation was observed between the georeferenced areas and the areas corroborated in the field, which is relevant for ecological restoration processes.
Keywords: Thorny broom (Ulex europaeus), remote sensing, satellite images, satellite spectroradiometry.
ÁREA AMBIENTAL
IDENTIFICACIÓN DE ZONAS CON PRESENCIA DE RETAMO ESPINOSO (Ulex europaeus), UTILIZANDO TELEDETECCIÓN ESPACIAL, EN EL MUNICIPIO DE SIBATÉ [CUNDINAMARCA]
IDENTIFICATION OF AREAS WITH PRESENCE OF THORNY BROOM (Ulex europaeus), USING SPATIAL REMOTE SENSING, IN THE MUNICIPALITY OF SIBATÉ [CUNDINAMARCA]
Edier Fernando Ávila Vélez edierf@hotmail.com
Recepción: 02 Junio 2022
Aprobación: 16 Septiembre 2022
Publicación: 17 Diciembre 2022
En la década de los años 50, debido a los procesos de colonización, en Colombia y en especial en su capital, Bogotá, se vivió un proceso avanzado de deforestación. Para contrarrestar este fenómeno, la administración local decidió realizar la plantación de la especie retamo espinoso (Ulex europaeus) (Empresa de Acueducto y Alcantarillado de Bogotá, 2021). Esta planta invasora es una especie que compite agresivamente con plantas nativas (Muñoz, 2009), causando pérdida de biodiversidad y el desplazamiento de algunas de ellas. A su vez, genera cambios en la estructura vegetal del ecosistema (Osorno, 2014).
Se han tenido reportes de la existencia y colonización de la mencionada planta en zonas de bosque alto andino y subpáramo (Sarmiento et al., 2016) de Colombia; igualmente, se ha documentado la presencia de tal planta invasora en otras partes del mundo, como Estados Unidos, Costa Rica, Sudamérica, India y Japón, entre otras. La presencia de retamo espinoso afecta drásticamente el sector agropecuario, disminuyendo las áreas de cultivo y desplazando la flora y fauna endémica (Sanguino, 2018).
El retamo espinoso crece en forma de matorral denso, conformando parches amplios de difícil penetración, y puede alcanzar hasta 7 metros de altura; esta planta acumula gran cantidad de ramas secas y en zonas del trópico florece todo el año (León y Vargas, 2011). Sumado a esto, la planta posee un alto poder calorífico, dada la cantidad de aceites en sus semillas, por lo cual, al quemarse, se propaga con mayor velocidad (Ocampo-Zuleta, 2019); infortunadamente, los habitantes y agricultores afectados por esta maleza [como medida de control] utilizan el corte, arranque y quema de la misma, sin considerar los daños ambientales que realizan al ecosistema y el beneficio que traen a su dispersión (Instituto Nacional de Investigación Agropecuaria [INIA], 2011). Sin embargo, se han implementado otras iniciativas en los territorios afectados por esta planta invasora, como su erradicación de los cerros orientales de Bogotá y la restauración ecológica (Barrera-Cataño, 2019; Sánchez et al., 2018; León y Vargas, 2011; Ocampo-Zuleta, 2019), mitigando en alguna medida su proliferación y afectación.
Algunos estudios acerca de la caracterización de zonas afectadas por la presencia de retamo espinoso se han realizado a través de visitas de campo y recorridos por predios y vías de los territorios (Sanguino, 2018), estas técnicas demandan mayor tiempo, costos económicos y pueden quedar áreas sin identificar. Debido a lo anterior, es necesario utilizar técnicas de mayor cobertura, como la teledetección espacial, (Beltrán, 2014; Everitt et al., 1987) para identificar las zonas con presencia de retamo espinoso y, así, contribuir a su control y posible erradicación. Por lo anterior, utilizar técnicas de teledetección espacial puede contribuir a mapear y delimitar áreas donde probablemente existe colonización de la planta invasora.
Estudios anteriores dan cuenta de la aplicación de técnicas de teledetección espacial en Colombia (Barreto, 2021) y en Uruguay (Sampayo y Vinçón, 2018), utilizando imágenes satelitales [Sentinel-2] de una resolución espacial media. También se encuentra lo realizado por Castro y Peñaloza (2015), que emplearon vehículos aéreos no tripulados en los cerros orientales de Bogotá [Colombia], pero con la limitante de menor cobertura espacial. Además, existen desarrollos informáticos con información relevante acerca de la identificación y cuantificación de zonas con presencia de retamo espinoso[1] (Secretaria Distrital de Ambiente de Bogotá, 2020), pero únicamente en la ciudad de Bogotá y en la zona rural de esta entidad territorial. No existe información de los municipios vecinos. Por eso, el objetivo de la presente investigación fue la identificación de zonas con presencia probable de dicha planta en el municipio de Sibaté [Cundinamarca, Colombia]. La presente investigación aporta una técnica innovadora, ya que se fundamenta en la utilización de sensores espaciales o teledetección de mayor resolución espacial y en el potencial de la resolución radiométrica [específicamente la longitud de onda del infrarrojo cercano] para la extracción de información temática de coberturas que pueden perjudicar y alterar el medio ambiente.
2.1. Zona de Estudio
La presente investigación se desarrolló en el municipio de Sibaté [Cundinamarca, Colombia], con latitud norte de 4°26’14’’ y 74°14’05’’ de longitud oeste, aproximadamente; posee una altura sobre el nivel del mar de 2700 metros y es considerada una despensa agrícola importante en la región. Esta zona está siendo afectada por retamo espinoso en áreas de conservación ecosistémica, como el Páramo de Sumapaz, y otras empleadas para la producción agrícola familiar y/o pastoreo.
2.2. Materiales
Como insumos satelitales, se adquirieron imágenes multiespectrales PlanetScope[2] (L3Harris Geospatial Solutions, 2022). de la zona de estudio, las cuales presentan los siguientes parámetros:
2.3. Metodología
Para la investigación se realizó un diseño experimental basado en cuatro variables independientes y una dependiente. Las cuatro variables independientes son las bandas o longitudes de onda que tiene la imagen satelital PlanetScope1:
- VI1 Banda Blue 455–515 nm
- VI2 Banda Green 500–590 nm
- VI3 Banda Red 590–670 nm
- VI4 Banda NIR 780–860 nm
VI1: Variable Independiente 1 y así sucesivamente
Y la variable dependiente es la presencia, o no, de retamo espinoso en la coordenada identificada en campo sobre el píxel o celda de la imagen satelital.
Hipótesis
Lo que se buscó fue una medida o un análisis estadístico-descriptivo que explicara el comportamiento de los valores radiométricos presentes en las longitudes de onda o bandas espectrales en las que se observa la presencia de retamo espinoso, teniendo en cuenta esta condición: si se modifican los valores o rangos radiométricos de las bandas o longitudes de onda, se esperaría que el píxel no tenga presencia de la planta invasora. Con estas premisas se formuló un análisis estadístico-descriptivo de los datos que pudiese trabajarse por rangos, debido a que la interacción de la energía electromagnética, capturada por las imágenes satelitales, está influenciadas por variables que no se pueden controlar, como la radiación solar, el estado atmosférico, el clima, la contaminación entre otras, y pueden variar de un lugar a otro en varios metros. En consecuencia, y según la experiencia de los investigadores, se decide trabajar por rangos radiométricos y no por medidas puntuales.
Análisis estadísticos
Con las características mencionadas anteriormente, se decidió implementar datos estadísticos-descriptivos de percentiles, los cuales son medidas de posición y se calculan con los datos de la tabla 2.
El primer paso de la investigación, después de tener los datos de los percentiles, fue conseguir las imágenes de los terrenos afectados por retamo espinoso. Se observa, en la figura 2, la distribución espacial de las tomas realizadas en campo, en los predios que estaban afectados por la planta invasora.
El segundo paso fue dividir las imágenes tomadas en campo, según criterios como el área afectada, la accesibilidad al terreno y la señal de GPS para su georreferenciación, en dos grupos. El primer grupo de imágenes se utilizó para analizar estadísticamente los valores de reflectancia [21 muestras]. Y el segundo grupo, compuesto por 6 muestras, se usó para verificar los datos arrojados por el procesamiento. A continuación, se obtuvieron los valores de los percentiles 25 y 75 con el fin de caracterizar los tamaños por medio de este dato estadístico (Cano et al., 2016) y asegurar, así, que los valores de reflectividad de retamo espinoso estuvieran presentes [tabla 3 y figura 3].
Al tomar el rango espectral entre los percentiles 25 y 75, estamos garantizado que el 75 % de los datos espectrales están por debajo de 70 niveles digitales y que el 75 % de los valores están por encima de 38 niveles digitales para la VI1 banda Blue [455 – 515 nm]. Por tal motivo, con estos valores tendríamos un rango de 38-70 niveles digitales para la banda VI1, y así con las demás variables independientes: VI2, VI3, VI4. Obsérvese la figura 3 y figura 4 (b) de este manuscrito.
El tercer paso fue aplicar los valores estadísticos de percentiles a la imagen satelital adquirida por cada banda o longitud de onda, según la tabla 3.
El cuarto paso consistió en aplicar un análisis estadístico espacial, denominado “Optimized Hot Spot Analysis”, al entorno del software ArcMap (Environmental Systems Research Institute [Esri], 2020). Esta técnica está muy bien documentada y validada, como una de las mejores para encontrar patrones espaciales, basados en la autocorrelación espacial de las entidades geográficas (Goovaerts, 2001; Liu et al., 2022).
Análisis Geoestadísticos
Del anterior análisis estadístico se obtuvieron 472 píxeles o celdas que cumplen la condición de la Tabla 3. El anterior proceso se realizó en el software ArcGIS 10.6.
Con la anterior capa de información temática se procedió a aplicar la técnica para detectar zonas calientes, denominada en el ámbito de los sistemas de información geográfica como “Optimized Hot Spot Analysis”. La anterior técnica se realizó con el software informático ArcGIS 10.6., usando la extensión de herramientas de estadística espacial. La técnica de zonas calientes analiza la dispersión o la generación de conglomerados estadísticamente significantes, por lo cual, cuando se aplica la herramienta, se genera un valor de “GiZScore”, el cual indica un conglomerado fuerte a nivel de entidades geográficas y con niveles de confianza significativos del 90 %, 95 % o 99 %. Los valores negativos se entienden como baja dependencia o autocorrelación para crear un conglomerado.
Como resultado principal, después de aplicar la metodología propuesta, se lograron identificar puntos probables o conglomerados de píxeles que son afectados por la presencia de retamo espinoso figura 5.
En la figura anterior se puede observar una zona con alta probabilidad de presencia de la planta invasora y ocho zonas con una probabilidad de presencia cercana al 95%. Se validó la metodología con los píxeles o muestras de verificación tomados en campo, encontrándose una gran correlación espacial entre ellas.
Para poder ubicar exactamente la distribución de retamo espinoso es necesario conocer ciertas características de esta planta, como su evolución fenológica y territorial. Teniendo en cuenta esos datos, se observaron algunos patrones espaciales de aparición de forma lineal, y a veces poligonal, muy cerca de las vías y cauces intermitentes de agua. Dada las características del sensor utilizado, no fue posible concluir explícitamente las zonas de aparición de la planta; por tal motivo se realizaron capturas de información, con una resolución espacial más detallada, utilizando vehículos aéreos no tripulados, los cuales pueden llegar a obtener resoluciones espaciales del orden de centímetros. Otra forma de ganar resolución espacial, sin perder resolución espectral, es utilizar la fusión de imágenes o bandas, como lo realizado por Wang et al. (2016), ya que con esta técnica se generan productos con gran capacidad en el detalle espacial.
Debido a lo anterior, se realizaron salidas para validar los datos calculados con la implementación de la metodología utilizada en este estudio, logrando corroborar varios puntos a través de vuelos realizados con aeronaves no tripuladas (UAV), como se evidencia en la figura 7. Por tal motivo, se puede inferir gran capacidad de adaptación del retamo espinoso a cualquier medio (Sanguino, 2018; Sarmiento et al., 2016), lo que genera una amenaza para el ecosistema, ya que las semillas pueden ser arrastradas por las corrientes superficiales de los cuerpo y cauces de agua y vientos, debido al tránsito de automotores.
Con la metodología utilizada en esta investigación se busca mostrar una manera práctica de clasificar píxeles, utilizando una técnica muy sencilla y menos costosa desde el punto de vista económico; debido a que, en muchas áreas del conocimiento, se desconocen las potencialidades de las imágenes satelitales y en general de las técnicas de la geomática (Pacino, 2018) para la extracción de información temática.
En este estudio se logró, a través de la teledetección espacial y radiometría satelital, caracterizar espectralmente el comportamiento de retamo espinoso en un territorio específico, como lo realizado por Ávila-Vélez et al. (2019) y Ormazábal et al. (2013), que caracterizaron cultivos de maíz utilizando imágenes satelitales y observaron la evolución espacial de los bosques utilizando imágenes Landsat-5, aplicando espectrorradiometría de campo [ esta técnica mide la energía reflejada de los elementos del paisaje]. Existen otras técnicas para clasificar elementos de un territorio usando imágenes satelitales (Haboudane et al., 2004; Perea et al., 2009; Stefanov et al., 2001; Yu et al., 2006); sin embargo, dependen en alto grado de la composición física y química del elemento en el territorio.
Se puede inferir, de acuerdo con lo hallado en este estudio, que el retamo espinoso posee una distribución espacial lineal, adaptada a los elementos del territorio y utiliza como agente multiplicador estos medios: las vías de comunicación, los vientos, las corrientes de aire y las escorrentías de agua que surgen a los lados de las vías. Este concepto es emitido a partir de los vuelos realizados con aeronaves no tripuladas, los cuales se efectuaron para corroborar los datos obtenidos con las técnicas geomáticas de este estudio.
CÓMO CITAR: Ávila-Vélez, E., Herrera-Martinez, L. y Celis-Forero, A. (2023). Identificación de
zonas con presencia de retamo espinoso (Ulex europaeus), utilizando teledetección espacial, en el
municipio de Sibaté [Cundinamarca]. Revista
de Investigación Agraria y Ambiental, 14(1), 201 - 219. https://doi.org/10.22490/21456453.5840
CONTRIBUCIÓN
DE LA AUTORÍA: Edier
Fernando Ávila-Vélez: Introducción, metodología, análisis
de datos, salidas de campo, resultados y escritura del borrador. Luis Alberto Herrera-Martínez:
metodología, análisis de datos, conclusiones, revisión y edición. Álvaro Celis-Forero: logística,
metodología, salidas de campo, revisión y edición.
CONFLICTO
DE INTERESES: Los autores declaran no tener ningún conflicto de intereses.
https://hemeroteca.unad.edu.co/index.php/riaa/article/view/5840 (html)
https://hemeroteca.unad.edu.co/index.php/riaa/article/view/5840/5961 (pdf)
Los autores agradecen a los habitantes y agricultores de la vereda “El Romeral” del municipio de Sibaté [Cundinamarca], a la Universidad de Cundinamarca, en alianza con el Ministerio de Ciencia, Tecnología e Innovación, por su apoyo para el desarrollo del proyecto 62851 “Diseño de estrategias para el control de retamo espinoso (Ulex Europeaus L.), en los municipios de Sibaté, Soacha y Pasca” y al grupo de investigación Geocartografía, perteneciente a la Universidad de Cundinamarca, por su apoyo en las salidas de campo.
