Introducción

En el proceso de la producción agrícola, con especial énfasis en la que está bajo riego, varios factores afectan su rendimiento, entre los cuales se pueden mencionar tres grupos que impactan la producción y los rendimientos de los cultivos.

El primer grupo lo constituyen los factores meteorológicos, como las heladas, lluvias torrenciales, cambios en las temperaturas, ciclones, etcétera; sin embargo, es muy difícil controlarlos y en la mayoría de los casos no es posible hacerlo; no obstante, se puede medir su efecto y estimar las posibles afectaciones, así como llevar a cabo algunas acciones para aminorar su impacto negativo.

Las estaciones meteorológicas automáticas son dispositivos con sensores que tienen la capacidad de registrar y recolectar datos de las variables ambientales, miden los principales factores que afectan la productividad de los cultivos, también incluyen un software para realizar cálculos automáticos que procesan la información y la vuelven disponible para el productor agrícola. Entre la información que proporcionan son de mucho interés las temperaturas máximas y mínimas, la precipitación y la evapotranspiración de referencia a nivel diario. En el Distrito de Riego Río Mayo se tiene una parte de la Red de Estaciones Meteorológicas Automáticas de Sonora (REMAS). Un par de estaciones situadas en el centroide del distrito de riego utilizan datos diarios como base de medición del efecto climático. Se pueden estimar los efectos térmicos en los cultivos y la demanda evapotranspirativa que afecta procesos como la fotosíntesis, base de la vida en el planeta.

El segundo grupo son los factores relacionados con el manejo de los cultivos y del suelo, que también son de importancia para la obtención de buenos rendimientos. Un buen agricultor sabe cómo y cuándo llevar a cabo las actividades que le permitan lograr buenos rendimientos, sobre todo si dispone de los medios económicos y mecánicos para realizarlas. No obstante, es posible que mejore la productividad del agua y el suelo si puede monitorear sus cultivos usando tecnologías relativamente nuevas, como el uso de sensores remotos, que es una herramienta eficiente para la detección de diferencias espaciales y temporales de los factores que impactan la productividad de los cultivos y, de ser posible, su corrección. También los sensores remotos mediante el uso de satélites son una herramienta de gran importancia para medir la variabilidad espacial del suelo en el desarrollo de los cultivos, con objeto de efectuar las correcciones necesarias para aminorar esta variación y sus efectos perniciosos.

Desde 2006, el Colegio de Postgraduados (Colpos) fue invitado a unirse al proyecto “Participatory Multi-Level EO-Assisted Tools for Irrigation Water Management and Agricultural Decision-Support” (PLEIADeS), financiado por la Comisión Europea. Dentro de él, se ha dado seguimiento al desarrollo de cultivos en los distritos de riego 038 Río Mayo, 041 Río Yaqui y 051 Costa de Hermosillo, utilizando sensores remotos. En especial, en el distrito del Río Mayo se ha monitoreado el desarrollo de los cultivos mediante índices de vegetación, principalmente al cultivo del trigo, que representa alrededor de un 70% de la superficie cultivada, apoyado con la participación de productores cooperantes. En el presente trabajo se hace un resumen de los resultados obtenidos, haciendo énfasis en años recientes cuando se han observado impactos importantes, como consecuencia de cambios en el clima.

Materiales y métodos

El Distrito de Riego 038 Río Mayo se encuentra ubicado en el noroeste del país, al sur del estado de Sonora, entre las latitudes norte 26.71 y 27.25 grados y longitudes oeste 109.8 y 109.37 grados, con una superficie regable de 98 520 ha; la principal fuente de abastecimiento de agua es la presa “Adolfo Ruiz Cortines”, con una capacidad del orden de los 1 100 hectómetros cúbicos, además de 130 pozos profundos.

Para darle seguimiento a los cultivos, se han utilizado principalmente imágenes de los satélites Landsat 5, 7 y 8, además de los satélites Spot 5 y 6, Deimos y RapidEye.

El distrito de riego es operado por la Sociedad de Responsabilidad Limitada (S. de R. L.), constituida por 16 módulos de riego y 11 582 usuarios registrados.

Como se ha indicado, para evaluar los aspectos meteorológicos, se utilizaron datos diarios de dos estaciones meteorológicas, SEMAY de la REMAS y Etchojoa.

El tercer grupo es la variabilidad espacial. Para el monitoreo espacial del suelo y temporal del desarrollo fenológico de los cultivos se usaron índices de vegetación, como: índice de vegetación de diferencias normalizadas (NDVI, por sus siglas en inglés) (^{Rouse, Hass, Schell, & Deering, 1974}) y el índice de estrés de humedad (MSI, por sus siglas en inglés) (^{Rock, Vugelmann, Williams, Voglemann, & Hoshizaki, 1986}).

El NDVI está constituido por las reflectancias en las bandas roja e infrarroja cercana de acuerdo con la siguiente fórmula:

Donde:

Este índice varía entre -1 y 1. Los valores negativos por lo general se presentan donde hay agua; valores menores que 0.2 casi siempre se muestran cuando hay suelo desnudo; para la vegetación verde, usualmente se tienen valores que varían entre 0.25 hasta valores cercanos a 1 cuando la vegetación está en pleno desarrollo.

El MSI es la relación entre las bandas del infrarrojo medio y el infrarrojo cercano, como se muestra a continuación:

Donde:

Los que varían de 0.25 hasta 0.5 indican que los cultivos no están estresados; conforme aumentan crece el estrés hasta llegar a 1 cuando los cultivos alcanzan el porcentaje de marchitamiento permanente.

Desde hace varios años se ha encontrado una relación lineal entre los valores del índice NDVI y el rendimiento del cultivo del trigo. Así, el Instituto Internacional para el Manejo del Agua (IWMI, por sus siglas en inglés) llevó a cabo un trabajo en la India que presentó en su Reporte 28 denominado “Performance Evaluation of the Bhakra Irrigation System, India, using Remote Sensing and GIS Techniques” (^{Sakthvadivel, Thiruvengadachari, Amerasinge, Bastiaanssen, & Molden, 1999}), donde se encontró la relación entre el NDVI y el rendimiento del trigo, utilizando 274 parcelas de prueba y donde se presenta la gráfica y función encontrada, con lo que se clasificó la productividad de las diferentes secciones de este distrito de un millón y medio de hectáreas, que se muestra en la Figura 1.

Figura 1 Función trigo y distrito de riego Bhakra. 

En México se ha encontrado una relación semejante en el Distrito de Riego 038 del Río Mayo; así, en el artículo “Earth Observation as a Support to Improve Water Use in Irrigated Agriculture” (^{Palacios, Palacios, Rodríguez, & Palacios, 2010}) se muestra una relación entre el rendimiento observado en 26 parcelas de productores de trigo y sus valores promedio del NDVI detectados durante el desarrollo de este cultivo, similar al obtenido en la India, como se muestra en la Figura 2.

Figura 2 Relación NDVI medio con el rendimiento del trigo. 

También en este artículo se muestra la forma en que se ha podido estimar la evapotranspiración del cultivo del trigo utilizando la relación ETc = Kc ETr, donde ETc es la evapotranspiración del cultivo, ETr es la evapotranspiración de referencia, que puede calcularse mediante la fórmula de Penman-Monteith y Kc coeficiente del cultivo, y puede estimarse en función de los valores del NDVI, como lo han demostrado varios autores, como ^{Choudhry, Ahmed, Idso, Reginato y Daughtry (1994)}. Para el cultivo del trigo, ^{Calera y González-Piqueras (2007)}, en particular, encontraron una relación lineal entre Kc y el valor del NDVI obtenido de cada imagen de los satélites Landsat, cuyo valor es:

En el mismo artículo de ^{Palacios y colaboradores (Palacios et al., 2010)} se muestra una relación similar obtenida en el Distrito de Riego 038 del Río Mayo, que es:

La evapotranspiración del cultivo estimada con esta función, para cada una de las 26 parcelas, se relacionó con el rendimiento obtenido y se obtuvo una función con un ajuste estadístico bueno, como se muestra en la Figura 3.

Finalmente, en el mencionado artículo se muestra una relación no lineal entre el MSI y el rendimiento del trigo, que es:

Con R ² = 0.794 y un error estándar de 0.443 t/ha. Como es de esperarse, la función es inversa; de manera que conforme aumenta el MSI, disminuye el rendimiento. Es de mencionarse que valores cercanos a la unidad de este índice indican que las plantas han llegado al porcentaje de marchitamiento permanente.

En los años agrícolas 2013-2014 y 2014-2015 se presentaron condiciones meteorológicas poco favorables para el cultivo del trigo; las horas frío disminuyeron, sobre todo en 2014-2015, y en consecuencia las unidades calor no fueron suficientes para lograr rendimientos buenos; en este año agrícola, el rendimiento tuvo una reducción cercana a 30%, como se mostrará.

Se ha observado que los cultivos requieren ciertas cantidades de energía para lograr cada etapa de su desarrollo; estas cantidades de energía se estiman de acuerdo con los grados-día acumulados. También se ha observado que hay una relación entre la variación del índice NDVI y los grados días acumulados (^{Miller, Lanier, & Brandt, 2001}), como se observa en la Figura 4, donde en las ordenadas se muestran los grados días acumulados (en inglés, Growing Degree Days GDD), y en las abscisas la variación de los índices NDVI (línea gruesa) y la relación simple:

Donde R y NIR ya fueron definidos.

Figura 4 Relación grados día-NDVI y relación simple. 

El acumulado de la curva de desarrollo del NDVI es una sigmoide, es decir, una curva en forma de una S, la cual se puede ajustar a una función matemática como la logística. Como ejemplo pueden verse en la Figura 5 los valores acumulados del NDVI de dos parcelas del módulo 15, donde se observa que en la parcela 9173 el valor acumulado del índice quedó abajo del de la parcela 9384; en efecto, en la primera se obtuvo un rendimiento de 6.8 t/ha de trigo y en la segunda su rendimiento fue de 7.1 t/ha. También en la Figura 5 se muestra en ambos casos en línea continua el ajuste a una función logística.

Figura 5 Valores acumulados del NDVI en dos parcelas. 

Resultados y discusión

Los efectos del manejo

Respecto a la posibilidad de detectar problemas durante el crecimiento de los cultivos, una ayuda importante para tomar decisiones que permitan prevenir reducciones en el rendimiento será utilizar el visor (^{Espinosa, 2013}), que está disponible en Internet, mediante el cual el productor puede apreciar cómo va el desarrollo de su cultivo al observar la coloración que genera en NDVI en las parcelas, observando por una parte la homogeneidad del crecimiento del cultivo, y el acumulado gráfico del NDVI, que puede compararse con el obtenido en el año anterior, mediante lo cual podrá tener una idea sobre el rendimiento esperado; así, como se muestra en la Figura 5, si la curva que se va generando es más baja que la del año anterior, bajo el supuesto de que se tenga un cultivo similar, se tendrá la sospecha de que algo no está bien; esto permitirá al productor buscar las razones de la reducción esperada, con objeto de corregir los posibles problemas, como se muestra en la Figura 6.

Para ver el Sistema de Monitoreo Satelital del Distrito de Riego 038 Río Mayo, hay que entrar a la página de Internet: http://hidro.colpos.mx:8080/sig-mon/

En esta página, el usuario no necesita registrarse sino poner donde se pregunta el Nombre: demo; el Password: demo; en tipo de usuario: user, y luego dar clic para entrar. Al ingresar, se despliega una lista de los números de los módulos; suponiendo que la parcela esté en el módulo 15, se busca en la lista este número y se le da un clic; en el centro de la pantalla aparecerán las parcelas de dicho módulo. Para encontrar la parcela es necesario dar el número de cuenta del padrón de usuarios o buscar la parcela por su posición dentro del módulo (^{Espinosa-Espinosa, Palacios-Vélez, Tijerina-Chávez, Flores-Magdaleno, & Quevedo-Nolasco, 2017}).

En general, hay una presentación que se ha hecho del conocimiento de los usuarios para buscar sus parcelas, pero lo interesante es que se puede ver con detalle el estado de las parcelas por el color del índice NDVI, no solamente para la fecha más reciente de cuando ha pasado el satélite, sino que además se pueden apreciar los valores acumulados del índice en años anteriores y el avance del año actual, como se muestra en la Figura 6.

Figura 6 Visor del sistema de monitoreo. 

Los efectos de la variabilidad

Otro factor que se ha estudiado ha sido el efecto de la variabilidad del suelo en el rendimiento de los cultivos, la que afecta los rendimientos, no solamente en el trigo sino en general en todos los cultivos establecidos en el distrito de riego.

Es común que exista una variabilidad espacial en las parcelas. A veces esta variabilidad es muy significativa y por lo mismo afecta en forma negativa al rendimiento de los cultivos que se establezcan; por esta razón se ha desarrollado la denominada agricultura de precisión. Sin embargo, es factible detectar esta variabilidad mediante las imágenes satelitales, observado la variabilidad del valor de los píxeles.

Por otra parte, es muy común que en los distritos de riego de México se compren los derechos de agua de parcelas pequeñas; es decir, un productor puede rentar varias parcelas y en estos casos la variabilidad espacial puede aumentar, en especial por la variabilidad que existe en colindancia entre estas parcelas y además la variabilidad dentro de cada parcela.

Como ejemplo se presentará un caso hipotético de un lote de 27 parcelas de unas 5 ha cada una dentro del módulo de riego número 15 del Distrito de Riego 038 del Río Mayo, comprendido por más de seis mil hectáreas, donde es muy común que sus parcelas ejidales sean de 5 ha en promedio y se renten en extensiones de consideración, que en algunos casos llegan a ser del orden de mil hectáreas.

Para el ejemplo se seleccionó un lote de parcelas en la parte central del módulo y utilizando la imagen del satélite Landsat 8, correspondiente al 14 de febrero de 2016, cuando se alcanzó el mayor valor promedio del índice NDVI, se evaluó la variabilidad espacial de esta superficie sembrada con trigo.

El valor promedio del índice NDVI en las 27 parcelas fue bastante aceptable, por lo que se logró un buen rendimiento; sin embargo, podría haberse obtenido un rendimiento mayor si no se hubiera tenido la variabilidad espacial, principalmente por las regiones límite entre parcelas, como se puede observar en la Figura 7, que del lado derecho muestra las curvas de isonivel del índice NDVI. En esta misma figura se presenta una gráfica donde se señala la superficie de 140 hectáreas dividida en 19 barras, mostrando cada una de ellas el área que le corresponde. El valor promedio de NDVI de todas las barras fue de 0.877, con una desviación estándar de 0.068, por lo que su coeficiente de variación fue de 7.78%, éste es un índice de variabilidad.

Figura 7 Variabilidad del índice NDVI en 27 parcelas. 

Conclusiones y recomendaciones

De los tres grupos de factores que influyen en el rendimiento de los cultivos, los componentes meteorológicos o del clima tienen una influencia notable en los rendimientos de los cultivos; sin embargo, es poco lo que puede hacerse para modificar su efecto negativo; no obstante, el seguimiento de su evolución mediante su medición con estaciones meteorológicas automáticas permite hacer algunas previsiones para moderar sus impactos económicos y evaluar el impacto negativo en la economía del productor, con objeto de que pueda tomar decisiones anticipadas para aminorar las posibles pérdidas.

Respecto a los factores de manejo, el conocimiento del monitoreo espacial y temporal del desarrollo fenológico del cultivo mediante el NDVI y el MSI le permitirá al productor conocer las condiciones en que se encuentran mediante la visión del estado de las parcelas, consultando al visor en la página Web, por la variación de la coloración generada por los cambios del NDVI, así como la posición relativa del valor obtenido de este índice cada vez que pasan los satélites sobre el distrito de riego, y que puede compararse con la posición debida al desarrollo en años anteriores; esto le permitirá al productor orientar sus decisiones, para tratar de mejorar, si es necesario, las condiciones del cultivo.

Una condición crítica en el desarrollo de los cultivos, generada por un estrés debido al ataque de una plaga, enfermedad o condición hídrica, se podrá notar por la posición relativa del valor del NDVI respecto a la del paso anterior del satélite, lo que permitirá llevar a cabo las acciones necesarias para tratar de aminorar los daños al cultivo.

El efecto de la variabilidad espacial tanto en una parcela como en un grupo de ellas se hace notar por el cambio del color del índice. Además, como se puede observar en la Figura 7, los efectos de orilla y de la condición diferente de algunas de las parcelas del grupo, manejado por un productor, puede afectar el futuro rendimiento del cultivo. Para aminorar este efecto, será conveniente, si la renta de las parcelas es por varios años, tratar de reducir esta variabilidad mediante varias acciones, como la nivelación de los terrenos, la aplicación diferenciada de fertilizantes y la forma de regar.

Los apoyos y recomendaciones que puede brindar el personal técnico que han contratado algunos módulos de los distritos de riego serán de mucha importancia en el asesoramiento a los productores para lograr mejores rendimientos. También este personal podrá asistirlos en el uso del visor.