1. Introducción
El monitoreo satelital inició a finales de la década de los cincuenta y ha permitido el entendimiento de los ecosistemas silvestres revelando aspectos importantes de su diversidad y complejidad (Ruiz, Savé, y Herrera Herrera, 2014). La recepción remota para monitorear en forma regular los cambios de la vegetación se realiza con algoritmos especializados, dando como respuesta los cambios ambientales (Leiva, 2003).
Por otro lado, para caracterizar la viabilidad de la agricultura de precisión se utilizan herramientas tecnológicas como son los Sistemas de Posicionamiento Global (GPS, por sus siglas en inglés), sensores planta-clima-suelo e imágenes multiespectrales obtenidas a partir de satélites, aviones o UAS/RPAS que permiten el manejo de la variabilidad espacial y temporal a nivel de subparcela de campo para mejorar el retorno económico y reducir el impacto ambiental (Díaz García-Cervigón, 2015).
La tuna es una planta de gran importancia en los sistemas agro-pastoriles de los valles interandinos que cuentan con las condiciones adecuadas para su asentamiento. Los frutos pueden ser consumidos en forma natural o en sus derivados como mermeladas y bebidas (López Orozco, Mercado Flores, Martínez Soto, y Magaña Ramírez, 2011).
Para el cultivo de la tuna los requerimientos agroclimáticos son los siguientes: los suelos deben tener una textura franca, franco arcilloso arenosa, arenosos, franco arenoso con un pH 6.5-8.5, aunque el cultivo de tuna es tolerante a la falta de agua sí se debería tener acceso a algún tipo de fuente de la misma con el fin de obtener mejores rendimientos. En cuanto a las condiciones climáticas apropiadas necesita de una temperatura anual de 16-28º C y una precipitación pluvial media anual de 150-1800 mm (Amaya Robles, 2009).
Por otro lado, en una investigación realizada por el Diario El Comercio el 15 de octubre de 2011, se señala que en el Ecuador poco o nada se conoce sobre el manejo de cultivo de tuna; a nivel nacional existen alrededor de 180 hectáreas de tuna sembradas y, puntualmente, indican que existen los siguientes tipos: amarilla con y sin espinas, la blanca y la silvestre repartidas en las provincias de Imbabura, Cotopaxi, Tungurahua, Chimborazo, Loja y Manabí (El Comercio, 2011).
La revista Líderes del Ecuador, señala que para producir tuna no se requiere de gran cantidad de agua ni de terreno húmedo, ya que existe una buena aceptación del producto en el mercado local por sus propiedades nutricionales y medicinales (Márquez, 2018).
El diario La Hora en octubre de 2011 publicó un artículo en el que señalaba que el cultivo de tuna ha ido ganando terreno en el norte del Ecuador, gracias a la creación del plan de producción e industrialización de la tuna en el año 2010, produciendo dos variedades la blanca y la amarilla ( La Hora, 2018)
El Diario El Telégrafo, el 24 de octubre de 2017, publicó que en el cantón Guano, al momento existen plantaciones de tuna ubicadas en las comunidades: Hela, Santa Rosa de Guano, Miraflores, Carrera Ambato, Chingazo Alto, Chingazo Bajo, Alacao, San José de Chocón, Pungal San Pedro y San José de Juntuz (Maggi, 2017).
Los Índices de Diferencias Normalizadas de Vegetación (NDVI, por sus siglas en inglés), índice de estrés hídrico (CWSI, por sus siglas en inglés) y el índice de vegetación ajustada al suelo (SAVI, por sus siglas en inglés) utilizan las propiedades de reflexión diferencial de la vegetación en el espectro visibles e infrarrojo para monitorear la vegetación tanto en tiempo como en espacio (Gilabert Navarro, García-Haro, y González Piqueras, 1997). Es útil para estimar el índice de área foliar, la cobertura y dinámica vegetal, la deforestación y fragmentación de hábitats, la evaluación del suelo y estado de la vegetación, así como también los cambios interanuales del NDVI relacionados con el fenómeno de El Niño en diferentes regiones y cambios a nivel global (Rodríguez y Bullock, 2018).
Estudios realizados sugieren que con un manejo adecuado de riego se influye de manera directa en la calidad del fruto, básicamente en la firmeza, ya que si el fruto es producido con un bajo déficit hídrico, este es más pequeño (Varela, Caldera, Zegbe, Mena, y Serna, 2014).
Esto ha motivado a que en la presente investigación se determinen los índices de NDVI, CWSI y SAVI del cultivo de la tuna en la zona antes señalada, para conocer el estado de salud de las plantaciones en una fecha determinada, aunque es un método simple, es un concepto importante para especificar el comportamiento de estos índices, combinados con el tipo y textura del suelo proporcionado por el Ministerio de Agricultura y Ganadería, se pueda implementar el sistema de riego más adecuado para la zona.
2. Metodología
En este estudio se analiza el comportamiento promedio del NDVI en las parcelas de tunas, cubre el área ubicada a 1°35′00″ S 78°38′00″ O, 2666 m.s..n. m., cerca de la ciudad de Riobamba se encuentran las parroquias Guano y San Andrés en el cantón Guano de la provincia de Chimborazo, en las comunidades de “Alacao, Carrera Ambato, Chingazo Alto y Bajo, Ela, La Merced, La Unión, Miraflores, San José de Chocón, San Andrés, San José de Juntuz y Santa Rosa” como se muestra la Figura 1, a partir de una imagen SENTINEL 2A que se obtuvo del Servicio Geológico de los Estados Unidos (USGS, por sus siglas en inglés) (https://earthexplorer.usgs.gov/) la información se encuentra en formato Geo TIFF, aplicando la corrección atmosférica se obtuvo una resolución espacial de 10 m con información de alta sensibilidad radiométrica en las 13 bandas espectrales (Casella, Barrionuevo, Pezzola, y Winschel, 2018).
Figura. 1:
Ubicación de las comunidades en el cantón Guano.
Las bandas en Sentinel 2A utilizadas para monitorear la vegetación son 12 (11), 8A, 4, donde la reflectancia puede variar dependiendo del contenido de la humedad en la región del infrarrojo de onda corta (SWIR, por sus siglas en inglés) en la hoja o en el suelo; lo cual indica que la vegetación vigorosa o irrigada y las áreas ribereñas muestran un color verde claro; es decir valores mayores del NDVI, mientras que la tierra seca y áreas naturales son de color verde opaco que indica valores menores del índice. (Rodríguez y Bullock, 2018). En general el cálculo del NDVI se define como:
[1]
Dónde: 
y 
son, respectivamente.
Los valores de reflectancia bidireccional de la superficie para las bandas del infrarrojo cercano y rojo, que son las bandas B2 y B1 de la imagen.(Salinas, Octavio, Martínez, y Morales, 2017).
El valor promedio del NDVI de cada parcela en el área de estudio se calculó con la siguiente formula:
[2]
Dónde: 
es el promedio de NDVI de las parcelas de cada comunidad, 
corresponde al NDVI de cada parcela que se encuentra dentro de la comunidad, y 
es el número de parcelas de la comunidad.(Salinas et al., 2017). En esta investigación se utilizó las bandas 8 y 4, obteniendo la ecuación:
[4]
Además, se analizó el CWSI, por ser considerado el índice verde más conocido y utilizado, en el que en forma general cuando se secan las hojas aumenta su reflectancia a lo largo de todo el espectro, produciendo un incremento en las bandas roja, en el infrarrojo medio y en las bandas de absorción de agua (Rock et al. 1985). Por lo que se puede definir como la determinación de la tasa de transpiración actual de un cultivo mediante la medición de la temperatura del dosel y el déficit de presión de vapor (López, Arteaga, Vásquez, López, y Sánchez, 2009).
El índice de estrés hídrico se calcula mediante la siguiente expresión:
[5]
Donde:
NIR= Infrarrojo cercano
SWIR = Infrarrojo medio
Aplicando las bandas correspondientes al infrarrojo cercano en infrarrojo medio que utiliza Sentinel 2. (Alderfasi y Nielsen, 2001), se expresa con la fórmula:
[6]
El SAVI, se aplica para que las alteraciones del valor del NDVI disminuyan en zonas áridas, debido a que puede corresponder a cubiertas vigorosas, pero poco densas o a cubiertas densas con poca vitalidad (Gonzaga, 2014). El SAVI se calcula mediante la siguiente expresión:
[7]
Donde:
IR = corresponde al infrarrojo cercano
R = corresponde al rojo
L = parámetro que varía según la densidad de la vegetación, para densidades
intermedias, para las imágenes de satélite, se toma 0.5.
El factor L (línea de suelo) es un parámetro que varía según la densidad de la vegetación y que está relacionado con la reflectividad del suelo. Para densidades intermedias, como en el caso de las imágenes de satélite, se toma L = 0.5. Este índice se utiliza con el fin de conseguir corregir algunos efectos externos al valor del NDVI en lugares donde el aporte de la reflectividad del suelo es elevada (Carvacho y Sánchez, 2010).
El SAVI se obtuvo con las bandas 8, 4 y con un valor L=0.5, aplicando la siguiente ecuación:
[8]
Finalmente, se analizaron los índices resultantes utilizando softwares libres y comerciales especializados.
La georreferenciación de las parcelas se realizó utilizando el GPS de la marca Garmin modelo GPSMAP 62sc con un receptor de alta sensibilidad que tiene un margen de error de ±4 metros, los datos tomados fueron verificados y validados antes del proceso cartográfico, combinando con la RGB (Red, Green, Blue) obtenida del Satélite Sentinel 2A.
3. Resultados y discusión
Como consecuencia de este estudio se obtuvo que el 10.73 % del total de cultivo de tuna que existe en el Ecuador está en la zona de estudio; esto es aproximadamente 19.31 ha repartidas en 66 parcelas a lo largo de las dos parroquias que se indica anteriormente de las 180 ha cultivadas en todo el territorio. Como muestra la Tabla 1.
Tabla. 1
Número de parcelas y extensión
Del muestreo obtenido en la imagen de las parcelas cultivadas de tuna se obtuvo los valores máximo y mínimo, es así que en comunidad la Unión se obtuvo un valor del NDVI de 0.40 y San Andrés con el 0.41, mientras que los valores más bajos están en las comunidades de Chingazo, Chocón y San José de Juntuz con un índice de 0.27; por lo tanto, el resultado obtenido se encuentra en el rango promedio de 0.27 a 0.41, evidentemente es un dato relativo de acuerdo a las circunstancias y la resolución de la imagen analizada.
Con estos antecedentes, el valor promedio obtenido del NDVI es de 0.31 de todas las parcelas muestreadas en el área. Como muestra la Figura 2.
Figura. 2:
Ejemplo de muestreo promedio de los valores de NDVI en las parcelas muestreadas
Se determinó la textura del suelo en la zona de cultivo de tuna, sobreponiendo las capas de cartografía del cantón Guano con la capa de textura de suelo proporcionado por el SNI del Ecuador, dando como resultado que el 99.9 % es gruesa, y un 0.1 % es fina. Como muestra la Figura 3.
Figura. 3:
Tipo de suelo en parcelas de tuna.
Además, se analizó el índice SAVI, determinando que para este tipo de textura (gruesa) la reflectancia recibida pertenece tanto al suelo como a la vegetación (tuna), mediante el cálculo de raster con las bandas 8 y 4 sumado el parámetro L con un valor de 0.5, se obtuvo valores referenciales aproximados bajos que pueden indicar que el cultivo de tuna tiene una cubierta de vegetación verde moderada como muestra la Figura 4:
[9]
Figura. 4:
Ejemplo de muestreo de valores SAVI.
En tanto que, los valores aproximados obtenidos de la tasa de transpiración después del cálculo raster entre las bandas 8 y 11 del CWSI muestran valores negativos, que pueden indicar que la pérdida de agua durante el día en la tuna es escasa. Como ilustra la Figura 5.
[10]
Figura. 5:
Ejemplo de muestreo de valores CWSI.
Por lo tanto, una de las posibles técnicas más adecuadas para incrementar la producción de la tuna es la de riego por goteo.
Para la georreferenciación de la zona, se tomaron los datos de las coordenadas (x, y, z), utilizando el GPS Garmin 62CS, con resolución de alta sensibilidad gracias a su antena Quadrifilar Helix, que da un margen de error de +/- 4 metros para la toma de cada una de las parcelas que varían desde 4 hasta 21 puntos dependiendo del tamaño del área cultivada. Los datos se importaron a tablas con la ayuda de una hoja de cálculo, en la que se registró la latitud, la longitud, la altura y la ubicación de la parcela en cada parroquia, para posteriormente importar en el software ArcMap 10.3 de ArcGis y con la ayuda de herramientas especializadas, se dibuja cada parcela sobre el mapa base de la textura de suelo.
El diseño de las tablas de registro de los datos contiene los siguientes campos: ID, X, Y, Z y UBICACIÓN, como muestra la Tabla 4.
Tabla. 4
Diseño de tablas de registro de datos.
Donde:
ID = número de punto
X, Y, Z = valores de latitud, longitud y altura
UBICACIÓN = nombre de la comunidad
Una vez importados los datos del GPS a la hoja de cálculo e importados a un programa comercial especializado, con la ayuda de las herramientas de importación y representación sobre el mapa base se obtuvo el resultado final como muestra la Figura 6.
Figura. 6:
Resultados de NDVI del AOI.
El mapa base utilizado es del servidor Maps de Microsoft, descargado utilizando el buscador Bing en el programa SasPlanet.
Sobre este tema existen publicaciones que han tratado comparaciones entre índices
NDVI, SAVI, TSAVI y NDII, en el contexto de los ecosistemas áridos y semiáridos en la parte central de la península de Baja California, demostrando que las variaciones espaciales y temporales no son paralelas, sin embargo, las diferencias entre NDVI y SAVI son relativamente menores (Rodríguez y Bullock, 2018).
En el artículo “Relaciones entre la precipitación, producción de biomasa e índices espectrales de la vegetación: alcances y limitaciones” publicado en marzo de 2018 modelan la relación entre el índice espectral NDVIcp de la biomasa con la precipitación total, usando modelos relativamente simples de procesos ecológicos y de manejo, desarrollando un modelo hibrido de aproximación (sensores remotos y modelos biofísicos) en el manejo de pastizales tropicales y bosques/selva (Paz y Díaz, 2018)
En el estudio “Multispectral aerial image processing system for precisión agriculture”, realizado en el año 2018, señalan que con la ayuda de ortofotos generaron los índices NDVI, GNDVI y SAVI sobre parcelas de caña de azúcar para identificar problemas de crecimiento en los surcos y tomar medidas en áreas puntuales de los cultivos.(Kharuf-Gutierrez, Orozco-Morales, Aday Díaz, y Pineda Ruiz, 2018). En esta investigación el estudio realizado de los cultivos de tuna se obtiene los índices NDVI, CWSI y SAVI con una imagen SENTINEL 2A, mejorando su resolución con correcciones radiométricas y atmosféricas, obteniendo resultados aproximados y referenciales, debido a que también existen otros factores importantes que influyen en el estado de salud del cultivo como son la fenología del cultivo, precipitación y temperatura, entre otros que no han sido considerados.
Una de las ventajas más importantes de este estudio es que se puede utilizar la reflectancia captada por los sensores de cualquier tipo de cultivo, en este caso de la tuna, para determinar los índices adecuados que permitan realizar un análisis profundo del comportamiento de la misma.
El procedimiento seguido en esta investigación deja abierta la posibilidad para que en el futuro se pruebe métodos alternos que conduzcan a nuevas conclusiones y sugerencias para el incremento de la producción de la tuna con el sistema de riego adecuado.
4. Conclusiones y recomendaciones
Los índices obtenidos en esta investigación son referenciales y aproximados y pueden ayudar al incremento de la producción de tuna de los agricultores del sector con la implementación de un sistema de riego de agua por goteo, considerando que en la actualidad se riega solo con agua de lluvia.
El índice de vegetación SAVI disminuye las alteraciones del valor del NDVI, influyendo en el estado vegetativo de un cultivo de tuna debido a que la zona es árida, considerando el valor de parámetro L para densidades intermedias.
El índice CWSI, obtenido mediante la medición de la temperatura del dosel y el déficit del valor de presión, indica que la tasa de transpiración es mínima en la plantación de tuna estudiada.
Para obtener valores más exactos de los índices en el cultivo de tuna es recomendable utilizar imágenes de menor resolución.
Las muestras han sido cultivadas en una parcela de en un área de 0.8 ha para realizar un seguimiento minucioso del crecimiento, numero de brotes, grosor, número de frutos, entre otros datos con la implementación de un sistema de riego por goteo.
