0% encontró este documento útil (0 votos)
75 vistas103 páginas

Universidad Mayor de San Andrés Facultad de Agronomía Carrera de Ingeniería Agronómica

Descargar como pdf o txt
Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1/ 103

INGENIERÍA Ingeniería Agronómica

UNIVERSIDAD MAYOR DE SAN ANDRÉS


FACULTAD DE AGRONOMÍA
CARRERA DE INGENIERÍA AGRONÓMICA

TESIS DE GRADO

VALIDACIÓN DE LA METODOLOGÍA DE LA FAO PARA EL CÁLCULO DE ETo CON


DATOS CLIMATICOS LIMITADOS PARA ZONAS DE ALTURA EN BOLIVIA

YUJRA ARUQUIPA RUDDY MILTON

LA PAZ – BOLIVIA
2011

1
INGENIERÍA Ingeniería Agronómica

UNIVERSIDAD MAYOR DE SAN ANDRÉS


FACULTAD DE AGRONOMÍA
CARRERA DE INGENIERÍA AGRONÓMICA

VALIDACIÓN DE LA METODOLOGÍA DE LA FAO PARA EL CÁLCULO DE ETo CON


DATOS CLIMATICOS LIMITADOS PARA ZONAS DE ALTURA EN BOLIVIA

Tesis de Grado presentado como requisito


Parcial para optar el Título de
Ingeniería Agronómica

RUDDY MILTON YUJRA ARUQUIPA

Asesores:

Ing. Ph. Dr. René Chipana R. …………………………………..

Ing. Gladys Tesoro Michel …………………………………..

Tribunal Examinador

Ing. Paulino Ruiz Huanca …………………………………..

Ing. Hugo Bosque Sanchez …………………………………..

Ing. Fernando Manzaneda …………………………………..

Aprobada

Presidente Tribunal Examinador …………………………………..

2011

2
INGENIERÍA Ingeniería Agronómica

PROLOGO

A lo largo de la historia, el hombre ha sabido aprender a dominar su medio


ambiente. Ha desarrollado plantas y variedades de plantas más perfectas, adaptadas a
sus necesidades. Ha concebido prácticas idóneas para utilizar el agua, los fertilizantes y
los plaguicidas del modo más eficaz posible con el objeto de aumentar la producción
agrícola. Pero no ha sabido dominar el clima, y sigue estando bajo la amenaza contante
de la sequia. Al disponer de agua en la cantidad limitada y al crecer la población y ser
necesaria una producción de alimentos mejor y más abundante, el agua a pasado a ser el
recurso natural mas valioso en casi todas las regiones del mundo y, por consiguiente,
resulta absolutamente imperativo llegar a una planificación realmente eficaz del
aprovechamiento del agua en la producción agrícola.

Se han ideado unas metodologías que permiten predecir el volumen exacto de agua
que se requiere para obtener una producción óptima de los distintos cultivos. Estos
métodos se refieren a las condiciones climáticas, agronómicas y edafológicas propias de
una zona dada. La transferencia de metodologías de una zona a otra muy distinta de
aquella en la que han sido concebidas sigue siendo problemática; a menudo se
necesitaran experimentos in situ que llevarán mucho tiempo y requerirán mucho trabajo, y
que a veces serán también muy onerosos (costosos), con objeto de someter a prueba y
calibrar los métodos en un conjunto nuevo de condiciones.

Hay que conocer exactamente la predicción cuantitativa de las necesidades de


riego con respecto a la producción agrícola a efectos de identificación y análisis de la
viabilidad de los proyectos de riego propuestos. Se requiere orientación sobre los métodos
de predicción más prometedores que quepa aplicar para determinar el aprovechamiento
más eficaz posible del agua disponible con fines de riego.

El método presentado en este trabajo fue formulado por el grupo consultivo de la


FAO sobre las necesidades de agua de los cultivos.

3
INGENIERÍA Ingeniería Agronómica

Esta Tesis la dedico con todo mi cariño a la


razón de mi vida mi Mariangel Yujra Cruz, a
mis papás Santiago y Betty por todo su
apoyo, cariño y comprensión. Y a mi familia
por brindarme su apoyo incondicional en todo
momento.

4
INGENIERÍA Ingeniería Agronómica

Quiero agradecer a Dios, por guiar mis pasos,


por darme la fuerza necesaria para seguir adelante y por
la bendición recibida en cada uno de los momentos que
así lo pedí.

Agradezco también al Dr. René Chipana,


quien me dio la oportunidad más grande en toda mi
carrera, poder terminarla.

De manera muy especial al Prof. Luis Santos


Pereira por su apoyo y valiosos aportes durante la
realización de la investigación.

A la Ing. Gladys Tesoro Michel y al


SENAMHI, en especial al Departamento de Archivos, un
sincero agradecimiento.

Agradezco de manera muy especial a Mario


Mendoza Ramos mi Tío quien, sin dudarlo, me ayudó
incondicionalmente en el transcurso de mi carrera.

5
INGENIERÍA Ingeniería Agronómica

ÍNDICE
Pág.
1. INTRODUCCION 1
1.1. Antecedentes 5
1.2. Justificación 6
2. OBJETIVOS 8
2.1. Objetivo General 8
2.2. Objetivos Específicos 8
3. REVISION BIBLIOGRAFICA 9
3.1. Proceso de evapotranspiración 9
3.1.1. Evaporación 9
3.1.2. Transpiración 10
3.2. Unidades 15
3.3. Factores que afectan la evapotranspiración 16
3.3.1. Variables climáticas 16
3.3.2. Factores de cultivo 19
3.3.3. Manejo y condiciones ambientales 20
3.4. ET calculada con datos meteorológicos 20
3.5. Formulación de la ecuación de Penman-Monteith 21
3.5.1. Ecuación de Penman-Monteith 21
3.5.2. Ecuación de FAO Penman-Monteith 21
4. LOCALIZACIÓN 22
5. MATERIALES Y MÉTODOS 23
5.1. Materiales y equipos 23
5.1.1. Datos Meteorológicos 23
5.1.2. Localización 24
5.1.3. Temperaturas 24
5.1.4. Humedad atmosférica 24
5.1.5. Radiación Solar 25
5.1.6. Velocidad del Viento 25
5.2. Equipos 25

6
INGENIERÍA Ingeniería Agronómica

5.3. Metodología 26
5.3.1. Factores meteorológicos que determinan la evapotranspiración 28
5.3.1.1. Temperatura del aire 28
5.3.1.2. Humedad del aire 29
5.3.1.2.1. Presión de vapor 30
5.3.1.2.2. Humedad relativa 30
5.3.1.2.3. Presión media de vapor de la saturación (es) 30
5.3.1.2.4. Pendiente de la curva de presión de
saturación de vapor (Δ) 32
5.3.1.2.5. Presión real de vapor (ea) derivada de
datos de humedad relativa 32
5.3.1.2.6. Déficit de presión de vapor (es - ea) 33
5.3.1.3. Radiación 33
5.3.1.3.1. Radiación extraterrestre (Ra) 33
5.3.1.3.2. Radiación extraterrestre para periodos diarios (Ra) 34
5.3.1.3.3. Duración solar relativa (n/N) 35
5.3.1.3.4. Duración máxima de la insolación (N) 36
5.3.1.3.5. Radiación solar o de onda corta (Rs) 36
5.3.1.3.6. Radiación solar en un día despejado (Rso) 37
5.3.1.3.7. Radiación neta solar o de onda corta (Rns) 38
5.3.1.3.8. Radiación neta de onda larga (Rnl) 38
5.3.1.3.9. Radiación neta (Rn) 39
5.3.1.4. Velocidad del viento 40
5.3.1.4.1. Relación del viento con la altura 40
5.4. Parámetros atmosféricos 41
5.4.1. Presión atmosférica 41
5.4.2. Calor latente de vaporización (λ) 42
5.4.3. Constante Psicrométrica (γ) 42
5.5. Estimación de datos climáticos faltantes 43
5.5.1. Estimación de datos faltantes de humedad 43
5.5.2. Estimación de los datos de radiación 43

7
INGENIERÍA Ingeniería Agronómica

5.5.3. Radiación solar derivada de las diferencias térmicas 44


5.6. Ecuación alternativa para el cálculo de la ETo 45
5.7. ANÁLISIS ESTADISTICO 46
5.7.1. Coeficiente de Determinación: 46
5.7.2. Coeficiente de Regresión 47
5.7.3. Error Cuadrado Medio (RMSE) 48
5.7.4. Promedio Absoluto del Error (AAE) 48
5.7.5. Promedio Relativo del error (ARE) 48
6. RESULTADOS Y DISCUCIONES 49
6.1. Estimación de la humedad a partir de la temperatura mínima 49
6.2. Estimación de la radiación solar a partir de datos diarios de temperatura
máxima y mínima para la estación meteorológica de Belén 51
6.3. Validación del procedimiento cuando no se cuenta con datos
de velocidad del viento 54
6.4. Estimación de la ETo cuando no se cuenta con datos de HR% y Rs. 57
6.5. Estimación de la ETo cuando no se cuenta con datos de HR% y U2. 60
6.6. Estimación de la ETo cuando no se cuenta con datos de Rs y U2. 62
6.7. Estimación de la ETo cuando no se cuenta con datos de HR, Rs y U2. 65
6.8. Cálculo de la ETo por la ecuación de Hargreaves para valores diarios 68
6.9. Hargreaves para periodos de cinco días 71
7. CONCLUSIONES 74
8. RECOMENDACIONES 75
9. BIBLIOGRAFÍA
10. ANEXOS

8
INGENIERÍA Ingeniería Agronómica

ÍNDICE DE ECUACIONES
Pág.
1. Ecuación FAO Penman-Monteith 27
2. Temperatura media del aire (Tmedia) 29
3. Humedad relativa (HR) 30
4. Presión de vapor a saturación en función de la temperatura (e°(T)) 31
5. Presión de vapor a saturación (es) 31
6. Pendiente e°(T) curva (Δ) 32
7. Presión real de vapor obtenida de HRmax y HRmin (ea) 32
8. Radiación extraterrestre para períodos diarios (Ra) 34
9. Conversión de grados decimales a radianes 34
10. Inverso de la distancia relativa Tierra-Sol (dr) 34
11. Declinación solar (δ) 34
12. Angulo horario al ocaso – función arco coseno (ωs) 35
13. Horas de insolación (N) 36
14. Radiación solar (Rs) 37
15. Radiación para días despejados cerca del nivel del mar (Rso) 37
16. Radiación para días despejados para altitudes mayores (Rso) 38
17. Radiación neta solar o de onda corta (Rns) 38
18. Radiación neta de onda larga (Rnl) 39
19. Radiación neta (Rn) 39
20. Ajuste de la velocidad del viento a la altura estándar (u2) 41
21. Presión atmosférica (P) 41
22. Constante psicrométrica (γ) 42
23. Estimación de la presión real de vapor a partir de Tmin (ea) 43
24. Estimación de la radiación solar a partir de diferencias de
temperatura (fórmula de Hargreaves) 44
25. Ecuación de evapotranspiración de referencia de Hargreaves 1985 45
26. Coeficiente de determinación R2 46
27. Coeficiente de Regresión b 47
28. Error cuadrado medio (RMSE) 48

9
INGENIERÍA Ingeniería Agronómica

29. Promedio absoluto del error (AAE) 48


30. Promedio relativo del error (ARE) 48

ÍNDICE DE FIGURAS
Pág.
Figura 1. Representación esquemática de un estoma 12
Figura 2. Estoma visto a través de un microscopio. 12
Figura 3. Repartición de la evapotranspiración en evaporación
y transpiración durante el periodo de crecimiento de un cultivo anual. 15
Figura 4. Factores que afectan la evapotranspiración con conceptos
relacionados de ET. 16
Figura 5. Ubicación Geográfica de las estaciones meteorológicas
utilizadas para el estudio. 23
Figura 6. Características del cultivo de referencia 26
Figura 7. Comparación de ETo – PM con datos completos y ETo
PM con valores estimados de humedad (HR). 49
Figura 8. Comparación de ETo – PM con datos completos y
ETo PM con valores estimados de Radiación Solar (Rs). 52
Figura 9. Comparación de ETo – PM con datos completos y
ETo PM calculado a partir de U2=2 m s-1. 55
Figura 10. Comparación de ETo – PM con datos completos y ETo PM calculado
con valores estimados de Radiación Solar (Rs) y humedad (HR) 58
Figura 11. Comparación de ETo – PM con datos completos y
ETo PM calculado con valores estimados de Humedad (HR)
y velocidad del viento = 2 m s-1. 60
Figura 12. Comparación de ETo – PM con datos completos y
ETo PM calculado con valores estimados de Radiación solar
y velocidad del viento = 2 m s-1 63
Figura 13. Comparación de ETo – PM con datos completos
y ETo PM calculado a partir de temperatura máxima y mínima. 65
Figura 14. Comparación de ETo – PM con datos completos y

10
INGENIERÍA Ingeniería Agronómica

ETo Hargreaves para valores diarios de T máx y mín. 68


Figura 15. Comparación de ETo – PM con datos completos y
ETo Hargreaves para periodos de cinco días 71

ÍNDICE DE CUADROS
Pág.
Cuadro 1. Características y localización de las estaciones meteorológicas
de las que se utilizaran los parámetros climáticos. 22
Cuadro 2. Nubosidad expresada en valores equivalentes de n/N 36

11
INGENIERÍA Ingeniería Agronómica

ÍNDICE DE TABLAS
Pág.
Tabla 1. Indicadores de desempeño para ETo calculados
con datos completos y con datos faltantes de HR. 51
Tabla 2. Indicadores de desempeño para ETo calculados
con datos completos y con datos faltantes de Rs. 53
Tabla 3. Indicadores de desempeño para ETo calculados
con datos completos y con datos faltantes de U2. 56
Tabla 4. Indicadores de desempeño para ETo calculados
con datos completos y con datos faltantes de HR y Rs. 59
Tabla 5. Indicadores de desempeño para ETo calculados
con datos completos y con datos faltantes de HR Y U2. 62
Tabla 6. Indicadores de desempeño para ETo calculados
con datos completos y con datos faltantes de Rs y U2. 64
Tabla 7. Indicadores de desempeño para ETo calculados
con datos completos y con datos faltantes de HR, R y U2. 67
Tabla 8. Indicadores de desempeño para ETo calculados
con datos completos y ETo por el método de Hargreaves. 70
Tabla 9. Indicadores de desempeño para ETo calculados
con datos completos y ETo por el método de Hargreaves
para periodos de cinco días. 73

12
INGENIERÍA Ingeniería Agronómica

Lista de los símbolos


Principales y acrónimos

apsi Coeficiente del psicrómetro [°C-1]


as Fracción de la radiación extraterrestre que llega a la tierra en un día nublado [-]
as+bs Fracción de la radiación extraterrestre que llega a la tierra en un día
despejado [-]
cp Calor específico [MJ kg-1 °C-1]
cs Capacidad calorífica del suelo [MJ m-3 °C-1]
dr Inverso de la distancia relativa Tierra-Sol [-]
E Evaporación [mm día-1]
e°(T) Presión de vapor a saturación a la temperatura del aire T [kPa]
es Presión de vapor a saturación para un período determinado [kPa]
ea Presión real de vapor [kPa]
es - ea Déficit de presión de vapor con respecto a la saturación
ET Evapotranspiración [mm día-1]
ETo Evapotranspiración del cultivo de referencia [mm día-1]
ETc Evapotranspiración del cultivo bajo condiciones estándar [mm día-1]
exp[x] 2,7183 (base del logaritmo natural) elevado a la potencia x
G Flujo de calor en el suelo [MJ m-2 día-1]
Gsc Constante solar [0,0820 MJ m-2 min-1]
HR Humedad relativa [%]
HRmáx Humedad relativa máxima diaria [%]
HRmedia Humedad relativa media diaria [%]
HRmín Humedad relativa mínima diaria [%]
h Altura del cultivo [m]
J Número del día en el año [-]
Kc Coeficiente del cultivo [-]
kRs Coeficiente de ajuste para la formula de radiación de Hargreaves [°C-0,5]
N Duración máxima de la insolación durante el día o el máximo número de
horas de sol durante el día [hora]

13
INGENIERÍA Ingeniería Agronómica

n Duración real de la insolación durante el día o el número real de horas de


sol durante el día [hora]
n/N Relación de la duración de la insolación o fracción de insolación [-]
P Precipitación [mm], presión atmosférica [kPa]
Ra Radiación extraterrestre [MJ m-2 día-1]
Rl Radiación de onda larga [MJ m-2 día-1]
Rn Radiación neta [MJ m-2 día-1]
Rnl Radiación neta de onda larga [MJ m-2 día-1]
Rns Radiación solar neta o radiación neta de onda corta [MJ m -2 día-1]
Rs Radiación solar o radiación de onda corta [MJ m-2 día-1]
Rso Radiación solar o de onda corta para cielo despejado [MJ m -2 día-1]
ra Resistencia aerodinámica [s m-1]
rs Resistencia del cultivo o resistencia superficial (total) [s m-1]
Rs/Rso Fracción de radiación solar o fracción de radiación de onda corta [-]
T Temperatura del aire [°C]
Tmáx Temperatura máxima diaria del aire [°C]
Tmedia Temperatura media diaria del aire [°C]
Tmín Temperatura mínima diaria del aire [°C]
u2 Velocidad del viento a 2 m sobre la superficie del suelo [m s-1]
uz Velocidad del viento a una altura z m sobre la superficie del suelo [m s-1]
z Altitud, altura sobre el nivel del mar [m]
α Albedo [-]
γ Constante psicrométrica [kPa °C-1]
Δ Pendiente de la curva de presión de vapor a saturación [kPa °C-1]
δ Declinación solar [rad]
ε Cociente entre los pesos moleculares del vapor de agua/aire seco
(= 0,622)
λ Calor latente de vaporización [MJ kg-1]
λET Flujo de calor latente [MJ m-2 día-1]
ω Ángulo del tiempo solar a la mitad del período horario o inferior [rad]
ωs Ángulo horario al ocaso [rad]

14
INGENIERÍA Ingeniería Agronómica

VALIDACIÓN DE LA METODOLOGÍA DE LA FAO PARA EL CÁLCULO DE ETo CON


DATOS CLIMATICOS LIMITADOS PARA ZONAS DE ALTURA EN BOLIVIA

1. INTRODUCCION

La evaporación y la transpiración ocurren simultáneamente y no hay una manera sencilla


de distinguir entre estos dos procesos. Combinados representan lo que se denomina
como evapotranspiración1.

Debido a la importancia de estos procesos para el uso en la planificación de proyectos,


producción de cultivos, manejo de recursos hídricos, asesoramiento ambiental y la
elaboración de calendarios de riego se consideró necesario realizar esta investigación que
nos permitirá contar con datos disponibles de evapotranspiración, a la falta de datos
climáticos en zonas de altura para un mejor aprovechamiento del agua; considerando que
en el momento de la siembra casi el 100% de la ET ocurre en forma de evaporación,
mientras que cuando la cobertura vegetal es completa, más del de 90% de la ET ocurre
como transpiración2.

El clima, las características del cultivo, el manejo y el medio de desarrollo son factores que
afectan la evaporación y la transpiración. Dentro los principales parámetros climáticos que
afectan la evapotranspiración están la radiación, la temperatura del aire, la humedad
atmosférica y la velocidad del viento3, estos datos, generalmente no se encuentran de
manera accesible, por lo que se plantea su cálculo mediante la metodología de la FAO
para datos meteorológicos faltantes.

1
ALLEN Richard G.; PEREIRA Luis S., DIRK Raes, MARTIN Smith, 2006. Evapotranspiración del Cultivo. Guías para
la determinación de los requerimientos de agua de los cultivos. Servicio de Recursos, Fomento y Aprovechamiento de
Aguas. FAO (Organización de las Naciones Unidas para la Agricultura y la Alimentación) Roma, estudio FAO riego y
drenaje Nº 56. Pág. 323.
2
BURMAN, R. y POCHOP, L.O. 1994.La evaporación, datos de evapotranspiración y Climática. Elsevier BV Ciencia,
Amsterdam.
3
ALLEN, R.G., SMITH, M., PERRIER, A., y PEREIRA, L.S. 1994a. Una actualización para la definición de
evapotranspiración de referencia. ICID Boletin. 43(2). Pág. 1-34.

15
INGENIERÍA Ingeniería Agronómica

El Altiplano de Bolivia, constituye una parte importante de la región andina en América del
Sur, las tierras agrícolas se encuentran entre los 3.400 a 4.000 m.s.n.m., donde las
precipitaciones se concentran entre los meses de diciembre a marzo, variando entre 80 a
800 mm por año, cerca al Salar de Uyuni en el Sur de esta región y entorno al Lago
Titicaca, respectivamente. El clima varía entre árido y semiárido, y la agricultura
predominantemente es asecano.

Aunque la mayoría de los cultivos locales están adaptados para producir en las duras
condiciones climáticas y extremas, los rendimientos son realmente bajos debido a que los
cultivos deben sacrificar la producción para desarrollar mecanismos de resistencia a la
sequía y las heladas4. Bajo estas consideraciones, los sistemas de riego serían necesarios
para mejorar la agricultura y la producción, así reducir las pérdidas debidas a las
condiciones climáticas o la baja disponibilidad del agua.

En este entorno, la escasez de agua es una limitación importante para la producción de


cultivos debido a los bajos niveles de precipitación, una tasa relativamente alta de
evapotranspiración y suelos con baja capacidad de retención de agua. Además, como
consecuencia de las bajas temperaturas registradas en invierno, el período de cultivo sólo
se abarca desde octubre a hasta abril.

Para este propósito, las estimaciones de la evapotranspiración fiable para la región, que
requiere el diseño y la gestión de sistemas de riego, en muchos casos son inexistentes
debido a la falta de datos climáticos para el Altiplano, como consecuencia, la
evapotranspiración de referencia (ETo) es frecuentemente estimada con la ayuda de
fórmulas como Thornthwaite (1954) y Hargreaves y Samani (1982, 1985).

Debido a ello, en las últimas décadas se han venido implementando varios proyectos de
riego con márgenes elevados de incertidumbre en el cálculo de la demanda de agua de
los cultivos, aspecto que se torna más significativo, puesto que los recursos hídricos
también son limitados.
4
BOSQUE, H., LEMEUR, R., VAN DAMME, P., 2000. Análisis Ecofisiológico del cultivo de la quinua (Chenopodium
quinoa Wild) en condiciones de estrés de sequía y salinidad. Tropicultura 18, 198–202.

16
INGENIERÍA Ingeniería Agronómica

Con las características ya mencionadas en las que se encuentra el altiplano boliviano,


donde los recursos hídricos son limitados y seriamente amenazada por la
sobreexplotación, es esencial la estimación de la necesidad de agua de los cultivos con
una mayor precisión. De esta manera, la buena gestión y planificación de los recursos
hídricos disponibles serán alcanzadas.

Las necesidades de agua deben ajustarse a la demanda atmosférica, que está


relacionada con las condiciones climáticas del medio en las que se desarrolla el cultivo, y
para hacer frente a estas cuestiones, la Organización de las Naciones Unidas para la
Agricultura y la Alimentación (FAO) propone una metodología para el cálculo de
evapotranspiración del cultivo basado en la utilización de la evapotranspiración de
referencia (ETo) y coeficientes de cultivo (Kc).

Estos coeficientes dependen de varios factores incluyendo tipo de cultivo, la fase del
crecimiento de los cultivos, la altura del dosel y su densidad5.

Para el cálculo diario de la evapotranspiración de referencia (ETo), el método de la FAO-


56 PM requiere datos diarios de temperatura del aire, máxima y mínima (Tmax y Tmin),
humedad relativa (HR), la radiación solar (R) y velocidad del viento (u). Lamentablemente,
para muchos lugares, como es el caso de Bolivia, las variables meteorológicas son a
menudo incompletas y/o no disponibles.

Desafortunadamente en los países en desarrollo como Bolivia existen pocas estaciones


meteorológicas, en las que además, las variables meteorológicas no son medidas
adecuadamente. Esta falta de datos meteorológicos fue resuelto parcialmente por
Hargreaves et al. (1985) desarrollando una relación sencilla para el cálculo de ETo. La
ecuación de Hargreaves-Samani6 requiere solamente datos medios de temperatura
máxima y mínima, usualmente disponible en la mayoría de las estaciones meteorológicas.

5
ALLEN, R.G., PEREIRA, L.S., RAES, D., SMITH, M., 1998. Crop Evapotranspiration-Guidelines for Computing Crop
Water Requirements. FAO Irrigation and drainage paper 56. United. Nations Food and Agriculture Organization, Rome.
6
HARGREAVES, G.H., SAMANI, Z.A., 1985. Reference crop evapotranspiration from temperature. Appl. Eng. Agric. 1
(2), 96–99.

17
INGENIERÍA Ingeniería Agronómica

Este método es adecuado para periodos semanales o mayores, aunque ha sido reportado
en literatura una cierta mejora para estimaciones diarias7 .

Gavilán et al. (2005), al efectuar la calibración regional de la ecuación de Hargreaves para


zonas semiáridas, encontraron que esta ecuación subestima la evapotranspiración de
referencia en zonas cercanas a las costas (con altas velocidades del viento y pequeña
variación de ΔT ) y sobrestima en las zonas interiores (con altas variaciones de ΔT y
bajas velocidades del viento), debido a que las medias anuales de la velocidad del viento
y el rango de variación de la temperatura influye en el comportamiento de dicha ecuación.
Es más adecuado para lugares ventosos con grande variación de la temperatura (ΔT) y
para las situaciones con condiciones de viento ligeras combinadas con bajos a moderados
valores de (ΔT).

Por otra parte, Allen et al. (1998), propone el uso de la ecuación de Hargreaves (HG)
como una alternativa de estimación para el cálculo de la ETo, ecuación que utiliza como
datos, la temperatura del aire (Máxima y mínima). Este método se comporta mejor para
datos semanales o más (Hargreaves y Allen, 2003). Por lo tanto, es importante evaluar la
exactitud de los procedimientos para calcular los datos faltantes a partir de ETo.

Aunque la FAO considera la ecuación Penman-Monteith como el único método de


referencia para determinar la ETo (Allen et al., 1998), este método rara vez se aplica
debido a su aparente complejidad y la necesidad excesiva de datos meteorológicos (en su
mayoría no está disponible). Sin embargo, la FAO, y el método de Penman-Monteith se
cree es la más fiable para la mayoría de las condiciones climáticas, gracias a la
incorporación de factores tanto físicos y biológicos y es ampliamente utilizado en el
mundo8.

7
HARGREAVES, G.H., ALLEN, R.G., 2003. History and evaluation of Hargreaves evapotranspiration equation. J. Irrig.
Drain. Eng. ASCE 129 (1), 53–63.
8
HESS, T., 1998. Tendencias en la evapotranspiración de referencia en el Norte Oriental de Zonas Áridas de Nigeria,
1961-1991. Journal of Arid Environments. 38, 99-115.

18
INGENIERÍA Ingeniería Agronómica

1.1 Antecedentes

FAO (2006), Indica que durante los últimos 50 años se han desarrollado una gran cantidad
de métodos más o menos empíricos por numerosos científicos y especialistas en todo el
mundo, con el fin de estimar la evapotranspiración a partir de diferentes variables
climáticas tales como los estudios comparativos de algunos métodos frecuentemente
usados para el cálculo de ETo, como el de Ivanov, la ecuación Hershkovich 1984, la
ecuación de Penman modificada para el césped y la referencia de la alfalfa y, más
recientemente, el método de FAO Penman-Monteith.

Según FAO (2006), los avances en la investigación y la valoración más precisa del uso de
agua del cultivo han revelado debilidades en las metodologías propuestas. Numerosos
investigadores analizaron el funcionamiento de los cuatro métodos para diversas
localidades. Aunque los resultados de tales análisis podrían haber sido influenciados por el
sitio, las condiciones de medida, el sesgo en la recopilación de los datos climáticos, llegó a
ser evidente, que los métodos propuestos no se comportan de la misma manera en
diversas zonas del mundo encontrando que:

 Los métodos basados en la ecuación de Penman podrían requerir una calibración


local de la función del viento para obtener resultados satisfactorios.
 Los métodos de la radiación demostraron buenos resultados en climas húmedos
donde el término aerodinámico es relativamente pequeño, pero su funcionamiento
en condiciones áridas fue errático y con tendencia a subestimar la
evapotranspiración.
 Los métodos basados en la temperatura del aire al ser empíricos, requieren
cuidadosas calibraciones locales para alcanzar resultados satisfactorios.
 Los métodos basados en la evaporación del tanque reflejan las limitaciones de
estimar la evapotranspiración del cultivo basándose en datos de evaporación libre
del agua. Estos métodos son susceptibles a las condiciones microclimáticas bajo
las cuales los tanques están funcionando y al rigor de mantenimiento de la
estación. Muchas veces su funcionamiento se ha presentado errático.

19
INGENIERÍA Ingeniería Agronómica

 La fórmula de FAO Penman-Monteith ha sido reportada por ambos equipos de


investigación (ASCE y europeo) como aquella que produce resultados
relativamente exactos y consistentes tanto en climas áridos como en húmedos.

Popova et al. (2006), validaron la metodología de la FAO en el sur de Bulgaria, para el


cálculo de ETo con datos limitados y encontraron que el procedimiento propuesto por la
FAO es adecuado proveyendo pequeños valores de error estándar, cuando comparado
con la ecuación de Hargreaves que tiende a sobrestimar la ETo9. En la misma línea
Jabloun y Salí (2008), evaluaron el método propuesto por la FAO para el cálculo de ETo
con datos limitados en Tunez. Encontraron que la diferencia entre los valores de ETo
calculados con la ecuación completa y con datos limitados son pequeñas, siendo mayor
cuando se usó la ecuación de Hargreaves10.

Para zonas de altura García et al. (2004), encontraron que la evapotranspiración de


referencia (ETo) calculada con la ecuación FAO-56 PM presenta similitud con los datos
observados en campo, en tanto que el método de Hargreaves es adecuado para
determinar la ETo en el altiplano norte y no así en el altiplano sur, debido a la exclusión del
factor aerodinámico.

1.2 Justificación

En Bolivia existe un número reducido de estaciones meteorológicas automáticas que


generan todos los parámetros climáticos necesarios para aplicar la ecuación FAO-56 PM,
la mayoría son estaciones convencionales o solo termopluviométricas, (miden solamente
la precipitación pluvial y temperatura), adjuntando a esto que, dichos datos no pasan por
un análisis de calidad que asegure su confiabilidad lo que limita aun mas la generación de
datos de evapotranspiración (ET).

9
POPOVA, Z., KERCHEVA, M., PEREIRA, L.S., 2006. Validation of the FAO methodology for computing ETo with
missing climatic data application to South Bulgaria. Irrig. Drain. 55, Pág. 201–215.
10
JABLOUN, M., SAHLI, A., 2008. Evaluation of FAO-56 methodology for estimating reference evapotranspiration
using limited climatic data Application to Tunisia. Agric. Water Manage. 95, 707–715.

20
INGENIERÍA Ingeniería Agronómica

Estos aspectos dificultan el cálculo preciso de la evapotranspiración, lo que lleva a una


subestimación o sobrestimación de las necesidades de riego. No obstante Allen et al.
(1998) propusieron métodos de estimación de la ET aplicando la ecuación FAO-56 PM, a
partir de datos limitados.

En estaciones de Bolivia probar la exactitud de los métodos bajo nuevas condiciones es


laborioso y necesita mucho tiempo y dinero. A pesar de ello, los datos de
vapotranspiración son necesarios con antelación al planeamiento de proyectos o para
programar calendarios de riego11.

El análisis del funcionamiento y la dispersión de resultados obtenidos por varios métodos


de cálculo revelan la necesidad de formular un único método estándar para el cómputo de
la ETo. Dado el mejor comportamiento y consistencia del método FAO Penman-Monteith
en el ámbito global, este es el único que se recomienda como método estándar de cálculo.
Esto se debe a que existe una elevada probabilidad de que este método prediga
correctamente los valores de ETo en una amplia gama geográfica y climática y cuenta con
previsiones para su uso en situaciones de falta de datos12.

11
PENMAN, H. L. 1948. La evaporación natural del agua abierta, tierra desnuda y césped. Proc. Roy. Soc. Londres,
A193, 120-146 p.
12
SMITH, M., ALLEN, R.G., MONTEITH, J.L., PERRIER, A., PEREIRA, L., y SEGEREN, A. 1992. Report of the expert
consultation on procedures for revision of FAO guidelines for prediction of crop water requirements. Organización de las
Naciones Unidas para la Agricultura y la Alimentación, Roma, 54 p.

21
INGENIERÍA Ingeniería Agronómica

11. OBJETIVOS

11.1. Objetivo General

 Validar la metodología de la FAO para el cálculo de la evapotranspiración (ETo)


con datos climáticos limitados para zonas de altura en Bolivia.

11.2. Objetivos Específicos

 Evaluar los resultados obtenidos con la ecuación de la FAO Penman-Monteith con


datos limitados y los obtenidos con datos completos.
 Efectuar un análisis de consistencia de la información meteorológica diaria de
precipitación, radiación solar, humedad relativa y velocidad del viento en zonas
altas de Bolivia.
 Evaluar la exactitud de los procedimientos de cálculo para las condiciones del
altiplano boliviano.
 Comparar la evapotranspiración en diversos períodos del año y en otras regiones
del altiplano boliviano

22
INGENIERÍA Ingeniería Agronómica

12. REVISION BIBLIOGRAFICA


12.1. La evapotranspiración

Se conoce como evapotranspiración (ET), a la combinación de dos procesos separados


por los que el agua se pierde, a través de la superficie del suelo por evaporación y por otra
parte mediante transpiración del cultivo13.

La evapotranspiración es la combinación de la evaporación desde la superficie del suelo y


la traspiración de la vegetación14.

Evapotranspiración es el resultado del proceso por el cual, el agua cambia de estado


líquido a gaseoso, y directamente, o a través de las plantas, vuelve a la atmósfera en
forma de vapor.

El término sólo es aplicable correctamente a una determinada área de terreno cubierta por
vegetación. Ante la ausencia de vegetación, sólo se puede hablar de evaporación.

La evapotranspiración (ET), es el proceso por el cual el agua es transferida desde la


superficie terrestre hacia la atmósfera. Incluye tanto la evaporación de agua en forma
sólida como líquida directamente del suelo o desde las superficies vegetales vivas o
muertas (rocío, escarcha, lluvia interceptada por la vegetación), como las pérdidas de
agua a través de las superficies vegetales, particularmente las hojas15.

12.1.1. Evaporación

La evaporación es el proceso por el cual el agua líquida se convierte en vapor de agua


(vaporización) y se retira de la superficie evaporante (remoción de vapor). El agua se
evapora de una variedad de superficies, tales como lagos, ríos, caminos, suelos y la
vegetación mojada16.

13
ALLEN Richard G. et al., 2006, Pág. 1
14
CHAVARRI VELARDE, Eduardo A. 2004, Universidad Nacional Agraria La Molina, Facultad de Ingeniería Agrícola,
Departamento de Recursos de Agua y Tierra. Hidrología Aplicada, Pág. 1-10.
15
http://www.miliarium.com/Proyectos/EstudiosHidrogeologicos/Memoria/Evapotranspiracion/evapotranspiracion.asp.
16
MONTEITH, J.L., 1965. La evaporación y Medio Ambiente. Simposios 19 de la Sociedad para la Biología
Experimental, University Press, Cambridge, 19:205-234.

23
INGENIERÍA Ingeniería Agronómica

Para cambiar el estado de las moléculas del agua de líquido a vapor se requiere energía.
La radiación solar directa y, en menor grado, la temperatura ambiente del aire,
proporcionan esta energía. A medida que ocurre la evaporación, el aire circundante se
satura gradualmente y el proceso se vuelve cada vez más lento hasta detenerse
completamente si el aire mojado circundante no se transfiere a la atmósfera o en otras
palabras no se retira de alrededor de la hoja.

El reemplazo del aire saturado por un aire más seco depende grandemente de la
velocidad del viento (Hanks, R.J. y Hill, R.W. 1980).

Por lo tanto, la radiación, la temperatura del aire, la humedad atmosférica y la velocidad


del viento son parámetros climatológicos a considerar al evaluar el proceso de la
evaporación17.

Cuando la superficie evaporante es la superficie del suelo, el grado de cobertura del suelo
por parte del cultivo y la cantidad de agua disponibles en la superficie evaporante son
otros factores que afectan el proceso de la evaporación. Lluvias frecuentes, el riego y el
ascenso capilar en un suelo con manto freático poco profundo, mantienen mojada la
superficie del suelo. En ausencia de cualquier fuente de reabastecimiento de agua a la
superficie del suelo, la evaporación disminuye rápidamente y puede cesar casi totalmente
en un corto lapso de tiempo18.

3.1.2. Transpiración

La transpiración es el proceso mediante el cual las plantas regresan el agua a la


atmósfera. Después de absorber agua del suelo, las plantas liberan agua a través de sus
hojas. La transpiración ayuda a la plantas a mantenerse frescas, de la misma forma que la
respiración mantiene frescos a los seres humanos y animales19.

17
MONTEITH, J.L., 1965.
18
PENMAN, H. L. 1948. 120-146 p.
19
http://www.windows.ucar.edu/tour/link=/earth/Water/transpiration.sp.html

24
INGENIERÍA Ingeniería Agronómica

La transpiración consiste en la vaporización del agua líquida contenida en los tejidos de la


planta y su posterior remoción hacia la atmósfera. Los cultivos pierden agua
predominantemente a través de los estomas. Estos son pequeñas aberturas en la hoja de
la planta a través de las cuales atraviesan los gases y el vapor de agua de la planta hacia
la atmósfera20.

El agua, junto con algunos nutrientes, es absorbida por las raíces y transportada a través
de la planta. La vaporización ocurre dentro de la hoja, en los espacios intercelulares, y el
intercambio del vapor con la atmósfera es controlado por la abertura estomática. Casi toda
el agua absorbida del suelo se pierde por transpiración y solamente una pequeña fracción
se convierte en parte de los tejidos vegetales.

La transpiración, igual que la evaporación directa, depende del aporte de energía, del
gradiente de presión del vapor y de la velocidad del viento. Por lo tanto, la radiación, la
temperatura del aire, la humedad atmosférica y el viento también deben ser considerados
en su determinación21.

El contenido de agua del suelo y la capacidad del suelo de conducir el agua a las raíces
también determinan la tasa de transpiración, así como la salinidad del suelo y del agua de
riego. La tasa de transpiración también es influenciada por las características del cultivo, el
medio donde se produce y las prácticas de cultivo22.

Diversas clases de plantas pueden tener diversas tasas de transpiración. Por otra parte,
no solamente el tipo de cultivo, sino también su estado de desarrollo, el medio donde se
produce y su manejo, deben ser considerados al evaluar la transpiración23.

20
Véase en la Figura 1
21
PEREIRA, L.S., PERRIER, A., ALLEN, R.G. y ALVES, I. 1996. Evapotranspiración: Revisión de conceptos y
tendencias de flujo. J. Riego y Drenaje. Ing. Agro. ASCE 25. (prensa).
22
PENMAN, H.L. 1963. La vegetación y la hidrología. Tech. Com. N º 53, de la Commonwealth Mesa de Suelos,
Harpenden, Reino Unido. 125 p.
23
PENMAN, H.L. 1963. 125 p.

25
INGENIERÍA Ingeniería Agronómica

FIGURA 1.
Representación esquemática de un estoma

Fuente: L. Pereira 1996

FIGURA 2.
Estoma visto a través de un microscopio.

Fuente: (David Sierra M., 2007, Articulo 240).

26
INGENIERÍA Ingeniería Agronómica

La transpiración es un traspaso de agua por parte de la planta hacia la atmósfera que


presenta un gradiente de humedad menor que el de la planta24.

En la Figura 2 se muestra las aberturas microscópicas (estomas) por los que la planta
realiza sus intercambios gaseosos y están formados por dos células con forma de judía
(células oclusivas) que rodean el agujero u ostiolo. Las células contiguas a las oclusivas
se llaman células accesorias y están relacionadas con el mecanismo de apertura del
estoma.

Los estomas son, pues, las "bocas" por las que respira la planta. La transpiración
considerada de forma normal es un mecanismo usado por la planta para realizar diversas
funciones:

Refrigerar las hojas: se pueden conseguir temperaturas en las hojas de hasta 15


ºC menos que en el aire que rodea la planta.

Dirigir y concentrar nutrientes absorbidos por las raíces.

La transpiración juega un papel muy importante en la subida de agua por el


xilema. Para que se produzca una transpiración acusada deben actuar distintos
factores externos, como el de la humedad atmosférica:

- La humedad en el suelo es otro factor importante ya que provoca una


mayor apertura de los estomas. Esto no debe tomarse como algo negativo
sino como lo normal en toda planta cuando tiene un sustrato húmedo, es decir,
no hay que dejar que se seque el sustrato pensando que se cierran más
estomas ya que el efecto es mucho peor. Hay que procurar aportar al suelo la
humedad suficiente para compensar la pérdida.

24
SIERRA, M. David, 2007, Apuntes de Botánica. Transpiración. Botanica Articulo 240 http://jardinactual.com/menu-
revista-articulos/240 LA_TRANSPIRACIÓN_DE_LAS_PLANTAS

27
INGENIERÍA Ingeniería Agronómica

- La temperatura: es evidente que las altas temperaturas provocan elevadas


pérdidas de agua y además de forma muy rápida.

- La velocidad del viento: el viento arrastra la capa de humedad que rodea el


estoma y deseca esa zona. La pérdida de agua es mayor cuanto mayor es la
velocidad del viento.

- La luz: provoca una mayor apertura de estomas y, por tanto, mayor transpiración.
Al ponerla en zonas sombreadas o con poca luz, la planta reacciona cerrando
estomas ya que sin luz no hay fotosíntesis ni intercambio de CO225.

Para disminuir los efectos de la transpiración algunas plantas presentan diferentes


sistemas de lucha como cutículas cerosas, pelos, hundimiento de estomas, espinas, hojas
enrolladas, etc.

La transpiración es, en definitiva, un mecanismo muy importante para la planta pero que
se torna peligroso en condiciones de calor y sequedad, llegando incluso a hacer inútiles
todos los esfuerzos llevados a cabo por la propia planta y por la ayuda que se ha podido
aportar desde el exterior26.

La evaporación y la transpiración ocurren simultáneamente y no hay una manera sencilla


de distinguir entre estos dos procesos. Aparte de la disponibilidad de agua en los
horizontes superficiales, la evaporación de un suelo cultivado es determinada
principalmente por la fracción de radiación solar que llega a la superficie del suelo27.

Esta fracción disminuye a lo largo del ciclo del cultivo a medida que el dosel del cultivo
proyecta más y más sombra sobre el suelo. En las primeras etapas del cultivo, el agua se
pierde principalmente por evaporación directa del suelo, pero con el desarrollo del cultivo y
finalmente cuando este cubre totalmente el suelo, la transpiración se convierte en el
proceso principal.

25
SIERRA, M. David, 2007.
26
SIERRA, M. David, 2007.
27
JENSEN, M.E., BURMAN, R.D., y ALLEN, R.G. (ed). 1990. Evapotranspiración y requerimientos de Agua para riego.
Manual de ASCE, Prácticas para diseñar Informes Nº 70, Es. Soc. Ing. civil., Nueva York, Estados Unidos de América,
360 p.

28
INGENIERÍA Ingeniería Agronómica

En la figura 3 se presenta la evapotranspiración dividida en sus dos componentes


(evaporación y transpiración) en relación con el área foliar por unidad de superficie de
suelo debajo de él. En el momento de la siembra, casi el 100% de la ET ocurre en forma
de evaporación, mientras que cuando la cobertura vegetal es completa, más del de 90%
de la ET ocurre como transpiración28.

FIGURA 3.
Repartición de la evapotranspiración en evaporación y transpiración durante el periodo de
crecimiento de un cultivo anual.

Fuente: L. Pereira 1996

12.2. Unidades
La evapotranspiración se expresa normalmente en milímetros (mm) por unidad de tiempo.
Esta unidad expresa la cantidad de agua perdida de una superficie cultivada en unidades
de altura de agua. La unidad de tiempo puede ser una hora, día, 10 días, mes o incluso un
período completo de cultivo o un año29.

28
PERRIER, A. 1978. Importance des définitions de l’évapotranspiration dans le domaine pratique de la
mesure, de l’estimation of de la notion de coefficients culturaux. XV’ Journal of Hydraulics, Société
Hydrotechnique de France, Question IV, Rapport 1:1-7 (en Francés).
29
PERRIER, A. 1978.

29
INGENIERÍA Ingeniería Agronómica

Como una hectárea tiene una superficie de 10 000 m 2 y 1 milímetro es igual a 0,001 m,
una pérdida de 1 mm de agua corresponde a una pérdida de 10 m 3 de agua por hectárea.
Es decir 1 mm día-1 es equivalente 10 m3 ha-1 día-1.

12.3. Factores que afectan la evapotranspiración

El clima, las características del cultivo, el manejo y el medio de desarrollo (figura 4) son
factores que afectan la evaporación y la transpiración30.

FIGURA 4.
Factores que afectan la evapotranspiración con conceptos relacionados de ET

Fuente: Elaboración Propia

3.3.1. Variables climáticas

Los principales parámetros climáticos que afectan la evapotranspiración son la radiación,


la temperatura del aire, la humedad atmosférica y la velocidad del viento31.

30
PERRIER, A. 1978.
31
ALLEN Richard G. et al, 2006. Pág. 5

30
INGENIERÍA Ingeniería Agronómica

a) Radiación solar. Es el factor determinante de la evapotranspiración ya que es la


fuente de energía de dicho proceso. El proceso de la evapotranspiración está
determinado por la cantidad de energía disponible para evaporar el agua. La
radiación solar es la más importante fuente de energía en el planeta y puede
cambiar grandes cantidades de agua líquida en vapor de agua. La cantidad
potencial de radiación que puede llegar a una superficie evaporante viene
determinada por su localización y época del año. Debido a las diferencias en la
posición del planeta y a su movimiento alrededor del sol, esta cantidad potencial de
radiación es diferente para cada latitud y para las diversas estaciones del año.

La radiación solar real que alcanza la superficie evaporante depende de la turbidez


de la atmósfera y de la presencia de nubes que reflejan y absorben cantidades
importantes de radiación. Cuando se determina el efecto de la radiación solar en la
evapotranspiración, se debe también considerar que no toda la energía disponible
se utiliza para evaporar el agua. Parte de la energía solar se utiliza también para
calentar la atmósfera y el suelo32.

b) Temperatura del aire. El aumento de temperatura en el aire facilita la


evapotranspiración ya que: En primer lugar, crea una convección térmica
ascendente, que facilita la aireación de la superficie del líquido; y por otra parte, la
presión de vapor de saturación es más alta.

La radiación solar absorbida por la atmósfera y el calor emitido por la tierra elevan
la temperatura del aire. El calor sensible del aire circundante transfiere energía al
cultivo y entonces ejerce un cierto control en la tasa de evapotranspiración. En un
día soleado y cálido, la pérdida de agua por evapotranspiración será mayor que en
un día nublado y fresco33.

32
ALLEN Richard G. et al, 2006. Pág. 29
33
ALLEN Richard G. et al, 2006. Pág. 32

31
INGENIERÍA Ingeniería Agronómica

c) Humedad atmosférica. Es un factor determinante en la evapotranspiración ya que


para que ésta se produzca, es necesario que el aire próximo a la superficie de
evaporación o transpiración no esté saturado (situación que es facilitada con
humedad atmosférica baja).

Mientras que el aporte de energía del sol y del aire circundante es la fuerza
impulsora principal para la evaporación del agua, la diferencia entre la presión de
vapor de agua en la superficie evapotranspirante y el aire circundante es el factor
determinante para la remoción de vapor. Áreas bien regadas en regiones áridas
secas y calientes, consumen grandes cantidades de agua debido a la gran
disponibilidad de energía y al poder de extracción de vapor de la atmósfera. En
cambio en regiones húmedas tropicales, a pesar de que el ingreso de energía es
elevado, la alta humedad del aire reducirá la demanda de evapotranspiración.

En este último caso, como el aire está ya cerca de saturación, puede absorber
menos agua adicional y por lo tanto la tasa de evapotranspiración es más baja que
en regiones áridas34.

d) El viento. Después de la radiación es el más importante, ya que renueva el aire


próximo a la superficie de evaporación o transpiración que está saturado. La
combinación de humedad atmosférica baja y viento resulta ser la que produce
mayor evapotranspiración.

El proceso de remoción de vapor depende en alto grado del viento y de la


turbulencia del aire, los cuales transfieren grandes cantidades de aire hacia la
superficie evaporante. Con la evaporación del agua, el aire sobre la superficie
evaporante se satura gradualmente con vapor. Si este aire no se substituye
continuamente por un aire más seco, disminuye la intensidad de remoción de vapor
de agua y la tasa de evapotranspiración disminuye.

34
ALLEN Richard G. et al, 2006. Pág. 33.

32
INGENIERÍA Ingeniería Agronómica

Bajo estas circunstancias, mucho vapor de agua puede ser almacenado en el aire
mientras que el viento puede promover el transporte del agua permitiendo que se
retire mayor cantidad de vapor de agua. Por otra parte, bajo condiciones
atmosféricas húmedas, la alta humedad del aire y la presencia de nubes hacen que
la tasa de evapotranspiración sea más baja.

Por tanto, bajo condiciones húmedas la velocidad del viento afecta la


evapotranspiración en un grado mucho menos importante que bajo climas áridos en
los que variaciones pequeñas en la velocidad del viento pueden dar lugar a
importantes variaciones en la evapotranspiración35.

El viento también produce un efecto secundario que es el enfriamiento de la superficie del


líquido y la consiguiente disminución de la evaporación36.

Se han desarrollado varios procedimientos para determinar la evaporación a partir de


estos parámetros. La fuerza evaporativa de la atmósfera puede ser expresada por la
evapotranspiración del cultivo de referencia (ETo). La evapotranspiración del cultivo de
referencia (ETo), representa la pérdida de agua de una superficie cultivada estándar37 .

3.3.2. Factores de cultivo

Según FAO (2006), El tipo de cultivo, la variedad y la etapa de desarrollo deben ser
considerados cuando se evalúa la evapotranspiración de cultivos que se desarrollan en
áreas grandes y bien manejadas. Las diferencias en resistencia a la transpiración, la altura
del cultivo, la rugosidad del cultivo, el reflejo, la cobertura del suelo y las características
radiculares del cultivo dan lugar a diferentes niveles de ET en diversos tipos de cultivos
aunque se encuentren bajo condiciones ambientales idénticas.

35
ALLEN Richard G. et al, 2006. Pág. 30
36
http://html.rincondelvago.com/evaporacion-y-evapotranspiracion.html
37
Burman, R.D., Jensen, M.E., y Allen, R.G. 1987. Los factores termodinámicos en la evapotranspiración. James, L.G.
y M.J. inglés (editores), Proc. Riego y Drenaje. La especificación. Conf., ASCE, Portland, Oregón, pág. 28-30.

33
INGENIERÍA Ingeniería Agronómica

La evapotranspiración del cultivo bajo condiciones estándar (ETc), se refiere a la demanda


evaporativa de la atmósfera sobre cultivos que crecen en áreas grandes bajo condiciones
óptimas de agua en el suelo, con características adecuadas tanto de manejo como
ambientales, y que alcanzan la producción potencial bajo las condiciones climáticas dadas
(FAO, 2006).

3.3.3. Manejo y condiciones ambientales

Los factores tales como salinidad o baja fertilidad del suelo, uso limitado de fertilizantes,
presencia de horizontes duros o impermeables del suelo, ausencia de control de
enfermedades y de parásitos y el mal manejo del suelo pueden limitar el desarrollo del
cultivo y reducir la evapotranspiración. Otros factores que se deben considerar al evaluar
la ET, son la cubierta del suelo, la densidad del cultivo y el contenido de agua del suelo.

Cuando se evalúa la tasa de ET, se debe considerar adicionalmente la gama de prácticas


locales de manejo que actúan sobre los factores climáticos y de cultivo afectando el
proceso de ET. Las prácticas del cultivo y el método de riego pueden alterar el microclima,
afectar las características del cultivo o afectar la capacidad de absorción de agua del suelo
y la superficie de cultivo (FAO, 2006).

12.4. ET calculada con datos meteorológicos

Debido a la dificultad de obtener mediciones de campo precisas (datos meteorológicos), la


ET se calcula comúnmente con datos meteorológicos. Una gran cantidad de ecuaciones
empíricas o semi-empíricas se han desarrollado para determinar la evapotranspiración del
cultivo o de referencia utilizando datos meteorológicos. Algunos de los métodos son
solamente válidos para condiciones climáticas y agronómicas específicas y no se pueden
aplicar bajo condiciones diferentes de las que fueron desarrolladas originalmente 38.

38
RIJTEMA, P.E., (1965). Analysis de la evapotranspiración. Agric. Res. Reporte. No. 69, Centre de Agric. Publ. y Doc.,
Wageningen, Países Bajos.

34
INGENIERÍA Ingeniería Agronómica

Numerosos investigadores han analizado el funcionamiento de los varios métodos del


cálculo para diversas localidades. Como resultado de una Consulta de expertos llevada a
cabo en mayo de 1990, el método FAO Penman-Monteith se recomienda actualmente
como el método estándar para la definición y el cálculo de la evapotranspiración de
referencia39.

12.5. Formulación de la ecuación de Penman-Monteith

12.5.1. Ecuación de Penman-Monteith

En 1948, Penman combinó el balance energético con el método de la transferencia de


masa y derivó una ecuación para calcular la evaporación de una superficie abierta de agua
a partir de datos climáticos estándar de horas sol, temperatura, humedad atmosférica y
velocidad de viento. Este método conocido como combinado fue desarrollado
posteriormente por muchos investigadores y ampliado a las superficies cultivadas por
medio de la introducción de factores de resistencia40.

12.5.2. Ecuación de FAO Penman-Monteith

Un panel de expertos e investigadores en riego fue organizado por la FAO en mayo de


1990, en colaboración con la Comisión Internacional para el Riego y Drenaje y con la
Organización Meteorológica Mundial, con el fin de revisar las metodologías previamente
propuestas por la FAO para el cálculo de los requerimientos de agua de los cultivos y para
elaborar recomendaciones sobre la revisión y la actualización de procedimientos a este
respecto (FAO, 2006).

El método reduce las imprecisiones del método anterior de FAO Penman y produce
globalmente valores más consistentes con datos reales de uso de agua de diversos
cultivos (FAO, 2006).

El método de FAO Penman-Monteith para estimar ETo, puede ser derivado de la ecuación
original de Penman-Monteith y las ecuaciones de la resistencia aerodinámica y superficial.

39
ALLEN Richard G. et al, 2006.
40
ALLEN et al, 1998.

35
INGENIERÍA Ingeniería Agronómica

13. LOCALIZACIÓN

Las estaciones meteorológicas de las que se tomaron los datos climáticos están ubicadas
en el departamento de la Paz (Viacha, Belén, Patacamaya), Oruro y Potosí, zonas
conocidas como andinas de Bolivia. En la figura 4 se muestra en un mapa la ubicación de
las estaciones meteorológicas.

En el cuadro 1 se muestra las estaciones meteorológicas y su ubicación, siendo estas


como ya se mencionaron lugares de altura en Bolivia y de las que se utilizaran los distintos
parámetros climáticos destinados a la investigación:

CUADRO 1
Características y localización de las estaciones meteorológicas de las que se utilizaron los
parámetros climáticos.

CARACTERISTICAS BELEN PATACAMAYA VIACHA ORURO POTOSÍ

Provincia Omasuyos Aroma Ingavi Cercado Tomas Frías


Departamento La Paz La Paz La Paz Oruro Potosí
Estación Belén Patacamaya Viacha Oruro Potosí (Los Pinos)

Latitud Sud: 16º 1' 0" 17º 12' 0" 16º 39' 0" 17º 58' 0" 19º 35' 0"

Longitud Oeste 68º 42' 0" 67º 55' 0" 68º 18' 0" 67º 7' 0" 65º 45' 0"

Altitud m/s/n/m 3820 3789 3850 3702 3950


Temperatura máxima
o 14.90 19.00 17.70 19.00 18.00
media ( C)
Temperatura mínima
o -0.40 1.00 -0.5 -0.13 -0.03
media ( C)
Humedad relativa
78.50 78.80 86.60 68.30 70.90
máxima media (%)
Humedad relativa
46.40 47.90 30.80 28.90 38.70
mínima media (%)
Velocidad media del
1.06 2.94 3.62 2.31 2.37
viento (m/s)
Radiación solar media
-2 -1 17.61 19.23 17.62 18.71 17.45
observada (MJ m d )
Radiación solar media
-2 -1 21.56 23.14 23.21 24.91 23.08
estimada (MJ m d )
Precipitación media
anual en el periodo 530.40 348.00 500.97 371.35 462.98
(mm)
Fuente: Elaboración propia

36
INGENIERÍA Ingeniería Agronómica

FIGURA 5.
Ubicación Geográfica de las estaciones meteorológicas utilizadas para el estudio.

Fuente: Elaboración propia

14. MATERIALES Y MÉTODOS


14.1. Materiales y equipos
14.1.1. Datos Meteorológicos

Se utilizo datos de localización; datos meteorológicos diarios de Temperatura Máxima y


Mínima, Radiación Solar, Humedad y Velocidad del Viento, de una serie de 5 años
calendario (2002 a 2006), correspondientes a cinco estaciones meteorológicas de zonas
de altura como lo son las estaciones de Belén, Viacha, Patacamaya, Oruro y Potosí (Los
Pinos); datos que fueron evaluados y validados por el Servicio Nacional de Meteorología e
Hidrología (SENAMHI).

37
INGENIERÍA Ingeniería Agronómica

14.1.2. Localización

Los datos de localización se refieren a la altura sobre el nivel del mar (m) de las zonas
para las que se determino la ETo, y su latitud (grados norte o sur), estos están
especificados en el Cuadro 1, presentado anteriormente. Estos datos son necesarios para
ajustar algunos parámetros climáticos al valor medio local de la presión atmosférica
(función de la elevación del sitio sobre nivel del mar) y para calcular la radiación
extraterrestre (Ra) y, en algunos casos, la duración máxima de la insolación (N). En los
procedimientos del cálculo para Ra y N, la latitud se expresa en radianes (grados
decimales π/180). Un valor positivo se utiliza para el hemisferio norte y un valor negativo
para el hemisferio sur.

14.1.3. Temperaturas

Para la aplicación de la fórmula FAO Penman-Monteith, se requiere información de


temperatura diaria (promedio) máxima y mínima en grados centígrados (°C). En los casos
en que solamente se tenga disponibilidad de temperaturas medias diarias, el cálculo de
ETo, todavía es válido pero se podría esperar una cierta subestimación de ETo, debido a
la relación no-linear de la presión de vapor de saturación en relación con la temperatura.
El uso de la temperatura media del aire en lugar de la temperatura máxima y mínima tiene
como resultado una presión de saturación más baja (es), y por lo tanto un déficit de presión
de vapor también más bajo (es - ea), resultante en una estimación más baja de la
evapotranspiración de referencia.

14.1.4. Humedad atmosférica

El valor de la presión real (promedio) diaria de vapor, (e a), en kilopascales (kPa), se


requiere como dato de entrada para la aplicación de la ecuación FAO Penman-Monteith.
En caso de que los valores de presión real de vapor no estén disponibles, estos pueden
ser derivados de la humedad relativa máxima y mínima (%), de los datos psicrométricos
(temperaturas de bulbo seco y mojado en °C) o de la temperatura del punto de rocío (°C).

38
INGENIERÍA Ingeniería Agronómica

14.1.5. Radiación solar

La radiación neta diaria (promedio), está expresada en megajoules por metro cuadrado
por día (MJ m-2 día-1). Estos datos no están disponibles comúnmente en forma directa pero
pueden derivarse de la radiación de onda corta (promedio), medida con un piranómetro o
de la duración real diaria (promedio) del sol brillante (horas por día) medida con el
heliógrafo.

14.1.6. Velocidad del viento

Para la aplicación de la ecuación FAO Penman-Monteith, se requiere contar con la


velocidad del viento (promedio) diaria en metros por segundo (m s-1) medida a 2 metros de
altura sobre el nivel del suelo. Es importante verificar la altura a la cual se mide la
velocidad del viento, pues velocidades del viento medidas a diversas alturas sobre la
superficie del suelo presentan con seguridad valores diferentes. En caso de no contar con
datos de velocidad de viento a 2 metros de altura, el procedimiento del cálculo para ajustar
la velocidad del viento a la altura estándar de 2 metros se presenta en la parte
metodológica.

14.2. Equipos

Para el procesamiento de la información de los parámetros climáticos (temperaturas,


radiación, humedad relativa y velocidad del viento, datos diarios de cinco estaciones, en
una serie de 5 años calendario, que corresponde 4 gestiones agrícolas, a partir del año
2002 al 2006), se utilizo una computadora, planillas en Excel, software (SISMET), material
recepcionado del SENAMHI, previamente validada.

Material de escritorio (papelería, bolígrafos, cuaderno, etc.); también fue utilizado para la
preparación del presente trabajo de investigación.

39
INGENIERÍA Ingeniería Agronómica

14.3. Metodología

Se utilizo la ecuación de FAO Penman-Monteith (Ecuación 1), para el cálculo de la


evapotranspiración de referencia diaria, que asume un cultivo de referencia como un
cultivo hipotético con una altura asumida de 0,12 m, con una resistencia superficial de 70 s
m-1 y un albedo de 0,23 y que representa a la evapotranspiración de una superficie
extensa de pasto verde de altura uniforme, creciendo activamente y adecuadamente
regado. (Figura 5).

El método planteado por la FAO Penmman-Monteith ha sido seleccionado como el método


por el cual la evapotranspiración de esta superficie de referencia (ETo), fue
inequívocamente determinada y con el cual se obtuvieron valores consistentes de ETo, en
las regiones altas de Bolivia mencionadas anteriormente.

FIGURA 6.
Características del cultivo de referencia

Fuente: ALLEN Richard G., 2006.

40
INGENIERÍA Ingeniería Agronómica

Ecuación. FAO Penmman-Monteith

(1)

Donde:
ETo Evapotranspiración de referencia (mm día-1),
Rn Radiación neta en la superficie del cultivo (MJ m -2 día-1)
Ra Radiación extraterrestre (mm día-1)
G Flujo del calor de suelo (MJ m-2 día-1)
T Temperatura media del aire a 2 m de altura (ºC)
u2 Velocidad del viento a 2 m de altura (m s-1)
es Presión de vapor de saturación (kPa)
ea Presión real de vapor (kPa)
es-ea Déficit de presión de vapor (kPa)
Δ Pendiente de la curva de presión de vapor (kPa ºC-1)
Constante psicrométrica (kPa °C-1)

La ecuación 1, utiliza datos climáticos de radiación solar, temperatura del aire, humedad y
velocidad del viento. Para asegurar la precisión del cálculo, Los datos climáticos en
especial la velocidad del viento debe ser medido a 2 metros de altura, sobre una superficie
extensa de pasto verde, cubriendo completamente el suelo y sin limitaciones de agua.

Para el cálculo de la ETo, con los datos meteorológicos completos se realizó el ajuste en
la altura de medición de la velocidad del viento a una altura de 2 m, que es la
recomendada por la FAO para evitar sobreestimaciones en los resultados.

La evapotranspiración de referencia (ETo), que se obtendrá proveerá un estándar de


comparación mediante el cual:

41
INGENIERÍA Ingeniería Agronómica

 Se podrá comparar la evapotranspiración en diversos períodos del año o en otras


regiones;
 Se podrá relacionar la evapotranspiración de otros cultivos.

Los datos hipotéticos utilizados, sobre los que se basa la ecuación FAO Penman-Monteith
serán tomados como un valor de comparación que homogeneizaran los resultados de tal
manera que los datos de diferentes zonas del mundo sean comparables entre sí.

14.3.1. FACTORES METEOROLÓGICOS QUE DETERMINAN LA


EVAPOTRANSPIRACÍON

Los factores meteorológicos que determinan la evapotranspiración son los componentes


del tiempo que proporcionan energía para la vaporización y extraen vapor de agua de una
superficie evaporante. Los principales parámetros meteorológicos que se deben
considerar se presentan a continuación.

14.3.1.1. TEMPERATURA DEL AIRE

El Servicio Nacional de Meteorología e Hidrología (SENAMHI), es la encarga de


recepcionar este tipo de información en Bolivia y de la cual se adquirió dicha información.
La temperatura del aire es medida con termómetros o termopares instalados en el abrigo.
Los termómetros de mínima y máxima registran la temperatura mínima y máxima del aire
en un período de 24 horas.

Debido a la relación no lineal de la humedad con la temperatura, ambas incluidas en la


ecuación FAO Penman-Monteith, la presión de vapor para cierto período se debe calcular
como la media entre la presión de vapor bajo la temperatura máxima y la presión de vapor
bajo la temperatura mínima del aire en ese periodo.

42
INGENIERÍA Ingeniería Agronómica

La temperatura máxima diaria del aire (Tmáx) y la temperatura mínima diaria del aire (Tmín)
son, respectivamente, la máxima y la mínima temperatura del aire observadas durante un
período de 24 horas, desde la medianoche. Tmáx y Tmín, por períodos más largos tales
como semanas, diez días o meses se obtienen dividiendo la suma de los valores diarios
respectivos por el número de días del período.

La temperatura media diaria del aire (Tmedia) se emplea en la ecuación FAO Penman-
Monteith solamente para calcular la pendiente de la curva de la presión de saturación de
vapor (Δ) y del impacto de la densidad media del aire (Pa), pues el efecto de las
variaciones de temperatura en el valor de estos parámetros climáticos es pequeño.

Para la estandarización de la Tmedia para períodos de 24 horas se define como el promedio


de las temperaturas máxima (Tmáx) y mínima diaria (Tmín) en lugar del promedio de las
mediciones horarias de temperatura.

(2)

Donde:
Tmedia : Temperatura media (ºC)
Tmáx : Temperatura Máx. (ºC)
Tmín : Temperatura Mínima (ºC)

14.3.1.2. HUMEDAD DEL AIRE

El contenido de agua del aire se puede expresar de diversas maneras. En


agrometeorología, la presión de vapor, la temperatura del punto de condensación (punto
de rocío) y la humedad relativa son expresiones comunes para indicar la humedad del
aire.

43
INGENIERÍA Ingeniería Agronómica

14.3.1.2.1. Presión de vapor

El vapor de agua es un gas y su presión contribuye a la presión atmosférica total. La


cantidad de vapor de agua en el aire se relaciona directamente con la presión parcial
ejercida por ese vapor de agua en el aire y esta es por lo tanto una medida directa del
contenido de vapor de agua del aire.

14.3.1.2.2. Humedad relativa

La humedad relativa (HR), expresa el grado de saturación del aire como el cociente entre
la presión real de vapor (ea) a una temperatura dada y la presión de saturación de vapor
(e°(T)) a la misma temperatura (T):

(3)

Donde:
ea : Presión real de vapor (kPa)
eº(T) : Presión de saturación T(kPa)

La humedad relativa es el cociente entre la cantidad de agua que el aire realmente


contiene a una determinada temperatura y la cantidad que podría contener si estuviera
saturado a la misma temperatura. Es adimensional y se expresa comúnmente como
porcentaje. Aunque la presión real de vapor puede ser relativamente constante a lo largo
del día, la humedad relativa fluctúa entre un máximo al amanecer y un mínimo a primeras
horas de la tarde. La variación de la humedad relativa se produce porque la presión de
saturación de vapor depende de la temperatura del aire. Como la temperatura del aire
cambia durante el día, la humedad relativa también cambia substancialmente.

14.3.1.2.3. Presión media de vapor de la saturación (es)

La presión de saturación de vapor puede ser calculada en función de la temperatura del


aire, pues depende de ella. La relación entre ambas variables se expresa como:

44
INGENIERÍA Ingeniería Agronómica

(4)

Donde
e°(T) : Presión de saturación de vapor a la temperatura del aire, T [kPa]
T : Temperatura del aire [°C]
exp [..] : 2,7183 (base del logaritmo natural) elevado a la potencia [..]

Los valores de la presión de saturación de vapor en función de la temperatura del aire se


presentan en el Anexo 4. Debido a la característica no-linear de la Ecuación 4, la presión
media de saturación de vapor para un día, fue calculada como el promedio de la presión
de saturación de vapor a la temperatura máxima media y la presión de saturación de vapor
a la temperatura mínima media del aire para ese período:

(5)

Donde:
es : Presión de saturación (kPa)
Tmax : Temperatura máxima (ºC)
Tmin : Temperatura mínima (ºC)

El uso de la temperatura media del aire en lugar de las temperaturas mínima y máxima da
lugar a subestimaciones de la presión media de saturación de vapor. El correspondiente
déficit de presión de vapor (un parámetro que expresa la energía evaporante de la
atmósfera), también será menor y el resultado será una cierta subestimación de la
evapotranspiración del cultivo de referencia.

Por lo tanto, es siempre recomendable el calcular la presión media de saturación de vapor


como el promedio de la presión de saturación de vapor a la temperatura máxima y a la
temperatura mínima del aire.

45
INGENIERÍA Ingeniería Agronómica

14.3.1.2.4. Pendiente de la curva de presión de saturación de vapor (Δ)

Para el cálculo de evapotranspiración, se requiere calcular la pendiente de la relación


entre la presión de saturación de vapor y la temperatura, (Δ). La pendiente de la curva a
una temperatura dada se da por:

(6)

Donde
Δ : Pendiente de la curva de la presión de saturación de vapor a la
temperatura del aire T [kPa °C-1]
T : Temperatura del aire [°C]
exp[.. ] : 2,7183 (base del logaritmo natural) elevado a la potencia [..]

Los valores de la pendiente Δ, para las estaciones meteorológicas y zonas en estudio


están en el Anexo 5. Para su aplicación en la ecuación FAO Penman-Monteith, la
pendiente de la curva de presión de vapor (Δ) se calcula usando la temperatura media del
aire (Ecuación 2).

14.3.1.2.5. Presión real de vapor (ea) derivada de datos de humedad relativa

La presión real de vapor (ea), se puede también derivar de la humedad relativa.


Dependiendo de la disponibilidad de datos de humedad, se pueden utilizar diversas
ecuaciones.

Disponibilidad de HRmáx y HRmín:

(7)

46
INGENIERÍA Ingeniería Agronómica

Donde:
ea : Presión real de vapor [kPa]
e°(Tmín) : Presión de saturación de vapor a la temperatura mínima diaria [kPa]
e°(Tmáx) : Presión de saturación de vapor a la temperatura máxima diaria [kPa]
HRmáx : Humedad relativa máxima [%].
HRmín : Humedad relativa mínima [%].

Para períodos de una semana, diez días o un mes, la HRmáx y la HRmín se obtienen
dividiendo la suma de los valores diarios por el número de días del período considerado.

14.3.1.2.6. Déficit de presión de vapor (es - ea)

El déficit de presión de vapor (es - ea), es la diferencia entre la presión de saturación de


vapor (es) y la presión real de vapor (ea) durante un determinado período. Para periodos
tales como una semana, diez días o un mes, e s se calcula de la Ecuación 5, usando el
promedio de Tmáx y Tmín en el período.

14.3.1.3. RADIACION
14.3.1.3.1. Radiación extraterrestre (Ra)

La radiación que choca a una superficie perpendicular a los rayos del sol en el extremo
superior de la atmósfera terrestre, se llama constante solar, y tiene un valor aproximado de
0,082 MJ m-2 min-1. La intensidad local de la radiación, sin embargo, está determinada por
el ángulo entre la dirección de los rayos solares y la superficie de la atmósfera. Este
ángulo cambia durante el día y es diferente en diversas latitudes y en diversas épocas del
año. La radiación solar recibida en la parte superior de la atmósfera terrestre sobre una
superficie horizontal se conoce como radiación (solar) extraterrestre, (Ra).

Si el sol se encuentra directamente encima de la cabeza, el ángulo de incidencia es cero y


la radiación extraterrestre es 0,082 MJ m-2 min-1. Así como las estaciones cambian, la
posición del sol, la longitud del día y la radiación extraterrestre también cambian. La
radiación extraterrestre es entonces una función de la latitud, la época del año y la hora
del día.

47
INGENIERÍA Ingeniería Agronómica

14.3.1.3.2. Radiación extraterrestre para periodos diarios (Ra)

La radiación extraterrestre, (Ra), para cada día del año y para diversas latitudes se puede
estimar a partir de la constante solar, la declinación solar y la época del año:

(8)

Donde
Ra : Radiación extraterrestre [MJ m-2 día-1]
Gsc : Constante solar = 0,082 MJ m-2 min-1,
dr : Distancia relativa inversa Tierra-Sol (Ecuación 10)
ωs : Ángulo de radiación a la puesta del sol (Ecuaciones 12) [rad]
: Latitud [rad] (Ecuación 9)

: Declinación solar (Ecuación 11) [rad].

La latitud ( ), expresada en radianes es positiva para el hemisferio norte y negativa para


el hemisferio sur. La conversión de grados decimales a radianes se da por:

(9)

La distancia relativa inversa Tierra-Sol, dr, y la declinación solar, δ, están dadas por:

(10)

(11)

Donde:
J : Es el número del día en el año entre 1 (1 de enero) y 365 (31 de diciembre).

48
INGENIERÍA Ingeniería Agronómica

Los valores de J para todos los días del año y una ecuación para estimar su valor se
presentan en el Anexo 6.

El ángulo de radiación a la hora de la puesta del sol, ωs, se da por:

(12)

14.3.1.3.3. Duración solar relativa (n/N)

La duración solar relativa (n/N), es otro cociente que expresa la nubosidad atmosférica. Es
el cociente de la duración real de la insolación, (n), y la duración máxima posible de la
insolación o de la luz del día (N).

Bajo condiciones de ausencia de nubosidad la duración real de la insolación es igual a la


insolación máxima posible (n = N) y el cociente es uno, mientras que en días muy
nublados (n) y el cociente puede tomar el valor de cero. En ausencia de datos directos de
Rs, la duración solar relativa, (n/N), se utiliza para derivar la radiación solar de la radiación
extraterrestre.

En la mayoría de las estaciones solo se realiza observaciones del cielo (lo sucedido en las
estaciones de las cuales se tomaron los datos de nubosidad), la nubosidad se expresa o
bien en decimos (0 a 10) o en octas (0 a 8).

Las observaciones están aproximadamente:


- Cielo despejado = 1 octa,
- Nubosidad parcial = 3 octas,
- Nubes = 6 octas,
- Cielo cubierto = 8 octas.

Las observaciones visuales de nubosidad se convirtieron en valores equivalentes de n/N.


con las relaciones que se presentan en el Cuadro 2.

49
INGENIERÍA Ingeniería Agronómica

CUADRO 2.
Nubosidad expresada en valores equivalentes de n/N

Nubosidad Nubosidad
0 .2 .4 .6 .8 0 .2 .4 .6 .8
(décimas) (décimas)*

0 .95 .90 .90 .90 .90 0 .90 .90 .90 .85 .85
1 .85 .85 .85 .85 .80 1 .85 .85 .80 .80 .80
2 .80 .80 .80 .75 .45 2 .75 .75 .75 .70 .70
3 .75 .70 .70 .70 .65 3 .65 .65 .65 .60 .60
4 .65 .65 .60 .60 .60 4 .55 .55 .50 .50 .45
5 .55 .55 .50 .50 .50 5 .45 .40 .40 .35 .35
6 .50 .45 .45 .45 .40 6 .30 .30 .25 .25 .20
7 .40 .35 .35 .35 .30 7 .15 .15 - - -
8 .30 .25 .25 .25 .20
9 .15 .15 .15 .15 -
* Octas: escala de 9 puntos en la cual 8 corresponde al cielo cubierto y 0 representa cielo despejado41.

14.3.1.3.4. Duración máxima de la insolación (N)

La duración máxima de la insolación N, está dada por:

(13)

Donde:
ωs : Es el ángulo de radiación a la hora de la puesta del sol calculado con
Ecuaciones 12.

14.3.1.3.5. Radiación solar o de onda corta (Rs)

Cuando la radiación atraviesa la atmósfera para llegar a la superficie terrestre, parte de


ella se dispersa, refleja o absorbe por los gases, las nubes y el polvo atmosféricos.

41
DOORENBOS J., PRUITT W. O., et.al. Las necesidades de agua de los cultivos, Estudio FAO Riego y
Drenaje 24 (Organización de las Naciones Unidas para la Agricultura y la alimentación), Roma 1984.

50
INGENIERÍA Ingeniería Agronómica

La cantidad de radiación que llega a un plano horizontal en la superficie terrestre se


conoce como radiación solar, (Rs). Como el sol emite energía por medio de ondas
electromagnéticas caracterizadas por longitudes de onda cortas, la radiación solar también
se conoce como radiación de onda corta.

Al no contar con mediciones directas de radiación solar (Rs), ésta se calculo a través de la
aplicación de la fórmula de Angstrom que relaciona la radiación solar con la radiación
extraterrestre y la duración relativa de la insolación:

(14)
Donde:
Rs: Radiación solar o de onda corta [MJ m-2 día-1],
n: Duración real de la insolación [horas],
N: Duración máxima posible de la insolación [horas],
n/N: Duración relativa de la insolación[-],
Ra: Radiación extraterrestre [MJ m-2 día-1],
as: Constante de regresión, que expresa la fracción radiación extraterrestre que
llega a la tierra en días muy nublados (n = 0),
as+ bs: Fracción de la radiación extraterrestre que llega a la tierra en días
despejados (n = N).

14.3.1.3.6. Radiación solar en un día despejado (Rso)

La radiación en días despejados (Rso), cuando n = N, puede ser calculada de la siguiente


forma:

Para localidades cerca al nivel del mar o cuando se dispone de valores calibrados
de as y bs:

(15)

51
INGENIERÍA Ingeniería Agronómica

Donde
Rso : Radiación solar en un día despejado [MJ m-2 día-1],
as+bs : Fracción de la radiación extraterrestre que llega a la superficie
terrestre en días de cielo despejado.

En casos en que los valores calibrados de as y bs no estén disponibles

(16)

Donde
z : Elevación de la estación sobre el nivel del mar [m].

14.3.1.3.7. Radiación neta solar o de onda corta (Rns)

La radiación neta de onda corta (Rns), resultante del equilibrio entre la radiación solar
entrante y la reflejada está dada por:

(17)

Donde
Rns : Radiación neta solar o de onda corta [MJ m-2 día-1],
α : Albedo o coeficiente de reflexión del cultivo, que es 0,23 para el cultivo
hipotético de referencia [adimensional],
Rs : Radiación solar entrante [MJ m-2 día-1].
Rns : En la Ecuación 38 está expresada en MJ m-2 día-1.

14.3.1.3.8. Radiación neta de onda larga (Rnl)

La cantidad de emisión de energía de onda larga es proporcional a la temperatura


absoluta de la superficie elevada a la cuarta potencia. Esta relación se expresa
cuantitativamente por la ley de Stefan-Boltzmann. Se debe tener en cuenta que el flujo de
energía neta que sale de la superficie terrestre es menor que la calculada y dada por la ley
de Stefan-Boltzmann debido a la absorción y radiación devuelta del cielo.

52
INGENIERÍA Ingeniería Agronómica

El vapor de agua, las nubes, el dióxido de carbono y el polvo absorben y emiten radiación
de onda larga. Por ello se deben conocer sus concentraciones para determinar el flujo
saliente neto. Como la humedad y la nubosidad tienen un papel importante, la ley de
Stefan-Boltzmann, se corrige por estos dos factores cuando se estima el flujo saliente neto
de la radiación de onda larga. De tal modo que se asume que las concentraciones de los
otros factores de absorción de radiación son constantes:

(18)

Donde

Rnl : Radiación neta de onda larga [MJ m-2 día-1],


Σ : Constante de Stefan-Boltzmann [4,903 x 10-9 MJ K-4 m-2 día-1],
Tmax,K : Temperatura máxima absoluta durante un periodo de 24 horas
[K = °C + 273,16],
Tmin,K : Temperatura mínima absoluta durante un periodo de 24 horas
[K = °C + 273,16],
ea : Presión de vapor real [kPa],
Rs/Rso : Radiación relativa de onda corta (valores ≤ 1,0),
Rs : Radiación solar medida o calculada (Ecuación 14) [MJ m-2 día-1],
Rso : Radiación en un día despejado (Ecuación 15 o 16) [MJ m-2 día-1].

14.3.1.3.9. Radiación neta (Rn)

La radiación neta (Rn) es la diferencia entre la radiación neta de onda corta (R ns) y la
radiación neta de onda larga (Rnl):

(19)

53
INGENIERÍA Ingeniería Agronómica

Donde:
Rn : Radiación neta
Rns : Radiación neta solar o de onda corta
Rnl : Radiación neta de onda larga

14.3.1.4. VELOCIDAD DEL VIENTO

El viento se caracteriza por su dirección y su velocidad. La dirección del viento se refiere a


la dirección de la cual el viento está soplando. Para el cómputo de la evapotranspiración,
la velocidad del viento es una variable importante. Como la velocidad del viento en una
localidad dada varía con el tiempo, es necesario expresarla como el promedio sobre un
intervalo determinado de tiempo. La velocidad del viento se mide en metros por segundo
(m s-1) o kilómetros por día (km día-1).

La velocidad del viento se mide con los anemómetros. Los anemómetros usados
comúnmente en las estaciones meteorológicas se componen de cazoletas o propulsores
que giran sobre un eje gracias a la fuerza del viento. El conteo del número de revoluciones
en un periodo dado, permite determinar la velocidad promedio del viento en dicho período.

14.3.1.4.1. Relación del viento con la altura

La velocidad del viento a diversas alturas sobre la superficie del suelo tienen valores
diferentes. La fricción superficial tiende a reducir la velocidad del viento que atraviesa la
superficie. La velocidad del viento es menor cerca de la superficie y aumenta con altura.

Por esta razón los anemómetros se colocan en una altura estándar elegida, 10 metros en
meteorología y 2 ó 3 metros en agrometeorología. Para el cálculo de la
evapotranspiración, se requiere la velocidad del viento medida a 2 metros sobre la
superficie.

54
INGENIERÍA Ingeniería Agronómica

Para ajustar los datos de velocidad del viento obtenido de instrumentos situados a
elevaciones diferentes a la altura estándar de 2m, se puede usar una relación logarítmica:

(20)

Donde
u2 : Velocidad del viento a 2 m sobre la superficie [m s-1],
uz : Velocidad del viento medida a z m sobre la superficie [m s-1],
z : Altura de medición sobre la superficie [m].

14.4. PARAMETROS ATMOSFERICOS

El efecto de los principales componentes del tiempo sobre la evapotranspiración se puede


determinó con la ayuda de las ecuaciones:

14.4.1. Presión atmosférica (P)

La presión atmosférica (P), es la presión ejercida por el peso de la atmósfera terrestre. La


evaporación en altitudes elevadas ocurre en parte gracias a la baja presión atmosférica
que se expresa con la constante psicrométrica. Este efecto es, sin embargo, pequeño y en
los procedimientos del cálculo, el valor medio para una localidad es suficiente.

Para calcular P, puede emplearse una simplificación de la ley de los gases ideales, a una
temperatura atmosférica estándar de 20°C:

(21)

Donde

P : Presión atmosférica [kPa]


Z : Elevación sobre el nivel del mar [m]

55
INGENIERÍA Ingeniería Agronómica

14.4.2. Calor latente de vaporización (λ)

El calor latente de vaporización (λ), expresa la energía requerida para cambiar una masa

de unidad de agua líquida a vapor de agua bajo presión y temperatura constantes. El valor
del calor latente de vaporización varía en función de la temperatura. Cuanto más elevada
sea la temperatura, menos energía será requerida.

Como λ, varía levemente dentro de rangos de temperaturas normales, se considera un

valor constante de 2,45 MJ kg-1 para la simplificación de la ecuación de FAO Penman-


Monteith. Este valor corresponde al calor latente de vaporización a una temperatura del
aire de alrededor de 20 °C.

14.4.3. Constante Psicrométrica (γ)

La constante psicrométrica, γ, se calcula por:

(22)

Donde

γ : Constante psicrométrica [ kPa °C-1],

P : Presión atmosférica [kPa],

λ : Calor latente de vaporización, 2,45 [ MJ kg-1],

Cp : Calor específico a presión constante, 1,013 x 10-3 [ MJ kg-1 °C-1],


ε : Cociente del peso molecular de vapor de agua /aire seco = 0,622.

El calor específico es la cantidad de energía requerida para aumentar la temperatura de


una unidad de masa de aire en 1°C, a presión constante. Su valor depende de la
composición del aire, fundamentalmente de su humedad.

56
INGENIERÍA Ingeniería Agronómica

Para condiciones atmosféricas medias se puede utilizar el valor Cp = 1,013 x 10-3 MJ kg-1
°C-1. Dado que se aplica un valor medio de presión atmosférica en cada localidad, la
constante psicrométrica también se mantiene constante para cada localidad.

14.5. ESTIMACIÓN DE DATOS CLIMÁTICOS FALTANTES

Si no se dispone de alguno de los datos climáticos requeridos o estos no pueden ser


calculados, se recomienda la estimación de estos datos con uno de los procedimientos
siguientes para con ello aplicar la ecuación de FAO Penman-Monteith, para el cálculo de
ETo.

14.5.1. Estimación de datos faltantes de humedad

En caso en que los datos de humedad no estén disponibles o sean de calidad


cuestionable, se puede llevar adelante una estimación de la presión real de vapor, e a, si se
asume que la temperatura del punto de rocío (T rocío) es similar a la temperatura mínima
diaria (Tmín).

Se asume implícitamente que a la salida del sol, cuando la temperatura del aire está
alrededor de Tmín, el aire está casi saturado con vapor de agua y que la humedad relativa
es de alrededor de 100%. Si Tmín se utiliza para representar Trocío entonces:

(23)

14.5.2. Estimación de los datos de radiación

Existen pocos instrumentos de medición de la radiación neta instalados en estaciones


agrometeorológicas. A falta de mediciones directas, la radiación de onda larga y la
radiación neta pueden ser derivadas de otros parámetros meteorológicos más comunes

57
INGENIERÍA Ingeniería Agronómica

p.e. horas de insolación, temperatura del aire y presión de vapor. En casos en que la
radiación solar no sea medida, esta puede ser estimada a partir de las horas de insolación.
Sin embargo, en casos en que no están disponibles las horas diarias de insolación (n), los
datos de radiación no pueden calcularse con los procedimientos de cálculo presentados
previamente.

14.5.3. Radiación solar derivada de las diferencias térmicas

La diferencia entre la temperatura máxima y mínima del aire se relaciona con el grado de
nubosidad en una localidad. Condiciones despejadas determinan altas temperaturas
durante el día (Tmáx), porque la atmósfera es transparente a la radiación solar entrante.

También se presentan bajas temperaturas durante la noche (T mín), porque menos


radiación saliente de onda larga es absorbida por la atmósfera. Por otra parte, en
condiciones cubiertas (nubladas), el valor de Tmáx, es relativamente más pequeño porque
una parte significativa de la radiación solar entrante nunca alcanza la superficie de la tierra
y es absorbida y reflejada por las nubes. Similarmente Tmin, será relativamente más alta
porque la cobertura de la nubosidad actúa como cobertor y disminuye la radiación saliente
neta de onda larga. Por lo tanto, la diferencia entre temperatura del aire máxima y mínima
(Tmax -Tmin), se puede utilizar como indicador de la fracción de la radiación extraterrestre
que alcanza la superficie de la tierra.

Este principio ha sido utilizado por Hargreaves y Samani, para desarrollar estimaciones de
ETo usando solamente datos de temperatura del aire.

La ecuación de radiación de Hargreaves, ajustada y validada en varias estaciones


meteorológicas en una amplia variedad de condiciones climáticas, está dada por:

(24)

58
INGENIERÍA Ingeniería Agronómica

Donde:
Ra : Radiación extraterrestre [MJ m-2 d-1],
Tmáx : Temperatura máxima del aire [ºC],
Tmín : Temperatura mínima del aire [ºC],
kRs : Coeficiente de ajuste (0,16 .. 0,19) [ºC-0,5].

La raíz cuadrada de la diferencia de temperaturas está muy relacionada con la radiación


solar diaria en una localidad dada. Los coeficientes de ajuste, kRs, son empíricos y se
diferencian para zonas del interior y las regiones costeras:

 Para las localidades «interiores», en donde la masa de tierra domina y las


masas de aire no están influenciadas fuertemente por un cuerpo grande del
agua, kRs =0,16;
 Para las localizaciones «costeras», situadas en la costa o cerca de una masa
grande de tierra y donde las masas de aire están influenciadas por un cuerpo
de agua cercano, kRs = 0,19.
 Recientes investigaciones para zonas de altura indican que en la practica estos
los resultados indican que el valor kRs = 0.14 se puede usar en todos lugares de
altura.

14.6. Ecuación alternativa para el cálculo de la ETo

Como una opción alternativa y comparativa, la ETo, se calculo usando la ecuación de


Hargreaves:

(25)
Donde:
ETo : Evapotranspiración de referencia (mm día-1),
Tmedia : Temperatura media (ºC),
Tmáx. : Temperatura media (ºC),
Tmín. : Temperatura mínima (ºC),
Ra : Radiación extraterrestre [MJ m-2 d-1].

59
INGENIERÍA Ingeniería Agronómica

Donde todos los parámetros han sido previamente definidos. Las unidades de ETo y Ra
(Radicación extraterrestre) son mm día-1. Antes de aplicar la Ecuación 23 en cada nueva
región, se debe verificar la validez de sus resultados comparándolos con estimaciones de
la ecuación de FAO Penman-Monteith (Ecuación 1) en estaciones meteorológicas donde
se dispone de datos de radiación solar, temperatura del aire, humedad, y velocidad del
viento.

También es importante recalcar que la Ecuación 23, tiene una tendencia a subestimar los
valores de ETo, bajo condiciones de viento fuerte (u2 > 3 m s-1) y a sobreestimar la ETo
bajo condiciones de elevada humedad relativa.

14.7. ANÁLISIS ESTADISTICO

El objetivo principal por el cual se utilizo el análisis de regresión fue proyectar el valor de la
variable dependiente conociendo o suponiendo valores para la variable independiente. La
confiabilidad de las proyecciones está dada por la confiabilidad de la ecuación, la cual se
mide a través del coeficiente de determinación y de los errores de los coeficientes de
regresión.

14.7.1. Coeficiente de Determinación:

El coeficiente de determinación (R2) nos muestra cuanto se ajusta la línea de regresión a


los datos y la misma está dada por:

(26)

Donde:

Oi y P i :(i = 1, 2, ..., n) son pares de observación


OyP : Son los respectivos valores medios

60
INGENIERÍA Ingeniería Agronómica

Como puede observarse, el coeficiente de determinación es la proporción de la variable


dependiente explicada por la variable independiente y por lo tanto está entre 0 y 1. Es
decir: 0 £ R2 £ 1.

A medida que el R2 se acerca a 1, la ecuación de regresión es más confiable y entre más


cercano esté el R2 de cero, la ecuación es menos confiable.

14.7.2. Coeficiente de Regresión

El coeficiente de regresión (b), está expresado en las mismas unidades de Y por cada
unidad de X. Indica el número de unidades en que varía Y cuando se produce un cambio,
en una unidad, en X (pendiente de la recta de regresión). Indica el número de unidades en
que se modifica la variable dependiente “Y” por efecto del cambio de la variable
independiente “X” o viceversa en una unidad de medida.

(27)

Donde:

Oi y P i : (i = 1, 2, ..., n) son pares de observación


OyP : Son los respectivos valores medios

El coeficiente de regresión puede ser: Positivo, Negativo y Nulo. Es positivo cuando las
variaciones de la variable independiente X son directamente proporcionales a las
variaciones de la variable dependiente “Y”. Es negativo, cuando las variaciones de la
variable independiente “X” son inversamente proporcionales a las variaciones de las
variables dependientes “Y”. Es nulo o cero, cuando entre las variables dependientes “Y” e
independientes “X” no existen relación alguna.

61
INGENIERÍA Ingeniería Agronómica

14.7.3. Error Cuadrado Medio (RMSE)

(28)

Donde:
n : Es el número de datos de las series Oi (medida) y Pi (simulada).

14.7.4. Promedio Absoluto del Error (AAE)

(29)

Donde:
Oi y P i :(i = 1, 2, ..., n) son pares de observación

14.7.5. Promedio Relativo del error (ARE)

(30)

Donde:
Oi y P i :(i = 1, 2, ..., n) son pares de observación

62
INGENIERÍA Ingeniería Agronómica

15. RESULTADOS Y DISCUSIONES

15.1. Estimación de la humedad a partir de la temperatura mínima

La estimación de la humedad para aquellas estaciones en las que estos datos no estén
disponibles se realizo a partir de la presión actual de vapor (ea), la que se calculó a partir
de la temperatura mínima (Tmín), por medio de la Ecuación 5. Asumiendo que a la salida
del sol, cuando la temperatura del aire está alrededor de Tmín, el aire está casi saturado
con vapor de agua y que la humedad relativa es de alrededor de 100%.

Una vez realizada la regresión lineal y obteniendo el coeficiente de determinación (R2), y el


coeficiente de regresión (b) mediante las ecuaciones 26 y 27, se comparó la ETo,
calculada mediante la Ecuación FAO-PM con datos observados o completos, con la
calculada a partir de la Tmín (ETo-Tmín) o datos limitados de humedad (figura 7).

FIGURA 7.
Comparación de ETo – PM con datos completos y ETo PM con valores estimados de humedad (HR) a
partir de datos diarios de temperatura mínima y de la presión actual de vapor ea, para estaciones de
Belén, Patacamaya, Viacha, Oruro y Potosí respectivamente.

63
INGENIERÍA Ingeniería Agronómica

Los valores de HR estimados a partir de la Ecuación 23 fueron similares a los observados


en las cinco localidades, lo que produjo como consecuencia valores de ETo simulados
similares a los valores de ETo – PM con datos completos, lo que es un indicador de que la
temperatura mínima es muy próxima a la temperatura del punto de roció, especialmente
en Belén.

64
INGENIERÍA Ingeniería Agronómica

TABLA 1.
Indicadores de desempeño para la comparación de resultados de ETo calculados con datos
completos y con datos faltantes de humedad relativa, para valores diarios.

Datos faltantes de HR
ESTACIÓN
RMSE* AAE** ARE***
METEOROLÓGICA
(mm d-1) (mm d-1) (%)
Belén 0.0700 0.0437 1.5671
Patacamaya 0.6229 0.2686 7.9882
Viacha 0.2393 0.1986 5.4555
Oruro 0.2333 0.1576 3.8722
Potosí 0.2542 0.1645 4.5467
* Error Cuadrado Medio, ** Promedio Absoluto del Error, *** Promedio Relativo del error

Los valores de RMSE variaron de 0,0700 a 0,6229 mm d-1 para Belén y Patacamaya,
respectivamente, dato que nos muestra que el modelo se comporta de manera muy buena
en el caso de Belén y buena en el caso de Patacamaya; AAE de 0,0437 a 0,2686 mm d-1
para Belén y Patacamaya, respectivamente; y ARE de 1,5671 a 7.9882 % porcentajes que
muestran la eficiencia y confiabilidad de los resultados para esta variable (HR), habiendo
obtenido Belén el valor más bajo y Patacamaya el valor más alto.

En tanto que el coeficiente de determinación R2 varió entre 0,99 a 0,85, para Belén y
Viacha, respectivamente reafirmando la poca variabilidad de las temperaturas mínimas
con relación a la temperatura del punto de roció.

15.2. Estimación de la radiación solar a partir de datos diarios de temperatura


máxima y mínima para la estación meteorológica de Belén

Cuando solamente se dispone de datos de temperatura, Rs se determinado a partir de la


Ecuación (24). Se asumió un valor para kRS igual a 0,16 ya que las estaciones
Meteorológicas no están influenciadas fuertemente por un cuerpo grande del agua.

65
INGENIERÍA Ingeniería Agronómica

Posteriormente se efectuó la regresión entre ETo-PM con datos completos y con datos
faltantes de Rs (figura 8).

FIGURA 8.
Comparación de ETo – PM con datos completos y ETo PM con valores estimados de Radiación Solar
(Rs) a partir de datos diarios de temperatura máxima y mínima, para estaciones de Belén,
Patacamaya, Viacha, Oruro y Potosí respectivamente.

66
INGENIERÍA Ingeniería Agronómica

En general, los valores de Rs estimados a partir de la Ecuación 24 fueron mayores que los
observados, lo que produjo como consecuencia valores de ETo simulados mayores que
los valores de ETo – PM con datos completos, debido probablemente a la gran diferencia
entre la temperatura máxima y mínima (siendo está en última en muchos casos negativa),
debido a la altura es una de las variables que interviene en el cálculo de la radiación solar.

Tabla 2.
Indicadores de desempeño para la comparación de resultados de ETo calculados con datos
completos y con datos faltantes de Radiación Solar, para valores diarios.

Datos faltantes de Rs
ESTACIÓN
RMSE* AAE** ARE***
METEOROLÓGICA
(mm d-1) (mm d-1) (%)
Belén 0.4430 0.3358 12.4542
Patacamaya 0.4574 0.3329 10.3533
Viacha 0.5484 0.4866 13.3368
Oruro 0.8093 0.6838 18.9276
Potosí 0.3952 0.2784 8.1232
* Error Cuadrado Medio, ** Promedio Absoluto del Error, *** Promedio Relativo del error

67
INGENIERÍA Ingeniería Agronómica

Los valores de RMSE variaron de 0,3952 a 0,8093 mm d-1 para Potosí y Oruro,
respectivamente, lo que indica que el modelo no tuvo un buen desempeño debido a la
variabilidad de los datos por la diferencia de altitudes que existe entre estaciones
meteorológicas; un AAE de 0,2784 a 0,6838 mm d-1 para Potosí y Oruro, respectivamente,
mostrando así la precisión de los cálculos para dicha variable; y una variación del 8,1232
al 18,9276% del ARE hacen que los datos no se distribuyan con simetría alrededor del
valor promedio obtenido con todos los datos completos para Potosí y Oruro,
respectivamente.

En tanto que el coeficiente de determinación R2, varió entre 0,74 a 0,86, para Belén y
Potosí, respectivamente, esto posiblemente a que la estación climática de Belén está
influenciada por una masa de agua como ser el lago Titicaca.

15.3. Validación del procedimiento cuando no se cuenta con datos de velocidad del
viento

Los datos de velocidad del viento (U2), no siempre están disponibles en las estaciones
agrometeorológicas lo que dificulta el cálculo de la ETo. En el presente estudio se asumió
el valor medio por defecto a nivel mundial, es decir, 2 m s-1. Los resultados se muestran en
la Figura 9.

Los valores de ETo – PM determinados empleando un valor para U2 igual a valor medio
mundial, fueron sobrestimados para Belén y subestimados para el caso de Viacha, puesto
que en la primera localidad la media de U2 para el periodo de tiempo considerado es 1,01
m s-1, en tanto que en la segunda es 3,62 m s -1. Lo que indica que están 1.00 abajo y 1.62
arriba de la media mundial propuesta (2 m s-1).

68
INGENIERÍA Ingeniería Agronómica

FIGURA 9.
-1
Comparación de ETo – PM con datos completos y ETo PM calculado a partir de U2=2 m s , para
estaciones de Belén, Patacamaya, Viacha, Oruro y Potosí respectivamente:

69
INGENIERÍA Ingeniería Agronómica

Los valores de ETo – PM determinados empleando un valor para U2 igual a valor medio
mundial, fueron sobrestimados para Belén y subestimados para el caso de Viacha, puesto
que en la primera localidad la media de U2 para el periodo de tiempo considerado es 1,01
m s-1, en tanto que en la segunda es 3,62 m s-1.

TABLA 3.
Indicadores de desempeño para la comparación de resultados de ETo calculados con datos
completos y con datos faltantes de Viento, para valores diarios.

Datos faltantes de U2
ESTACIÓN
RMSE* AAE** ARE***
METEOROLÓGICA
(mm d-1) (mm d-1) (%)
Belén 0.3355 0.2624 10.1922
Patacamaya 0.3923 0.2852 8.3058
Viacha 0.2742 0.2492 6.6148
Oruro 0.4424 0.3499 9.8079
Potosí 0.3263 0.2277 6.2966
* Error Cuadrado Medio, ** Promedio Absoluto del Error, *** Promedio Relativo del error

70
INGENIERÍA Ingeniería Agronómica

Los valores de RMSE, variaron de 0,2742 a 0,4424 mm d-1 para Viacha y Oruro,
respectivamente lo que indica que el modelo tuvo un buen desempeño tomando en cuenta
el promedio de la velocidad del viento propuesto por FAO - PM; un AAE de 0,2277 a
0,3499 mm d-1 para Potosí y Oruro respectivamente, mostrando así la precisión de los
cálculos para dicha variable; y una variación del 6,2966 al 10,1922% del ARE hacen que
los datos no se distribuyan con simetría alrededor del valor promedio recomendado por el
método FAO-PM ya que se encontró una variación de +/- 1 m s-1 en estaciones como
Potosí y Belén, respectivamente.

En tanto que el coeficiente de determinación R2, varió entre 0,80 a 0,96, para Oruro y
Viacha, respectivamente. Mostrándose así la confiabilidad de los datos y resultados
expuestos en los caculos.

15.4. Estimación de la ETo cuando no se cuenta con datos de HR% y Rs.

El procedimiento empleado calcula la ETo, mediante la Ecuación de FAO-PM con datos


faltantes de HR% y Rs. Donde Rs fue estimada a partir de la Tmax y Tmin (Ecuación 24) y
ea estimada a partir de la Tmin (Ecuación 23). Por tanto las variables serán la velocidad del
viento y la temperatura.

De igual forma se realizo una regresión lineal para realizar la comparación entre los
cálculos de ETo con datos completos y los calculados con datos faltantes de Humedad y
Radiación solar, en la figura 10 se muestran los resultados de estos cálculos.

71
INGENIERÍA Ingeniería Agronómica

FIGURA 10.
Comparación de ETo – PM con datos completos y ETo PM calculado con valores estimados de
Radiación Solar (Rs) a partir de datos diarios de temperatura máxima y mínima y con valores
estimados de humedad (HR) a partir de datos diarios de temperatura mínima para estaciones de
Belén, Patacamaya, Viacha, Oruro y Potosí respectivamente.

72
INGENIERÍA Ingeniería Agronómica

TABLA 4.
Indicadores de desempeño para la comparación de resultados de ETo calculados con datos
completos y con datos faltantes de Humedad y Radiación solar, para valores diarios.

Datos faltantes de HR y Rs
ESTACIÓN
RMSE* AAE** ARE***
METEOROLÓGICA
(mm d-1) (mm d-1) (%)
Belén 0.4691 0.3616 13.3926
Patacamaya 0.7261 0.5958 18.3963
Viacha 0.5248 0.4319 12.4743
Oruro 0.7013 0.5530 15.9264
Potosí 0.5061 0.3673 10.7953
* Error Cuadrado Medio, ** Promedio Absoluto del Error, *** Promedio Relativo del error

Tomando en cuenta la falta de dos datos como lo son la Rs y la HR, los valores de RMSE,
variaron de 0,4691 a 0,7261 mm d-1 para Belén y Patacamaya respectivamente,
ratificando lo obtenido para cada variable independientemente y demostrando que el
modelo tiene un comportamiento aceptable cuanto estas variables están ausentes; el AAE
de 0,3616 a 0,5958 mm d-1 para Belén y Patacamaya respectivamente, demostrando así la
precisión de los cálculos cuando no se cuentan con estas dos variables; y un ARE de

73
INGENIERÍA Ingeniería Agronómica

10,7953 a 18,3963 % para Potosí y Patacamaya respectivamente, porcentaje que hace


que los resultados no se distribuyan con simetría en relación con los valores obtenidos
cuando se cuentan con todos los datos requeridos.

En tanto que el coeficiente de determinación R2 varió entre 0,66 a 0,78, para Viacha y
Oruro, respectivamente, brindando la confiabilidad requerida para el uso de estos datos en
sus distintas aplicaciones.

15.5. Estimación de la ETo cuando no se cuenta con datos de HR% y U2.

Cálculo realizado a partir de la Tmin para la Humedad y un valor medio por defecto a nivel
mundial de 2 m s-1, para la velocidad del viento. Por lo que las variables serán la radiación
solar y la temperatura.

En la figura 11, se muestra los resultados de la regresión lineal conjuntamente la


comparación entre los cálculos de ETo con datos completos y los calculados con datos
faltantes de Humedad y velocidad del viento.

FIGURA 11.
Comparación de ETo – PM con datos completos y ETo PM calculado con valores estimados de
-1
Humedad (HR) a partir de datos diarios de temperatura mínima y velocidad del viento = 2 m s
(promedio de 2000 estaciones meteorológicas en todo el mundo), para estaciones de Belén,
Patacamaya, Viacha, Oruro y Potosí respectivamente.

74
INGENIERÍA Ingeniería Agronómica

75
INGENIERÍA Ingeniería Agronómica

TABLA 5.
Indicadores de desempeño para la comparación de resultados de ETo calculados con datos
completos y con datos faltantes de Humedad y Velocidad del viento, para valores diarios.

Datos faltantes de HR y U2
ESTACIÓN
RMSE* AAE** ARE***
METEOROLÓGICA
(mm d-1) (mm d-1) (%)
Belén 0.3822 0.3205 12.2592
Patacamaya 0.5434 0.4343 13.3754
Viacha 0.4021 0.3501 9.2161
Oruro 0.4939 0.3702 9.7919
Potosí 0.4596 0.3607 10.0894
* Error Cuadrado Medio, ** Promedio Absoluto del Error, *** Promedio Relativo del error

Cuanto los datos faltantes son la HR y la U2 los valores de RMSE variaron de 0,3822 a
0,5434 mm d-1 para Belén y Patacamaya respectivamente, indicando que el modelo
aplicado se adecua y comporta de manera aceptable a la falta de estas dos variables; el
AAE de 0,3205 a 0,4343 mm d-1 para Belén y Patacamaya respectivamente, resultados
que muestran la precisión de los datos obtenidos a la falta de estas dos variables; y un
ARE de 9,2161 a 13,3754 % para Viacha y Patacamaya respectivamente, resultados que
muestran en valores porcentuales que los datos no se comportan con simetría a
comparación de los datos obtenidos cuando se cuenta con todos los datos necesarios.

En tanto que el coeficiente de determinación (R2), varió entre 0,65 a 0,87, para
Patacamaya y Belén, respectivamente mostrando así la confiabilidad de los cálculos y los
resultados obtenidos mediante el método FAO-PM.

15.6. Estimación de la ETo cuando no se cuenta con datos de Rs y U2.

Donde Rs fue estimada a partir de la Tmax y Tmin y un valor medio por defecto a nivel
mundial de 2 m s-1, para la velocidad del viento. Por lo que las variables serán la radiación
solar y la velocidad del viento.

76
INGENIERÍA Ingeniería Agronómica

La figura 12, nos muestra los resultados de la regresión lineal conjuntamente la


comparación entre los cálculos de ETo con datos completos y los calculados con datos
faltantes de velocidad del viento y Radiación solar.

FIGURA 12.
Comparación de ETo – PM con datos completos y ETo PM calculado con valores estimados de
-1
Radiación solar y velocidad del viento = 2 m s (promedio de 2000 estaciones meteorológicas en todo
el mundo), para estaciones de Belén, Patacamaya, Viacha, Oruro y Potosí respectivamente.

77
INGENIERÍA Ingeniería Agronómica

TABLA 6.
Indicadores de desempeño para la comparación de resultados de ETo calculados con datos
completos y con datos faltantes de Radiación Solar y Velocidad del viento, para valores diarios.

Datos faltantes de Rs y U2
ESTACIÓN
RMSE* AAE** ARE***
METEOROLÓGICA
(mm d-1) (mm d-1) (%)
Belén 0.6245 0.5361 20.4271
Patacamaya 0.5855 0.4529 14.3306
Viacha 0.4018 0.3148 8.8812
Oruro 0.9442 0.8021 23.4756
Potosí 0.5127 0.4095 11.9976
* Error Cuadrado Medio, ** Promedio Absoluto del Error, *** Promedio Relativo del error

Cuando no se cuentan con datos de Rs y U2 los valores de RMSE, variaron de 0,4018 a


0,9442 mm d-1 para Viacha y Oruro, respectivamente, para el caso de Oruro el modelo no
muestra un comportamiento muy bueno debido a los cabios bruscos en la velocidad del
viento y a la variabilidad de los daros de radiación, no así en las demás estaciones donde
el comportamiento del modelo es bueno; el AAE de 0,3148 a 0,8021 mm d-1 para Viacha y
Oruro respectivamente, de igual forma en la estación de Oruro no se muestra una buen
resultado, pero para las demás estaciones la precisión de los resultados los esperados y
reflejados anteriormente; y un ARE de 8,8812 a 23,4756 % para Viacha y Oruro,

78
INGENIERÍA Ingeniería Agronómica

respectivamente, porcentajes que hacen que los datos no sean simétricos a los obtenidos
cuando se cuenta con todos los datos.

En tanto que el coeficiente de determinación R2 varió entre 0,64 a 0,77, para Patacamaya
y Viacha, respectivamente, resultados me ratifican la confiabilidad de los valores
obtenidos.

15.7. Estimación de la ETo cuando no se cuenta con datos de HR, Rs y U2.

El cálculo de la ETo con datos estimados de HR, Rs y U2 se realizo a partir de las


ecuaciones 23, 24 y un valor medio por defecto a nivel mundial de 2 m s-1, para la
velocidad del viento.

FIGURA 13.
Comparación de ETo – PM con datos completos y ETo PM calculado a partir de temperatura máxima
y mínima, para estaciones de Belén, Patacamaya, Viacha, Oruro y Potosí respectivamente.

79
INGENIERÍA Ingeniería Agronómica

80
INGENIERÍA Ingeniería Agronómica

TABLA 7.
Indicadores de desempeño para la comparación de resultados de ETo calculados con datos
completos y con datos faltantes de Humedad, Radiación Solar y Velocidad del viento, para valores
diarios.

Datos faltantes de HR, Rs y U2


ESTACIÓN
RMSE* AAE** ARE***
METEOROLÓGICA
(mm d-1) (mm d-1) (%)
Belén 0.6751 0.5869 22.2705
Patacamaya 0.7884 0.6486 20.7511
Viacha 0.4126 0.3104 9.0625
Oruro 0.8328 0.6939 20.5649
Potosí 0.6233 0.4952 14.5737
* Error Cuadrado Medio, ** Promedio Absoluto del Error, *** Promedio Relativo del error

Para el caso en el que no se cuenta con valores de HR, Rs y U 2 los valores de RMSE,
variaron de 0,4126 a 0,8328 mm d-1 para Viacha y Oruro, respectivamente, valores que
reflejan lo obtenido para cada variable independientemente y a la falta de las tres el
comportamiento es aceptable tomando en cuenta la variabilidad de los datos de (Rs), los
promedios de U2 distintos a la media recomendada y la diferencia de altitudes de las
estaciones meteorológicas; el AAE de 0,3104 a 0,6939 mm d-1 para Viacha y Oruro
respectivamente obteniéndose así la precisión de los resultados; y un ARE de 9,0625 a
22.2705 % para Viacha y Belén, respectivamente, porcentaje que hace que los resultados
no se comporten de manera simétrica en relación con los resultados obtenidos cuando se
tiene los datos completos.

En tanto que el coeficiente de determinación R2 varió entre 0,52 a 0,67, para Patacamaya
y Belén, respectivamente, y a la falta de estas tres variables aun se ratifica la confiabilidad
de los resultados obtenidos.

81
INGENIERÍA Ingeniería Agronómica

15.8. Cálculo de la ETo por la ecuación de Hargreaves para valores diarios

Se determinó la ETo con valores diarios por el método de Hargreaves (Ecuación 25)
utilizando únicamente datos de temperatura máxima y mínima, lo que fue comparado con
ETo-PM con datos completos. Los resultados se muestran en la Figura 14.

Se realizo una regresión lineal para realizar la comparación entre los cálculos de ETo con
datos completos y ETo calculado con Hargreaves para valores diarios.

FIGURA 14.

Comparación de ETo – PM con datos completos y ETo Hargreaves para valores diarios de
temperatura máxima y mínima, para estaciones de: Belén, Patacamaya, Viacha, Oruro y Potosí.

82
INGENIERÍA Ingeniería Agronómica

Para el caso del cálculo de ETo por el método de Hargreaves para valores diarios en
todas las localidades fue sobreestimada la ETo, a excepción de Viacha; principalmente en
Belén, debido a las bajas velocidades predominantes del viento, a la gran diferencia de
Tmax- Tmin y a la humedad relativa que fue superior a las demás conjuntamente
Patacamaya. En la localidad de Viacha hubo una ligera subestimación debido a la
predominancia de vientos moderados a fuertes (U2 = 3,62 m s-1). Se encontró un mejor

83
INGENIERÍA Ingeniería Agronómica

ajuste para la localidad de Belén, puesto que la diferencia entre Tmax- Tmin fue menor que
en otras localidades, debido a la influencia del lago Titicaca, que actúa como
termorregulador.

TABLA 8.
Indicadores de desempeño para la comparación de resultados de ETo calculados con datos
completos y ETo por el método de Hargreaves para valores diarios.

ETo PM – Hargreaves
ESTACIÓN (datos diarios)
METEOROLÓGICA RMSE* AAE** ARE***
(mm d-1) (mm d-1) (%)
Belén 0.4247 0.3321 12.4293
Patacamaya 0.6229 0.4851 14.8489
Viacha 0.4287 0.3576 10.0328
Oruro 0.6104 0.4717 13.1533
Potosí 0.6014 0.4773 13.1607
* Error Cuadrado Medio, ** Promedio Absoluto del Error, *** Promedio Relativo del error

Haciendo una comparación entre el método Hargreaves y el método FAO- PM los valores
de RMSE variaron de 0,4247 a 0,6229 mm d-1 para Belén y Patacamaya, respectivamente,
mostrando que método FAO-PM muestra un comportamiento bueno a comparación con el
método Hargreaves que llega a sobre estimar los resultados especialmente en Oruro,
Potosí y Patacamaya; el AAE de 0,3321 a 0,4851 mm d-1 para Belén y Patacamaya,
respectivamente, reflejando así la precisión de los cálculos y por ende los resultados ; y un
ARE de 10,0328 a 14,8489 % para Viacha y Patacamaya, respectivamente, porcentajes
que hacen que los resultados obtenidos no se distribuyan con simetría.

En tanto que el coeficiente de determinación R2 varió entre 0,52 a 0,70, para Patacamaya
y Belén, respectivamente, ratificando la confiabilidad de los resultados.

84
INGENIERÍA Ingeniería Agronómica

15.9. Hargreaves para periodos de cinco días

Se determinó la ETo para periodos de cinco días por el método de Hargreaves (Ecuación
25) utilizando únicamente datos de temperatura máxima y mínima, lo que fue comparado
con ETo-PM con datos completos. Los resultados se muestran en la Figura 15.

Al igual que en los anteriores casos, se realizo la regresión lineal para comparar los
resultados de los calculados por la ecuación de Penmman Montheit (ecuación 1) para
valores de cinco días.

FIGURA 15.

Comparación de ETo – PM con datos completos y ETo Hargreaves para periodos de cinco días, para
las estaciones de: Belén, Patacamaya, Viacha, Oruro y Potosí.

85
INGENIERÍA Ingeniería Agronómica

Para valores de ETo calculados en periodos de cinco días al igual que en el anterior caso,
a excepción de Viacha, en todas las localidades fue sobreestimada la ETo, principalmente
en Belén, debido a las bajas velocidades predominantes del viento, a la gran diferencia de
Tmax- Tmin y a la humedad relativa que fue superior a las demás conjuntamente
Patacamaya. En la localidad de Viacha hubo una ligera subestimación debido a la
predominancia de vientos moderados a fuertes (U2 = 3,62 m s-1).

86
INGENIERÍA Ingeniería Agronómica

También se encontró un mejor ajuste para la localidad de Belén, puesto que la diferencia
entre Tmax- Tmin fue menor que en otras localidades, debido a la influencia del lago
Titicaca, que actúa como termorregulador.

TABLA 9.
Indicadores de desempeño para la comparación de resultados de ETo calculados con datos
completos y ETo por el método de Hargreaves para periodos de cinco días.

ETo PM – Hargreaves
ESTACIÓN (periodos de 5 días)
METEOROLÓGICA RMSE* AAE** ARE***
(mm d-1) (mm d-1) (%)
Belén 0.3269 0.2599 9.4236
Patacamaya 0.5720 0.4504 14.0658
Viacha 0.4350 0.3636 9.8836
Oruro 0.7004 0.4921 16.7912
Potosí 0.4660 0.3689 10.0210
* Error Cuadrado Medio, ** Promedio Absoluto del Error, *** Promedio Relativo del error

Para periodos de 5 días los valores de RMSE, variaron de 0,3269 a 0,7004 mm d-1 para
Belén y Oruro, respectivamente, lo que indica que FAO presentan un comportamiento
aceptable y no así el método de Hargreaves, para periodos de 5 días; el AAE de 0,2599 a
0,4921 mm d-1 para Belén y Oruro, respectivamente, mostrando que el método Hargreaves
llega a subestimar y en algunos casos sobrestimar la precisión de los resultados; y un
ARE de 9,4236 a 16,7912% para Belén y Oruro, respectivamente, porcentaje de valores
que hacen que no exista una simetría en los resultados.

En tanto que el coeficiente de determinación (R2) varió entre 0,32 a 0,84, para Oruro y
Belén, respectivamente, ratificando la confiabilidad de los resultados, recomendando se
haga otras cálculos con otras estaciones en Oruro.

87
INGENIERÍA Ingeniería Agronómica

16. CONCLUSIONES

Los resultados de ETo estimados fueron obtenidos mediante los procedimientos


mencionados y fueron comparados con datos completos, que fueron tomados como
referencia. Las comparaciones fueron hechas mediante regresiones a partir del origen
entre los datos estimados y los de referencia, utilizando parámetros estadísticos como el
coeficiente de determinación (R2), el cuadrado medio del error (RMSE), el promedio
absoluto del error (AAE) y el promedio relativo del error (ARE).

Por los resultados obtenidos podemos concluir que:

Los valores de HR estimados a partir de la Ecuación 23 fueron similares a los observados


en las cinco localidades (Belén, Viacha, Patacamaya, Oruro y Potosí), lo que produjo
como consecuencia valores de ETo simulados similares a los valores de ETo – PM con
datos completos, lo que es un indicador de que la temperatura mínima es muy próxima a
la temperatura del punto de roció, especialmente en Belén.

Visto que hay pocos estudios de radiación en lugares de altitud y no hay dados de
observaciones en Bolivia que permitan una calibración de la ecuación de Angström, fue
necesario partir de la hipótesis que, siendo menos larga la espesura de la atmosfera en
altitud, se reduce la radiación difusa y aumenta la radiación directa. Así se pudo estimar la
radiación global en altitud.

Se puede advertir que las variables que más influyeron en el cálculo de ETo-PM con datos
limitados (ETo-PM Temp) fueron la radiación solar y el viento, debido a la sobreestimación
en el primer caso y la gran variabilidad en el segundo.

La ecuación de Hargreaves fue superior a la ETo-PM Temp en la localidad de Belén


debido a la menor variación entre la temperatura máxima y mínima por efecto de
atenuación del lago Titicaca, así como las bajas velocidades del viento. Este
comportamiento también fue observado por García et al. (2004).

88
INGENIERÍA Ingeniería Agronómica

En la localidad de Potosí se obtuvieron los valores más bajos de desempeño de las


diferentes ecuaciones, quizá porque estaba ubicada a mayor altura que las demás. En
zonas del Altiplano boliviano debido a las considerables alturas que se encuentran existe
gran variación de la temperatura diurna y nocturna, así como de la velocidad del viento y la
radiación solar presenta valores altos, sumado las bajas humedades del aire debido a las
condiciones de aridez y semiaridez, donde también estaría implicada la advección. Estos
aspectos dificultan la determinación de la ETo con datos limitados.

Los errores de observación para condiciones aerodinámicas variadas en las estaciones y


variaciones en las series de datos originan resultados y errores diferentes en la ETo
debido a la utilización de estimadores para la tensión de vapor actual a partir de la
temperatura mínima y del viento por un promedio de 2 m s -1. Esos errores influencian los
resultados de ETo cuando se estima con la FAO-PM solo con datos de temperatura.
Resulta así que no hay diferencias significativas cuando se compara con la ecuación
empírica de Hargreaves-Samani. Estudios más profundados son necesarios.

17. RECOMENDACIONES

Por los resultados obtenidos se recomienda el uso de la metodología de la FAO Penman –


Monteith para el cálculo de ETo en las localidades donde se realizo esta investigación
(zonas altas de Bolivia); ya que la misma toma en cuenta todos los parámetros climáticos
a comparación de la ecuación de Hargreaves.

Después de realizar la presente investigación y haber obtenido los resultados que se


mostraron anteriormente se recomienda para posteriores investigaciones, realizar
validaciones para otras estaciones meteorológicas con datos posteriores al de esta
investigación y para otras zonas del altiplano boliviano.

Se recomienda realizar como propuesta de investigación, la calibración de los valores de


as y bs para el cálculo de radiación solar (ecuación 14) puesto que estos son necesarios
para el cálculo de la ETo y para próximas validaciones.

89
INGENIERÍA Ingeniería Agronómica

Se recomienda el uso de la ecuación de Hargreaves para el cálculo de la ETo solamente


para periodos de cinco días, mensuales y anuales ya que el mismo sobre estima la ETo
para datos diarios como es el caso de esta investigación.

El método de Hargreaves se recomienda para localidades donde no es apropiado importar


datos de radiación de una estación cercana, ya sea porque no ocurren condiciones de
clima homogéneas, o porque no se cuenta con datos para la región.

Por último recomendar a nuestro Servicio Nacional de Hidrología y Meteorología


(SENAMHI) realizar los esfuerzos necesarios para brindarnos una información completa y
calidad para realizar investigaciones con una mayor exactitud que nos permitan mejorar de
una manera eficiente y eficaz el uso de los recursos hídricos y por consiguiente mejorar la
producción de los alimentos.

90
INGENIERÍA Ingeniería Agronómica

9. BIBLIOGRAFIA

 ALLEN, R.G., PEREIRA, L.S., RAES, D., SMITH, M., 1998. Crop
Evapotranspiration-Guidelines for Computing Crop Water Requirements. FAO
Irrigation and drainage paper 56. United. Nations Food and Agriculture
Organization, Rome.

 Allen, R.G., Smith, M., Perrier, A., y Pereira, L.S. 1994a. Una actualización para
la definición de evapotranspiración de referencia. ICID Boletin. 43(2). 1-34 p.

 Burman, R. y Pochop, L.O. 1994.La evaporación, datos de


evapotranspiración y Climática. Elsevier BV Ciencia, Amsterdam.

 Burman, R.D., Jensen, M.E., y Allen, R.G. 1987. Los factores termodinámicos en
la evapotranspiración. James, L.G. y M.J. inglés (editores), Proc. Riego y Drenaje.
La especificación. Conf., ASCE, Portland, Oregón, pág. 28-30.

 Bosque, H., Lemeur, R., Van Damme, P., 2000. Análisis ecofisiológico del cultivo
de la quinoa (Chenopodium quinoa Wild) en condiciones de estres de sequ´ıa y
salinidad. Tropicultura 18, 198–202.

 Bevan, K. 1979. Sistematización del análisis de Penman-Monteith para la


estimación de la evapotranspiración. J. Hydrol. 44:169-190 p.

 CHAVARRI VELARDE, Eduardo A. 2004, Universidad Nacional Agraria La Molina,


Facultad de Ingeniería Agrícola, Departamento de Recursos de Agua y Tierra.
Hidrología Aplicada, Pág. 1-10.

 David Sierra M., 2007, Apuntes de Botánica. Transpiración. Botanica Articulo 240
http://jardinactual.com/menu-revista-articulos/240 LA_TRANSPIRACIÓN_DE_LAS_PLANTAS,

 DOORENBOS J., PRUITT W. O., et.al. Las necesidades de agua de los cultivos,
Estudio FAO Riego y Drenaje 24 (Organización de las Naciones Unidas para la
Agricultura y la alimentación), Roma 1984.

91
INGENIERÍA Ingeniería Agronómica

 GARCIA, M., RAES, D., ALLEN, R., HERBAS, C., 2004. Dynamics of reference
evapotranspiration in the Bolivian highlands (Altiplano). Agric. Forest Meteorol. 125,
67–82.

 Hargreaves GL, Hargreaves GH, Riley JP. 1985. El periódico de Riego y Drenaje
Los beneficios agrícolas para la Senegal Río cubeta. 111(2): 113-124 p.

 HARGREAVES, G.H., SAMANI, Z.A., 1985. Reference crop evapotranspiration


from temperature. Appl. Eng. Agric. 1 (2), 96–99.

 HARGREAVES, G.H., ALLEN, R.G., 2003. History and evaluation of Hargreaves


evapotranspiration equation. J. Irrig. Drain. Eng. ASCE 129 (1), 53–63.

 Hanks, R.J. y Hill, R.W. 1980. Respuesta del modelo de riego para la cosecha
para suelos salinos y diferencias climáticas. Centro de Información de Irrigación
internacional, No. 6, Prensa de Pergamon, Elmsford, New York, Estados Unidos de
América, 63 pág.

 Hess, T., 1998. Tendencias en la evapotranspiración de referencia en el Norte


Oriental de Zonas Áridas de Nigeria, 1961-1991. Journal of Arid Environments. 38,
99-115.

 JENSEN, M.E., BURMAN, R.D., y ALLEN, R.G. (ed). 1990. Evapotranspiración y


requerimientos de Agua para riego. Manual de ASCE, Prácticas para diseñar
Informes Nº 70, Es. Soc. Ing. civil., Nueva York, Estados Unidos de América, 360 p.

 JABLOUN, M., SAHLI, A., 2008. Evaluation of FAO-56 methodology for estimating
reference evapotranspiration using limited climatic data Application to Tunisia. Agric.
Water Manage. 95, 707–715.

 MONTEITH, J.L., 1965. La evaporación y Medio Ambiente. Simposios 19 de la


Sociedad para la Biología Experimental, University Press, Cambridge, 19:205-234.

92
INGENIERÍA Ingeniería Agronómica

 PEREIRA, L.S., PERRIER, A., ALLEN, R.G. y ALVES, I. 1996. Evapotranspiración:


Revisión de conceptos y tendencias de flujo. J. Riego y Drenaje. Ing. Agro. ASCE
25. (prensa).

 Penman, H. L. 1948. La evaporación natural del agua abierta, tierra desnuda y


césped. Proc. Roy. Soc. Londres, A193, 120-146 p.

 Penman, H.L. 1963. La vegetación y la hidrología. Tech. Com. N º 53, de la


Commonwealth Mesa de Suelos, Harpenden, Reino Unido. 125 p.

 PERRIER, A. 1978. Importance des définitions de l’évapotranspiration dans le


domaine pratique de la mesure, de l’estimation of de la notion de coefficients
culturaux. XV’ Journal of Hydraulics, Société Hydrotechnique de France, Question
IV, Rapport 1:1-7 (en Francés).

 ALLEN Richard G., PEREIRA Luis S., RAES Dirk, SMITH Martin, 2006.
Evapotranspiración del cultivo. Guías para la determinación de los requerimientos
de agua de los cultivos. Servicio de Recursos, Fomento y Aprovechamiento de
Aguas. FAO (organización de las naciones unidas para la agricultura y la
alimentación) roma, estudio FAO riego y drenaje Nº 56. pp 323.

 RIJTEMA, P.E., (1965). Analysis de la evapotranspiración. Agric. Res. Reporte. No.


69, Centre de Agric. Publ. y Doc., Wageningen, Países Bajos.

 SMITH, M., ALLEN, R.G., MONTEITH, J.L., PERRIER, A., PEREIRA, L., y
SEGEREN, A. 1992. Report of the expert consultation on procedures for revision of
FAO guidelines for prediction of crop water requirements. Organización de las
Naciones Unidas para la Agricultura y la Alimentación, Roma, 54 p.

93
INGENIERÍA Ingeniería Agronómica

ANEXO 1

UNIDADES Y SÍMBOLOS
Prefijos
Las unidades pueden ser usadas como tales o como múltiplos

Tera (T) que significa 1012 Centi (c) 10-2


Giga (G) 109 Mili (m) 10-3
Mega (M) 106 Micro (μ) 10-6
Kilo (k) 103 Nano (n) 10-9
Hecto (h) 102 Pico (p) 10-12
Deca (da) 101 Femto (f) 10-15
Deci (d) 10-1 Atto (a) 10-18

Temperatura
Unidad estándar: grados Celsius (°C)
_________________________________________________________
Grados Fahrenheit (°F) (°C) = (°F-32) 5/9
Kelvin (K) 1 K = (°C) + 273,16

PRESION
(Presión del aire, presión de vapor)
Unidad estándar: kilopascal (kPa)
_______________________________________________________________
milibar (mbar) 1 mbar = 0,1 kPa
bar 1 bar = 100 kPa
centímetro de agua (cm) 1 cm de agua = 0,09807 kPa
milímetro de mercurio (mm Hg) 1 mm Hg = 0,1333 kPa
atmósferas (atm) 1 atm = 101,325 kPa
libras por pulgada cuadrada (psi) 1 psi = 6,896 kPa

94
INGENIERÍA Ingeniería Agronómica

VELOCIDAD DEL VIENTO


Unidad estándar: metro por segundo (m s-1)
________________________________________________________________
kilómetro por día (km día-1) 1 km día-1 = 0,01157 m s-1
milla náutica /hora (nudo) 1 nudo = 0,5144 m s-1
pies por segundo (ft s-1) 1 ft/s = 0,3048 m s-1

RADIACIÓN
Unidad estándar: mega joule por metro cuadrado por día (MJ m-2 día-1)
o como evaporación equivalente en mm por día (mm día-1)
____________________________________________________________________
evaporación equivalente (mm día-1) 1 mm día-1 = 2,45 MJ m-2 día-1
joule por cm2 por día (J cm-2 día-1) 1 J cm-2 día-1 = 0,01 MJ m-2 día-1
caloría por cm2 por día (cal cm-2 día-1) 1 cal = 4,1868 J = 4,1868 10-6 MJ
1 cal cm-2 día-1 = 4,1868 10-2 MJ m-2 día-1
watt por m2 (W m-2) 1 W = 1 J s-1
1 W m-2 = 0,0864 MJ m-2 día-1

EVAPOTRANSPIRACIÓN
Unidad estándar: milímetro por día (mm día-1)
__________________________________________________________________
m3 por hectárea por día (m3 ha-1 día-1) 1 m3 ha-1 día-1 = 0,1 mm día-1
litro por segundo por hectárea (l s-1 ha-1) 1 l s-1 ha-1 = 8,640 mm día-1
radiación equivalente en mega joules por
metro cuadrado por día (MJ m-2 día-1) 1 MJ m-2 día-1 = 0,408 mm día-1

95
INGENIERÍA Ingeniería Agronómica

ANEXO 2

Presión atmosférica (P) para diferentes altitudes (z)

96
INGENIERÍA Ingeniería Agronómica

ANEXO 3

Constante psicrométrica (γ) para diferentes altitudes (z)

97
INGENIERÍA Ingeniería Agronómica

ANEXO 4

Presión de saturación de vapor (eº(T)) para diferentes temperaturas (T)

98
INGENIERÍA Ingeniería Agronómica

ANEXO 5

Pendiente de la curva de presión de vapor (Δ) para diferentes temperaturas (T)

99
INGENIERÍA Ingeniería Agronómica

ANEXO 6
Número de día en el año (J)

Continua…….

J puede ser determinado para cada día (D) del mes (M) usando
J = NUMERO ENTERO (275 M/9 – 30 + D) – 2
SI (M < 3) ENTONCES J = J + 2
ó también, SI (año bisiesto y (M > 2)) ENTONCES J = J + 1

Para cálculos de 10 días, use J para día D = 5, 15 y 25


Para cálculos mensuales, J en el día 15 del mes viene dado aproximadamente por

J = NUMERO ENTERO (30,4 M – 15)

100
INGENIERÍA Ingeniería Agronómica

101
ANEXO 7
Radiación extraterrestre diaria (Ra) para diferentes latitudes para el día 15vo del mes1
INGENIERÍA
Ingeniería Agronómica

102
ANEXO 8
Insolación máxima diaria (N) para diferentes latitudes para el día 15vo del mes1
INGENIERÍA
Ingeniería Agronómica

103

También podría gustarte

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy