ETP Metodos
ETP Metodos
ETP Metodos
LISÍ
LISÍMETROS
El cálculo empírico de la Evapotranspiración potencial se realiza mediante lisímetros. El
problema es que no siempre se dispone de estos aparatos cuando se van a estimar los
requerimientos hídricos de los cultivos, por lo que es necesario recurrir a métodos indirectos
basados en modelos teóricos.
MÉTODOS TEÓ
TEÓRICOS
La necesidad de desarrollar procedimientos prácticos estandarizados ha llevado a la
formulación de diversos métodos teóricos para calcular la Evapotranspiración potencial a partir
de las variables climáticas y/o algunas características de la vegetación.
BASADOS EN LA RADIACIÓ
RADIACIÓN
BASADOS EN LA TEMPERATURA
BASADOS EN LA EVAPORACIÓ
EVAPORACIÓN
COMBINADOS
1
CÁLCULO DE LA EVAPOTRANSPIRACIÓ
EVAPOTRANSPIRACIÓN POTENCIAL MEDIANTE
Cuerpo cilí
LISÍ
LISÍMETROS Colchó
cilíndrico Colchón hidrá
hidráulico
57cm
CULTIVO ESCALA
CILINDRO
CÁMARA
TUBO DEL
SUELO NEUMÁTICA
MANÓMETRO
Capacidad: 200 l 70cm
88cm MANÓMETRO
GRAVA VÁLVULA
DRENAJE
ASFALTO CUBIERTA
COLCHON TABLONES DE CAUCHO
DE AGUA DE MADERA
Un lisímetro hidráulico consta de un cuerpo cilíndrico relleno de suelo en el que se cultiva pasto.
En su base tiene un sistema de drenaje y reposa sobre un colchón hidráulico conectado a un
manómetro hidráulico que nos permite controlar la presión del líquido que hay en el interior del
colchón. Conocida la cantidad de agua que recibe el cultivo y la cantidad de agua que es drenada,
se puede calcular la cantidad de agua que se evapora a partir de los cambios de presión detectados
por el manómetro.
2
CÁLCULO DE LA EVAPOTRANSPIRACIÓ
EVAPOTRANSPIRACIÓN POTENCIAL MEDIANTE
Cuerpo cilí
LISÍ
LISÍMETROS Colchó
cilíndrico Colchón hidrá
hidráulico
57cm
CULTIVO ESCALA
CILINDRO
CÁMARA
TUBO DEL
SUELO NEUMÁTICA
MANÓMETRO
Capacidad: 200 l 70cm
88cm MANÓMETRO
GRAVA VÁLVULA
DRENAJE
ASFALTO CUBIERTA
COLCHON TABLONES DE CAUCHO
DE AGUA DE MADERA
CÁLCULO DE LA EVAPOTRANSPIRACIÓ
EVAPOTRANSPIRACIÓN POTENCIAL MEDIANTE
Cuerpo cilí
LISÍ
LISÍMETROS Colchó
cilíndrico Colchón hidrá
hidráulico
57cm
CULTIVO ESCALA
CILINDRO
CÁMARA
TUBO DEL
SUELO NEUMÁTICA
MANÓMETRO
Capacidad: 200 l 70cm
88cm MANÓMETRO
GRAVA VÁLVULA
DRENAJE
ASFALTO CUBIERTA
COLCHON TABLONES DE CAUCHO
DE AGUA DE MADERA
Los litros de vapor de agua evapotranspirados (L) durante un periodo de 24h se calculan como la diferencia
entre la lámina consumida entre un riego y otro (Lc) y las pérdidas por drenaje (D):
L = Lc - D
Estimación de la evapotranspiración potencial (ETP): se calcula como el conciente entre los litros de vapor
de agua evapotranspirados y la superficie (A) del lisímetro (2.550 cm2):
ETP = L · 10.000 / A
3
CÁLCULO DE LA EVAPOTRANSPIRACIÓN POTENCIAL
MEDIANTE MÉTODOS TEÓRICOS
Métodos basados en la radiació
radiación
Fueron desarrollados originariamente por Makkink para cultivos en condiciones húmedas en
Holanda.
Su validez puede extenderse a un amplio rango de condiciones climáticas introduciendo un
coeficiente de corrección en función de las condiciones locales de viento y humedad.
4
COMPARACIÓN DE DIVERSOS MÉTODOS TEÓRICOS CON LOS
RESULTADOS OBTENIDOS MEDIANTE MÉTODOS EMPÍRICOS
La FAO (1990) encargó a un panel de expertos un estudio para calibrar diversos métodos de
cálculo teórico de la evapotranspiración potencial.
La calibración se realizó calculando la bondad de ajuste entre los resultados empíricos obtenidos
mediante lisímetros y los obtenidos mediante fórmulas teóricas basadas en datos climáticos y en
datos sobre el estado de la vegetación.
Dado que algunos métodos sobreestiman la evapotranspiración potencial en función del grado de
humedad atmosférica, los resultados se han diferenciado dependiendo de que el clima sea
húmedo o seco.
1
Sobre o infraestimación (%) calculada a partir de los datos de 11 lisimetros situados en cada localización
2
Error estandard de las estimaciones en mm/día
5
COMPARACIÓN DE DIVERSOS MÉTODOS TEÓRICOS CON LOS
RESULTADOS OBTENIDOS MEDIANTE MÉTODOS EMPÍRICOS
MÉTODOS COMBINADOS
Penman-Monteith 4% 0.32 -1% 0.49
FAO-24 Penman (c=1) 29% 0.93 12% 0.69
FAO-24 Penman (corrected) 35% 1.14 18% 1.1
FAO -PPP-17 Penman 16% 0.67 6% 0.68
Penman (1963) 14% 0.6 -2% 0.7
Penman 1963, VPD #3 20% 0.69 6% 0.67
1972 Kimberley Penman 18% 0.71 6% 0.73
1982 Kimberley Penman 10% 0.69 3% 0.54
Businger-van Bavel 32% 1.03 11% 1.12
1
Sobre o infraestimación (%) calculada a partir de los datos de 11 lisimetros situados en cada localización
2
Error estandard de las estimaciones en mm/día
El panel de expertos consultados por la FAO recomendó la adopción del método de Penman-
Monteith como estándar para calcular la evapotranspiración potencial y especificó los
procedimientos para calcular los diversos parámetros climáticos y fisiológicos que intervienen en
su cálculo.
( Rn - G) + c p ( ea - ed ) / r a
ET o =
+ (1+ r s / r a )
6
CÁLCULO DE LA EVAPOTRANSPIRACIÓN POTENCIAL
POR EL MÉTODO DE PENMAN-MONTEITH (propuesta inicial)
NUMERADOR: (kJ m-2 s-1 kPa ºC-1)
Flujo de calor latente de la atmósfera
( Rn - G) + c p ( ea - ed ) / r a libre y de la planta por unidad de
ET o = superficie, teniendo en cuenta la
+ (1+ r s / r a ) resistencia aerodinámica de la planta.
Todo ello se pondera por la intensidad de
la presión de vapor de agua en función
λETo : Flujo de calor latente de evaporación [kJ m-2 s-1]
de la temperatura.
Rn : Flujo de radiación neta en superficie [kJ m-2 s-1]
G : Flujo de calor del suelo [kJ m-2 s-1]
ρ : Densidad de la atmósfera [kg m-3] DENOMINADOR: (kPa ºC-1)
cp : Calor específico del aire húmedo [kJ kg-1 ºC-1] Intensidad de la presión de vapor de
(ea-ed): Déficit de presión de vapor [kPa] agua en función de la temperatura, tanto
en la atmósfera libre como en la planta,
rs : Resistencia superficial [s m-1]
teniendo en cuenta la superficie de suelo
ra : Resistencia aerodinámica [s m-1] recubierta de vegetación y su resistencia
Δ : Pendiente de la curva de presión de vapor [kPa ºC-1] aerodinámica.
γ : Constante psichrométrica [kPa ºC-1]
λ : Calor latente de evaporación [MJ kg-1]
7
EL BALANCE HÍDRICO DE TENERIFE
INFLITRACIÓN ESCORRENTÍA
120
Precipitación (mm)
100
80
60
40
20
0
0 500 1000 1500 2000
Altura (m)
8
VARIACIÓN ALTITUDINAL DE LA PRECIPITACIÓN EN
TENERIFE
El relieve ejerce una enorme influencia sobre la distribución espacial de la precipitación. Los
valles y laderas expuestos a los vientos del sector Norte son relativamente más lluviosos que los
expuestos al Sur. Las zonas más lluviosas se encuentran en la fachada Norte entre los 750 y los
1500 metros.
120
Precipitación (mm)
100
80
60
40
20
0
0 500 1000 1500 2000
Altura (m)
… en el Monte
Las precipitaciones de la Esperanza
más altas de enero (130-140 mm)
se registran en el
Macizo de Teno
(160-170 mm) y …
9
VARIACIÓN ESPACIAL DE LA ESCORRENTÍA EN TENERIFE
Macizo
de Anaga
Macizo de
Teno
Altos de
Vilaflor
Los mayores porcentajes de escorrentía se producen en las zonas recubiertas por emisiones
volcánicas antiguas debido a que los suelos poseen una permeabilidad moderada o baja.
Por el contrario, los suelos de las zonas cubiertas por emisiones volcánicas recientes presentan
una elevada permeabilidad y el porcentaje de escorrentía es bajo.
Las zonas donde se produce la mayor recarga de lluvia coincide con las zonas de mayor
pluviosidad que están recubiertas por materiales volcánicos recientes.
El agua de recarga se refiere al volumen de agua procedente de las precipitaciones que se
infiltra en el terreno y ocupa total o parcialmente los poros del suelo y del subsuelo.
El agua de infiltración alcanza el subsuelo y conecta con el sistema acuífero general.
10
EVOLUCIÓN DE LA SUPERFICIE FREÁTICA EN TENERIFE
500 500
2000 2000
Desde 1925 a 1997 las entradas de agua por precipitación no han compensado su extracción a
través de galerías y pozos.
250
250 250
Isopiezas cada 50 m.
11
EVOLUCIÓN DE LA SUPERFICIE FREÁTICA EN TENERIFE
CAPA SUPERFICIAL
CAPA FREÁTICA
(límite superior
de saturación)
3000 CAPA SATURADA
ZÓCALO
Altitud (m)
2000
1000
En la capa de saturación se acumula la mayor parte del agua del subsuelo. Está limitada por la
capa freática, que representa el límite superior de la zona de saturación; y el zócalo,
prácticamente impermeable.
La morfología de la capa freática reproduce la topografía de la isla. Su altura es máxima en la
zona de Las Cañadas y desde ahí desciende gradualmente hasta alcanzar los 0 m en la costa.
Posiblemente la escorrentía sea mayor debido a que tan sólo se ha tenido en cuenta la
permeabilidad del terreno en función del tipo de sustrato (coladas recientes o antiguas) pero hay
otro factores relacionados con el uso del suelo que pueden aumentar estos valores de escorrentía.
12
CÁLCULO DE LA EVAPOTRANSPIRACIÓN POTENCIAL POR EL
MÉTODO DE THORNWAITE
a
T
ETP 16 10
I
I = Índice calorífico anual = suma de los valores del índice calorífico mensual (i), que se calcula
como: T
1.514
i
5
T = Temperatura media mensual en ºC.
a = Función polinómica del índice calorífico anual (I), que se calcula como:
a 0.00000067 5 I 3 0.0000771 I 2 0.01792 I 0.49239
Dado que todos los meses no tienen el mismo número de horas de sol y este valor varía con la
latitud es necesario corregir el valor obtenido de evapotranspiración mensual en función del
número de horas de sol reales:
Hs
ETPc ETP
360
El valor obtenido debe ser calibrado según se trate de un clima seco (+ 37%) o húmedo (+ 4%)
Se calcula la evapotranspiración potencial (ETP) mediante la fórmula de Thornwaite, a partir de los datos
mensuales de temperatura media (T) y número de horas de sol. El resultado obtenido se calibra en función
del tipo de clima.
13
CÁLCULO DEL DÉFICIT HÍDRICO
Mes Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic Annual
Temperatura media mensual (T) 14 14 14 15 17 19 21 22 21 19 17 15
Precipitación media mensual (P) 140 112 110 68 47 40 28 33 81 119 131 111 1020
Nº horas de sol 328 314 372 385 420 417 425 406 370 356 323 323
1.514
Índice calorífico mensual i = (T/5) 4.9 4.7 4.9 5.3 6.2 7.4 8.7 9.3 8.9 7.5 6.3 5.4
a
Evapotrasp. Potencial mensual (mm) ETP = 16·(10·T/I) ·(NºHsol/360) 41 37 47 53 69 84 104 108 93 74 54 45 809
Evapotrasp. Potencial mensual (mm) Calibrada para zonas húmedas (+4%) 43 39 49 55 71 88 108 112 97 77 56 47 841
Var. mensual en la reserva de agua G = P-ETP 97 73 62 13 -24 -48 -80 -79 -15 42 75 64 179
Reserva de agua mensual CR = (P-ETP) + CRo si CR < 0 -> CR = 0 (max100) 100 100 100 100 76 28 0 0 0 42 100 100 746
Evapotransp. Real mensual ETR = CR + P si CR+P>ETP-> ETP ó CR+P 43 39 49 55 71 67 28 33 81 77 56 47 646
Déficit hídrico Dh = ETP - ETR 0 0 0 0 0 20 80 79 15 0 0 0 195
Índice calorífico anual I = ∑ i 79
a =0.016·I+0.5
a 0.000000675 I 3 0.0000771 I 2 0.01792 I 0.49239 1.8
mm P = ETR + G + Dh
140
1020 = 646 + 179 + 195
120
100
80
60
40
20
0
Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic
P ETP CR ETR Dh
mm P = ETR + G + Dh mm P = ETR + G + Dh
180 180
363 = 374.2 - 600 + 588.8 127 = 127 - 1337 + 1337
160 160
140 140
120 120
100 100
80 80
60 60
40 40
20 20
0 0
Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic
P ETP CR ETR Dh P ETP CR ETR Dh
14
VARIACIÓN ALTITUDINAL DEL DÉFICIT HÍDRICO EN
TENERIFE
La Laguna 798 m. Izañ
Izaña 2367 m.
mm P = ETR + G + Dh mm P = ETR + G + Dh
140 160
774 = 424 + 45 + 305 504 = 302 - 320 + 522
120 140
120
100
100
80
80
60
60
40
40
20 20
0 0
Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic
Mes Ene Feb Mar Abr May Jun* Jul* Ago* Sep Oct Nov Dic
Nº horas de sol 328 314 372 385 420 417 425 406 370 356 323 323
Reducción del 10% 296 282 335 347 378 292 298 285 333 321 291 290
* Se ha considerado una reducción del 30% debido a la mayor frecuencia del mar de nubes
120 120
100 100
80 80
60 60
40 40
20 20
0 0
Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic
15