Espanhol (pdf)
Artigo em XML
Referências do artigo
Tradução automática
Enviar este artigo por email
Citado por SciELO 
Similares em
SciELO 
Permalink
1. Introducción
El arroz es uno de los principales cultivos de cerea les en el mundo (Fornasiero et al., 2022; Rahman et al., 2022), que mantiene la demanda nutricional y la seguridad alimentaria del 50 % de la población mundial (Jiang et al., 2019; Rahman et al., 2022). Anualmente son sembradas aproximadamente 160 millones de hectáreas que producen 750 millones de toneladas de arroz, previéndose un aumento del 28% de la demanda al 2050 (Jiang et al., 2019). Recientemente, el incremento de la población mundial ocasionó un aumento en la demanda de agua para la producción de alimento (Ouda & Zohry, 2022). Además, la Organización de las Naciones Unidas para Alimentación y la Agricultura (FAO), estimó la escasez del suministro de agua que presentan al menos 4 mil millones de personas en el mundo durante un mes al año (Hussain et al., 2022). Para 2025, la FAO ha proyectado que el 20% de los campos de arroz irrigados presentarán escasez de agua y además el 29% puede padecer de escasez económica de agua (Ishfaq et al., 2020). Por otro lado, el cambio climático está ligado directamente al agua, provocando alteraciones en el ciclo hidrológico, causando una distribución irregular de las precipitaciones y el aumento de la temperatura es diferentes partes del mundo, ocasionando mayores pérdidas en evaporación y evapotranspiración de los cultivos, reflejando un aumento de la necesidad de los recursos hídricos, siendo esto critico ante escenarios futuros de sequía (Roy et al., 2022; Ouda & Zohry, 2022; Zobeidi et al., 2022).
El Perú se encuentra entre uno de los países que presenta escasez de agua media a alta (20% - 40%), según la relación de extracción y disponibilidad del recurso hídrico en el mundo (Hussain et al., 2022). Sin embargo, uno de sus principales cultivos demandantes de agua es el arroz, con un consumo de 12000 m3 ha-1 y 20000 m3 ha-1 en la costa del Perú (Neira et al., 2020). Además, el consumo per cápita anual por persona es de 63,5 kg de arroz, se siembran 4175970 m2 en las que se producen 3436637 t, con un rendimiento promedio de 0,00082 t m-2 (FAOSTAT, 2023). En este sentido, la costa norte de Perú (Piura, Lambayeque y La Libertad), dominadas por planicies costeras áridas, aproximadamente se siembran 715000 m2 con una producción de 1017267 t para el periodo del 2021 (FAOSTAT, 2023).
El cambio climático está afectando la producción del arroz, debido al aumento de calor (alta temperatura y humedad), de los periodos de sequía, y de la salinidad (Ouda & Zohry, 2022; Zobeidi et al., 2022). La necesidad de ahorro de agua en arrozales ha fomentado que se investiguen variedades que necesiten un aporte hídrico menor y además optando planificaciones prácticas para un uso más eficiente del agua, como la técnica de riego, alternancia de humedecimiento y secado (AWD), evapotranspiración (ET), trasplante en suelo no encharcado, cultivo en suelo saturado, sistema de camas elevadas, sistema de arroz aeróbico, siembra directa en seco (Ishfaq et al., 2020; Cheng et al., 2022). La técnica AWD se basa en la alternancia de fases de inundación y secado del suelo a partir las etapas de macollamiento hasta de llenado de grano (Mote et al., 2022; Cheng et al., 2022) e implica un constante monitoreo del estado del cultivo durante la fase de secado, ya que la falta de riego puede provocar estrés hídrico en la planta, reduciendo de rendimiento (De la Cruz et al., 2022; Mote et al., 2022; Gao et al., 2023). Se han realizado muchos esfuerzos para promover la AWD. Adoptado principalmente por los países asiáticos. Por ejemplo, China posee el 41% de sus campos de arrozales irrigados por AWD, luego se encuentra Tailandia, Bangladesh, Vietnam y Filipinas (Chu et al., 2018; Ishfaq et al., 2020; Cheng et al., 2022; Yu et al., 2023).
La ET estima el requerimiento de agua en los cultivos; puede obtenerse utilizando ecuaciones empíricas o algoritmos de detección remota, además de métodos tradicionales de medición como los lisímetros, sensores de humedad y torres de flujo turbulento. Sin embargo, poseen un elevado costo, adquisición e instalación y ofrecen información puntual a escala local y para un cultivo en específico (Niu et al., 2020; Saha et al., 2022). Con la implementación de los sensores remotos, se dispone de información continua de diversas imágenes dentro del rango del espectro electromagnético, que permite caracterizar la fenología del cultivo, suelo, agua y entre otros (Saha et al., 2022). El uso de sensores remotos y el balance de energía en la superficie permite estimar la ET.
El modelo METRICTM (Mapping Evapotranspiration at High Resolution using Internalized Calibration) (Allen et al., 2007), calcula la ET como un residual en el balance de la energía de la superficie y es una variante del SEBAL (Surface Energy Balance Algorithm for Land), que es un algoritmo basado en un proceso de balance de energía (Bastiaanssen et al., 1998). Inicialmente utilizado con información de satélite de bajo costo; sin embargo, la principal limitación es la falta de re solución espacial para estudios a nivel parcelario, la presencia de nubosidad atmosférica e insuficientes datos en tiempo real (Niu et al., 2020). Por lo tanto, genera incertidumbre en la cuantificación de ciertas variables de la ecuación del balance energético.
En consecuencia, se ha desarrollado una nueva tecnología como la utilización de Vehículos Aéreos No Tripulados (VANT) a la necesidad de obtener imágenes de alta resolución, sin errores atmosféricos y adquisición de imágenes en tiempo real; permitiendo la monitorización de cultivos homogéneos; sin embargo, el procesamiento de las imágenes requiere una mayor capacidad computacional que requiere ser evaluada (Pintér et al., 2022). Con esta tecnología se han realizado investigaciones en diferentes cultivos alrededor del mundo, por ejemplo: viñedos en Estados Unidos (Nassar et al., 2021), cebada en Dinamarca (Hoffmann et al., 2017), olivo en Chile (Ortega-Farías et al., 2016), olivo en España (Ramírez-Cuesta et al., 2019) y entre otros.
Esta investigación tuvo como objetivo pro poner un sistema del riego con menor consumo de agua en arrozales a partir de la estimación de la evapotranspiración obtenida de imágenes deriva das de un VANT, para proporcionar una guía en la gestión adecuada del riego en la costa norte del Perú.
2. Metodología
2.1 Área de estudio
El estudio se enmarca en un área total de 12400 m2 con arrozales de la variedad INIA 508- Tinajones con un marco de plantación de 0,20 m x 0,20 m y tres condiciones de riego: (i) riego por inundación continua (IC) (5900 m2), ii) riego con alternancia de humedecimiento y secado (AWD) (2600 m2), y (iii) riego por filtraciones laterales (FL) (3900 m2), ubicado en el sector de riego Overazal, distrito de Chongoyape, región Lambayeque, Perú, con una altitud de 211 m.s.n.m. y un clima árido y templado (Figura 1).
Se realizó en un periodo de seis meses, entre agosto del 2021 a enero del 2022, presentando temperatura mínima y máxima de 19 ºC y 31 ºC, respectivamente, humedad relativa máxima de 94% y precipitación de 16,44 mm, obtenidos de la estación meteorológica Tinajones del SENAMHI.
Figura 1 Ubicación del área de estudio. (a) Distribución de parcelas según sistema de riego, (b) Ubicación del área de estudio en la costa norte del Perú, (c) Imagen satelital del distrito de Chongoyape. El sistema de referencia es EPSG:32171. La imagen se obtuvo a partir de vuelos de VANT.
Según Köppen registra un clima seco árido cálido. El suelo es franco arcilloso (28% arcilla, 48% arena, 24% de limo), conductividad eléctrica (CE) de 0,406 dS/m, potencial de hidrógeno (pH) de 7,51, capacidad de intercambio catiónico (CIC) de 12 meq/100 g, rela ción carbono/ nitrógeno (C/N) < 0,08, materia or gánica de 2,26%, nitrógeno total de 1287 mg kg-1, fosforo disponible de 20,7 mg∙kg-1, boro < 0,50 mg∙kg-1, densidad aparente de 1,65 g∙cm-3, capacidad de campo de 0,35 cm3∙cm-3 y porosidad del 45,28%.
2.2 Sistema de riego
Se tuvieron tres sistemas de riego (IC, AWD, y FL). El AWD se realizó entre los estadios de macollamiento y botoneo (Figura 2). En IC se realizaron 15 riegos distribuidos: uno en almacigo, tres en fase vegetativa, siete en fase reproductiva y cuatro en fase de maduración, con una duración de 4 a 5 horas por riego; para AWD se aplicó nueve riegos, reduciéndose a: uno en fase vegetativa y cuatro en la fase reproductiva; y en FL no se realizó riego directamente, sino recibió filtración lateral (FL) de la poza 4 adyacentes (Rahman et al., 2022).
Se identificó la fenología de acuerdo a los grados días de calor (GDD) en base a las ecuaciones (1) y (2), índice útil para la identificación del crecimiento del cultivo (Porras et al., 2020):
()1
()2
Donde, T a , es la temperatura media del aire; T c-min , es el umbral mínimo de temperatura del cultivo, a partir del cual no se produce desarrollo vegetal.
La fuente de agua de riego proviene del canal que abastece al reservorio Tinajones; y tiene las siguientes características químicas: CE de 221 µS cm-1, pH de 7,6, cationes (Ca2+ de 1,96; Mg2+ < 0,31; Na+ de 0,18; K+ de 0,0,07) meq l-1, aniones (Cl-1 < 0,3; NO3 2- < 0,16; SO4 2- de 0,48) meq l-1. Según U.S. Soild Salinity Laboratory la clasificación del agua de riego es de tipo C1-S1, debido al bajo contenido de sodio y escasa salinidad del agua, con una relación de absorción de sodio (RAS) de 0,169.
El volumen de agua aplicado se estimó con el caudal y tiempo de riego, medido en el canal de distribución, con el método de la velocidad media (flotador) y área transversal al flujo. Para analizar los efectos del sistema del riego, se estimó la variación en: el consumo de agua (ΔCA) (Ecuación 3), el rendimiento (ΔY) (Ecuación 4), y la eficiencia de uso de agua (ΔEUA) (Ecuación 5), entre AWD y IC (Cheng et al., 2022).
()3
()4
()5
Donde, las unidades de medida para CA son m3∙ha- 1, para Y son kg∙ha- 1 y para EUA son kg∙m- 3.
2.3 Estimación de la evapotranspiración
La Figura 3 representa el procedimiento secuencial para la estimación de la evapotranspiración del cultivo (ETc). Se emplearon dos VANTs (1) para colectar imágenes multiespectrales (2), RGB (3) y térmicas (4), a partir del cual se estimaron los albedos (5), e índices de vegetación como, el Índice de Vegetación de Diferencia Normalizada (NDVI) (6,7), el Índice de Área Foliar (IAF) (8), previamente ajustados con información de campo.
El modelo digital de terreno (MDT) (9) se estimó del ortomosaico RGB. Se colectó información climática (10) como humedad relativa (HR) (11) y otras variables climáticas (16) como temperatura del aire (Ta), velocidad del viento (V), altura de césped (h), evapotranspiración de referencia (ETo). Se estimó la temperatura superficial (12), ajustada con información de campo (13). Con la información del MDT se calculó la presión atmosférica (P)(17), y luego la densidad del aire (ρair) (18). La información climática permitió calcular la resistencia aerodinámica (rah) (19); y en el caso del flujo de calor sensible del aire (H)(22), se calculó el diferencial de temperatura (dT) (21) a partir de la identificación de valores umbrales llamados pixel frío y caliente (20). Con los datos previos, se estimaron los componentes del balance de energía como la radiación neta (Rn) (14), el flujo de calor del suelo (G)(15), (H)(22) y el flujo de calor latente (LE) (23). En última instancia, se estimó la fracción de evapotranspiración de referencia (ETrF) (24) y la ETc (25).
2.3.1. Vehículo aéreo no tripulado (VANT)
Se realizaron los monitoreos aéreos con dos cuadricópteros (i) Matrice 210 con cámara multiespectral Parrot Sequoia (Parrot S.A., París, France) y (ii) Matrice 300 RTK (DJI, Shenzhen, China) con cámara térmica DJI Zenmuse H20T (DJI, Shenzhen, China); con cinco vuelos por VANT, haciendo un total de 10 vuelos. La cámara multiespectral capturó 1 000 imágenes aproximadamente entre todos los sensores, con calibración sobre un panel de reflectancia llamada panel de reflectancia (AIRNOV). La cámara térmica capturó 694 imágenes térmicas y ópticas (RGB, TIF radiométrico). Mayores detalles de las características de los sensores se indican en la Tabla 1, y características del plan de vuelo en la Tabla 2.
Figura 3 Procedimiento secuencial para la estimación es la evapotranspiración. Siendo temperatura del aire (Ta), velocidad del viento (V), altura de césped (h), humedad relativa (HR), evapotranspiración de referencia (ETo).
Tabla 1 Característica de los sensores de las cámaras VANT
| VANT | Cámara | v (m s-1) | Traslape frontal y lateral | Tiempo de captura (s) | Tiempo de vuelo (min) | Tamaño del píxel (cm/píxel) |
| Matrice 210 | Parrot Sequia | 7,4 | 85/80 (multiespectral) | 1,6 | 10 | 12,4 |
| Matrice 300 - RTK | Zenmuse H20T | 3,4 | 75/75 (térmica) 90/60 (RGB) | 1 | 30 | 8,9 3,2 |
Tabla 2 Características del plan de vuelo. Plan de vuelo (altura de 100m, hora de vuelo: entre las 11:00 am y 13:00 pm*)
* Para minimizar el sombreado, bajo condiciones de cielo despejado y con un viento ligero.
Los valores de temperatura de las imágenes térmicas se calibraron con información de un radiómetro Apogee MI-210 (Apogee Instruments, MI-210, Utah, Estados unidos) sobre nueve coberturas conocidas (aluminio, vegetación muerta, vegetación viva, poliestireno expandido, tela amarrilla, tela negra, tela roja, tela verde y suelo desnudo) enmarcadas en un cuadrado de PVC de 1 m x 1 m (largo y ancho). La medición se realizó a una altura de 40 cm, colectándose un total de 81 lecturas por día de vuelo, haciendo un total de 405 lecturas.
En la calibración de las imágenes térmicas, se obtuvo una ecuación de corrección por día de vuelo (Figura 4), además se generó una ecuación general (Y radiómetro = 0,8061∙X H20T + 14,667) cuyo coeficiente de determinación (R2), coeficiente de Pearson (r) y error cuadrático medio (RMSE) fueron de 0,953, 0,976 y 12,77, respectivamente. Estos re sultados indican una buena correlación entre las temperaturas de la cámara térmica H20T y el radiómetro Apogee MI-210, mejor al obtenido por Torres-Rua (2017), quien relacionó imágenes térmicas colectadas con un VANT y datos registrados en campo con un radiómetro.
Se georreferenciaron ocho puntos de control en suelo, empleando el equipo de prospección GNSS del D-RTK 2, con una precisión de ±1 cm + 1 ppm y ±2 cm + 1 ppm, en horizontal y vertical, respectivamente (ppm = por cada 1 km de aumento en la distancia, la precisión será 1 mm menos).
Figura 4 Relación entre las temperaturas registradas por el radiómetro Apogee MI-210 y cámara térmica H20T.
2.3.2. Medición en campo de NDVI y IAF
Se instalaron nueve puntos de muestreo distribuidos en seis pozas, identificadas con varillas de acero de 1,5 m y en la parte superior, un cilindro de poliestireno expandido de 0,50 m de diámetro y 0,10 m de altura cubiertas de papel aluminio. En estos nueve puntos se realizaron mediciones específicas sobre el NDVI, IAF y rendimiento del cultivo.
La medición del NDVI se realizó con el reflectómetro portátil GreenSeaker™, marca Trimble, colectando la lectura a una altura de 60 cm sobre las plantas de arroz. La medición del IAF se realizó con un ceptómetro portátil ACCUPAR LP-80 (METER Group, Pullman, WA, USA) que consta de una matriz lineal de 80 sensores independientes de radiación fotosintéticamente activa (PAR) y un sensor externo, con lecturas en el rango de longitud de onda de 400 a 700 nm en unidades de μmol m-2 s-1 y precisión ± 5 %; mide la cantidad de luz interceptada en el dosel del cultivo. El sensor PAR se mantuvo debajo de las plantas de arroz y el sensor externo a una altura de 60 cm sobre las plantas de arroz. Asimismo, se registró cuatro mediciones por punto de muestreo, con un total de 36 datos por fecha de registro, haciendo un total de 180 datos para el sensor NDVI y IAF.
Se obtuvo el ajuste lineal del NDVI medido en campo y estimado con imágenes multiespectrales a todas las fechas de monitoreo, generándose ecuaciones de correlación para cada día de vuelo para la relación NDVI (Y NDVI green seeker = m∙X NDVI CAMARA H20T + n) (Figura 5a). En la relación IAF y NDVI, se obtuvo una relación lineal para cada día de vuelo (X IAF ceptómetro = m∙Y NDVI corregido + n) (Figura 5b), relaciones específicas para las condi ciones climáticas locales y del dosel del cultivo que se encontraba entre macollamiento y estado de huso.
Al respecto, Sisheber et al. (2022) obtienen una relación lineal del IAF (ceptómetro) con NDVI estimados de imágenes LANDSAT (R2 = 0,45). Por otro lado, Mendoza et al. (2017) indican que las mediciones del IAF con ceptómetro en el cultivo de ají (Capsicum annuum) en invernadero y IAF obtenido por método destructivo, pueden llegar a tener un R2 de 0,99.
2.3.3. Preprocesamiento de las imágenes VANT
El preprocesamiento de las imágenes se realizó con el software fotogramétrico Pix4Dmapper Pro (Pix4D S.A., Prilly, Switzerland), principalmente tres procedimientos estándares (alineación de fotos, generación de nube de puntos y generación de ortomosaicos con modelos digitales de superficie (MDS), la secuencia de levantamiento y generación de mosaicos de las imágenes de un VANT) (Río et al., 2019). El software presenta plantillas para generar mosaicos con imágenes RGB, imágenes multiespectrales, imágenes térmicas y modelos en 3D; por ello el programa tiene un registro de modelo de cámaras, de cuales reconoce los formatos de las imágenes capturadas permitiéndole procesarlas.
Se generaron ortomosaicos de las imágenes RGB, multiespectrales y térmicas, respectivamente; luego se homogenizo la resolución a 14 cm/píxel por el método resample, y se realizó una clasificación supervisada (máxima verosimilitud) de las imágenes mediante el software ArcGIs 10,4 para obtener la máscara de vegetación por cada ortomosaico.
Figura 5 Relación entre NDVI estimados de imágenes multiespectrales y medido con GreenSeaker en campo (a), relación entre NDVI corregido y LAI con ceptómetro ambas medidas en campo (b).
Se calculó el NDVI empleando las bandas espectrales RED y NIR y se comparó con el NDVI medido en campo con el GreenSeakerTM, obteniéndose una ecuación lineal que permitió estimar los valores de NDVI calculados. La fórmula empírica para hallar el NDVI se presenta a continuación en la ecuación (6):
()6
El IAF, que representa la cantidad de radiación que las plantas utilizan en la fotosíntesis, se obtuvo de la relación entre el IAF medido en campo por el ceptómetro y NDVI corregido con datos recolectados en campo, acorde a lo realizado por Fan et al. (2009).
2.3.4. Modelo de balance de energía
Según Allen et al. (2007), Mapping EvapoTranspiration at high Resolution with Internalized Calibration (METRIC), es originalmente un modelo basado en imágenes de satélite que estima la ET mediante la aplicación de balance de energía superficial. La ET se calcula como un residuo de la ecuación de balance de energía superficial mediante el flujo de calor latente (LE). La fórmula empleada para el cálculo del LE es presentada en la ecuación (7).
()7
Donde, R n es la radiación neta (W∙m-2); G es el flujo de calor del suelo (W∙m−2); y H es el flujo de calor sensible (W∙m−2) (Allen et al., 2011).
El cálculo de R n requiere de manera previa la obtención de imágenes de NDVI, e imágenes IAF, ajustadas a datos de campo; además el MDT calculado se aplica en la radiación de onda corta incidente y radiación de onda larga incidente, la T s y el IAF en el cálculo de la radiación de onda larga saliente; esta secuencia es detallada por Bastiaanssen et al. (1998).
El G representa el calor conducido al suelo, su cálculo depende del valor del IAF obtenido, de la T s obtenidos de las imágenes capturadas por el VANT y la R n , y se calcula mediante las siguientes ecuaciones (8a y 8b). (Allen et al., 2011).
() ()8
Donde T S es la temperatura de la superficie (K) y IAF es adimensional (Allen et al., 2011).
El flujo de calor sensible (H) es la tasa de pérdida de calor por el aire a través de la convección y conducción, se calcula a partir de la rugosidad de la superficie, velocidad del viento y la diferencia de temperatura de la superficie y del aire, extrayendo datos de la estación meteorológica se calcula la densidad del aire (ρ air , kg∙m-3), Cp es el calor específico del aire a presión constante (J∙kg-1∙K-1), el dT (K) es la diferencia de la temperatura entre las alturas que definen la capa de intercambio calórico entre el aire y la superficie y r ah es la resistencia aerodinámica (s∙m-1) correspondiente a dT (Allen et al., 2011; Tasumi, 2019). El flujo de calor sensible (H) es la transferen cia de calor y vapor que se da entre la cobertura vegetal y la atmósfera ocurre entre Z1(superficie evaporante) y Z2 (nivel de referencia) (Allen et al., 2011).
Se calcula el gradiente de temperatura (dT) para cada pixel basado en la relación lineal entre el dT y la temperatura de la superficie (Ts) para los píxeles de anclaje frío y caliente (Allen et al., 2011), aplicando cada resultado a la ecuación (9) (Acharya & Sharma, 2021).
()9
En la selección del píxel frío y caliente se definieron condiciones límites para el equilibrio energético. Los píxeles están en función de la homogeneidad entre sus píxeles vecinos y la distancia a la estación meteorológica. El píxel frío se identificó en campo agrícola, con selección de un subconjunto de píxeles con valores de temperatura en el rango ± 0,2 K del promedio del 20% más frío, extraído del área del 5% de los valores más altos del NDVI. Finalmente, se seleccionó el 2% del 5% de los valores más altos del NDVI, siendo el píxel candidato final el promedio de este 2%. Por lo contrario, el píxel caliente se identificó en suelos desnudos sin vegetación, con selección de subconjunto de píxeles con valores de temperatura en el rango del ± 0,2 K del promedio del 20% más caliente, extraído del área del 10% de los valores más bajos de NDVI. La selección de los píxeles finales se da en función de la proximidad del valor promedio y la homogeneidad entre sus píxeles vecinos (Morton et al., 2013; Bhattarai et al., 2017).
3. Resultados y discusión
3.1 Rendimiento y consumo de agua
Se identificaron las fases fenológicas del cultivo, a 66, 98 y 130 DDS para vegetativa, reproductiva y de maduración, destacando los estadios de macollamiento, estado de huso y floración, a 66, 76 y 89 DDS respectivamente (Porras et al. 2020).
En la Tabla 3, se presenta el consumo de agua y número de riegos aplicados por fase del cultivo de acuerdo a la condición de riego (IC y AWD) en sus pozas correspondientes, observando un ahorro de agua del 41,90% con AWD respecto a IC; y una reducción del 13,66% del rendimiento, pero un incremento de la EUA en 48,94%. Al respecto, los cortes de agua en la etapa reproductiva (floración) con AWD han impactado negativamente en el rendimiento, ya que debieron estar entre macollamiento y estado de huso (Ishfaq et al., 2020). La implementación de AWD en la etapa reproductiva causó una reducción del rendimiento en comparación con la implementación de AWD en la fase vegetativa. Además, la textura franco-arenosa pudo influir en un mayor estrés hídrico, por el poco almacenamiento del agua; por lo tanto, se debe explorar la gestión óptima del AWD sin penalizar el rendimiento en condiciones del suelo actual.
3.2 Selección de píxel frio y caliente
La selección de los píxeles extremos caliente y frío para las fechas de campo (Tabla 4, muestra que los píxeles fríos corresponden a una temperatura superficial de (T s ) mínima de 283,22 K, y valores máximos de NDVI, IAF y albedo de 0,83, 3,03 y 0,33, respectivamente. Los píxeles calientes corresponden a una temperatura superficial (T s ) máxima de 331,23 K, valores mínimos de NDVI, IAF y albedo de 0,06, 0,005 y 0,21, respectivamente. Las características anteriormente señaladas corroboran lo expuesto por Allen et al. (2011) y Morton et al. (2013), quienes afirman que el proceso de verificación de los píxeles frio y caliente se realizan en áreas separadas, es decir que el píxel frío se debe seleccionar en coberturas vegetales, mientras que el píxel caliente se debe seleccionar en suelos desnudos sin vegetación.
3.3 Balance de energía
En la Figura 6 se presenta la variación espacial de los componentes de balance de energía Rn, G, H y LE en macollamiento (23 y 28 de octubre), en primordio floral (14 y 17 de noviembre) y estado de huso (20 de noviembre) del cultivo de arroz. En la Tabla 4 se muestra los resultados del NDVI, IAF y componentes del balance de energía por sistema de riego. Para el riego con IC los valores varían en el rango de 540,63 a 621,09 W∙m-2, de 47,99 a 56,69 W∙m-2, y 159,68 a 206,31 W∙m-2 para Rn, G y H, respectivamente. Para AWD, los valores varían de 553,76 a 655,75 W∙m-2, 52,70 a 58,11 W∙m-2, y 139,36 a 206,39 W∙m-2 para Rn, G y H, respectivamente. Para FL, los valores varían entre 527,88 a 591,21 W∙m-2, 46,71 a 52,51 W∙m-2, y 160,00 a 208,76 W∙m-2 para Rn, G y H, respectivamente.
Según Allen et al. (2011) y Nassar et al. (2021), la Rn puede encontrarse entre 100 a 800 W∙m-2 dependiendo de la variación diurna, alcanzando un mayor valor al medio día. Otros autores indican que G varía de 0 a 180 W∙m-2 (Nassar et al., 2021), y H varía de 50 a 400 W∙m-2 (Acharya & Sharma, 2021). Los resultados obtenidos en este estudio están dentro de estos rangos. Sin embargo, Nassar et al. (2021) indican que los componentes del balance de energía alcanzan valores más altos al medio día, con variación horaria según estación del año (verano, primavera, otoño e invierno) y valores más altos en verano.
Figura 6 Variación espacial de Radiación Neta (Rn), Flujo de calor del suelo (G), flujo de calor del sensible (H) y flujo de calor latente (LE) en macollamiento (49 DDS, 54 DDS), en primordio floral (71 DDS, 74 DDS) y estado de huso (77 DDS) del cultivo de arroz.
Tabla 4 Selección del píxel frío y caliente con valores de NDVI, albedo, temperatura y las constantes de calibración [dT = a + b∙T S ; (Allen et al., 2011)]
Tabla 5 Resumen de los valores del índice de vegetación bajo riego de inundación continua (IC), alternancia de humedecimiento y secado (AWD), y filtración lateral (FL)
La variación espacial del H (Figura 6) está directamente influenciada por la temperatura superficial. La variaron térmica del día 54 DDS, ocurre por el apagado y encendido del VANT para realizar el cambio de batería. La cámara H20T es una cámara no refrigerada, por lo que el propio calor interno del sensor influye en la captura de datos de temperatura. Esto podría generar una variación en las lecturas, por lo que se sugiere no apagar el equipo para cambio de batería.
3.4 Evapotranspiración del cultivo diario (ETc)
En la Figura 7 se aprecia la variación espacial de la ETc diaria, con valores (x, S) para las fases de macollamiento, primordio floral y estado de huso, según condición de riego.
En la Tabla 5 se presentan los valores de la ETc, según estadio y sistema de riego, obteniéndose valores (X±S) de 7,28, 8,21 y 6,29 mm∙d-1 (macollamiento), 5,62, 6,36 y 5,33 mm∙d-1 (primordio floral) y 5,68, 5,85 y 5,41 mm∙d-1 (estado de huso), con valores máximos de 13,72, 7,60 y 7,35 mm∙d-1 para IC, y 13,17, 7,59 y 7,35 mm∙d-1 para AWD.
En la Tabla 6 se muestra los promedios y desviación estándar de los componentes de balance de ener gía, bajo diferentes sistemas de riego. Se observa que en un riego AWD hay una reducción de 6,56% y 32,43% para NDVI e IAF, 21,28% para al bedo y 11,58%, 8,09% y 2,45% para Rn, G y H, res pectivamente, respecto a FL. Se obtienen valores promedios de ETc de 7,07, 6,79 y cd 6,54 mm∙d-1 para IC, FL y AWD, respectivamente. Además, se obtiene con el riego AWD un ahorra en el consumo de agua de 42% (ΔCA), pero una reducción del rendimiento del 14% (ΔY); esto se traduce en un incremento del 49% (ΔEUA) en la eficiencia del uso de agua.
Figura 7 Variación espacial de la evapotranspiración (ET) en macollamiento (49 DDS, 54 DDS), primordio floral (71 DDS, 74 DDS) y estado de huso (77 DDS) del cultivo de arroz.
Tabla 6 Resultados finales de los componentes de balance de energía, Evapotranspiración del cultivo (ETc) y otros resultados, según el riego
| Detalle | Inundación continua (IC) | Filtración Lateral (FL) | Alternancia de humedecimiento y secado (AWD) |
|---|---|---|---|
| Albedo | 0,18 ± 0,01 | 0,18 ± 0,01 | 0,20 ± 0,01 |
| Térmica (ºC) | 22,19 ± 0,68 | 22,71 ± 0,89 | 22,34 ± 0,78 |
| NDVI | 0,78 ± 0,02 | 0,80 ± 0,01 | 0,80 ± 0,02 |
| IAF (m2∙m-2) | 1,8 ± 0,27 | 2,07 ± 0,17 | 2,07 ± 0,30 |
| Rn (W∙m-2) | 622,17 ± 13,65 | 618,12 ± 12,34 | 604,06 ± 14,86 |
| G (W∙m-2) | 44,14 ± 2,85 | 43,11 ± 1,84 | 43,00 ± 2,76 |
| H (W∙m-2) | 198,14 ± 10,62 | 206,45 ± 12,64 | 199,23 ± 13,40 |
| LE (W∙m-2) | 402,25 ± 21,87 | 376,93 ± 23,67 | 367,03 ± 26,79 |
| ETc (mm∙d-1) | 7,07 ± 0,43 | 6,79 ± 0,46 | 6,54 ± 0,53 |
| Rendimiento (t∙ha-1) | 10,54 ± 1,24 | 10,31 ± 0,35 | 9,10 ± 0,98 |
| Consumo de agua (m-3) | 22287,23 | - | 12950,1 |
| Eficiencia en el consumo de agua (EUA) (kg∙m-3) | 0,47 | - | 0,7 |
| Variación del consumo del agua (∆IW) | -41,90% | ||
| Variación del rendimiento (∆Y) | -13,66% | ||
| Variación de la eficiencia en el consumo del agua (∆EUA) | -48,94% |
Los valores de IAF en los arrozales con AWD, son los menores en comparación a los otros métodos de riego, esto por el estrés hídrico ocasionado por riego de alternancia, humedecimiento y secado en la etapa de macollamiento. La Rn están dentro del rango obtenido por Lee & Kim, (2016), quien monitoreo cultivo de arroz con información de imágenes del satélite MODIS, obteniendo valores de Rn de 400 a 800 W∙m-2.; con valores de G, de 28 a 78 W∙m-2 y H hasta 241 W∙m-2; siendo valores cercanos a los obtenidos por Lee & Kim (2016).
Los valores de NDVI registrados están en aumento, según aumente los días del cultivo, esto como respuesta al desarrollo del cultivo y fertilización. Sin embargo, los valores promedios son menores a 0,8 lo que podría deberse a una dosis menor de fertilización por los altos costos en el mercado, además puede haberse afectado por una geminación no uniforme y leve disminución de la densidad poblacional (Gonzalez & Mayorga, 2018). Sin embargo, estos valores son mayores a los obtenidos por Garcia & Martinez, (2009) que registraron valores de 0,36 a 0,71 en arrozales en Colombia, tomadas con un VANT.
En la Tabla 7, se presentan valores de ETc estimados para diferentes condiciones climáticas, sistema de riego y método de ETc; con valores máximos de ETc se encuentran en el clima tropical, templado mediterráneo y árido (Senegal, Estados Unidos y Irán) con 10,50, 8,40 y 8,40 mm∙d-1. En otros estudios presenta valores superiores a lo estimado en la zona de estudio, además presentan mayor afinidad en las condiciones climáticas. Sin embargo, en las condiciones climáticas de templado (Lee & Kim, 2016), húmedo tropical (Sawadogo et al., 2020), monzónico tropical (Saha et al., 2022) y monzónico subtropical (Liu et al., 2018) y continental cálido (Wei et al., 2023), presentan valores inferiores de 5,30, 6,17, 5,64, 4,03 y 5,10 mm∙d-1. Aunque, en el clima subtropical húmedo de China (Wei et al., 2023), se asemeja al resultado obtenido de 7,07 mm d-1. Además, de las similitudes en algunas condiciones climáticas, también presenta diferentes ETc debido a las variedades de arroz, morfología de suelo y a las condiciones antrópicas.
La fase de cultivo que registro mayor promedio de ETc fue en la fase de macollamiento, por presentar riego unos días antes del vuelo del VANT, como fue el 28 de octubre (54 DDS), presentando esta fecha el mayor ETc, además ese día, se registró la temperatura máxima de los días de vuelo, condición que favoreció un incremento de la ETc.
Tabla 7 Valores de evapotranspiración del cultivo (ETc) de diversos estudios
* Alternancia de humedecimiento y secado (AWD).
Además, es posible que la temperatura superficial captada por la cámara térmica esté alterada por el cambio de batería realizada durante el monitoreo del VANT.
La ETc en el arrozal, fue menor en las fases de pri mordio floral y estado de huso, debido a un des perfecto con la bomba de riego, por lo cual se tuvo cinco días sin riego. La reducción del rendimiento puede ser causada por la reducción de la dosis de fertilización empleada y control parcial del ataque de plagas (Pyricularia oryzae y Tagosodes orizicolus), debido al alto costo en el mercado.
4. Conclusiones
Este estudio evalúa los mosaicos multiespectrales y térmicos de alta resolución espacial mediante la calibración de los parámetros del modelo de balance de energía para la estimación del modelo METRIC en arrozales. El ajuste de las imágenes térmicas es importante para la implementación del modelo METRIC, por lo tanto, el uso de un radiómetro con lecturas en diferentes coberturas obtuvo un RMSE de 12,5 ºC, mostrando resultados confiables para la estimación de los componentes del modelo METRIC.
La calibración de los índices de vegetación NDVI e IAF con mediciones de campo, resultaron en un RMSE de 0,099 y 1,219, respectivamente. Los componentes de balance de energía Rn, G, H y LE, mostraron valores aceptables dentro de los rangos establecidos, para el área de estudio en las diferentes fases de desarrollo del cultivo.
Se obtienen valores medios de 6,34 y 5,84 mm∙d-1 para riego con inundación continua y alternancia de humedecimiento y secado, respectivamente. Además, se obtiene un ahorro en el consumo de agua del 42%, pero una reducción del rendimiento del 14%, que se traduce en un incremento del 49% en la eficiencia del uso de agua, resultados iniciales que permitirán optimizar el uso del agua y mejorar la gestión del riego. Para futuras investigaciones se sugiere realizar un estudio con mayor margen de área, con diferentes variedades que sean resistente al estrés hídrico y la utilización en conjunta de imágenes satelitales de Sentinel y VANT, y obtener un plan de gestión hídrica para afrontar el cambio climático.
