UNIVERSIDAD NACIONAL DE COLOMBIA

Sede Bogotá

HYDROLOGY
2024-2S

Group # 05

Submitted to
Professor Nestor Mancipe

Prepared by:

Santiago Cortes Delgado Andres Camilo Amado Chacon

Jorge Esteban Cantor Bernal Jean Carlo Vargas Avila

Due date: 22 January 2025

Score:

Note:
My signature on this document indicates that: [1] I have done my fair share of this
assignment; [2] I have discussed in detail all problems and solutions of my teammates; [3] I
have checked my teammates work for soundness and accuracy; and [4] I am competent to
solve any problems handled by my teammates.
 Estudio hidrológico de la cuenca del río Saldaña en las en las subzonas hidrográficas
de Alto y Medio Saldaña y río Ataco para estimar proyectos de captación de agua
para su potabilización en el municipio de Ataco-Tolima.
1. Objetivos
1.1 Objetivo General
Realizar un estudio hidrológico de la cuenca del río Saldaña evaluando la posibilidad de un
proyecto de captación de agua para su potabilización y distribución como estrategia para
suplir las necesidades en las veredas aledañas y cabecera municipal del Ataco principalmente
en épocas de sequía.
1.2 Objetivos Específicos
• Caracterizar la cuenca del rio Saldaña identificando parámetros morfológicos, físicos,
geográficos e hidrológicos y climáticos que permiten evaluar la disponibilidad del
recurso hídrico.
• Calcular la oferta hídrica de la cuenca delimitada disponible mediante análisis de
precipitación y evaporación, para proyectos de captación y potabilización para las
comunidades del municipio de ataco; estableciendo caudales máximos para las
temporadas de sequía, en las que se afecta el suministro del líquido a la población.
• Estimar la demanda de agua potable por parte de la población del municipio de Ataco
y su proyección a futuro garantizando la sostenibilidad del recurso.
2. Delimitación de la cuenca
2.1 Delimitación Digital
La delimitación digital se realizo mediante el software ArcGIS Pro haciendo uso del
modelo de elevación digital (DEM) y una capa de recursos hídricos en formato SHP, el DEM
se obtuvo del satélite ASTER por la página EARTHDATA de la NASA, mientras que la capa
de ríos se obtuvo de HYDROSHEDS.
Inicialmente, se procesó el DEM para poder trabajarlo de manera adecuada
reacondicionándolo y rellenando espacios sin información con la herramienta “Fill” para asi,
proceder a utilizar las herramientas de dirección de flujo y acumulación de flujo necesarias
para ver hacia donde se dirige el flujo de la cuenca. Con este proceso aplicado, se realizó un
análisis del DEM que permitió obtener los drenajes de la cuenca necesarios para los cálculos
morfométricos al igual que encontrar el drenaje principal. La delimitación de la cuenca se
puede encontrar en el Anexo XX al final del documento.
2.2 Delimitación manual:
Para realizar esta delimitación de la cuenca seleccionada, usamos un mapa base con
elementos cartográficos que permitieran hacer el cálculo de los parámetros solicitados área
y perímetro, un punto de drenaje con datos conocidos que nos permitió empezar a delimitar
la cuenca basándonos en las líneas divisoras de agua, esta delimitación empezó en el punto
de drenaje y se fue mirando en sentido antihorario el camino optimo hasta volver a cerrar la
figura obtenida, una vez establecida la figura por medio de elementos geométricamente
conocidos y con dimensiones medibles por la escala en la que se imprimió el mapa se calculó
el área de la cuenca, y usando el método de la cuerda se calculó el perímetro, los resultados
de este trabajo se entregan físicamente.
 3. Análisis Morfométrico
Para la realización del análisis morfométrico, se utilizaron tablas y medidas obtenidas en
el proceso mencionado de la delimitación digital obteniendo longitudes, coeficientes, entre
otros valores que permiten analizar la forma de la cuenca, señalados en la siguiente tabla 1.

Tabla 1. Parámetros Morfométricos de la cuenca

Referente al área y perímetro, se puede afirmar que es una cuenca debido a que su área es
mayor a 300 𝑘𝑚2 según Trujillo, 2010. La longitud axial representa la distancia que tiene el
cauce principal desde su inicio en la cuenca hasta el punto más lejano de la misma, indicando
por su gran longitud que pueden tener redes más simples que se distribuyen a lo largo de la
cuenca. El ancho promedio, si se compara con la longitud axial es considerablemente menor;
lo cual representa que la cuenca es de forma alargada. Si se estudia el ancho máximo con el
ancho promedio, hay una diferencia de igual manera considerable que en este caso indica una
cuenca con forma irregular haciendo referencia de que hay zonas mas estrechas u otras zonas
mas amplias.
Para el estudio de la forma y la figura de la cuenca, se hace uso de los coeficientes de
Gravelius, Miller, Elongación y el factor de forma de Horton. Según los valores obtenidos,
la cuenca es de forma oval – oblonga a rectangular -oblonga, es alargada, moderadamente
achatada.
Al analizar la red de drenaje de la cuenca se puede afirmar que la red drenaje no es muy
ramificada debido a que el valor es menor a 3 indicando un flujo menos disperso, presenta
una densidad de drenaje densa al tener valores bajos, la escorrentía se distribuye en el cause
más extenso, por lo cual no hay muchos retardos en la cuenca. Respecto al tiempo de
concentración, se puede decir que no es un tiempo muy elevado para el área de la cuenca
debido a que una gota de agua tardaría aproximadamente 166 minutos en llegar hasta la zona
más alejada del cauce principal.
Respecto a el relieve de la cuenca, se obtuvo una elevación media de 1140.92 metros sobre
el nivel del mar teniendo máximos de 3032 metros y mínimos de 503 metros, mientras que
se observa una pendiente del 21.86% indicando que la cuenca se encuentra en una pendiente
moderada favoreciendo a un escurrimiento medianamente rápido y favoreciendo procesos de
erosión. Si se observa la curva hipsométrica de la figura 1, se observa que tiene una forma
equilibrada entre las zonas altas y bajas donde hay erosión en las partes altas y puede generar
un riesgo moderado de inundaciones. Las zonas altas pueden retener agua mientras que las
zonas bajas favorecen la captación hacia los causes principales.

Elevación (msnm)

% Área Acumulada

Figura 1. Curva Hipsométrica. Fuente: Autor

4. Precipitación
Se presenta el análisis de 30 años de precipitaciones con datos mensuales del municipio
de Ataco, Tolima. Se uso valores mensuales en vez de diarios, debido a la disponibilidad de
la información extraída del IDEAM. Además, se descargó datos de NASA POWER, para
tener una zona completa con información. Aunque el uso de datos mensuales puede limitar
la precisión de los resultados esperados, puede brindar suficiente información para realizar
el proyecto.
4.1. Tratamiento de datos
Después de descargar los datos de 30 años de precipitación en (mm/mes) 10 estaciones
repartidas en el área seleccionada, se trató la información para mejorar su calidad. Se uso el
código propuesto en R por García-Echeverri (2023), en donde se identifica los datos atípicos
con el método de Boxplot, se los elimina y se completan los datos mediante el método
racional deductivo.
El primer paso de este método es hacer una caracterización, en donde se calcula el
porcentaje de datos faltantes, y se crea un espacio en la matriz que pueda contener los datos
vacíos. Posteriormente, se identifican los datos atípicos por medio del diagrama de cajas y
bigotes, en donde los datos se toman a partir de los meses del respectivo año. Al finalizar se
obtiene un diagrama de bigotes de los 12 meses del año con todos los datos, en donde se
identifican los datos atípicos. Al finalizar se cuenta el número de datos atípico y se los elimina
de las series mensuales. Para el llenado de datos faltantes, se identifica cual es el porcentaje
de datos a rellenar, y se usa el método racional deductivo. Este método consiste en calcular
la precipitación media anual de los años sin los datos faltantes, y se calcula los porcentajes
de su respectivo registro. Después, se procede a calcular el valor medio mensual de los
porcentajes y se procede a rellenar los datos faltantes usando la siguiente ecuación:
Σ(𝑃)
𝑃𝑖 = ⋅𝑆
1200 − Σ(𝑆𝑖 ) 𝑖
Ecuación 1: Cálculo de valor medio mensual.
 Donde Si corresponde al porcentaje promedio de los meses desconocidos, Pi a la
precipitación mensual desconocida en cada año incompleto y P a las precipitaciones
conocidas en años incompletos.

4.2. Especialización de precipitaciones

Para estimar la precipitación media mensual en la cuenca en Ataco, Tolima, se realizaron
interpolaciones mensuales IDW con los datos de 10 estaciones para el año de 1990,
considerando los cambios de altura de la zona, que varían entre 500 m y 3000 m. Además,
se decide extraer los valores de las precipitaciones
entregadas por los IDW en las 5 principales áreas de
drenaje, para obtener la precipitación del año de
1990 de la cuenca, como se muestra en la ilustración
1.
Ilustración 1. Ejemplo de IDW y subcuencas.

4.3. Avance de resultados

En las figuras 2 y 3, sugiera una distribución
trimodal de precipitación que varía 212,73 mm/mes
en los periodos húmedos (Febrero, Abril, Mayo,) y
124,24 mm/mes en los periodos secos (Mayo, Junio, Julio-Agosto-Septiembre, Noviembre).

Figura 2. Promedios mensuales de P y ETP.

Figura 3. Coeficiente pluviométrico (1990).

Por último, la precipitación media anual obtenida es de 2253,10 mm/año. Este valor se
acerca al valor reportado del IDEAM en el periodo de 1981-2010 (ver valor en el anexo), el
cual tiene un rango 2000-2500 mm/año.
5. Evapotranspiración
Para este caso de estudio, se optó por calcular la evapotranspiración potencial en mm/mes,
puesto que representa un escenario más crítico para determinar el balance hídrico de la
cuenca, determinando la cantidad de agua que podría ser evaporada y transpirada desde la
superficie incluyendo la vegetación, estimando valores por encima de la evapotranspiración
real en concordancia a los objetivos del proyecto. Con los datos obtenidos mensualmente se
obtuvieron valores anuales y promedios diarios.
5.1 Métodos para calcular la evapotranspiración
Para el cálculo de la evapotranspiración potencial en mm/mes para un periodo de 1990-
2020 con base a la disponibilidad de datos hidrometeorológicos se empleo el método de Turc
Modificado (1962) en el que interviene la radiación global incidente (Rg) en Cal/día/cm2, la
temperatura media mensual en °C y una constante (K) que depende del mes considerando
(0.4 para los meses de 30 y 31 días y 0.37 para febrero). La ecuación esta expresada de la
siguiente manera dado que se tienen humedades relativas superiores al 50 %.
𝑇
𝐸𝑇𝑃 = 𝑘 ( ) (𝑅𝑔 + 50)
𝑇 + 15
Ecuación 2: Ecuación método de Turc Modificado para humedades relativas >50%

Los datos de temperatura media mensual (T) se tomaron de la estación “Mesa de Pole” y
Los valores de la radiación global se tomaron de la página de NASA-POWER para la latitud
y longitud de la “Estación mesa de Pole”. El resultado general indica una evapotranspiración
promedio mensual de 114.26 mm, lo que indica una evapotranspiración diaria de 3.75 mm.
Otro método para analizar la evapotranspiración fue el método de Hargreaves el cual
calcula valores promedios diarios, donde involucra parámetros de temperatura media
mensual, mínima y máxima; y radiación solar extraterrestre en mm/día (𝑅𝑔 ). La ecuación se
expresa de la siguiente manera.
𝐸𝑇𝑃 = 0.0023 × (𝑇𝑚𝑒𝑑𝑖𝑎 + 17.78) × 𝑅𝑔 × (𝑇𝑚á𝑥 − 𝑇𝑚í𝑛 )0.5
Ecuación 3: Ecuación método de Hargreaves simplificada

Sin embargo, para este estudio se usó la ecuación de Hargreaves ajustada por el IDEAM,
la cual ajusta los valores de los coeficientes al contexto colombiano. Esta ecuación esta
expresada de la siguiente manera.
𝐸𝑇𝑃 = 0.00216 × (𝑇𝑚𝑒𝑑𝑖𝑎 + 17.78) × 𝑅𝑔 × (𝑇𝑚á𝑥 − 𝑇𝑚í𝑛 )0.47
Ecuación 4: Ecuación método de Hargreaves ajustado por el IDEAM

Los datos de temperatura media mensual, mínima y máxima se tomaron de la estación

“Mesa de Pole” y los valores de radiación solar extraterrestre para los 12 meses se obtuvieron
a partir de la interpolación en la tabla N°12 de Castro 1985 para la latitud de la estación
“Mesa de Pole”. El resultado obtenido indica una evapotranspiración promedio mensual de
130.43 mm, lo que indica una evapotranspiración diaria de 4.28 mm. Sin embargo, se tiene
que este método puede presentar sobrestimaciones en los resultados por lo cual se recomienda
comparar con datos oficiales. Frente a la variación presentada un tercer método para el
cálculo de la evapotranspiración potencial promedio diaria fue a través de la ecuación de
cenicafé (válida para la zona colombiana). Donde ℎ representa la elevación media de la
cuenca analizada la cual fue calculada en los parámetros morfométricos y correspondieron a
un valor de 1140 msnm.
𝐸𝑇𝑃 = 4.658𝑒 −0.0002ℎ
Ecuación 3: Ecuación Cenicafé
 El resultado obtenido fue de 3.71 mm/día, lo que representa un valor promedio mensual
de 113.06 mm; resultado muy similar al de Turc modificado.
Finalmente, para justificar los resultados obtenidos estos se comparan con datos oficiales
del IDEAM con el mapa de distribución de evapotranspiración potencial, que muestra el
promedio multianual de ETP por una paleta de colores; analizando la zona correspondiente
a la cuenca delimitada se encontró que esta se posicionada sobre la franja de los 1200-1400
mm/anual (Ver anexo 8.4); lo cual ratifica el resultado obtenido por el método de Turc
modificado donde la ETP promedio anual fue de 1370 mm/anual.
6. Conclusiones Preliminares

• Los parámetros morfométricos calculados de la cuenca de estudio nos muestran un

área de forma irregular y también alargada, por lo tanto, esperamos obtener un
hidrograma uniforme con un pico de caudal claramente identificado, al analizar la red
de drenaje vemos que no hay muchos retardos a lo largo de ella entonces se espera
un escurrimiento rápido que genera más susceptibilidad a procesos erosivos en el
terreno.
• Frente a la precipitación se puede decir que tiene 3 picos, y su variación mensual está
en un rango de 212,73 mm/mes y 124,24 mm/mes. En cuanto a la evapotranspiración
potencial, se tiene un comportamiento más estable sin variaciones significativas
oscilando entre los 3.6 y 3.9 mm/día. A manera general se puede apreciar que el
balance hídrico presenta valores superiores de precipitación frente a la
evapotranspiración, sobre todo en los meses de abril y Octubre; exceptuando los
meses de junio, julio, agosto y septiembre correspondientes a temporadas de sequía.
• El trabajo realizado permitió entender la dificultad a la hora de encontrar datos
hidrometeorológicos para las zonas de estudio y la calidad de estos, por lo cual se
necesario el ajuste para su debido procesamiento.
• Para efectuar el balance hídrico de la cuenca se es necesario calcular otros factores
como caudales e infiltración, los cuales son indispensables para en estudio
hidrológico.
 7. Referencias

Alcaldía municipal de Ataco. (2002). ESTUDIO CLIMÁTICO, HIDROLÓGICO,

CARACTERIZACIÓN DE AGUAS SUPERFICIALES Y EVALUACIÓN DE LA
CALIDAD FÍSICO- QUÍMICA DEL AGUA. Ibagué.
Chitiva, L. P. (2013). ESTUDIO DE IMPACTO AMBIENTAL PARA EL PROYECTO DE
EXPLOTACIÓN DE ORO DE ALUVIÓN DEL CONTRATO DE CONCESIÓN
No.HKU-08051 . Ingeominas.
Chow, V. T. (1988). Applied hydrology. McGraw-Hill.
COLCIENCIAS-Unidad de Planeación Minero Energéticas (UPME). (1999). Balances
Hídricos en Colombia-Informe Final. Medellín.
ddd. (dddd de ddddd de ddd). Obtenido de ddd: ddd
ECOS INTERACTIVA. (22 de Noviembre de 2024). Ecos del COMBEIMA . Obtenido de
https://www.ecosdelcombeima.com/region/nota-138481-pobladores-de-ataco-
continuan-consumiendo-agua-no-tratada
Gupta, H. V. (2008). Reconciling theory with observations: Elements of a diagnostic
approach to model evaluation. Hydrological Processes.
Instituto de Hidrología, Meteorología y Estudios Ambientales (IDEAM). (s.f.). IDEAM.
Obtenido de http://www.ideam.gov.co/
Ministerio de Vivienda Ciudad y Territorio. (8 de Junio de 2017). Resolución 330 RAS. 182.
N., L. G. (2012). Metodología para el cálculo de la humedad del suelo usando parámetros
topográficos (MDE), climáticos y edáficos en un sector del piedemonte
depositacional del municipio de Villavicencio. Bogotá, Colombia.
Trujillo, A. U. (2010). Guía básica para la caracterización morfométrica de cuencas
hidrográficas (1st ed.). Universidad del Valle.
García-Echeverri, C. (2023). Identificación de atípicos y llenado de faltantes. RPubs.
https://rpubs.com/cagarciae/atipicos_faltantes

Wajarai Chuim, J. L., & Jaramillo Veliz, J. J. (2023). Evaluación del balance hídrico
superficial de la subcuenca hidrográfica del río Jipijapa. Revista Científica Arbitrada
Multidisciplinaria PENTACIENCIAS, 5(3).
https://doi.org/10.59169/pentaciencias.v5i3.493
 8. Anexos
8.1 Resultado Delimitación Digital

8.2 Resultados de procesamiento de datos de precipitación por estación.

CONDER
EL
(22020070)

Precipitacion mm/mes de 1990-2020

800

600

400

200

0
SANTIAG
O PEREZ
(22020050)

Precipitacion mm/mes de 1990-2020

800

600

400

200

PAN DE
AZUCAR
(22050080)

Precipitacion mm/mes de 1990-2020

800

600

400

200

0
SAN
PEDRO
(22050070)

Precipitacion mm/mes de 1990-2020

1000

800

600

400

200

MESA DE
POLE
(22055020)

Precipitacion mm/mes de 1990-2020

1000

800

600

400

200

0
CASA DE
ZINC
(22020040)

Precipitacion mm/mes de 1990-2020

800

600

400

200

GAITAN
(22010070)

Precipitacion mm/mes de 1990-2020

800

600

400

200

0
ATACO
(22050040)

Precipitacion mm/mes de 1990-2020

800

600

400

200

HERRER
A
(22010010)

Precipitacion mm/mes de 1990-2020

400

300

200

100

0
PUNTO
DE NASA
POWER

Precipitacion mm/mes de 1990-2020

400

300

200

100

0
8.3 Atlas Precipitación media anual 1981-2010
8.4 Mapa Distribución Evaporación Potencial Total