Pontificia Universidad Javeriana

Ingeniería civil

Proyecto Final

Asignatura:
Hidrología

Estudiantes:
Alejandro Peñaranda
Laura Catalina Balaguera
Laura Catalina Gutiérrez
Luis Carlos Rocha

Profesora:
Rosa Angelica Rodríguez Fernández

Bogotá D.C
2022

1
Contenido
Introducción: ..................................................................................................................................... 3
Caracterización de cuencas hidrográficas: ..................................................................................... 3
Puente Quebrada Negra................................................................................................................... 3
• Caracterización de la cuenca hidrográfica........................................................................... 3
• Caracterización de la red de drenaje.................................................................................... 6
• Relieve de la cuenca .......................................................................................................... 10
• Análisis de la amenaza ante crecientes súbitas en los puentes vehiculares¡Error! Marcador
no definido.
Puente Únete ................................................................................................................................. 17
• Caracterización de la cuenca hidrográfica......................................................................... 27
• Caracterización de la red de drenaje.................................................................................. 28
• Relieve de la cuenca .......................................................................................................... 30
Puente Blanco................................................................................................................................ 31
• Caracterización de la cuenca hidrográfica......................................................................... 31
• Caracterización de la red de drenaje.................................................................................. 33
• Relieve de la cuenca .......................................................................................................... 34
• Análisis de la amenaza ante crecientes súbitas en los puentes vehiculares ....................... 35
Análisis de precipitación:................................................................................................................ 35
Información de precipitación media mensual Anual ..................................................................... 35
Precipitación media / Ciclo Interanual .......................................................................................... 36
Precipitación promedio multianual anual y mensual al interior de la cuenca ............................... 37
Evaporación:.................................................................................................................................... 38
Temperatura promedio multianual anual en cada estación ........................................................... 38
Evaporación potencial ................................................................................................................... 38
Índice de aridez y Evaporación Real ............................................................................................. 39
Balance Hídrico: .............................................................................................................................. 42
Balance hidrológico de largo plazo en la cuenca .......................................................................... 42
Conclusión: ...................................................................................................................................... 43
Bibliografía: ..................................................................................................................................... 43

2
Introducción:
Durante el desarrollo del proyecto el cual consistió en la caracterización de tres cuencas
hidrográficas haciendo en una primera instancia una caracterización de los parámetros
geométricos de la cuenca, como su área o su perímetro llegando a hacer la caracterización de
la red de drenaje que compone esta cuenca, su respectivo relieve. Por último, para esta
sección realizado una revisión de posible amenaza que este pudiera presentar a estructuras
como puentes que están en el cause principal por un evento de creciente súbita.
Para un segundo momento se hace un análisis de la precipitación, evaporación y Balance
hídrico de la cuenca de Puente blanco con lo cual se obtienen datos de esta cuenca para
múltiples periodos de tiempo y se llega a una mejor comparación de cómo se mueve el agua
dentro de esta teniendo en cuenta lo que se precipita y lo que se evapora por causa de las
temperaturas que se presentan en este sector.

Caracterización de cuencas hidrográficas:

Se comenzó la caracterización de las tres cuencas usando los modelos digitales de terreno
con las respectivas elevaciones además se realizó mediante el software geográfico
GRASSGIS haciendo exactamente el mismo proceso para las tres cuencas obteniendo la red
de drenaje los parámetros geométricos de las tres cuencas.
Para esto se procedió a generar la región en la cual se procede a trabajar, realizando la re-
projeccion del mapa del sistema coordenado de WGS84 a Origen Nacional para lo cual se
realizó la conversión de las coordenadas generando la nueva región de trabajo y realizando
la re-projeccion de tanto el modelo de elevación y el punto de concentración de cada una de
las cuencas.
Para posteriormente proceder a la cauterización de cada cuenca generando un modelo de
direcciones de flujo y con esto lograr delimitar la cuenca para poder obtener su área, longitud,
pendiente, ordenamiento de Horton, Strhler, o la curva hipsométrica de cada una de las
cuencas.
Puente Quebrada Negra
• Caracterización de la cuenca hidrográfica
Área y perímetro de la cuenca

Se encuentran usando el mapa que contiene la línea vectorial del límite de la cuenca, y
con el comando de [Link], se obtiene un reporte del área y del perímetro en kilómetros
cuadrados y kilómetros, respectivamente.
𝐴 = 3.77626735731336𝑘𝑚2
𝑃 = 11.5885974128079 𝑘𝑚
Longitud y ancho de la cuenca
Para medir la longitud de la cuenca se utilizan las herramientas de análisis del mapa, para
medir la distancia longitudinal. Por otro lado, se encuentra el ancho (W) de la cuenca, por

3
medio de la relación entre área y longitud. El ancho de la cuenca es como ese ancho promedio
de la cuenca se relaciona con el factor de forma y la longitud la distancia desde el punto de
concentración hasta el punto más lejano de la cuenca
𝐿𝑜𝑛𝑔𝑖𝑡𝑢𝑑 𝑑𝑒 𝑙𝑎 𝑐𝑢𝑒𝑛𝑐𝑎 = 2.588 𝑘𝑚
𝐴 3.77626735731336𝑘𝑚2
𝑊= = = 1.45915 𝑘𝑚
𝐿 2.588 𝑘𝑚
Parámetros de forma
En cuanto a los parámetros de forma de la cuenca, se calculó el índice de compacidad de
Gravelius (k) y el factor de forma (kf). El primero, es la relación entre el perímetro de la
cuenca y el perímetro de un círculo con la misma área, nos dice que tan compacta o alargadas
es la cuenca dependiendo del valor de k. Si k>1.75, esta es alargada, si k˂1.25, la cuenca es
compacta.
El valor se obtiene mediante:
𝑃 11.5885974128079𝑘𝑚
𝐾 = 0.282 = 0.282 = 1.6817
√𝐴 √3.77626735731336 𝑘𝑚2

𝐴 𝑊 1.45915𝑘𝑚
𝐾𝑓 = = = = 0.563814
𝐿2 𝐿 2.588𝑘𝑚
Debido a que el índice de Gravelius indica que 1.25<K<1.75 se determina de que a pesar de
que la cuenca no se encuentra en uno de estos valores, si tiende por cercanía a 1.75 a ser
alargada, sin embargo, su forma sigue siendo muy equilibrada entre su ancho y su largo. Esta
tiene cauces más cortos respecto al cauce principal, lo cual se traduce en que las
precipitaciones se amortiguan por el río, por cual el río crece paulatinamente. Sus tiempos de
concentración son altos debido a su forma y terreno pues al no poseer zonas de elevadas
pendientes el agua tiende a estancarse más y, por lo tanto, sus crecidas son atenuadas. Para
el segundo parámetro, se hace la relación entre área y la longitud al cuadrado de la cuenca;
si este es mayor a 0.75, la cuenca tiende a ser cuadrada, y, si es menor a 0.25, la cuenca tiende
a ser alargada, esto debido a que el factor de forma compara la geometría de la cuenca con
un cuadrado. Igualmente, el factor de forma muestra que Kf=0.563814 es mayor a 0.25 y más
cercano al 0.75 por lo cual la cuenca tiende a ser cuadrada.

Kc Forma de la cuenca
1.00-1.25 Redonda a oval oblonga
1.25-1.50 De oval redonda a oval oblonga
1.50-1.75 De oval oblonga a rectangular oblonga
>1.75 Rectangular
Tabla 1: Índice de Gravelius

4
Su forma según la tabla anterior se puede definir como, una cuenca de oval oblonga a
rectangular oblonga.
Longitud cauce principal
Por otra parte, para lograr obtener el valor de la longitud del cauce principal, se creó una
máscara de recorte para que Grass solo incluya la información dentro del límite de la cuenca.
Como se observa en la ilustración, y a partir del comando [Link], se recortó el mapa raster
de colores, y luego, con el comando de [Link], se recortaron las corrientes a partir del límite
de cuenca.

Ilustración 1 Cuenca de Análisis

Este cauce principal se puede interpretar como la longitud de la corriente perteneciente al

máximo orden ya que es como el “soporte-tronco”, si se compara con un árbol, de toda la
cuenca. Con este mapa intervenido, se obtuvo la longitud del cauce principal, que para la
cuenca de Puente Quebrada Negra dio 2.9480 km

Ilustración 2 Longitud del Cauce principal SIG

5
 • Caracterización de la red de drenaje
Dentro de la cuenca, el agua que cae en esta se escurre mediante un sistema de corrientes, o
una red de drenaje que lleva toda el agua hacia un mismo punto de salida (ArcGIS, 2019).
Este sistema tiene forma de árbol, debido a que sus corrientes se van ramificando a partir de
una corriente principal. Para cuantificar esta red, Horton definió relaciones para estas
corrientes, también conocidas como índices de Horton.
Orden de la cuenca
A las corrientes se les da un orden acumulativo según su tamaño, siendo el cauce principal
el mayor. Las corrientes comienzan como orden 1, y, si se juntan dos corrientes del mismo
orden, estas forman una corriente de orden mayor. En cambio, si se juntan dos corrientes de
distinto orden, la corriente continuará con el orden mayor. Esto permite el análisis fractal de
la red.

Ilustración 3 Ordenamiento de una cuenca orden 4. (Salgar, 2021)

Las corrientes de menor orden, consideradas como corrientes jóvenes, se relacionan con
mayores pendientes en la parte alta de la cuenca, y, por lo tanto, las corrientes de mayor orden
se consideran con menores pendientes debido a que se encuentran en las zonas bajas. Esta
jerarquización influye en gran medida en procesos de erosión, transporte de sedimentos y
deposición. Estos procesos deben analizarse en conjunto con la topografía y las poblaciones
cercanas para evitar posibles desastres y desarrollar sistemas de evacuación o de alerta
temprana, ya que, dependiendo de la bifurcación de las corrientes, la red tendrá una respuesta
bien sea lenta o rápida ante un pulso de precipitación. Para obtener esta información de la

6
cuenca analizada, y a partir del comando con el comando [Link], se obtuvieron estas
estadísticas de la red de drenaje de Strahler. Se utilizaron los mapas de la red de drenaje de
Strahler, el mapa raster de direcciones y el mapa de elevaciones para obtener las estadísticas
de las corrientes y de la cuenca. El orden mayor obtenido es de 5.

Ilustración 4 Orden de la cuenca SIG

Densidad de drenaje
Es la longitud de los cauces dentro de la cuenca a razón del área total de drenaje (A), es uno
de los parámetros más importantes de la morfometría de una cuenca dando una idea de la
erosión de esta, relaciona la suma de las longitudes de las corrientes con el área de la cuenca.
Como se observa, la longitud total que se tiene de las corrientes es de 33.1139 km y con el
área obtenida anteriormente, se calcula la densidad:
𝑙𝑜𝑛𝑔 𝑡𝑜𝑡𝑎𝑙 𝑐𝑜𝑟𝑟𝑖𝑒𝑛𝑡𝑒𝑠 33.1139 𝑘𝑚
𝐷𝑒𝑛𝑠𝑖𝑑𝑎𝑑 𝑑𝑒 𝑑𝑟𝑒𝑛𝑎𝑗𝑒 = = = 8.7689𝑘𝑚
á𝑟𝑒𝑎 3.77626735731336 𝑘𝑚2

Este valor representa cada cuántos kilómetros cuadrados hay un kilómetro de corriente. La
densidad es baja, lo que significa que cada 8.7689 kilómetros cuadrados hay una corriente.
Los valores de densidades u ordenes de corrientes pequeñas se observan donde los suelos son
muy resistentes a la erosión o muy permeables, es un área pobremente drenada, con respuesta
hidrológica lenta. Por la manera en que la cuenca se desarrolla en su interior y los parámetros
encontrados se puede decir que la cuenca posee un sistema de drenaje dendrítico que se
identifica por las siguientes características:
• Sistema de canales que asemejan la ramificación de un árbol.
• Se da sobre pendientes suaves donde los materiales responden homogéneamente a la
erosión.
• Las zonas planas con rocas sedimentarias presentan un patrón dendrítico
Relación de bifurcación
Asimismo, se calcularon los índices de Horton. El primero, la relación de bifurcación (Rb), la cual
explica que tantas veces se divide, o se ramifica, una corriente. Para esto, se utilizaron los datos de
orden de corrientes y el número de corrientes por cada uno. A partir de esta información se hizo el
respectivo ajuste de mínimos cuadrados para sacar el valor de Rb en Excel. Para ello, se graficó el

7
orden i de las corrientes versus el logaritmo del número de corrientes por orden. A estos valores se
les sacó una línea de tendencia con su respectiva ecuación y, utilizando el valor de la pendiente m, se
obtuvo Rb:

Cálculo de Rb (ajuste)
i Ni Ln(Ni)
1 91 4,51085951
2 15 2,7080502
3 5 1,60943791
4 1 0
5 1 0
m= -ln(Rb) -1,173
Rb 3.23

Tabla 2 Cálculo de relación de bifurcación

𝑦 = 𝑚𝑥 + 𝑏 → 𝑦 = −1,173 + 5,2846; 𝑚 = −1,173 = ln(𝑅𝑏) ; 𝑅𝑏 = 𝑒 𝑚 = 𝑒 1,173 = 3.23

BIFURCACIÓN
5
y = -1,173x + 5,2846
4 R² = 0,937
3

0
1 2 3 4 5
-1

Ilustración 5 Relación de Bifurcación

Este valor relaciona el número de canales i-1 con el de i, para identificar la correlación entre
cuantas veces se dividen las corrientes. El rango normal de este índice está entre 3 a 5, por lo
que el valor se encuentra dentro de este rango (Cardona, s.f)

Relación de longitudes
Este índice relaciona las longitudes promedio por cada orden i, es decir, que tanto desarrollo
tienen las corrientes de orden uno en comparación con los demás órdenes. Como la cuenca,
según los factores de forma y lo visto en los planos del SIG, su forma tiende a ser alargada
ya que se aleja a la semejanza de un círculo y se acerca al de una elipse. Esto se relaciona
con el desarrollo del cauce principal, el cual debería ser mayor a los cauces de orden menor.

8
 Cálculo de Rl (ajuste)
i Ni Ln(Ni)
1 11,5911 2,45023756
2 8,1592 2,09914612
3 7,7253 2,04450066
4 2,6903 0,98965271
5 2,948 1,08112697
m= -ln(Rl) -0,3848
Rl 0,680587

Tabla 2Cálculo de relación de longitudes

Para calcular este índice se utilizaron los órdenes i, y las longitudes promedio obtenidas del
SIG de todos los cauces de cada orden. Al obtener la gráfica de los órdenes versus el
logaritmo de las longitudes, se hizo el respectivo ajuste y se obtuvo la ecuación:
y = mx + b ￫ y = −0.3848x + 2.8872; m = −0.3848 = −ln(RL); RL = 𝑒 𝑚 = 𝑒 0.3848
= 0.680587

Ilustración 6 Relación de longitudes

Relación de áreas
Este índice relaciona las áreas promedio por cada orden i, es decir, que tanto desarrollo tienen
las corrientes de orden uno en comparación con los demás órdenes. Como la cuenca, según
los factores de forma y lo visto en los planos del SIG, su forma tiende a ser alargada ya que
se aleja a la semejanza de un círculo y se acerca al de una elipse. Esto se relaciona con el
desarrollo del cauce principal, el cual debería ser mayor a los cauces de orden menor.
TIEMPO DE CONCENTRACIÓN
Es el tiempo mínimo necesario para que todos los puntos de la superficie de la cuenca
contribuyan simultáneamente al caudal recibido en la salida. Para encontrar este tiempo,

9
distintos autores encontraron distintas fórmulas para el cálculo de este tiempo de
concentración (tc), para distintas regiones en el mundo. Por esta razón, muchos de los
cálculos de este tiempo varían su valor y unidades, ya que cada ecuación fue diseñada para
lugares distintos. En este caso se escogieron dos ecuaciones principales:
Ecuación de Kirpich: Ecuación de Temez:
𝐿 𝐿
𝑇𝑐 = 0.066 ∗ ( )0.77 𝑇𝑐 = 0.3 ∗ ( 0.25 )0.77
√𝑆 𝑆
La L es la longitud del canal desde aguas L como la longitud del cauce principal en
arriba hasta la salida, en unidades de pies. Y kilómetros, el área de la cuenca y S es la
S es la pendiente en pies/pie calculado a pendiente del cauce principal, es decir, de
partir del SIG. Esta ecuación se obtuvo los datos de la corriente de mayor orden de
convirtiendo la longitud de la cuenca de km la cuenca de puente quebrada negra.
a pies, y con la pendiente, pasar este tiempo
a horas.
Tc(min)=74.1032min=1.23horas Tc(min)=207.138min=3.4523horas

Extensión media de escorrentía superficial

Este parámetro explica que tanta área disponible que tiene una cuenca para la escorrentía.
Esto depende de la pendiente, ya que, a partir de esta, la capacidad que tiene una cuenca de
crear escorrentía superficial varía. Para esta cuenca, el parámetro dio de 0.02851Este valor
explica la distancia media que el agua de la precipitación tendrá que transportarse hasta un
cauce de agua cercano.

Se calcula a partir de la siguiente ecuación:

𝐴 3.77626735731336𝑘𝑚2
𝐸𝑠 = = = 0.02851 𝑘𝑚
4𝐿𝑖 4 ∗ 33.1139 𝑘𝑚
Es: Extensión de escurrimiento superficial
A= Superficie de la cuenca en 𝑘𝑚2
Li=Largo total de los cursos de agua en Km
• Relieve de la cuenca
Pendiente de la cuenca
La pendiente de la cuenca es de suma importancia, ya que, según la pendiente de una cuenca,
el tiempo de concentración varía, además de la rapidez de evacuación y escorrentía del agua.
Para calcularla se utilizaron 2 métodos. El primero fue por medio del SIG, a partir del
comando [Link], y utilizando el mapa ráster de elevación (MDT), se generó un mapa
10
de salida de pendiente ráster en unidad de porcentaje. Todo esto se realiza siempre con la
máscara activa, solo interesa ver lo que pasa dentro de la cuenca no fuera de ella, para una
mejor interpretación se adjunta el mapa de pendientes únicamente de la cuenca

Ilustración 7Mapa de pendientes SIG

Como se observa en el mapa, las pendientes altas son las de colores rojo, negro, mientras que
las pendientes bajas de colores verdes. Es posible observar las partes planas de la cuenca con
facilidad y aún más teniendo en cuenta la Ilustración 1., como algunas partes de pendientes
muy bajas hacia Puente Quebrada negra. Luego de tener este mapa, se utilizó el comando
[Link] para calcular los momentos estadísticos de un mapa, en este caso de pendientes, saca
el promedio, el mínimo, el máximo, y muchas más estadísticas principales de cada pixel.

11
 Ilustración 8 Datos SIG

La imagen anterior muestra los resultados arrojados por el SIG para la pendiente promedio. Este
valor es de 47.7398%, lo que indica que es un terreno fuerte.

El segundo método utilizado para calcular este mismo parámetro fue a partir del criterio de Alvord,
que se caracteriza por medio de la siguiente ecuación:

∑ 𝐿𝑖 ∗ 𝐸𝑞
𝑆𝑐𝑢𝑒𝑛𝑐𝑎 =
𝐴
La pendiente de la cuenca (S) es igual a la sumatoria de las longitudes (Li) de las curvas de
nivel de la cuenca, multiplicado por la equidistancia vertical (Eq) entre curvas, y dividido
sobre el área (A). Para obtener los datos necesarios primero se comenzó utilizando el
comando [Link], deponiendo como mapa de entrada del modelo digital del terreno en
formato ráster, y como salida, un mapa tipo vectorial debido a que las curvas de nivel que se
obtienen son líneas. Para esto se sacaron tipos de curvas con información, una para
incremento de niveles de contorno de 50m.
Sin embargo, para obtener el criterio de Alvord se necesitan las longitudes de las curvas, por
lo que a partir del comando [Link], y, con los planos de curvas ya obtenidos, se utilizaron
para sacar las longitudes en metros de cada uno.

12
 PARA 50
cat level length
5 950 1.35475181
6 1000 2.3027588
7 1050 2.48474264
8 1100 2.24371723
9 1150 3.90520777
10 1200 4.32112413
11 1250 2.21905686
12 1300 2.39116301
13 1350 4.95343009
14 1400 5.23171644
15 1450 2.5431793
16 1500 2.58737335
17 1550 2.63686041
18 1600 5.20887661
19 1650 4.79071298
20 1700 4.37717033
21 1750 4.75259479

Tabla 3 Categorías longitudes curva de nivel 50m

Para las curvas que van cada 50 metros se obtuvieron 17 categorías, cada una con el total de
longitudes por cada una. Para obtener el criterio de Alvord (50m), se hizo la sumatoria de las
longitudes, con un total de 58.3044366km = 58304.4366 m. Cabe aclarar que las medidas de
tomaron en metros, por lo que el área se pasó de km2 a m.
58.3044366 𝑘𝑚 ∗ 50
𝑆 𝑐𝑢𝑒𝑛𝑐𝑎 50 𝑚 = = 77.19%
3.77626735731336𝑘𝑚2 ∗ 1000 ∗ 1000

Para las categorías obtenidas de las curvas de nivel a cada 10m se hizo el mismo
procedimiento. Cabe aclarar que todas las medidas se tomaron en metro por lo que el área se
pasó de km2 a m

13
 PARA 10
cat level length
18 920 0.25068896
19 930 0.44611707
20 940 1.21304559
21 950 1.35475181
22 960 1.56801297
23 970 1.69676866
24 980 1.77653352
25 990 1.86386521
26 1000 2.3027588
27 1010 2.52336688
28 1020 2.58035803
29 1030 2.68564305
30 1040 2.29867053
31 1050 2.48474264
32 1060 2.67424655
33 1070 2.68464387
34 1080 2.85495019
35 1090 2.97598309
36 1100 2.24371723
37 1110 3.31984572
38 1120 3.58152092
39 1130 3.76970753
40 1140 3.85880952
41 1150 3.90520777
42 1160 4.0070271
43 1170 2.09931962
44 1180 2.13289596
45 1190 2.12929384
46 1200 4.32112413
47 1210 4.28216706
48 1220 4.4145916
49 1230 4.52735625
50 1240 2.19708468
51 1250 2.21905686
52 1260 2.26737347
53 1270 4.39179162
54 1280 4.4815767
55 1290 4.56997422
56 1300 2.39116301
57 1310 4.63955766
58 1320 4.66686887
59 1330 4.72701673
60 1340 4.88253727
61 1350 4.95343009
62 1360 5.11292115
63 1370 2.61998735
64 1380 2.64819603
65 1390 2.61789102
66 1400 5.23171644
67 1410 5.33315573
68 1420 5.41459647
69 1430 2.69337291
70 1440 2.63194297
71 1450 2.5431793
72 1460 2.53563188
73 1470 2.61249061
74 1480 2.59173352
75 1490 2.58435598
76 1500 2.58737335
77 1510 2.53114219
78 1520 4.96348612
79 1530 4.96864157
80 1540 2.6773709
81 1550 2.63686041
82 1560 2.57817239
83 1570 5.24740892
84 1580 2.62902569
85 1590 5.26382791
86 1600 5.20887661
87 1610 5.42015039
88 1620 5.3313681
89 1630 5.28523298
90 1640 5.04241979
91 1650 4.79071298
92 1660 4.67355153
93 1670 2.30392062
94 1680 4.57356446
95 1690 2.21292577
96 1700 4.37717033
97 1710 4.53046762
98 1720 4.59715583
99 1730 4.53079403
100 1740 2.31014734
101 1750 4.75259479
102 1760 4.76111191
103 1770 3.47453341
104 1780 3.41097291
105 1790 0.64914375

14
 296.706459 𝑘𝑚 ∗ 10
𝑆 𝑐𝑢𝑒𝑛𝑐𝑎 10𝑚 = = 78.57%
3.77626735731336𝑘𝑚2 ∗ 1000 ∗ 1000
Es evidente que por cada método empleado las pendientes tuvieron una variación muy
amplia, esto en parte se debe debido a la definición a la cual se trabajó cada método. SIG
tiene una mejor precisión puesto que trabaja al pixel por lo cual posiblemente su pendiente
sea muy similar a la real, pero si bien se calculan por el criterio de Alvord se encuentra que
la exactitud de este está dada por la definición de las curvas de nivel. Entre menos elevación
exista entre ellas más certero será el resultado. De acuerdo con las pendientes encontradas y
la siguiente tabla se determina que el terreno es de tipo escarpado en su mayoría.
PENDIENTE MEDIA % TERRENO
0-2 Llano
2-5 Suave
5-10 Moderadamente accidentado
15-25 Muy accidentado
25-20 Escarpado
>50 Muy escarpado

Ilustración 9 Tabla clasificación de terrenos mediante pendientes

Elevación media de la cuenca

Ilustración 10 Datos obtenidos del SIG- estadísticas del MDT para la cuenca

Para obtener esta elevación promedio se utilizó el comando [Link], y, utilizando el modelo
digital del modelo, se sacaron las diferentes estadísticas. Como se observa en la imagen de
los resultados obtenidos e GRASS, la media obtenida fue de 1370.51 metros. La altura
mínima de 917 metros se espera que esté en las zonas más cercanas al punto de concentración
que es el más bajo ya que si hasta ahí desencadena, es porque la pendiente es menor, entonces
esta más bajo que el punto más alto de 1797 metros.

Curva hipsométrica de la cuenca

Esta curva se utiliza para representar la variación de elevación de una cuenca, es un parámetro
de relieve que permite describir la cuenca y se representa mediante un gráfico que relaciona
el área de drenaje que existe por encima o por debajo de un nivel. Además, esta curva también
15
funciona como indicativo para la edad de la cuenca y para hacer una comparación con otras
cuencas; si la curva tiene una derivada negativa significa que el cambio de sus pendientes es
alto, mientras que, si la curva genera una derivada positiva, significa que no tiene mayores
cambios en sus pendientes. Para crear esta curva, se comenzó utilizando el comando de
[Link], y, utilizando el modelo digital del terreno que contiene las elevaciones se obtuvieron
los valores de las elevaciones versus el percentil, lo que indica cuanto porcentaje hay de
alturas en la cuenca. Las mayores alturas se muestran cercanas al 0, mientras que las menores
se acercan a 1.

CURVA HIPSOMÉTRICA CURVA HIPSOGRAFICA

Ilustración 11 Curva hipsométrica e hipsografica

Pendiente del cauce principal

La pendiente, en primera instancia, se calculó en el SIG como se muestra en la imagen. Este
valor indica que en la parte más baja de la cuenca, en donde ya se puede considerar un terreno
plano pues es el punto más bajo de la cuenca y hasta ahí llega el agua que fluye por toda la
red de drenaje. El agua en este punto se dirige ya hacia un nuevo proceso. La pendiente del
cauce principal es de 0.299%.

Ilustración 12 Pendiente del cauce principal SIG

16
Puente Únete

A continuación, se presenta el procedimiento de cómo se realizó la caracterización de la

cuenca hidrográfica de Puente Únete (ubicado en Casanare) realizando la extracción de la
red de drenaje, caracterización geométrica de área, perímetro y como se obtienen los
ordenamientos de esta cuenca. Y el paso a paso de la caracterización de cada uno de los
distintos parámetros de la cuenca con sus respectivos datos mostrando el paso a paso de lo
que se hizo al obtener información de esta cuenca, el cual se muestra a continuación:

Figura 1 (cargar el DEM)

Como se muestra en la figura anterior se puede apreciar la carga del modelo digital de terreno
obtenido para poder utilizarlo en el procedimiento, pero antes de esto fue necesario realizar
la creación del directorio de mapas con la creación de la ubicación de WGS84 para poder
cargar el mapa que se descargó de Hydroshed, e importar el archivo tif correspondiente con
el comando [Link], con el propósito de poder trabajar sobre este tif el área de interés. Cabe
resaltar que es importante definir de forma de adecuada la región de trabajo con las
coordenadas correctas, para que no quede una zona muy grande sin utilizar ni se corte la
cuenca a la hora del análisis. Después se procede a pasar al entorno de trabajo con origen
nacional, lo cual requiere que se transformen las coordenadas de un sistema a otro, para esto
se utilizó la página del servicio geológico colombiano para poder realizar la conversión de
las coordenadas a origen nacional de forma correcta .

17
 Ilustración 13: DEM con región de trabajo en WGS84

Se realizo la creación de la región de trabajo con el comando de [Link] se delimita la región

de trabajo con las coordenadas transformadas para empezar a establecer y delimitar la cuenca
identificando el punto de concentración como un archivo vectorial.

Ilustración 14: Cambio de coordenadas

N 2154718.16
S 2127074.66
E 5060907.97
O 5033235.47

Tabla 4 coordenadas origen nacional

18
 Ilustración 15: Definición de región de trabajo a Origen Nacional

Ilustración 16: Mapa origen nacional

Luego de realizar la re-projeccion del DEM dentro de la región de trabajo y del punto de
concentración usando los comandos [Link] para el mapa raster y [Link] para el punto de la
estación, para así poder comenzar con la delimitación y caracterización de la cuenca.

19
 Ilustración 17: Comando [Link]

Empezando para obtener los arroyos y las direcciones de flujo de la cuenca y de la red de
drenaje, a partir del mapa de elevación (DEM), junto con el flujo mínimo para la acumulación
que hay en la cuenca y obtenemos las corrientes en formato ráster y en formato vectorial,
cabe resaltar que según el flujo mínimo escogido para la acumulación de flujo para arroyos
los mapas pueden verse más o menos saturados, por lo cual hay que escoger este flujo con
atención.

Ilustración 18: Mapas de salida comando [Link]

Al ejecutar el comando se obtienen estos mapas donde se puede ver las direcciones de flujo
de la cuenca y de todo el sector que está dentro de la región de trabajo especificada en origen
nacional.
20
Usando el comando [Link] se obtienen el mapa vectorial de red de flujo del sector
especificado según el flujo mínimo elegido.

Ilustración 19:Mapa de flujo vectorial delimitando la cuenca

Ilustración 20: Máscara cuenca delimitada

Se genera una máscara del espacio de la cuenca ya delimitado permitiendo identificar la red
de drenaje de la cuenca hasta el punto de concentración generado desde el inicio del trabajo
con la cuenca. Luego esto se usa para encontrar datos como el área de la cuenca entre otros
parámetros.
Se hace el comando para poder delimitar el área de la cuenca generando también un mapa
vectorial a partir de la máscara generada anteriormente.

21
 Ilustración 21: Mascara Vectorial

Se realiza la modificación a la máscara para poder identificar los demás parámetros

geométricos de la cuenca dejando ver toda la red de drenaje de la cuenca en cuestión usando
el comando [Link].

Ilustración 22: capas generadas hasta el momento y cuenca delimitada

Este es el mapa que se obtiene después de realizar todo el procedimiento descrito

anteriormente y se delimita la cuenca.

22
 Ilustración 23: comando [Link]

Usando el comando [Link] se obtiene el dato de área de la cuenca definiendo el área de la

zona vectorial definida. Luego de convierte de vectorial a un archivo de líneas para realizar
la medida del perímetro de la cuenca. Esto usando el comando [Link] y luego volviendo a
usar al comando [Link] se medirá la longitud de la línea asi se mide el perímetro de la
cuenca.

23
 Ilustración 24: Corrientes de la cuenca

Usando el comando [Link] se puede obtener las diferentes corrientes del mapa de la cuenca
analizada.
Usando el comando [Link] para encontrar los datos de ordenamiento de las cuencas
se obtiene el ordenamiento de la cuenca tanto por el método de Horton y el método de
Strahler, que para el caso de la cuenca de puente únete es de 5, también se obtiene la pendiente
y longitud del cauce principal las cuales son 0.0617 y 21.7123km respectivamente.

Ilustración 25: Resultado de ordenamiento de la cuenca

24
 Ilustración 26:Resultados [Link]

Ilustración 27:Resultados [Link] 2

Utilizando el comando [Link] se puede sacar la longitud del cause principal el orden
de la cuenca y la pendeinte promedio de esta generando apatrir de datos estadisticos que
genera el programa.
Para generar las pendientes las pendientes y orientaciones para la cuenca que se esté
analizando a partir del mapa ráster y formato en el cual va a salir el valor de la pendiente
del cauce principal.

25
 Ilustración 28:Resultados [Link]

Con el comando [Link] se obtienen las curvas de nivel con las cuales se puede identificar
el relieve de la cuenca el cual se hizo para el caso de puente únete con una separación de 50m
dando una mejor aproximación a el relieve del lugar.

Ilustración 29:Comando [Link]

26
Por último se realiza por medio del comando [Link] el programa realiza el cálculo de la
gráfica hipsométrica de las cuencas a partir de lo obtenido con
• Caracterización de la cuenca hidrográfica
Área y perímetro de la cuenca.:

Ilustración 30: Mapa cuenca Puente Únete

Para la cuenca de Puente únete se obtuvo el área y de perímetro de esta cuenca usando el
mapa obtenido del mapa que contiene la línea vectorial de longitud de la cuenca en
kilómetros.
𝐴 = 138.94502569 𝑘𝑚2
𝑃 = 70.269 𝑘𝑚
Longitud y Ancho de la cuenca:
Usando las herramientas de medición longitudinal del programa es posible realizar la
medición de la longitud de la cuenca y utilizando la relación entre área y longitud de la cuenca
es posible encontrar el acho de esta.
𝐿 = 15.79 𝑘𝑚
𝐴
𝑊= = 8.799558308 𝑘𝑚
𝐿
Parámetros de forma de la cuenca:
Para los los parámetros de forma de la cuenca, se calculó el índice de compacidad de
Gravelius (k) y el factor de forma (kf). El primero, es la relación entre el perímetro de la
cuenca y el perímetro de un círculo con la misma área, nos dice que tan compacta o alargadas
es la cuenca dependiendo del valor de k.
El valor se obtiene mediante:

27
 𝑃 70.269 𝑘𝑚
𝐾 = 0.282 = 0.282 = 1.683820661
√𝐴 √138.495 𝑘𝑚2

𝐴 𝑊 8.799𝑘𝑚
𝐾𝑓 = = = = 0.557251425
𝐿2 𝐿 15.79𝑘𝑚
Debido a que el índice de Gravelius indica que 1.25<K<1.75 se determina de que a pesar de
que la cuenca no se encuentra en uno de estos valores, si tiende por cercanía a 1.75 a ser
alargada, sin embargo, su forma sigue siendo muy equilibrada entre su ancho y su largo. Esta
tiene cauces más cortos respecto al cauce principal, lo cual se traduce en que las
precipitaciones se amortiguan por el río, por cual el río crece paulatinamente. Sus tiempos de
concentración son altos debido a su forma y terreno pues al no poseer zonas de elevadas
pendientes el agua tiende a estancarse más y, por lo tanto, sus crecidas son atenuadas. Para
el segundo parámetro, se hace la relación entre área y la longitud al cuadrado de la cuenca;
si este es mayor a 0.75, la cuenca tiende a ser cuadrada, y, si es menor a 0.25, la cuenca tiende
a ser alargada, esto debido a que el factor de forma compara la geometría de la cuenca con
un cuadrado. Igualmente, el factor de forma muestra que Kf=0.557 es mayor a 0.25 y más
cercano al 0.75 por lo cual la cuenca tiende a ser cuadrada.

• Caracterización de la red de drenaje

Dentro de la cuenca, el agua que cae en esta se escurre mediante un sistema de corrientes, o
una red de drenaje que lleva toda el agua hacia un mismo punto de salida (ArcGIS, 2019).
Este sistema tiene forma de árbol, debido a que sus corrientes se van ramificando a partir de
una corriente principal. Para cuantificar esta red, Horton definió relaciones para estas
corrientes, también conocidas como índices de Horton.
Orden de la cuenca
A las corrientes se les da un orden acumulativo según su tamaño, siendo el cauce principal
el mayor. Las corrientes comienzan como orden 1, y, si se juntan dos corrientes del mismo
orden, estas forman una corriente de orden mayor. En cambio, si se juntan dos corrientes de
distinto orden, la corriente continuará con el orden mayor.

28
 Ilustración 31: Ordenamiento de cuenca de orden 5

Las corrientes de menor orden, consideradas como corrientes jóvenes, se relacionan con
mayores pendientes en la parte alta de la cuenca, y, por lo tanto, las corrientes de mayor orden
se consideran con menores pendientes debido a que se encuentran en las zonas bajas. Esta
jerarquización influye en gran medida en procesos de erosión, transporte de sedimentos y
deposición. Estos procesos deben analizarse en conjunto con la topografía y las poblaciones
cercanas para evitar posibles desastres y desarrollar sistemas de evacuación o de alerta
temprana, ya que, dependiendo de la bifurcación de las corrientes, la red tendrá una respuesta
bien sea lenta o rápida ante un pulso de precipitación. Para obtener esta información de la
cuenca analizada, y a partir del comando con el comando [Link], se obtuvieron estas
estadísticas de la red de drenaje de Strahler. Se utilizaron los mapas de la red de drenaje de
Strahler, el mapa raster de direcciones y el mapa de elevaciones para obtener las estadísticas
de las corrientes y de la cuenca. El orden mayor obtenido es de 5.

29
 • Relieve de la cuenca

Ilustración 32:Mapa relieve cuenca

Se puede ver a partir del relieve de la cuenca que se puede identificar los puntos más altos de
la región desde donde empiezan las líneas de corriente y bajan en rumbo al punto de
concentración de la cuenca esto se identifica de los puntos que se ven como las partes altas
según las curvas de nivel generadas y coincide con las líneas de corriente que se generan de
la red de drenaje de la cuenca, en esto se puede identificar los puntos desde donde escurre el
agua en la cuenca con rumbo al cause principal.

30
Puente Blanco

• Caracterización de la cuenca hidrográfica.

Por último, se hará la caracterización de la cuenca hidrográfica que se encuentra en el

Puente Blanco el cual se encuentra en Cartago, Colombia (N: 2100000; S: 1990180.77; O:
4650000; E: 4750000, Origen Nacional). Es importante decir que se realizó el mismo
procedimiento que ya se describió en las cuencas anteriores, para obtener los diferentes
parámetros de la cuenca como el Área (A), Longitud de la Cuenca y Cauce Principal (L, Lcp),
Perímetro (P), Pendiente, ancho (S), índice de Capacidad de Gravelius (Kc), factor de forma
(Kf) y el orden de la cuenca.

Ilustración 33: Ubicación puente Blanco

31
 Ilustración 34: Mapa elevación de corriente puente Blanco

Área y perímetro de la cuenca.:

Para la cuenca de Puente Blanco se obtuvo el área y el perímetro, usando el mapa obtenido
del mapa que contiene la línea vectorial de longitud de la cuenca en kilómetros.
𝐴 = 2781.26 𝑘𝑚2
𝑃 = 371.38 𝑘𝑚
Longitud y Ancho de la Cuenca:
Usando las herramientas de medición longitudinal del programa es posible realizar la
medición de la longitud de la cuenca y haciendo uso de la relación entre área y longitud de
la cuenca es posible encontrar el acho de esta
𝐿 = 85.78 𝑘𝑚
𝐴
𝑊= = 32.42 𝑘𝑚
𝐿
Parámetros de forma de la cuenca:
Para los parámetros de forma de la cuenca, se calculó el índice de compacidad de Gravelius
(kc) y el factor de forma (kf), Estos dos valores logra obtener mediante las siguientes
ecuaciones:
𝑃
𝐾𝑐 = 0.282 = 1.99
√𝐴

𝐴 𝑊
𝐾𝑓 = = = 0.378
𝐿2 𝐿

32
Teniendo en cuenta que el Índice de Gravelius es mayor a 1.75 se puede decir que la cuenca
tiene una forma Rectangular. Este parámetro indica la forma en la que la red de drenaje de la
cuenca está organizada, para este caso la red de drenaje es de forma rectangular

Kc Forma de la cuenca
1.00-1.25 Redonda a oval oblonga
1.25-1.50 De oval redonda a oval oblonga
1.50-1.75 De oval oblonga a rectangular oblonga
>1.75 Rectangular
Tabla 5: Índice de Gravelius

Ilustración 35: Formas de drenaje, tomado de: [Link]

• Caracterización de la red de drenaje

Orden de la cuenca
Para poder obtener el orden de la cuenca Puente Blanco se utilizó el comando [Link],
el cual arroja toda la información de la cuenca, obteniendo el orden ( 7 ) , la medida del cauce
principal ( Lc = 146.299 Km ) y otros parámetros importantes como la pendiente promedio
del cauce principal ( Sc = 0.0248 ) , la relación de área, longitud, la densidad de drenaje,
entre otros.

33
 Ilustración 36: Información red de drenaje

Para poder terminar el análisis de la red de drenaje de la cuenca de puente Blanco, es

importante realizar la Curva Hipsométrica la cual indica como va variando la altura con
respecto al área de la cuenca, se puede observar que la cuenca tiene varias pendientes altas,
lo cual puede generar Caudales altos que deben ser tenidos en cuenta en el momento de
diseñar alguna estructura por esta zona.

Ilustración 37: curva hipsométrica Puente Blanco

• Relieve de la cuenca

Haciendo uso del mapa de pendientes generado en el programa, fue posible calcular
la pendiente promedio de la cuenca, para esto se utilizó el comando [Link] el cual

34
 arroja los parámetros estadísticos del mapa de pendientes. El valor medio sería el
valor de la pendiente promedio (S=27.67)

Ilustración 38: Valores estadísticos Mapa de pendientes de elevación de la cuenca

Puente Blanco
Geometria de la cuenca Parametros de forma
2 Indice de
Area (km ) 2781.26 1.99
Gravelius
Perimetro Facot de
371.383 0.378
(km) forma
Longitid Caracterizacion red de
85.78
Cuenca (km) drenaje
Ancho (Km) 32.4232 Orden 7
Longitud
Pendiente
Cauce
Promedio 27.6701 146.299
Principal
Cuenca
(km)
Pendiente
Promedio
0.0367
Cauce
Principal

Ilustración 39: Datos Cuenca Puente Blanco

• Análisis de la amenaza ante crecientes súbitas en los puentes vehiculares

Análisis de precipitación:
Información de precipitación media mensual Anual

Para poder obtener la información de precipitación media mensual y anual de la zona a

analizar, fue necesario tener los datos de los instrumentos de medición de la precipitación
que se encuentran en esta zona, los cuales se encuentran en la página del IDEAM. Al obtener
los diferentes resultados de precipitación para cada mes y para cada año de las diferentes
estaciones, se procedió a subir estos datos al Grass para así poder obtener los diferentes
mapas y la precipitación media dentro de la cuenca. Es importante resaltar que primero se
tiene que subir los datos en WGS84 y luego si Re proyectarlo a Origen Nacional, puesto que
las coordenadas de latitud y longitud esta WGS84.

35
 Precipitacion Media Multianual
latitud longitud jan feb mar apr may jun jul aug sep oct nov dec Anual
4.521 -75.278 93.694 112.274 193.697 233.606 230.948 185.026 157.245 160.790 226.906 220.700 159.932 113.423 2088.242
4.384 -75.497 46.848 58.338 86.616 130.346 141.148 105.619 85.406 81.334 94.307 104.429 87.089 61.068 1082.548
4.007 -75.608 59.465 84.894 120.713 172.174 172.677 126.094 106.535 94.442 120.835 135.765 120.577 83.062 1397.233
4.125 -75.500 95.445 95.429 141.594 189.397 188.658 137.003 154.900 138.445 129.248 197.932 209.374 129.403 1806.829
4.224 -76.147 62.065 65.613 108.645 162.129 129.226 73.452 68.710 52.613 100.387 143.016 136.771 87.677 1190.303
4.589 -75.966 72.823 93.613 148.032 192.323 183.290 109.903 85.903 88.323 133.655 160.748 141.142 90.507 1500.262
4.054 -75.995 150.539 135.968 166.981 238.323 169.716 105.752 55.570 61.613 116.200 289.597 258.777 183.659 1932.693
4.310 -76.020 73.126 81.168 142.165 194.784 175.755 116.968 98.803 82.677 126.961 179.226 168.748 121.148 1561.530
4.172 -76.056 86.816 102.119 170.148 209.397 156.423 70.092 72.837 73.355 123.022 194.339 193.361 130.187 1582.097
4.087 -76.101 104.474 108.442 188.471 215.294 172.274 109.404 80.339 90.795 144.639 218.616 201.652 154.381 1788.779
4.516 -76.041 59.161 68.194 111.839 171.129 142.710 101.194 81.258 71.903 109.313 131.670 111.969 80.777 1241.116
4.114 -76.043 106.484 124.787 185.858 212.571 156.323 88.587 73.004 71.065 127.429 257.545 222.429 167.252 1793.333
4.383 -76.074 49.677 64.884 115.710 161.581 140.552 88.774 85.010 76.984 114.141 138.887 112.639 82.723 1231.560
4.057 -75.806 80.633 76.355 94.774 108.774 118.061 62.200 45.493 36.013 73.000 143.416 152.145 99.667 1090.531
4.688 -75.962 84.574 93.161 131.561 185.000 174.258 107.935 100.639 98.974 130.493 156.952 141.497 91.655 1496.699
4.033 -75.888 73.089 71.629 98.787 118.139 96.710 49.870 38.537 38.625 60.768 147.416 133.977 93.307 1020.854
4.090 -76.224 81.962 88.985 144.623 164.510 138.329 81.233 62.385 64.382 107.886 158.661 141.700 93.181 1327.836
4.028 -76.169 95.877 106.035 149.713 168.812 110.471 64.815 53.128 58.864 124.110 165.635 154.777 107.643 1359.881
4.343 -76.190 47.169 57.774 92.677 152.455 132.935 70.903 64.323 56.282 97.935 125.954 123.694 82.863 1104.964
4.971 -76.039 149.516 157.126 217.687 254.723 252.258 182.952 177.730 171.748 176.813 302.345 292.416 217.981 2553.295
4.794 -75.987 76.616 90.385 157.284 180.077 188.197 105.855 85.877 90.945 138.029 169.413 157.274 102.855 1542.808
4.685 -76.036 85.419 77.371 136.665 173.419 145.323 98.065 97.258 94.161 137.897 144.332 135.074 82.742 1407.726
4.419 -76.101 36.438 61.248 85.677 138.510 126.626 71.617 81.186 63.964 99.139 111.510 102.459 67.080 1045.452
4.532 -76.094 41.129 68.968 105.903 164.810 133.032 77.774 74.452 70.387 102.452 125.495 128.874 66.210 1159.485
4.616 -76.034 40.323 72.862 92.419 145.097 118.781 74.387 85.839 75.645 102.416 107.267 93.832 51.978 1060.846
4.530 -76.060 49.794 67.274 96.198 147.935 133.277 86.687 80.259 74.289 103.981 119.955 106.995 60.464 1127.107
4.481 -76.080 45.900 63.298 93.362 145.475 118.718 69.945 70.210 56.953 96.658 121.927 99.394 65.530 1047.369
4.332 -75.835 122.506 124.303 187.923 226.081 179.571 94.232 69.239 74.177 132.326 249.974 242.161 173.677 1876.171
4.403 -75.915 100.026 117.145 186.655 225.016 191.732 124.697 96.732 88.955 157.313 226.000 187.248 122.545 1824.065
4.677 -75.786 104.037 117.755 191.377 235.674 221.097 157.316 105.910 108.094 168.019 217.600 228.326 148.190 2003.395
4.409 -75.875 99.423 123.645 188.955 229.058 185.903 119.806 96.883 88.645 152.977 240.171 191.355 121.306 1838.128
4.410 -76.882 139.455 147.145 191.884 264.400 223.858 155.326 108.516 96.387 180.197 255.132 271.632 185.429 2219.361
4.185 -75.832 156.432 158.471 233.271 293.239 249.335 143.329 103.216 96.106 187.881 339.306 289.748 219.006 2469.342
4.844 -75.558 172.558 171.290 245.765 307.065 267.745 184.890 132.813 152.781 212.800 331.400 300.755 229.748 2709.610
4.886 -75.558 161.684 161.519 220.748 261.948 216.910 148.874 121.810 137.358 169.274 321.565 305.313 213.971 2440.974
4.801 -75.858 102.687 126.484 171.355 211.226 193.806 132.516 116.645 106.516 146.387 206.597 203.771 135.758 1853.749
4.962 -75.742 152.706 156.816 226.216 289.016 275.045 235.077 206.803 181.810 220.226 252.132 273.029 187.539 2656.416
4.812 -75.843 116.161 134.852 194.084 228.894 210.606 136.823 121.548 119.065 164.139 218.326 224.955 157.539 2026.990
4.816 -75.737 133.268 145.648 237.658 250.832 248.265 189.813 134.800 129.787 193.261 289.719 289.816 204.632 2447.500
4.893 -75.883 98.641 116.032 181.942 217.142 243.419 164.581 129.323 116.355 148.581 195.171 216.687 140.077 1967.950
4.952 -75.490 108.452 124.677 166.903 213.648 206.968 144.845 97.468 99.132 151.794 251.023 240.881 139.881 1945.671
4.898 -76.677 309.287 279.971 306.897 355.855 365.568 266.997 240.635 278.590 310.071 493.674 511.374 400.619 4119.539
4.573 -76.202 68.645 92.484 116.903 175.290 133.871 92.710 71.871 81.806 132.902 176.845 160.406 89.217 1392.951

Tabla 6: Datos de precipitación media multianual mensual y anual

Precipitación media / Ciclo Interanual

Teniendo en cuenta los datos anteriores se procedió a realizar el cálculo de la precipitación
media mensual y anual en cada una de las estaciones presentadas anteriormente y a partir de
estos datos se obtuvo la precipitación promedio/media mensual de la zona a analizar la cual
es de 144 mm/mes
Para poder realizar el ciclo interanual, se tomaron los datos la estación con coordenadas
4.6884, -75.9617 y se procedió a realizar el grafico correspondiente, donde se encontró que
entre los meses marzo, abril y mayo es donde la precipitación media es mayor que los otros
meses, esto tiene validez, puesto que en estos meses normalmente se presenta el fenómeno
de la niña. Los meses que menos precipitación media hay en esta zona son en los meses de
junio, julio, agosto los cuales tiene sentido puesto que en esa época está presente el fenómeno
del niño

36
 Ilustración 40: ciclo interanual

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

Teniendo los datos de precipitación de la zona y la limitación de la cuenca a analizar, es

posible obtener la precipitación media multianual mensual y anual, haciendo una
interpolación del mapa raster de la cuenca y el mapa vectorial de las diferentes estaciones.
Haciendo uso del comando [Link]. se encuentra los parámetros estadísticos de los diferentes
meses y de la precipitación media anual dentro de la cuenca.

Precipitación Media
dentro de la cuenca
(mm/mes)
Enero 104.804
Febrero 116.281
Marzo 176.718
Abril 218.251
Mayo 195.019
Junio 127.480
Julio 98.949
Agosto 96.549
Septiembre 149.082
Octubre 222.926
Noviembre 206.212
Diciembre 141.949
Anual 1854.220
Tabla 7: Precipitación Media dentro de la cuenca

37
 Ilustración 41: Mapa Precipitación enero

Evaporación:
Temperatura promedio multianual anual en cada estación

De la misma forma que se procesaron los datos de la precipitación se realizó con los datos de
temperaturas obtenidos por la página del IDEAM. Para este caso, es importante buscar las
estaciones que también tengan los datos de precipitación, para poder realizar los cálculos de
la evaporación real, potencial y el índice de aridez de la zona y dentro de la cuenca. Es
importante resaltar que el cálculo para la temperatura media anual se obtuvo a partir de la
media de las diferentes temperaturas registradas a lo largo del año.

Temperatura Media multianual mensual

latitud longitud jan feb mar apr may jun jul aug sep oct nov dec Anual

4.03 -75.89 11.75 11.81 11.94 12.02 12.26 12.14 11.78 11.88 11.84 11.46 11.38 11.56 11.82
4.24 -76.03 23.30 23.43 23.49 23.41 23.36 23.39 23.41 23.68 23.59 23.08 22.91 23.15 23.35
4.09 -76.22 24.03 24.28 24.24 23.96 23.79 23.72 24.04 24.15 23.97 23.40 23.19 23.54 23.86
4.03 -76.17 22.56 22.77 22.70 22.56 22.51 22.62 22.68 22.98 22.72 22.31 22.28 22.32 22.58
4.53 -76.06 24.40 24.77 24.80 24.41 24.20 24.18 24.41 24.69 24.37 23.82 23.79 24.00 24.32
4.45 -75.77 22.22 22.37 22.25 22.03 22.05 22.21 22.46 22.59 21.97 21.48 21.51 21.85 22.08
4.19 -75.83 18.96 18.96 18.87 18.78 18.89 19.06 19.29 19.46 19.02 18.48 18.44 18.68 18.91
4.82 -75.74 21.89 22.06 21.92 21.59 21.51 21.58 21.83 21.99 21.58 20.98 20.99 21.42 21.61

Tabla 8: Temperatura media multianual mensual y anual

Evaporación potencial
Para poder hacer el cálculo de la evaporación potencial se utilizó la ecuación de thornthwaite,
la cual requiere del cálculo del índice calórico, la constante empírica y tener en cuenta el
parámetro de ajuste por número de horas de luz. Haciendo uso de estos parámetros y de la

38
temperatura promedio en cada mes, se procedió a realizar el cálculo de la evaporación
potencial para cada mes.

Ilustración 42: Método de Thornthwaite

Evapotranspiración Potencia ( Etpr ) (mm/mes) / Evaporación Potencial Mensual
Latitud longitud Ene Febr Mar Abr May Jun Jul Ago Sept Oct Nov Dic Anual
4.03305556 -75.888 55.865 50.793 56.913 55.679 58.665 56.309 56.074 56.581 54.745 54.349 52.342 54.854 663.169
4.23891667 -76.033 96.883 89.257 100.331 97.665 100.641 98.541 101.276 104.225 99.346 95.156 90.126 94.778 1168.225
4.09008333 -76.224 104.346 97.416 108.265 102.938 104.638 101.257 107.826 108.966 102.615 97.253 91.547 97.786 1224.854
4.02833333 -76.169 91.010 84.448 93.628 90.484 93.327 92.179 95.159 98.147 91.766 89.044 85.807 88.119 1093.118
4.53 -76.060 107.998 102.505 114.784 107.452 108.670 105.709 111.501 115.028 106.604 101.192 97.387 102.338 1281.168
4.45472222 -75.766 88.794 81.894 90.339 86.677 90.166 89.631 94.494 95.496 85.750 82.245 79.762 84.813 1050.060
4.18505556 -75.832 71.246 64.543 71.501 69.643 73.137 72.675 76.274 77.285 71.135 68.296 65.696 68.880 850.311
4.81586111 -75.737 87.100 80.452 88.485 83.823 86.193 84.825 89.403 90.657 83.417 79.258 76.718 82.304 1012.635

Tabla 9: Evaporación Potencial

Índice de aridez y Evaporación Real

Haciendo uso de la evaporación potencial obtenido anteriormente, se utilizaron las

ecuaciones de budyco, el cual indica que el índice de aridez es la división entre la evaporación
potencial sobre la precipitación media multianual mensual o anual. Al obtener los datos de
los índices de aridez mensual y anual, se procedió a realizar el cálculo de la evaporación real,
haciendo uso de la formula del método a utilizar.
Indice de Aridez promedio multianual mensual
Latitud longitud Ene Febr Mar Abr May Jun Jul Ago Sept Oct Nov Dic Anual
4.03305556 -75.888 0.76 0.71 0.58 0.47 0.61 1.13 1.46 1.46 0.90 0.37 0.39 0.59 0.65
4.23891667 -76.033 1.18 1.00 0.69 0.59 0.73 1.21 1.62 1.62 0.92 0.60 0.64 1.02 0.88
4.09008333 -76.224 1.09 0.92 0.72 0.61 0.95 1.56 2.03 1.85 0.83 0.59 0.59 0.91 0.90
4.02833333 -76.169 1.83 1.26 0.97 0.61 0.70 1.06 1.19 1.32 0.88 0.74 0.80 1.46 0.97
4.53 -76.060 0.77 0.70 0.60 0.41 0.49 0.68 1.03 1.19 0.59 0.40 0.36 0.55 0.58
4.45472222 -75.766 0.57 0.52 0.39 0.30 0.36 0.63 0.92 0.99 0.46 0.24 0.28 0.39 0.43
4.18505556 -75.832 0.53 0.44 0.30 0.28 0.29 0.38 0.57 0.60 0.37 0.24 0.23 0.34 0.35
4.81586111 -75.737 0.28 0.29 0.29 0.24 0.24 0.32 0.37 0.33 0.27 0.16 0.15 0.21 0.25

Tabla 10: Índice de Aridez

Evapotranspiración Real ( Er ) (mm/mes) / Evaporacion Real
Latitud longitud Ene Febr Mar Abr May Jun Jul Ago Sept Oct Nov Dic Anual
4.03305556 -75.888 43.414 40.498 48.101 49.011 48.951 36.711 31.428 31.573 39.854 49.486 47.339 46.135 543.123
4.23891667 -76.033 61.590 61.839 80.559 81.941 79.565 61.730 52.736 54.378 71.628 79.632 74.254 65.217 859.174
4.09008333 -76.224 69.374 70.309 85.766 85.777 74.486 54.134 47.653 51.637 77.416 81.819 76.879 70.928 891.765
4.02833333 -76.169 43.535 51.862 65.810 75.339 74.714 62.027 60.390 58.457 67.400 69.915 65.507 49.336 769.626
4.53 -76.060 83.530 82.197 96.121 96.698 95.165 85.386 76.340 72.732 89.519 91.348 88.939 87.352 1082.293
4.45472222 -75.766 75.309 70.875 81.790 80.556 82.270 74.190 68.306 66.464 75.881 77.484 74.523 76.787 939.162
4.18505556 -75.832 61.233 57.378 66.359 65.029 67.991 65.887 64.738 64.793 64.781 64.449 62.134 63.304 779.068
4.81586111 -75.737 81.246 74.933 82.394 79.104 81.336 78.369 81.333 83.579 78.063 76.180 73.928 78.248 953.199

Tabla 11: Evaporación Real

Teniendo estos datos organizados, se subieron al Grass, para poder obtener los diferentes
mapas de la zona y tener la información multianual mensual y anual de evaporación

39
potencial, real y el índice de aridez dentro de la cuenca. Esto fue posible usando la
interpolación del mapa de la zona y los diferentes datos obtenidos anteriormente.

Ilustración 43:Ecuación de Budyco

Evaporacion Real Indice Aridez Medio Evaporacion Potencial

Media mensual mensual mutianual Media mensual
mutianual (mm/mes) (mm/mes) mutianual (mm/mes)
Enero 71.746 Enero 0.567 Enero 85.147
Febrero 67.663 Febrero 0.510 Febrero 78.512
Marzo 78.007 Marzo 0.397 Marzo 86.923
Abril 76.746 Abril 0.318 Abril 83.307
Mayo 78.183 Mayo 0.369 Mayo 86.304
Junio 71.382 Junio 0.582 Junio 85.194
Julio 67.955 Julio 0.814 Julio 89.383
Agosto 67.594 Agosto 0.853 Agosto 90.748
Septiembre 72.909 Septiembre 0.465 Septiembre 83.153
Octubre 74.270 Octubre 0.267 Octubre 79.615
Noviembre 71.523 Noviembre 0.276 Noviembre 76.847
Diciembre 72.759 Diciembre 0.404 Diciembre 81.466
Anual 894.011 Anual 0.426 Anual 81.466

Tabla 12 Índice de aridez, evaporación real y potencial dentro de la cuenca

40
 Ilustración 44:Evaporación Real Enero

Ilustración 45: Evaporación Potencial Enero

41
 Ilustración 46: índice de largo plazo en la cuenca

Balance Hídrico:

Balance hidrológico de largo plazo en la cuenca

Teniendo en cuenta los datos encontrados anteriormente, es posible mediante un balance

hídrico encontrados los caudales modelados de la cuenca y así poder camparlos con los
caudales observados.
Para este caso se hizo un balance hídrico el cual consiste en restar la precipitación obtenida
con la evaporación real que hay dentro de la cuenca y multiplicarla por el área total de la
cuenca, esto sería igual al caudal modelado promedio en cada mes.
Donde se puede observar que en los meses donde la precipitación es mayor que la
evaporación se tienen caudales más altos que otros meses, por lo que genera que exista un
aumento en el caudal promedio del rio, el cual es de suma importancia en el momento de
diseñar una infraestructura en esta zona pues que al no tenerse en cuenta esta puede provocar
consecuencias negativas como pérdidas económicas, fallos en la estructura o perdidas de
vida.

42
 Evaporacion Real Caudal Modelado
Precipitacón (m/s) 3
(m/s) (m /S)

Enero 3.99E-08 Enero 2.73E-08 Enero 34.99

Febrero 4.42E-08 Febrero 2.57E-08 Febrero 51.45
Marzo 6.72E-08 Marzo 2.97E-08 Marzo 104.47
Abril 8.30E-08 Abril 2.92E-08 Abril 149.76
Mayo 7.42E-08 Mayo 2.97E-08 Mayo 123.65
Junio 4.85E-08 Junio 2.72E-08 Junio 59.37
Julio 3.77E-08 Julio 2.59E-08 Julio 32.80
Agosto 3.67E-08 Agosto 2.57E-08 Agosto 30.64
Septiembre 5.67E-08 Septiembre 2.77E-08 Septiembre 80.62
Octubre 8.48E-08 Octubre 2.83E-08 Octubre 157.33
Noviembre 7.85E-08 Noviembre 2.72E-08 Noviembre 142.54
Diciembre 5.40E-08 Diciembre 2.77E-08 Diciembre 73.23
Anual 7.06E-07 Anual 3.40E-07 Anual 1016.21
Tabla 13:Balance Hídrico

Conclusión:

Para finalizar, es importante concluir diferentes aspectos. El primer aspecto para concluir es
que es de suma importancia conocer la geomorfológica de la cuenca, puesto que esto da un
indicio de cómo va a hacer la respuesta hidrología de la cuenca a diversos eventos climáticos.
Además, al conocer la geomorfológica de la cuenta se pueden prevenir varias catástrofes
ambientales como la inundación de la cuenca y por ende de las zonas aledañas. Es importante
resaltar que al tener una cuenca de mayor tamaño, se tiene un orden mayor de la cuenca y
por ende los caudales y los diferentes aspectos al ser mayores pueden tener una influencia en
la zona mucho mayor que las cuencas de menor tamaño.
En segunda instancia es importante resaltar que es de suma importancia hacer un análisis de
precipitación y evaporación real media multianual mensual y anual para la zona, puesto que
a partir de estos datos se pueden obtener caudales modelados parecidos a los que se observan
y poder saber en qué periodos de tiempo se van a tener mayores caudales, para así poder
diseñar y realizar una obra de infraestructura como puente y no tenga inconvenientes de
inundación u otros fenómenos que se pueden genera a partir de los grandes caudales.

Bibliografía:
(Salgar, 2021) Cuenca del río magdalena, Proyecto hidrología 2021-30.

43