texto en 
Inglés (pdf)
Artículo en XML
Referencias del artículo
Traducción automática
Enviar artículo por email
Citado por SciELO 
Similares en
SciELO 
Permalink
Introducción
Según la Comisión Nacional del Agua (Conagua), el uso consuntivo del recurso hídrico por el sector agrícola es de 76% (Conagua, 2018). El sistema de producción agrícola es el usuario principal del agua, pero no necesariamente la utiliza de manera eficiente (Reyes-González, Kjaersgaard, Trooien, Hay, & Ahiablame, 2018). Los productores aplican más agua de riego de la que requiere el cultivo. Ante la perspectiva de mejorar el aprovechamiento del agua de riego y eficientar su uso, es esencial implementar medidas que estimen con mayor eficiencia la demanda hídrica de los cultivos.
El requerimiento hídrico de los cultivos es equivalente al valor de la evapotranspiración real del cultivo (ETc). Para su cálculo deben incluirse datos de las condiciones meteorológicas, de la capacidad de retención de humedad del suelo y de la etapa fenológica del cultivo. Un método útil para estimar la ETc es multiplicar la evapotranspiración de referencia (ETo) por un coeficiente de cultivo (Kc) (Allen, Pereira, Raes, & Smith, 1998). La ETo se estima con base en datos meteorológicos de temperatura del aire, velocidad del viento, radiación solar y humedad relativa. Dada la naturaleza cambiante de los factores que intervienen en su cálculo, el Kc puede ser tomado como referencia en la literatura (Doorenbos & Pruitt, 1977; Allen et al., 1998) o bien, calculado directamente en el terreno a través de sensores; estas observaciones se asociarían de forma directa con la etapa fenológica del cultivo.
Factores bióticos y abióticos que afectan el crecimiento y desarrollo del cultivo, como estrés hídrico, deficiencia-exceso de nutrientes, daño por plagas y enfermedades, entre otros, afectan el valor de Kc y lo desvían del valor de referencia publicado. Un número importante de autores sugieren que es mejor utilizar valores de Kc provenientes de índices de vegetación, en lugar de Kc de tablas, ya que éste representa la situación actual del cultivo (Bausch, 1993; Gontia & Tiwari, 2010; Lei & Yang, 2014; Adamala, Rajwade, & Reddy, 2016), y su uso como referencia puede mejorar la programación y empleo eficiente del agua (Kar, Kumar, & Martha, 2007).
El índice de vegetación de diferencia normalizada (Normalized Difference Vegetation Index-NDVI; Rouse, Haas, Schell, & Deering, 1974) se ha utilizado para el monitoreo de la vegetación, estimación de rendimiento de cultivos y detección de estrés (Kamble, Kilic, & Hubbard, 2013). En diversos trabajos de investigación se ha demostrado que los índices de vegetación (estimados con satélite o a través de dispositivos en campo) son empleados como base confiable para pronosticar con precisión los Kc de diferentes cultivos (Bausch & Neale, 1987; Jayanthi, Neale, & Wright, 2007; Irmak et al., 2011; Reyes-González et al., 2019). Se han incluido cultivos como alfalfa (Wright, 1982); trigo (Garatuza-Payán, Tamayo, Watts, & Rodríguez, 2003); frijol (Laike, Tilahum, K., & Hordofa, 2006); caña de azúcar (González-Dugo & Mateos, 2008); cebolla (Trout, 2008); sorgo (Singh & Irmak, 2009); uva (Campos, Neale, Calera, Balbontín, & Gonzalez-Piqueiras, 2010), y maíz (Reyes-González et al., 2019). Sin embargo, la relación entre el NDVIcampo proveniente del GreenSeeker con el Kc tomado del manual de la FAO-56 (Allen et al., 1998) en cártamo forrajero no se ha explorado.
El cártamo es una planta que crece cada año con un tallo principal, ramas primarias que pueden producir ramas secundarias y terciarias. Cada rama termina con una cabeza de flor (capítulo) provista de espinas (Dobrin & Marin, 2015). La altura de planta, número de ramas y presencia de espinas en las hojas depende de la variedad y del medio ambiente (Dajue & Mundel, 1996). El cártamo es tolerante a altas temperaturas, salinidad y sequía (Dwiedi, Upadhyaya, & Hegde, 2005; Tayebi, Afshari, Farahvash, Sinki, & Nezarat, 2012). Estas propiedades son importantes en el desarrollo de la planta y su rendimiento (Saruhan, Saglam, Demiralay, & Kadioğlu, 2012; Tayebi et al., 2012).
Investigaciones realizadas en Irán, Turquía, India, Rumania, Italia y Brasil, entre otros, se han enfocado en conocer las necesidades de diferentes niveles de nitrógeno (Arslan & Bayraktar, 2015; Bonfim-Silva Paludo, Sousa, de Freitas Sousa, & Da Silva, 2015); densidad de plantación (Mohammadi, Ghassemi-Golezani, Zehtab-Salmasi, & Nasrollahzaade, 2013; Vaghar, Shamsi, Kobraee, & Behrooz, 2014); niveles de riego (Istanbulluoglu, 2009; Jalali, Salehi, & Bahrani, 2011), y eficiencia en el uso del agua en cártamo (Kar et al., 2007; Lovelli, Perniola, Ferrara, & Di Tommaso, 2007). Sin embargo, en México existe poca o nula información que haya documentado obtener un Kc en cártamo forrajero para mejorar la programación del riego. Debido a lo anterior, los objetivos de la investigación fueron monitorear la altura y el índice de vegetación NDVI del cártamo durante dos ciclos de cultivo, desarrollar un nuevo Kc con base en el índice de vegetación para cártamo forrajero en la Comarca Lagunera, México, y calcular la ETc tomando en cuenta los nuevos valores de Kc.
Materiales y métodos
Área de estudio
La investigación se llevó a cabo durante los ciclos de otoño-invierno 2016-2017 y 2017-2018 en el establo E. y C. Asociados, en el Ejido Purísima, municipio de Matamoros, Coahuila, México (Figura 1). La preparación del terreno consistió en barbecho, rastreo y nivelación con láser. Se evaluaron tres cultivares de cártamo: uno con espinas (Guayalejo) y dos sin espinas (Selkino y CD-868). La fertilización se realizó sobre el terreno barbechado, utilizando una dosis de 150-50-00 unidades de N, P2O5 y K2O, respectivamente. Como fuente de N se usó un fertilizante nitrogenado con inhibidor de la nitrificación (3.4-Dimetilpirazol fosfato) y como fuente de P2O5 se empleó fosfato monoamónico. Los fertilizantes se incorporaron con un paso de rastra, para después nivelar el terreno. La parcela experimental fue de dos hectáreas.
Establecimiento del cultivo
La siembra se llevó a cabo con una sembradora de granos pequeños en un suelo seco el 10 y 12 de diciembre de 2016 y 2017, respectivamente, utilizando una densidad de siembra de 50 kg ha-1; tres días después, se aplicó el riego de siembra. En el ciclo 2016-2017 se hicieron tres riegos de auxilio a los 20, 37 y 58 días después de la siembra (DDS); y en el ciclo 2017-2018, dos riegos de auxilio a los 30 y 48 DDS. El manejo del cultivo estuvo a cargo del productor.
Medición del NDVI
El NDVI se empezó a medir a los 37 DDS en el ciclo 2016-2017 y a los 25 en el ciclo 2017-2018, a una altura de 0.6 m sobre el dosel del cultivo cada ocho días, pues es la frecuencia en que los satélites Landsat 7 y 8 pasan por el área de estudio. Se efectuaron cinco mediciones por día por cultivar en 10 m lineales en el centro de la parcela. Las mediciones fueron a las 10:00 a.m., con orientación de este a oeste, con el equipo GreenSeeker, diseñado por la Universidad Estatal de Oklahoma y comercializado por la compañía NTech Industries (Ukiah, CA, EUA). El GreenSeeker funciona dirigiendo un haz de luz roja (0.660 μm), así como un- infrarrojo cercano (0.780 μm) al dosel del cultivo. El campo de visión es un óvalo que va de 25 a 50 cm, dependiendo de la altura a la que se encuentre el GreenSeeker con respecto al dosel del cultivo. La altura de medición que se usó fue de 0.6 m, teniendo un campo de visión de 40 cm. El NDVIcampo obtenido del GreenSeeker para los tres cultivares de cártamo se comparó con el NDVIsatelital estimado por sensores remotos.
Imágenes de satélite
Se utilizaron cinco imágenes de satélite en el ciclo 2016-2017 y siete en el 2017-2018 para estimar el NDVI. Las imágenes se colectaron de los satélites Landsat 7 Enhanced Thematic Mapper plus (ETM+) y Landsat 8 Operational Land Imager (OLI) y Thermal Infrared Sensor (TIRS) (path = 30, row = 42) (Tabla 1). Cada imagen fue tratada por corrección radiométrica, atmosférica y superficie iluminada (Kjaersgaard & Allen, 2010). Las imágenes se seleccionaron con base en la cobertura temporal y las condiciones sin nubes. Las imágenes con nubes ubicadas a más de 10 km del área de interés se consideraron aceptables y fueron procesadas siguiendo las recomendaciones de Allen, Tasumi y Trezza (2007), y Kjaersgaard y Allen (2010), utilizando el modelo METRIC, a través del software ERDAS Imagine versión 14.0.
Tabla 1 La tabla muestra fecha, DDS y número de Landsat de las imágenes utilizadas en el estudio.
| Fecha de imagen | DDS | # Landsat |
|---|---|---|
| 19/01/17 | 37 | 8 |
| 27/01/17 | 45 | 7 |
| 04/02/17 | 53 | 8 |
| 20/02/17 | 69 | 8 |
| 08/03/17 | 85 | 8 |
| 06/01/18 | 25 | 8 |
| 14/01/18 | 33 | 7 |
| 22/01/18 | 41 | 8 |
| 30/01/18 | 49 | 7 |
| 23/02/18 | 73 | 8 |
| 03/03/18 | 81 | 7 |
| 11/03/18 | 89 | 8 |
El NDVI se estimó con el valor de reflectancia de las longitudes de onda rojo e infrarrojo cercano. Para el producto de Landsat 7 ETM+ se usaron las bandas 3 (0.63 a 0.69 μm) y 4 (0.77 a 0.90 μm), y para Landsat 8 OLI las bandas 4 (0.636 a 0.673 μm) y 5 (0.851 a 0.879 μm). Para evitar la mezcla de pixeles de NDVIsatelital se utilizó el pixel ubicado en el centro de la parcela, al igual que las mediciones con GreenSeeker. En último lugar se obtuvo un modelo de datos asimilados entre el NDVIsatelital y el NDVIcampo.
Altura de planta
Se midió con una cinta métrica desde la superficie del suelo hasta la última rama del cultivo, tomando 10 lecturas por cultivar, por día, y fueron promediadas. Las lecturas se registraron en la misma fecha en la que se tomó el NDVIcampo (ocho días), coincidiendo con el paso del satélite Landsat.
Relación entre NDVIcampo y Kc FAO-56
Se estableció la relación entre el NDVIcampo y el Kc real con base en el manual FAO-56 (Allen et al., 1998) para ambos ciclos. La ecuación general obtenida en los dos ciclos se usó para proponer un valor de Kc como una función de NDVI. Al final, los valores del nuevo Kc se multiplicaron por ETo para generar la evapotranspiración real del cultivo (ETc) y compararla con la evapotranspiración de referencia, la cual se tomó de la estación meteorológica del INIFAP, campo experimental La Laguna, al igual que la precipitación y las unidades calor.
Resultados y discusión
Precipitación y evapotranspiración
La precipitación acumulada fue de 25.6 mm para el ciclo 2016-2017 y de 38.4 mm para 2017-2018. La precipitación se observó con mejor distribución en el primer año de estudio, ya que los mayores eventos sucedieron al final del ciclo, después de los riegos de auxilio (Figura 2). Los valores diarios de evapotranspiración de referencia (ETo) variaron de 1.0 a 6.8 mm día-1 en el ciclo 2016-2017 y de 0.4 a 5.5 mm día-1 en el ciclo 2017-2018. Sin embargo, los valores acumulados durante el ciclo del cultivo (de siembra a cosecha) fueron similares (267 mm para el ciclo 2016-2017 y 266 mm para el ciclo 2017-2018). Como era de esperarse, el valor mínimo de ETo se registró al inicio de cada ciclo de producción y el valor máximo al final del mismo. Estos valores coincidieron con las temperaturas promedio anotadas en la Comarca Lagunera, donde las temperaturas frías (10 °C) son características del mes de enero y las templadas (20 °C) en marzo. Por otro lado, la Figura 3 muestra la acumulación de unidades calor en los dos ciclos de cultivo, registrándose 1 075 unidades calor en el primero y 1 050 en el segundo.
Figura 2 Precipitación, riego y ETo en dos ciclos de cultivo de cártamo forrajero en Matamoros, Coahuila.
Curvas de NDVI
Las curvas de NDVI mostraron estar asociadas con las etapas fenológicas del cultivo y sólo al final del ciclo de producción se notaron divergentes. Esto debido a que el último riego para el ciclo 2016-2017 se aplicó el 9 de febrero de 2017 (58 DDS), además de recibir una precipitación de 14.6 mm antes de la cosecha; mientras que en 2017-2018 el último riego fue el primero de febrero de 2018 (51 DDS) y recibió una lluvia de 7 mm antes de la cosecha. Así, los valores de NDVI en 2016-2017 se mantuvieron en 0.8 y los valores en 2017-2018 disminuyeron a 0.6 al final del ciclo (Figura 4). Sin embargo, Tasumi, Allen, Trezza y Wright (2005) reportaron que los valores de NDVI no son afectados por la humedad del suelo en un estudio que realizó en 383 campos agrícolas al sur de Idaho, EUA. En general, los valores de NDVI tomados en campo con el GreenSeeker para ambos ciclos fueron en promedio 0.2, 0.8 y 0.6 para la etapa inicial, madurez y senescencia, respectivamente. Valores similares de NDVI fueron reportados por Hajare, Patil y Verma (2008) durante dos ciclos en el cultivo de cártamo sembrado en tres tipos de suelo en Nagpur, India, aunque al final del ciclo caen hasta un valor de 0.2 cuando el cultivo se encontraba seco.
Figura 4 Curvas de NDVI tomados con el GreenSeeker para tres cultivares de cártamo. Los círculos rojos representan el NDVI derivado del satélite en dos ciclos de cultivo.
La discrepancia entre valores de NDVI, derivado del satélite y los valores tomados en campo con el GreenSeeker fueron mayores al inicio y al final del ciclo (Figura 4). Esta diferencia se debió a que las mediciones en campo son directamente al dosel del cultivo, cubriendo un área de alrededor de 0.4 m; mientras que con el satélite, las estimaciones son de un área de 30 x 30 m, midiendo follaje y suelo desnudo. Las mediciones en campo se realizaron el mismo día y a la misma hora en que el satélite pasa por el área de estudio. Las imágenes utilizadas fueron de días sin nubes.
Curvas similares de NDVI en trigo, caña de azúcar, algodón y maíz fueron reportadas por Garatuza-Payán et al. (2003); Zhang, Anderson, y Wang (2015); French, Hunsaker y Thorp (2015), y Reyes-González et al. (2019), respectivamente. Estos autores encontraron discrepancia en la primera y última etapa del cultivo entre las imágenes de satélite y las mediciones en campo. Los resultados sugieren que el periodo de mayor diferencia entre el satélite y el GreenSeeker se presenta cuando los valores de NDVI son menores a 0.6, pues con el satélite se recomienda utilizar el índice de vegetación ajustado al suelo (SAVI) (Huete, 1988), ya que éste minimiza la reflectancia del suelo en las primeras y últimas etapas del cultivo cuando hay escasa cobertura vegetal. Por otro lado, Reyes-González et al. (2018) documentaron que al inicio y término de la estación de crecimiento del maíz, el valor obtenido de ET a partir del NDVI fue menos preciso. Sin embargo, para propósitos de calendarios de riego, este método puede ser aceptable. En otro estudio, Anderson, Allen, Morse y Kustas (2012) reportaron que la ETc calculada a partir de índices de vegetación siempre fue subestimada en áreas irrigadas de Idaho, EUA.
En la Figura 5 se muestra el desarrollo de la cobertura vegetal del cártamo variedad CD-868 durante el ciclo 2017-2018, el cual empezó con valores bajos; después se va incrementando y al final disminuye ligeramente. Dicha figura presenta la estrecha relación entre la cobertura vegetal y el índice de vegetación.
Relación entre NDVIcampo y NDVIsatelital
Con las imágenes del satélite se obtuvieron cinco valores de NDVI en 2017-2018 y siete en 2017-2018. La Figura 6 expone la relación entre el NDVI adquirido en campo y el NDVI obtenido con las imágenes. Se encontró una buena relación entre lo medido y estimado, con un alto coeficiente de determinación (r 2 = 0.96). Este grado de confiabilidad muestra que el satélite es una herramienta útil, no sólo para estimar el NDVI a escala regional sino también a escala de campo. Altas correlaciones (r 2 = 90) de NDVI entre lo medido en campo con el radiómetro (MSR87) y estimado con sensores remotos fueron analizados por Duchemin et al. (2006) en 13 campos agrícolas en Marruecos, África.
Altura de planta
La Figura 7 indica la dinámica de crecimiento en altura de planta de los tres cultivares de cártamo evaluados durante los dos ciclos de estudio. A los 25 DDS, el cultivo mostró una altura de ~ 5 cm y en la etapa de roseta una altura de 10 cm (40 DDS). A continuación, a los 44 DDS, se inició la elongación del tallo hasta alcanzar la etapa reproductiva con el comienzo de la formación del capítulo a los 80 DDS, con una altura de planta entre los 90 y 100 cm. El tallo y ramas continúan creciendo junto con la diferenciación y desarrollo de las partes del capítulo que se hicieron distinguibles a los 93 DDS, aunque el capítulo aún no estaba abierto. En esta fecha, las plantas alcanzaron entre 111 y 151 cm de altura. Los cultivares Guayalejo y CD-868 presentaron la mayor altura de planta, y el cultivar Selkino mostró la menor en ambos ciclos. En trabajos anteriores se reportó que el cultivo presentó un ciclo de crecimiento entre 82 y 95 DDS, alcanzando alturas entre 82 y 93 cm (Reta et al., 2014; Reta et al., 2017). Esta comparación indica que las alturas de planta fueron superiores a las observadas en trabajos previos en la región. También fueron superiores a las alturas reportadas por Tayebi et al. (2012), quienes obtuvieron alturas máximas de 80 cm en cártamo establecido con diferentes niveles de evapotranspiración en Tabriz, Irán. Sin embargo, alturas similares (de 104 a 149 cm) fueron reportadas por Kizil, Cakmak, Kirici y Inan (2008) en cultivares de cártamo establecidos en Eskishehir, Turquía. No obstante, Bagheri, Saeedi y Zaeh (2006) mencionan que la humedad disponible en el suelo afectó la altura de planta en un estudio realizado en Isfahán, Irán.
Por otro lado, la Figura 8 señala una relación polinomial de segundo orden entre la altura de planta y el NDVI de los cultivares de cártamo, lo que explica entre 64 y 92% la variación total observada del NDVI en el primer ciclo (2016-2017), y entre 83 y 87% en el segundo ciclo (2017-2018).
Relación entre NDVIcampo y Kc FAO-56
Los valores de NDVI se tomaron de tres cultivares de cártamo evaluados durante dos ciclos, mientras que los valores de Kc real se obtuvieron del manual de la FAO-56. La Figura 9 muestra la relación entre el NDVIcampo de tres cultivares de cártamo y el Kc tabulado para dos estaciones de crecimiento. En el primer ciclo de cultivo la relación entre el NDVIcampo y el Kc mostró una r 2 = 0.93; para el segundo ciclo de 0.89 indicó una buena relación entre NDVIcampo y Kc. Sin embargo, en el segundo ciclo, el valor de “r” disminuyó ligeramente, esto se debió a que el cultivo presentó condiciones de estrés por causa de la mala programación de los riegos, pues sólo se aplicaron dos riegos de auxilio a la mala distribución de la lluvia durante el ciclo de cultivo. Por tal motivo, la relación entre NDVIcampo y Kc fue ligeramente baja, ya que los Kc de la FAO están desarrollados para cultivos bien regados y sin restricciones de humedad y nutrimentos (Allen et al., 1998). No obstante, los valores de Kc FAO-56 después de los 50 DDS fue de 1.0, mientras que los valores de NDVIcampo oscilaron entre 0.65 y 0.8.
Figura 9 Relación entre el NDVIcampo de tres cultivares de cártamo forrajero y el Kc de la FAO-56 en dos estaciones de crecimiento.
Se han reportado altas correlaciones entre NDVIcampo y Kc, tomado de tablas de la FAO para diversos cultivos. Rocha, Perdigão, Melo y Henriques (2012), así como Reyes-González et al. (2016) reportaron altos coeficientes de determinación (r 2 = 0.99), para el cultivo de maíz y alfalfa. Consoli, Russo y Snyder (2006); Singh e Irmak (2009); Gontia y Tiwari (2010); Kukal, Irmak y Kilic (2017), y Palacios-Vélez, Palacios-Sánchez y Espinosa-Espinosa (2018) reportaron buenas relaciones entre NDVI y Kc para el cultivo de naranja (r 2 = 0.80), sorgo (r 2 = 0.87), trigo (r 2 = 0.80), soya (r 2 = 0.81) y uva (r 2 = 0.83), respectivamente. Las diferencias entre altas y bajas relaciones se dan principalmente por el manejo de cultivo.
Desarrollo de la ecuación general de Kc
El NDVI tomado en campo y el Kc obtenido de la FAO-56 se usó para desarrollar ecuaciones de regresión lineal para dos estaciones de crecimiento (Figura 9). Las ecuaciones obtenidas para ambos ciclos fueron comparadas utilizando el método de prueba de “t” para probar la diferencia estadística entre dos regresiones independientes. Una vez realizada la prueba, el valor de “t” calculada fue menor (0.76) que el valor de “t” tabulada (1.68), lo cual indica que no hubo diferencia estadística significativa entre las dos ecuaciones. Con base en este resultado, todos los datos (dos ciclos) se agruparon y se generó una nueva ecuación para el cultivo de cártamo en la Comarca Lagunera (Figura 10).
Se han generado diversas ecuaciones a partir de índices de vegetación para diferentes cultivos. Por ejemplo, Consoli et al. (2006) desarrollaron una ecuación utilizando el NDVI en el cultivo de naranja en Italia (Kc = 1.2562 * NDVI + 0.287); Singh e Irmak (2009), para soya en EUA (Kc = 1.217 * NDVI - 0.034); Campos et al. (2010), para uva en España (Kc = 1.44 * NDVI - 0.10); Gontia y Tiwari (2010) para trigo en la India (Kc = 2.7109 * NDVI - 0.424); Rocha et al. (2012) para maíz en Portugal (Kc = 1.37 * NDVI - 0.017), y Reyes-González, Hay, Kjaersgaard, y Neale (2015) para alfalfa en México (Kc = 2.112 * NDVI - 0.4989). Estas ecuaciones se generaron para desarrollar nuevos valores de Kc en condiciones locales, como lo recomiendan Doorenbos y Pruitt (1977), y Allen et al. (1998), con el objetivo de mejorar el manejo y uso del agua en la agricultura.
La ecuación generada (Kc = 1.0644 * NDVI + 0.2128), en la Figura 10 se empleó para desarrollar una nueva curva de Kc para el cártamo forrajero, tomando en consideración los días después de la siembra (Figura 11). Valores similares de Kc para cártamo las reportaron Kar et al. (2007) en Dhenkanal, India, quienes informaron valores de 0.37 en la etapa inicial, 0.94 en desarrollo, 1.29 en madurez y 0.35 al final de la estación de crecimiento. Sin embargo, en las dos últimas etapas (madurez y final) hubo una discrepancia con los encontrados en este estudio, lo cual fue más notable al final del ciclo, ya que el cártamo evaluado en la Comarca Lagunera se utiliza para forraje y se cosecha en verde, por lo que no alcanza la madurez fisiológica, como el cártamo que se utiliza para semilla.
Comparación entre ETo y ETc
Los valores de ETo se tomaron de la estación meteorológica del campo experimental La Laguna, localizado a 15 km del área de estudio, mientras que los valores de ETc son resultado de multiplicar los valores nuevos de Kc por la ETo. La Figura 12 indica la comparación de ETo y ETc para el ciclo 2017-2018. Al inicio del ciclo (1-20 DDS), los valores de ETo fueron alrededor de 2 mm día-1, mientras que para la ETc fueron 0.23 mm día-1. En la etapa de desarrollo (21-57 DDS), los valores de ETo fluctuaron en 2.5 mm día-1 y los valores de ETc fueron de 0.3 a 2.5 mm día-1. En la etapa de madurez (58 a 85 DDS), los valores entre ETo y ETc fueron similares (3.5 mm día-1). Al final del ciclo (86-95 DDS), la ETo se mantuvo en 4.5 mm día-1 y la ETc disminuyó a 3.0 mm día-1. En general, la figura ratifica que, a principios del ciclo de cultivo, la ETo diaria fue mayor a la ETc y a mediados de ciclo los valores fueron similares; esto indica que con la ETc, la demanda hídrica del cultivo es acorde con las etapas fenológicas del cultivo. Valores de ETo mayores (6.0 mm día-1) a los encontrados en el estudio realizado por los autores fueron reportados por Lovelli et al. (2007) y Kar et al. (2007). Sin embargo, la ETo de diferentes variedades de cártamo, bajo diferentes condiciones de clima y suelo pueden variar de 200 a 1 000 mm por ciclo (Istanbulluoglu, 2009).
Entendiendo las diferentes etapas de desarrollo del cultivo y aplicando la cantidad exacta de volumen de agua, los productores podrán mejorar el calendario de riego, el manejo del agua e incrementar la sustentabilidad del recurso hídrico en los cultivos. Reyes-González et al. (2016) informaron que los agricultores deberían usar la ETc en lugar de la ETo para la programación de los riegos en regiones áridas y semiáridas donde el agua es escasa.
Conclusiones
Los valores de NDVI y de altura fueron afectados por la calendarización de los riegos y por la precipitación recibida durante las dos estaciones de crecimiento. Por otro lado, la relación entre el NDVIcampo y el Kc de FAO-56 fue buena, aunque en el segundo ciclo de crecimiento bajo, a consecuencia de la mala programación de los riegos y la distribución de la lluvia, pues los Kc de la FAO se desarrollaron para cultivos sin restricciones de humedad.
La ecuación generada entre NDVIcampo y Kc FAO-56 para ambos ciclos de cultivo se utilizó para desarrollar nuevos valores de Kc para cártamo forrajero. Ello significó que los índices de vegetación son una herramienta útil para generar Kc en el nivel local. Los valores nuevos de Kc obtenidos provocarán que se aplique el agua de acuerdo con las diferentes etapas fenológicas y/o requerimientos hídricos del cultivo, ayudando a mejorar la programación del riego y haciendo más eficiente el uso del agua en la agricultura.
