Saltar al contenido

Evapotranspiración: conceptos, FAO-56 y métodos de cálculo

Pablo RojasActualizado el 13 min de lectura

En este artículo
  1. Evaporación, transpiración y evapotranspiración
  2. Factores que controlan la evapotranspiración
  3. ET0, ETc y ETc ajustada según FAO-56
  4. Ecuación de Penman-Monteith FAO-56
  5. Método de Hargreaves-Samani
  6. Tanque evaporímetro clase A
  7. Coeficiente de cultivo Kc
  8. Medición directa y estimación satelital
  9. Relación con el balance hídrico y el riego
  10. Ejemplo resuelto
  11. Errores frecuentes
  12. Preguntas frecuentes
  13. Fuentes
¿Qué es la evapotranspiración?

La evapotranspiración (ET) es el agua que pasa a la atmósfera como vapor desde el suelo, por evaporación, y desde las plantas, por transpiración. Se expresa como lámina de agua, normalmente en mm/d. El estándar para calcularla es el Estudio FAO Riego y Drenaje n.º 56 (FAO-56): la ecuación de Penman-Monteith entrega la evapotranspiración de referencia y un coeficiente de cultivo la traslada a cada cultivo.

Evaporación, transpiración y evapotranspiración

La evaporación es el paso del agua líquida a vapor y su remoción desde una superficie húmeda: suelo, follaje mojado o agua libre. La transpiración es la vaporización del agua contenida en los tejidos de la planta, que sale por los estomas de las hojas. En un campo ambos procesos ocurren a la vez y se estudian juntos. Según FAO-56, al sembrar casi el 100 % de la ET es evaporación y, con el cultivo cubriendo el suelo, más del 90 % es transpiración. Es la vía por la que el ciclo del agua devuelve a la atmósfera parte de la precipitación continental.

Evaporar agua consume unos 2,45 MJ/kg. Por eso 1 mm/d equivale a 2,45 MJ m⁻² d⁻¹ y a 10 m³/ha/d, y el factor 0,408 de las ecuaciones FAO-56 convierte energía en lámina de agua.

Factores que controlan la evapotranspiración

  • Radiación solar: es la principal fuente de energía; depende de la latitud, la época del año, la nubosidad y el albedo.
  • Temperatura del aire: aporta calor sensible y eleva la presión de vapor de saturación.
  • Humedad del aire: el déficit de presión de vapor mide cuánto vapor más admite el aire; es mayor en climas áridos.
  • Viento: renueva el aire húmedo que queda sobre la superficie; su efecto crece con aire seco y cálido.
  • Agua disponible: cuando baja la humedad del suelo, los estomas se cierran y la ET real cae bajo la demanda del clima.
  • Cultivo y manejo: especie, etapa, altura y cobertura, además de salinidad, densidad y método de riego.

ET0, ETc y ETc ajustada según FAO-56

FAO-56 separa el efecto del clima, del cultivo y del manejo en tres magnitudes:

  • Evapotranspiración de referencia (ET0ET_0): la de un cultivo hipotético de 0,12 m de altura, resistencia superficial de 70 s/m y albedo de 0,23, similar a una pradera de pasto verde, uniforme, en crecimiento activo y sin falta de agua. Depende solo del clima.
  • Evapotranspiración del cultivo en condiciones estándar (ETcET_c): la de un cultivo sano, bien fertilizado, en un campo extenso, con humedad óptima y producción plena. Se calcula como ETc=Kc ET0ET_c = K_c\,ET_0 y es el límite superior del consumo.
  • Evapotranspiración ajustada (ETc ajET_{c\,aj}): la ET real cuando el campo se aparta de esas condiciones por déficit hídrico, salinidad, baja densidad o plagas. Para el déficit hídrico, el capítulo 8 introduce el coeficiente de estrés KsK_s, entre 0 y 1.
ETc aj=Ks Kc ET0ET_{c\,aj} = K_s\,K_c\,ET_0

FAO-56 desaconseja el término «evapotranspiración potencial» por la ambigüedad de sus definiciones.

Ecuación de Penman-Monteith FAO-56

FAO-56 adopta esta ecuación (ecuación 6) como único método estándar para calcular ET0ET_0 con datos meteorológicos. Suma un término de radiación y uno aerodinámico:

ET0=0,408 Δ (Rn−G)+γ 900T+273 u2 (es−ea)Δ+γ (1+0,34 u2)ET_0 = \frac{0{,}408\,\Delta\,(R_n - G) + \gamma\,\dfrac{900}{T + 273}\,u_2\,(e_s - e_a)}{\Delta + \gamma\,(1 + 0{,}34\,u_2)}

donde ET0ET_0 está en mm/d; RnR_n es la radiación neta y GG el flujo de calor en el suelo, ambos en MJ m⁻² d⁻¹; TT, la temperatura media diaria a 2 m (°C); u2u_2, el viento a 2 m (m/s); ese_s y eae_a, las presiones de vapor de saturación y real (kPa); Δ\Delta, la pendiente de la curva de presión de vapor, y γ\gamma, la constante psicrométrica, ambas en kPa/°C. El capítulo 3 de FAO-56 detalla cada término.

Presión de vapor y constante psicrométrica

e∘(T)=0,6108 exp⁡(17,27 TT+237,3),Δ=4098 e∘(T)(T+237,3)2e^{\circ}(T) = 0{,}6108\,\exp\left(\frac{17{,}27\,T}{T + 237{,}3}\right), \qquad \Delta = \frac{4098\,e^{\circ}(T)}{(T + 237{,}3)^2}

Δ\Delta se evalúa con la temperatura media, pero ese_s es el promedio de e∘(Tmax)e^{\circ}(T_{max}) y e∘(Tmin)e^{\circ}(T_{min}). La presión real sale de la temperatura de rocío, ea=e∘(Trocıˊo)e_a = e^{\circ}(T_{\text{rocío}}), o de las humedades relativas extremas:

ea=12[e∘(Tmin) HRmax100+e∘(Tmax) HRmin100],γ=0,665×10−3 Pe_a = \frac{1}{2}\left[e^{\circ}(T_{min})\,\frac{HR_{max}}{100} + e^{\circ}(T_{max})\,\frac{HR_{min}}{100}\right], \qquad \gamma = 0{,}665 \times 10^{-3}\,P

La presión atmosférica PP (kPa) se estima con la altitud zz (m) como P=101,3 [(293−0,0065 z)/293]5,26P = 101{,}3\,[(293 - 0{,}0065\,z)/293]^{5{,}26}.

Radiación neta

La radiación extraterrestre RaR_a depende solo de la latitud y del día del año (ecuaciones 21 a 25). La radiación solar RsR_s se mide con piranómetro o se estima con la fórmula de Ångström, Rs=(0,25+0,50 n/N) RaR_s = (0{,}25 + 0{,}50\,n/N)\,R_a, donde nn son las horas de sol medidas y NN la duración máxima del día. Luego:

Rn=(1−α) Rs−RnlR_n = (1 - \alpha)\,R_s - R_{nl}

con albedo α\alpha = 0,23. La onda larga neta RnlR_{nl} se calcula con las temperaturas extremas, eae_a y la nubosidad relativa Rs/RsoR_s/R_{so}, con RsoR_{so} como radiación de cielo despejado (ecuación 39). En pasos diarios, GG se toma igual a cero.

Viento y datos faltantes

El viento medido a otra altura zz se lleva a 2 m con u2=uz 4,87/ln⁡(67,8 z−5,42)u_2 = u_z\,4{,}87/\ln(67{,}8\,z - 5{,}42), y no conviene usar u2u_2 menor que 0,5 m/s. Si falta la radiación, Rs=kRsTmax−Tmin RaR_s = k_{Rs}\sqrt{T_{max} - T_{min}}\,R_a, con kRsk_{Rs} ≈ 0,16 en el interior y ≈ 0,19 en la costa. Si falta la humedad, ea≈e∘(Tmin)e_a \approx e^{\circ}(T_{min}) en estaciones rodeadas de vegetación regada.

Método de Hargreaves-Samani

Hargreaves y Samani (1985) propusieron una ecuación que usa solo las temperaturas extremas y la radiación extraterrestre, calibrada con ocho años de datos diarios de lisímetros con festuca en Davis, California. FAO-56 la recoge como ecuación 52:

ET0=0,0023 (Tmedia+17,8) Tmax−Tmin  RaET_0 = 0{,}0023\,(T_{media} + 17{,}8)\,\sqrt{T_{max} - T_{min}}\;R_a

con ET0ET_0 y RaR_a en mm/d: RaR_a en MJ m⁻² d⁻¹ se multiplica por 0,408. La amplitud térmica actúa como indicador de nubosidad y humedad. Conviene usarla cuando la estación registra solo temperaturas, para rellenar series o en pasos de diez días o mensuales. FAO-56 advierte que subestima con vientos fuertes (u2>3u_2 \gt 3 m/s) y sobrestima con humedad relativa alta. También pide verificarla en cada región contra Penman-Monteith y señala que Penman-Monteith con RsR_s y eae_a estimados suele ser algo más precisa.

Tanque evaporímetro clase A

El tanque clase A es un recipiente circular de 120,7 cm de diámetro y 25 cm de profundidad, montado sobre una plataforma a 15 cm del suelo. Como el agua libre no evapora igual que un cultivo, su evaporación EpanE_{pan} (mm/d) se corrige con un coeficiente:

ET0=Kp EpanET_0 = K_p\,E_{pan}

La tabla 5 de FAO-56 da KpK_p entre 0,35 y 0,85 según si el tanque está rodeado de cultivo verde (caso A) o de barbecho seco (caso B), la distancia de esa cubierta a barlovento, el viento y la humedad relativa media. En el caso A, con 10 m de cubierta, viento de 2 a 5 m/s y humedad de 40 a 70 %, KpK_p = 0,70: una evaporación media de 5,5 mm/d da ET0ET_0 = 3,9 mm/d. El método se recomienda para períodos de diez días o más, y un nivel de agua 10 cm por debajo de la norma (5 a 7,5 cm bajo el borde) puede generar errores de hasta 15 %.

Coeficiente de cultivo Kc

El coeficiente de cultivo simple KcK_c resume las diferencias entre el cultivo y el pasto de referencia: altura, resistencia de las hojas, albedo y evaporación del suelo. FAO-56 divide la temporada en cuatro etapas: inicial, de la siembra a cerca de 10 % de cobertura; desarrollo, hasta la cobertura completa efectiva; media, hasta el inicio de la madurez, y final, hasta la cosecha o la senescencia completa. KcK_c es constante en las etapas inicial (Kc iniK_{c\,ini}) y media (Kc medK_{c\,med}), y varía linealmente en las de desarrollo y final hasta Kc finK_{c\,fin}. La tabla 12 entrega valores para clima subhúmedo, con HRminHR_{min} ≈ 45 % y u2u_2 ≈ 2 m/s:

CultivoKc iniKc medKc finAltura máxima (m)
Maíz para grano0,31,200,60 a 0,352
Papa0,51,150,750,6
Tomate0,61,150,70 a 0,900,6
Poroto o frijol seco0,41,150,350,4
Vid de mesa o para pasas0,300,850,452
Palto (aguacate), sin cubierta vegetal0,600,850,753

El Kc finK_{c\,fin} del maíz es 0,60 si se cosecha con grano húmedo y 0,35 si el grano se seca en el campo. Los Kc iniK_{c\,ini} suponen riegos de frecuencia habitual; con aspersión muy frecuente o lluvia diaria pueden llegar a 1,0 o 1,2. En climas más secos o ventosos, Kc medK_{c\,med} (y Kc finK_{c\,fin} cuando es de 0,45 o más) se corrige según la altura del cultivo hh (m):

Kc med=Kc med(tab)+[0,04 (u2−2)−0,004 (HRmin−45)](h3)0,3K_{c\,med} = K_{c\,med(tab)} + \left[0{,}04\,(u_2 - 2) - 0{,}004\,(HR_{min} - 45)\right]\left(\frac{h}{3}\right)^{0{,}3}

La corrección vale para u2u_2 entre 1 y 6 m/s y HRminHR_{min} entre 20 y 80 %. Para riegos de alta frecuencia, FAO-56 propone el coeficiente dual, que separa la transpiración (KcbK_{cb}) de la evaporación del suelo (KeK_e).

Medición directa y estimación satelital

Los lisímetros aíslan un volumen de suelo con cultivo. Los de pesada miden la ET por cambio de masa, con precisión de centésimas de milímetro; los de drenaje, por diferencia entre agua aplicada y drenada. Son costosos y se reservan para investigación y calibración. Los métodos micrometeorológicos derivan el calor latente del balance de energía Rn−G−λET−H=0R_n - G - \lambda ET - H = 0: la razón de Bowen usa gradientes de temperatura y humedad, y la covarianza de remolinos correlaciona las fluctuaciones rápidas del viento vertical y del vapor. En 22 sitios de la red FLUXNET, Wilson et al. (2002) hallaron un desbalance de energía medio cercano a 20 %, atribuible a flujos turbulentos subestimados o a energía disponible sobrestimada.

Los modelos satelitales calculan la ET como residuo del balance de energía a partir de la temperatura de superficie. METRIC (Allen, Tasumi y Trezza, 2007), derivado de SEBAL, usa la banda térmica de Landsat, píxeles de anclaje frío y caliente y la ET de referencia de estaciones terrestres. OpenET combina seis modelos con píxeles de 30 m; contra torres de flujo en 24 sitios agrícolas, sus errores medios absolutos fueron de 13,6 a 21,6 mm/mes (Melton et al., 2022).

Relación con el balance hídrico y el riego

En el balance hídrico de una cuenca, la ET real cierra la ecuación entre precipitación, escorrentía y cambio de almacenamiento. En la zona de raíces, FAO-56 lleva un balance diario del agotamiento del suelo. Mientras este no supera el agua fácilmente aprovechable, una fracción pp (cercana a 0,5 en muchos cultivos) del agua disponible entre capacidad de campo y punto de marchitez, el cultivo evapotranspira a la tasa ETcET_c; más allá, KsK_s baja linealmente hasta cero.

La ETcET_c de la etapa media fija el caudal de punta de los sistemas de riego y de los canales de riego, una vez sumadas las pérdidas. La diferencia entre precipitación y evapotranspiración es la base del SPEI, uno de los índices de sequía, y la ET del agua de lluvia corresponde al agua verde de la huella hídrica. En Chile, la Red Agrometeorológica del INIA publica la ET0ET_0 de sus estaciones.

Ejemplo resuelto

Se reproduce el ejemplo 18 de FAO-56: 6 de julio en Uccle (Bruselas), a 50°48' N y 100 m de altitud, con TmaxT_{max} = 21,5 °C, TminT_{min} = 12,3 °C, HRmaxHR_{max} = 84 %, HRminHR_{min} = 63 %, viento de 2,78 m/s medido a 10 m y 9,25 h de sol. Los cálculos se hicieron con python3.

1. Penman-Monteith FAO-56

MagnitudEcuaciónResultado
TT; PP; γ\gamma9; 7; 816,9 °C; 100,1 kPa; 0,0666 kPa/°C
Δ\Delta130,122 kPa/°C
u2u_2472,078 m/s
ese_s; eae_a; es−eae_s - e_a11, 12 y 171,997; 1,409; 0,589 kPa
RaR_a; NN (día 187)21 y 3441,09 MJ m⁻² d⁻¹; 16,1 h
RsR_s; RsoR_{so}35 y 3722,07; 30,90 MJ m⁻² d⁻¹
RnsR_{ns}; RnlR_{nl}; RnR_n38 a 4017,00; 3,71; 13,28 MJ m⁻² d⁻¹
Término de radiación62,81 mm/d
Término aerodinámico61,07 mm/d
ET0ET_063,88 mm/d

Los valores intermedios coinciden con los de FAO-56, que informa 3,9 mm/d.

2. Hargreaves-Samani

RaR_a = 41,09 × 0,408 = 16,76 mm/d, Tmedia+17,8T_{media} + 17{,}8 = 34,7 °C y 9,2\sqrt{9{,}2} = 3,033. Entonces ET0=0,0023×34,7×3,033×16,76ET_0 = 0{,}0023 \times 34{,}7 \times 3{,}033 \times 16{,}76 = 4,06 mm/d, 4,6 % sobre Penman-Monteith, coherente con su tendencia a sobrestimar con aire húmedo (humedad media de 73,5 %).

3. ETc de maíz en la etapa media

PasoCálculoResultado
Kc med de tablaMaíz para grano, h = 2 m1,20
Ajuste climático[0,04 × 0,078 − 0,004 × 18] × (2/3)^0,3−0,061
Kc med ajustado1,20 − 0,0611,14
ETc1,14 × 3,88 mm/d4,4 mm/d
Volumen semanal4,42 mm/d × 7 d × 10 m³/(ha·mm)309 m³/ha

Sin ajuste climático sería de 4,66 mm/d. Con déficit hídrico, la ET real sería KsK_s veces ese valor.

Errores frecuentes

  • Usar RaR_a en MJ m⁻² d⁻¹ en Hargreaves: sin el factor 0,408, la ET0ET_0 sale 2,45 veces mayor.
  • Calcular ese_s con la temperatura media: en el ejemplo, e∘(16,9)e^{\circ}(16{,}9) = 1,925 kPa y el déficit de presión de vapor baja 12 %.
  • No llevar el viento a 2 m: los 2,78 m/s medidos a 10 m superan en 34 % al u2u_2 correcto.
  • Tomar ETcET_c como ET real: con déficit hídrico, salinidad o baja densidad hay que aplicar KsK_s u otros ajustes.
  • Usar KcK_c de tabla sin corregir: Kc medK_{c\,med} cambia con el viento y la humedad, y Kc iniK_{c\,ini} con la frecuencia de mojado.
  • Aplicar el tanque con KpK_p = 1 o a datos diarios: KpK_p va de 0,35 a 0,85 y el método es confiable en períodos de diez días o más.

Preguntas frecuentes

¿Cuánta agua representa 1 mm de evapotranspiración?

Un milímetro equivale a 1 litro por metro cuadrado, o 10 m³ por hectárea. Un cultivo que evapotranspira 6 mm/d consume 60 m³/ha cada día.

¿Qué método usar si la estación solo registra temperaturas?

Estimar RsR_s con la amplitud térmica y eae_a con la temperatura mínima y aplicar Penman-Monteith, como indica FAO-56. Hargreaves-Samani sirve para pasos de diez días o mensuales si se verifica localmente.

¿Por qué la ETc puede ser mayor que la ET0?

Los cultivos más altos y rugosos que el pasto de referencia, con más área foliar, intercambian vapor con la atmósfera con mayor eficiencia. Por eso el maíz o la papa superan KcK_c = 1 en la etapa media.

¿Se puede medir la evapotranspiración con satélites?

Se estima, no se mide: los modelos resuelven el balance de energía con la temperatura de superficie y datos meteorológicos. Entregan mapas a escala de predio, pero requieren validación con torres de flujo o lisímetros.

Fuentes

  • Allen, R. G., Pereira, L. S., Raes, D. y Smith, M. (1998). Crop evapotranspiration: Guidelines for computing crop water requirements (Estudio FAO Riego y Drenaje n.º 56). FAO. fao.org
  • Hargreaves, G. H. y Samani, Z. A. (1985). Reference crop evapotranspiration from temperature. Applied Engineering in Agriculture, 1(2), 96-99. doi.org
  • Allen, R. G., Tasumi, M. y Trezza, R. (2007). Satellite-based energy balance for mapping evapotranspiration with internalized calibration (METRIC): Model. Journal of Irrigation and Drainage Engineering, 133(4), 380-394. ascelibrary.org
  • Wilson, K., Goldstein, A., Falge, E. et al. (2002). Energy balance closure at FLUXNET sites. Agricultural and Forest Meteorology, 113, 223-243. ameriflux.lbl.gov
  • Melton, F. S., Huntington, J., Grimm, R. et al. (2022). OpenET: Filling a critical data gap in water management for the western United States. Journal of the American Water Resources Association, 58(6), 971-994. usgs.gov
  • Instituto de Investigaciones Agropecuarias (INIA). Agrometeorología: evapotranspiración de referencia. agrometeorologia.cl