Escuela de Postgrado Universidad Nacional Agraria La Molina

Maestría en Recursos Hídricos

Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde

CLASE IV
EVAPORACIÓN

1. Introducción

La evaporación es un proceso esencial del ciclo hidrológico, pues se estima

que el 75% aproximadamente de la precipitación total anual que ocurre sobre
los continentes retorna a la atmósfera en forma de vapor, directamente por
evaporación o a través de las plantas por transpiración. Sin embargo, la forma
más importante de evaporación es probablemente aquella que se produce
desde el mar, ya que es la fuente primordial de toda el agua que precipita en
los continentes.

La cantidad de agua evaporada constituye una pérdida directa desde los

almacenamientos superficiales (embalses) y los subterráneos (acuíferos), de
manera que es necesaria su estimación en todos los análisis hidrológicos
relativos a la distribución del agua y en los estudios de balance y operación de
embalses.

2. Evaporación

La evaporación tiene dos fases. La primera fase ocurre cuando las moléculas
de agua adquieren suficiente energía a través de la superficie de agua y
escapan a la atmósfera. Esta energía la provee la radiación solar mediante el
calor latente de vaporización principalmente. La segunda fase consiste en el
transporte de las moléculas de vapor de agua desde la superficie de agua
hacia la atmósfera. Esto es controlado por la diferencia entre la presión de
vapor del cuerpo de agua y por la velocidad del viento debido a convección.

Existe por lo menos existen ocho factores de los que depende la tasa de
evaporación:
 Radiación solar
 Temperatura del aire
 Presión atmosférica
 Humedad relativa
 Temperatura del agua
 Velocidad del viento
 Calidad del agua
 Geometría de la superficie evaporante

La evaporación desde la superficie terrestre comprende la evaporación directa

desde la superficie del suelo y desde la superficie de la vegetación. La
transpiración a través de las hojas de las plantas consiste en extraer agua
desde las raíces de éstas, transportada hacia arriba a lo largo de sus tallos y
difundirla a la atmósfera a través de pequeñas aberturas en las hojas llamadas
estomas. Los procesos de evaporación desde la superficie terrestre y de
transpiración de la vegetación se conocen con el nombre de
evapotranspiración. La evapotranspiración depende de los dos factores
descritos anteriormente para la evaporación de una superficie abierta de agua,

1 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde
y de un tercer factor, dado por el suministro de humedad hacia la superficie de
evaporación. La evapotranspiración potencial es la evapotranspiración que
puede ocurrir desde una superficie bien cubierta por vegetación cuando el
suministro de humedad es ilimitado, y se calcula de una forma similar a la que
se aplica para la evaporación sobre una superficie abierta de agua. La
evapotranspiración real cae por debajo de su nivel potencial a medida que el
suelo se seca.

2.1 Definiciones

a. Calor Latente de vaporización del agua (Iv)

Las moléculas del agua en estado líquido se encuentran sujetas entre
ellas debido a las fuerzas de atracción intermoleculares. En el vapor de
agua las moléculas se encuentran 10 veces más separadas que en el
estado líquido, debido a que las fuerzas intermoleculares son mucho
más pequeñas. Durante la evaporación, la separación entre las
moléculas se incrementa, el trabajo producido es realizado en contra las
fuerzas de atracción intermolecular y la energía es absorbida. La energía
requerida se llama calor latente de vaporización del agua (Iv), y
disminuye ligeramente su valor con el incremento de la temperatura del
agua.
Si ‘T’ es la temperatura del agua en grados centígrados, el calor latente
de vaporización se puede calcular mediante la siguiente ecuación:

lv  2.501x106  2370T , en (J/kg)

Esto significa que 2.5 millones de joules se requieren para evaporar un
kilogramo de agua, donde joules es la unidad estándar de energía en el
Sistema Internacional (SI).

b. Calor Sensible (Hs)

Existe una porción de energía proveniente de la radiación entrante a la
superficie de la tierra que no es utilizada por la evaporación; sin embargo
ésta calienta la atmósfera en contacto con el suelo y por tanto el aire
caliente se desplaza hacia arriba. Ello produce un flujo de energía
asociado a un flujo de calor 'sensible', debido a que éste cambia la
temperatura del aire, la cual puede ser medida.

El cambio de temperatura es proporcional al producto C p  a , donde a 
es la densidad del aire y Cp es el calor específico del aire a presión
constante, tomado como: 1.01 kJ kg-1 K-1.
La densidad del aire húmedo puede ser estimada desde la siguiente
ecuación.
P
 a  3.486 ; kg m-3
275  T
Donde P es la presión atmosférica en kPa y T la temperatura del aire en
grados Centigrados.

El flujo de calor sensible o flujo de calor por unidad de área (Hs),

comúnmente asciende desde el suelo durante el día y desciende durante

2 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde
la noche, soportando la pérdida de energía de radiación desde la
superficie de la tierra. Localmente el flujo de calor sensible puede
descender durante el día si la vegetación está húmeda y la evaporación
es rápida.

c. Radiación neta
La radiación neta1 , es la diferencia entre la radiación entrante y saliente
de longitudes de onda cortas y largas. Es el equilibrio entre la energía
absorbida, reflejada y emitida por la superficie terrestre o la diferencia de
la radiación neta de onda corta entrante (Rns) y la radiación neta de
onda larga saliente (Rnl). Rn es normalmente positiva durante el día y
negativa durante la noche. El valor diario total para Rn es casi siempre
positivo para las 24 horas, excepto en condiciones extremas de latitudes
elevadas.
Hay dos métodos para estimar la radiación neta. El primero consiste en
medir la radiación neta de onda corta y de onda larga. Para ello se utiliza
un instrumento denominado radiómetro, que dispone de cuatro sensores
que miden la radiación de onda larga incidente, la de onda larga
reflejada, la de onda corta incidente y la de onda larga reflejada.
El segundo método recomendado por la FAO (1989), es a partir de los
valores estimados de la radiación neta de onda corta entrante Rns y la
radiación neta de onda larga Rnl.
Asimismo, actualmente la radiación neta (Rn) puede ser estimada
mediante técnicas de teledetección, en conjunto con la implementación
de un modelo denominado SEBAL. La metodología del algoritmo SEBAL
consiste en determinar los valores de albedo, los índices de vegetación
(IVDN, IAF), temperatura de superficie, radiación de onda corta incidente
y radiación de onda larga emitida. Ello se puede hacer mediante el
procesamiento digital de imágenes satelitales como la Landsat 5 TM.

1
Richard G. Allen, en el libro publicado por la FAO “Crop evapotranspiration - Guidelines for
computing crop water requirements - FAO Irrigation and drainage paper 56” (1998).

3 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde
2.2 Métodos conceptuales para calcular la evaporación

2.2.1 Método del Balance de Energía

Para desarrollar las ecuaciones de continuidad y energía aplicables a la

evaporación desde un tanque de evaporación como el que se muestra en la
siguiente figura, se supone un tanque circular que contiene agua, en el cual la
tasa de evaporación se calcula midiendo la tasa de disminución de la altura de
la superficie del agua.
Tasa de flujo de
vapor
Radiación neta
Calor sensible al aire
Hs Rn Mv= ρw. A.E

ρa

ρw h
E=-dh/dt
Area A

Calor conducido
hacia el suelo
G
Considerando un área unitaria en la superficie de agua, la fuente de
energía calórica es el campo de flujo neto de radiación Rn, medida en watts por
metro cuadrado; el agua suministra un campo de flujo de calor sensible Hs a la
corriente de aire y un campo de flujo de calor de suelo G a la superficie del
suelo, luego el cambio de calor en el tiempo dH/dt = Rn -Hs – G.

Si se supone que la temperatura del agua dentro del volumen de control es

constante en el tiempo, el único cambio en el calor que se almacena dentro del
volumen de control es el que se produce en la energía interna del agua
evaporada, el cual es igual a lv.Mv, donde lv es el calor latente de
vaporización y Mv la tasa de flujo de vapor. Luego se puede escribir la
siguiente ecuación:

Rn – Hs – G = lv . Mv

Sustituyendo para Mv y considerando A = 1 m2 , se puede resolver para Er :

Er 
1
Rn  H s  G 
lv  w

4 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde
Si suponemos que el campo de flujo de calor sensible Hs y el campo de flujo de
calor de suelo G equivalen a cero, entonces la tasa de evaporación Er puede
calcularse como la tasa a la cual toda la radiación neta de entrada se absorbe
por la evaporación:

Er 
1
R n 
lv  w

Ejemplo 1

Utilizando el método de balance de energía, calcule la tasa de evaporación

desde una superficie abierta de agua, si la radiación neta es 200 W/m 2 y la
temperatura del aire es 25°C, suponiendo que no existen campos de flujo de
calor sensible o de calor de suelo.

A partir de la siguiente ecuación: lv  2.501x10  2370T , en (J/kg)

6

El calor latente de vaporización a 25°C es lv = 2 441 kJ/kg.

De acuerdo a la ecuación de Thiesen-Scheel-Diesselhorst para la densidad del

agua:


 w  10001 
T  288.9414 T  3.98632 
 508929.2(T  68.12963) 

La densidad del agua será ρw = 997 kg/m3, así que sustituyendo en la ecuación
anterior, resulta:

200
Er 
2,441x10 3 x997
 8.22 x10 8 m / s
 8.22 x10 8 x1,000 x86,400mm / día
 7.10mm / día

5 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde
2.2.2 Método Aerodinámico

Además del suministro de energía calórica, existe un segundo factor que

controla la tasa de evaporación desde una superficie abierta de agua y
depende de la habilidad para transportar el vapor fuera 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, y estos dos
procesos pueden analizarse utilizando simultáneamente las ecuaciones de
transporte de masa y de energía en el aire.

En el volumen de control que se muestra en la siguiente figura se puede

considerar un plano horizontal de área unitaria localizado a una altura z por
encima de la superficie.

Radiación neta Rn
Elevación z

Elevación z

Elevación z
Flujo del aire

Evaporación E

Velocidad del viento u Temperatura del aire T Humedad Específica

qv

0.622k 2  a eas  ea u2

Mv 
 w pln z 2 / z0 2
donde Mv es la tasa de flujo de vapor, k es la Constante de Von Karman = 0.4,
zo la altura de la superficie, eas la presión de vapor en la superficie, ea la
presión de vapor a la altura z2, u2 la velocidad del viento a la altura z2, a la
densidad del aire y p es la presión atmosférica del aire igual en zo y z2.

Recordando que Mv se define para una superficie de área unitaria, como una
tasa de evaporación equivalente Ea , que se expresa en dimensiones de (L/T),
puede establecerse haciendo Mv = w.Ea en la ecuación anterior:

E a  Be as  e a  , donde B se conoce como coeficiente de transferencia de

vapor.
0.622k 2  au2
B
p w lnz2 / z0 2

El coeficiente de transferencia de vapor ´B´, varía de un lugar a otro.

Este tipo de ecuación fue propuesto por primera vez por Dalton en 1802.

6 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde
Ejemplo 2

Calcule la tasa de evaporación de una superficie abierta de agua utilizando el

método aerodinámico con una temperatura de 25°C, una humedad relativa del
40% , una presión de aire de 101.3 kPa y una velocidad de viento de 3 m/s ,
todas medidas a una altura de 2.0 m por encima de la superficie del agua . Se
supone una altura de rugosidad de z0 = 0.03 cm.

El coeficiente de transferencia de vapor B está dado por la ecuación anterior,

usando k = 0.4, ρa = 1.19 kg/m3 para aire a 25 °C y ρw = 997 kg/m3 . Luego

0.622k 2  au2
B
p w lnz2 / z0 2
0.622 x0.42 x1.19 x3


101.3x103 x997ln 2 /(3x10 4 ) 2

= 4.54 x 10-11 m/Pa·s

La tasa de evaporación utilizando eas =3167 Pa a 25°C según tabla, ea =HR .

eas = 0.4 x 3167 = 1267 Pa. :

Ea = B(eas - ea)
= 4.54 x 10-11(3167 – 1267)
= 8.62 x 10 –8 m/s

= 7.45 mm/día

Se puede mencionar que para el cálculo de eas se puede utilizar la ecuación:

 17.27T 
eas  611exp  , en pascales y T = temperatura del aire (°C )
 237.3  T 

7 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde
2.2.3 Método combinado (balance de energía y aerodinámico)

La evaporación puede calcularse utilizando el método aerodinámico cuando el

transporte de vapor es significativo o aplicar el método de balance de energía
cuando el suministro de energía es también significativo. Pero normalmente,
estos dos factores no son limitantes, luego es conveniente una combinación de
los dos métodos. En el método del balance de energía es difícil calcular el
campo de flujo de calor sensibles Hs. Sin embargo, como el calor se transfiere
por convección a través del aire que se localiza encima de la superficie del
agua, y el vapor de agua se transfiere por convección, en forma similar puede
suponerse que el campo de flujo de calor del vapor lv.Mv y el campo de flujo de
calor sensible Hs son proporcionales, en donde la constante de
proporcionalidad se conoce como la relación de Bowen β (Bowen,1926):
Hs

lvmv
Por otra parte, la ecuación de balance de energía con campo de flujo de calor
de suelo G = 0 puede escribirse como :
Rn = lv mv(1+ β)

La relación de Bowen se calcula utilizando simultáneamente las ecuaciones

para transporte de vapor y para calor, lo cual es similar a la unión de las
ecuaciones de transporte de vapor y de momentum que se utilizan en el
desarrollo de la ecuación de Thornthwaite-Holzman.

Si consideramos como nivel de referencia la superficie de evaporación y la

corriente de aire por encima de ésta , respectivamente, puede demostrarse que
la tasa de evaporación (Er )que se calcula de la tasa de radiación neta y la tasa
de evaporación que se establece utilizando los métodos aerodinámicos (E a ) se
combinan para dar un valor estimado ponderado de evaporación E por :
 
E Er  Ea
 

donde  es la constante psicrométrica y  es el gradiente de la curva de

presión de saturación del vapor a una temperatura de aire ‘T’, los factores  /(
+) y  /(  +  ), suman la unidad.
4,098eas

237.3T 2 , en Pascal/°C.

  66.8 , en Pascal/°C. (ó 0.27 mm Hg/°F)

La ecuación anterior del método de combinación para el cálculo de evaporación

fue desarrollada inicialmente por Penman (1948).

8 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde
El método combinado es más preciso para el cálculo de la evaporación. Las
principales suposiciones del balance de energía son que prevalezca un flujo de
energía de estado permanente y que los cambios en el almacenamiento de
calor en el tiempo o en cuerpo de agua no sean significativos. Estas
suposiciones limitan la aplicación del método a intervalos de tiempo diarios o
mayores, y a situaciones que no involucran grandes capacidades de
almacenamiento de calor, como las que posee un gran lago o embalse. La
principal suposición del método aerodinámico esta asociada con la forma del
coeficiente de transferencia de vapor ‘B’. Se han propuesto muchas formas
empíricas para ‘B’, ajustadas localmente con información de viento y otras
observaciones meteorológicas.

El método combinado es apropiado para aplicarse a áreas pequeñas con

información climatológica detallada. La información requerida incluye la
radiación neta, la temperatura del aire, la humedad, la velocidad del viento, y la
presión del aire. Cuando parte de esta información no es disponible, deben
utilizarse las ecuaciones de evaporación más simples que requieren menos
variables (American Society of Civil Engineers,1973; Doorenbos y Pruitt, 1977).

En el caso de la evaporación desde áreas grandes, las consideraciones de

balance de energía dominan la tasa de evaporación. Para tales casos Priestley
y Taylor (1972), determinaron que el segundo término de la ecuación de
combinación es aproximadamente el 30% del primero, luego puede rescribirse
como la ecuación de evaporación de Priestley – Taylor :

E Er


donde  = 1.3. Otros investigadores han confirmado la validez de esta

aproximación, variando levemente el valor de  de una localidad a otra.

9 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde
Ejemplo 3

Aplique el método de combinación para calcular la tasa de evaporación desde

una superficie de agua abierta sujeta a una radiación neta de 200 W/m2, una
temperatura de 25 °C ,una humedad relativa del 40 % y una velocidad de
viento de 3 m/s, todas registradas a una altura de dos metros, y a una presión
atmosférica de 101.3 kPa.

La tasa de evaporación correspondiente a una radiación neta de 200 W/m 2 es

Er =7.10 mm/día, desde el ejemplo anterior y el método aerodinámico arroja E a
=7.45 mm/día, para la temperatura del aire, humedad y condiciones de viento
dadas.

4098eas

237.3  T 2
4098 x3167

237.3  252
=188.7 Pa/°C
 
E Er  Ea
 
= 0.738 x 7.10 + 0.262 x 7.45

= 7.2 mm/día

Fuente Principal: Ven Te Chow, Maidment, Mays 'Hidrología Aplicada'.

10 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde
2.3 Algoritmo para estimar la evaporación

Existen diversos métodos para estimar la evaporación desde el punto de

vista de la hidrología práctica hasta la investigación hidrológica rigurosa.

Los investigadores interesados en conseguir modelos suficientemente

realistas utilizan frecuentemente cuencas experimentales instrumentadas y
crean complejas descripciones del movimiento de la energía y agua en la
interfase: suelo - vegetación - atmósfera. El uso de dichos modelos
complejos requiere una gran cantidad de datos de variables meteorológicas
y muestreo frecuente, así como extensa información de los parámetros de
suelo y vegetación difícilmente disponibles aún en la literatura.

De acuerdo a las líneas de investigación convencionales existen dos

procesos de evaporación, una para superficies libres o abiertas de agua
llamadas evaporación potencial Eo y otra para algún cultivo específico
llamada evaporación potencial para el cultivo de referencia Erc.
Cabe indicar que existe una excepción al usar estos dos tipos de procesos
y se trata del procedimiento para estimar la evaporación en el caso de
vegetación forestal. Se ha demostrado que en dicho caso, la diferencia
entre la tasa de evaporación cuando la vegetación está húmeda y cuando
está seca es muy marcada. La hidrología forestal trata este caso de
manera más detallada.

2.3.1 Estimación de la energía disponible para la evaporación (A)

La energía disponible en la superficie terrestre es normalmente el control

primario de la evaporación desde superficies libres de agua así como para
la evaporación desde un campo cubierto de algún cultivo de referencia.

Generalmente para el uso de métodos o técnicas para estimar la cantidad

de evaporación, es necesaria la estimación preliminar de la energía
disponible.

Seguidamente se presenta el diagrama de flujo de los estimadores

preliminares de la técnica de cálculo de la energía disponible para el
proceso de la evaporación.

11 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde

ALGORITMO
Estimación de la energía disponible para el proceso de evaporación NO
La máxima horas de luz posible (N) puede ser calculada por:
24 Se dispone de
N Ws (Ecuación 1)
 mediciones locales de
Donde Ws es el 'ángulo de horas de atardecer' ó 'sunset hour la radiación solar (St
angle' en radianes dados por: en MJ m-2día-1) ?
Ws  arccos  tan  tan   Se dispone de registros SI Calcular (n/N);
donde N desde
Donde  es la latitud del lugar (positivo en el hemisferio norte y de horas de sol (n)? SI
la ecuación (1)
negativo en el hemisferio sur) y  es la declinación solar (en
NO
radianes) dados por:
 2 
  0.4093sen J  1.405  Estimar la evaporación desde
 365  tanque evaporímetro
Donde 'J' es el número de día Juliano. Dividir por lv  2.501  0.002361T
Su cálculo puede ser realizado utilizando calculadoras en (MJ/kg) para obtener St en Si existen valores de los coeficientes
disponibles en:http://www.nr.com/julian.html equivalentes de agua (mm/día) de angstrom , seleccionar as y
http://quasar.as.utexas.edu/Billinfo/JulianDateCalc.html Seleccionar el valor del albedo bs ; sino usar: as = 0.25; bs =
() desde la siguiente tabla 0.50
La máxima radiación solar extraterrestre (So ) puede ser Cobertura Albedo ()
calculada por: Agua libre 0.08 - Calcular la radiación extraterrestre
Foresta 0.11-0.16

So  15.392d r Ws sensen  cos  cos senWs  Cultivos tallo largo (p.e. Caña azúcar) 0.15.0.20 So (Ecuación 2)
Cereales 0.20-0.26
Cultivos talla corto 0.20-0.26 - So = So * lv (MJm-2dia-1)
(Ecuación 2) ; en mm/día Grass y pastos 0.20-0.26
Suelo desnudo 0.10(Húmedo)-0.35(Seco)
Nieve y Hielo 0.20(Viejo)-0.1(Nuevo) - Calcular la radiación solar incidente
Donde dr es la distancia relativa entre la tierra y el sol dada por: de onda corta
 2   n
d r  1  0.033 cos J St   as  bs  So
Algoritmo  365  Calcular la radiación neta de onda  N
corta. S n  St (1   ) ; en MJ m-2día-1 ; en MJ m-2 día-1
Se dispone de Obtener la presión de vapor 'e' con el dato
mediciones locales NO de la temperatura a punto de rocío
de la radiación neta Se dispone de datos mediante tabla o ecuación proporcionada
SI
(Rn en MJ m-2día-1)? de temperatura a en la Clase III Precipitación.Calcular la
punto de rocío ? emisividad neta e'  a  b e
SI  e e

Dividir Rn entre lv  2.501  0.002361T NO  ae = 0.34 y be = -0.14 (Coef.Emisividad)

en (MJ/kg) para obtener Rn en mm/día
 Utilizar la siguiente ecuación
en MJm-2dia-1 Ah   wc p qiTi  qoTo  PT p   4.19 x103 qiTi  qoTo  PTp 
Se requiere estimar la energía SI
por advección (Ah) ? Donde qi, qo : Caudales de entrada y salida por unidad de área de espejo de agua (mm), P:
Precipitación (mm), Ti, To, Tp: Temperatura en °C del caudal de entrada, salida y de
NO
precipitación. La energía disponible para la evaporación A = Rn + Ah ; mm/día.
A = Rn (mm/día)

12 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde

Calcular la emisividad neta mediante la

Se dispone de
NO siguiente ecuación en función de la
medición de la
temperatura media del aire (°C)
temperatura
mínima del aire ?

e'  0.02  0.261exp  7.77 x104 T 2 

SI

Reemplazar la temperatura al punto de

rocío con la temperatura mínima del aire
y obtener la presión de vapor 'e'
mediante tabla o ecuación proporcionada
en la Clase III Precipitación.

Calcular la emisividad neta e'
mediante e'  ae  be e

Seleccionar los coeficientes de nubosidad ac y bc,

Usar ac = 1.35 ; bc = -0.35 para áreas áridas ó ac = 1.00; bc = 0.00
para áreas húmedas. Entonces:
- Calcular la radiación solar para cielo claro (Sto ) Sto  as So
; MJm-2 día-1 St
f  ac  bc
- Calcular el factor de nubosidad (f) Sto

- Calcular la radiación neta de onda larga (Ln)
Ln   fe ' T  273.2
4
; MJm-2día-1
Donde  : Constante de Stefan-Boltzmann (4.903x10-9 MJ m-2 °K-
4dia-1); T:Temperatura media del aire (°C).

- Calcular la Radiación neta

Rn=Sn + Ln ; MJ m-2 día-1

13 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde
Ejemplos de website para el cálculo del número juliano:

http://www.geocities.com/xgarciaf/java/juliano.htm,
http://www.astrogea.org/foed/efemerides/dia_julia.htm

2.3.2 Cálculos de otras variables de entrada para estimar la evaporación

2.3.2.1 Déficit de Presión de Vapor (D)
Algunos métodos para la estimación de la evaporación requieren el valor del
déficit de la presión de vapor (D), que es la diferencia entre la presión de vapor a
saturación 'es' y la presión de vapor actual o ambiente 'e', promedio en un
periodo de tiempo determinado.
Como la presión de vapor a saturación no es una función lineal de la
temperatura, este procedimiento para calcular el déficit de presión de vapor a
partir de registros climatológicos puede sufrir un error estimado en 30%,
particularmente en zonas áridas donde el ciclo de temperatura diaria es
frecuentemente grande.
En términos generales, es posible calcular el déficit de presión de vapor,
calculando individualmente la presión de vapor para la temperatura máxima y
para la temperatura mínima y luego tomando el promedio de los valores, lo que
da una mejor estimación del déficit promedio diario que realizando el promedio
de las temperaturas (o humedades relativas) y luego calcular el déficit.
Como se explicó anteriormente en la Clase III, los valores de humedad
atmosférica son calculados utilizando:
- La Hrmáx y Hrmín en porcentaje,
- El promedio diario de las temperaturas del bulbo húmedo y seco (Ts y Th) en °C.
- El promedio diario de la temperatura del punto de rocío (td) en °C.

El tiempo de toma de la medición es importante pero frecuentemente no se

encuentra registrado.
Afortunadamente la presión de vapor actual del aire tiende a ser constante y una
medida por día puede ser suficiente.
Si no se tiene el verdadero valor del déficit de presión de vapor (D), éste podrá ser
estimado mediante datos climatológicos, siguiendo el siguiente métodos de cálculo.

a. Humedad relativa
Datos de entrada: Tmáx, Tmín, HR
 e (T )  e s (Tmín )  (100  HR )
D   s máx  ; en kPa
 2  100
 17.27T 
Donde es  0.6108 exp ; en kPa.
 237.3  T 

b. Temperaturas del bulbo húmedo y del bulbo seco

Datos de entrada: Tmáx, Tmín, Ts, Th
es (Tmáx )  es (Tmín )
D e
2
Donde e: Presión de Vapor (Kpa) (Tabla 4.4.2) para psicrómetros de aspiración y no
aspiración.

14 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde

Fuente: Maidment David, 'Handbook of Hydrology', McGraw-Hill, 1992.(Pág 4.34, 4.35)

15 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde
c. Temperatura al punto de rocío
Datos de entrada: Tmáx, Tmín, Td
es (Tmáx )  es (Tmín )
D  es (Td )
2
Finalmente, debe tenerse en cuenta la siguiente corrección por altura de medición,
debido a que la técnica descrita anteriormente es apropiada para calcular el déficit de
presión de vapor y la velocidad del viento a 2.0 m de altura desde el suelo.
U 2' 34.9648

U2  z  0.08   zu  0.08 
ln e  ln 
 0.001476   0.01476 
Donde: U'2 es la velocidad del viento corregida a 2.0 m; ze y zu son las alturas de registro
de la humedad y velocidad del viento respectivamente.

2.3.3 Estimación de la tasa de evaporación desde una superficie libre de agua

La ecuación es la siguiente:
E p  Fp1 A  Fp2 D ; en mm/día
Donde A es la energía disponible para evaporación(mm/día) dado por A=(Rn + Ah).
1
El coeficiente Fp es función de la temperatura y de la elevación del sitio.

Fp1 
 
Donde  es el gradiente de la curva de presión de saturación del vapor a una
temperatura de aire T (°C) y es (kPa) es la presión de vapor a saturación.
4098es
 -1
237.3  T 2 ; kPa °C

P
 es la constante psicrométrica = 0.0016286 ; kPa°C-1
lv
5.256
 293  0.0065Z 
P es la Presión Atmosférica = P  101.3  ; kPa, Z es la elevación en
 293 
metros.

Elevación Temperatura Vel.Viento FP1 FP2 Frc1 Frc2 Fpthúmedo Fptárido

(m) (°C) (m/s)
0 10 3 0.553 3.028 0.383 2.937 0.696 0.962
0 10 6 0.553 4.895 0.293 4.495 0.696 0.962
0 30 3 0.781 1.505 0.643 1.588 0.985 1.360
0 30 6 0.781 2.433 0.546 2.697 0.985 1.360
1000 10 3 0.582 2.832 0.411 2.803 0.733 1.012
1000 10 6 0.582 4.578 0.318 4.336 0.733 1.012
1000 30 3 0.801 1.371 0.670 1.470 1.010 1.394
1000 30 6 0.801 2.216 0.575 2.524 1.010 1.394
Fuente: Maidment, 'Handbook of Hydrology', pág. 4.37

16 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde

El coeficiente FP2 está en función de la temperatura, velocidad del viento y altitud del
lugar.
   6.431  0.536U 2 
FP2   
    Iv
Donde U2 es la velocidad del viento medida a 2.0 m.
é Iv el calor latente de vaporización = 2.501-0.002361 Ts (MJkg-1)

Por ejemplo

Para la siguiente información:

- Radiación neta (MJm-2día-1) = 17.0

- Tmáx.aire (°C) = 35
- Tmín.aire (°C) = 4
- T del agua (°C) = T.punto rocío (°C) = 5
- HR(%) = 80
- Vel.Viento a 2m. (m/s) = 3.0
- Horas de sol (hr) = 9
- Altitud (msnm) = 1000
- Latitud Sur (°) = 12
- Albedo = 0.08
- Q1(mm) = 4, Q0 (mm) = 2, Pp (mm) = 3
- T1(°C) = 10, T0(°C) = 9, Tp(°C) = 8

La evaporación calculada por el método es: 3.053 mm/día

17 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde
2.4 Métodos para estimar la evaporación a partir de información registrada en tanques
de evaporación y evaporímetros

La información del tanque de evaporación provee la mejor indicación de evaporación en

superficie de aguas abiertas.

Básicamente existen dos tipos de aparatos para medir la evaporación: los tanques de
evaporación y los evaporímetros que son superficies húmedas de porcelana o papel. Las
evaporaciones medidas en tales aparatos pueden ser consideradas como máximas o
potenciales y son indicativas del poder evaporante de la atmósfera, de manera que
aplicando coeficientes de reducción se podrán deducir los valores más probables de la
evaporación en la superficie de interés.

Los tanques de evaporación utilizados actualmente en los diversos países son de

dimensiones, formas y características diferentes.

En el Perú se utiliza principalmente el denominado Tanque clase 'A' del US Weather

Bureau, que es un tanque montado sobre un enrejado de madera de forma que su base
está a 5 o 10 cm, por encima del suelo, permitiendo la circulación del aire por debajo.
Está construido de fierro galvanizado y es llenado hasta 5.1 cm bajo su borde. Las
mediciones se realizan apoyando en un tubo un tornillo micrométrico cuya punta se
enrasa con el nivel del agua.

Tanque de evaporación Clase A

18 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde

Fuente: Meteorological Stations, Munro.

Para el caso de estimar la evaporación en presas, lagos o cualquier tipo de gran

almacenamiento, se debe tener en cuenta que la evaporación es mayor en depósitos
pequeños que en los grandes por lo que se deben corregir los datos registrados por
algún factor que varía entre 0.6 y 0.8. En general 0.7 es un buen valor del factor de
corrección.

Kohler y Parmelei (1967), teniendo en consideración que para un cuerpo de agua se

debe involucrar el almacenamiento de calor y la energía por advección del agua que
entre y sale de dicho cuerpo de agua, propusieron la siguiente ecuación que tiene en
cuenta la presión de vapor a saturación del cuerpo de agua y del tanque evaporímetro.

esL  eZ
E  K' ETanque
esp  eZ

19 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde
Donde:
esL : Presión de vapor a saturación para la máxima temperatura que resulte sobre el
espejo de agua (kPa).
esp : Presión de vapor a saturación para la máxima temperatura en el tanque de
evaporación (kPa).
eZ : Presión de vapor media del aire a una altura Z sobre el espejo de agua (kPa)
K’ : Coeficiente en función al tipo de tanque de evaporación. Para el tanque de
evaporación clase A igual a 1.5 para Z = 4 m. (Webb, 1966).

Nota: 01 kpa = 10 mb

Ejemplo:

 Temperatura media del aire a 4.0 m = 30oC

 Temperatura del lago
o Máxima = 20 oC
o Mínima = 18 oC
 Temperatura del tanque de evaporación
o Máxima = 28 oC
o Mínima = 25 oC
 Humedad Relativa = 25%
 Velocidad del viento a 4.0 m = 6m/s
 Evaporación de tanque = 9 mm.
Solución:
 Desde tablas o fórmulas:
o esL = 2.337 kpa
o esp = 3.781 kpa
o es = 4.243 kpa
25
 Presión de vapor del aire = es  (4.243)  1.061kpa
100
 2.337  1.061 
 E L  1.5 (9.0)  6.33mm
 3.781  1.061 

Los evaporímetros son aparatos que miden el poder evaporante de la atmósfera, se

utilizan esferas o placas de porcelana porosa saturada de agua, éstos dispositivos son
empleados por meteorólogos y agrónomos en investigaciones y estudios de
transpiración.

En el Perú se utiliza con frecuencia el evaporímetro o atmómetro Piche, el cual está

constituido por un tubo cilíndrico de vidrio de 25 cm de diámetro, el tubo está graduado y
cerrado por su parte superior y abierto por su parte inferior, la cual es obstruida por un
papel filtro normalizado de 30 mm de diámetro y 0.5 mm de espesor. Llenado el aparato
con agua destilada, ésta se evapora progresivamente a través del papel, el proceso de
evaporación en este instrumento está ligado al déficit hidrométrico y no toma en cuenta
la influencia de la insolación, pues el atmómetro se coloca bajo abrigo y tampoco es
afectado por la lluvia.

20 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde

Ejemplo : Evaporación desde el embalse Tinajones

21 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde

2.5 Método de balance de agua

Este es un método indirecto para calcular la evaporación. Se basa en la ecuación de

continuidad que para un gran almacenamiento es:
E = I - O - V
Donde
E : Volumen de evaporación en el intervalo de tiempo  t considerado.
I : Volumen de entrada al vaso en el  t considerado (Precipitación directa y
escurrimiento).
O : Volumen de salidas del vaso en el  t considerado (infiltración y escurrimiento,
en el caso de presas el último se forma por las salidas de la obra de toma y el
vertedero de excedencias).
V : Cambio en el volumen almacenado en el  t considerado.

La ecuación anteriormente señalada se puede usar con fines estadísticos para estimar la
evaporación que podría presentarse en un vaso de almacenamiento dado o en vasos
cercanos a él. Sin embargo, para que sus resultados sean confiables, I, O y  V deben
estar medidos con precisión; dado que E es en general un término relativamente
pequeño en comparación con los demás, errores leves en la medición de I, O ó  V
conducirán a errores considerables en la estimación de E.

Debe señalarse que existe otra aplicación de la ecuación de balance hídrico muy
utilizada para calcular de manera indirecta el volumen de entrada a un embalse, ésta es:
I = O + E +  V.

22 18/04/13
 Escuela de Postgrado Universidad Nacional Agraria La Molina
Maestría en Recursos Hídricos
Curso: Métodos de Análisis en Ingeniería de Recursos Hídricos
Dr. Eduardo A. Chávarri Velarde
Ejemplo: Balance Hídrico del Embalse Pasto Grande (PEPG -INADE).

VARIACION DE NIVELES Y VOLUMENES DEL EMBALSE "PASTO GRANDE"

AÑO DESCRIPCION ENE FEB MAR ABR MAY JUN JUL AGO SET OCT NOV DIC
Volumen (MMC) 3.20 3.30 3.40 3.45 3.10
Area (ha) 700.00 730.00 750.00 770.00 690.00
Altitud (m.s.n.m.) 4,517.40 4,517.43 4,517.44 4,517.45 4,517.42
1989
Evaporación (m³/s) 0.13 0.32 0.38 0.42 0.38
Variac.Volumen (m³/s) 0.039 0.384 0.019 (0.131)
Ingreso (m³/s) 0.362 0.763 0.442 0.254
Volumen (MMC) 6.00 7.80 10.10 10.20 12.00 15.40 16.00 16.49 15.70 15.00 15.15 22.70
Area (ha) 1,200.00 1,400.00 1,580.00 1,610.00 1,620.00 1,700.00 1,730.00 1,740.00 1,724.00 1,705.00 1,720.00 1,859.00
Altitud (m.s.n.m.) 4,517.75 4,517.87 4,518.08 4,518.14 4,518.19 4,518.41 4,518.47 4,518.49 4,518.45 4,518.41 4,518.44 4,518.82
1990
Evaporación (m³/s) 0.28 0.50 0.45 0.64 0.56 0.30 0.60 0.64 0.78 0.81 0.95 0.94
Variac.Volumen (m³/s) 1.080 0.740 0.860 0.039 0.300 1.310 0.220 0.295 (0.305) (0.261) 0.193 2.680
Ingreso (m³/s) 1.355 1.239 1.312 0.684 0.862 1.613 0.825 0.932 0.478 0.550 1.145 3.623
Volumen (MMC) 51.40 66.60 97.00 102.40 102.50 104.00 104.00 102.60 101.30 98.60 96.30 95.70
Area (ha) 2,665.00 2,920.00 3,480.00 3,575.00 3,580.00 3,610.00 3,610.00 3,585.00 3,558.00 3,504.00 3,465.00 3,457.00
Altitud (m.s.n.m.) 4,520.09 4,520.65 4,521.59 4,521.73 4,521.74 4,521.78 4,521.78 4,521.75 4,521.71 4,521.63 4,521.57 4,521.55
1991
Evaporación (m³/s) (0.46) 0.37 (0.49) 0.66 0.99 0.86 1.31 1.17 1.74 1.44 1.06 1.05
Variac.Volumen (m³/s) 10.715 6.280 11.350 2.083 0.037 0.579 - (0.523) (0.502) (1.008) (0.887) (0.224)
Ingreso (m³/s) 10.254 6.652 10.858 2.746 1.031 1.441 1.314 0.650 1.235 0.432 0.177 0.826
Volumen (MMC) 106.80 106.80 103.80 100.70 99.00 98.50 97.60 95.80 93.40 90.70 90.40 89.60
Area (ha) 3,660.00 3,660.00 3,605.00 3,555.00 3,515.00 3,500.00 3,482.00 3,458.00 3,410.00 3,355.00 3,345.00 3,330.00
Altitud (m.s.n.m.) 4,521.86 4,521.86 4,521.77 4,521.70 4,521.64 4,521.62 4,521.60 4,521.56 4,521.48 4,521.41 4,521.39 4,521.37
1992
Evaporación (m³/s) 0.18 1.06 1.96 1.39 1.20 0.84 1.23 1.27 1.55 1.60 1.86 1.76
Variac.Volumen (m³/s) 4.144 - (1.120) (1.200) (0.630) (0.193) (0.336) (0.672) (0.926) (1.008) (0.116) (0.299)
Ingreso (m³/s) 4.323 1.056 0.839 0.189 0.569 0.647 0.895 0.603 0.626 0.592 1.747 1.460
Volumen (MMC) 121.50 124.50 136.70 136.70 135.70 135.30 133.60 134.90 131.80 129.40 128.80 136.70
Area (ha) 3,650.00 3,705.00 3,940.00 3,940.00 3,915.00 3,905.00 3,880.00 3,898.00 3,845.00 3,795.00 3,787.00 3,940.00
Altitud (m.s.n.m.) 4,526.40 4,526.48 4,526.81 4,526.81 4,526.78 4,526.76 4,526.73 4,526.75 4,526.67 4,526.61 4,526.59 4,526.81
1,993
Evaporación (m³/s) (0.47) 1.41 0.16 1.04 1.23 1.13 1.34 0.83 1.58 1.23 1.58 0.60
Variac.Volumen (m³/s) 9.894 1.240 4.555 - (0.373) (0.154) (0.635) 0.485 (1.196) (0.896) (0.231) 2.949
Ingreso (m³/s) 9.423 2.654 4.715 1.042 0.855 0.971 0.704 1.312 0.387 0.334 1.346 3.551
Volumen (MMC) 153.80 188.30 178.00 172.10 169.00 168.30 168.30 165.60 164.50 163.00 156.10 158.60
Area (ha) 4,215.00 4,650.00 4,543.00 4,485.00 4,445.00 4,440.00 4,440.00 4,405.00 4,385.00 4,365.00 4,265.00 4,280.00
Altitud (m.s.n.m.) 4,527.24 4,527.97 4,527.78 4,527.69 4,527.62 4,527.61 4,527.61 4,527.55 4,527.52 4,527.49 4,527.33 4,527.35
1,994 Evaporación (m³/s) (0.01) (0.31) 0.94 0.60 1.53 1.50 1.58 1.77 2.54 2.55 2.96 0.57
Variac.Volumen (m³/s) 5.862 14.137 (3.696) (1.775) (1.157) (0.270) - (1.008) (0.424) (0.560) (2.662) 0.299
Ingreso (m³/s) 5.849 22.029 4.416 2.336 0.371 1.232 1.585 0.759 2.112 1.988 0.295 0.872
Rebose (m³/s) 8.198 7.177 3.509
Volumen (MMC) 162.00 162.00 180.00 171.00 165.80 163.00 160.00 155.40 150.60 145.80 138.30 136.60
Area (ha) 4,250.00 4,295.00 4,540.00 4,470.00 4,410.00 4,365.00 4,330.00 4,255.00 4,185.00 4,100.00 3,970.00 3,860.00
Altitud (m.s.n.m.) 4,527.30 4,527.37 4,527.77 4,527.67 4,527.56 4,527.49 4,527.43 4,527.31 4,527.19 4,527.06 4,526.85 4,526.71
1,995 Evaporación (m³/s) 1.05 1.68 (0.34) 1.40 1.54 1.40 1.42 1.95 1.78 2.17 2.03 1.35
Variac.Volumen (m³/s) (0.709) 1.075 6.869 (1.929) (1.941) (1.045) (1.120) (1.717) (1.852) (1.792) (2.894) (1.904)
Ingreso (m³/s) 0.344 2.751 6.533 0.564 0.727 1.354 1.283 1.164 0.821 1.365 0.712 0.737
Salida (m³/s) 1.095 1.132 0.997 0.984 0.927 0.895 0.984 1.580 1.290
Volumen (MMC) 142.60 157.40 163.80 167.50 164.50 160.60 158.10 154.60 148.80 140.40 133.30 132.70
Area (ha) 4,040.00 4,285.00 4,380.00 4,430.00 4,385.00 4,335.00 4,300.00 4,248.00 4,155.00 4,005.00 3,875.00 3,865.00
Altitud (m.s.n.m.) 4,526.96 4,527.37 4,527.52 4,527.59 4,527.52 4,527.44 4,527.38 4,527.30 4,527.14 4,526.90 4,526.72 4,526.70
1,996 Evaporación (m³/s) (0.28) 0.06 1.14 0.89 1.36 1.31 1.29 1.31 1.70 2.16 1.24 0.29
Variac.Volumen (m³/s) 3.500 5.900 2.389 1.427 (1.120) (1.504) (0.746) (1.306) (2.237) 3.136 (2.739) (0.224)
Ingreso (m³/s) 4.626 6.991 4.257 3.107 1.005 0.373 1.159 0.713 0.928 6.812 0.334 1.960
Salida (m³/s) 1.404 1.035 0.728 0.790 0.767 0.567 0.611 0.712 1.464 1.518 1.831 1.898

Fuente Bibliográfica:
- Maidment David 'Handbook of Hydrology', McGraw Hill, 1992.
- Ven Te Chow 'Hidrología Aplicada', McGraw Hill, 1994.
- Gupta Ram ‘Hydrology and Hydraulic Systems’, 2008.

Trabajo encargado N°03

1. Correlacione los cálculos realizados de Er, Ea y E con la evaporación calculada

mediante el método de balance de agua utilizando datos reales de algún embalse.
2. Correlacione los cálculos de evaporación realizados mediante la aplicación del
algoritmo expuesto en el punto 2.3 con datos reales de algún embalse.

23 18/04/13