UNIVERSIDAD NACIONAL DE COLOMBIA

Facultad de Ingeniería.
Dpto. de Ingeniería Civil y Agrícola

ASIGNATURA: HIDROLOGÍA

Andrés Vargas Luna

PROYECTO : CUENCA RÍO LA VIEJA

Presentado por:

Nicolás Ochoa Valencia

Nicolás Daniel Cifuentes Cruz
Jose Alejandro Gonzalez Ortiz
Wendy Lizette Martinez Garcia

BOGOTÁ D.C.

23 de enero de 2020
 CONTENIDO
1. Objetivos
1.1 Objetivo General
1.2 Objetivos Específicos
2. Datos a Utilizar
2.1 Descripción de la Cuenca.
2.2 Modelo Digital de Terreno -MDT
2.2.1 MDT a utilizar
2.2.2 Definición de la cuenca sobre el MDT
2.2.3 Importación del límite de la cuenca y la red de drenaje al SIG
2.2.4 Procesamiento del MDT
2.2.5 Mapa en el informe
2.3 Comparación de Modelo Digital de Terreno con Planos Físicos
3. Análisis y Resultados
3.1 Caracterización Geométrica de la Cuenca.
3.2 Caracterización de la Red de Drenaje
3.3 Relieve de la Cuenca
3.4 Análisis los Datos de Precipitación
3.4.1 Llenado y recortado de la información
3.4.2 Información de precipitación anual
3.4.3 Información de precipitación mensual
3.4.4 Mapas de precipitación media
3.4.5 Precipitación promedio multianual anual y mensual al interior de la cuenca
3.5 Evaporación
3.5.1 Temperatura promedio multianual anual en cada estación
3.5.2 Índice de aridez
3.5.3 Evaporación potencial
3.5.4 Evaporación real
3.6 Balance Hidrológico
3.6.1 Datos de caudal de la estación
3.6.2 Balance hidrológico de largo plazo en la cuenca
3.6.3 Balance de agua en el suelo y cálculo de caudales
3.7 Conclusiones
4. Bibliografía
LISTA DE ANEXOS
1. Objetivos

1.1 Objetivo General

Realizar un estudio geográfico e hidrológico de la cuenca del río La Vieja para poder calcular
su régimen de precipitación, evaporación, temperatura y caudal, con el fin de comparar los
resultados teóricos con los experimentales.

1.2 Objetivos Específicos

● Desarrollar y analizar todos los componentes de caracterización geomorfológica de la

cuenca por medios virtuales y con planos físicos.
● Desarrollar cálculos y análisis necesarios para determinar el ciclo anual de lluvias y de
evaporación en la cuenca a partir de datos de precipitación y temperatura.
● Comparar los resultados del ciclo de caudal mensual multianual calculado (obtenido
con el modelo de balance), con los datos reales de caudales multianuales mensuales
registrados.

2. Datos a Utilizar

2.1 Descripción de la Cuenca.

La cuenca asignada cierra en la estación limnigráfica de Cartago cuyas coordenadas son:

4,757611° latitud; -75,89925° Longitud y tiene una altitud de 937 msnm.

La extensión de la cuenca La Vieja abarca la mayoría del departamento del Quindío, además
de algunos municipios y veredas aledañas en el Valle del Cauca y Risaralda. Los centros
urbanos más importantes que componen la cuenca son Armenia, Buenavista, Calarcá
(Quindío), Taguales, Corozal (Valle del Cauca) y Pereira (Risaralda).

Las elevaciones de la cuenca se enmarcan entre 926 y 4795 msnm; la mayor elevación se
encuentra entre los picos de las montañas andinas y la menor se ubica en el encuentro del río
Cauca con el río La Vieja. Se aumenta la altura al ir del noroeste al sureste.

Las variaciones geológicas de la cuenca se dividen en:

UNIDADES DE SUELOS TRANSPORTADOS: Aquí se encuentran depósitos volcánicos (flujos

piroclásticos) y depósitos fluviales aluviales (debido a glaciares y cauces antiguos y viejos).
Los depósitos piroclásticos son predominantes en la parte noroeste de la cuenca, mientras que
Los suelos aluviales son predominantes en el costado oeste de la cuenca.
ROCAS DURAS: Suelos compuestos por secuencia interestratificada de rocas metamórficas y
cuerpos Ígneos intrusivos. Tienen poca frecuencia en la cuenca, imponiéndose básicamente de
pequeños depósitos en el centro de esta

ROCAS INTERMEDIAS TIPO 1: Suelos compuestos por secuencia interestratificada de rocas

de origen volcano-sedimentario, con efectos de metamorfismo dinámico. Estos suelos están
posicionados dentro de la cuenca en una grande hilera diagonal de este a suroeste.

ROCAS INTERMEDIAS TIPO 2: Suelos compuestos por secuencia interestratificada de rocas

metamórficas (esquistos micáceos y anfibolitas). Estos suelos están posicionados dentro de la
cuenca en una grande hilera diagonal en sentido suroeste.

ROCAS INTERMEDIAS TIPO 3: Suelos compuestos por secuencia interestratificada de rocas

metamórficas (esquistos cuarzo sericíticos, grafíticos, filitas verdes). Estos suelos están
posicionados dentro de la cuenca en una grande hilera diagonal en sentido sureste.

ROCAS BLANDAS: Secuencia interestratificada de areniscas, conglomerados y arcillolita. Esta

secuencia geológica se ubica con poca presencia al costado izquierdo de la cuenca.

SUELOS RESIDUALES Y SAPROLITO: (pequeños depósitos en la parte central de la cuenca)

(INGEOMINAS, 2005)
 Imagen 1. Cuenca enmarcada en Google Earth.Se aprecia los departamentos y el relieve del lugar

2.2 Modelo Digital de Terreno -MDT

Para abarcar la cuenca completa en una escala adecuada se adquirieron 5 planchas tanto de
forma física como de MDTs (por la página asf.alaska​). Las planchas utilizadas son: la 224, 225,
243, 244 y 262, entregadas por el IGAC. Todas cuentan con escala 1:100.000 y tienen
coordenadas Magna Sirgas 4686.

Para adquirir los MDTs se siguió paso a paso la guía brindada en clase. Primero se accede a
la página ASF Data Research Vertex y se introducen dos coordenadas que aproximan el área
de la cuenca, a continuación se busca y descarga la imagen satelital que se ajuste mejor a la
cuenca. Se guardan los archivos descargados, se le agrega el tipo de coordenadas y se abre el
DEM en ArcMap. Se genera un polígono en Excel con los puntos antes utilizados, el cual va a
recortar el DEM. El siguiente paso a es generar el flujo con las herramientas “Spatial analyst
tool” e “Hydrology”, y finalmente se delimita la cuenca seleccionando en la lista “Hydrology” la
opción “watershed”.
Una aclaración para la realización del proyecto es que se aumentó el detalle de la red hídrica
(Value >200.000), pues así el trabajo de Arcgis se ajustó correctamente a la información
contenida en las Planchas..

Imagen 2. Cuenca Enmarcada en Arcgis.Se aprecia los departamentos y el relieve del lugar

2.3 Comparación de Modelo Digital de Terreno con Planos Físicos

 Imagen 3. Cuenca Enmarcada en Arcgis y Planos de la cuenca del IGAC

Claramente se indica en las imágenes que el detalle del MDT es bastante acertado con
respecto al mapa físico, se puede ver que la silueta enmascara describe correctamente la
cuenca, la red hídrica no es exactamente la misma dado que el MDT hace una aproximación
con los datos de relieve que recibe para generar esta red, mientras que en el mapa físico
podemos apreciar la red real de ríos y quebradas que se encuentra dentro de la cuenca. Por lo
tanto se hace un ajuste de la red hídrica en el dem para que se parezca a la red hídrica del
mapa físico, sin embargo por más aproximación que se haga, siempre hay variación en las
líneas de flujo.

3. Análisis y Resultados

3.1 Caracterización Geométrica de la Cuenca.

El estudio de la geomorfología de la cuenca del río La Vieja se basó en el análisis del DEM
(​Digital Elevation Model)​ de la zona, usando el software ArcMap.

A continuación se presentan los datos y análisis de resultados de la longitud del cauce

principal, cálculo de la pendiente, densidad de drenaje, número del orden del cauce principal,
cálculo de la relación de bifurcación y la elaboración de la curva hipsométrica.
 Tabla 1. Características geométricas de la cuenca.

El área y perímetro de la cuenca fueron calculados al enmascarar la cuenca en Arcgis

Longitud axial de la cuenca

Este dato se calculó tomando la distancia entre la estación limnigráfica y el punto más lejano.
Cabe destacar que la longitud axial es menor a la longitud del cauce principal, ya que este
empieza desde el punto de concentración, baja hacia el centro de la cuenca y sube hacia la
parte nororiental de la misma.

Ancho de la cuenca

Se calculó dividiendo el área entre la longitud de la cuenca.

Índice de compactación de Gravelius

Índice utilizado para obtener la “circularidad” de la cuenca, un valor igual a 1 o cercano a este
indica una cuenca circular, en donde los tiempos de concentración desde cualquier punto de la
cuenca serán muy similares. La cuenca del río La Vieja cuenta con un Kc de 1,65 lo cual indica
que es alargada.
 Imagen 4. Forma de cuencas según el índice de Gravelius. Obtenida de
https://www.prontubeam.com/articulos/articulos.php?Id_articulo=26

Factor de forma

Se determina para saber si la cuenca es cuadrada, alargada o achatada. Una cuenca achatada
es más propensa a sufrir inundaciones, ya que con un cauce corto en caso de una creciente, el
cauce tiende a desbordarse ya que el tiempo de respuesta de la misma es muy bajo, el factor
de forma es igual al área de la cuenca dividida entre la longitud de la cuenca al cuadrado.

En la cuenca del río la vieja se obtuvo un valor de 0,446 que indica que la cuenca es alargada,
coincidiendo con la información suministrada por el índice de compacidad, por lo tanto se
esperaría que el tiempo de concentración de la misma sea alto.

3.2 Caracterización de la Red de Drenaje

Orden de la cuenca:

En la tabla 2 se encuentran los datos determinados a partir del esquema de ordenamiento de

corrientes de Strahler-Horton.Teniendo en cuenta la jerarquía de órdenes propuesta por
Horton, en cual los tributarios sin ramificaciones son de orden uno, el orden dos es una
corriente formada por corrientes de orden uno. En el caso de que se encuentren dos corrientes
de diferente orden la corriente generada tendrá el orden mayor entre la corrientes
generadoras.​Francisco.J.(1989)​.
 Imagen 5. Distribución órdenes de cauce.

Se tomaron los números de orden del cauce de la cuenca y se contaron, teniendo en cuenta
que dos cauces con el mismo número de orden consecutivos se cuentan como un solo cauce.
Además, se calculó la longitud y se consignaron los resultados en la siguiente tabla:
 Tabla 2. Ordenamiento usando DEM en software ArcMap

De esta forma, se determinó que el número de orden del cauce principal es de 7.

Este resultado nos muestra que la red de drenaje de la cuenca del río La Vieja es bastante
ramificada, aunque es necesario mencionar que este esquema de ordenamiento depende en
gran medida de la escala en la cual se encuentre representada la red de drenaje, siendo que a
mayor escala (en este caso manejamos una cartografía en escala 1:100000) se podrán
determinar mayor número de afluentes, caso que sería contrario si se trabajara en una escala
mucho menor.

Relación de bifurcación:

Para el cálculo de la relación de bifurcación de la cuenca se utilizó la siguiente relación

geométrica:

En este caso la relación de bifurcación es el promedio hallado para cada orden:

Tabla 3. Relación de bifurcación.

Además, se utilizó las leyes de similitud de Horton en donde la relación de bifurcación es
directamente proporcional a:

Este valor se confirmó gráficamente mediante la extracción de la pendiente del gráfico de

Relación de bifurcación, para así obtener el valor final de Rb. Lo anterior se presenta a
continuación:

Gráfico 1. Línea de tendencia relación de bifurcación.

Con la ecuación de la línea que mejor se ajusta a los datos, extraemos la pendiente de la
misma, la cual es el valor que acompaña a la x en la ecuación con un valor de -0,5386.

En donde la relación de bifurcación está dada por:

𝑅𝑏 = 𝑙𝑜𝑔−1(𝑏)

De acuerdo con los datos obtenidos, se determina una relación de bifurcación Rb de 3,46
 𝑅𝑏 = 𝑙𝑜𝑔−1(0,5386) = 3,46

Los resultados no presentan una diferencia significativa por lo que se confirma que el
procedimiento está bien hecho.

La relación de bifurcación es un parámetro que está influenciado por la forma de la cuenca y

que generalmente varía entre 3 y 5. El valor obtenido de 3,46 corresponde al hecho de que la
forma de la cuenca es algo alargada y a la vez ensanchada. En este sentido, este es un factor
que limita el tiempo de respuesta de la cuenca ante eventos de precipitación.

Relación de longitudes:

Para el cálculo de la relación de longitudes de la cuenca se utilizó la siguiente relación

geométrica:

En este caso la relación de longitudes es el promedio hallado para cada orden:

Tabla 4. Relación de longitudes.

Utilizando las leyes de similitud de Horton para la relación de longitudes se obtiene:

Este valor se confirmó gráficamente mediante la extracción de la pendiente del gráfico de
relación de longitudes, para así obtener el valor final de RL. Lo anterior se presenta a
continuación:

Gráfico 2. Línea de tendencia relación de longitudes.

Con la ecuación de la línea que mejor se ajusta a los datos, extraemos la pendiente de la
misma, la cual es el valor que acompaña a la x en la ecuación con un valor de 0,2525.

En donde la relación de longitudes está dada por:

𝑅l = 𝑙𝑜𝑔−1(𝑏)

De acuerdo con los datos obtenidos, se determina una relación de longitudes RL de 1,79

𝑅l = 𝑙𝑜𝑔−1(0,2525) = 1,79

La relación de longitudes es un parámetro que generalmente varía entre 1,5 y 3,5. El valor
obtenido varía entre 2 y 2,7 por los diferentes métodos hallado, corresponde al hecho de que la
forma de la cuenca es alargada.
Densidad de drenaje:

Una vez calculada la red de drenaje,se midió la longitud de cada afluente perteneciente a esta.
Para realizar la medición la cuenca se subdividió en zonas de las cuales se obtuvo la
respectiva longitud de la red de drenaje.

Donde Lu es la longitud total de los cauces de la cuenca y A es el área total de la cuenca, con
2
la medición se obtuvo un valor de Lu = 3106, 708 km y A = 2777, 01 km , obteniendo una
2
densidad de drenaje de 1, 12 km/km .

La densidad de drenaje de ​1,12 km/km​2 obtenida para la cuenca del río La Vieja, es un valor
de baja densidad de drenaje, generando volúmenes de escorrentía bajos y baja susceptibilidad
a las avenida torrenciales.

Los valores que indican una alta densidad de drenaje se presentan en áreas con presencia de
suelos relativamente impermeables, con pendientes fuertes y escasa cobertura vegetal;
mientras que, los valores que indican una baja densidad de drenaje son aquellos que se
presentan en áreas con suelos compuestos por materiales resistentes a la erosión o muy
permeables y su relieve es bajo (Linsley R. K., 1977).

Tiempo de concentración:

El tiempo de concentración es considerado como el tiempo de viaje de una gota de agua que
escurre desde el lugar más lejano de la cuenca hasta el punto de interés. Para su cálculo se
emplean diferentes fórmulas empíricas que relacionan otros parámetros de la cuenca ya
conocidos.( Ven Te Chow, 1994).

A continuación se presentan dos de los ecuaciones utilizadas para hallar el tiempo de

concentración de la cuenca del río La Vieja.

Fórmula de Kirpich

Donde:
L= Longitud del cauce principal (m)
S= Pendiente promedio del recorrido (adm)
Tc= Tiempo de concentración (min)
Fórmula de Passini

Donde:
Lc= Longitud del cauce principal (km)
So= Pendiente promedio del cauce principal (adm)
Tc= Tiempo de concentración (horas)
A= Área de la cuenca (km​2​)

Como se puede ver, las siguientes ecuaciones se correlacionan bastante en sus resultados
generando un promedio de tiempo de concentración de 36 horas. A pesar de esto, al comparar
estos resultados con un documento oficial de la CAR CVC sobre el cauce del río la Vieja se
aprecia que el tiempo de concentración varía notablemente ( según el método Kirpich 46,58
horas y según Passini 53,73 horas). Sin embargo, se concluye que la diferencia de tiempo se
debe a la diferencia del punto de salida tomado, pues en este documento se tomó la salida en
la estación de Cartago y en el documento oficial toma la salida en el encuentro con el río Cauca
(Corporación regional del valle del cauca, 2014)

3.3 Relieve de la Cuenca

Pendiente de la cuenca (utilizando el valor promedio del mapa de pendientes)

Imagen 6. Pendientes reclasificadas obtenida del DEM.
 Tabla 5. Cálculo de pendiente promedio de la cuenca.

Se generó un mapa de pendientes por medio de la herramienta Surface -> Slope, teniendo en
cuenta el archivo raster de las alturas en la cuenca. Luego se clasifica el mapa y se le crean 10
intervalos iguales. Para obtener la pendiente promedio se toma la información brindada en la
tabla de atributos y se genera un excel.

Curva hipsométrica

Para la obtención de la curva hipsométrica se usó información de las alturas en la cuenca y

usando herramientas de ARCGIS, se logró calcular un histograma de áreas por encima de las
determinadas alturas y dividirlas entre el área de la cuenca para expresarlas en forma de
porcentaje y de esta manera obtener la curva.

Se reclasifico la altura en Arcgis mediante el archivo raster que brindaba la elevación con la
herramienta Reclass. Luego calculamos el área con la herramienta de estadística, zonal.

Se generaron 12 intervalos iguales para el cálculo de las áreas.

 Gráfica 3. Curva hipsométrica.

Imagen 7. Tipos de curvas hipsométricas. Obtenida de

https://ingeciv.com/que-es-una-curva-hipsometrica/

Comparando la curva obtenida con la imagen 6, se puede observar que se trata de una curva
tipo C, característica de cuencas sedimentarias, de avanzada edad, por lo que se puede inferir
que la cuenca del río la vieja lleva una gran cantidad de tiempo siendo afectada por erosión y
otros factores.

Elevación media de la cuenca

Para hallar este dato se utilizó un promedio ponderado en el que se toma una cota de altura
dentro de la cuenca multiplicada por su correspondiente área .
Con este proceso la altitud promedio de la cuenca arroja un valor de 1817,9 msnm. Valor que
se confirma con el punto correspondiente al 50% en el eje x de la curva hipsométrica de la
cuenca.

Pendiente del cauce principal

Al conocer la red de drenaje del río La Vieja, se estableció el cauce principal como el afluente
de mayor longitud, como puede apreciarse en la ilustración 8, es el afluente resaltado por su
color azul claro. Luego, se midió la longitud del cauce principal obteniendo: 𝐿𝐶𝑃= 146 𝑘𝑚.

Imagen 8.Cauce principal y curvas de nivel.

Usando la fórmula de Taylor-Schwarz, se tomaron 25 tramos de igual longitud y se interpolo los
valores entre curvas de nivel para el cálculo de la pendiente y la gráfica de perfil del cauce
principal.

Gráfica 4. Perfil del cauce principal.

 Tabla 6. Pendiente del cauce principal por Taylor - Schwarz.

El valor de la pendiente es de 4,49% y se debe a que el cauce nace en la zona montañosa de

la cuenca y se suaviza hacia la región del valle del cauca, teniendo una zona de altas
pendientes y una zona de llanura.

3.4 Análisis los Datos de Precipitación

3.4.1 Llenado y recortado de la información

Tanto para hacer el llenado de datos de precipitación como de temperatura se eligió una
ventana temporal de 1980 a 2015, ya que es un periodo en común en el que las estaciones
escogidas para el análisis tienen suficientes datos. También se realizó el método de polígonos
de Thiessen para determinar cuáles estaciones tienen influencia sobre la cuenca; las
estaciones que no tenían influencia fueron obviadas, de modo que no se llenaron sus datos ni
se usaron para interpolar.

Para completar los datos se realizó una Interpolación con el método IDW, pues brinda unos
datos muy acertados con pocas estaciones. Este método depende de la distancia entre
estaciones y el valor de precipitación (o dado el caso, temperatura), así que es necesario usar
las coordenadas de las estaciones. Ante la insuficiencia de datos, en este llenado de datos se
tuvo que completar toda la información de un año en diferentes en estaciones.
(Ver Anexo llenado_datos)

3.4.2 y 3.4.3 Información de precipitación mensual y anual

En la siguiente tabla se puede observar la precipitación mensual multianual (mm) y la

precipitación anual multianual (mm) que registraron la estaciones con influencia en la cuenca.

Tabla 7. Promedio de precipitación mensual y anual multianual para cada estación.

Para obtener esta información se usaron el llenado de datos de las estaciones, luego se calculó
el promedio de los datos de precipitación mensual desde los años 1980 a 2015, generando así
el promedio mensual multianual. luego se suman las precipitaciones anteriores y se genera el
promedio anual multianual de cada estación

3.4.3 Información de precipitación mensual

Al analizar los datos de las estaciones y las posiciones de las mismas, se llega a la conclusión
de que la estación que mejor describe la precipitación dentro de la cuenca es Alcalá, pues está
posicionada en un punto central entre las todas las estaciones. La información correspondiente
a su promedio anual multianual y su promedio mensual multianual se encuentra dentro de la
tabla 7.

El régimen de precipitación en la cuenca es un régimen bimodal como se puede apreciar en la

gráfica 5, ya que hay altos registros de precipitación en 2 épocas diferentes del año. El régimen
dentro de la cuenca cuenta con 2 picos en abril y octubre; y 2 valles en enero y agosto.

Gráfica 5. Régimen de precipitación de la cuenca.

3.4.4 Mapas de precipitación media ​ (VER ANEXO: MAPAS DE

PRECIPITACIÓN)

Para realizar los mapas de interpolación en Arcgis se usó la herramienta Toolbox-> Spacial
analisis->interpolate->IDW. En este programa se indica que en la interpolación se deben tener
en cuenta todas la estaciones que tengan influencia por los polígonos de Thiessen, luego se
enmascara la cuenca para solo observar la interpolación dentro de la cuenca.

En los mapas de precipitación, se observa por la tonalidad azul más fuerte, que los meses de
octubre y noviembre tuvieron los mayores valores mensuales multianuales de precipitación y en
menor manera el mes de abril también tuvo altas precipitaciones, en contraste con los meses
de junio, julio y agosto que tienen una tonalidad más verdosa representando precipitaciones de
menor magnitud.

El mapa de precipitación anual multianual nos permite apreciar que la precipitación está más
concentrada en la parte nororiental de la cuenca, donde se ve un tono azul muy fuerte, la
estación ubicada en este punto, Salento, fue la que registró los mayores valores de
precipitación en general, hay otros dos puntos en donde se concentran los registros altos de
precipitación más hacia el sur de la cuenca en la ubicación de dos estaciones, y hacia el sur de
la cuenca se observan los menores valores.

Al estar el punto de mayor concentración de lluvia relativamente alejado del punto de cierre de
la cuenca, esto podría influir en el tiempo de respuesta de la cuenca a eventos de precipitación,
ya que el mayor volumen de precipitación que cae sobre la cuenca debe recorrer una distancia
importante para llegar al punto de cierre.

3.4.5 Precipitación promedio multianual anual y mensual al interior de la cuenca

Al generar los mapas raster a partir de los datos de precipitación, en las estadísticas de estos
mapas se pueden leer los valores promedio que saca ArcGis de manera automática.

(VER ANEXO: Balance Hídrico)

3.5 Evaporación

3.5.1 Temperatura promedio multianual anual en cada estación.

Para cada una de las 8 estaciones con influencia en la cuenca con información de temperatura
disponible, se calculó el promedio de la temperatura para cada mes entre los años 1980 a
2015, así como el promedio anual multianual. Esta información se ordena en la siguiente tabla.

Tabla 8. Datos de temperatura multianual de la cuenca.

3.5.2 Índice de aridez

El índice de aridez se calculó usando una relación entre las ecuaciones de Budyko y Turc.
Usando la ecuación de Turc, se reemplazaron los datos de temperatura media mensual y anual
multianual para encontrar el valor del parámetro L en cada estación, el cual sería reemplazado
en la fórmula donde se relaciona con los datos de precipitación mensual y anual multianual.

🇷
En este apartado cabe aclarar que para el análisis mensual, se tuvo que multiplicar los valores
de precipitación mensuales registrados por 12 y al final se dividió el valor de E /P encontrado
entre 12 para ajustar los valores a la ecuación, ya que esta está pensada para el cálculo anual.

🇷
🇷
Con el valor del factor E /P encontrado, se usa la herramienta buscar objetivo de EXCEL para
encontrar un valor de índice de aridez Φ que ajuste la ecuación de Budyko al valor de E /P de
la ecuación de Turc.

Los resultados se consignan en la siguiente tabla.

Tabla 9. Índice de aridez calculado.

Como se puede observar en la tabla, el valor del índice de aridez en todas las estaciones se
muestra en su máximo valor en dos temporadas: de Diciembre a Enero y de Julio a Agosto. Así
se muestra nuevamente el régimen bimodal de precipitación de la zona, donde los meses con
alto índice de aridez, tienen una alta temperatura y baja precipitación generalmente.
Con esta información, se realizaron los 12 mapas mensuales y el mapa anual de índice de
aridez en el periodo entre 1980 y 2015 para la cuenca usando ARCGIS.

3.5.2.1 Mapas de índice de aridez ​(VER ANEXO: MAPAS ARIDEZ)

Los mapas generados con los valores del índice de aridez, indican que los meses de junio, julio
y agosto son los meses más áridos en promedio durante el año, meses en los cuales los
registros de precipitación muestran los valores más bajos. El mapa anual multianual nos da una
mejor idea de que la concentración de los valores altos del índice, se da hacia la parte
occidental de la cuenca, contrastando así, con los mapas de precipitación en los que la
concentración se daba hacia la parte nororiental de la misma.

Estas claras diferencias dentro de la cuenca podrían estar relacionadas a la geografía de la

cuenca y a las actividades económicas enmarcadas en las zonas, ya que la cuenca está en el
límite entre varios departamentos, y en las zonas en las que se analiza en este contexto. La
parte más árida se da hacia la región del Valle del Cauca y la parte más lluviosa hacia la región
del Quindío. Esto podría estar relacionado a que parte del Quindío tiene región montañosa
donde las temperaturas son más bajas y en a parte de la cuenca que hace parte de El Valle del
Cauca, como su nombre lo indica, la región es muy plana, dando lugar a temperaturas más
altas.

3.5.3 Evaporación potencial

Una vez obtenido el valor del índice de aridez para cada estación a nivel mensual y anual, se
multiplica este valor por la respectiva precipitación registrada por la estación. Este valor es un
indicativo de la evaporación máxima teórica que se podría dar en la cuenca bajo condiciones
óptimas.

Los resultados de los cálculos se muestran en la siguiente tabla.

Tabla 10. Evaporación potencial calculada

Cabe resaltar que en varios casos la evaporación potencial calculada, excede por mucho el
valor de precipitación registrado en la estación, esto podría explicarse ya que las estaciones
registraron altas temperaturas y muy bajas precipitaciones, haciendo que el valor del índice de
aridez aumenta bastante.

3.5.3.1 Mapas de evaporación potencial ​(VER ANEXO: MAPAS

EVAPORACIÓN POTENCIAL)

En los mapas de evaporación potencial, se puede ver un fenómeno similar a lo observado en

los mapas de índice de aridez, donde los valores más altos de evaporación potencial se dan
hacia el occidente de la cuenca (en la región del Valle del Cauca). Esto, concuerda
perfectamente con lo observado en los mapas de índice de aridez, ya que en las zonas más
áridas, se dan las temperaturas más altas y por lo tanto hay más probabilidad de evaporación
hacia estas zonas. En los mapas se muestra un fenómeno interesante, ya que hacia el sur de
la cuenca se observan los menores valores de evaporación potencial; sin embargo, esto puede
deberse a que las estaciones de registro de temperatura están ubicadas hacia el norte de la
cuenca, Esto crean así una carga hacia esa dirección, por lo que para describir de una manera
más acertada las condiciones de la cuenca, se hacen necesarias estaciones ubicadas hacia el
sur de la misma.

Aunque no es tan claro como en los demás mapas, hacia los meses de junio, julio, agosto y
enero se ven los valores más altos de evaporación potencial, concordando con las épocas de
menor precipitación del régimen de la cuenca.

3.5.4 Evaporación real

El valor de evaporación real se halló en el proceso para calcular el índice de aridez aplicando
las ecuaciones de Budyko y Turc vistas con anterioridad, para estos cálculos se usaron los
datos de precipitación y temperatura multianual mensual y anual.
Se consignan los resultados en la siguiente tabla.

Tabla 13. Evaporación real calculada

Si se comparan estos valores con los de precipitación registrados por estas estaciones se
puede observar que la evaporación real siempre es menor a la precipitación. En los casos en el
que se registró una temperatura alta y una baja precipitación el valor de evaporación real
calculado es casi igual al de precipitación.

Tabla 13. Precipitación mensual multianual y anual multianual registrada en las estaciones con datos de
temperatura.

3.5.4.1 Mapas de evaporacipon real

Con los mapas de evaporación real de la cuenca se tiene información transpuesta a la

suministrada por mapas de evaporación potencial, donde, los meses con valores más altos
corresponden a los más áridos de la cuenca, es decir, enero, junio, julio y agosto. Sin embargo,
la evaporación real calculada muestra que en estos meses se dieron los valores más bajos de
evaporación, y los meses pico en el régimen de precipitaciones, noviembre y abril, presentaron
los valores más altos de evaporación.

Esto puede deberse a que en estos meses al haber más agua disponible, lo que también
aumenta el ancho de los ríos y el área superficial de los mismos, se disponía de más agua para
evaporar que en los meses más áridos. Es decir, aunque los mapas de evaporación potencial
indican que la temperatura era el mayor factor para la evaporación, los valores reales
calculados muestran que en los meses en los que hay más precipitaciones y por lo tanto más
agua disponible se presenta la mayor evaporación.

Esto también se puede explicar viendo la comparación entre los valores de precipitación
registrados y de evaporación real calculada, que son muy similares en los meses menos
lluviosos. Aquí un gran porcentaje del agua que cae en forma de precipitación es evaporada, y
aunque el potencial de la cuenca para evaporar el agua sea muy grande, no se tiene agua
disponible para hacerlo.

El mapa anual multianual, muestra una concentración de evaporación muy marcada hacia el
norte de la cuenca, disminuyendo gradualmente hacia el punto de concentración en el sur, pero
al igual que la región occidental de la cuenca, el norte de la misma carece de zonas
montañosas. Por lo que la suposición de las altas temperaturas debido a los valles sigue siendo
válida.
 3.6 Balance Hidrológico

3.6.1 Datos de caudal de la estación

Para el estudio de la cuenca del río la Vieja, se usó como estación de cierre, la estación
Limnigráfica Cartago, y al igual que con las estaciones de medición de precipitación, se escogió
una ventana temporal de 1980 a 2015.

A estos datos se les calculó el caudal promedio multianual mensual y anual, datos que se
usaron para analizar el régimen de la cuenca.

Tabla 14. Registro de caudales de la estación Cartago

Con los datos de caudal promedio mensual multianual, se procede a realizar la gráfica 6, donde
se puede observar un claro régimen bimodal de caudal, con picos en los meses de mayo y
noviembre, y cuyos valores más bajos se dan en el mes de agosto.

Si se compara este régimen, con los valores de índice de aridez de la cuenca, se evidencia que
el mes de agosto tiene unos altos valores de aridez. lo anterior se contrasta con los bajos
valores del mismo índice en los meses de mayo y noviembre donde se dan los picos de caudal
de la cuenca. Indicándonos que la incidencia de la temperatura y precipitación se relacionan
directamente con la respuesta de la cuenca.

Gráfica 6. Caudal promedio mensual multianual.

Los datos de caudal fueron comparados con el histórico de eventos de temperatura (Climate
Prediction Center Internet Team, 2017) ​donde se resaltan los meses con temperaturas
anormalmente altas o bajas. Se superpuso estos eventos con los registros de caudal, y así se
da una buena idea la ocurrencia de los fenómenos microclimáticos conocidos como NIÑO y
NIÑA.
 Tabla 15. Comparación de fenómenos extremos de temperatura y registro de caudales.

En la tabla 15, se resaltaron con azul claro los meses en los que fue registrada una
temperatura anormalmente baja, y en naranja los meses en los que se registró una temperatura
anormalmente alta. En la última columna a la derecha, correspondiente al caudal promedio
anual, se usó un formato de color condicional para resaltar en azul los caudales más altos y en
rojo los más bajos.

El formato de colores usado permite evidenciar que en los años en los que se registraron
temperaturas muy altas, los caudales son generalmente bajos; en contraste, los años en los
que la temperatura fue baja, coinciden con caudales altos. Se puede observar que en cierta
forma que los fenómenos de temperaturas extremas se dan en periodos de aproximadamente
dos o tres años en esta región, y que en cierta forma se alternan temporadas de altas
temperaturas con temporadas frías.
Para analizar esta información de mejor manera, se discriminaron los años con temperaturas
predominantemente altas (NIÑO), los de temperaturas bajas (NIÑA) y los que comparten
ambas características en proporciones similares o simplemente no predomina ninguna. El valor
de caudales altos y bajos se usó como referencia para discriminar los años en los que se
presentaban dudas.

Tabla 16. Discriminación de años NIÑO, NIÑA y Normales.

Con la información de caudales mensuales multianuales promedio en los años discriminados,
se procedió a graficar estos valores para compararlos.

Gráfico 7. Caudales promedio mensuales en años NIÑO, NIÑA y Normales.

Este gráfico es un indicativo del efecto de los fenómenos del NIÑO y la NIÑA en Colombia,
donde evidentemente, los años NIÑA, se caracterizan por unos altos caudales y unas bajas
temperaturas, lo que se podría asociar a una alta precipitación en la zona. Por el contrario, los
años NIÑO, tienen unas altas temperaturas, lo que podría indicar muy bajas precipitaciones y
terminando así en caudales bajos. Se puede destacar que independientemente del fenómeno,
el régimen de la cuenca sigue manteniendo su forma bimodal, con picos en los meses de
Abril-Mayo y Noviembre-Diciembre.

3.6.2 Balance hidrológico de largo plazo en la cuenca

Tomando los valores de precipitación y evaporación real medias a nivel mensual multianual y
anual multianual calculadas mediante la interpolación de los mapas en ArcGIS, se procede a
realizar un balance hídrico simple para estimar los valores de caudal esperados en la cuenca,
asumiendo que el almacenamiento dentro de la cuenca permanece constante, e ignorando la
infiltración e intercambios con cuencas vecinas se tiene:

P-E=Q

Los datos de precipitación y evaporación calculados mediante la interpolación de los mapas, se

expresan en las tablas con unidades de m³/s para poder compararlas con los registros que se
tienen de la estación limnigráfica de cierre de la cuenca.
 Tabla 17. Balance hídrico a nivel mensual y anual y comparación con registros.

En la tabla 17, se usó la fórmula mencionada previamente para calcular el caudal mensual
multianual y anual multianual, y a su vez se comparó con los datos de caudales medios
mensuales multianuales extraídos de los registros de la estación Cartago, mostrando en la
última fila el error entre ambos caudales, para facilitar el análisis, se grafican ambos caudales.

Gráfica 8. Comparación de caudal calculado y real.

Como evidencian la tabla 17, y la gráfica 8, en algunos meses se tuvo un error de menos del
30% pero en otros, principalmente octubre, se observan grandes diferencias en los valores
calculados y reales. Se analiza la gráfica hasta el mes de mayo y se tienen valores muy
cercanos, pero cuando empieza la época de bajos caudales, el caudal esperado es mucho
menor que el real, ocurriendo lo contrario en la temporada inmediatamente posterior, donde el
caudal esperado supera por mucho el caudal real.

Los variaciones entre los caudales esperado y real, pueden deberse a varios factores, pero
principalmente se podría explicar debido a que en el balance hídrico se ignoró por completo la
infiltración en la cuenca. Si se observa la gráfica, en el primer pico del año se tiene una buena
aproximación de los datos, lo que podría deberse a que la capacidad de infiltración del suelo se
superó y casi en su totalidad la precipitación se convirtió en escorrentía superficial, ahora bien,
en los meses valle de la gráfica, se calculó un valor mucho menor al registrado, pero siguiendo
con la idea de la infiltración se puede explicar perfectamente.

Al no haber una gran cantidad de precipitación, el agua infiltrada en el suelo de la cuenca

empieza a hacer aportes al caudal de la misma, explicándolo de otra forma, el agua que fue
almacenada en el suelo durante la época más lluviosa, sale unos meses después debido a que
su tiempo de concentración es mucho menor al flujo superficial, dando como resultado un
retraso en la respuesta de la cuenca, así que el caudal real excede al esperado ya que está
siendo alimentado por agua almacenada en la cuenca de periodos anteriores de precipitación.

Justo después de esta temporada “seca”, en el segundo pico del año, el acuífero que ha
perdido el agua contenida en él empieza otra vez a recargarse en la temporada de lluvias, por
lo que en el segundo pico del año, el caudal real es menor al calculado, ya que una parte de la
precipitación se pierde en infiltración y no se convierte directamente en escorrentía superficial.

Esta explicación gana aún más fuerza analizando los valores anuales multianuales de caudal,
ya que solo difieren en un 4.55%, esto es debido a que en la época de bajas precipitaciones se
subestima el caudal ya que se ignora el aporte del agua subterránea, y en el segundo pico se
sobrestima el caudal ignorando la infiltración, por lo que se compensa de alguna manera el
caudal, dando como resultado una gran similitud entre el caudal esperado y real.

3.7 Conclusiones

La cuenca estudiada, por su magnitud, presenta condiciones muy diferenciadas en varias

zonas, como lo pueden indicar los mapas del índice de aridez de la misma, teniendo en cuenta
que la cuenca está entre varios departamentos, se recomienda su estudio en zonas más
pequeñas para hacer una caracterización más precisa de la misma.

Por la geometría de la cuenca, la respuesta de la misma a los eventos de precipitación no es

inmediata pero sí es bastante rápida, esto es evidenciado por los regímenes de precipitación y
caudales, los cuales tienen los picos en las mismas épocas del año, mostrando que la
respuesta de la cuenca a estos eventos puede ser cuestión de horas en los eventos más
extremos.

El análisis a nivel anual, revela que la cuenca se ve afectada por los eventos macroclimaticos
de una gran manera durante ciclos de años, en los que al ser afectada por el fenómeno del
NIÑO disminuye en gran proporción el caudal de la misma y contrariamente en el fenómeno de
la NIÑA se pueden dar crecientes considerables en el cauce principal del río la Vieja.

Los datos analizados de caudal de la estación son una clara muestra de la importancia de la
infiltración de la misma y de los suelos que la componen, ya que los tipos y usos que se dan al
suelo afectan muy directamente la respuesta de la cuenca a eventos de precipitación,
indicando en este caso que en las épocas más secas del año el aporte de agua subterránea se
hace importante en el estudio del caudal.

4. Bibliografía

Chow, V. T. (1977). ​Applied Hydrology.​ McGraw-Hill. 1977

Climate Prediction Center Internet Team. (2017). ​National Weather Service Center.
Obtenido de National Weather Service Center:
https://origin.cpc.ncep.noaa.gov/products/analysis_monitoring/ensostuff/ONI_v5.
php

Córdova, M. (30 de abril de 2016). ​prontubeam. Obtenido de prontubeam:

https://www.prontubeam.com/articulos/articulos.php?Id_articulo=26

Corporación regional del valle del cauca (2014). La Vieja POMCA en Ajuste Fase
Diagnóstico. Capitulo 7 Morfometría

ingeciv. (9 de febrero de 2015). ​ingeciv. Obtenido de ingeciv:

https://ingeciv.com/que-es-una-curva-hipsometrica/

INGEOMINAS. (2005). ​Resumen de la Geología del Departamento del Quindío. Bogotá

D.C.

IGAC (2019). planchas: 224, 225, 243, 244, 262. Escala 1:100.000

​ insley, R. K. (1977). Hidrología para Ingenieros. Bogotá: McGraw Hill Latinoamericana

L
S.A.
 LISTA DE ANEXOS

● ·​ ​Geometría de la cuenca
■ o ​Caracterización geométrica.xlsx

Archivo de Excel en donde se encuentran los cálculos de las

características geométricas de la cuenca (área, pendientes, curva
hipsométrica, etc.).

● ·​ ​llenado de datos.xls
■ Archivo de Excel que muestra cómo se realizó la interpolación para poder
completar datos de temperatura y precipitación dentro de las estaciones.
También tiene información sobre precipitación mensual y anual
multianual e información sobre los fenómenos del niño y la niña
● ·​ ​Mapas interpolados

■ o​ ​Mapas de precipitación
■ o​ ​Mapas de evaporación potencial
■ o​ ​Mapas de evaporación real
■ o​ ​Mapas de índice de aridez
● ·​ ​Maps temp

■ o​ ​Evaporación.xls
■ Archivo de Excel con las tablas de resultados de los cálculos de
evaporación real y potencial.
■ o​ ​Aridez.xls
■ Archivo de Excel que muestra los cálculos para obtener el índice de
aridez de cada cuenca
● ·​ ​Balance hídrico.xls

○ Archivo de Excel en donde se obtiene el valor de precipitación, índice de aridez y

evaporación mensual multianual brindado por Arcmap. En este archivo también
se hace una comparación entre el caudal real y el caudal teórico