Universidad Mayor de San Andrés Facultad de Agronomía Carrera de Ingeniería Agronómica
Universidad Mayor de San Andrés Facultad de Agronomía Carrera de Ingeniería Agronómica
Universidad Mayor de San Andrés Facultad de Agronomía Carrera de Ingeniería Agronómica
TESIS DE GRADO
LA PAZ – BOLIVIA
2011
1
INGENIERÍA Ingeniería Agronómica
Asesores:
Tribunal Examinador
Aprobada
2011
2
INGENIERÍA Ingeniería Agronómica
PROLOGO
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.
3
INGENIERÍA Ingeniería Agronómica
4
INGENIERÍA Ingeniería Agronómica
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
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
Í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
Í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
13
INGENIERÍA Ingeniería Agronómica
14
INGENIERÍA Ingeniería Agronómica
1. INTRODUCCION
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.
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
Estos coeficientes dependen de varios factores incluyendo tipo de cultivo, la fase del
crecimiento de los cultivos, la altura del dosel y su densidad5.
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 .
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.
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:
19
INGENIERÍA Ingeniería Agronómica
1.2 Justificación
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
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
22
INGENIERÍA Ingeniería Agronómica
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.
12.1.1. Evaporación
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).
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
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
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
FIGURA 2.
Estoma visto a través de un microscopio.
26
INGENIERÍA Ingeniería Agronómica
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:
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 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.
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.
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
FIGURA 3.
Repartición de la evapotranspiración en evaporación y transpiración durante el periodo de
crecimiento de un cultivo anual.
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.
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
30
PERRIER, A. 1978.
31
ALLEN Richard G. et al, 2006. Pág. 5
30
INGENIERÍA Ingeniería Agronómica
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
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.
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.
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
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.
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
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.
CUADRO 1
Características y localización de las estaciones meteorológicas de las que se utilizaron los
parámetros climáticos.
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"
36
INGENIERÍA Ingeniería Agronómica
FIGURA 5.
Ubicación Geográfica de las estaciones meteorológicas utilizadas para el estudio.
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
38
INGENIERÍA Ingeniería Agronómica
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.2. Equipos
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
FIGURA 6.
Características del cultivo de referencia
40
INGENIERÍA Ingeniería Agronómica
(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.
41
INGENIERÍA Ingeniería Agronómica
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í.
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.
(2)
Donde:
Tmedia : Temperatura media (ºC)
Tmáx : Temperatura Máx. (ºC)
Tmín : Temperatura Mínima (ºC)
43
INGENIERÍA Ingeniería Agronómica
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)
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 [..]
(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.
45
INGENIERÍA Ingeniería Agronómica
(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 [..]
(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.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).
47
INGENIERÍA Ingeniería Agronómica
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)
(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.
(12)
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).
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).
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.
(13)
Donde:
ωs : Es el ángulo de radiación a la hora de la puesta del sol calculado con
Ecuaciones 12.
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
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).
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.
(16)
Donde
z : Elevación de la estación sobre el nivel del mar [m].
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.
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
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
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.
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].
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
55
INGENIERÍA Ingeniería Agronómica
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.
(22)
Donde
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.
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)
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.
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.
Este principio ha sido utilizado por Hargreaves y Samani, para desarrollar estimaciones de
ETo usando solamente datos de temperatura del aire.
(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].
(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.
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.
(26)
Donde:
60
INGENIERÍA Ingeniería Agronómica
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:
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
(28)
Donde:
n : Es el número de datos de las series Oi (medida) y Pi (simulada).
(29)
Donde:
Oi y P i :(i = 1, 2, ..., n) son pares de observación
(30)
Donde:
Oi y P i :(i = 1, 2, ..., n) son pares de observación
62
INGENIERÍA Ingeniería Agronómica
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%.
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
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ó.
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.
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
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.
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.
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.
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
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
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.
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.
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
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
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
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.
88
INGENIERÍA Ingeniería Agronómica
17. RECOMENDACIONES
89
INGENIERÍA Ingeniería Agronómica
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.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.
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.
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.
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.
92
INGENIERÍA Ingeniería Agronómica
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.
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
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
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
96
INGENIERÍA Ingeniería Agronómica
ANEXO 3
97
INGENIERÍA Ingeniería Agronómica
ANEXO 4
98
INGENIERÍA Ingeniería Agronómica
ANEXO 5
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
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