INTRODUCCIÓN
La práctica de riego es uno de los procesos más complejos que debe realizar el agricultor debido a la multitud de factores implicados en el manejo del agua: se requiere de información técnica para el conocimiento preciso entre el agua usada y el rendimiento de los cultivos (García et al., 2009; Domínguez et al., 2012).
En la investigación y planificación agrícola, el desarrollo de modelos para simular diferentes procesos relacionados con su eficiencia, se ha convertido en una práctica común que, sustentada con la información científica disponible, es útil para pronosticar resultados en situaciones y condiciones específicas, lo que permite plantear nuevas hipótesis y orientar la investigación o el manejo hacia los puntos más críticos (Gormley y Sinclair, 2003).
Entre los modelos existentes se encuentra el AquaCrop-FAORaes et al.(2009); Steduto et al.(2009), modelo de aplicación general de cultivo, el cual se ha utilizado en varias zonas del mundo bajo diferentes condiciones ambientales. El mismo es desarrollado por la FAO con la ventaja de ser relativamente sencillo en comparación con otros modelos de simulación de rendimientos, por lo que con solo calibrar algunos parámetros se obtienen resultados similares a la realidad. Esto facilita el uso del modelo en tanto la información necesaria para adaptarlo a las condiciones de nuestro país no sea tan grande en comparación a otros modelos y la plataforma del software resulte ser de fácil manejo.
Recientemente, González et al. (2018) y Gonzaléz et al. (2019) calibraron y validaron el modelo AquaCrop para maíz, sorgo y soya en suelo Ferralítico Rojo compactado en Cuba (municipio Alquízar, Provincia Artemisa). Los resultados de la calibración del modelo permitieron optimizar los parámetros fundamentales de suelo y cultivo para su aplicación en las condiciones de estudio, destacándose la duración de las fases fenológicas, el índice de cosecha, los límites de agotamiento de la reserva de agua para definir el estrés hídrico y la curva número para la determinación del escurrimiento superficial.
La validación de un modelo es la aprobación, a través de procedimientos estadísticos adecuados, del ajuste entre lo simulado y observado de tal modo que las interpretaciones del sistema real, sean correctas. Sin embargo, la sensibilidad de los modelos al rango de condiciones agroecológicas en que se desarrollan los sistemas a evaluar es un requisito indispensable para definir su utilidad e implica, de manera estricta, el estudio de la variación de sus resultados en función de diferenciaciones graduales en el valor de sus parámetros (Tittonell, 2007, citado por Giayetto y Morla (2013).
Varios autores han realizado análisis de sensibilidad al modelo AquacropcomoRaes et al. (2006); Cortés et al. (2013); Giayetto y Morla (2013), variando datos de entrada de parámetros climáticos, parámetros del cultivo, profundidad de raíces y el factor de agotamiento para el estrés hídrico y del suelo; encontrando por ejemplo, que alteraciones del coeficiente del cultivo pueden provocar sobre y subestimación de la transpiración, el estrés hídrico y la producción del cultivo.
El frijol común (Phaseolus vulgaris L.) forma parte del grupo de las leguminosas comestibles. Este grano se considera estratégico, no solo por sus propiedades nutricionales y culinarias, sino además por su presencia en los cinco continentes del mundo y su importancia para el desarrollo rural y social de muchas economías (Thornto, 2014). La mayor contribución del frijol común en Cuba está asociada a la sustitución de importaciones y la seguridad alimenta, los granos contienen una abundante gama de vitaminas, fibra vegetal y minerales, entre los que se destaca el hierro. Sin embargo, su mayor valor nutricional radica, básicamente, en un alto contenido proteico que oscila entre el 12% y el 25% del peso de las semillas secas, es decir, 2,5 veces mayor al de los cereales (IIG-MINAG, 2013).
El objetivo de este trabajo fue calibrar el modelo AquaCrop-FAO para el cultivo de frijol en suelo Ferralítico Rojo compactado para las condiciones climáticas de la región de Alquizar, Artemisa y realizar un análisis de sensibilidad del modelo a los cambios medioambientales.
MÉTODOS
Localización del área de estudio
El modelo se calibró con datos provenientes de experimentos de campo realizados durante la época de invierno (noviembre-marzo), en la Estación Experimental del Instituto de Investigaciones de Riego y Drenaje, en el municipio de Alquízar, coordenadas: Latitud 22º 46' N y Longitud 82º 36' W, altura sobre el nivel medio del mar 6 m, a 12 km de la costa. Estos experimentos se llevaron a cabo con el objetivo de estudiar el efecto de la escasez de agua sobre el desarrollo del cultivo y el rendimiento del frijol negro de la variedad Cueto 25-9 (Giralt, 1979). La fecha de siembra fue el 17 de noviembre, la madurez se alcanzó a los 109 días después de la siembra, la densidad de población fue 222 000 plantas ha-1, el rendimiento máximo fue de 2,63 t ha-1 y el método de riego fue aspersión semiestacionario de media intensidad.
Según Giralt (1979), la temperatura media mínima y máxima variaron entre los 20 y 31°C respectivamente, la lluvia total para los meses en estudio fue de 167,7 mm. El mayor valor de precipitación caída correspondió al mes de diciembre con un total de 64,1 mm, mientras que el mes de menor valor fue noviembre con 12,8 mm mensual. La evapotranspiración de referencia (ETo) tuvo un valor total para los meses en estudio de 371,9 mm. Los valores máximos de evapotranspiración diaria se observaron en las fases de floración-fructificación (3,34-3,98 mm). Los promedios mensuales máximos y mínimos corresponden a el mes de enero (132,9 mm) y noviembre (30,7 mm), respectivamente.
En la Figura 1 se muestra el balance hídrico decenal de precipitación y ETo para la zona y meses en estudio, donde se refleja la variación de ambos parámetros. Este balance hídrico es negativo para 7 decenas (58%) de un total de 12 existentes en el periodo de estudio, por lo que el riego es imprescindible para la obtención de buenos rendimientos agrícolas.
FIGURA 1.
Balance hídrico decenal para los meses de noviembre-marzo en la Estación Experimental del Instituto de Investigaciones de Riego y Drenaje, en Alquízar.
El tipo de suelo en estudio es Ferralítico Rojo compactado según la segunda clasificación genética de los suelos de Cuba por Instituto de Suelos (1979); Hernández et al. (1999), de la zona sur de La Habana. Esta zona forma parte de la llanura cársica Habana-Matanzas, muy plana, con una altura inferior a 30 metros del nivel del mar y donde la agricultura es la actividad económica fundamental. Las propiedades físicas fundamentales para cada capa del perfil del suelo Ferralítico Rojo compactado han sido reportadas y actualizadas por Cid (2006) y Cid et al.(2011) y se resumen en la Tabla 1 y conforman el fichero de suelo que se utilizó para todas las corridas del modelo.
TABLA 1.
Análisis granulométrico y otras propiedades y parámetros físicos e hidrofísicos del perfil del suelo Ferralítico Rojo compactado
Nota: da a θ FC densidad volumétrica o aparente del suelo determinada; θ FC contenido de agua en el suelo correspondiente a capacidad de campo; θ wp contenido de agua en el suelo correspondiente al punto de marchitez permanente; θ sat contenido de agua en el suelo a la saturación; K sat conductividad hidráulica del suelo saturado; Pt porosidad total del suelo; Vi nf.bas velocidad de infiltración básica en el suelo.
Según estos autores de acuerdo a su funcionamiento hídrico, el suelo presenta un comportamiento diferenciado en la capa de 0-20 cm con respecto a los valores de densidad aparente, el proceso de retención de humedad y de la conductividad hidráulica. A partir de los 40 cm hay una capa más compactada que diferencia estos procesos hidrodinámicos.
Este suelo pertenece al grupo A de los suelos cubanos por clases hidrológicas (cultivo en hileras, surcos rectos y condiciones hidrológicas malas, debido a la compactación del suelo) según López et al. (1998), por lo que se adoptó la curva número (CN) 71 para representar el escurrimiento superficial y se consideró que no existe estancamiento de agua en superficie. Para la calibración del modelo se optó por no evaluar de manera particular ninguna reducción del crecimiento asociada al nivel de fertilización, evitando la mayor complejidad que supone analizar la interacción entre la disponibilidad de agua y la nutrición de los cultivos.
Modelo AquaCrop
De acuerdo a lo reportado por Raes et al. (2009), el modelo AquaCrop consta de varias ecuaciones que, con datos de clima, densidad de población, características genéticas, tipo de suelo, nivel de fertilización y nivel de déficit hídrico, simulan el crecimiento y rendimiento del cultivo. En la Tabla 2 se muestran las ecuaciones que AquaCrop utiliza para el cálculo de biomasa y producción de grano.
Tabla 2.
Ecuaciones usadas por Aquacrop para la estimación de biomasa y producción de grano (tomado de Ruiz et al., 2018)
* Hi- índice de cosecha según Raes et al.(2009) relación entre el rendimiento en grano y la biomasa total producida en el cultivo del frijol.
Respetando el concepto original de un vínculo directo entre el uso del agua en el cultivo y su rendimiento, el modelo AquaCrop evolucionó a partir del enfoque del Estudio FAO: Riego y Drenaje Nº 33, al separar la evaporación del suelo no productiva (E) de la transpiración del cultivo productiva (Tr) y al estimar la producción de biomasa directamente a partir de la transpiración del cultivo real mediante un parámetro de productividad del agua. En la mayoría de los cultivos, solo una parte de la biomasa producida se destina a los órganos cosechables para producir rendimiento (Y), y la relación entre el rendimiento y la biomasa se conoce como índice de cosecha (Hi), por tanto, los procesos subyacentes que culminan en B y Hi se diferencian claramente entre sí. En consecuencia, la separación de Y en B y Hi permite considerar los efectos de las condiciones ambientales y el estrés en B y Hi por separado (Raes et al., 2009).
Una descripción minuciosa de las principales características de AquaCrop aparece en el Boletín de la FAO: Riego y Drenaje Nº 66 por Steduto et al. (2012); donde se explica detalladamente los conceptos básicos y los procedimientos fundamentales de cálculo, acompañados de descripciones adicionales relacionadas con los requerimientos de datos de entrada, la interfaz del usuario y los resultados del modelo.
Se debe especificar el contenido de agua de las capas del perfil del suelo al comienzo del período de simulación, a menos que esté en capacidad de campo. Se consideró un contenido de agua en el suelo de 0,37 cm3·cm-3, que representa el 90% del contenido de agua en el suelo a capacidad de campo (0,39 cm3·cm-3) hasta una profundidad de 0,40 m.
Análisis de sensibilidad
Para entender la sensibilidad del modelo a los cambios en las características medioambientales y cómo esa susceptibilidad puede afectar las salidas en biomasa y rendimientos, se seleccionaron algunas de las principales entradas del modelo que fueron modificadas dentro de un rango específico para observar cómo se veían alterados los rendimientos del cultivo, según procedimiento propuesto por Cortés et al.(2013). Todas las comparaciones de variación porcentuales que se realizaron fueron respecto a los valores de los rendimientos finales en la etapa de ajuste del modelo. El análisis de sensibilidad se llevó a cabo sobre las siguientes entradas:
Diferente profundidad del Nivel freático
Diferente nivel de fertilidad del suelo
Diferentes láminas rápidamente aprovechables (RAW).
Considerando diferentes texturas de suelo
Respecto a las simulaciones de los valores de rendimiento para diferentes texturas en este análisis se tomaron las características hidrofísicas que se encuentran por defecto en el módulo de suelos de AquaCrop (Tabla 3).
TABLA 3.
Características hidrofísicas de las texturas evaluadas con el modelo AquaCrop.
PMP: contenido de agua del suelo en % de volumen a punto de marchitez permanente
CC: contenido de agua del suelo en % de volumena capacidad de campo
SAT contenido de agua del suelo en % de volumen a saturación
TAW: Agua Disponible Total en el Suelo (Total Available Soil Water)
Ksat: Conductividad hidráulica del suelo saturado
Validación del modelo
Después de realizar la calibración del modelo AquaCrop con el experimento de Giralt (1979), la validación se hizo mediante análisis estadísticos de los datos de biomasa, rendimiento y cobertura del dosel, reales y los simulados con los parámetros de entrada.
La bondad del ajuste se valoró a través del coeficiente de determinación R2 y los estadígrafos: raíz cuadrada del cuadrado medio del error (RMSE), la eficiencia de modelación (EF) y el índice de Willmott (d), calculados según las expresiones propuestas por Araya et al.(2016).
De acuerdo con Heng et al. (2009), es factible el uso de estas expresiones, ya que la RMSE representa una medida global entre los valores observados y simulados, es decir, un indicador de la incertidumbre, debido a que toma las mismas unidades de la variable simulada y, por consiguiente, el valor más cercano a cero indica un desempeño bueno en la simulación. EF tiene como valor máximo y óptimo 1. Si EF es menor que 0 el modelo predice peor que utilizando el valor medio de los datos observados. El valor d varia de 0 a 1 y en ambos análisis se aplicó secuencialmente a los datos observados y simulados.
RESULTADOS Y DISCUSIÓN
Calibración del modelo
Para el ajuste del modelo AquaCrop se realizaron simulaciones, inicialmente sin modificar los parámetros conservativos del cultivo, y se compararon las salidas en cuanto a rendimiento respecto a la información de campo. De acuerdo a los resultados encontrados en el proceso de ajuste se modificaron algunos parámetros conservativos con relación a los que propone el modelo en su base de datos, como se pueden observar en la Tabla 4.
TABLA 4.
Parámetros utilizados en la calibración de AquaCrop en frijol para suelo Ferralítico Rojo compactado
La cobertura inicial y máxima del dosel vegetal fueron del 2,22 y 92%, respectivamente, mientras que el coeficiente de decline del dosel por día fue estimado como 10,9%, valores que difieren de los originales calibrados en el modelo para frijol. Aquellos parámetros relacionados con la duración de las diferentes etapas de desarrollo del cultivo sufrieron también variaciones con respecto a los usados como valores de partida para comenzar la calibración; destacándose la cantidad de días después de plantado el frijol para lograr la máxima cobertura, el comienzo de la floración y la senescencia; 65, 48 y 80 días, respectivamente.
El índice de cosecha (Hi) considerado fue de 0,33, por ser el valor más común obtenido bajo condiciones óptimas de desarrollo durante la calibración del modelo. Este valor se encuentra en el rango (0,3 a 0,5) que propone Wani et al., (2012).
El modelo tiene predeterminado el valor de la productividad del agua normalizada para frijol de 15 g m-2 (Hsiao et al., 2009; Raes et al., 2012). Los ajustes de AquaCrop normalizan la productividad del agua según la concentración de CO2 atmosférico, año tras año, en el caso particular de este trabajo este valor fue de 14 g m-2, inferior al considerado por Hsiao et al.(2009), y aspecto a tener en cuenta en otras simulaciones con este modelo.
En la Tabla 5, se muestran los resultados de la calibración del modelo AquaCrop para frijol en suelo Ferralítico Rojo compactado en Alquízar. Los valores de rendimiento y productividad del agua están en el intervalo publicado en trabajos precedentes de Gonzaléz (2013). En el caso particular del rendimiento en experimentos de campo varió entre 0,63 - 2,67 t·ha-1 y la productividad en base al rendimiento alcanzó valores de 0,97kg·m-3.
TABLA 5.
Resultados de la corrida del modelo para la calibración del frijol
Resultados muy similares a los encontrados en este trabajo fueron publicados por Ruiz et al.(2018). Estos autores, para analizar el desempeño agronómico del frijol con un ciclo de 128 días en diferentes escenarios climáticos en Colombia, se apoyaron en el modelo AquaCrop para simular productividad del agua en base a la lámina bruta e índice de cosecha, encontrando valores entre 0,87-0,98 kg m-3 y 35%, respectivamente.
En la Figura 2 se presenta la salida del modelo AquaCrop para el cultivo del frijol. En la misma, no se aprecia reducción en el desarrollo de la biomasa potencial y el suelo estuvo a capacidad de campo en casi todo el ciclo de cultivo, no presentando estrés hídrico. Solo se apreció un agotamiento de la humedad del suelo en la zona radicular (reducción promedio en la expansión foliar en un 5%) en la segunda decena de diciembre y enero, causado por la ausencia de lluvia y riego en varios días consecutivos, lo que no provocó reducción en los rendimientos finales.
Se encontró así que el valor de cobertura vegetal final (CC), determinado con datos experimentales fue de un 87%, mientras que el simulado fue de 89%, en todo el ciclo se sobrestimó la cobertura por el modelo, presentándose una diferencia de solo el 2%. Es importante que la calibración de valores de cobertura foliar sea correcto ya que para el modelo Aquacrop este parámetro afecta la tasa transpiratoria y en consecuencia la acumulación de biomasa (Farahani et al., 2009).
FIGURA 2.
Salidas del modelo para frijol.
En la Figura 3 se muestra el seguimiento de la evaporación y transpiración a lo largo del ciclo del cultivo. Se observa una transpiración máxima de 4,3 mm entre los 66-70 días de plantado el frijol, período que coincide con la floración y generalmente considerado como de mayor consumo y además, según Merino (2006) crítico para alcanzar altos rendimiento en los granos.
FIGURA 3.
Seguimiento de la evaporación y la transpiración del frijol.
Por su parte, Ergo (2013) ha señalado que no todo el ciclo del cultivo resulta igualmente crítico, sino que este periodo comienza en el inicio de la floración y finaliza 10 a 12 días después del llenado de granos, o puede prolongarse hasta el comienzo del máximo tamaño de los granos, alrededor de los 90 días de plantado, con valores de transpiración que pueden alcanzar los 3,4 mm día-1.
Análisis de sensibilidad
Comportamiento del modelo AquaCrop a diferente profundidad del nivel freático
Tomando como condición de entrada en el modelo una profundidad del nivel freático de 8 m, se puede observar en la Figura 4 que el comportamiento del rendimiento en el cultivo de frijol ante variaciones en la profundidad del nivel freático (de 1 a 2 m) no presenta afectaciones importantes para profundidades mayores a un metro; para este mismo rango de profundidad la tasa de cierre estomático no se incrementa.
FIGURA 4.Variación.
del rendimiento y cierre de estomas simulados con el modelo AquaCrop bajo diferentes profundidades del nivel freático.
Mientras que para niveles inferiores a un metro se mostrarían reducciones considerables en el rendimiento del grano, esto ocasionado, según el modelo, por un cierre de estomas dada la diminución de oxígeno en el suelo, la cual reduce la respiración de la planta y afecta el proceso de transpiración. Resultados similares fueron encontrados por Herrera et al. (1988) para frijol en Cuba y Cortés et al. (2013) en zonas representativas de producción de maíz en Colombia.
Adicionalmente, para niveles freáticos elevados, cercanos a la máxima profundidad de la raíz (0,4 m), la disminución del rendimiento se encuentra entre el 79 y 87%, presentándose un cierre de estomas de entre el 65% y 82%. Para niveles más elevados, se presentaría pérdida total del grano, ya que la planta no alcanzaría un desarrollo y producción del grano por asfixia.
Comportamiento del modelo AquaCrop a diferente nivel de fertilidad del suelo
En el modelo AquaCrop la reducción de la fertilidad está directamente relacionada con la modificación de los parámetros de desarrollo del follaje, y en consecuencia con la producción de biomasa y rendimiento. En la Figura 5, se evaluó el comportamiento del rendimiento respecto a la variación del nivel de fertilidad y se observa que a medida que disminuye este nivel, la producción lo hace gradualmente. Teniendo como base el ajuste de fertilidad de 80% para frijol en suelo ferralítico, se encontró que bajo condiciones de fertilidad de suelo de 100% se observan incrementos de 25,4% en el rendimiento. Es importante resaltar que el modelo AquaCrop asume nivel de fertilidad constante durante todo el ciclo de cultivo.
Por otro lado, con condiciones de fertilidad baja o pobre (alrededor del 40%), el rendimiento del cultivo se reduce en cerca del 49,8%, respecto al valor de referencia; y en condiciones de muy pobre puede alcanzar el 68%.
FIGURA 5.
Variación del rendimiento agronómico simulado con el modelo AquaCrop bajo diferentes niveles de fertilidad del cultivo.
Comportamiento del modelo AquaCrop a diferentes porcentajes de agua fácilmente aprovechable (RAW)
AquaCrop realiza balances hídricos internos con el fin de determinar las necesidades de agua para la planta y con base en el contenido de humedad, variables hidrofísicas del suelo, entre otras; simula el desarrollo, crecimiento y producción. En la Tabla 6, se puede observar que la disponibilidad de agua durante el ciclo del cultivo es vital para el buen desarrollo de la planta. La aplicación de riego en el cultivo permitiría aumentar los rendimientos en un 25% y reducir el porciento en la inhibición de la expansión del dosel vegetal de 39 a 0%.
TABLA 6.
Variación del rendimiento, lámina de riego y porciento en la inhibición de la expansión del dosel vegetal para diferentes porcentajes de agua fácilmente aprovechable (RAW)
| RAW | Lamina (mm) | Rendimiento | % expansión |
|---|---|---|---|
| 0% RAW | 247,9 | 2,64 | 0 |
| 20% RAW | 228,7 | 2,64 | 0 |
| 40% RAW | 216 | 2,63 | 1 |
| 50% RAW | 211 | 2,63 | 4 |
| 60% RAW | 205 | 2,56 | 7 |
| 80% RAW | 195 | 2,48 | 18 |
| 100% RAW | 185 | 1,96 | 39 |
Cuando se mantiene el límite de la humedad del suelo entre el 0% (Capacidad de Campo) y el 50% de RAW como criterio de aplicación de riego, no hay incrementos considerables en las salidas del modelo AquaCrop para el rendimiento, por lo que se estaría sobre aplicando agua (63 mm), perdiéndose esta lámina de riego en percolación profunda y escorrentía superficial. Lo anterior indica que para el cultivo del frijol no es necesario llegar a condiciones de capacidad de campo para obtener los mejores rendimientos.
De Faria et al. (1997) encontraron una disminución importante en el rendimiento cuando los niveles de humedad del suelo se redujeron a menos del 40% del agua total disponible y recomendaron que el nivel no debe descender del 60%.
En un ambiente de crecimiento húmedo Heinemann et al.(2016), compararon varias cantidades de agua (5, 6 y 7 mm·día-1) y concluyeron que los mayores rendimientos se obtuvieron cuando se aplicaron 6 mm·día-1 y la humedad disponible del suelo era el 70% en los primeros 30 cm en el momento de la aplicación. El frijol es un cultivo de raíces poco profundas extrae el 85% del agua de los primeros 40 cm del perfil del suelo (Muñoz et al., 2007).
Según Fernández (2013), el modelo AquaCrop es sumamente sensible al agua que ingresa al perfil del suelo, dado que esta se convierte en transpiración y luego por el parámetro de productividad de agua se convierte en biomasa, por lo que su variación puede disminuir considerablemente la biomasa y los rendimientos esperados.
Comportamiento del modelo AquaCrop para diferentes texturas de suelo
En la Figura 6 se evaluó el comportamiento del rendimiento respecto a la variación de las características texturales del suelo y se observa que la textura franco (F) permite un mejor desarrollo del cultivo, ya que la cantidad de agua almacenada y las tasas de infiltración permiten que las raíces absorban con mayor facilidad la cantidad de agua requerida para su desarrollo y a su vez se puedan anclar al terreno adecuadamente.
FIGURA 6.
Variación del rendimiento agronómico simulado con el modelo AquaCrop considerando diferentes texturas del suelo.
Para esta textura se observó un aumento en rendimiento de 14,5 % con respecto al suelo de referencia utilizado. Para las demás texturas las salidas simuladas fueron siempre menores a las halladas en la etapa de ajuste.
Según resultados de Cortés et al. (2013), los suelos con fracciones de arenas disminuyen la productividad agrícola, debido a su baja capacidad de retención de humedad, siendo menor la cantidad de agua disponible para ser absorbida y posteriormente transpirada por la planta.
Raes et al. (2006) consideran que en los suelos más permeables las modificaciones de la conductividad hidráulica producen cambios menores en los rendimientos, mientras que en suelos con mayor contenido de arcilla una disminución drástica de la conductividad hidráulica da lugar a alteraciones importantes de los valores simulados del escurrimiento superficial, la cantidad de agua infiltrada en el perfil del suelo, la tasa de transpiración y la correspondiente producción. Variaciones de la cantidad máxima de agua que se puede conservar en la zona de la raíz solamente produjo efectos moderados.
En el presente estudio las mayores disminuciones se presentaron para la textura arcillosa con valores mayores al 42% en rendimiento. En el caso de las texturas arcillosas el parámetro que afecta las salidas de rendimiento es la conductividad hidráulica saturada, la cual al ser baja impide el libre movimiento del agua en el perfil del suelo.
El valor de referencia utilizado fue un suelo Ferralítico Rojo compactado, 57-61% arcilla, que presentan buenas condiciones de drenaje interno pero ha perdido propiedades físicas debido a su sobrexplotación (compactación).
Los suelos Ferralíticos Rojos son generalmente arcillosos (hasta más de 90% de partículas menores de 0,001 mm) en el subtipo típico, lo que unido a la casi ausencia de minerales alterables y poco contenido de limo, constituye un índice del alto grado de evolución de los mismos. Predominan los minerales arcillosos no dilatables.
Estos suelos tienen un alto contenido de arcilla (igual o más de 90%) y un pequeño contenido de materia orgánica, sin embargo muestran uniformidad física y química en profundidad. Presentan además buen estado de agregación y no son dilatables. Trabajos anteriores en el mismo suelo y lugar de Cid (1995), señalan la existencia de una capa compacta entre los 30 y los 50 cm. Esta capa endurecida es típica en este tipo de suelo debido fundamentalmente a la cantidad y calidad de arcilla predominante y a factores antrópicos.
En este análisis de sensibilidad del modelo AquaCrop se deben considerar como variables importantes la disponibilidad de agua, el grado de fertilidad del suelo y las características hidrofísicas del suelo. Las láminas de riego y los posibles aportes por ascenso capilar modifican de manera importante los valores de producción. De igual forma las disminuciones en el porcentaje de fertilidad del suelo tienen una relación directa con la biomasa y el rendimiento del frijol.
Validación del modelo
En la Tabla 7 se exponen los estadígrafos seleccionados en la validación del modelo en el cultivo de frijol para las variables de salida: biomasa aérea, rendimiento y cobertura vegetal. El seguimiento de la biomasa medida en campo y simulada con el modelo, muestra valores para esta variedad de 7,5 t ha-1 y de 7,83 t ha-1, respectivamente, presentando una diferencia de 330 kg ha-1, lo que representa en términos porcentuales menos del 5%. Los resultados indican que las predicciones de la biomasa del modelo fueron buenas (RMSE=0,82 t ha-1, EF=0,88 y d=0,94).
TABLA 7.
Estadígrafos seleccionados en la validación del modelo
| Estadígrafo | Variables de salida | ||
|---|---|---|---|
| Biomasa | Rendimiento | Cobertura | |
| RMSE | 0,82 | 0,21 | 0,72 |
| EF | 0,88 | 0,94 | 0,90 |
| d | 0,94 | 0,99 | 0,96 |
| r | 0,99 | 0,99 | 0,95 |
| pendiente | 0,87 | 0,90 | 1,03 |
| R2 | 0,96 | 0,97 | 0,92 |
Para las tres variables de salida cobertura del dosel, producción de biomasa y rendimiento, la raíz cuadrada del cuadrado medio del error (RMSE) indica un desempeño bueno en la simulación. La eficiencia está cercana a 1, que es el valor óptimo y d está en el rango establecido (entre 0 y 1). El anterior comportamiento se corrobora a través del coeficiente de correlación con valores de 0,99, 0,99 y 0,95, respectivamente, una relación fuertemente significativa (valores de probabilidad menores a 0,05) entre el seguimiento de las variables simuladas y las variables determinadas a nivel de campo.
Estos resultados demuestran que el modelo AquaCrop explica adecuadamente el desarrollo del crecimiento de frijol, evaluados bajo las condiciones ambientales y edáficas de suelo Ferralítico Rojo compactado en Alquízar, provincia Artemisa.
CONCLUSIONES
Los resultados de la calibración del modelo permitieron optimizar los parámetros fundamentales de suelo y del cultivo de frijol para su aplicación en las condiciones de estudio, destacándose la duración de las fases fenológicas, el índice de cosecha y la curva número.
Se evidencia la posibilidad de alcanzar rendimientos de 2,63 t ha-1, si las precipitaciones son suficiente en este periodo poco lluvioso y el cultivo no se somete a ningún tipo de limitación salvo: la genética vegetal, la radiación solar y la temperatura.
El análisis de sensibilidad del modelo AquaCrop a los cambios en las características medioambientales demostró cómo se pueden afectar las salidas en rendimientos del frijol y que se deben considerar como variables importantes la disponibilidad de agua, la profundidad del nivel freático, el grado de fertilidad del suelo y las características hidrofísicas del suelo.
La validación del modelo demostró la precisión de los parámetros calibrados para la simulación de biomasa, cobertura vegetal y rendimiento, con valores de eficiencia de simulación superiores a 0,88, errores cuadráticos medios menores al 10% y coeficientes de determinación por encima de 0,92.