Estimacin de la Relacin Adimensional de Budyko en Colombia

Paola Andrea Arias Gmez

[email protected]
Germn Poveda Jaramillo
[email protected]
Posgrado en Aprovechamiento de Recursos Hidrulicos
Escuela de Geociencias y Medio Ambiente
Universidad Nacional de Colombia
Sede Medelln

Preparado para presentacin en el

XVI Seminario Nacional de Hidrulica e Hidrologa
Sociedad Colombiana de Ingenieros
Sociedad de Ingenieros del Quindo
Universidad del Quindo
Corporacin Autnoma Regional del Quindo
Armenia 29, 30 y 31 de octubre de 2004

RESUMEN
Los balances de agua y energa sobre cuencas hidrogrficas condicionan la formacin del terreno, la
distribucin y productividad de la vegetacin, etc. Se define el nmero adimensional de Budyko, B, como la
relacin entre la precipitacin media anual (P) y la evapotranspiracin potencial media anual (PET),
B=P/PET, en cuencas hidrogrficas. Este nmero B resume el balance entre la disponibilidad hdrica (P) y la
disponibilidad energtica (PET) y ha sido empleado para caracterizar la disponibilidad hdrica de las plantas
y la aridez y la productividad primaria neta de los ecosistemas. Estudios clsicos han encontrado que para
B=1, que denominamos nmero de Budyko crtico, Bc, se presentan condiciones climatolgicas,
geomorfolgicas y de biodiversidad nicas que hacen que esta relacin adimensional sea de especial inters
en hidroclimatologa y ecologa. En este trabajo se presentan los mapas de la estimacin del nmero de
Budyko para Colombia a escala de 5 minutos, y para Antioquia a escala de 30 segundos de arco,
implementando diversos mtodos para la estimacin de la evaporacin potencial anual. Se concluye que en
Colombia predomina la vegetacin limitada por energa.

INTRODUCCIN
En un artculo clsico, Melton (1958) presenta una relacin emprica entre la densidad de drenaje de las

cuencas hidrogrficas, Dd, un ndice de la geometra de la red de drenaje que refleja la erosin del terreno
por el agua, y el ndice P-E de Thornthwaite para 80 cuencas en el suroeste rido de Estados Unidos. Ms
adelante Madduma-Bandara (1974) extendi este hallazgo al clima tropical de Sri Lanka al demostrar que Dd
exhibe una relacin directa con el ndice P-E para 24 cuencas (Abrahams 1984). El ndice P-E se define
como diez veces la suma de los 12 valores mensuales de la relacin entre la precipitacin ( Pi ) y la
evaporacin potencial ( PETi ), i=1,2,12, (Thornthwaite 1933). El entendimiento de esta relacin no lineal
sigue siendo un tema de investigacin sobre las relaciones entre el terreno, la vegetacin, el agua y la
energa (Milne et al. 2002). Para superar las limitaciones de la variabilidad estacional, se propone un ndice
climtico adimensional como la relacin entre la precipitacin promedio anual precipitacin ( P ), y la
evaporacin promedio anual ( PET ), ambas con dimensiones [L/T], (Milne et al., 2002),

B = P / PET = Pl / R
en donde

PET

hmeda, PET

se asume como el equivalente en agua de la radiacin neta,

= R /l ,

(1)

para una superficie

y donde l es el calor latente de vaporizacin. Al parmetro adimensional B lo

denominaremos el Nmero de Budyko en honor de M. Budyko quien hizo contribuciones fundamentales al

entendimiento de las relaciones entre el clima y la vida (Budyko 1974). El ndice P-E y B estn relacionados
por un factor del orden de 100. Por tanto, cuando se quieren inferir relaciones entre
relacin original entre

Dd

Dd

y B, el eje x de la

y el ndice P-E se debe reescalar adecuadamente por tal factor. De cualquier

forma, la relacin no cambiar cualitativamente por esta transformacin.

El Nmero de Budyko constituye una componente fundamental de los balances de agua y energa en
cuencas hidrogrficas, y tiene aplicaciones fundamentales como el ndice de Aridez (UNESCO, 1979), la
prediccin de productividad primaria neta de ecosistemas, en la descomposicin y en la estimacin de
contenido de Carbono en el suelo (Ojima et al. 1996), en el efecto del cambio climtico en los procesos de
la hidrologa superficial (Koster et al. 2001), y en la definicin de vegetacin limitada por agua, por energa o
por ambas (Nemani et al., 2003), y en la regulacin natural de la biodiversidad (Olff et al. 2002). Hay una
transicin de fase entre vegetacin limitada por agua a vegetacin limitada por energa en el Nmero Crtico
de Budyko, Bc=1.0. En este trabajo se estima el nmero adimensional de Budyko para Colombia, y se
discuten distintas aplicaciones alrededor de su estimacin.

INFORMACIN EMPLEADA
Para la realizacin de este trabajo se emplearon mapas de diferentes variables hidroclimatolgicas

obtenidos durante los proyectos Balances Hidrolgicos de Colombia (UNALMED-UPME (1999)), Atlas
Hidrolgico de Colombia (UNALMED-UPME (2000)) y Balance Hidrolgico y Atlas Digital de Antioquia
(UNALMED-CTA (2001)), realizados por el grupo de investigacin del Posgrado en Aprovechamiento de
Recursos Hidrulicos de la Universidad Nacional de Colombia Sede Medelln.
Los mapas obtenidos en el proyecto Balances Hidrolgicos de Colombia (UNALMED-UPME (1999))
corresponden a todo el territorio colombiano y tienen una resolucin espacial de 5 minutos de arco. Por su
parte, los mapas obtenidos en el proyecto Balance Hidrolgico y Atlas Digital de Antioquia (UNALMED-CTA
(2001)) corresponden a todo el departamento de Antioquia y tienen una resolucin espacial de 30 segundos
de arco.

MARCO TERICO
La vegetacin es una componente determinante en el comportamiento climtico y energtico de los

ecosistemas terrestres. Es uno de los factores ms representativos del ciclo global del carbono, ya que
define las tasas de almacenamiento y de liberacin del CO2. Su presencia modifica por completo las
caractersticas de los balances hdrico y energtico, por lo que es necesario considerar su interaccin con
los factores climticos, de manera que el entendimiento de las dinmicas y de los procesos fsicos
relacionados con la productividad primaria neta sea completo. Esto constituye un gran aporte para un mejor
discernimiento en un tema de tanta importancia en la vida de los ecosistemas terrestres como lo es el
calentamiento global.

3.1

Evapotranspiracin

La evapotranspiracin es una componente fundamental en el balance hidrolgico, aunque es, a su vez, la

variable con mayores incertidumbres en su medicin y estimacin. Est compuesta por la evaporacin desde
superficie y por la transpiracin de las plantas. La dinmica de la evapotranspiracin est condicionada por
tres factores fundamentales: disponibilidad de agua, disponibilidad de energa para cambio de estado del
agua y condiciones aerodinmicas propicias para asimilacin y transporte de vapor de agua. En sntesis, la
evapotranspiracin resume la interaccin entre la disponibilidad hdrica y la disponibilidad energtica,
interaccin que se ve afectada por las condiciones ambientales existentes, siendo la vegetacin uno de los
factores determinantes en su comportamiento.
A pesar de la existencia de una gran variedad de propuestas para la estimacin de la evapotranspiracin
(Brutsaert (1982)), no es claro cul es la ms adecuada, ms an si se considera que el conocimiento que
se tiene acerca de esta variable es reducido, adems de las claras dificultades que existen para medirla y
estimarla. De esta variedad de metodologas, algunas de las cuales han sido ampliamente empleadas para

la estimacin de la evapotranspiracin (tanto potencial como real), muchas son basadas en regresiones
obtenidas de registros empricos correspondientes a regiones con caractersticas climticas, energticas y
de vegetacin particulares, lo que limita la validez de los resultados obtenidos al aplicarlas en regiones con
caractersticas diferentes. Otras de estas metodologas parten de suposiciones y de desarrollos tericos que
no siempre consideran todas las variables implicadas o que no dan cuenta adecuadamente de los procesos
fsicos inherentes a este fenmeno. De hecho, la gran mayora de las metodologas clsicas para estimar la
evapotranspiracin no incluyen el factor vegetacin, totalmente determinante en el comportamiento de esta
variable.
Aunque el objetivo central de este trabajo no es ahondar en el tema de la evapotranspiracin, realizando
la revisin bibliogrfica para desarrollarlo se encontr una nueva metodologa para la estimacin de la
evapotranspiracin real que considera la vegetacin como una de las variables requeridas para su
aplicacin. Dadas las consideraciones tericas de esta nueva formulacin, puede convertirse en una
herramienta de fcil aplicacin para la estimacin de una variable tan compleja como la evapotranspiracin.
Esta metodologa fue propuesta por Zhang et al (Zhang et al (2001)). A continuacin se describe brevemente
su fundamentacin. En la Seccin 4.1 se evalan los resultados de su aplicacin en las regiones estudiadas
en este trabajo (Colombia y Antioquia).

3.1.1 Mtodo de Zhang et al (2001) para la estimacin de la evapotranspiracin real

Una manera de aproximarse a la evapotranspiracin es asumir que la evapotranspiracin desde
superficie est controlada por la disponibilidad de agua y por la demanda de la atmsfera. La disponibilidad
de agua puede ser aproximadamente igual a la precipitacin; la demanda atmosfrica representa la mxima
evapotranspiracin posible (evapotranspiracin potencial).
Bajo condiciones muy secas, la evapotranspiracin potencial excede la precipitacin, y la
evapotranspiracin real es igual a la precipitacin. Bajo condiciones muy hmedas, la disponibilidad de agua
excede la evapotranspiracin potencial y la evapotranspiracin real se aproximar asintticamente a la
evapotranspiracin potencial. Basndose en estas consideraciones, Budyko (1958) postul las siguientes
relaciones bajo condiciones secas

ETP 1

Rn / P

(1)

donde P es la precipitacin, ET es la evapotranspiracin real y Rn es la radiacin neta. Para condiciones

muy hmedas se tiene

ET R n

Rn / P 0

(2)

Zhang et al (2001) encontraron que la siguiente expresin satisface las condiciones (1) y (2)

PET
ET
P
=
1
P
PET PET
+
1+ w

P
P
1+ w

(3)

donde w es el coeficiente de disponibilidad de agua en las plantas. En las expresiones (1) y (2)
planteadas por Budyko, la radiacin neta Rn es un subrogado de la evapotranspiracin potencial. Sin
embargo, Zhang et al introdujeron el trmino PET (evapotranspiracin potencial) en lugar de Rn. Por esto
aparece la variable PET en la expresin (3).
El efecto de w en la evapotranspiracin es mnimo bajo condiciones muy secas o muy hmedas. La
mxima diferencia en la relacin entre evapotranspiracin real y precipitacin entre rboles y plantas
herbceas ocurre cuando la precipitacin anual es igual a la demanda atmosfrica (PET/P=0). Bajo esta
condicin, la habilidad de los rboles de extraer agua de zonas ms profundas les permite que usen el agua
que ha sido almacenada en los perodos de inactividad, mientras que las plantas herbceas con races ms
profundas pueden permitir que el agua se escape de su zona de races.
De manera que la expresin (3) sea vlida para regiones con proporciones variables de bosques y
cultivos, Zhang et al propusieron un modelo generalizado. Siguiendo el trabajo de Eagleson (1982), estos
autores asumieron que la evapotranspiracin anual de la vegetacin herbcea de la regin en estudio
(incluyendo la evaporacin desde el suelo) y la de los bosques est ponderada linealmente de acuerdo con
sus reas. De esta manera

ET = fET f + (1 f ) ETh

(4)

donde ET es la evapotranspiracin real total anual en mm, f es la fraccin de cobertura vegetal, ETf es la
evapotranspiracin real anual en los bosques (mm) y ETh es la evapotranspiracin real anual de las plantas
herbceas (en mm).
A partir de los anlisis realizados por Zhang et al, stos consideraron que la evapotranspiracin potencial
PET puede ser asumida como constante para las zonas boscosas y para las zonas de plantas herbceas,
donde la constante depende del tipo de vegetacin. As mismo, asumieron que la disponibilidad de agua en
la planta es una constante que depende del tipo de vegetacin. De esta manera, los valores encontrados por
Zhang et al fueron PET=1410 mm y w=2.0 para bosques y PET=1100 mm y w=0.5 para zonas herbceas.
De esta manera, la expresin (4) puede ser expresada como

1410
1100

1+ 2*
1 + 0.5 *

P
P

P
+ (1 f )
ET = f
1410
1100
P
P

+
1 + 0.5 *
+
1+ 2*

1410
1100
P
P

(5)

Como puede observarse de la expresin (5), el planteamiento generalizado propuesto por Zhang et al
para la estimacin de la evapotranspiracin real incluye tres factores fundamentales en el comportamiento
de esta variable: la disponibilidad de agua, representada por P; la disponibilidad de energa, representada
por PET; y las condiciones vegetales representadas por w y f. Esta expresin es de fcil aplicacin, aunque
parte del supuesto de que para regiones con igual tipo de cobertura, la evapotranspiracin potencial es
constante, as como la disponibilidad de agua en las plantas. Es necesario ser cuidadosos con estas
suposiciones, ya que en algunos casos pueden no ser vlidas, por tanto, la aplicacin de esta metodologa
requiere su previa verificacin.

3.1.2 La radiacin neta como un estimativo de la evapotranspiracin potencial

Como se puede ver de las expresiones (1) y (2), Budyko (1958) emplea la radiacin neta Rn como una
aproximacin a la evapotranspiracin potencial PET. Esta aproximacin es tal que PETRn/Lv, donde Lv es
el calor latente de vaporizacin del agua.
Budyko (1974) argument que para el caso de promedios anuales, el lmite superior de la
evapotranspiracin es igual a la radiacin neta ya que la sumas anuales del flujo de calor latente no
proporciona una entrada significativa de energa a la superficie terrestre; por esta razn, el flujo de calor
latente puede estar determinado por la radiacin neta.

3.2

La relacin adimensional de Budyko

Una manera muy aceptada de resumir la hidroclimatologa de largo plazo en una regin es graficando la
relacin entre la evapotranspiracin real media y la precipitacin media (ET/P) vs la relacin entre la
evapotranspiracin potencial media y la precipitacin media (PET/P) (Sankarasubramanian y Vogel (2002)).
Los ndices ET/P y PET/P son denominados relaciones de evapotranspiracin y de aridez, respectivamente.
Budyko (1974), lanz la hiptesis de una relacin funcional generalizada entre la relacin de
evapotranspiracin y el ndice de aridez, tal que

ET
PET
=

P
P
Ms exactamente,

(6)


PET
P
ET = P * 1 exp
* PET * Tanh

P
PET

0.5

(7)

La expresin (7) es una herramienta acertada para la estimacin de la evapotranspiracin real a partir de
la evapotranspiracin potencial y que ha sido ampliamente empleada en la literatura.

RESULTADOS Y CONCLUSIONES

4.1

Aplicacin del mtodo de Zhang et al (2001) en Colombia

Para verificar qu tan vlido es el modelo propuesto por Zhang et al para la estimacin de la
evapotranspiracin real, ste fue aplicado en las regiones de estudio: Colombia y Antioquia.
Para ello se emplearon los mapas de precipitacin obtenidos en los trabajos Balances Hidrolgicos de
Colombia (UNALMED-UPME (1999)) para Colombia (con una resolucin espacial de 5 minutos de arco) y
Balance Hidrolgico y Atlas Digital de Antioquia (UNALMED-CTA (2001)) para Antioquia (con una resolucin
espacial de 30 segundos de arco). Estos mapas fueron obtenidos por interpolacin mediante Kriging con
deriva externa, donde la deriva empleada fue el mapa de precipitacin del Estudio Nacional del Agua (ENA),
para el caso de Colombia, y por interpolacin con triangulacin con deriva externa por subregiones, donde la
deriva fue el mapa de precipitacin del ENA reescalado, para Antioquia.
Con respecto a la fraccin de cobertura vegetal, sta fue estimada empleando el mapa de cobertura
vegetal obtenido en el trabajo Atlas Hidrolgico de Colombia (UNALMED-UPME (2000)). As, se consider
que f corresponde al porcentaje de pixeles en los cuales se tiene cobertura vegetal boscosa con respecto al
nmero total de pixeles del mapa de cobertura en consideracin. El valor obtenido es de 0.49 para Colombia
y de 0.29 para Antioquia.
De esta manera se estim la evapotranspiracin real en Colombia mediante la expresin (5). Como ya se
coment, la expresin (5) tiene implcita la consideracin de evapotranspiracin potencial y disponibilidad de
agua en las plantas constantes para cada tipo de cobertura, cuyos valores fueron obtenidos con informacin
de 68 cuencas de Estados Unidos de Norteamrica. Por esta razn, se consider evapotranspiracin
potencial constante para cada tipo de cobertura, pero los valores se obtuvieron a partir de diferentes mapas
de evapotranspiracin potencial elaborados en los trabajos Balances Hidrolgicos de Colombia para
Colombia (con una resolucin espacial de 5 minutos de arco) y Balance Hidrolgico y Atlas Digital de
Antioquia para Antioquia (con una resolucin espacial de 30 segundos de arco). Los mapas de
evapotranspiracin potencial considerados fueron los obtenidos por las metodologas de Cenicaf, Morton,
Penman, Thornwaite, Turc y Turc Modificado.

Con respecto a la disponibilidad de agua en las plantas, ya que en nuestro medio no se cuenta con
informacin significativa de esta variable, se consideraron los valores propuestos por estos autores: 2.0 para
bosques y 0.5 para plantas herbceas. De esta manera, se obtuvieron seis mapas de evapotranspiracin
real mediante el modelo de Zhang et al, sustituyendo los valores de evapotranspiracin potencial (1410 mm
y 1100 mm) por los encontrados a partir de cada uno de los mapas anteriores.
Para validar los resultados obtenidos, se estimaron los caudales medios de largo plazo en diferentes
cuencas instrumentadas de la regin en estudio (108 cuencas en Colombia y 56 cuencas en Antioquia) va el
balance hdrico de largo plazo (P-ET=R, donde R es la escorrenta superficial) y se compararon con los
caudales reales en las cuencas. Para la extraccin de las cuencas se emple un MDT de 5 minutos de arco
para Colombia y otro de 30 segundos de arco para Antioquia, ambos suministrados por el USGS.
Los resultados obtenidos para Colombia muestran que la metodologa de Zhang et al presenta errores
bastante altos para Colombia (del orden del 80%). Estos errores pueden deberse, en gran medida, a la
calidad del MDT empleado para la extraccin de las cuencas, ya que en muchos casos, stas no coinciden
con las cuencas reales, por lo que los resultados se ven distorsionados. Vale la pena aclarar que en un
principio se contaba con 303 cuencas para Colombia, pero se debi reducir a 108 cuencas debido a la mala
calidad del MDT para dichas cuencas.
Con respecto a la aplicacin del modelo de Zhang et al en Antioquia, los errores obtenidos son mucho
menores que los obtenidos para Colombia (del orden del 20%). Estos errores en el balance son bastante
satisfactorios. La disminucin sustancial en el error se debe a la mejor calidad del MDT empleado y a la
mejor resolucin de los mapas de las variables empleadas.
En general, se obtuvo que los menores errores en la estimacin de los caudales medios anuales, tanto
para Colombia como para Antioquia, se presentan al emplear la evapotranspiracin potencial de Cenicaf y
de Morton, que son dos de las metodologas que han mostrado mejores estimativos de la evapotranspiracin
real en Colombia (UNALMED-UPME (1999), UNALMED-UPME (2000)), UNALMED-CTA (2001)). Por su
parte, los mayores errores se presentan al emplear la evapotranspiracin potencial de Penman y de Turc
Modificado. Tanto para Colombia como para Antioquia, se encuentra que al emplear los valores de PET
dados por Zhang et al se obtienen errores intermedios a los obtenidos por las seis metodologas para
estimar PET consideradas.
De esta manera, se encuentra que la metodologa de Zhang et al es una herramienta que entrega
resultados aceptables. Debe avanzarse en el estudio del parmetro w, que resume la disponibilidad de agua
en las plantas, de manera que se obtengan mejores resultados.

4.2

Verificacin de la aproximacin PETRn/Lv en Colombia

Para verificar que tan vlida es la aproximacin PETRn/Lv para las condiciones hdricas, energticas y
de vegetacin en Colombia, se estim el error de Rn/Lv con respecto a PET, tanto para todas las regiones en
estudio, como para las cuencas seleccionadas, para los seis mtodos de evapotranspiracin potencial
considerados. Los mapas de radiacin neta empleados corresponden a los obtenidos en los trabajos
Balances Hidrolgicos de Colombia (UNALMED-UPME (1999)) y Balance Hidrolgico y Atlas Digital de
Antioquia (UNALMED-CTA (2001)). Estos mapas se estimaron mediante la metodologa propuesta por
Morton (1982) y estn a una resolucin de 5 minutos y 30 segundos de arco, respectivamente.
Los resultados obtenidos muestran que en promedio el error en esta aproximacin es del 40% para
Colombia y del 20% para Antioquia. Nuevamente, estas diferencias en los errores se deben a la mejor
resolucin y calidad de los mapas y del MDT de Antioquia. La aproximacin considerada presenta menores
diferencias con las metodologas de Cenicaf y de Morton, para ambas resoluciones.
Algo que debe considerarse en este anlisis, es que la radiacin neta empleada no corresponde a
mediciones sino que fue estimada por una metodologa iterativa. Por lo tanto, existen incertidumbres no slo
en los mapas de evapotranspiracin potencial sino tambin en el mapa de radiacin neta.
En general, es posible afirmar que la aproximacin propuesta por Budyko (1958), PETRn/Lv, se cumple
aceptablemente en la regin de estudio, por lo que se puede convertir en un estimativo de la
evapotranspiracin potencial en regiones en las que no se tenga informacin de sta.

4.3

Estimacin de la relacin adimensional de Budyko en Colombia

Considerando la relacin adimensional de Budyko dada por la expresin (7), se analiz la

hidroclimatologa de largo plazo para las cuencas en consideracin (108 para Colombia y 56 para Antioquia).
Para ello se emplearon los mapas de precipitacin media anual ya mencionados y los mapas de
evapotranspiracin real y potencial obtenidos por las metodologas de Cenicaf, Morton, Penman,
Thornwaite, Turc y Turc Modificado (UNALMED-UPME (1999) y UNALMED-CTA (2001)). En la Figura 1 se
observan los resultados obtenidos para las cuencas en Colombia con la evapotranspiracin potencial de
Cenicaf y la evapotranspiracin real estimada por cuatro de las seis metodologas consideradas.
Las asntotas A y B en la Figura 1 corresponden a los lmites superiores en la relacin de evaporacin o a
los lmites inferiores en la escorrenta. Para regiones con abastecimiento de humedad ilimitado, la
evapotranspiracin real se acerca a la evapotranspiracin potencial, lo que est ilustrado por la lnea A.
Similarmente, cuando la evapotranspiracin potencial excede la precipitacin, la evapotranspiracin real se
acerca a la precipitacin, que se muestra por la lnea B.

Figura 1.

Relacin ET/P vs PET/P para las 108 cuencas en consideracin en Colombia. PET es estimada mediante
la metodologa de Cenicaf. ET es estimada por los mtodos de Cenicaf, Morton, Penman y Thornwaite.
La precipitacin es interpolada por Kriging con deriva externa (UNALMED-UPME (1999)).

Los resultados obtenidos tanto para Colombia como para Antioquia son similares. Se observan valores
del ndice de aridez PET/P inferiores a 2.0, lo que no permite observar el comportamiento asinttico con la
lnea B, esto es, regiones en las que la evapotranspiracin potencial excede la precipitacin. As, los
resultados muestran que la mayora de nuestras cuencas estn en ambientes limitados por energa, esto es,
existe una gran cantidad de agua para disponible, pero sta no es evaporada debido a limitaciones
energticas.
En algunas de las metodologas se observan grandes dispersiones en los resultados de las diferentes
cuencas, sin embargo, la metodologa de Cenicaf, tanto para estimar ET como PET, es la que muestra un
comportamiento ms uniforme para las distintas cuencas.
Al inverso del ndice de aridez PET/P lo denominaremos B=P/PET (Nmero climtico de Budyko). Este
nmero B resume el balance entre la disponibilidad hdrica (P) y la disponibilidad energtica (PET) y ha sido
empleado para caracterizar la disponibilidad hdrica de las plantas y la aridez y la productividad primaria neta

de los ecosistemas. De estudios realizados por diversos autores se ha encontrado que para B=1 (que
denominaremos Bc: Budyko Crtico) se presentan condiciones climatolgicas, geomorfolgicas y de
biodiversidad nicas que hacen que esta relacin adimensional sea de especial inters en la
hidroclimatologa y otras ramas de la ciencia. Conociendo el gran inters que despierta esta relacin
climtica, sta es estimada tanto para las diferentes cuencas como para Colombia y Antioquia (pixel a pixel).
En la Figura 2 se observa el mapa del nmero climtico de Budyko para Colombia estimado con el mapa de
PET de Cenicaf y el mapa de P interpolado por Kriging con deriva ENA (UNALMED-UPME (1999)).

Figura 2.

Nmero climtico de Budyko (P/PET) estimado con la evapotranspiracin potencial de Cenicaf y la

precipitacin de Kriging con Deriva ENA (UNALMED-UPME (1999)). La resolucin del mapa es de 5
minutos de arco.

Se observa que las regiones con menores valores de B corresponden a la costa Caribe (a excepcin de
la Sierra Nevada de Santa Marta) y algunas regiones ubicadas en las partes altas de las cordilleras. Los
mayores valores de B se observan en el litoral Pacfico, donde la precipitacin es alta, y en la cuenca
Amaznica, donde tambin se registran valores altos de pluviosidad. En general predominan los valores de

B superiores a 1 (Bc), indicando, como ya se haba encontrado, que nuestras cuencas son, en su mayora,
zonas limitadas por energa.
En general, son muy pocas las regiones con valores de B cercano al crtico (esto se obtuvo para la
mayora de las metodologas de estimacin de PET). Las regiones en las que se presentan valores de B
(estimado con PET de Cenicaf y P de Kriging con deriva ENA) cercanos al crtico (entre 0.8 y 1.2) se
encuentran concentradas en la costa Atlntica en cercana a las estribaciones de las cordilleras occidental y
central. Se observan algunas zonas de poca extensin en las partes altas de la cordillera central. Para las
metodologas de Penman y Morton se observa una regin de gran extensin en la Orinoqua con valores de
B cercanos a 1.
Un estudio bien interesante sobre este tema consistira en analizar las zonas de vida de estas regiones y
definir su estructura y caractersticas (flora, fauna, geomorfologa, etc) para encontrar relaciones con B, que,
como ya se ha dicho, es una sntesis del balance hdrico y energtico. Desafortunadamente, en nuestro pas
no se cuenta con informacin suficiente, detallada y de calidad al respecto, por lo que en este trabajo no se
pudo ahondar en la identificacin de tales relaciones.

4.4

Conclusiones

Se estim la evapotranspiracin real mediante la metodologa de Zhang et al (2001). Los resultados en la

estimacin de los caudales medios anuales muestran que es una metodologa aceptable para la estimacin
de ET. Esta metodologa es interesante ya que involucra el papel de la vegetacin en el proceso de la
evapotranspiracin. Sin embargo, es necesario corroborar mejor las hiptesis de esta metodologa para las
condiciones colombianas.
Se verific la aproximacin PETRn/Lv para las condiciones climticas colombianas. Se encontraron
errores entre el 20% y el 40%, lo que indica que esta aproximacin puede ser empleada aceptablemente en
regiones en las que no se conozca la evapotranspiracin potencial.
Se estim la relacin adimensional de Budyko. Se observa que en Colombia predomina la vegetacin
limitada por energa (valores de B superiores a 1).

REFERENCIAS
ABRAHAMS, A. D. Channel networks: a geomorphological perspective. Water Resources Research,
20:161-188, 1984.
BRUTSAERT, W. (1982). Evaporation into the atmosphere, theory, history and applications. Boston:
Klower Academic Publishers.
BUDYKO, M. (1958). The heat balance of the earths surface. Translate from russian by N.A Stepanova,
U.S Department of Commerce., Washington D.C, 259 pp.

BUDYKO, M. (1974). Climate and life. Translate from russian by D.H Miller. Academic. San Diego,
California.
EAGLESON, P. (1982). Ecologycal optimally in water-limited natural soil vegetation systems: theory and
hypotesis. Water Resources and Research, 18, 325-340.
KOSTER, R.D., P.A. DIRMEYER, P.C.D. MILLY, Y G.L. RUSSELL. Comparing GCM-generated land surface

waster budgets using a simple common framework. In Land Surface Hydrology, Meteorology, and Climate:
Observations and Modeling, Water Sciences and Application, Vol. 3, pp. 95-105. American Geophysical

Union, 2001.
MADDUMA-BANDARA. Drainage density and effective precipitation, Journal of Hydrology, 21,187-190,
1974.
MELTON, M. A. Correlation structure of morphometric properties of drainage systems and their controlling
agents. Journal of Geology, 66, 442-460. Chicago, 1958.
MILNE, B.T, V.K GUPTA Y C. RESTREPO. A Scale Invariant Coupling of Plants, Water, Energy and
Terrain. Ecoscience, 9, 191-199, 2002.
MORTON, F. (1982). Operational estimates of areal evapotranspiration and their significance to the
science and practice of hydrology. En N. H. R. Institute (Ed.), Journal of Hydrology, chapter 3, Otawa:
Environmental Canada, 1st edition, 237-263.
NEMANI, R, C.D. KEELING, H. HASHIMOTO, W.M. JOLLY, S.C PIPER, TUCKER, C.J. Climate-driven
increases in global terrestrial net primary production from 1982 to 1999. Science, 300, 1560-1563, 2003.
OJIMA, D.S., PARTON, W.J, COUGHENOUR, M.B Y SCURLOCK, J.M. impact of climate and
atmospheric carbon dioxide changes on grasslands of the world. Global Change: Effects on Coniferous
Forests and Grasslands. Scientific Committee on Problems of the Environment (SCOPE), Chichester, NY,
USA, Wiley, 271309, 1996.
OLFF, H., M. E. RITCHIE, Y H. H. T. PRINS. Global

environmental controls of diversity in large

herbivores. Nature, 415,901-904, 2002.

SANKARASUBRAMANIAN, A. Y VOGEL, R. (2002). Annual climatology of the united states. Water
Resources and Research, 38, 19.1-19.12.
THORNTHWAITE, C.W, The climates of the Earth. Geogr. Rev, 23:433-440, 1933.
UNALMED-CTA (2001). Balance Hidrolgico y Atlas digital de Antioquia. Universidad Nacional de
Colombia Sede Medelln-Convenio Interinstitucional Ctedra del Agua. Informe final.
UNALMED-UPME (1999). Balances Hidrolgicos de Colombia. Universidad Nacional de Colombia Sede
Medelln-Unidad de Planeacin Minero Energtica. Informe final.
UNALMED-UPME (2000). Atlas Hidrolgico de Colombia. Universidad Nacional de Colombia Sede
Medelln-Unidad de Planeacin Minero-Energtica. Informe final.
UNESCO. Mapa Mundial de Zonas ridas. UNESCO, 1979.
ZHANG, L., DAWES, W., Y WALKER, G. (2001). Response of mean annual evapotranspiration to
vegetation changes at catchment scale. Water Resources and Research, 37, 701-708.