Paola A.

Arias
MÉTODOS DE ESTIMACIÓN DE Adriana Cuartas
LA EVAPORACIÓN EN COLOMBIA Hernán A. Moreno
Janeth Barco
(RESÚMEN) Jaime I. Vélez
Grupo Atlas Hidrológico de Colombia Oscar J. Mesa
Germán Poveda

EVAPOTRANSPIRACIÓN Y SU ESTIMACIÓN

La evapotranspiración es una componente fundamental en el balance hidrológico, aunque

es, a su vez, la variable con mayores incertidumbres en su medición y estimación. Está
compuesta por la evaporación desde superficie y por la transpiración de las plantas. Debido
a que los aportes de la evaporación y de la transpiración son difíciles de cuantificar, surge
la necesidad de estimarlos conjuntamente empleando variables que influyen de alguna
manera en el comportamiento de ambas componentes.

La dinámica de la evapotranspiración está condicionada por tres factores fundamentales:

disponibilidad de agua, disponibilidad de energía para cambio de estado del agua y
condiciones aerodinámicas propicias para asimilación y transporte de vapor de agua.

La disponibilidad de agua en una zona es determinada por la precipitación. A largo plazo,

ésta es la responsable del flujo base de las corrientes, del suministro de agua a los acuíferos
por infiltración y en la zona vadosa para suministro de agua a las plantas, las cuales, a su
vez, interceptan gran parte del agua precipitada.
La radiación neta es el remanente del balance energético global, el cual está afectado por
factores astronómicos y por condiciones físicas y climáticas. La disponibilidad de energía
para cambio de estado del agua es suministrada por la radiación neta proveniente del sol.
De esta radiación, una parte es responsable del calentamiento de la superficie en forma de
calor sensible, que es reflejada por la temperatura superficial, mientras que la otra parte es
la responsable de proveer calor latente de vaporización.
Las condiciones aerodinámicas propicias para la evaporación son influenciadas por la
velocidad del viento, la cantidad de vapor de agua en el aire, la resistencia aerodinámica de
la superficie y por el tipo y cantidad de vegetación presentes.
Se distinguen dos tipos de evapotranspiración: la evapotranspiración potencial y la
evapotranspiración real. La evapotranspiración potencial es la máxima evapotranspiración
posible considerando que las únicas condiciones limitantes en el fenómeno son las
radiativas, esto es, la única limitante para la evaporación es la disponibilidad de energía. La
evapotranspiración real es la evapotranspiración presente considerando las condiciones,
tanto de disponibilidad de energía como de disponibilidad hídrica y características
aerodinámicas.

Debido a la falta de mediciones directas de evapotranspiración es necesario recurrir a

técnicas de estimación de esta variable a diferentes escalas temporales y espaciales.
 2

Este escrito realiza una recopilación de metodologías empleadas para la estimación de la

evapotranspiración, tanto potencial como real. La mayoría de estos procedimientos se
emplearon en el proyecto Balances Hidrológicos de Colombia (UNALMED-
COLCIENCIAS-UPME, 1999) y en Barco y Cuartas (1999), aunque se ha trabajado en el
desarrollo de nuevas metodologías con mayores consideraciones acerca de la física de este
fenómeno.

EVAPOTRANSPIRACIÓN POTENCIAL

Ecuación de Cenicafé
Esta ecuación fue obtenida por el Centro Nacional de Investigaciones del Café
(CENICAFÉ) mediante regresión a los valores obtenidos de aplicar el método de Penman a
los datos de sus estaciones climáticas en Colombia para la obtención de evaporación
potencial (Jaramillo, 1989 y 1997 comunicación personal). Su cálculo es sencillo, ya que
sólo depende de la altura sobre el nivel del mar, así

ETP =4.658 exp (− 0.0002 h ) (1)

ETP es la evaporación potencial en mm/día y h es la cota sobre el nivel del mar en m.

Ecuación de García y López

Por medio de esta ecuación se calcula la evaporación potencial. Está basada en las
correlaciones hechas entre la evaporación medida en seis estaciones netamente tropicales
dentro del rango latitudinal 15°N y 15°S y la temperatura y la humedad relativa, obteniendo
la ecuación

ETP = 1.21* 10 N ( 1 − 0.01HR ) + 0.21T − 2.3 (2)

la cual es válida sólo para temperaturas medias mayores o iguales a 10 °C. El valor de HR
se estima mediante
7.45T
n= (3)
(234.7 + T)

HR8:00 HORAS + HR14 HORAS

HR = (4)
2
Esta expresión estima ETP en mm/día.

Ecuación de Thornwaite
La ecuación de Thornwaite (1944), da una estimación de la evaporación potencial basada
en numerosos experimentos efectuados con lisímetros. Ésta fue modificada por el mismo
 3

autor en 1948, cuando relacionó datos de temperatura y evapotranspiración potencial en

regiones irrigadas del oeste de Estados Unidos.

Los factores tomados dentro de este cálculo son la temperatura media mensual y la
temperatura media anual
a
⎛ T⎞
ETP = 1.6 ⎜ 10 ⎟ (5)
⎝ I⎠
donde ETP es la evaporación potencial en cm/mes, T es la temperatura media mensual en
°C, I es el índice calórico anual dado por
⎛ ⎛ Tanual ⎞1.514 ⎞
I = 12 ⎜ ⎜ ⎟ ⎟ (6)
⎜⎝ 5 ⎠ ⎟
⎝ ⎠
y a es un exponente dado en función de I, así

( ) ( ) ( )
a = 675 ⋅ 10 −9 I 3 − 771 ⋅ 10 −7 I 2 + 179 ⋅ 10 −4 I + 0.492 (7)

Turc Modificado
Una vez obtenida la ecuación para el cálculo de la evaporación potencial anual de una
cuenca vertiente, Turc (1962) encontró una expresión basada en las observaciones de las
cajas lisimétricas de Versailles, Rothamsted, Groningen y otras estaciones situadas en
Escocia, África del Norte y las Indias. La expresión obtenida considera la evaporación real
en función de la radiación global incidente, la temperatura, la humedad relativa y una
constante que depende del mes o período considerado.
La ecuación está dada como
⎛ T ⎞
ETP = K ⎜ (
⎟ R g + 50 ) si HR>50% (8)
⎝ T + 15 ⎠

⎛ 50 − H R ⎞
⎟ (R g + 50 ) ⎜ 1 +
⎛ T ⎞
ETP = K ⎜ ⎟ si HR <50% (9)
⎝ T + 15 ⎠ ⎝ 70 ⎠

donde HR es la humedad relativa media mensual.

⎛ 50 − H R ⎞
En la práctica, el término de corrección ⎜ 1 + ⎟ interviene sólo en el caso de
⎝ 70 ⎠
climas desérticos o subdesérticos.
ETP es la evapotranspiración potencial en mm/ mes, K es una constante igual a 0.40 para
meses de 30 y 31 días, 0.37 para el mes de febrero y 0.13 para 10 días, T es la temperatura
media mensual en grados centígrados, Rg es la radiación solar global incidente del mes
considerado expresada en cal/ cm2 /día. (Remenieras, 1974).
 4

Método combinado aerodinámico y de balance de energía (Penman)

En el caso de la estimación de la evaporación potencial, las dos limitantes principales son el

suministro de energía y el transporte del vapor. Bajo esta condición, es necesario la
combinación de estos dos aspectos para dar un valor estimado de la evaporación potencial.
Esta ecuación fue desarrollada inicialmente por Penman en 1948 y está dada por

∆ γ
E= Er + Ea (10)
∆ +γ ∆ +γ

En la ecuación para el cálculo de la evaporación γ es la constante psicrométrica, la cual se

expresa como
CpKhP
γ= (11)
0.622 Lv K w

Cp es el calor específico a presión constante; este valor es 1004 J/Kg.K; P es la presión

atmosférica en Pa; Lv es el calor latente de vaporización dado por

Lv = 2.501* 10 6 − 2.370T (12)

expresado en J/Kg y T es la temperatura del aire en °C. Kh/Kw es la relación entre las
difusividades de calor y el vapor que se toma comúnmente como 1 (Priestley y Taylor,
1972). ∆ es la pendiente de la curva de presión de saturación del vapor a la temperatura del
aire en Pa/°C expresada como
es
∆ = 4098 (13)
( 237.3 + T ) 2

es es la presión de vapor de saturación, en Pa, a la temperatura del aire y está definido

como
⎛ αT ⎞
e s = 611 exp⎜⎜ ⎟⎟ (14)
⎝ β +T ⎠

con α = 17.27 y β = 237.3 si T ≥ 0, ó , α = 21.88 y β = 265.5 si T ≤ 0. T está en grados

Celsius.
Er (m/s) representa la componente radiativa a la evaporación, dada por

Rn
Er = (15)
Lv ρ w

donde Rn es la radiación neta en w/m2; Lv es el calor latente de vaporización del agua que
se calcula como se dijo anteriormente y ρw es la densidad del agua tomada como 1000
Kg/m3.
 5

Además del suministro de energía calórica, la tasa de evaporación desde una superficie
abierta de agua está controlada también por la habilidad para transportar el vapor lejos de la
superficie. La tasa de transporte se determina por el gradiente de humedad en el aire
cercano a la superficie y la velocidad del viento a través de dicha superficie. Mediante un
análisis matemático, aplicando las ecuaciones de transporte de masa y momentum en el
aire, se obtiene la expresión para la evaporación en su componente aerodinámica

E a = B( e s − e ) (16)

Ea está en m/s, es es la presión de vapor de saturación correspondiente a la temperatura del

aire en Pa, T es la temperatura en grados Celsius y e es la presión de vapor, en Pa,
expresada en términos de la temperatura del punto de rocío Td en grados celsius

⎛ α .Td ⎞
e = 611 exp⎜⎜ ⎟⎟ (17)
⎝ β + Td ⎠

B es el coeficiente de transferencia de vapor y está dado por

0.622 k 2 ρ a u 2
B= 2
(18)
⎡ ⎛z ⎞⎤
Pρ w ⎢ Ln⎜⎜ 2 ⎟⎟⎥
⎣⎢ ⎝ z o ⎠⎦⎥

k es la constante de Von Karman, aproximadamente 0.4, ρa es la densidad del aire en kg/m3

u2 es la velocidad del viento en m/s medida a una altura z2 en cm, P es la presión
atmosférica a una altura dada, en Pa, ρw es la densidad del agua, en kg/m3 y zo es la altura
de la rugosidad en cm, la cual puede ser expresada como Winkler (1993)

z o = 0.13h (19)

donde h, es la altura de la vegetación que se puede estimar mediante

h = 0.4 NDVI + 0.04 (20)

El NDVI es el índice normalizado de vegetación, tomado de información de satélite.

Método combinado aerodinámico y de balance de energía (Penman) con aproximación

de Priestley y Taylor

Priestley y Taylor, determinaron, para muchas regiones, que el segundo término de la

ecuación podría representar una fracción del primero, de la siguiente forma
∆
E = 1.3 Er (21)
∆+ γ
 6

Donde la componente aerodinámica ha sido aproximada como el 30% de la componente

radiativa.

Ecuación de Meyer
Está ecuación, que calcula la evaporación potencial, fue propuesta en 1915 y está basada en
la ley de Dalton.
⎛ u ⎞
ETP = C 1 ( e s − e a )⎜ + 1⎟ (22)
⎝ 16.09 ⎠

donde ETP es la evaporación potencial en cm; es, la presión de vapor de saturación media
en pulgadas de mercurio; ea la presión de vapor media mensual en pulgadas de mercurio y
u, la velocidad media mensual del viento medida a 10 m de la superficie en km/h. C1 es un
coeficiente empírico que puede tomarse como 38 para depósitos pequeños y 28 para
grandes depósitos (Remenieras, 1974).

EVAPOTRANSPIRACIÓN REAL

Método de Penman-Monteith

Esta ecuación está basada en el método combinado aerodinámico y de balance de energía o

Penman (1948 y 1946). El método permite calcular la evaporación con variables
meteorológicas como radiación neta, temperatura y humedad relativa, además de estimar
otras relacionadas con las características de los estomas y la aerodinámica del cultivo que
son calculadas con la altura media del cultivo, la velocidad del viento y con el índice de
área de hojas. La ecuación tiene la forma:

⎛ ρC P ⎞
∆(R N − G ) + ⎜⎜ ⎟⎟(e − e d )
⎝ ra ⎠
ETR = (23)
⎛ ⎛ r ⎞⎞
Lv ⎜⎜ ∆ + γ⎜⎜ 1 + s ⎟⎟ ⎟⎟
⎝ ⎝ ra ⎠ ⎠

donde ETR es la evapotranspiración real en mm/dia, Lv es el calor latente de vaporización

definido anteriormente, ∆ es el gradiente de la presión de vapor de saturación con respecto
a la temperatura, expresado en Pa/C°, Rn es la radiación neta, G es la densidad de flujo de
calor del suelo, ra es la resistencia aerodinámica en s/m evaluada como

2
⎛ ⎛ 2 ⎞⎞
4.72⎜⎜ Ln⎜⎜ ⎟⎟ ⎟
⎟
⎝ ⎝ Z0 ⎠⎠
ra = (24)
(1 + 0.54.u )
 7

Z0 es la altura de la rugosidad de la superficie expresada en m; u es la velocidad del viento

a 2 m de altura de la superficie en m/s; rs es la resistencia de la superficie en s/m expresada
como,

200
rs = (25)
LAI

Usualmente en el trópico G puede ser despreciada.

Ecuación de Turc
Esta ecuación estima la evapotranspiración real basándose en un balance de masas, en
función de elementos meteorológicos simples como la temperatura y la precipitación de la
cuenca, aplicadas a medidas de largo plazo. Turc en 1954, adaptó la familia de curvas de la
forma ETR = f (P,T), a partir de observaciones hechas en 254 cuencas situadas en todos los
climas de la tierra (Remenieras, 1974; Brutsaert, 1982). De esta manera obtuvo

P
ETR = para (P / L) > 0.316 (26)
2
P
0 .9 +
L2

ETR = P para (P / L) <0.316 (27)

Donde ETR es la evapotranspiración real en mm/año, P es la precipitación del año

considerado en mm, T es la temperatura media anual en grados centígrados y L = 300 + 25
T + 0.05 T3. Este último parámetro fue hallado por medio de regresiones hechas en las 254
cuencas y corresponde a la evapotranspiración potencial, ETP, en mm/año.

Ecuación de Coutagne

Esta ecuación se basa en las características meteorológicas de la cuenca. Da una estimación

de la evaporación real en función de la precipitación y la temperatura así,
ETR = P -λP2 (28)
Donde ETR es la evaporación real en m / año, P es la precipitación en m, y
1
λ= (29)
0.8 + 0.14T
en donde T es la temperatura media anual en °C. La ecuación (28) sólo es aplicable para P
comprendido entre (8λ)-1 y (2λ)-1. Si las precipitaciones son inferiores a (8λ)-1 entonces
ETR=P. Si P > (2λ)-1, entonces ETP = (4λ ) = 0.20 + 0.035T (Remenieras, 1974).
−1
 8

Método de Morton
Morton (1983), basado en la ecuación de balance de energía postula con base en la
hipótesis de Bouchet (Bouchet, 1983), una relación complementaria entre la evaporación
potencial ETP, la evaporación real ETR y la evaporación de un ambiente húmedo ETW.

En términos de ecuaciones esta relación se puede escribir,

ETR = 2 ETW − ETP (30)

donde ETR es la evaporación real, ETW es la evaporación de un ambiente húmedo y ETP

es la evaporación potencial, dadas en W/m2. La relación complementaria también puede
expresarse gráficamente así,

Figura 1 Representación esquemática de la relación complementaria entre

evaporación potencial y real con suministro de energía de radiación
constante.

La Figura 1 es una representación esquemática de la ecuación (30) bajo condiciones de

suministro de radiación constantes. La abscisa representa el suministro de agua a la
superficie suelo-planta de un área, una cantidad que es usualmente desconocida. Cuando
no hay disponibilidad de agua para la evaporación real (extremo izquierdo de la Figura 1)
resulta que ETR es igual a cero, el aire es muy caliente y seco y ETP alcanza el máximo
valor (2ETW, la evaporación potencial de un ambiente seco). Como el suministro de agua
para la superficie suelo-planta del área se incrementa (moviéndose a la derecha de la Figura
1 el resultado equivalente en el incremento de ETR causa que el aire que está pasando por
encima se vuelve frío y más húmedo, el cual a su vez produce un decremento equivalente
en ETP. Finalmente cuando el suministro de agua para la superficie suelo-planta del área
se ha incrementado suficientemente los valores de ETR y ETP convergen al mismo punto.
La definición convencional de la evaporación potencial es la misma definición que la
evaporación real de un ambiente húmedo.
 9

Como puede verse la relación complementaria de Morton (1982) es muy similar a la

propuesta por Bouchet en 1963.
ETW se obtiene de la expresión
⎛ ⎞
⎜ ⎟
⎜ R np ⎟ 1
ETW = b1 + b2 ⎜ ETP ≤ ETW ≤ ETP (31)
γ p⎟
;
⎜ 1+ ⎟
2
⎜ ∆ p ⎟⎠
⎝
b1 y b2 son constantes; b1 = 14 W/m2, b2 = 1.2, ∆P es la pendiente de la presión de vapor de
saturación a temperatura de equilibrio, Rnp es la radiación neta de la superficie del suelo a
temperatura de equilibrio dada por,

(
Rnp = ETP + γ p f T T p − T ) (32)

En la ecuación (32) γ es la constante psicrométrica ≈ 0.66 hPa/oC, p es la presión

atmosférica en hPa, TP es la temperatura de equilibrio en °C; esta es la temperatura a la cual
la ecuación de balance de energía y la ecuación de transferencia de vapor para una
superficie húmeda da el mismo resultado, T es la temperatura promedio del aire en °C, ETP
es la evapotranspiración potencial, en W/m2 expresada como,

(
ETP = Rn − λf T T p − T ) (33)
donde Rn es la radiación neta, en W/m2 para la superficie del suelo a temperatura del aire
dada por,
Rn = (1 − α ) G − B p (34)

En la ecuación (34) α es el albedo, G es la radiación global incidente; esta es trasmitida a

través de las nubes, Bp es la pérdida de radiación neta de onda larga por la superficie de
suelo a la temperatura del aire dada como
⎡ ⎛ p ⎞ ⎤
B p = εσ (T + 273)4 ⎢1 − ⎜⎜ 0.71 + 0.070V D ⎟⎟ (1 + φ )⎥ (35)
⎢⎣ ⎝ pS ⎠ ⎥⎦

⎛ p ⎞
El término ⎜⎜ 0.71 + 0.070V D ⎟ es la relación de la radiación de la atmósfera a cielo
⎝ p S ⎟⎠
descubierto, φ es el incremento proporcional de la radiación en la atmósfera debido a las
nubes, expresada como

[
φ = 0.18 (1 − c 2 )(1 − S )2 + c 2 (1 − S )0.5 ] pp
S
(36)
 10

⎛e ⎞
c 2 = 10⎜ d − S − 0.42 ⎟ 0 ≤ c 2 ≤ 1.0 , S es la relación de horas observadas de brillo solar y
⎝ e ⎠
el máximo posible, εσ = 5.22*10-8W/m2K4, ε es la emisividad de la tierra, tal que ε =0.92,
pS es la presión atmosférica a nivel del mar. e d es la presión de vapor de saturación en hPa.

⎛ 17.27Td ⎞
e d = 6.11 exp⎜⎜ ⎟⎟ (37)
⎝ Td + 273 ⎠
Td es la temperatura promedio a punto de rocío en oC.

En la ecuación (34), en el trópico, Bp ≥ 0.05εσ (T + 273)4 mientras que en la ecuación (33)

λ es el coeficiente de transferencia de calor, estimado como
⎛ 4εσ (T + 273 )3 ⎞
λ =γ p+⎜ ⎟ (38)
⎜ fT ⎟
⎝ ⎠
γ p = γ pS (p/pS ) = (0.66 hPa/oC)(p/pS)., fT es el coeficiente de transferencia de vapor
estimado como
0.5
⎛p ⎞ ⎛ fZ ⎞
f T = ⎜⎜ S ⎟⎟ ⎜⎜ ⎟⎟ (39)
⎝ p ⎠ ⎝ ξ ⎠
fZ = 28 W/m2mbar, ξ es el factor de estabilidad adimensional dado como

1 ⎛ e ⎞ ∆R nc
= 0.28⎜ 1 + d ⎟ + (40)
ξ ⎝ e ⎠ ⎛p ⎞
0.5
γ p⎜⎜ S ⎟⎟ bO f Z (e − e d )
⎝ p ⎠
bO = 1, Rnc = Rn con Rnc = > 0

Para estimar Tp , VP, ∆P, es necesario hacer un proceso iterativo como sigue:
Se debe escoger inicialmente TP’, VP’, y ∆P’ igual TP, Vp, ∆P e iterar
 11

⎡ Rn ⎤
⎢ + ed − e P ' +λ (T − TP ' )⎥
f
δTp' = ⎣ T ⎦
(∆ P ' +λ )
TP = TP ' +δTp'
(41)
⎡ α T ⎤
e P = 6.11 exp ⎢ 1 P ⎥
⎣ (T P + β )⎦
α βV
∆P = 1 P 2
(TP + β )

Se itera hasta obtener⎮δTP⎮ ≤ 0.01 oC.

Método de Jiang – Islam

El método propuesto por Jiang y Islam (1999) es un esquema simple que estima
evapotranspiración real a partir de información obtenida de sensores remotos. Esta técnica
se basa en la relación existente entre algunos parámetros superficiales que son de fácil
medición (temperatura, radiación neta y vegetación) y un parámetro índice de la resistencia
efectiva de la superficie a la evaporación.
Dentro de las metodologías para la estimación del flujo de calor latente en superficie, una
de las más comúnmente empleadas es el método residual, el cual está basado en el flujo de
energía en la superficie, dado por

R n − G = H + λE (42)

donde Rn es la radiación neta, G es el flujo de calor en el suelo, H es el flujo de calor

sensible y λE es el flujo de calor latente.
Algunas de las formulaciones empleadas para la estimación de λE pueden ser expresadas
de la forma
⎡ ∆ γ ⎤
λE = β ⎢ A( R n − G ) +B f ( u )( e * a −e a )⎥ (43)
⎣ ∆ +γ ∆ +γ ⎦

donde ea es la presión de vapor en la atmósfera, e*a es la presión de vapor saturado en la

atmósfera, ∆ es la pendiente de la presión de vapor saturado a la temperatura del aire Ta, γ
es la constante psicrométrica, f(u) es una función del viento en la que los parámetros A y B
dependen del modelo particular y β es el parámetro de Budyko-Thornwaite-Mather (ver
Brutsaert, 1982; Parlange et al., 1995).

Nótese que la expresión (43) puede aproximarse a la expresión de Priestley y Taylor

empleando A=α, β=1.0 y B=0.
 12

Observaciones provenientes de FIFE, HAPEX-MOBILHY y ARME indican la existencia

de una fuerte relación entre el parámetro α y características físicas como humedad,
temperatura y conductancia de la superficie. Basados en esta evidencia y en la simplicidad
de la expresión de Piestley y Taylor, Jiang y Islam propusieron un modelo de evaporación
dado por
⎡ ∆ ⎤
λE = φ ⎢( R n − G )
∆ + γ ⎥⎦
(44)
⎣
con λE en w/m2.

La información de radiación neta (Rn) puede ser obtenida de sensores remotos. Aunque en
ocasiones se pueden tener estimaciones de G, estos autores proponen estimar esta variable
como
G = 0.583 exp( −2.13 NDVI )R n si NDVI>0

G = 0.583 R n si NDVI<0

En la expresión (44), φ representa un complejo parámetro que resume los efectos

combinados de los parámetros de Priestley y Taylor, α y β. En ausencia de advección y
convección significativas, el flujo de calor latente no puede exceder Rn-G, y por tanto, φ
∆ +γ
está en el rango comprendido entre 0 y .
∆
En la Figura 2 se observa el comportamiento de la temperatura superficial sensada por
satélite (To) con el índice normalizado de vegetación (NDVI) obtenido por Jiang-Islam.
Vale la pena aclarar que la información con la se obtuvieron estos resultados corresponde a
imágenes procesadas de un paso del satélite AVHRR NOAA-14 sobre la región de trabajo
del proyecto Southern Great Plains 1997 (SGP97), por lo que se observa una alta
variabilidad de los valores de temperatura.

Figura 2 Temperatura derivada (To) vs NDVI. Información requerida para la

interpolación del parámetro φ. Tomada de Jiang-Islam (1999).
 13

Puede observarse que para valores de NDVI negativos se presentan las más bajas
temperaturas, mientras que para valores positivos, los valores de temperatura se ven
acotados por un triángulo, de manera que a medida que aumenta el NDVI se observan
menores temperaturas superficiales.
La interpolación del parámetro φ se basa en el hecho de que para un valor de NDVI
cualquiera, su valor mínimo, φmin=0, está asociado a una temperatura máxima, Tmax,
∆ +γ
mientras que su valor máximo, φmax= , está asociado a una temperatura mínima, Tmin.
∆
Así, una vez definidos los valores límites de este parámetro, interpola linealmente en todos
los valores de temperatura, Ti, para cada valor de NDVI.
Jiang-Islam han encontrado que el valor de φ para cada clase de NDVI es muy cercano a
1.26.

Ecuación de Budyko
Esta ecuación fue obtenida por Budyko (1974) y está basada en un balance de masas, a
partir de estudios realizados por Schreiber en 1904 y por Ule y Penck. Budyko observó que
cuando la precipitación decrece, la escorrentía superficial también decrece, pero cuando la
precipitación aumenta, el caudal tiende a algún valor límite. Schreiber propone una
ecuación de interpolación para la escorrentía anual en cuencas de Europa central, dada por

⎛ a⎞
R = P exp ⎜ − ⎟ (45)
⎝ P⎠

donde a es una constante para una cuenca dada. Ahora, una forma para la
evapotranspiración anual es
⎡ ⎛ a ⎞⎤
ETR = P ⎢1 − exp ⎜ − ⎟⎥ (46)
⎣ ⎝ P ⎠⎦

En 1911, Dekop propone la ecuación

⎛ P ⎞
ETR = ETP tanh ⎜ ⎟ (47)
⎝ ETP ⎠

Budyko (1948 y 1974) postula que las siguientes condiciones son validas, en el límite de
condiciones muy secas,

R ETR ETP
→0 ó →1 cuando →∞
P P P
y en el límite de condiciones muy húmedas,
ETP
ETR → ETP cuando →0
P

Budyko interpola las ecuaciones (46) y (47), proponiendo la media geométrica.

 14

1
⎡ ⎛ P ⎞⎛ ⎛ ETP ⎞ ⎛ ETP ⎞ ⎞⎤ 2
ETR = ⎢ ETP ⋅ P ⋅ tanh ⎜ ⎟ ⎜⎜ 1 − cosh ⎜ ⎟ + senh ⎜ ⎟ ⎟⎟⎥ (48)
⎣ ⎝ ETP ⎠ ⎝ ⎝ P ⎠ ⎝ P ⎠ ⎠⎦

donde ETR es la evaporación real, P es la precipitación, ETP es la evaporación potencial

en mm por el período considerado (Budyko, 1974 y Brutsaert, 1982).
En la ecuaciòn (48) se observa un comportamiento casi igual al de las ecuaciones
propuestas por Turc y Budyko para el cálculo de la evaporación real, ya que estas tuvieron
un desarrollo teórico semejante. Ambas ecuaciones están expresadas en función de la
precipitación y la evaporación potencial, donde la evaporación potencial en Turc es el
parámetro L, el cual está en función de la temperatura.
1.2

0.8
ETR/P

0.6

0.4

TURC
0.2 BUDYKO

0
0 1 2 3 4 5 6

ETP/P

Figura 3 Relación adimensional entre la evaporación real y potencial para la

ecuación de Turc y Budyko.

Ecuación del Factor Regional

Para esta metodología el estimativo de la evaporación real anual se realiza a través de una
ecuación empírica de la forma
P
E= 1
(49)
⎧⎪ ⎛ P ⎞ α ⎫⎪ α
⎨1 + ⎜⎜ ⎟⎟ ⎬
⎪⎩ ⎝ R n ⎠ ⎪⎭

Donde P es la precipitación anual (mm), Rn es el equivalente de agua de la radiación neta en

(mm) y α es un parámetro ajustable que depende de la escala de estimación de la variable
(Choudhury, 1999). Para un experimento realizado en 10 grandes cuencas del mundo con
áreas del orden de 106 km2 este autor encontró que α = 1.8 con un error medio absoluto de
36 mm/año. El efecto de la variabilidad espacial de P y Rn en la determinación de la
evaporación desde la ecuación empírica explica las posibles diferencias del valor de alfa
dependiendo de la escala regional de trabajo. El método asume que la evaporación anual no
 15

excede los valores correspondientes a P ó Rn y dispone, para las cuencas de trabajo, de la

relación µ = P/Rn en el rango 0.23 < µ < 1.97, típico de regiones poco húmedas con altos
valores de radiación neta. Alrededor de esta metodología, se propone el cálculo del factor
alfa para la región comprendida entre 5 y 9° de laltitud Norte y 73°30’ y 77°30’ W que
abarca a Antioquia y zonas de departamentos vecinos, que minimiza el error medio
absoluto entre las evaporaciones estimadas y las obtenidas al realizar el balance de largo
plazo en 52 cuencas con áreas entre 25 y 5300 km2 con condiciones climáticas diferentes.
Se encuentra una valor de α = 1.91 con mejores resultados en el rango de 0.85 < µ < 1.9,
pero aplicable a la zona de trabajo con 0.85 < µ < 6.37 con valores altos de µ típicos de
regiones tropicales húmedas.

Método de Kazama-Sahoo
Kazama y Sahoo (1999) estudiaron la relación entre la evapotranspiración real y el índice
normalizado de vegetación (NDVI). Para ello estudiaron el comportamiento de ambas
variables en cinco cuencas con similitudes hidrológicas y climáticas pero con diferencias
geológicas y topográficas, dos de ellas ubicadas en Vietnam y tres ubicadas en Tailandia,
para finalmente aplicar los resultados obtenidos en la cuenca del río Mekong. Este estudio
se efectuó para un mes típicamente seco y un mes típicamente lluvioso, de manera que se
pudieran encontrar las diferencias estacionales existentes. La evapotranspiración fue
obtenida a partir del balance hidrológico de largo plazo para cada una de las cuencas. Estos
autores encontraron que para su región de estudio, la evapotranspiración y el NDVI siguen
una relación logarítmica de la forma

ETR = 141.9 ln( NDVI ) + 337.7 (50)

Referencias

Baldion, J. V., La Evapotranspiración Potencial, Instituto Colombiano de Hidrología,

Meteorología y Adecuación de Tierras - HIMAT, Bogotá, D.E, 1986.

Barco, O.J, L.A. Cuartas., Estimación de la evaporación en Colombia (Tesis de Pregrado,

Ingeniería Civil), Medellín, Universidad Nacional de Colombia, 1998.

Brutsaert, W. Evaporation into the atmosphere, theory, history and applications. Klower
Academic Publishers, Boston, 1982.

Cassiraga, F., Incorporación de información blanda para la cuantificación de la

incertidumbre: aplicación a la hidrogeología. [Tesis de Doctorado] Universidad
Politécnica de Valencia, Escuela Superior de Ingenieros de Caminos, Canales y
Puertos, 1999.

Chaves, B., A. Jaramillo, Regionalización de la distribución del brillo solar en Colombia

por métodos de conglomeración estadística. Atmósfera 25, 1996.
 16

Choudhury, B., Evaluation of an empirical equation for annual evaporation using field
observations and results from a biophysical model. Journal of Hydrology, 216, 99-
110, 1999.

Chow, V. T., D.R. Maidment, L.W. Mays. Hidrología Aplicada. Mc Graw Hill, 1994.

Eslava R. J. A. Régimen de la Presión Atmosférica en Colombia. Academia Colombiana

de Ciencias Exactas, Físicas y Naturales. Colección Jorge Alvarez Lleras. N°8.
Santa Fe de Bogotá, D.C. 1995.

Frére, M., Rijks, J. Q., Rea, S., Evaporación, En: Estudio Agroclimatológico de la zona
Andina, Nota técnica 161, OMM – N° 506, FAO-UNESCO-OMM, 81-297, 1978.

García, B.J., D.J López, Fórmula para el cálculo de evapotranspiración potencial adaptada
al trópico. Agronomía Tropical. Venezuela, 1970.

Holdrige, L.R. Life Zone Ecology. Tropical Science Center, San José de Costa Rica, 1978.

Ianko’s Gis Page: http://www.ian-ko.com/triangulated_irregular_network.htm [citado 3 de

Noviembre de 2001].

Jaramillo R. A, Relación entre la Evaporación y los elementos Climáticos. (nota técnica).

Cenicafé , vol. 40 N° 3, 1989.

Jiang, L., I. Shafiqul, A methodology for estimation of surface evapotranspiration over

large areas using remote sensing observations.Geophysical Research Letters, 26, 17,
2773-2776, 1999.

Kazama, S., G.B. Sahoo, Estimation of evapotranspiration in the Mekong river basin.
INCEDE Report, 19, 103-111, 2000.

Linsley, R.K., M.A Kohler, J.L Paulhus, Hidrología para Ingenieros. Mc Graw-Hill.
Segunda edición. 1986.

Lippmann, P., An introduction to computing with neural nets. IEEE Acoustics, Speech, and
Signal Processing Magazine, 1-22, 1987.

Morton, F.I, Operational estimates of areal evapotranspiration and their Significance to the
science and practice of hydrology. In : Journal of Hydrology. National Hydrology
Research Institute, Enviromental Canada, Ottawa, Agosto, 1982.

Poveda, G. Y O.J Mesa, Efectos Hidrológicos de la Deforestación. ENERGÉTICA.

Facultad de Minas Universidad Nacional de Colombia. Sede Medellín. N°16, pág.
91 – 102. Diciembre, 1995.

Remenieras, G. Tratado de Hidrología Aplicada. Técnicos Asociados S.A. Segunda

edición. Barcelona. 1974.
 17

Rodríguez, H., F. González. Manual de radiación solar en Colombia. H. Rodríguez y F.

González Editores, Bogotá, 1992.

Stewart, J.B., J.W. Finch, Application of remote sensing to forest hydrology. Journal of
Hydrology, 150, 701-716, 1993.

Vélez, J. I., G. Poveda, O. Mesa, Balances hidrólogicos de Colombia, Universidad

Nacional de Colombia, Sede Medellín, Primera Edición, 2000.