1.

ESTUDIOS HIDROLÓGICOS E HIDRAULICOS

Y DE GEOMORFOLOGÍA – AVENIDAS
TORRENCIALES PARA LA PTAR EL “PARAÍSO”

1-1

1
Tabla de Contenido

1. ESTUDIOS HIDROLÓGICOS E HIDRAULICOS Y DE

GEOMORFOLOGÍA – AVENIDAS TORRENCIALES PARA LA PTAR
EL “PARAÍSO” ...................................................................................................... 1

1.1 Estudio hidrológico y de crecientes pata la construcción de la PTAR ................. 12

1.1.1 INFORMACIÓN BASE ................................................................................................. 12

1.1.1.1 Información secundaria ............................................................................. 12

1.1.1.2 Información primaria .................................................................................. 16

1.2 Climatología ................................................................................................................ 18

1.3 Morfometría de las cuencas ...................................................................................... 18

1.4 Análisis temporal y espacial de las precipitaciones ............................................... 20

1.4.1 Precipitaciones medias .............................................................................................. 20

1.4.1.2 Detección datos anómalos ..................................................................... 1-23

1.4.1.3 Análisis de homogeneidad y tendencia .................................................. 1-24

1.4.2 Precipitaciones máximas ........................................................................................ 1-27

1.4.3 Elaboración de curvas Intensidad Frecuencia Duración..................................... 1-28

1.4.3.2 Curvas IFD estación Termales ............................................................... 1-32

1.4.3.3 Curvas IFD estación Boquerón .............................................................. 1-33

1.4.3.4 Curvas IFD estación Playa Rica ............................................................ 1-34

1.4.3.5 Curvas IFD estación Pez Fresco ........................................................... 1-35

1.4.3.6 Curvas IFD estación Aeropuerto Matecaña ........................................... 1-36

1.4.3.7 Curvas IFD estación La Laguna Otún .................................................... 1-37

1.5 Análisis temporal y espacial de la Temperatura .................................................. 1-38

1.6 Análisis temporal y espacial de la Evaporación ................................................... 1-41

1.7 Caracterizar régimen hidrológico y Oferta hídrica ............................................... 1-46

1.8 Régimen de Caudales del Río Otún ....................................................................... 1-50

1-2

2
1.8.1 Caudales medios río Otún ...................................................................................... 1-50

1.8.2 Caudales mínimos río Otún .................................................................................... 1-52

1.8.3 Caudal ecológico río Otún ...................................................................................... 1-56

1.8.4 Caudales máximos río Otún ................................................................................... 1-57

1.8.4.1 Análisis de frecuencia ............................................................................ 1-57

1.8.4.2 Modelo lluvia escorrentía en HEC - HMS .............................................. 1-60

1.8.5 Análisis eventos extremos ..................................................................................... 1-73

1.8.5.1 Eventos torrenciales ............................................................................... 1-73

1.8.5.2 Eventos mínimos .................................................................................... 1-77

1.9 Modelamiento Hidráulico ........................................................................................ 1-78

1.9.1 Metodología para estimar los perfiles de flujo ..................................................... 1-78

1.9.2 Consideraciones para el modelado ....................................................................... 1-78

1.9.2.1 Rugosidad .............................................................................................. 1-78

1.9.2.2 Coeficientes de contracción y expansión ............................................... 1-82

1.9.2.3 Condiciones de frontera y Leeves .......................................................... 1-82

1.9.3 Geometría del modelo ............................................................................................. 1-83

1.9.4 Casos de estudio ..................................................................................................... 1-85

1.9.5 Análisis de resultados hidráulicos ........................................................................ 1-86

1.9.6 Conclusiones ........................................................................................................... 1-94

1.10 Geomorfología. ........................................................................................................ 1-97

1.10.2 Subunidades geomorfológicas – Geomorfología local. .................................... 1-102

1.10.3 Abanico de Pereira (Dap). ..................................................................................... 1-103

1.10.4 Abanico incisado (Dapi). ....................................................................................... 1-104

1.10.5 Depósito de vertiente (Qdv). ................................................................................. 1-105

1.10.5.1 Depósito de vertiente (Caída de bloques) del lote San Cayetano ....... 1-106

1-3

3
1.10.5.2 Depósito Coluvial del lote El Paraíso ................................................... 1-107

1.10.6 Terrazas aluviales. ................................................................................................. 1-107

1.10.6.1 Terraza cero (T0) .................................................................................. 1-107

1.10.6.2 Terraza uno (T1) .................................................................................. 1-108

1.10.6.3 Terraza dos (T2) ................................................................................... 1-110

1.10.6.4 Terraza tres (T3). .................................................................................. 1-111

1.10.6.5 Terraza cuatro (T4). .............................................................................. 1-111

1.11 Avenidas Torrenciales. ......................................................................................... 1-113

1.11.2 Torrencialidad histórica. ....................................................................................... 1-116

1.11.3 Estudios anteriores. .............................................................................................. 1-116

1.11.4 Antecedentes históricos ....................................................................................... 1-117

1.11.5 Otros eventos recientes. ....................................................................................... 1-118

1.11.6 Análisis de torrencialidad de la cuenca. ............................................................. 1-119

1.11.6.2 Otún alto. .............................................................................................. 1-120

1.11.6.3 Río Bárbaro. ......................................................................................... 1-122

1.11.6.4 Río San Juan. ....................................................................................... 1-124

1.11.6.5 Otún medio. .......................................................................................... 1-125

1.11.6.6 Río San José. ....................................................................................... 1-127

1.11.6.7 Quebrada Dosquebradas. .................................................................... 1-128

1.11.6.8 Otún Bajo.............................................................................................. 1-129

1.11.6.9 Quebrada la arenosa. ........................................................................... 1-130

1.12 Conclusiones Ronda Hídrica ................................................................................ 1-131

BIBLIOGRAFÍA 1-132

1-4

4
Lista de Tablas

Tabla 1-1 Relación de estaciones de medición climática utilizada en el estudio

POMCA del río Otún. ............................................................................................... 13

Tabla 1-2 Lista de estaciones suministradas por IDEAM en el 2019 para el estudio ............. 16

Tabla 1-3 Resumen de parámetros morfométricos de las cuencas de interés ....................... 20

Tabla 1-4 Lista de las estaciones de medición de precipitación total empleadas en el

estudio ..................................................................................................................... 21

Tabla 1-5 Análisis de anómalos para las estaciones utilizando la metodología de

DIXON y GRUBBS ............................................................................................... 1-24

Tabla 1-6 Análisis de homogeneidad y cambio en la media ................................................ 1-25

Tabla 1-7 Análisis de tendencia de las estaciones ............................................................... 1-25

Tabla 1-8 Precipitación máxima en 24 horas en distintos períodos de retorno para las
estaciones solicitadas ........................................................................................... 1-28

Tabla 1-9 Valores de los coeficientes a, b, c y d para el cálculo de las curvas

intensidad-duración-frecuencia, IDF, para Colombia ........................................... 1-29

Tabla 1-10 Intensidad Frecuencia Duración de la estación Termales ................................... 1-32

Tabla 1-11 Intensidad Frecuencia Duración de la estación Boquerón ................................... 1-33

Tabla 1-12 Intensidad Frecuencia Duración de la estación Playa Rica ................................. 1-34

Tabla 1-13 Intensidad Frecuencia Duración de la estación Pez Fresco ................................ 1-35

Tabla 1-14 Intensidad Frecuencia Duración de la estación Aeropuerto Matecaña ................ 1-36

Tabla 1-15 Intensidad Frecuencia Duración de la estación La Laguna Otún......................... 1-37

Tabla 1-16 Magnitud de las temperaturas promedios en la estación Aero. Matecaña

[26135040] ............................................................................................................ 1-39

Tabla 1-17 Gradientes de variación de temperatura por cada 100 m para la región
andina ................................................................................................................... 1-40

Tabla 1-18 Comparativo entre los métodos de estimación de la evaporación y la

estación Apto Matecaña ....................................................................................... 1-45

Tabla 1-19 Balance entre precipitación y evaporación en la estación aeropuerto

Matecaña .............................................................................................................. 1-47

1-5

5
Tabla 1-20 Oferta y balance hídrico multianual de las cuencas de la zona ........................... 1-49

Tabla 1-21 Caudales mínimos en la estación La Bananera ................................................... 1-55

Tabla 1-22 Caudales mínimos estimados en la PTAR ........................................................... 1-56

Tabla 1-23 Caudales medios mensuales multianuales registrados en la estación La

Bananera .............................................................................................................. 1-57

Tabla 1-24 Pruebas teóricas de ajuste en orden de jerarquía para caudales ........................ 1-58

Tabla 1-25 Caudales máximos en la estación La Bananera para diferentes periodos

de retorno.............................................................................................................. 1-58

Tabla 1-26 Caudales máximos (m 3/s) en la estación La Bananera de estudios previos

y actual.................................................................................................................. 1-59

Tabla 1-27 Caudales máximos (m 3/s) del río Otún en la desembocadura con el río
Cauca.................................................................................................................... 1-59

Tabla 1-28 Métodos utilizados para el modelo HEC-HMS ..................................................... 1-60

Tabla 1-29 Parámetros morfométricos de las subcuencas .................................................... 1-61

Tabla 1-30 Tiempos de concentración para cada una de las subcuencas ............................ 1-64

Tabla 1-31 Número de curva para diferentes usos del suelo según condiciones de
humedad antecedente Tipo II ............................................................................... 1-66

Tabla 1-32 Numero de curva de las cuencas de análisis ....................................................... 1-69

Tabla 1-33 Precipitaciones máximas y % de área aferente ................................................... 1-70

Tabla 1-34 Precipitación máximas en mm para diferentes duraciones y periodos de

retorno................................................................................................................... 1-70

Tabla 1-35 Caudales máximos estimados con el HEC-HMS ................................................. 1-72

Tabla 1-36 Caudales máximos para diferentes periodos de retorno en la PTAR .................. 1-72

Tabla 1-37 Resumen de eventos hidráulicos, precipitaciones, desbordamientos

cubiertos por la prensa ......................................................................................... 1-74

Tabla 1-38 Coeficiente n1 en función del grado de irregularidad de la sección ..................... 1-80

Tabla 1-39 Coeficiente n2 en función de la variación en la sección transversal del canal

principal................................................................................................................. 1-80

1-6

6
Tabla 1-40 Coeficiente n3 en función del grado de obstrucciones presentes en el canal
principal................................................................................................................. 1-81

Tabla 1-41 Coeficiente n4 en función del grado de obstrucciones presentes en el canal

principal debido a la vegetación ........................................................................... 1-81

Tabla 1-42 Coeficiente m de corrección por meandros en el tramo de estudio ..................... 1-81

Tabla 1-43 Rugosidad seleccionada para el canal principal .................................................. 1-82

Tabla 1-44 Niveles máximos estimados en la PTAR .............................................................. 1-89

Tabla 1-45 Diferencias entre depósitos aluviales y torrenciales .......................................... 1-113

Tabla 1-46 Eventos de crecientes más relevantes ............................................................... 1-118

1-7

7
Lista de Figuras

Figura 1-1. Distribución de las subcuencas del río Otún analizadas en el plan de
manejo de la cuenca (POMCA) ............................................................................... 15

Figura 1-2. Distribución espacial de las cuencas de los afluentes más importantes del
río Otún .................................................................................................................... 19

Figura 1-3. Análisis de anómalos por metodología de Grubbs (izq.) y Dixon (der.) para
la estación Pez fresco ........................................................................................... 1-23

Figura 1-4. Análisis de anómalos por metodología de Grubbs (izq.) y Dixon (der.) para
la estación Termales ............................................................................................. 1-23

Figura 1-5. Análisis de anómalos por metodología de Grubbs (izq.) y Dixon (der.) para
la estación Apto Matecaña ................................................................................... 1-24

Figura 1-6. Promedio de la precipitación total mensual de las estaciones

representativas de la cuenca del río Otún ............................................................ 1-26

Figura 1-7. Campo de precipitación promedio multianual sobre la cuenca del río Otún ........ 1-27

Figura 1-8. Curvas IDF para 100 años de periodo de retorno por diferentes
metodologías ........................................................................................................ 1-31

Figura 1-9. Polígonos de Thiessen ......................................................................................... 1-32

Figura 1-10. Curvas Intensidad Frecuencia Duración de la estación Termales ....................... 1-33

Figura 1-11. Curvas Intensidad Frecuencia Duración de la estación Boquerón ...................... 1-34

Figura 1-12. Curvas Intensidad Frecuencia Duración de la estación Playa Rica ..................... 1-35

Figura 1-13. Curvas Intensidad Frecuencia Duración de la estación Pez Fresco .................... 1-36

Figura 1-14. Curvas Intensidad Frecuencia Duración de la estación Aeropuerto

Matecaña .............................................................................................................. 1-37

Figura 1-15. Curvas Intensidad Frecuencia Duración de la estación La Laguna Otún ............ 1-38

Figura 1-16. Variación mensual de la temperatura promedio mínima, media y máxima ......... 1-39

Figura 1-17. Variación espacial de la temperatura promedio multianual .................................. 1-41

Figura 1-18. Variación de la evaporación promedio mensual en la estación Apto.

Matecaña .............................................................................................................. 1-44

Figura 1-19. Evapotranspiración real estimada con el método de Turc ................................... 1-46

1-8

8
Figura 1-20. Variación mensual promedio de la precipitación y la evaporación de tanque
en la estación Aeropuerto Matecaña .................................................................... 1-47

Figura 1-21. Balance entre Precipitación-Evapotranspiración real total anual ......................... 1-48

Figura 1-22. Caudal medio diario registrado en la estación La Bananera ................................ 1-50

Figura 1-23. Caudal medio mensual registrado en la estación La Bananera ........................... 1-51

Figura 1-24. Curva de duración de caudales en la estación La Bananera ............................... 1-52

Figura 1-25. Caudales mínimos estimados en La Bananera .................................................... 1-55

Figura 1-26. Caudales máximos para diferentes ajustes .......................................................... 1-58

Figura 1-27. Subcuencas del modelo HEC-HMS ..................................................................... 1-61

Figura 1-28. Coberturas CLC para las cuencas de análisis ..................................................... 1-68

Figura 1-29. Coberturas simplificadas CLC para las cuencas de análisis ................................ 1-69

Figura 1-30. Cuencas hidrográficas en HEC-HMS ................................................................... 1-71

Figura 1-31. Hidrógrafas para diferentes periodos de retorno .................................................. 1-72

Figura 1-32. Crecientes del río Otún a la altura del viaducto César Gaviria Trujillo.
Tomada diario El Mundo, 15 noviembre de 2017. ............................................... 1-74

Figura 1-33. Eventos torrenciales del río Otún y Consota en noviembre de 2017 ................... 1-77

Figura 1-34. Material que compone el lecho del cauce del río Otún ........................................ 1-79

Figura 1-35. Localización general de las secciones y del tramo de estudio ............................. 1-83

Figura 1-36. Esquema de detalle del modelo en la zona de la PTAR ...................................... 1-84

Figura 1-37. Isométrico del modelo realizado en HEC-RAS..................................................... 1-85

Figura 1-38. Perfil de flujo de línea base para la creciente de 2,33 años de periodo de
retorno................................................................................................................... 1-86

Figura 1-39. Perfiles de profundidad del flujo para diferentes crecientes ................................. 1-87

Figura 1-40. Velocidades del flujo para diferentes crecientes .................................................. 1-88

Figura 1-41. Ancho superficial para diferentes crecientes ........................................................ 1-89

Figura 1-42. Mancha de inundación para la creciente de 2,33 años de periodo de retorno
.............................................................................................................................. 1-90

1-9

9
Figura 1-43. Mancha de inundación para la creciente de 100 años de periodo de retorno
.............................................................................................................................. 1-91

Figura 1-44. Mancha de inundación para la creciente de 100 años de periodo de retorno
mayorada .............................................................................................................. 1-92

Figura 1-45. Manchas de inundación para otros escenarios poco probables .......................... 1-94

Figura 1-46. Geomorfología regional ........................................................................................ 1-98

Figura 1-47. Modelo digital de terreno del sector occidental del río Otún .............................. 1-100

Figura 1-48. Subunidades geomorfológicas ........................................................................... 1-103

Figura 1-49. Perfil de terraza Abanico de Pereira (Qap) margen derecha del río Otún ......... 1-104

Figura 1-50. Perfil de terrzas QT1 y QT2 margen izquierda del río Otún ............................... 1-109

Figura 1-51. Perfil de terrza QT4 margen derecha del río Otún .............................................. 1-112

Figura 1-52. Geoformas área de estudio ................................................................................ 1-112

Figura 1-53. Amenaza por avenida torrencial ......................................................................... 1-115

Figura 1-54. Delimitación de la cuenca y sitios estudiados sobre base IGAC 1:200.000 ...... 1-120

Figura 1-55. Morrenas y cicatrices de deslizamientos de la parte alta del Otún – Barbo ...... 1-121

Figura 1-56. Rasgos morfodinámicos y de depositación, Otún, sitio La Pastora ................... 1-122

Figura 1-57. Huella del flujo torrencial o creciente de mayo de 2007 ..................................... 1-123

1-10

10
Lista de Fotografías

Fotografía 3-1. Vista parcial del lote Valparaíso con el Abanico de Pereira con escarpe
y Terrazas aluviales ............................................................................................ 1-105

Fotografía 3-2. Depósito de vertiente por caída de bloques en el lote San Cayetano ........ 1-106

Fotografía 3-3. Aspecto de la terraza T0, en el borde del río .............................................. 1-108

Fotografía 3-4. Composición de los depósitos aluviales en el sitio del proyecto ................ 1-109

Fotografía 3-5. Superficie de la Terraza T1, con auger que muestra el espesor del
horizonte A > 25 cm ............................................................................................ 1-110

Fotografía 3-6. Disposición paralela al río de la Terraza T2, lote San Cayetano ................ 1-111

Fotografía 10-1. Cauce del río Barbo en Pezfresco luego de la creciente del 14 de
octubre de 2011 .................................................................................................. 1-124

Fotografía 10-2. Huellas de creciente San Juan el 14/10/2011. B) el Otún urbano el

mismo evento ..................................................................................................... 1-125

Fotografía 10-3. a) Terraza T2 típica del cauce medio (urbano) del Otún. b) Marca
de erosión reciente en la margen del Otún, sector La Bananera ....................... 1-126

Fotografía 10-4. Aspecto de la confluencia del río San José al Otún................................ 1-127

Fotografía 10-5. Depósito que conforma la Terraza-abanico de la Q. Dosquebradas ...... 1-128

Fotografía 10-6. Sistema de terrazas aluviales del Otún en Estación Pereira .................. 1-129

1-11

11
1.1 Estudio hidrológico y de crecientes pata la construcción
de la PTAR
1.1.1 INFORMACIÓN BASE

Este numeral está discriminado en información secundaria y primaria, a continuación, se presenta

una breve descripción para cada ítem.

1.1.1.1 Información secundaria

El estudio principal tomado como referencia corresponde al Plan de Ordenamiento y Manejo de

la Cuenca (POMCA) del río Otún, elaborado por el consorcio Ordenamiento cuenca del río Otún
en el año 2007.

También se consideró, pero en menor medida el componente hídrico e hidráulico del estudio de
riesgo del departamento de Risaralda, adelantado por la Corporación Autónoma Regional de
Risaralda (CARDER,1998).

Es de resaltar que, para la elaboración de los componentes hidrológicos y morfométricos del

POMCA, el estudio consideró la base de datos de registros de estaciones de medición
presentados en la fase de aprestamiento. A continuación, se presenta en forma de resumen la
información de las estaciones consideradas de gran relevancia en el POMCA:

1-12

12
 Tabla 1-1 Relación de estaciones de medición climática utilizada en el estudio POMCA del río Otún.
Entidad Código Estación Tipo Estado Longitud Latitud Inicio Fin Años
CHEC 26120310 BREMEN PG SUSPENDIDA 4,68333 -75,53333 1977 1999 29
SUB ESTA
CHEC 26130060 4,88333 -75,63333 1960 33
ROSA PM ACTIVA
CHEC 26130070 SAN RAMON PM ACTIVA 4,85000 -75,58333 1970 13
SUB ESTA
CHEC 26130120 4,83333 -75,70000 1975 17
ROSA PM ACTIVA
CHEC 26130130 LA FLORIDA PM ACTIVA 4,90000 -75,66667 1970 17
FNC 26125240 LA CATALINA CP ACTIVA 4,75000 -75,75000 1986 9
LOS
FNC 26130210 4,81667 -75,83333 1968 32
CAMBULOS PM ACTIVA
P
FNC 26130230 4,80000 -75,66667 1970 28
TRATAMIENTO PG ACTIVA
FNC 26130260 EL BOSQUE PM ACTIVA 4,83333 -75,68333 1968 18
FNC 26130310 LA MORENITA PM ACTIVA 4,90000 -75,76667 1991 11
FNC 26130510 COMBIA PM ACTIVA 4,85000 -75,78333 1991 15
FNC 26135060 EL JAZMIN CO ACTIVA 4,91667 -75,63333 1960 37
FNC 26135070 EL CEDRAL CO ACTIVA 4,78333 -75,53333 1960 30
IDEAM 21210180 TOCHE PG ACTIVA 4,51950 -75,40983 1986 27
IDEAM 26100830 PALMASOLA PM ACTIVA 4,68842 -75,96167 1980 32
LA
IDEAM 26110150 4,79383 -75,98650 1946 44
INMACULADA PM MANT.
IDEAM 26115040 LA UNION CP ACTIVA 4,53122 -76,06217 1967 24
IDEAM 26120150 ALCALÁ PM ACTIVA 4,67647 -75,78631 1972 41
IDEAM 26120160 SALENTO PM ACTIVA 4,63628 -75,56917 1974 39
IDEAM 26130020 TERMALES PM ACTIVA 4,84417 -75,55806 1971 41
IDEAM 26130170 POTREROS PM ACTIVA 4,88556 -75,55806 1970 43
IDEAM 26130180 SAN ISIDRO PM ACTIVA 4,80097 -75,85833 1970 43
IDEAM 26130200 EL RECUERDO PM ACTIVA 4,96242 -75,74167 1970 25
IDEAM 26130220 MARACAY PM ACTIVA 4,81217 -75,84306 1978 25
IDEAM 26130540 PLAYA RICA PM ACTIVA 4,75672 -75,59719 1993 20

1-13

13
 Entidad Código Estación Tipo Estado Longitud Latitud Inicio Fin Años
HDA LAS
IDEAM 26130550 4,73333 -75,51667 1993 1998 3
BRISAS PM SUSP.
IDEAM 26130560 PEZ FRESCO PM ACTIVA 4,73272 -75,57556 1993 20
IDEAM 26130570 NUEVO LIBARE PM ACTIVA 4,78206 -75,64583 1999 9
APTO
IDEAM 26135040 4,81586 -75,73722 1947 43
MATECANA SP ACTIVA
UN PEREIRA
IDEAM 26135080 4,80000 -75,70000 1963 1971 2
UTP CO SUSP.
IDEAM 26135100 LA BOHEMIA CO ACTIVA 4,87803 -75,87756 1963 25
IDEAM 26135110 VERACRUZ CO MANT. 4,86764 -75,62781 1977 25
IDEAM 26135140 LA LAGUNA CO ACTIVA 4,77642 -75,41250 1981 19
IDEAM 26135160 SAN REMO PM ACTIVA 4,84167 -75,57997 1981 26
IDEAM 26140180 LA VIRGINIA PM ACTIVA 4,89250 -75,88278 1979 34
IDEAM 26145030 LA CAMELIA CO ACTIVA 5,07936 -75,96706 1963 25
IDEAM 26150150 PAPAYAL PM ACTIVA 4,95167 -75,48972 1970 25
IDEAM 26155110 APTO LA NUBIA SP ACTIVA 5,02978 -75,46992 1968 25
IDEAM 26155150 LAS BRISAS CP ACTIVA 4,93453 -75,35039 1981 25

Fuente: Adaptado de POMCA Río Otún (2017)

1-14

14
Dentro del componente morfológico del POMCA también denominado hidrográfico, fueron
identificados los principales afluentes que conforman las subcuencas del río Otún y se mencionan
a continuación: quebradas Barbo, San Juan, San José, Dos Quebradas, La Combia, Azul, Grande
y Hato Viejo, siendo los dos últimos, afluentes directos al río cauca.

En la siguiente figura se presenta la distribución espacial de las subcuencas de interés del

POMCA del río Otún:

Figura 1-1. Distribución de las subcuencas del río Otún analizadas en el plan de
manejo de la cuenca (POMCA)

Fuente: POMCA Río Otún (CONSORCIO, 2017)

Por otra parte, el estudio de diagnóstico del riesgo adelantado por la corporación regional del
departamento de Risaralda (CARDER, 2014) dentro de su componente climatológico también
realizó un análisis climático a escala de cuenca y de región de variables como la temperatura,
precipitación y la humedad relativa, detallando de manera particular las cuencas de los ríos Otún,
Consota, la Quebrada Cestillal y el río Barbas.

El estudio concluye que el régimen hidrológico característico de la región es bimodal y presenta

afectaciones en los caudales según la fase del fenómeno ENSO que se encuentra activa. La fase

1-15

15
cálida refleja grandes cambios en el periodo enero-marzo, mostrando una disminución en los
caudales medios en este periodo. Por el contrario, la fase fría del fenómeno ENSO se ve reflejada
mayormente en el período diciembre-marzo, traduciéndose en un aumento considerable en los
caudales medios en este periodo de tiempo.

A pesar de la antigüedad y escala de este diagnóstico, algunos resultados obtenidos en el análisis

espacial de las variables climatológicas, en particular de la precipitación, fueron utilizadas como
referencia cualitativa en este estudio.

1.1.1.2 Información primaria

A continuación, se presenta una breve descripción de la información básica tenida en cuenta para
el presente estudio

I. Información hidroclimatológica

Se solicitaron los parámetros de: precipitación diaria total, máxima en 24 horas, evaporación
promedio diaria, temperatura promedio mínima, media y máxima mensual, caudales medios,
mínimos y máximos a nivel diario de la red climatológica de la red nacional de estaciones IDEAM
en la zona de influencia de la cuenca del río Otún.

El listado de las estaciones climatológicas recibidas por parte de la entidad se encuentra a

continuación:

Tabla 1-2 Lista de estaciones suministradas por IDEAM en el 2019 para el estudio
Código Estación Tipo Estado Latitud Longitud Inicio Fin Años
26130020 TERMALES PM ACTIVA 4.84417 -75.55806 1971 2017 47
26130180 SAN ISIDRO PM ACTIVA 4.80097 -75.85833 1970 2017 48
26130190 BOQUERON PM SUSP. 4.83333 -75.66667 1970 1980 11
26130540 PLAYA RICA PM ACTIVA 4.75672 -75.59719 1993 2017 25
26130560 PEZ FRESCO PM ACTIVA 4.73272 -75.57556 1993 2017 25
26135040 APTO MATECANA SO ACTIVA 4.81586 -75.73722 1947 2017 71
26135140 LA LAGUNA CO ACTIVA 4.77642 -75.41250 1994 2017 24
26137240 LA LAGUNA DEL OTUN LM ACTIVA 4.77667 -75.41222 1996 2015 20
26137110 LA BANANERA 609 LM ACTIVA 4.77339 -75.63306 1971 2013 43
26130550 LAS BRISAS HDA PM SUSP. 4.73333 -75.51667 1993 1998 6
26135080 UNIV PEREIRA UTP CO SUSP. 4.80000 -75.70000 1963 1971 9

Fuente: IDEAM

Es importarte anotar que se complementaron estos registros con información de promedios

climatológicos emitidos por IDEAM e información secundaria de estudios previos de la cuenca,
mencionados a continuación.

II. Información cartográfica

 Registro fotográfico extraída de periódicos y videos de registro de inundaciones en el sitio

de estudio.

1-16

16
 Información de Corine Land Cover Colombia (IDEAM, 2007), para definir el número de
curva de la cuenca.
 Información Lidar con batimetría en el tramo de estudio incluida para construcción de la
geometría del modelo hidráulico, suministrada por Conhydra S.A.

1-17

17
1.2 Climatología
A continuación, se presenta el análisis climático de la zona partiendo de una caracterización
morfométrica de la cuenca principal sobre el río Otún y las subcuencas asociadas a los afluentes
principales, un análisis temporal y espacial de variables como la precipitación, temperatura y la
evaporación; caracterizando, de esta manera, el régimen climatológico en la zona de interés.

1.3 Morfometría de las cuencas

La cuenca del río Otún tiene correspondencia con las características físicas presentes en la
Cordillera central, por lo que tiene la peculiaridad de una topografía abrupta. De lo anterior, el río
presenta características típicas a las de una corriente de montaña: alta pendiente y lecho rocoso.
La dirección del río es de Oriente-occidente y tiene confluencia con el río Cauca, atravesando el
casco urbano del municipio de Pereira e influenciando de manera directa las condiciones
hidráulicas de la cabecera.

Dentro de la caracterización morfológica de la zona de interés a nivel de cuenca, se tienen

afluentes como las quebradas: Dosquebradas, San José, San Juan, Barbo, Bosque y Azul. Para
cada una de las cuencas se midió el área, perímetro, longitud del cauce, las cotas superior e
inferior, ancho promedio de la cuenca, pendiente promedia del cauce y con la ayuda de sistemas
de información geográfica se estimaron otros parámetros morfométricos. A continuación, se
presentan las ecuaciones utilizadas para su estimación:

 Pendiente media de la cuenca (P) (%): basada en la relación entre el desnivel máximo
de la cuenca y la longitud del cauce más largo
Cota max  Cota min
P
Lcauce
Ecuación 1.1 Pendiente media de la cueca o cauce principal

 Ancho de la cuenca (W) (km): relación entre el área y la longitud de la cuenca

A
W
L

Ecuación 1.2 Ancho medio de la cuenca

 El Índice de compacidad de Gravelius (Kc): es la relación entre el perímetro de la cuenca

y el perímetro de un círculo de área igual a la de la cuenca.

P
Kc  0.28
A
Ecuación 1.3 Índice de compacidad de Gravelius

 Factor de forma (kf): relaciona el ancho de promedio de la cuenca y la longitud de la

misma

1-18

18
 A
Kf 
L2
Ecuación 1.4 factor o coeficiente de forma de la cuenca

 Indice de alargamiento de la cuenca (Ia): relaciona la longitud máxima de la cuenca y el

ancho perpendicular

Lm
Ia 
L
Ecuación 1.5 Índice de alargamiento de la cuenca

 Densidad de drenaje de las cuencas (km/km 2): es la relación de la sumatoria de todas las
longitudes de los cauces y el área de la cuenca

D
L i

A
Ecuación 1.6 Densidad de drenaje de la cuenca

A continuación, se muestra la distribución espacial de las cuencas de los afluentes principales del
río Otún.

Figura 1-2. Distribución espacial de las cuencas de los afluentes más importantes del
río Otún

1-19

19
Fuente: Hidrased, 2019.

Tabla 1-3 Resumen de parámetros morfométricos de las cuencas de interés

Q. Dosquebradas

Q. San Juan
Q. San José
Q. Bosque
Parámetro

Q. Barbo

Q. Azul
PTAR
Área (km2) 409 61 32 59 27 22 12
Perímetro (km) 189 47 40 48 44 33 24
Cota máxima (msnm) 4937 4479 4715 2248 3156 3581 4714
Cota mínima (msnm) 1124 1837 3208 1289 1489 1676 3195
Longitud del cauce principal (km) 72 19 13 15 17 13 9
Sumatoria drenajes (km) 841 128 63 135 55 35 22
Pendiente (%) 0,1 0,1 0,1 0,1 0,1 0,2 0,2
Ancho (km) 5,7 3,1 2,6 3,8 1,6 1,7 1,4
Factor forma (kf) 0,1 0,2 0,2 0,3 0,1 0,1 0,2
Coeficiente compacidad (kc) 2,6 1,7 2 1,8 2,4 2 2
Alargamiento (ia) 12,8 6,2 4,9 4 11 7,5 6,3
Densidad de drenaje (km/km2) 2,1 2,1 2 2,3 2,1 1,6 1,8

Fuente: Hidrased, 2019.

La cuenca del río Otún a la altura de la PTAR es considerada en términos de área como
“intermedia pequeña” (400<A<600 km2), con un índice de alargamiento alto (Ia>4,9), por lo que
en términos de forma se considera “muy alargada”, y con un coeficiente de compacidad que le
permite ser caracterizada como “ovalada rectangular”; posee una densidad de drenaje que oscila
en 2,1 clasificándose como “bien drenada”, es decir con presencia significativa de drenajes en
toda su extensión.

1.4 Análisis temporal y espacial de las precipitaciones

A continuación, se detalla el análisis de la precipitación de los valores medios y máximos de la
zona y el proceso de elaboración de las curvas de intensidad-frecuencia-duración (IFD) para las
estaciones que representan espacialmente la zona de interés.

1.4.1 Precipitaciones medias

Para el análisis temporal y espacial de las precipitaciones medias de la zona se utilizaron los
registros de las estaciones meteorológicas del IDEAM solicitadas por Hidrased S.A.S., la base de
datos hidrológica del POMCA 2017, y como información complementaria, los promedios
climatológicos del IDEAM entre el período 1980 y 2010.

En la Tabla 1-4 se muestra un resumen de las estaciones de medición utilizadas para el análisis
de las precipitaciones medias de la zona.

1-20

20
 Tabla 1-4 Lista de las estaciones de medición de precipitación total empleadas en el estudio
Estudio Entidad Código Estación Tipo Estado Inicio Fin Años
Hidrased 2019 IDEAM 26130020 TERMALES PM ACTIVA 1971 2017 47
Hidrased 2019 IDEAM 26130180 SAN ISIDRO PM ACTIVA 1970 2017 48
Hidrased 2019 IDEAM 26130190 BOQUERON PM SUSPENDIDA 1970 1980 11
Hidrased 2019 IDEAM 26130540 PLAYA RICA PM ACTIVA 1993 2017 25
Hidrased 2019 IDEAM 26130560 PEZ FRESCO PM ACTIVA 1993 2017 25
Hidrased 2019 IDEAM 26135040 APTO MATECANA SO ACTIVA 1947 2017 71
Hidrased 2019 IDEAM 26135140 LA LAGUNA CO ACTIVA 1994 2017 24
POMCA 2017 CHEC 26120310 BREMEN PG SUSPENDIDA 1977 1999 23
POMCA 2017 CHEC 26130060 SUB ESTA ROSA PM ACTIVA 1978 1995 18
POMCA 2017 CHEC 26130070 SAN RAMON PM ACTIVA 1968 1991 24
POMCA 2017 CHEC 26130120 SUB ESTA ROSA PM ACTIVA 1978 1995 18
POMCA 2017 CHEC 26130130 LA FLORIDA PM ACTIVA 1978 1995 18
POMCA 2017 FNC 26130210 LOS CAMBULOS PM ACTIVA 1964 1996 33
POMCA 2017 FNC 26130230 P TRATAMIENTO PM ACTIVA 1969 1997 29
POMCA 2017 FNC 26130260 EL BOSQUE PM ACTIVA 1978 1996 19
POMCA 2017 FNC 26130510 COMBIA PM ACTIVA 1981 1996 16
POMCA 2017 FNC 26135060 EL JAZMIN CO ACTIVA 1960 1997 38
POMCA 2017 FNC 26135070 EL CEDRAL CO ACTIVA 1972 2002 31
POMCA 2017 IDEAM 21210180 TOCHE PG ACTIVA 1986 2013 28
POMCA 2017 IDEAM 26100830 PALMASOLA PM ACTIVA 1981 2013 33
POMCA 2017 IDEAM 26110150 LA INMACULADA PM MANT. 1969 2013 45
POMCA 2017 IDEAM 26115040 LA UNION CP ACTIVA 1988 2012 25
POMCA 2017 IDEAM 26120150 ALCALÁ PM ACTIVA 1972 2011 40
POMCA 2017 IDEAM 26120160 SALENTO PM ACTIVA 1974 2013 40
POMCA 2017 IDEAM 26130170 POTREROS PM ACTIVA 1970 2013 44
POMCA 2017 IDEAM 26130200 EL RECUERDO PM ACTIVA 1988 2013 26
POMCA 2017 IDEAM 26130220 MARACAY PM ACTIVA 1978 2013 36
POMCA 2017 IDEAM 26135100 LA BOHEMIA CO ACTIVA 1988 2013 26
POMCA 2017 IDEAM 26135110 VERACRUZ CO MANT. 1988 2013 26
POMCA 2017 IDEAM 26135160 SAN REMO PM ACTIVA 1987 2013 27

1-21

21
 Estudio Entidad Código Estación Tipo Estado Inicio Fin Años
POMCA 2017 IDEAM 26140180 LA VIRGINIA PM ACTIVA 1979 2013 35
POMCA 2017 IDEAM 26145030 LA CAMELIA CO ACTIVA 1988 2013 26
POMCA 2017 IDEAM 26150150 PAPAYAL PG ACTIVA 1988 2013 26
POMCA 2017 IDEAM 26155110 APTO LA NUBIA SP ACTIVA 1988 2013 26
POMCA 2017 IDEAM 26155150 LAS BRISAS SP ACTIVA 1988 2013 26
IDEAM PROM IDEAM 26125060 APTO EL EDEN SP ND 1980 2010 ND
IDEAM PROM IDEAM 26125070 BELLA LA CO ND 1980 2010 ND
IDEAM PROM IDEAM 26140200 TRIBUNA LA PM ND 1980 2010 ND
IDEAM PROM IDEAM 21240070 ANZOATEGUI PM ND 1980 2010 ND
IDEAM PROM IDEAM 21210160 DARIEN EL PG ND 1980 2010 ND
IDEAM PROM IDEAM 21210120 ESMERALDA LA PG ND 1980 2010 ND
IDEAM PROM IDEAM 21210020 JUNTAS LAS PG ND 1980 2010 ND
IDEAM PROM IDEAM 21210220 PALMAR EL PG ND 1980 2010 ND
IDEAM PROM IDEAM 21210030 PASTALES PG ND 1980 2010 ND
IDEAM PROM IDEAM 21210110 PLACER EL PG ND 1980 2010 ND
IDEAM PROM IDEAM 21210100 RANCHO EL PG ND 1980 2010 ND
IDEAM PROM IDEAM 21240030 SAN JUAN DE CHINA PM ND 1980 2010 ND
IDEAM PROM IDEAM 21210080 SECRETO EL PG ND 1980 2010 ND
IDEAM PROM IDEAM 21250500 LIBANO PM ND 1980 2010 ND
IDEAM PROM IDEAM 21250050 ALTO DEL OSO PM ND 1980 2010 ND
IDEAM PROM IDEAM 21250110 MURILLO PM ND 1980 2010 ND
IDEAM PROM IDEAM 21245090 ERMITA LA ME ND 1980 2010 ND
IDEAM PROM IDEAM 21245140 STA ISABEL ME ND 1980 2010 ND
IDEAM PROM IDEAM 21250170 BOQUERON PM ND 1980 2010 ND
IDEAM PROM IDEAM 26120600 PIEDRAS DE MOLER PM ND 1980 2010 ND
IDEAM PROM IDEAM 26100330 SOCORRO EL PG ND 1980 2010 ND
IDEAM PROM IDEAM 26100290 CRUCES PM ND 1980 2010 ND
IDEAM PROM IDEAM 26100300 OBANDO PM ND 1980 2010 ND
IDEAM PROM IDEAM 26100790 ITALIA LA PM ND 1980 2010 ND
IDEAM PROM IDEAM 26100350 PTO FRAZADAS PM ND 1980 2010 ND

Fuente: Hidrased, 2019. con base en catálogo nacional de estaciones IDEAM 2019.

1-22

22
1.4.1.2 Detección datos anómalos

Con el fin de analizar la consistencia de los registros de precipitación de las siete estaciones solicitados al
IDEAM para el estudio, se utilizaron las pruebas de datos anómalos de Grubbs (1969) y de Dixon (1983).
Los resultados obtenidos para algunas de las estaciones se encuentran desde en las figuras mostradas a
continuación y el resumen de los mismos en Tabla 1-5. Es de resaltar que los círculos rojos representan
los datos anómalos y las líneas azules la serie de datos.

140 140
120 120
100 100
Precipitación (mm)

Precipitación (mm)
80 80
60 60
40 40
20 20
0 0
1993 2002 2011 1993 2002 2011

Figura 1-3. Análisis de anómalos por metodología de Grubbs (izq.) y Dixon (der.) para la
estación Pez fresco

140 140
120 120
100 100
Precipitación (mm)

Precipitación (mm)

80 80
60 60
40 40
20 20
0 0
1971 1980 1989 1998 2007 2016 1971 1980 1989 1998 2007 2016

Figura 1-4. Análisis de anómalos por metodología de Grubbs (izq.) y Dixon (der.) para la
estación Termales

1-23

1-23
 140 140
120 120
100 100
Precipitación (mm)

Precipitación (mm)
80 80
60 60
40 40
20 20
0 0
1947 1956 1965 1974 1983 1992 2001 2010 1947 1956 1965 1974 1983 1992 2001 2010

Figura 1-5. Análisis de anómalos por metodología de Grubbs (izq.) y Dixon (der.) para la
estación Apto Matecaña

Fuente: Hidrased, 2019.

Tabla 1-5 Análisis de anómalos para las estaciones utilizando la metodología de DIXON y
GRUBBS
Número de anómalos
Código Nombre
DIXON GRUBBS
26130020 TERMALES 810 809
26130180 SAN ISIDRO 931 931
26130190 BOQUERON 108 195
26130540 PLAYA RICA 472 472
26130560 PEZ FRESCO 458 458
26135040 APTO MATECANA 1449 1449
26135140 LAGUNA LA 325 325

Fuente: Hidrased, 2019.

A pesar de que las pruebas mencionadas encontraron para cada estación numerosos datos anómalos
desde el punto de vista estadístico, luego de analizar individualmente cada uno de los registros, y teniendo
en cuenta la alta variabilidad de las series de tiempo de las precipitaciones de la zona, ninguno de ellos
fue descartado para el análisis. En general, los valores marcados corresponden a eventos aislados o
“micro-picos” dentro de la serie, que incluso se encuentran dentro de los valores esperados de la
precipitación acumulada de la zona.

De la misma forma, no se consideró necesario rellenar la información de los años faltantes en las series
debido a que gran parte de éstas tienen entre 15 y 45 años de información.

1.4.1.3 Análisis de homogeneidad y tendencia

Las series de Pettit (1979) y Von Neumman (Kang & Yusof, 2012) fueron empleadas para el análisis de
homogeneidad, mientras que para el análisis tendencia se utilizó la prueba de Mann-Kendall (Mann, 1945)
(Kendall, 1955) y el estimativo de la Pendiente de Sen. A continuación, se presentan los resultados
obtenidos de ambos análisis.

1-24

1-24
 Tabla 1-6 Análisis de homogeneidad y cambio en la media
Von Neumann Pettit
Código Nombre
Cambio N Cambio Fecha media1 media2
26130020 TERMALES Si 1,51 Si 11/05/1990 9,42 7,26
26130180 SAN ISIDRO Si 1,81 Si 21/11/1975 6,84 4,49
26130190 BOQUERON Si 1,97 Si 25/12/1975 12,73 9,41
26130540 PLAYA RICA Si 1,61 Si 21/09/2005 7,33 8,65
26130560 PEZ FRESCO Si 1,54 Si 24/07/1998 6,65 8,36
26135040 APTO MATECANA Si 1,77 Si 21/03/1998 5,94 7,09
26135140 LAGUNA LA Si 1,60 Si 06/07/2002 3,29 2,92

Fuente: Hidrased, 2019.

Tabla 1-7 Análisis de tendencia de las estaciones

Mann-kendall
Código Nombre
Tendencia Tau P Sen
26130020 TERMALES SI -0,064 0,003 0,00000
26130180 SAN ISIDRO No -0,035 1,000 0,00011
26130190 BOQUERON No -0,037 0,056 0,00000
26130540 PLAYA RICA No 0,006 0,385 0,00029
26130560 PEZ FRESCO No 0,009 1,000 0,00019
26135040 APTO MATECANA No 0,051 1,000 0,00000
26135140 LAGUNA LA SI -0,044 0,004 0,00000

Fuente: Hidrased, 2019.

El análisis de la homogeneidad de las estaciones muestra que existe un leve cambio de media en algunas
de las estaciones con significancia estadística, presente en los primeros 10 años de registro de algunas
estaciones. Este cambio en la media puede ser debido a posibles ajustes en los equipos de medición o
procesamiento de la información reportada por el IDEAM.

Por otro lado, los resultados de la prueba de Mann-Kendall de las estaciones y el estimador de Sen
muestran que la mayoría de las estaciones no tienen una tendencia con significancia estadística y las
estaciones que si se presume de tendencia (Termales y La Laguna), la pendiente de Sen estimada es muy
cercana a 0, desestimando lo mostrado por la prueba. En general ninguna de las estaciones presenta una
tendencia clara en sus datos.

La variabilidad del total de la precipitación mensual promedio en las estaciones más representativas de la
zona, tales como: Termales, San Isidro, Boquerón, Playa Rica, Pez Fresco, La Laguna y Aeropuerto
Matecaña se muestra a continuación.

1-25

1-25
 450
400
350
Precipitación total (mm)

300
250
200
150
100
50
0
ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC
Mes
TERMALES SAN ISIDRO BOQUERON PLAYA RICA
PEZ FRESCO APTO MATECANA LA LAGUNA

Figura 1-6. Promedio de la precipitación total mensual de las estaciones representativas de la

cuenca del río Otún

Fuente: Hidrased, 2019.

En la anterior figura es posible observar un comportamiento bimodal de la precipitación que caracteriza a

esta zona, donde se evidencia un aumento en la magnitud promedio en el periodo marzo-mayo y
septiembre-noviembre separados por periodos de menos precipitación. Esta variación corresponde con el
paso de la zona de convergencia intertropical (ZCIT), definida como un cinturón de bajas presiones que
recorre latitudinalmente la zona tropical y que en la zona Andina colombiana influye en un aumento tanto
en la magnitud como en la frecuencia de las lluvias.

En particular, los meses de octubre y noviembre registran los mayores promedios de precipitación,
mientras julio y agosto los menores. La estación Boquerón presenta para el mes agosto un promedio
atípico respecto al resto de las estaciones.

Las precipitaciones promedio de la región de interés tienen una variabilidad que va desde los 80 mm en el
mes de julio hasta los 370 mm en los meses de octubre y noviembre, con un promedio mensual multianual
cercano a los 190 mm y un promedio anual total cercano a los 2400 mm.

Con el fin de caracterizar espacialmente la precipitación total anual de la zona, se utilizaron todos los
registros promedios multianuales de las estaciones y mediante la interpolación geo-espacial, utilizando el
estimador de Kriging Ordinario (KO), se interpoló un campo de precipitación. A continuación, se muestra
en la Figura 1-7 la variación espacial de la precipitación promedio multianual de la zona.

1-26

1-26
 Figura 1-7. Campo de precipitación promedio multianual sobre la cuenca del río Otún

Fuente: Hidrased, 2019.

En la anterior figura es posible observar que existen al menos tres regiones espacialmente marcadas en
la distribución de la precipitación total anual a lo largo de la cuenca del río Otún, la cual se encuentra
orientada en sentido oriente-occidente. La zona de la cuenca alta (1) presenta las menores precipitaciones
con variaciones entre 1000-1400 mm/año, por otra parte, la cuenca media (2) presenta las mayores
precipitaciones con una variación desde 2200 hasta 3000 mm/año y la cuenca baja (3), donde se encuentra
ubicada la zona aledaña al lugar de la construcción de la PTAR, presenta variaciones entre 1700-2200
mm/año. En la margen izquierda del río Otún (zona sur) se presentan menores magnitudes de precipitación
que en la margen derecha (zona norte), presentándose variaciones promedio de hasta 400 mm/año entre
ellas.

En el mapa M-PTAR-PER-HHS-CLIMA-01 anexo a este estudio es posible encontrar, con más detalle, el
resultado del campo de precipitación y sus correspondientes isoyetas calculadas, así como la ubicación
espacial de las estaciones empleadas.

1.4.2 Precipitaciones máximas

Para el análisis de las lluvias extremas se utilizaron los registros de precipitaciones máximas en 24 horas
de las estaciones relevantes a la zona de estudio: Termales, Boquerón, Playa Rica, Pez Fresco,
Aeropuerto Matecaña.

1-27

1-27
El análisis de frecuencia y los ajustes con las funciones de distribución de probabilidad Gumbel, Log
Gumbel, Pearson, Log Pearson, Normal y Log Normal, fueron utilizados para estimar los valores totales
acumulados de las máximas precipitaciones en 24 horas con períodos de retorno desde 2,33 hasta 100
años para las estaciones mencionadas con anterioridad.

Los valores asociados a las lluvias máximas en24 horas para distintos periodos de retorno fueron
determinados por medio de las pruebas estadísticas Chi Cuadrado y Smirnov - Kolmogorov, con las cuales
es posible estimar la bondad del ajuste de las distribuciones de probabilidad utilizadas sobre los datos de
cada estación, y los criterios espaciales de la distribución de la lluvia. Los resultados se presentan en la
Tabla 1-8.

Tabla 1-8 Precipitación máxima en 24 horas en distintos períodos de retorno para las
estaciones solicitadas
Estación Termales Boquerón Playa rica Pez fresco Aero Matecaña La Laguna
Tr (años) Gumbel Pearson Gumbel Pearson Pearson Log Pearson
1,01 46 11 56 45 48 17
2,33 79 101 86 96 80 34
5 94 115 99 107 91 41
10 106 122 109 114 98 46
20 117 127 120 120 105 51
25 121 128 123 121 107 53
30 124 129 125 123 109 54
50 132 132 133 126 113 57
100 143 134 143 130 119 61
500 169 138 166 138 131 71
1000 180 139 176 140 136 75

Fuente: Hidrased, 2019.

1.4.3 Elaboración de curvas Intensidad Frecuencia Duración

Uno de los principales objetivos del presente estudio consiste en determinar las curvas de Intensidad
Frecuencia Duración (IFD) de cada una de las estaciones con influencia en la zona, y realizar un análisis
espacial, es decir, definir el área de influencia y aplicabilidad de cada curva.

Se realizó un análisis de las IFD estimadas en el POMCA para las mismas estaciones y se comparó con
otras metodologías ampliamente usadas en el medio tales como: (INVIAS, 2009) , (Pulgarin, 2009) ,
(Hershfield, 1961). Como punto de observación para esta comparación se tomó la estación del aeropuerto
Matecaña, debido a su longitud de registro que superan los 40 años y la confiablidad que se tiene en la
calidad de los registros de esta estación. A continuación, se mencionan las diferentes metodologías
utilizadas para la elaboración de las curvas Duración – Intensidad – Frecuencia (DIF)

Hershfield (1961)

𝑃𝐷 = 𝑃𝑇𝑟 ∗ (𝑎 + 𝑏 𝐿𝑛(𝐷))
Ecuación 1.7 Precipitación para cada duración según Hershfield

1-28

1-28
Donde,

PD: precipitación para cada duración (mm)

PTr: precipitación asociada a cada periodo de retorno (mm)

a: coeficiente equivalente a -0,1

b: coeficiente equivalente a 0,17

D: curación (horas)

Manual INVIAS (2009)

𝑎 × 𝑇 𝑏 × 𝑀𝑑
𝑖=
𝑡 𝑐
( )
60
Ecuación 1.8 Intensidad para cada duración según INVIAS

Donde,

i: intensidad de precipitación, en milímetros por hora (mm/h).

T: periodo de retorno en años.

M: precipitación máxima promedio anual en 24 h a nivel multianual

t: duración de la lluvia, en minutos (min).

a, b, c, d: parámetros de ajuste de la regresión y sus valores se presentan en la siguiente tabla

Tabla 1-9 Valores de los coeficientes a, b, c y d para el cálculo de las curvas intensidad-
duración-frecuencia, IDF, para Colombia
Región A B C D
Andina 0,94 0,18 0,66 0,83
Caribe 24,85 0,22 0,5 0,1
Pacífico 13,92 0,19 0,58 0,2
Orinoquía 5,53 0,17 0,63 0,42

Fuente: (INVIAS, 2009)

Pulgarín (2009)

1-29

1-29
 0,12 𝑚 𝑡 𝜃
𝑖 = ([0,88 𝑚 − 0,004] + [ × 𝑦]) ( )
0,5772 1440

Ecuación 1.9 Intensidad para cada duración según Pulgarín

Donde,

i: intensidad para un periodo de retorno T y una duración t, expresada en mm/hora

t: duración de la tormenta en minutos

T: periodo de retorno en años

M: intensidad promedia diaria anual (mm/hora)

Q: exponente de escalamiento -0.829

Y: −ln[− ln(1− (1/T ))] es la transformación del periodo de retorno T

Si se desea obtener valores de precipitación o de intensidad para rangos menores de una hora, se sugiere
encontrar el valor de intensidad de 1 hora mediante la ecuación anterior y transformarla en precipitación
de 1 hora y usar los estimativos empíricos de Bell descritos en la Ecuación 1.10.

𝑃𝑇𝑑
= 0,54 𝑑 0,25 − 0,5 𝑝𝑎𝑟𝑎 5 ≤ 𝑑 ≤ 120
𝑃𝑇60
Ecuación 1.10 Escalamiento con base en Bell

Donde,

PTd: precipitación máxima en un intervalo de tiempo d (en minutos), con un

período de retorno de T años.

PT60: es la precipitación máxima en un intervalo de tiempo de 60 minutos con un período de retorno de T

años.

Aplicando cada una de las metodologías anteriormente expuestas, se estimó la curva IDF de la estación
Matecaña para 100 años de periodo de retorno. Donde se observa que la curva presentada en el POMCA
tiende a sobrestimar las intensidades en duraciones cortas, mientras que la metodología de Hershfiel y de
INVIAS son muy similares para todas las duraciones. La curva generada con Pulgarín (Pulgarin, 2009)
para duraciones mayores a 90 min es muy similar a la presentada en el POMCA, pero en duraciones
menores a 30 min es menor la intensidad generada por esta. Con base en lo anterior se seleccionó la
metodología de Pulgarín para estimar la IDF del presente estudio, debido a que es conservadora con

1-30

1-30
respecto a Hershfiel y INVIAS, pero no sobrestima los caudales para pequeñas duraciones comparándola
con la curva presentada en el POMCA de 2017.

400
Hershfiel (1961)
350 POMCA -Vargas y Díaz Granados (1998)
300 Pulgarin (2009)
Intensidad (mm/h)

Manual Invias (2011)

250

200

150

100

50

0
0 20 40 60 80 100 120 140 160 180
Duración (min)

Figura 1-8. Curvas IDF para 100 años de periodo de retorno por diferentes metodologías

Fuente: Hidrased, 2019.

A continuación, se presenta las curvas IDF generadas para cada una de las estaciones de análisis y en la
Figura 1-9 se muestra el área de influencia de éstas, es importante mencionar que, el área aferente a cada
estación se estimó mediante los polígonos de Thiessen y herramientas SIG.

1-31

1-31
 Figura 1-9. Polígonos de Thiessen

Fuente: Hidrased, 2019.

1.4.3.2 Curvas IFD estación Termales

En la Tabla 1-10 se presentan los valores de intensidad para diferentes duraciones y periodos de retorno
de la estación Temarles y en la Figura 1-10 se muestra de manera gráfica.

Tabla 1-10 Intensidad Frecuencia Duración de la estación Termales

Precipitación Tr1 Intensidad (mm/hora)
diaria (años) 5 8 10 12 15 20 25 30 45 60 90 120 180 360
79 2,33 111 100 93 88 82 73 67 62 52 46 33 26 18 10
94 5 156 141 132 125 115 104 95 88 74 65 46 37 26 15
106 10 200 180 168 159 147 132 121 112 95 83 59 47 33 19
121 25 262 235 221 208 193 173 159 147 124 109 78 61 44 25
132 50 313 281 264 249 230 207 190 176 148 130 93 73 52 29
143 100 368 331 311 293 271 244 223 207 174 153 109 86 61 35

Fuente: Hidrased, 2019.

1 Tr: Periodo de retorno

1-32

1-32
 400
2.33 años
350 5 años
10 años
300 25 años
intensidad (mm)

50 años
250
100 años

200

150

100

50

0
0 20 40 60 80 100 120 140 160 180
Duración (min)
Figura 1-10. Curvas Intensidad Frecuencia Duración de la estación Termales

Fuente: Hidrased, 2019.

1.4.3.3 Curvas IFD estación Boquerón

En la Tabla 1-11 se presentan los valores de intensidad para diferentes duraciones y periodos de retorno
de la estación Boquerón y en la Figura 1-11 se muestran gráficamente.

Tabla 1-11 Intensidad Frecuencia Duración de la estación Boquerón

Precipitación Intensidad(mm/hora)
Tr (Años)
diaria 5 8 10 12 15 20 25 30 45 60 90 120 180 360
79 2,33 141 127 119 112 104 93 86 79 67 59 42 33 24 13
94 5 191 172 161 152 141 127 116 108 91 79 57 45 32 18
106 10 230 207 194 183 169 152 140 130 109 95 68 54 38 22
121 25 277 249 234 221 204 184 168 156 131 115 82 65 46 26
132 50 312 280 263 248 230 206 189 175 147 129 92 73 52 29
143 100 345 310 291 274 254 228 209 194 163 143 102 81 58 32

Fuente: Hidrased, 2019.

1-33

1-33
 400
2.33 años
350 5 años
10 años
300 25 años
intensidad (mm)

50 años
250
100 años

200

150

100

50

0
0 20 40 60 80 100 120 140 160 180
Duración (min)
Figura 1-11. Curvas Intensidad Frecuencia Duración de la estación Boquerón

Fuente: Hidrased, 2019.

1.4.3.4 Curvas IFD estación Playa Rica

En la Tabla 1-12 se presentan los valores de intensidad para diferentes duraciones y periodos de retorno
de la estación Playa Rica y en la Figura 1-12 se muestran gráficamente.

Tabla 1-12 Intensidad Frecuencia Duración de la estación Playa Rica

Precipitación Intensidad(mm/hora)
Tr (Años)
diaria 5 8 10 12 15 20 25 30 45 60 90 120 180 360
79 2,33 120 108 101 95 88 79 73 67 57 50 35 28 20 11
94 5 164 148 139 131 121 109 100 93 78 68 49 38 27 15
106 10 206 185 174 164 152 136 125 116 98 86 61 48 34 19
121 25 265 239 224 211 196 176 161 149 126 110 79 62 44 25
132 50 314 283 265 250 231 208 191 177 149 130 93 73 52 30
143 100 367 330 309 292 270 243 222 206 174 152 109 86 61 34

Fuente: Hidrased, 2019.

1-34

1-34
 400
2.33 años
350 5 años
10 años
300 25 años
intensidad (mm)

50 años
250
100 años

200

150

100

50

0
0 20 40 60 80 100 120 140 160 180
Duración (min)
Figura 1-12. Curvas Intensidad Frecuencia Duración de la estación Playa Rica

Fuente: Hidrased, 2019.

1.4.3.5 Curvas IFD estación Pez Fresco

En la Tabla 1-13 se presentan los valores de intensidad para diferentes duraciones y periodos de retorno
de la estación Pez Fresco y en la Figura 1-13 se muestran gráficamente.

Tabla 1-13 Intensidad Frecuencia Duración de la estación Pez Fresco

Precipitación Intensidad(mm/hora)
Tr (Años)
diaria 5 8 10 12 15 20 25 30 45 60 90 120 180 360
79 2,33 134 120 113 106 98 89 81 75 63 56 40 31 22 13
94 5 178 160 150 142 131 118 108 100 84 74 53 42 30 17
106 10 215 194 181 171 159 143 131 121 102 89 64 50 36 20
121 25 262 236 221 209 193 174 159 148 124 109 78 61 44 25
132 50 298 268 251 237 220 197 181 168 141 124 88 70 50 28
143 100 334 300 281 265 246 221 202 188 158 138 99 78 56 31

Fuente: Hidrased, 2019.

1-35

1-35
 400
2.33 años
350 5 años
10 años
300 25 años
intensidad (mm)

50 años
250
100 años

200

150

100

50

0
0 20 40 60 80 100 120 140 160 180
Duración (min)
Figura 1-13. Curvas Intensidad Frecuencia Duración de la estación Pez Fresco

Fuente: Hidrased, 2019.

1.4.3.6 Curvas IFD estación Aeropuerto Matecaña

En la Tabla 1-14 se presentan los valores de intensidad para diferentes duraciones y periodos de retorno
de la estación Aeropuerto Matecaña y en la Figura 1-14 se muestran gráficamente.

Tabla 1-14 Intensidad Frecuencia Duración de la estación Aeropuerto Matecaña

Precipitación Tr Intensidad(mm/hora)
diaria (Años) 5 8 10 12 15 20 25 30 45 60 90 120 180 360
79 2,33 112 100 94 89 82 74 68 63 53 46 33 26 19 11
94 5 151 136 127 120 111 100 92 85 72 63 45 35 25 14
106 10 185 167 156 148 137 123 113 104 88 77 55 43 31 17
121 25 232 208 195 184 171 153 141 130 110 96 69 54 39 22
132 50 268 241 226 213 197 177 163 151 127 111 79 63 45 25
143 100 305 275 258 243 225 202 185 172 145 127 91 71 51 29

Fuente: Hidrased, 2019.

1-36

1-36
 350
2.33 años
300 5 años
10 años
250 25 años
intensidad (mm)

50 años
100 años
200

150

100

50

0
0 20 40 60 80 100 120 140 160 180
Duración (min)
Figura 1-14. Curvas Intensidad Frecuencia Duración de la estación Aeropuerto Matecaña

Fuente: Hidrased, 2019.

1.4.3.7 Curvas IFD estación La Laguna Otún

En la Tabla 1-15 se presentan los valores de intensidad para diferentes duraciones y periodos de retorno
de la estación La Laguna Otún y en la Figura 1-15 se muestran gráficamente.

Tabla 1-15 Intensidad Frecuencia Duración de la estación La Laguna Otún

Precipitación Intensidad(mm/hora)
Tr (Años)
diaria 5 8 10 12 15 20 25 30 45 60 90 120 180 360
79 2,33 48 43 40 38 35 32 29 27 23 20 14 11 8 5
94 5 68 62 58 54 50 45 42 39 32 28 20 16 11 6
106 10 87 79 74 69 64 58 53 49 41 36 26 20 15 8
121 25 114 102 96 90 84 75 69 64 54 47 34 27 19 11
132 50 135 121 114 107 99 89 82 76 64 56 40 32 23 13
143 100 158 142 133 125 116 104 96 89 75 65 47 37 26 15

Fuente: Hidrased, 2019.

1-37

1-37
 180
2.33 años
160
5 años
140 10 años
25 años
intensidad (mm)

120 50 años
100 años
100

80

60

40

20

0
0 20 40 60 80 100 120 140 160 180
Duración (min)
Figura 1-15. Curvas Intensidad Frecuencia Duración de la estación La Laguna Otún

Fuente: Hidrased, 2019.

1.5 Análisis temporal y espacial de la Temperatura

Para el análisis de la temperatura de la zona se utilizó de forma representativa los registros de la estación
Aeropuerto Matecaña (26135040) ubicada aproximadamente en la cota 1342 msnm, la cual cuenta con
cerca de 70 años de información de este parámetro (1948-2017) de forma consistente y sin vacíos de
información.

A continuación, en la Figura 1-16 y de forma tabular en la Tabla 1-16; se muestran tanto los promedios
mensuales de la estación, así como el promedio de los mínimos y el promedio de los máximos.

1-38

1-38
 30
temperatura promedio (°c)

25

20

15

10

0
ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC
Mes

Mínimo Media Máximo

Figura 1-16. Variación mensual de la temperatura promedio mínima, media y máxima

Fuente: Hidrased, 2019.

Tabla 1-16 Magnitud de las temperaturas promedios en la estación Aero. Matecaña [26135040]
Temperatura Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic Prom.
Mínimo 15,1 15,2 15,3 15,2 15,2 14,9 14,9 14,9 14,8 14,6 14,9 14,9 15,0
Media 22,2 22,3 22,2 21,8 21,6 21,7 22,0 22,2 21,8 21,2 21,2 21,7 21,8
Máximo 27,6 27,8 27,6 27,0 26,6 26,6 27,1 27,5 27,2 26,4 26,3 26,7 27,0

Fuente: Hidrased, 2019. basada en los datos IDEAM

Las temperaturas medias promedio de cada mes presentan variaciones entre 21.2 °C (octubre-noviembre)
y 22.3 °C (febrero), con una diferencia cercana a 1.1°c. La máxima variabilidad entre la temperatura
promedio mínima y máxima se detectó para los meses de febrero y agosto, alcanzando valores de 12,7 y
de 12,6 °C respectivamente. En general, la temperatura media promedio multianual de la estación de 21.8
°C se considera como referente de la zona de interés.

Debido a la estrecha relación que existe entre la temperatura promedio y la elevación, es posible estimar
espacialmente esta variable con una alta confiabilidad, utilizando como pivote una estación con registros
consistentes en el tiempo.

De esta manera, siguiendo la metodología propuesta por Barón (2003) para el IDEAM, donde se
regionalizó los gradientes de temperatura por cada 100 metros utilizando más 1200 estaciones de registro
climático (Tabla 1-17) para toda Colombia, es posible extraer una relación lineal de la temperatura
promedio multianual de una zona utilizando un modelo de elevación digital (DEM).

1-39

1-39
 Tabla 1-17 Gradientes de variación de temperatura por cada 100 m para la región andina
Mes Gradiente Coeficiente de ajuste (R2)
Enero 0,66 0,973
Febrero 0,66 0,968
Marzo 0,64 0,967
Abril 0,61 0,971
Mayo 0,61 0,974
Junio 0,65 0,977
Julio 0,67 0,971
Agosto 0,69 0,970
Septiembre 0,67 0,972
Octubre 0,62 0,973
Noviembre 0,60 0,977
Diciembre 0,63 0,977

Fuente: adaptado de Barón (2003)

Con la temperatura promedio multianual en función de la elevación obtenida de un modelo digital del
terreno y los gradientes de la tabla anterior, se obtiene la distribución espacial de este parámetro. Tal y
como se muestra en la Figura 1-17

1-40

1-40
 Figura 1-17. Variación espacial de la temperatura promedio multianual

Fuente: Hidrased, 2019.

La cuenca del río Otún presenta variaciones de la temperatura promedio multianual desde los 0°C (o
inferiores) hasta los 24°C en la parte baja de la misma. Es posible encontrar con mayor detalle esta
variación en el M-PTAR-PER-HHS-CLIMA-03 a este informe.

1.6 Análisis temporal y espacial de la Evaporación

El volumen de evapotranspiración en la cuenca se calculó con los métodos de CENICAFÉ (1989), Budyko
(1974), Turc (1954) , Countagne (Rojo, 2000) y Thornthwaite (1948)

Para aplicar el método de CENICAFÉ (1989) es necesario estimar inicialmente la evapotranspiración

potencial, que corresponde a la máxima evapotranspiración posible, considerando que la única condición
limitante es la disponibilidad de energía. La evapotranspiración potencial multianual se evaluó (utilizando
herramientas de algebra de mapas mediante el Sistema de Información Geográfica Arc GIS) en función
de la altitud de las cuencas, aplicando la siguiente expresión:

ETP = 365× 4,658× e -0,0002h

Ecuación 1.11 Evapotranspiración potencial según CENICAFÉ

1-41

1-41
Donde,

ETP: evapotranspiración potencial (mm/año)

h: elevación en cada punto de la cuenca sobre el nivel del mar (msnm)

Posteriormente se evaluó la evapotranspiración real que considera la disponibilidad energética, hídrica y

condiciones aerodinámicas. Para este fin se utilizó la ecuación de Budyco (1974), obtenida por medio de
balance de masas:

0 ,5
  P    ETP   ETP   
ETR=  ETP×P× tanh   ×1- cosh  +senh    
  ETP     P   P  

Ecuación 1.12 Evapotranspiración real según Budyco

Donde,

ETR: evapotranspiración real (mm/año)

ETP: evapotranspiración potencial (mm/año)

P: precipitación media anual (mm)

La ecuación de Turc (1954) permite estimar la evapotranspiración anual en mm en función de la

precipitación y de la temperatura de la cuenca mediante la siguiente expresión:

P P
ETR= ,para  0,316
P2 L
( 0,9+ )
L2

P
ETR=P,para  0,316
L
Ecuación 1.13 Evapotranspiración real de Turc

1-42

1-42
Donde,

ETR: evapotranspiración real media anual (mm)

P: precipitación media anual (mm)

L: está una función de la temperatura (t), según la siguiente expresión:

L = 300 + 25T + 0.05T 3

Ecuación 1.14 Coeficiente L de Turc

La metodología de Coutagne (Rojo, 2000)se basa en las condiciones meteorológicas de la cuenca. Es una
estimación de evapotranspiración real en función de la precipitación y de la temperatura como se muestra
a continuación:

𝐸𝑇𝑅 = 𝑃 − 𝜆𝑃2
Ecuación 1.15 ETR según Coutagne

1
𝜆= 𝑠𝑖 (8𝜆)−1 < 𝑃 < (2𝜆)−1
0,8 + 0,14𝑇

Si P < (8𝜆)−1 𝑒𝑛𝑡𝑜𝑛𝑐𝑒𝑠 ETR = P

Si P > (2𝜆)−1 𝑒𝑛𝑡𝑜𝑛𝑐𝑒𝑠 ETR = (4𝜆)−1 = 0,2 + 0,035𝑇

Ecuación 1.16 coeficiente de Coutagne

Donde,

ETR: evapotranspiración real en mm/año

P: precipitación en m

T: temperatura media anual

La ecuación de Thornthwaite (1948) es bastante antigua, y se basa en numerosos experimentos con

lisímetros y de una estimación de evapotranspiración potencial en función de la temperatura media
mensual como se muestra a continuación:

𝑇 𝑎
𝐸𝑇𝑃 = 1,6 (10 )
𝐼
Ecuación 1.17 Evapotranspiración de Thornthwaite

Donde,

1-43

1-43
ETP: Evapotranspiración Potencial (cm/mes)

T: temperatura media mensual en °C

I: índice calórico anual que se define como:

𝑇𝑎𝑛𝑢𝑎𝑙 1,514
𝐼 = 12 [( ) ]
5
Ecuación 1.18 Índice calórico de Thornthwaite

𝑎 = (675 × 10−9 )𝐼 3 − (771 × 10−7 )𝐼 2 + (179 × 10−4 )𝐼 + 0,492

Ecuación 1.19 Alfa de la metodología de Thornthwaite

De esta manera, con los resultados obtenidos de las precipitaciones medias multianuales (numeral 1.4.1)
y la temperatura promedio multianual (numeral 1.5) fueron estimados de forma espacial los valores de la
evaporación potencial (ETP) y evaportranspiración o evaporación real (ETR).

La estación aeropuerto Matecaña, como se ha mencioado anteriormente para otras variables climáticas,
posee una longitud de registro de información y consistencia apropiada en la caracterización climática y
también en particular en la medición de la evaporación de tanque (ET), cuya variación mensual se presenta
en la siguiente figura.

140

120
Evaporación (mm)

100

80

60

40

20

0
ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC
MES
Figura 1-18. Variación de la evaporación promedio mensual en la estación Apto. Matecaña

Fuente: Hidrased, 2019.

1-44

1-44
A continuación, en la siguiente tabla, se presentan los resultados obtenidos de las metodologías para la
ubicación de la estación Aeropuerto Matecaña así como la comparación con el promedio multianual de los
registros de tanque de evaporación (ET) de la misma.

Tabla 1-18 Comparativo entre los métodos de estimación de la evaporación y la estación Apto
Matecaña
Método Evaporación promedio (mm/año) Error relativo (%)
Estación A. Matecaña (ET) 1252 -
CENICAFÉ (ETP) 776 38%
Thornthwaithe (ETP) 833 33%
Budyko-CENICAFÉ (ETR) 717 43%
Budyko-Thornthwaithe (ETR) 765 39%
Coutagne (ETR) 962 23%
Turc (ETR) 1202 4%

Fuente: Hidrased, 2019.

A partir de la tabla comparativa anterior, es posible observar que las metodologías de CENICAFÉ,
Thornthwaithe y la formulación de Budyko, subestiman la evaporación de la zona cercana a la estación
Aeropuerto Matecaña, mientras que la metodlogías de Coutagne y Turc se acercan en su orden de
magnitud, siendo esta última la que mejor ajuste presenta en cuanto a los registros multianuales promedio
de la misma.

De esta manera, en la Figura 1-19 se muestra de forma espacial la variación de la evapotranspiración real
(ETR) obtenida a partir de la metodología de Turc (1954)

1-45

1-45
 Figura 1-19. Evapotranspiración real estimada con el método de Turc

Fuente: Hidrased, 2019.

1.7 Caracterizar régimen hidrológico y Oferta hídrica

Con el fin de caracterizar la oferta hídrica de la zona, tanto de forma local como regional, se realizó una
comparación entre los valores de precipitación promedio y evaporación promedio multianual de la estación
Aeropuerto Matecaña, así como también un balance hídrico a nivel espacial para la cuenca del río Otún.

A continuación, en la Figura 1-20 y en la Tabla 1-19 para la estación Aeropuerto Matecaña, se muestra la
variación promedio mensual de la precipitación, junto a la evaporación media registrada y su
correspondiente diferencia.

1-46

1-46
 350

300
Lámina promedio (mm)

250

200

150

100

50

0
ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC
Mes
Pr (mm) ET (mm)

Figura 1-20. Variación mensual promedio de la precipitación y la evaporación de tanque en la

estación Aeropuerto Matecaña

Fuente: Hidrased, 2019.

Tabla 1-19 Balance entre precipitación y evaporación en la estación aeropuerto Matecaña

Parámetro Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic
Pr (mm) 130,1 140,9 227,5 254,2 254,0 182,1 127,8 125,4 188,5 298,9 288,6 203,2
ET (mm) 128,4 118,9 129,9 111,8 103,5 102,1 116,0 123,2 114,1 113,9 107,7 115,4
P-ET (mm) 1,7 22,0 97,6 142,5 150,5 80,0 11,8 2,1 74,5 185,0 180,8 87,9

Fuente: Hidrased, 2019.

A nivel de cuenca, el balance hídrico promedio de una región a lo largo de un periodo largo consta de la
relación entre la precipitación, la evaporación y el área de drenaje aferente a cada punto, por lo que se
asume que la perdida de volumen por recarga o infiltración tiende a cero, resumiendo la oferta volumétrica
a la siguiente expresión:

Q  ( P  ETR) Acuenca

Ecuación 1.20 Ecuación de balance hídrico a largo plazo

Donde,

Q: caudal

1-47

1-47
P: precipitación

ETR: evapotranspiración real

A: área de la cuenca

Utilizando los resultados obtenidos de la precipitación total promedio multianual (numeral 1.4.1) y los
valores de la evapotranspiración real de Turc (numeral 1.6), se calculó el balance hídrico a nivel puntual
(pixel) y los resultados se muestran a continuación.

Figura 1-21. Balance entre Precipitación-Evapotranspiración real total anual

Fuente: Hidrased, 2019.

Los valores de balance entre la precipitación promedio total anual y la evapotranspiración total anual
estimada fueron acumulados siguiendo los patrones de drenaje de la cuenca y el área de cada pixel, y
realizando una equivalencia entre unidades, es posible obtener del balance un caudal promedio de balance
hídrico.

La Tabla 1-20 muestra la magnitud promedio del aporte volumétrico de caudal del río Otún, por medio de
balance hídrico, a la altura de la ubicación de la PTAR, así como los aportes individuales de cada
subcuenca de los afluentes principales.

1-48

1-48
 Tabla 1-20 Oferta y balance hídrico multianual de las cuencas de la zona

Rendimiento hídrico (l/s

Cuenca Área (km2) Q. Balance hídrico (m3/s)
/km2)

Río Otún 484,92 22,88 47,2

Rio Otún (PTAR) 409,41 20,40 49,8
Q. Barbo 60,95 3,78 62,1
Q. Dosquebradas 58,52 3,00 51,2
Q. Azul 85,75 2,02 23,6
Q. San José 26,99 1,56 57,8
Q. San Juan 21,67 1,39 64,0
Q. Bosque 32,22 0,74 22,9
R. Azul 11,70 0,28 24,2

Fuente: Hidrased, 2019.

De los resultados mostrados anteriormente se observa que los valores de caudal promedio a largo plazo
obtenidos utilizando el balance hídrico tienen el mismo orden de los caudales promedios diarios
encontrados en el análisis de caudal medio del río Otún y que son detallados en el siguiente numeral.

1-49

1-49
1.8 Régimen de Caudales del Río Otún
A continuación, se realiza el análisis de caudales característicos del río Otún.

1.8.1 Caudales medios río Otún

Los caudales medios del río Otún se analizaron con la estación La Bananera, la cual, representa más del
50% del área de la cuenca aferente a la zona de estudio, por lo que se previó que esta estación podría
brindar una muy buena aproximación al comportamiento de los caudales de la planta de tratamiento de
aguas residuales (PTAR). Los caudales medios a nivel diario se muestran en la Figura 1-22 y a nivel
mensual en la Figura 1-23

120
Estación La Bananera Caudal Medio = 12.6 m3/s
100

80
Cadual (m 3/s)

60

40

20

0
1971 1977 1983 1989 1995 2001 2007 2013
Fecha

Figura 1-22. Caudal medio diario registrado en la estación La Bananera

Fuente: Hidrased, 2019.

En la Figura 1-23 se observa que el comportamiento de caudales presenta una bi-modalidad, es decir, dos
periodos secos y dos periodos lluviosos. Los periodos secos se identificaron entre enero-febrero y julio-
septiembre, mientras que los lluviosos se dan entre abril-mayo y noviembre-diciembre. El caudal medio de
la bananera es de aproximadamente 12.6 m3/s. En la Figura 1-24 se observa la curva de duración de
caudales en la estación La Bananera.

1-50

1-50
 20

18

16

14
Caudal (m 3/s)

12

10

0
ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC
Mes

Figura 1-23. Caudal medio mensual registrado en la estación La Bananera

Fuente: Hidrased, 2019.

Finalmente, con el caudal medio en la estación La Bananera se estima el caudal medio en la PTAR,
mediante transposición con respecto al área (ver Ecuación 1.21), el cual da un valor de caudal medio de
18,6 m3/s.

Area PTAR
𝑄 𝑃𝑇𝐴𝑅 = 𝑄 𝐵𝑎𝑛𝑒𝑛𝑒𝑟𝑎 ( )
Area bananera
Ecuación 1.21 Transposición con base en el área

1-51

1-51
 100%
90%
Fracción del tiempo en cual el
caudal es mayor al indicado

80% Serie Total

70% Periodo seco

60% Periodo Lluvioso

50% Caudal medio

40%
30%
20%
10%
0%
0 10 20 30 40 50 60 70 80 90 100
Caudal (m3/s)
Figura 1-24. Curva de duración de caudales en la estación La Bananera

Fuente: Hidrased, 2019.

1.8.2 Caudales mínimos río Otún

Para la estimación de caudales mínimos en el río Otún a la altura de la planta de tratamiento de aguas
residuales (PTAR) se realizó inicialmente un análisis de frecuencia de los caudales mínimos registrados
en la estación La Bananera operada por el IDEAM y ubicada sobre el río Otún. El área aferente a esta
estación corresponde a más del 50 % del área total de la cuenca aferente a la PTAR, por lo tanto, el análisis
de frecuencia de sus registros puede entregar muy buenas aproximaciones del comportamiento de los
caudales mínimos en la cuenca.

Con base en lo anterior se realizó un análisis de frecuencia considerando funciones de distribución de

probabilidad Gumbel y LogGumbel para la serie de caudales mínimos diarios anuales de la estación La
Bananera operada por el IDEAM, de las cuales se seleccionó la función de distribución que presentó el
mejor ajuste sobre el río Otún para diferentes períodos de retorno.

Adicionalmente, para complementar los estimativos de caudales mínimos por análisis de frecuencia se
aplicaron ecuaciones de regionalización. Las ecuaciones utilizadas fueron EPM (2005: Método I, Método
III, IIIA y IIIB) y Bolaños (1995) -Método III y Método IV-, las cuales fueron generadas para la aplicación en
todo el departamento de Antioquia, sin embargo, se ha analizado que son funcionales en otras regiones.

Las ecuaciones de regionalización mencionadas recurren a parámetros propios de la cuenca tales como
el área, la precipitación promedia multianual, la longitud del cauce principal y la pendiente media de la
cuenca. Con base en tales valores se obtienen expresiones para la estimación de valores mínimos medios

1-52

1-52
y de desviación estándar, los cuales, en la ecuación general de frecuencia, junto con el cálculo del factor
de frecuencia, permiten estimar caudales mínimos para varios períodos de retorno, como se indica a
continuación:

QTR  ˆ  KTR * ˆ
Ecuación 1.22 Ecuación general de frecuencia

Donde,

KTR: factor de frecuencia

̂ : media estimada de los caudales (m 3/s)

ˆ : desviación estándar estimada de los caudales (m 3/s)

   1  
K TR   0,45  0,779 * ln   ln   
   TR  
Ecuación 1.23 Factor de frecuencia

Donde,

TR: periodo de retorno (años)

Las expresiones para la estimación de los valores medios y las desviaciones estándar son las siguientes:
 EPM (2005):

Método I:

Media  0,01114 * A1,087

Desv _ est .  0,002948 * A1,087
Ecuación 1.24 Método I EPM (2005)

Método III, IIIA y IIIB:

Media  0,01114 * A1,087

Desv _ est .  0,002948 * A1,087
Ecuación 1.25 Método III EPM (2005)

1-53

1-53
 Media  0,01114 * A1,087
Desv _ est .  0,002948 * A1,087
Ecuación 1.26 Método IIIA EPM (2005)

Media  10 5,281 * A 0,903 * Pma 1,086

Desv _ est .  10 5,437 * A 0,683 * Long 0,516 * Pend 0,322 * Pma 0,868
Ecuación 1.27 Método IIIB EPM (2005)
 Bolaños (1995):

Método III:

Media  10 5,155 * A 0,942 * Pma 1,025

Desv _ est .  10 5,464 * A 0,997 * Pend 0,218 * Pma 0,894
Ecuación 1.28 Método III Bolaños (1995)

Método IV: regionalización para todo el departamento

Media  10 5,155 * A 0,942 * Pma 1,025

Desv _ est .  10 5,301 * A 0,991 * Pma 0,889
Ecuación 1.29 Método IV Bolaños (1995)

Donde,

A: área de la cuenca (km)

Pma: precipitación media multianual (mm/año)

Pend: pendiente media del cauce principal (%)

Long: longitud del cauce principal (km).

Realizando el análisis de frecuencia y las metodologías de regionalización anteriormente mencionadas se

estimaron los caudales mínimos en la estación La Bananera para diferentes períodos de retorno, tal y
como se presentan en la Tabla 1-21, donde se aprecia que los órdenes de magnitud son similares, Sin
embargo, el ajuste de Gumbel para periodo de retorno altos genera valores negativos y las dos
metodologías de regionalización aunque se ajustan muy bien al comportamiento de los datos, estiman

1-54

1-54
valores menores a los registrados para periodos de retorno bajos, por lo anterior la que mejor ajuste
presenta es la distribución de frecuencia de Log Gumbel ( línea verde de la Figura 1-25)

Tabla 1-21 Caudales mínimos en la estación La Bananera

Periodo de retorno Caudal mínimo (m3/s) en la estación La Bananera
(años) Gumbel Log Gumbel EPM Bolaños
1,01 8,2 11,8 8,3 7,2
2,33 4,6 4,1 3,6 3,0
5 3,1 2,6 2,9 2,4
10 1,9 1,8 2,6 2,1
25 0,3 1,1 2,2 1,8
50 -0,9 0,8 2,0 1,6
100 -2,1 0,5 1,8 1,4

Fuente: Hidrased, 2019.

14 Gumbel (La Bananera)

Log Gumbel (La Bananera)
12 BOLAÑOS
EPM
10 La Bananera
Caudal(m 3 /s)

0
1.01 1.1 1.25 1.5 2.33 5 10 25 50 100
Periodo de retorno ( Años) Papel probabilístico Gumbel

Figura 1-25. Caudales mínimos estimados en La Bananera

Fuente: Hidrased, 2019.

Finalmente, los caudales mínimos para diferentes periodos de retorno se requieren en el punto de la PTAR,
por lo que se aplica una transposición de los caudales mínimos desde la estación La Bananera hasta el
punto de interés, para los análisis de frecuencias esta transposición se realiza con base en el área de las
cuencas como se indica en la Ecuación 1.30, mientras que para la regionalización no es necesaria la
transposición dado que con la información de la cuenca en la PTAR se aplican los métodos. Los valores
de caudales mínimos para diferentes periodos de retorno en la PTAR se presentan en la Tabla 1-22.

1-55

1-55
 Area PTAR
𝑄 𝑇𝑟 𝑃𝑇𝐴𝑅 = 𝑄 𝑇𝑟 𝐵𝑎𝑛𝑒𝑛𝑒𝑟𝑎 ( )
Area bananera
Ecuación 1.30 Transposición de análisis de frecuencia con base en el área

Tabla 1-22 Caudales mínimos estimados en la PTAR

Caudal mínimo (m3/s) en PTAR
Periodo de retorno (años)
Gumbel Log Gumbel EPM Bolaños
1,01 12,0 17,4 13,5 12,7
2,33 6,8 6,1 5,8 5,5
5 4,6 3,8 4,7 4,4
10 2,7 2,6 4,1 3,8
25 0,4 1,6 3,5 3,3
50 -1,3 1,1 3,2 3,0
100 -3,1 0,8 2,9 2,7

Fuente: Hidrased, 2019.

Con base en lo anterior se concluye que el caudal mínimo para el río Otún en el sitio de PTAR para un
periodo de retorno de 2,33 años es de aproximadamente 6 m 3/s.

1.8.3 Caudal ecológico río Otún

Para estimar el caudal mínimo ecológico o caudal mínimo remanente que se debe conservar en los casos
de desviación o aprovechamiento del recurso hídrico de una corriente de agua, y que corresponde al
caudal requerido para el sostenimiento del ecosistema, la flora y la fauna de la misma, existen diversas
metodologías según la Resolución 0865 de 2004 expedida por el Ministerio de Ambiente, Vivienda y
Desarrollo Territorial (MAVDT) y según el informe sobre la Metodología para la estimación del caudal
ambiental en proyectos licenciados, elaborado por la Universidad Nacional de Colombia, sede Bogotá,
para el MAVDT en 2008. Estas metodologías pueden clasificarse de la siguiente forma:

 Metodologías hidrológicas.
 Metodologías con enfoque ecológico: hidrológicas con enfoque ecológico, hidráulicas con
enfoque ecológico, de simulación de hábitat fluvial y holísticas o funcionales.
 Metodologías que consideran explícitamente la calidad del agua.

Dentro de las metodologías hidrológicas ampliamente utilizadas se encuentra el porcentaje de descuento,

adoptada por el IDEAM, que define como caudal mínimo ecológico el valor correspondiente al 25% del
caudal medio mensual multianual más bajo de la corriente en estudio (ver Ecuación 1.31).

25
𝑄 𝑒𝑐𝑜𝑙𝑜𝑔𝑖𝑐𝑜 = 𝑄 𝑀𝑖𝑛 𝑀𝑒𝑛𝑠𝑢𝑎𝑙 ×
100
Ecuación 1.31 Caudal ecológico

Donde,

1-56

1-56
Qecologico: caudal ecológico (m 3/s)

RMinMesual: caudal medio mensual multianual más bajo (m 3/s)

De lo anterior y con base en los registros de caudales en la estación La Bananera, se estimaron los
caudales medios mensuales multianuales (ver Tabla 1-23). El caudal medio mensual multianual más bajo
corresponde al valor de agosto con 7,6 m3/s y por tanto el caudal mínimo ecológico en el río Otún a la
altura de la estación La Bananera es de 1,86 m 3/s (25% del caudal medio mensual multianual más bajo de
la corriente).

Tabla 1-23 Caudales medios mensuales multianuales registrados en la estación La Bananera

Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic
11,5 10,9 12,2 15,0 15,5 11,2 8,3 7,6 9,0 15,2 18,8 16,4

Fuente: Hidrased, 2019.

Finalmente, el objetivo es estimar el caudal ecológico en el río Otún a la altura de la planta de tratamiento
de aguas residuales, por lo que se transpone el caudal hallado de la estación La Bananera a la PTAR con
base en el área de las cuencas. El caudal ecológico en la PTAR es de 2,78 m 3/s

1.8.4 Caudales máximos río Otún

En este numeral se presentan los caudales máximos para diferentes periodos de retorno estimados tanto
por análisis de frecuencia como por modelo lluvia escorrentía (HEC-HMS). Estos valores son comparados
con los caudales máximos obtenidos en estudio previos.

1.8.4.1 Análisis de frecuencia

Con el objetivo de estimar los crecientes el río Otún se realizó un análisis de frecuencia a la serie de
caudales máximos de la estación La Bananera operada por el IDEAM, con las funciones de distribución
de probabilidad Gumbel, Log Gumbel, Pearson, Log Pearson, Normal y Log Normal. En la Figura 1-26 se
presenta los registros de la estación como las distribuciones de frecuencia anteriormente mencionadas.

1-57

1-57
 300
LA BANANERA Normal
Pearson Gumbel
Caudal maximos instantaneo (m 3/s)

250 Log Gumbel Log Normal

Log Pearson

200

150

100

50

0
1.01 1.1 1.5 2.33 5 10 25 50 100
Período de retorno en años
Papel probabilístico Normal

Figura 1-26. Caudales máximos para diferentes ajustes

Fuente: Hidrased, 2019.

Adicionalmente se realizaron pruebas teóricas de ajuste, para determinar cuál función representa de
manera teórica la distribución de los caudales de la estación La Bananera, este análisis dio como resultado
que la distribución Pearson es la que mejor representa los caudales, como SE muestra en la Tabla 1-24.

Tabla 1-24 Pruebas teóricas de ajuste en orden de jerarquía para caudales

2
Jerarquía Smirnov-Kolmogorov
1 LogNormal Pearson
2 LogPearson Normal
3 Gumbel LogNormal
4 Pearson LogPearson
5 LogGumbel Gumbel
6 Normal LogGumbel

Fuente: Hidrased, 2019.

En la Tabla 1-25 se presenta los caudales máximos estimados para diferentes periodos de retorno en la
estación La Bananera.

Tabla 1-25 Caudales máximos en la estación La Bananera para diferentes periodos de retorno

1-58

1-58
 Periodo de retorno (años) 2,33 5 10 25 50 100
Caudal máximo (m3/s) 140 170 190 212 227 241

Fuente: Hidrased, 2019.

Comparando estos estimativos de caudales máximos en la estación la Bananera con los obtenidos en
estudios previos (ver Tabla 1-26), se tiene un valor de un orden de magnitud similar al presentado en el
POMCA del 2016. Lo anterior confirma que el valor estimado en el presente estudio es acorde con otros
estudios recientes.

Tabla 1-26 Caudales máximos (m3/s) en la estación La Bananera de estudios previos y actual
Periodo de retorno
HIDRASED 2019 CARDER 1986 POMCA 2016
(años)
2,33 140
5 170 183 167
10 190 213 191
25 212 255
50 227 288 238
100 241 325 257

Fuente: Hidrased, 2019.

Adicionalmente, en 2014 la CARDER presentó un estudio donde se estimaron los caudales característicos
de las principales cuencas de todo el departamento y se determinaron los caudales tanto mínimos como
máximos del río Otún en la desembocadura con el río Cauca. Así las cosas, realizando la transposición
con base en el área de los caudales estimados con análisis de frecuencia por el estudio actual y
comparándolos con los estimados por la CARDER (ver Tabla 1-27) se observa que los caudales
entregados por este último son menores.

Tabla 1-27 Caudales máximos (m3/s) del río Otún en la desembocadura con el río Cauca
Periodo de retorno (años) HIDRASED 2019 CARDER 2014
Transposición LOGNORMAL GUMBEL
2,33 244 178 178
5 295 222 222
10 330 255 258
25 370 291 303
50 396 314 336
100 420 332 370

Fuente: Hidrased, 2019.

Finalmente, teniendo el análisis de caudales máximos del río Otún en la estación la Bananera y en la
desembocadura con el río Cauca, se procede a estimarlos en la ubicación de la PTAR, pero primero se
construyó el modelo lluvia escorrentía como se muestra a continuación.

1-59

1-59
1.8.4.2 Modelo lluvia escorrentía en HEC - HMS

Para la estimación de caudales máximos mediante lluvia escorrentía se construyó un modelo hidrológico
en el software HEC-HMS (Hydrologic Modeling System) desarrollado por el cuerpo de ingeniero de Estados
Unidos (United States Army Corps of Engineers -USACE) para el modelamiento hidrológico y está
diseñado para simular los procesos hidrológicos completos del sistema de cuencas hidrográficas
dendríticas. Es un modelo ampliamente usado en el medio.

Para este caso se utilizaron las configuraciones y métodos presentadas en la Tabla 1-28 para el análisis
de crecientes del río Otún a la altura de la zona de la PTAR. Con el fin de tener un mayor detalle sobre
cada uno de los métodos, se puede consultar el Manual de Usuario del HEC-HMS en su versión 4.2
(Brunner, 2016)

Tabla 1-28 Métodos utilizados para el modelo HEC-HMS

Método de transformación Clark Unit Hydrograph
Método de Perdida SCS Curve Number (CN)
Caudal Base No
Modelo de precipitación Frecuency Storm

Fuente: Hidrased, 2019.

I. Parámetros de las subcuencas

Inicialmente se definieron las subcuencas en las cuales se fraccionó el modelo hidrológico, para este caso
se realizó la partición de la cuenca en 7 subcuencas, 4 de ellas de alta pendiente, 2 de pendiente media y
1 baja pendiente. En la Figura 1-27 se puede observar las subcuencas en la que fue dividida el área de
estudio hidrológica, y en la Tabla 1-29 se observan los parámetros morfométricos tales como: área,
perímetro, pendiente, entre otras.

1-60

1-60
 Figura 1-27. Subcuencas del modelo HEC-HMS

Fuente: Hidrased, 2019.

Tabla 1-29 Parámetros morfométricos de las subcuencas

Área Cota min Cota Max Delta cota Longitud De Pendiente
Subcuenca
(km2) (msnm) (msnm) (m) Cauce (m) Media (%)
Q. SAN JUAN 21,7 1676,0 3581,0 1905,0 12702 15,0%
Q. DOSQUEBRADAS 58,5 1289,0 2248,0 959,0 15300 6,3%
Q. SAN JOSE 27,0 1489,0 3156,0 1667,0 17262 9,7%
Q. AZUL 85,7 3035,0 4937,0 1902,0 18360 10,4%
Q. BARBO 61,0 1837,0 4479,0 2642,0 19447 13,6%
EST. LA BANANERA 109,7 1568,0 3035,0 1467,0 29888 4,9%
PTAR 45,8 1124,0 1568,0 444,0 24139 1,8%

Fuente: Hidrased, 2019.

II. Duración

La duración de las tormentas máximas para diferentes períodos de retorno se consideró igual al tiempo de
concentración de cada subcuenca, ya que se ha observado que esta duración es la más crítica para la
generación de las crecientes, porque en este tiempo toda la escorrentía producida en cada sitio de la

1-61

1-61
cuenca alcanza a sumarse en el punto de interés, por lo que teóricamente, duraciones mayores no generan
un aumento en el caudal pico de la hidrógrafa.

El tiempo de concentración de las cuencas estudiadas fue calculado como el promedio de los valores
obtenidos mediante expresiones propuestas por diferentes autores como Kirpich (1940), Ventura-Heron
(MATA-LIMA, y otros, 2007), Giandotti (1940) , U.S.C.E (Feldman, 2000), Bransby (MOTH, 1998) , y Chow
(1998), las cuales son función de los parámetros morfométricos de las cuencas, como son el área, la
longitud del cauce principal y la diferencia máxima de altura. Las expresiones utilizadas para estimar este
parámetro fueron las siguientes:

Kirpich (1990):

0,77
 
 L 
T  0,066  
c  S 
 o 
Ecuación 1.32 Tiempo de concentración según Kirpich

Donde,

Tc: tiempo de concentración (horas)

L: longitud desde la estación de aforo hasta la divisoria, siguiendo el cauce principal (km)

So: relación entre la diferencia de cotas entre los puntos extremos de la cuenca y la longitud L -
Pendiente media del cauce-(m/m)-
Ventura-Heron:

0,75
 L 
T  0,30   
c  S 0,25 
 o 
Ecuación 1.33 Tiempo de concentración según Ventura-Heron

Donde,

Tc: tiempo de concentración (horas)

L: longitud del cauce principal (km)

So: pendiente media del cauce (m/m)

Giandotti (1990):

1-62

1-62
 4  A  1,5  L
T 
c 25,3  L  S
o
Ecuación 1.34 Tiempo de concentración según Giandotti

Donde,

Tc: tiempo de concentración (horas)

A: área de la cuenca (km 2)

L: longitud del cauce principal (km)

So: pendiente media del cauce (m/m)

U.S.C.E:

0,76
 L 
T  1,67  0,126   
c  S 0,25 
 o 
Ecuación 1.35 Tiempo de concentración según U.S.C.E.

Donde,

Tc: tiempo de concentración, en horas.

L: longitud del cauce principal, en kilómetros.

So: pendiente media del cauce, en m/m.

Bransby:

14,7  L
T  0,1
c A  S 0,2
o
Ecuación 1.36 Tiempo de concentración según Bransby

Donde,

Tc: tiempo de concentración, en horas.

A: área de la cuenca, en km 2.

L: longitud del cauce principal, en kilómetros.

1-63

1-63
So: pendiente media del cauce, en m/m.

Chow (1994):

0,64
 L 
T  1,67  0,123   
c  S 0,5 
 o 
Ecuación 1.37 Tiempo de concentración según Chow

Donde,

Tc: tiempo de concentración (horas)

L: longitud del cauce principal (km)

So: pendiente media del cauce (m/m)

Con base en los parámetros morfométricos de las subcuencas se estimaron los tiempos de concentración
de cada una y el Lag time. El Lag Time el tiempo entre el centro de masa del hietograma y el pico de la
hidrógrafa, y se define como 0.6 veces el tiempo de concentración. Los resultados se presentan en la Tabla
1-30.

Tabla 1-30 Tiempos de concentración para cada una de las subcuencas

Tiempo de concentración (min) Tiempo Lag
Subcuenca concentración Time
B2 C3 G4 K5 U6 V7 Asumido
(horas) (horas
Q. San Juan 201 115 65 58 125 173 173 2,88 1,73
Q. Dosquebradas 260 171 130 94 170 234 234 3,90 2,34
Q. San Jose 291 161 86 87 172 236 236 3,94 2,36
Q. Azul 272 164 111 89 177 244 244 4,07 2,44
Q. Barbo 282 156 88 84 176 242 242 4,04 2,42
La Bananera 502 284 170 172 296 405 405 6,75 4,05
PTAR 538 340 225 213 303 415 415 6,91 4,15

Fuente: Hidrased, 2019.

De la tabla anterior, los tiempos de concentración seleccionados son los estimados por medio de la
ecuación propuesta por Ventura, la cual fue desarrollada con base en el estudio de cuencas intermedias

2 Bransby
3 Chow
4 Giandotti
5 Kirpich
6 U.S.C.E.
7 Ventura

1-64

1-64
rurales, y que contrasta con otras formulaciones, como la Giandotti cuyo estudio compila cuencas muy
variables en extensión (170-70.000 km2), o Kirpich la cual se deriva de cuencas muy pequeñas (0.004-
0.453 km2)

III. Parámetros de perdida

Para estimar las pérdidas se utilizó el método propuesto por el U. S. Soil Conservation Service, Numero
de Curva (CN), el cual es función de la cobertura vegetal, del tipo y uso del suelo, de la impermeabilidad
del terreno y su humedad antecedente.

El número de curva depende del tipo del suelo, del tipo de tratamiento superficial que éste presente, de su
condición hidrológica y de su humedad antecedente que puede relacionarse con la precipitación en los
cinco días anteriores. En relación con la humedad antecedente del suelo se diferencian tres grupos:

AMC-I: Para suelos secos.

AMC-II: Para suelo intermedios.

AMC-III: Para suelos húmedos.

Para tener en cuenta las condiciones de humedad antecedente a la lluvia máxima de cada período de
retorno, se consideró la condición Tipo II, que corresponde a una humedad promedia.

Además, se requiere ubicar la zona dentro de un grupo hidrológico del suelo, el cual está definido como el
potencial de infiltración después de una lluvia prolongada. El tipo de suelo se clasifica en cuatro categorías
hidrológicas, según el método del SCS. En orden decreciente de permeabilidad las categorías hidrológicas
son:
 Grupo de Suelo A: alta infiltración (baja escorrentía). Suelos derivados de rocas metamórficas
cubiertos con vegetación (bosque o rastrojo alto) con grado de meteorización 30/50 según Brand
(1988) y con discontinuidades en la matriz de suelo producto del proceso de descomposición de la
roca. Las estructuras heredadas funcionan como canales de flujo principales. La pendiente en este
grupo de suelo debe ser inferior a un 7%. También se incluyen los suelos que presentan un efecto
geológico marcado por diaclasas, bandeamientos, y cuyo buzamiento se da en ángulos mayores a
60°.
 Grupo de Suelo B: infiltración moderada (escorrentía moderada). Suelos derivados de rocas
metamórficas, ígneas o sedimentarias poco denudados con grado de meteorización 0/30 según Brand
(1988) y con discontinuidades difícilmente cartografiables en campo. Las estructuras heredadas ya
no son tan importantes en la permeabilidad del conjunto, por lo cual la permeabilidad primaria es la
que controla el flujo de agua en el suelo. La pendiente en este grupo de suelo debe ser inferior a un
10%.
 Grupo de Suelo C: infiltración baja (escorrentía de moderada a alta). Suelos residuales derivados de
cualquier tipo de roca, y cuyo grado de meteorización se clasifique como residuales maduros (SR,
Brand – 1988), en los cuales las propiedades y minerales del material parental no son fácilmente

1-65

1-65
 identificables. Esta formación tiene como característica principal que está cubierta por pastos
manejados y sus pendientes varían entre un 5 y 15%.
 Grupo de Suelo D: muy poca infiltración (alta escorrentía). Suelos derivados de rocas metamórficas,
ígneas y sedimentarias, cuyo grado de meteorización varía entre 50/90 y en donde la pendiente media
excede 15%. El tipo de vegetación en este grupo no es de mayor importancia ya que el fuerte gradiente
topográfico impide la recarga de los acuíferos en la zona, generando por el contrario un flujo
hipodérmico hasta aflorar nuevamente en superficie y hacer parte del agua que se aporta al canal
principal de la cuenca.

El número de curva fue definido ponderando en el área los usos del suelo de cada cuenca utilizando
herramientas SIG, y tomando como información base el Corine Land Cover Colombia (IDEAM, 2007) que
permite describir, caracterizar, clasificar y comparar las características de la cobertura de la tierra,
interpretadas a partir de la utilización de imágenes de satélite de resolución media (Landsat), para la
construcción de mapas de cobertura a diferentes escalas. En la Tabla 1-31 se presentan los valores de
referencia de los números de curva, según uso del suelo y grupo hidrológico para la condición de humedad
antecedente Tipo II.

Tabla 1-31 Número de curva para diferentes usos del suelo según condiciones de humedad
antecedente Tipo II

1-66

1-66
 Grupo Hidrológico del
Descripción de usos del suelo Suelo
A B C D
Zona de parqueo pavimentada, Techada. 98 98 98 98
Calles y carreteras:
Pavimentos con bordillos y alcantarillas 98 98 98 98
Afirmado 76 85 89 91
Suelo compactado 72 82 87 89
Cultivos (Cosecha agrícola) del suelo:
Sin tratamientos de conservación (Sin terrazas) 72 81 88 91
Con tratamientos de conservación (Terrazas, Contornos) 62 71 78 81
Pastos o gama de suelos
Pobre (< 50% cobertura de los suelos por pastos muy continuos) 68 79 86 89
Buena (50-75% Cobertura del suelo por pastos poco continuos) 39 61 74 80
Prado (pasto sin pastar y corto para follaje) 30 58 71 78
Rastrojo (bueno, >75% cobertura de suelo) 30 48 65 73
Bosques y selvas
Pobre (árboles pequeños – rastrojo destruido después de quema o
45 66 77 83
rozada)
Regular (rozado pero sin quema, algunos rastrojos) 36 60 73 79
Bueno (sin rozado, cobertura de suelo rastrojo) 30 55 70 77
Espacios abiertos (pastos, parques, canchas de golf, cementerios, etc.):
Regular (cobertura de pastos 50-75% del área) 49 69 79 84
Bueno (cobertura de pastos >75% del área) 39 61 74 80
Zonas comerciales y de negocios (85% impermeable) 89 92 94 95
Zonas industriales (72% impermeable) 81 88 91 93
Áreas residenciales:
Lotes de 505 m2, cerca de 65% impermeable 77 85 90 92
Lotes de 1.011 m2, cerca de 38% impermeable 61 75 83 87
Lotes de 2.023 m2, cerca de 25% impermeable 54 70 80 85
Lotes de 4.046 m2, cerca de 20% impermeable 51 68 79 84

Fuente: Chow (1994)

En las siguientes figuras se observan las coberturas del CLC para las cuencas de análisis y en la Figura
1-28 se presenta el número de curva CN estimado para cada una de las cuencas que componen el modelo
hidrológico.

1-67

1-67
 Figura 1-28. Coberturas CLC para las cuencas de análisis

Fuente: Hidrased, 2019.

En la figura anterior se presentaron todas las coberturas según el CLC, donde se puede apreciar más de
30 categorías. Para entender de una mejor manera la distribución espacial de la cobertura en la cuenca,
se realizó una simplificación gráfica, únicamente para efectos de presentación (ver Figura 1-29), donde se
puede observar la cobertura categorizada en 5 clases principales.

1-68

1-68
 Figura 1-29. Coberturas simplificadas CLC para las cuencas de análisis

Fuente: Hidrased, 2019.

Tabla 1-32 Numero de curva de las cuencas de análisis

Porcentaje
Numero de curva
Nombre Cobertura domínate cobertura
(CN)
dominante
Q. San Juan Bosque denso 72% 65,2
Q. Dosquebradas Cultivos permanentes arbustivos 32% 71,8
Q. San José Plantación forestal 37% 74,4
Q. Azul Herbazal 69% 72,6
Q. Barbo Bosque denso 71% 64,9
La Bananera Bosque denso 44% 67,5
PTAR Cultivos permanentes arbustivos 20% 74,9

Fuente: Hidrased, 2019.

IV. Parámetros de precipitación

El método para el análisis de la precipitación corresponde a - frecuency storm - el cual requiere la magnitud
de la lluvia en milímetros para diferentes duraciones de cada uno de los periodos de retorno analizados.

1-69

1-69
Inicialmente, se realizó el análisis de frecuencia de las estaciones pluviométricas de la zona, como se
indicó en el capítulo de climatología. Con base en esta precipitación y en el porcentaje de cuenca aferente
de cada estación definido a partir de los polígonos de Thiessen, se estimó la lluvia máxima en 24 horas
para diferentes periodos de retorno que representa la climatología de toda la cuenca y los resultados se
presentan en la Tabla 1-33.

Tabla 1-33 Precipitaciones máximas y % de área aferente

Aero
Estación Termales Boquerón Playa rica Pez fresco La Laguna
Matecaña
% Área 5% 18% 15% 25% 4% 34%
Tr (años) Gumbel Pearson Gumbel Pearson Pearson Log Pearson
1,01 46 11 56 45 48 17
2,33 79 101 86 96 80 34
5 94 115 99 107 91 41
10 106 122 109 114 98 46
20 117 127 120 120 105 51
25 121 128 123 121 107 53
30 124 129 125 123 109 54
50 132 132 133 126 113 57
100 143 134 143 130 119 61
500 169 138 166 138 131 71
1000 180 139 176 140 136 75

Fuente: Hidrased, 2019.

Con los valores máximos de precipitación para diferente periodo de retorno, se procede a estimar la
magnitud de la precipitación para diferentes duraciones, en la Tabla 1-34 se presenta estos valores que
son ingresados al modelo climatológico del HEC-HMS.

Tabla 1-34 Precipitación máximas en mm para diferentes duraciones y periodos de retorno

Tiempo 2,33 años 5 años 10 años 25 años 50 años 100 años
5 min 10 12 13 14 15 16
15 min 23 26 28 31 32 34
1 hora 38 43 47 51 54 57
2horas 46 52 57 62 65 68
3 horas 50 57 62 68 72 75
6 horas 58 66 72 78 83 87
12 horas 65 75 82 89 94 98
24 horas 73 84 91 99 105 110

Fuente: Hidrased, 2019.

Contando con los parámetros morfométricos, de transformación, de perdida de cada una de las
subcuencas y la caracterización de la precipitación máxima en la zona, se procede a la construcción del
modelo en el HEC HMS, en la Figura 1-30 se puede observar las conexiones entra cada una de las
subcuencas.

1-70

1-70
 Figura 1-30. Cuencas hidrográficas en HEC-HMS

Fuente: Hidrased, 2019.

El modelo hidrológico se simuló durante un día y medio, sin embargo, la precipitación máxima es en 24
horas, y por lo tanto las hidrógrafas resultantes también lo son. En la Figura 1-31 se presentan las
hidrógrafas de salida de HEC-HMS y en la Tabla 1-35 se presentan los valores de caudales máximos
obtenidos para la estación La Bananera y en la PTAR.

1-71

1-71
 500
2.33 Años
450 5 Años
10 Años
400
25 Años
350 50 Años
100 Años
Cadual (m 3/s)

300

250

200

150

100

50

0
00:00 03:00 06:00 09:00 12:00 15:00 18:00 21:00 00:00 03:00 06:00 09:00 12:00
Tiempo
Figura 1-31. Hidrógrafas para diferentes periodos de retorno

Fuente: Hidrased, 2019.

Tabla 1-35 Caudales máximos estimados con el HEC-HMS

Periodo de retorno (años) 2,33 5 10 25 50 100
Caudal máximo (m3/s) La Bananera 117 170 210 256 293 325
Caudal máximo (m3/s) PTAR 179 260 320 390 446 494

Fuente: Hidrased, 2019.

Los resultados obtenidos por el modelo lluvia escorrentía construido en HEC-HMS fueron comparados con
las estimaciones de caudales máximos a partir de la transposición del análisis de frecuencia (Ecuación
1.30) tal y como se muestra en la Tabla 1-36. A partir de los resultados se evidenció que el modelo lluvia
escorrentía tiende a sobrestimar un poco los caudales, aunque para el periodo de retorno de 2,33 años el
valor de caudal máximo es menor. De acuerdo a lo anterior los caudales seleccionados e hidrógrafas para
la modelación hidráulica son los entregados por el modelo HEC-HMS dado que es la metodología más
conservadora.

Tabla 1-36 Caudales máximos para diferentes periodos de retorno en la PTAR

1-72

1-72
 HIDRASED 2019 CARDER 2014
Periodo de
Transposición
retorno Transposición Transposición
análisis HEC-HMS
(años) LOGNORMAL GUMBEL
frecuencia
2,33 207 179 151 150
5 250 260 188 188
10 280 320 215 218
25 313 390 246 256
50 335 446 265 284
100 355 494 281 313

Fuente: Hidrased, 2019.

1.8.5 Análisis eventos extremos

A continuación, en los siguientes numerales, se complementa el análisis de eventos extremos tanto para
eventos máximos (torrenciales) como para mínimos (sequias). Esto se hizo con base en información
secundaria, como noticias o reportes.

1.8.5.1 Eventos torrenciales

Dentro de los registros de inundaciones históricas de la zona, es posible encontrar recuentos de noticias
relacionados con eventos puntuales donde se han dado altas precipitaciones en la cuenca del río Otún o
la del río Consota, presentándose desbordamientos, deslizamientos de ladera e inundaciones en zonas
tanto rurales como urbanas de la ciudad de Pereira.

1-73

1-73
 Figura 1-32. Crecientes del río Otún a la altura del viaducto César Gaviria Trujillo. Tomada
diario El Mundo, 15 noviembre de 2017.

En la Tabla 1-37 se muestra un resumen de algunos eventos históricos de precipitación, desbordamiento

de cauces, deslizamientos en temporada invernal o atención de desastres por riesgos hidráulicos que han
sido cubiertos por algunos diarios o cadenas radiales:

Tabla 1-37 Resumen de eventos hidráulicos, precipitaciones, desbordamientos cubiertos por

la prensa
Año Época Evento Afectaciones Fuente
2019 Abril Lluvias e inundaciones Dosquebradas Caracol Radio
Deslizamientos y 80 personas afectadas
2018 Noviembre desbordamiento de la zona rural. 20 viviendas Caracol Cadena
Quebrada Florida afectadas
Intervención de 17
2018 Enero puntos críticos ante Caracol Cadena
temporada invernal
Alerta naranja Río Otún.,
2018 Julio Altas Precipitaciones 16 viviendas afectadas Caracol Cadena
por inundaciones

1-74

1-74
Año Época Evento Afectaciones Fuente
Villa del Prado, la
2017 Mayo Lluvias e inundación avenida Sur, la avenida El Tiempo
del Río
Desbordamiento del río EL Mundo, El
2017 Noviembre 300 damnificados
Otún y Consota tiempo
Inundaciones y colapso inundaciones en 12
2017 Mayo de sistema de barrios de Pereira y Caracol Cadena
alcantarillado colapso del alcantarillado
Lluvias y
Colapso de puente
2017 Enero desbordamiento rio Caracol Cadena
peatonal Dos Quebradas
Otún
Deslizamientos e
2016 Septiembre inundaciones laderas Caracol Cadena
Delmiro Otún
Reporte de riesgo El Nuevo siglo
2016 Agosto
cuenca rio Otún Periódico
Inundaciones en
2015 Noviembre viviendas afectadas RCN Radio
Dosquebradas
Alerta de
2013 Marzo precipitaciones El Tiempo
históricas 40 años
Crecientes rio Otún y Evacuación preventiva
2011 Octubre Caracol Cadena
Consota por creciente
Alerta amarilla por ríos
2009 Marzo 46 familias El Tiempo
Otún, Consota, osos
Sector de El Caracol-La
Curva, en la vía a
Armenia, y en los barrios
Lluvias torrenciales. Las Brisas y Villa
2008 Junio Declaración de urgencia Santana resultaron El Tiempo
manifiesta destechadas como
consecuencia de lluvias
acompañadas por
intensos vientos.
60% de las vías
Daño en vías de
2005 Junio afectadas por temporada El Tiempo
Risaralda por invierno
invernal
63 viviendas de Pereira
afectadas el jueves
pasado por el
Desbordamiento del rio desbordamiento del río
2002 Marzo El Tiempo
Otún Otún. Incendio en
subestación
Dosquebradas por el
desbordamiento

1-75

1-75
 Año Época Evento Afectaciones Fuente
Desbordamiento 80 familias, barrio
1996 Febrero quebradas Oso y rio Monserrate oriente de El tiempo
Consota Pereira
Desbordamiento rio Destrucción de 80
1979 Defensa Civil
Otún viviendas

Fuente: Adaptado de El Tiempo, RCN radio, Caracol cadena

Cronológicamente, de la tabla anterior, el evento más antiguo reportado data de 1979 y corresponde al
desbordamiento de río Otún y es además el evento de origen hidráulico más grande del afluente atendido
por la defensa civil de Risaralda. A raíz de este evento y de las problemáticas de riesgos hidráulicos
ligados al río Otún, en 1981 se crea la corporación autónoma regional de Risaralda (CARDER) con el fin
de prevenir y gestionar aspectos de los recursos hídricos, ecológicos y desarrollo de actividades
económicas del departamento.

En el año 1995 la misma entidad, junto la defensa Civil, reportan que los ríos Consota, Dosquebradas y la
quebrada El Oso se desbordan, ocasionando la afectación de 2500 viviendas, 831 personas evacuadas,
2 casos fatales, y 6 heridos.

Recientemente, la emergencia invernal de 2011 y de 2017 en los periodos de Abril-Mayo y Octubre-

Noviembre mostraron que todos los afluentes de la zona, y en particular el río Otún y el río Consota
provocaron la afectación de cientos de familias, directamente por inundaciones o por desestabilización de
las laderas de su valle, sin embargo, estos periodos invernales han sido reportados como recurrentes
durante los últimos 15 años, afectando directamente tanto la zona urbana como rural del municipio de
Pereira.

1-76

1-76
 Figura 1-33. Eventos torrenciales del río Otún y Consota en noviembre de 2017

Fuente: Tomadas periódicos el mundo, y Caracol Emisora. 15/11/2017

1.8.5.2 Eventos mínimos

Dentro de los períodos de estiaje, se reporta en artículos de periódicos, un incremento en los casos de
enfermedades como el Dengue (enero de 2010), disminución en el caudal del río Otún en agosto de 2015,
asociado a la oleada de calor de la temporada inicial de este año, y en 2016 una alerta en los niveles de
la laguna del Otún, ubicada en la parte alta de la cuenca y alimentada directamente por los páramos de la
zona.

1-77

1-77
1.9 Modelamiento Hidráulico
En este capítulo se presenta la metodología, los casos de estudio y los resultados estimados mediante
simulación hidráulica, para el tramo de estudio del río Otún.

1.9.1 Metodología para estimar los perfiles de flujo

La modelación hidráulica requerida para estimar las características de un tramo del río Otún en la zona
donde se plantea localizar el proyecto de la PTAR contó con información base de campo, la cual consistió
en información LiDAR y levantamiento topográfico de detalle.

Para el procesamiento de la información topográfica y para determinar los perfiles de flujo, así como los
demás parámetros hidráulicos asociados a los diferentes caudales, se utilizaron dos programas de
computación: El primero es el HEC-RAS, desarrollado por el U.S. Army Corps of Engineers, (Brunner,
2016) que resuelve las ecuaciones de energía y momentum, por el método estándar de pasos para
determinar los perfiles de flujo. El segundo programa usado para mapear la geometría del cauce y las
secciones fue el RAS Mapper, que es una plataforma de información geográfica integrada al modelo
hidráulico. Esto garantiza mayor precisión en los datos geométricos del modelo, tanto en el trazado
horizontal como en las secciones transversales, haciendo que el modelo hidráulico definitivo represente
de mejor manera las condiciones reales.

1.9.2 Consideraciones para el modelado

Los perfiles de flujo para cada uno de los escenarios estudiados se determinaron con base en las
siguientes consideraciones:

1.9.2.1 Rugosidad

La rugosidad de las llanuras de inundación se determinó teniendo en cuenta parámetros tales como:
material que compone el lecho (nb), estabilidad de las secciones del cauce (n 1), efecto relativo de las
obstrucciones (n3), Vegetación (n4) y efectos de meandro (m 5)

En relación con la rugosidad del canal y de las llanuras de inundación, se siguió la metodología propuesta
por el U.S. Geological Survey (Aldridge & Garrett, 1973) y Chow (1983) que tienen en cuenta que para las
llanuras de inundación se pueden tener rugosidades bastante diferentes a las del canal o cauce principal
de la corriente y por lo tanto deben evaluarse de manera independiente. Al igual que para el canal, en la
llanura de inundación, se establece una rugosidad base que se ajusta a través de diversos factores para
obtener la rugosidad total.

La rugosidad base se estimó a partir del diámetro medio, en campo se observaron rocas de hasta 0,30 m
(ver Figura 1-34), que según la ecuación de Strickler (Yen, 1992), da un valor de rugosidad del orden de
n=0,028.

1-78

1-78
 Figura 1-34. Material que compone el lecho del cauce del río Otún

Los factores que afectan la selección del coeficiente de rugosidad, n, son el tipo y tamaño de los materiales
que componen el lecho y las bancas del canal, así como su forma. Para estimar el efecto de estos factores,
Cowan (1956) propone la siguiente expresión para determinar el valor de n de un canal:

n  (nb  n1  n 2  n3  n 4 )m

Ecuación 1.38 Cálculo de la rugosidad de Manning

Dónde,

nb: rugosidad base para la sección del canal, según el material presente en el lecho del cauce.

n1: factor de corrección para el efecto de las irregularidades de la superficie.

n2: factor que tiene en cuenta las variaciones en forma y tamaño de la sección transversal del canal.

n3: factor para las obstrucciones en la sección.

n4: valor que considera la vegetación existente sobre el canal principal y las condiciones del flujo cuando
estas se encuentran sumergidas debido a crecientes extremas.

1-79

1-79
m: corrección por meandros en el canal principal, que depende de la relación entre la longitud total de
meandros en el tramo de estudio, respecto a la longitud de tramos rectos.

A continuación, se presenta en detalle la descripción de los factores de corrección definidos con

anterioridad para estimar el coeficiente n de rugosidad de Manning y se indican los valores sugeridos para
cada uno de ellos:

 n1: la corrección por irregularidades se debe tener en cuenta por la condición de bancas erodadas que
exponen aristas de material o incluso raíces de vegetación circundante. En la Tabla 1-38 se presentan
los valores de este factor, en donde suave corresponde a una sección completamente estable, menor
a una sección con bancas ligeramente socavadas o afectadas por procesos erosivos, moderado
representa procesos de erosión lateral frecuentes en el tramo de análisis, y severo indica que la sección
es completamente inestable a lo largo de todo el tramo de estudio.

Tabla 1-38 Coeficiente n1 en función del grado de irregularidad de la sección

Condiciones del canal Valor
Suave 0
Menor 0,005
Grado de irregularidad n1
Moderado 0,010
Severo 0,020

Fuente: modificada del U.S. Geological Survey (Aldridge & Garrett, 1973)

 n2: la corrección por las variaciones de la sección en cuanto a la forma y el tamaño se debe
principalmente a cambios bruscos por reducción en la capacidad hidráulica o cambios bruscos en la
dirección del flujo, pero en última instancia depende de la magnitud del cambio entre secciones (Benson
& Dalrymple, 1967). En la Tabla 1-39 se presentan los valores de este factor.

Tabla 1-39 Coeficiente n2 en función de la variación en la sección transversal del canal

principal
Valor
Condiciones del canal
Inferior Superior Medio
Variación Gradual 0 0,05 0,025
sección n2 Ocasionalmente alternante 0,050 0,01 0,030
transversal Frecuentemente alternante 0,010 0,015 0,0125

Fuente: modificada del U.S. Geological Survey (Aldridge & Garrett, 1973)

 n3: mediante este factor se consideran obstrucciones en la sección como la presencia de raíces, rocas
grandes, pilas, escombros e incluso patrones de flujo como cauces trenzados, los cuales incrementan
la rugosidad. Los niveles de obstrucción se clasifican como despreciables, menores, apreciables o
severos (Chow V. T., 1983) cómo se indica en la Tabla 1-40.

1-80

1-80
 Tabla 1-40 Coeficiente n3 en función del grado de obstrucciones presentes en el canal
principal
Valor
Condiciones del canal
Inferior Superior Medio
Insignificante 0 N/A 0
Menor 0,01 0,015 0,013
Efecto relativo de las obstrucciones n3
Apreciable 0,020 0,030 0,025
Severo 0,040 0,060 0,050

Fuente: modificada del U.S. Geological Survey (Aldridge & Garrett, 1973)

 n4: el incremento de la rugosidad por la vegetación existente depende del porcentaje del perímetro
mojado que se encuentra cubierto por la misma, de la densidad de la vegetación por encima de la línea
máxima de la creciente, y de su localización con respecto al sentido del flujo. En la Tabla 1-41 se
presentan los valores típicos de este coeficiente.

Tabla 1-41 Coeficiente n4 en función del grado de obstrucciones presentes en el canal

principal debido a la vegetación
Valor
Condiciones del canal
Inferior Superior Medio
Baja 0,005 0,01 0,0075
Media 0,010 0,025 0,018
Vegetación n4
Alta 0,025 0,050 0,038
Muy alta 0,050 0,100 0,075

Fuente: modificada del U.S. Geological Survey (Aldridge & Garrett, 1973)

m: la corrección por meandros es menor para el rango de 1,0 a 1,20, apreciable si está entre 1,20 y
1,50, y severa si es mayor a 1,50 (Chow V. T., 1983). En la Tabla 1-42 se presentan los valores típicos
para corrección del coeficiente de rugosidad por la presencia de meandros.

Tabla 1-42 Coeficiente m de corrección por meandros en el tramo de estudio

Condiciones del canal Valor
Menor 1,0-1,2
Efecto de los meandros del cauce m Apreciable 1,2-1,5
Severo >1,5

Fuente: modificada del U.S. Geological Survey (Aldridge & Garrett, 1973)

De acuerdo a lo anterior se procesó la información y se obtuvo la rugosidad, tal y como se presentan a

continuación:

1-81

1-81
 Tabla 1-43 Rugosidad seleccionada para el canal principal
nb Rugosidad base 0,028
n1 Grado de irregularidad 0,005
n2 Variación sección transversal 0,025
n3 Efecto relativo de las obstrucciones 0,013
n4 Vegetación 0,005
m5 Efecto de los meandros del cauce 1,0
ncanal Rugosidad general 0,076

Fuente: Hidrased, 2019.

Finalmente, la rugosidad para el lecho del canal se estimó en 0,076 y de 0,100 para las llanuras de
inundación que corresponde a un terreno con árboles, matorrales (Chow V. T., 1983).

1.9.2.2 Coeficientes de contracción y expansión

De otro lado, la selección de los coeficientes de contracción y expansión del flujo debido a variaciones
entre la geometría de dos secciones, las pérdidas locales debido a descargas o cambios de sección se
determinaron mediante las cartas del USACE ( (Hydrologic Engineering System , 1991). Los valores
seleccionados se muestran a continuación:

 Coeficiente de contracción: 0,10

 Coeficiente de expansión: 0,30

1.9.2.3 Condiciones de frontera y Leeves

Adicionalmente en la modelación, se emplearon “leeves” en los sitios en donde éstos generan restricciones
al nivel del flujo. Sin embargo, en este caso, dadas las características del tramo de estudio, en donde la
sección tiene una alta capacidad hidráulica, no fueron determinantes a la hora de definir la extensión de la
inundación.

Para la simulación hidráulica, se verificó que la condición de frontera no generara falsos forzamientos en
las cercanías de las mismas. En este caso se utilizó la profundidad normal aguas arriba y aguas abajo del
tramo de estudio con un valor medio de pendiente del lecho del orden de 0,018 m/m.

Adicional a lo anterior, los perfiles de flujo (para las crecientes de 2,33 hasta 100 años de periodo de
retorno) se calcularon para flujo permanente para las condiciones sub y súper críticas para identificar
adecuadamente el tipo de régimen hidráulico, y evitar falsos forzamientos.

1-82

1-82
1.9.3 Geometría del modelo

A partir de la información topográfica, se construyó la geometría del río Otún con ayuda del software HEC-
RAS. Como punto de partida, es recomendable que, para cauces con pendientes fuertes y caracterizados
por altas velocidades, se su utilicen espaciamientos del orden de 30 m, mientras que, para ríos con
pendientes suaves y uniformes, el espaciamiento puede ser de hasta 1500 m (USACE, 2010a)8.

Adicionalmente, al tener en cuenta para el espaciamiento la pendiente media del lecho (0,018 m/m) y de
profundidad la del canal formador (del orden de 3,0 m) se tiene que el espaciamiento es de 25 m. (Samuels,
1990). Sin embargo, las secciones del modelo quedaron con un espaciamiento variable, con el fin de tener
en cuenta detalles de cambios morfológicos del canal, por lo que se tienen espaciamientos mínimos de
hasta 13,5 m. El tramo de estudio quedó conformado por una longitud del orden de 2580 m, y se muestra
en la Figura 1-35.

Figura 1-35. Localización general de las secciones y del tramo de estudio

Fuente: Hidrased, 2019.

De forma complementaria, en la Figura 1-36 se muestra la disposición de las secciones en el tramo de la

PTAR, mientras la Figura 1-37 presenta el isométrico del modelo para todo el tramo de estudio, en donde
se observa las características de un canal definido en una sección con geometría encañonada y con pocas
llanuras inundables.

8En este sentido es importante mencionar que el tener secciones muy próximas, puede resultar en problemas de
inestabilidad debido a que la derivada con respecto a la distancia tiende a ser sobrestimada, especialmente cuando
el flujo está sobre el ramal ascendente de las hidrógrafas.

1-83

1-83
 Figura 1-36. Esquema de detalle del modelo en la zona de la PTAR

Fuente: Hidrased, 2019.

1-84

1-84
 Figura 1-37. Isométrico del modelo realizado en HEC-RAS

Fuente: Hidrased, 2019.

1.9.4 Casos de estudio

Este escenario corresponde al análisis de los niveles y velocidades para la geometría actual y las
crecientes máximas que fueron estimadas en el aparte hidrológico. Es de resaltar que para la creciente de
diseño (100 años de periodo de retorno) se realizó una mayoración del caudal para tener en cuenta el
efecto de flujos hiperconcentrados.

1-85

1-85
Algunos autores como Pierson (Pierson, 2005) definen el flujo hiperconcentrado como una etapa
intermedia entre la inundación (1) y el flujo de lodos y escombros (2). Para la creciente de análisis se
consideran aportes de sedimentos que no superan el 4% en volumen, mientras que, en el flujo de lodos y
escombros, este valor puede ser mayor al 60%. Algunos otros estudios consideran que un 40% de
mayoración al volumen de la hidrógrafa para casos extremos en regiones semi-áridas, está asociado a un
flujo hiperconcentrado. Finalmente, en un estudio más reciente (Rodríguez Espinosa, Córdoba Guerrero,
& Delgado-Granados, 2017) se hace mención a una concentración que varía entre el 35-40% en volumen
para los flujos hiperconcentradas. Para este caso se usó la creciente de 100 años mayorada en un 40%.

1.9.5 Análisis de resultados hidráulicos

En la Figura 1-38, se muestra el perfil de flujo simulado para una creciente de 2,33 años de período de
retorno, esto a manera de ejemplo, puesto que los resultados procesados se indican seguidamente. De
este resultado se destaca que para las crecientes analizadas el flujo tiende a ser de naturaleza subcritica.
En algunos sectores muy puntuales se generan algunos controles hidráulicos, debido principalmente a
estrechamientos del flujo.

Figura 1-38. Perfil de flujo de línea base para la creciente de 2,33 años de periodo de retorno

Fuente: Hidrased, 2019.

1-86

1-86
En la Figura 1-39 se muestran las profundidades máximas estimadas para diferentes periodos de
referencia (2,33 a 100 años de periodo de retorno). Las profundidades tienden a ser muy uniformes para
cada periodo de retorno, es decir que no tienen mucha variabilidad. Es de anotar que para la creciente de
2,33 años la profundidad máxima promedio en el sitio de la PTAR es del orden de 3,0 m, mientras que
para la creciente de 100 años de periodo de retorno este valor es de aproximadamente 5,0 m.

8.0
2.33 10 25 50 100 100 mayorada
7.0

6.0
Profundidad (m)

5.0

4.0

3.0
Inicio PTAR

2.0
Final PTAR
Cruce

1.0

0.0
2600 2400 2200 2000 1800 1600 1400 1200 1000 800 600 400 200 0
Absicisa por el eje del río Otún (m)

Figura 1-39. Perfiles de profundidad del flujo para diferentes crecientes

Fuente: Hidrased, 2019.

Asimismo, se evaluaron las velocidades de flujo a lo largo del tramo de estudio y se obtuvieron los
siguientes perfiles. En el tramo de la PTAR las velocidades alcanzaron valores hasta de 3,0 m/s para la
creciente de 2,33 años, y valores de hasta 5,0 m/s para la creciente de 100 años de periodo de recurrencia.
Para la creciente mayorada, los valores de velocidad máxima son cercanos a 6,0 m/s. En cualquier caso,
los valores de velocidad son elevados y pueden tener efectos directos sobre la capacidad erosiva del flujo.

1-87

1-87
 9.0
2.33 10 25 50 100 100 mayorada
8.0

7.0
Velocidad (m/s)

6.0

5.0

4.0

3.0
Inicio PTAR

Final PTAR
2.0
Cruce

1.0

0.0
2600 2400 2200 2000 1800 1600 1400 1200 1000 800 600 400 200 0

Absicisa por el eje del río Otún (m)

Figura 1-40. Velocidades del flujo para diferentes crecientes

Fuente: Elaboración propia

De otro lado, en la Figura 1-41 se presenta un perfil con la medida del ancho superficial del flujo. Esto
para ilustrar que los mismos tienden a tener poca variabilidad respecto a cada creciente analizada. Lo
anterior se debe a que el flujo queda confinado casi por completo por la sección natural del río Otún. Por
ejemplo, para la creciente de 2,33 años en el tramo de la PTAR el ancho superficial es del orden de 35 m,
mientras que para la de 100 años es del orden de 50 m. En todos los casos estos valores se incrementan
debido a que se generan ensanchamientos naturales, y para la creciente de 100 años desbordamiento del
flujo sobre la margen izquierda.

1-88

1-88
 Final PTAR
Inicio PTAR
Cruce
140

120
Ancho superficial (m)

100

80

60

40

20
2.33 10 25 50 100 100 mayorada
0
2600 2400 2200 2000 1800 1600 1400 1200 1000 800 600 400 200 0
Absicisa por el eje del río Otún (m)

Figura 1-41. Ancho superficial para diferentes crecientes

Fuente: Hidrased, 2019.

En la Figura 1-42, Figura 1-43 y Figura 1-44 se muestra la mancha de inundación para la creciente de
2,33, 100 y 100 años mayorada de periodo de retorno en el tramo de la PTAR. De acuerdo a los resultados
obtenidos se observa cómo se generan leves desbordamientos, pero en ningún caso se ve afectada la
Planta de Tratamiento de aguas residuales. Adicionalmente en la Tabla 1-44 se presentan los niveles
generados para las crecientes mencionadas.

Tabla 1-44 Niveles máximos estimados en la PTAR

Niveles (msnm) Tr 2,33 años Tr 100 años Tr 100 años Mayorada
Cruce 1125.98 1127.56 1128.28
PTAR 1118.84 1121.25 1121.88

Fuente: Hidrased, 2019.

1-89

1-89
 Figura 1-42. Mancha de inundación para la creciente de 2,33 años de periodo de retorno

Fuente: Hidrased, 2019.

1-90

1-90
 Figura 1-43. Mancha de inundación para la creciente de 100 años de periodo de retorno

Fuente: Hidrased, 2019.

1-91

1-91
 Figura 1-44. Mancha de inundación para la creciente de 100 años de periodo de retorno
mayorada

Fuente: Elaboración propia

Finalmente, para efectos de conocer la afectación de las inundaciones ante eventos extremos se corrieron
algunos caudales adicionales para determinar cuáles eran los que ocasionaban desbordamientos hacia la
infraestructura de la PTAR. Las manchas de inundación para estos escenarios poco probables se
muestran a manera de resumen en la Figura 1-45.

Para un caudal de 1.000 m3/s el flujo queda confinado dentro del canal principal. No existe afectación
desde el punto de vista de extensión de la inundación a la infraestructura proyectada. Para el caudal de
2.000 m3/s, de más baja probabilidad de ocurrencia, hay contacto del flujo sobre la pata del lleno
proyectado. Finalmente, para el caso de 3.500 m3/s el flujo alcanzaría a tocar dos de los tanques
propuestos para el esquema actual.

1-92

1-92
Caudal de 1.000 m3/s

Caudal de 2.000 m3/s

1-93

1-93
 Caudal de 3.500 m3/s
Figura 1-45. Manchas de inundación para otros escenarios poco probables

Fuente: Hidrased, 2019.

Estos últimos casos, superan más de un millón de veces el caudal de 100 años de periodo de retorno (en
términos estadísticos), por lo tanto, desde el punto de vista netamente hidrológico e hidráulico es muy
improbable que el proyecto sea afectado con la dinámica de inundaciones del río Otún en este tramo.

En todo caso, se debe validar con un análisis geomorfológico la posible divagación del lecho, y de ser
aplicable, proponer las respectivas obras de mitigación para mantener el equilibrio del lecho y la estabilidad
de la PTAR, pues las velocidades máximas para las crecientes (2,33 a 100 años) son considerables, y con
alto potencial erosivo.

1.9.6 Conclusiones

La información topográfica y batimétrica a disposición, levantada mediante LiDAR en 2019 para la zona
de construcción de la PTAR, se consideró completa y apropiada para el modelamiento hidráulico de los
niveles y mancha de inundación de la zona de construcción de la planta de tratamiento de aguas
residuales. Esta información fue complementada con la información satelital gratuita del ALOS-PALSAR,
que tiene cubrimiento mundial en resolución de pixel de 12.5 m, apropiada para el análisis hidrológico y
morfométrico realizado.

1-94

1-94
Los estudios previos, en particular los realizados por CARDER y el plan de ordenamiento de la cuenca del
río Otún (POMCA), fueron utilizados como referencia en los componentes de morfometría y en la
recolección y validación de información climatológica. Se utilizaron en este estudio cerca de 28 estaciones
con información histórica de precipitación presentes en la base de datos hidrológica del POMCA del río
Otún como complemento a las solicitadas en este estudio.

Se analizaron los registros mediáticos de los eventos de inundación o de riesgo relacionado con el clima
en la ciudad de Pereira y sus localidades aledañas, encontrándose al menos 18 eventos torrenciales
cubiertos por la prensa o con registros digitales, desde 1979 y abril de 2019, donde se evidencian algunos
efectos directos de afectaciones por deslizamientos producto de lluvias, inundaciones y desbordes, así
como fallas en los sistemas de recolección de aguas pluviales y evacuaciones. La temporada donde más
se encontraron registros y cubrimiento de los eventos se encuentra en los meses de octubre-noviembre,
coincidente con los meses con mayor precipitación promedio multianual.

La cuenca del río Otún a la altura de la ubicación de la PTAR tiene un área de drenaje cerca a los 409 km 2,
con un perímetro de 189 km y con una diferencia entre la cota máxima y mínima de cerca de 3.800 m. El
río Otún hasta el punto de interés tiene una longitud principal cercana a los 72 km; por área la cuenca es
considerada intermedia pequeña, con un índice de alargamiento alto (Ia=12,8), de forma ovalada-
rectangular según si índice de compacidad (Kc=2.6) y con una densidad de drenaje significativa (2,1
km/km2) con un patrón detrítico en sus afluentes. Algunas subcuencas de importancia al río Otún son: la
Q. Barbo (61 km2), la Q. Bosque (32 km 2), Q. Dosquebradas (59 km 2), Q. San José (27 km 2), Q. San Juan
(22 km2) y la Q. Azul (12 km 2).

La información base utilizada para la caracterzación climática de la zona fue compilada a partir de registros
de la red nacional de medición del IDEAM, así como también información proveniente de el plan de
ordenamiento y manejo del río Otún. En total se utilizó información de precipitaciones, temperatura,
evaporación, niveles y caudales según el caso, de las estaciones: Termales, San Isidro, Boquerón, Playa
Rica, Pez Freso, La Laguna y en particular los registros de las estaciones Aeropuerto Matecaña y La
bananera 609. Estas dos últimas estaciones cuentan con información completa de 70 y 47 años de
registros respectivamente.

Se analizaron estadísticamente los registros de las estaciones solicitadas y, utilizando las pruebas de
Dixon y Grubbs, se encontraron valores anómalos en las series que no fueron descartados en ningún caso
debido a que luego de analizarse de forma individual se consideraron registros válidos dentro de los valores
esperados para las variables. La prueba de tendencia Mann-Kendall arrojó que, en general, las estaciones
no presentan una tendencia significativa en el tiempo. Así mismo, las pruebas de Von Neumann fueron
empleados para evaluar la homogeneidad de los registros, encontrándose leves cambios en las medias
de las estaciones, principalmente en los años 1975, 1990, 1998 y 2002, se infiere que este comportamiento
está relacionado con posibles cambios de las curvas de calibración, ajustes de equipos o cambios en la
metodología de pre procesamiento y no afectan de manera significativa los análisis hidrológicos de los
mismos.

El análisis temporal de la precipitación realizado con las estaciones solicitadas muestra un comportamiento
bimodal, con un aumento promedio de los valores en los meses marzo-abril-mayo y en los meses
septiembre-octubre-noviembre intercalándose por periodos de menos precipitación en los meses
restantes. La variabilidad de los registros va desde los 80 mm en Julio y agosto, hasta los 370 mm en los
meses octubre-noviembre, con un promedio mensual de 190mm y un total anual promedio cercano a los
2400 mm. Se utilizó el método de Kriging Ordinario (KO) para la interpolación espacial de la precipitación

1-95

1-95
total anual promedio, observándose 3 zonas de distribución de la precipitación total: (1) cuenca alta con
1000-1400 mm/año, (2) cuenca media: 2200-3000 mm/año, y (3) cuenca baja, zona aledaña a la PTAR,
con variaciones de 1700-2200 mm/año.

El análisis temporal de la temperatura fue realizado a través de los registros de la estación Aeropuerto
Matecaña representativa de la zona de estudio, donde es posible encontrar variaciones de temperatura
promedio en los 21.2°C (octubre-noviembre) hasta 22.3 °C (febrero) y con una diferencia entre las
temperaturas máximas y mínimas promedio de los 12.6 °C en los meses de agosto y febrero. De forma
espacial, la temperatura promedio anual fue estimada utilizando la relación de gradientes con respecto a
la elevación (Barón, 2008), encontrándose valores promedio desde los 0°C (o inferiores) en la parte alta
de la cuenca, hasta los 24°C en la parte baja.

En cuanto al análisis de la evaporación, se evaluaron los métodos de estimación de la evaporación

potencial (ETP) propuestos por Cenicafé, Thornwaithe y la evapotranspiración real (ETR) mediante los
métodos de Budyko, Turc y Coutagne. Los métodos se basan en correlaciones con la temperatura,
precipitación y la elevación. Comparándose con los registros de evaporación de tanque de la estación
Aeropuerto Matecaña, el método de Turc ofreció los mejores resultados con 9% de error relativo promedio
en el valor promedio anual reportado por la estación.

Con el fin de caracterizar la oferta hídrica promedio anual (P-ETR) se utilizaron los registros de la estación
aeropuerto matecaña, encontrándose que los meses de febrero y agosto ofrecen la menor cantidad del
recurso hídrico mientras que los meses de noviembre y octubre los mayores. Se utilizaron además los
resultados espaciales de los campos de precipitación y evapotranspiración para estimar a nivel de cuenca
la oferta hídrica, encontrándose valores promedios multianuales del mismo orden de magnitud que los
ofrecidos por las estaciones de caudal a disposición.

Realizando un análisis comparativo de metodología de estimación de curvas Intensidad-Duración-

Frecuencia se seleccionó la metodología de Pulgarín para estimar la IDF de la zona de estudio, debido a
que es conservadora con respecto a Hershfield y INVIAS, pero no sobrestima los caudales para pequeñas
duraciones comparándola con la curva presentada en el POMCA de 2017.

El comportamiento de caudales medios en la zona de estudio presenta una bimodalidad, es decir, dos
periodos secos y dos periodos lluviosos. Los periodos secos se identificaron entre enero-febrero y julio-
septiembre, mientras que los lluviosos se dan entre abril-mayo y noviembre-diciembre. El caudal medio de
la bananera es de aproximadamente 12.6 m3/s y en la PTAR es de 18,6 m3/s.

Los caudales máximos estimados con el modelo lluvia-escorrentía se seleccionaron para la construcción
del modelo, dado que son mas conservadores y su valor, aunque es mayor a los estimados en análisis de
frecuencia, no generan diferencias significativas en los niveles estimados.

Las profundidades estimadas mediante el modelo hidráulico tienden a ser muy uniformes para cada
periodo de retorno, es decir que no tienen mucha variabilidad.

Las velocidades de flujo alcanzaron valores hasta de 3,0 m/s para la creciente de 2,33 años, y valores de
hasta 5,0 m/s para la creciente de 100 años de periodo de recurrencia. Para la creciente mayorada, los
valores de velocidad máxima son cercanos a 6,0 m/s. En cualquier caso, los valores de velocidad son
elevados y pueden tener efectos directos sobre la capacidad erosiva del flujo.

1-96

1-96
El flujo queda confinado casi por completo por la sección natural del río Otún para todas las crecientes
evaluadas, únicamente se generan leves desbordamientos, pero en ningún caso se ve afectada la Planta
de Tratamiento de aguas residuales.

Realizando un análisis de desbordamiento se determinó que para un caudal de 1.000 m3/s el flujo queda
confinado dentro del canal principal. Este caudal supera superan en términos estadísticos más de un millón
de veces el caudal de 100 años de periodo de retorno, por lo tanto, desde el punto de vista netamente
hidrológico e hidráulico es muy improbable que el proyecto sea afectado con la dinámica de inundaciones
del río Otún en este tramo.

Se recomienda realizar un análisis geomorfológico la posible divagación del lecho, y de ser aplicable,
proponer las respectivas obras de mitigación para mantener el equilibrio del lecho y la estabilidad de la
PTAR, pues las velocidades máximas para las crecientes son considerables, y con alto potencial erosivo.

Las manchas de inundación presentadas cumplen con los requisitos hidrológicos e hidráulicos de definición
de ronda hídrica indicados por la CARDER, sin tener en cuenta componentes de sedimentos o nutrientes.

1.10 Geomorfología.
El área de estudio se localiza hacia la vertiente Occidental de la Cordillera Central en un ambiente fluvio-
volcánico y denudacional que predomina sobre un ambiente estructural. La acumulación de productos
volcánicos como lavas, cenizas, bloques y detritos, los cuales han sido transportados, mezclados y re
trabajados por ríos como el Otún y en otros casos provienen del derretimiento del hielo o nieve durante la
erupción volcánica, que han suavizado la topografía del paisaje.

En buena medida, la geomorfología de la zona está definida desde la cartografía geológica realizada por
INGEOMINAS, desde 1984, sin embargo, con posteridad se han trazado algunas directrices para la
elaboración de los mapas geomorfológicos y por esta razón, en este trabajo se asimila en lo posible a la
subdivisión del paisaje contenida en el POMCA actualizado (CONSORCIO ORDENAMIENTO CUENCA
RÍO OTÚN, 2017)

En el aparte de Geomorfología del POMCA del río Otún en escala 1:25.000 (Figura 1-46) y en el perfil de
la Ilustración 3-1, se aprecia que el Abanico de Pereira (clasificado como Glacis de Erosión), abarca toda
la zona de influencia del proyecto. Solamente en las márgenes inmediatas del río, se aprecia una franja
estrecha de color naranja clasificada como Llanura Aluvial, rodeada de algunas áreas de forma alargada
identificadas con color rojo y descritas como Terrazas Fluviales, que corresponde justamente con las zonas
“planas” dentro de las cuales se localizarán algunas de las estructuras de la PTAR. Las unidades
planteadas para la geomorfología regional de la zona de estudio, se especifican a continuación de acuerdo
con el estudio citado, al igual que su descripción particular, de nuevo, se advierte que este es un estudio
a escala 1:25.000 y necesariamente presenta diferencias con el 1:5.000 que corresponde a este trabajo,
comenzando por la precisión de la base topográfica, lo que incide hasta en la forma del cauce del río.

1-97

1-97
 Figura 1-46. Geomorfología regional
Fuente: Elaboración propia, 2019. Modificado (CONSORCIO ORDENAMIENTO CUENCA RÍO OTÚN,
2017).

Ilustración 3-1. Perfil geomorfológico de la cuenca del río Otún

Fuente: CONSORCIO ORDENAMIENTO CUENCA RÍO OTÚN, 2017.

1-98

1-98
Glacis de Erosión. Comprende todo el abanico de Pereira como una unidad de topografía muy suave en
la que se encuentran montículos a lomeríos de baja altura, en general menores de 50 m entre topes y
valles, con lomos alargados en dirección W, de topes redondeados y valles estrechos en forma de artesa.
Las vertientes tienen pendientes bajas a moderadas, entre 10° y 20° de inclinación. El paisaje erosional
original existente antes del Holoceno se encuentra suavizado por depósitos de ceniza, algunos de ellos
más antiguos de esta época, tan viejos que se encuentran altamente meteorizados, lo que puede indicar
hasta centenares de miles de años en los cuales se habría incisado el Abanico.

Llanura de inundación. Este nombre corresponde a las zonas donde se ha denudado el Abanico y
comprende no sólo las vertientes del mismo hacia el río Otún, sino también depósitos de vertiente de
diferentes tipos localizados a lo largo de ellas. Se caracteriza por pendientes moderadas a altas, del orden
de 20 a 45°, con vertientes rectas, en general cortas, de menos de 600 m y se pueden extender desde el
borde superior del escarpe hasta el río. También contiene escarpes subverticales en los cuales se puede
apreciar la composición del abanico, como en el escarpe del norte del área urbana de Pereira (Figura 1-
47).

Terrazas Fluviales. En esta categoría se enuncian en el POMCA, depósitos de diversos tamaños y

composición, de origen aluvial que se han depositado en diferentes épocas en las márgenes del río Otún.
El estudio referido no separa niveles ni determina edades para ellas.

“Como aspectos prominentes a nivel fisiográfico, se destaca en la parte alta de la cuenca, la presencia de
relieves contrastantes, derivados de la evolución de la actividad volcánica, una franja de relevancia a nivel
de evolución y rasgos tectónicos entre las Fallas Palestina al oriente y parte alta de la cuenca, la parte
central cruzada por trazas asociadas a la Falla Silvia – Pijao, con cobertura y enmascaramiento de flujos
de lodo y cenizas volcánicas, mientras en la porción distal de la Cuenca, se presenta un prominente
piedemonte de compleja interpretación evolutiva, pero asociada en general a grandes aportes derivados
de actividad explosiva de las fuentes volcánicas, durante el Plio-Pleistoceno inferior y medio, que
igualmente representa una cubeta sedimentaria delimitada por fallas, donde la Falla Quebradanueva
marca su límite oeste .”

Una visualización aproximada de la geomorfología la proporciona el modelo digital de terreno (Figura 1-

47) y en la imagen tomada de Google Earth (Ilustración 3-2), donde se destaca el área dónde se localiza
el proyecto de la PTAR El Paraíso. Se aprecia claramente la ubicación dentro de una topografía suave que
es la delimitada por el Abanico de Pereira, el cual, en el sitio, se encuentra incisado por el trabajo erosivo
del río Otún en las últimas decenas de miles de años.

1-99

1-99
 Figura 1-47. Modelo digital de terreno del sector occidental del río Otún
Fuente: Elaboración propia, 2019.

De acuerdo con CONSORCIO ORDENAMIENTO CUENCA RÍO OTÚN, 2017, debido a la forma de la
cuenca, esta presenta una probabilidad media en cuanto a crecientes repentinas o avenidas torrenciales,
aunque este tipo de eventos se relaciona más a la parte alta de la cuenca, como se presentará en el
capítulo de amenazas de este informe.

1-100

1-100
 Ilustración 3-2. Fisiografía del contexto del área de influencia

Elaboración propia, con base en imagen de Google Earth, 2019.

Ilustración 3-3. Área de Influencia

Elaboración propia, con base en imagen de Google Earth, 2019.

1-101

1-101
1.10.2 Subunidades geomorfológicas – Geomorfología local.

Mediante el análisis de información secundaria como las memorias y planchas del Servicio Geológico
Colombiano (SGC) a escala 1:100.000, el modelo de elevación digital del terreno, fotointerpretación,
puntos de control en campo y análisis de la información, para el objetivo de este proyecto se llegó a la
definición de las subunidades geomorfológicas y a la identificación de los procesos morfodinámicos a
escala 1:5.000 (Figura 1-48).

La geomorfología se presenta de acuerdo con la Propuesta de Estandarización de la Cartografía La

geomorfología se presenta de acuerdo con la Propuesta de Estandarización de la Cartografía
Geomorfológica en Colombia (Servicio Geológico Colombiano, 2011), alcanzando el nivel de detalle de
subunidades que permite desarrollar información necesaria para el proyecto y el análisis de la Gestión del
Riesgo.

Las subunidades geomorfológicas identificadas durante este trabajo corresponden a una subdivisión de
las Unidades Geomorfológicas enunciadas por el POMCA, definiéndose esta última como una geoforma
individual genéticamente homogénea. Las subunidades están determinadas fundamentalmente por los
contrastes morfológicos y morfométricos, que relacionan el tipo de material o la disposición estructural de
los mismos, con la correspondiente topografía del terreno.

1-102

1-102
 Figura 1-48. Subunidades geomorfológicas
Fuente: Elaboración propia, 2019. Modificado (CONSORCIO ORDENAMIENTO CUENCA RÍO OTÚN,
2017).

En el área de influencia las subunidades que configuran el relieve son producto de un ambiente
denudacional y fluvial, el ambiente morfogenético hace referencia a las condiciones físicas, químicas,
bióticas y climáticas bajo las cuales se generaron las geoformas, identificando su origen ya sea endógeno
o exógeno que dan lugar al desarrollo y evolución de estas geoformas.

1.10.3 Abanico de Pereira (Dap).

Esta unidad se clasifica con un ambiente morfogenético Denudacional (erosional), aunque su origen es
fluvio-volcánico porque en el presente su superficie se encuentra sometida a procesos tales como
meteorización química, erosión por escorrentía y movimientos en masa. Estos procesos han modelado
sobre el abanico original un relieve de colinas suaves que se presentan en el sector norte del lote El Paraíso
y en la parte sur de la vía Pereira – Cartago del San Cayetano.

1-103

1-103
Consiste en una superficie suave, con montículos largos, hasta de 2 km de longitud, muy bajos, de menos
de 10 m de altura, orientados con dirección general SW y pendientes menores de 5°, que sólo son
excedidas en las depresiones formadas por algunas cañadas como resultado de menos de 10 m. de
incisión que presenta la geoforma actual (Figura 1-49). En toda su extensión, la topografía ha sido
suavizada por numerosos episodios de caída de cenizas volcánicas que conforman también la parte
superior de la secuencia fluvio-volcánica original del Abanico.

Figura 1-49. Perfil de terraza Abanico de Pereira (Qap) margen derecha del río Otún

Fuente: Elaboración propia, 2019.

Al interior de la unidad son prácticamente inexistentes los procesos morfodinámicos, bien sea por
movimientos en masa o profundización o socavación de cauces del orden 1 ó 2. En cauces mayores sí se
presentan estos fenómenos, intensificados por el efecto antrópico de urbanización o mal manejo de
cultivos, pero éstos no se presentan en el área de influencia directa del Proyecto

1.10.4 Abanico incisado (Dapi).

La profundización y erosión lateral del río Otún, activa desde el final de la formación del Abanico, ha
entallado en el sitio un escarpe de aproximadamente 80 a 100 m de profundidad que contiene vertientes
rectas, cortas, de menos de 500 m de longitud con pendientes desde moderadas, a partir del 7% como el
filo principal del lote Paraíso, hasta sub-verticales como en los escarpes del extremo noroccidental del lote
San Cayetano y en la vertiente superior del mismo. Las vertientes muestran escasos rasgos de incisión
interna (Fotografía 3-1).

En varios sectores del escarpe el proceso de formación no ha sido exclusivamente erosivo, sino que han
ocurrido pulsos temporales de depositación aluvial, que se podrían agrupar como componentes
(elementos) de esta geoforma pero que por claridad para el lector se describen en los numerales siguientes
como geoformas de depósitos aluviales y de vertiente.

1-104

1-104
 Fotografía 3-1. Vista parcial del lote Valparaíso con el Abanico de Pereira con escarpe y
Terrazas aluviales
Fuente: Elaboración propia, 2019.

Esta geoforma se encuentra cubierta superficialmente por unos 5 m de depósitos de cenizas de caída
piroclástica en estado de meteorización avanzada, lo que parece indicar varios miles a decenas de miles
de años de formación. Tomando como criterio su morfología, la edad del escarpe es avanzada, desde
hace más de 40.000 años como lo afirman algunos estudios realizados en el abanico, que indican más de
estos 40.000 años para algunos de los depósitos más superficiales que lo cubren y puedan ser coetáneos
con el principio de la incisión.

1.10.5 Depósito de vertiente (Qdv).

En el área de influencia del proyecto se destacan dos depósitos de vertiente que han suavizado la
topografía del escarpe y que no sólo se presentan en este sitio, sino que son típicos del cañón del Otún
en muchos tramos a lo largo del río. En la margen derecha se encuentra un depósito de carácter coluvial,
mientras que en la margen izquierda se encuentra un depósito de caída de bloques. Ambos depósitos se
formaron durante el proceso de profundización del río, pero su edad es difícil de estimar por la carencia de
estudios sistemáticos a lo largo del río y de la dificultad de obtener materiales que permitan establecer su
edad por métodos radiométricos.

1-105

1-105
1.10.5.1 Depósito de vertiente (Caída de bloques) del lote San Cayetano

En la mitad inferior de la vertiente sur del escarpe en el lote San Cayetano, se puede encontrar un depósito
con una pendiente aproximada de 20 a 35°, localizado aproximadamente entre las cotas 1140 y 1170,
compuesto en superficie por bloques redondeados de tamaño métrico a decimétrico de andesitas,
dispersos sobre la vertiente, lo que indica un proceso de exhumación por erosión superficial prolongada
por lo menos a lo largo de cientos de años, como lo evidencia su meteorización superficial que le da un
aspecto rugoso. El único afloramiento encontrado permite ver una serie de bloques y guijarros
redondeados de andesitas, flotantes en una matriz areno arcillosa de color pardo, cohesiva que le da
soporte al conjunto del depósito. Para un mayor detalle textural y composicionalmente (ver
ANEXOS_MEDIO ABIÓTICO_ 1 Estudio Geotécnico).

La génesis del depósito es exclusivamente por lavado de la matriz de un estrato superior compuesto de
bloques rocosos que paulatinamente se sueltan y ruedan pendiente abajo para alcanzar un sitio de reposo,
mecanismo que forma un depósito de caída que se conoce con el nombre de “Talus”.

Fotografía 3-2. Depósito de vertiente por caída de bloques en el lote San Cayetano
Fuente: Elaboración propia, 2019.

Sobre esta unidad se proyecta la construcción de los sistemas sedimentadores iniciales de la PTAR, por
lo que la investigación geotécnica debe tener un cuidado especial en las marcadas diferencias que
sepresentan en las propiedades físicas de estos suelos en distancias cortas, a veces de metros, además
de la dificultad constructiva que representan las fundaciones en este tipo de depósitos.

1-106

1-106
1.10.5.2 Depósito Coluvial del lote El Paraíso

En este sector del Proyecto PTAR, hacia la parte baja de la vertiente se encuentra una superficie de
inclinación suave y recta, de menos de 7° de inclinación, donde a lo largo de miles de años se han
depositado materiales finos arrastrados por aguas superficiales de escorrentía y depositadas en la parte
baja del cañón del río, específicamente sobre la terraza T1 que se encuentra bordeando el río en el sitio,
en la Fotografía 3-1 se delimita como depósito de vertiente.

Un afloramiento pobre visible en el corte de la vía de acceso a la parte baja del lote muestra una masa
limoarcillosa de color pardo claro compacta y cohesiva con guijarros dispersos sobre ella.

1.10.6 Terrazas aluviales.

La incisión permanente del río Otún sobre el Abanico ha logrado a través del tiempo formar un cañón
relativamente estrecho, de menos de 1 km de ancho, proceso que ha permitido la acumulación de
materiales de su lecho. Al ocurrir cambios en su trayectoria y nivel a lo largo del tiempo, el río abandona
el lecho anteriormente activo, dejando en sus bordes, superficies planas compuestas por bloques y
guijarros, que responden al estilo particular de la corriente en los diversos tramos y épocas de depositación
que se conocen con el nombre genérico de terrazas aluviales.

Con el proceso de incisión se han formado varios niveles de terrazas, depositados en las márgenes del
Otún que se encuentran a diferentes alturas sobre el lecho actual, siendo las más altas de mayor edad.
Este sistema es más o menos constante en los diferentes tramos del río y corresponde a la dinámica fluvial
particular en los diferentes tramos del mismo. En la zona de influencia del lote de la PTAR se encontraron
cinco niveles de terrazas que numeradas de más baja a más alta se describen a continuación.

1.10.6.1 Terraza cero (T0)

Este elemento geomorfológico se encuentra en la margen del río como acumulaciones de materiales
recientes que alcanzan una altura hasta de 1.5 – 1.7 m sobre el nivel de aguas medias, con una
composición similar en tamaño y composición que el lecho actual. Se caracteriza por estar constituida por
bloques métricos y decimétricos con cantidades subordinadas de guijarros, en mayor proporción de
composición andesítica, casi todos con alta redondez y mediana a alta esfericidad, con excepción de una
menor proporción de bloques de rocas metamórficas que se encuentran subredondeados, tal como se
aprecia en la Fotografía 3-3.

Una característica particular de este nivel es la acumulación de basuras especialmente plásticos en su

borde que muestran el nivel de máximos caudales, al igual que bolsas y trapos que durante los episodios
de crecientes quedan colgados de las ramas de los árboles indicando así los niveles máximos alcanzados
en los últimos años. Este nivel corresponde de manera aproximada a la creciente de 2.33 años calculada
por métodos hidrológicos, aunque a veces puede indicar recurrencias mayores hasta de 5 a 7 años, dada
la prolongada duración de muchos de éstos elementos plásticos en la superficie del terreno y la variabilidad
climática en la zona andina.

El límite entre la T0 y la T1 es un escarpe erosivo vertical, lo que indica que las crecientes máximas recientes
a subrecientes del río (de los últimos cientos a menos de mil años), han erosionado el depósito de la T 1,
pero no han alcanzado a depositar materiales encima de él, lo que permite concluir que cerca de la cota

1-107

1-107
de T1 está el valor del caudal máximo esperado en un período de un poco menos de 1.000 años que es el
tiempo necesario para redondear un escarpe en el medio andino colombiano (Parra, 2016).

Fotografía 3-3. Aspecto de la terraza T0, en el borde del río

Fuente: Elaboración propia, 2019. (Se aprecian plásticos atorados en las rocas y el escarpe vertical del
depósito de la T1)

1.10.6.2 Terraza uno (T1)

Este nivel se encuentra aproximadamente a 1.5 m de altura de la terraza T 0 en el sitio del proyecto, es
decir, a 3.2 m de altura sobre el nivel de aguas medias del río. Aflora sobre ambos lotes del proyecto
formado un sistema pareado entre ellos (depósitos de igual altura a ambos lados del cauce ver . Se
compone de bloques, cantos y guijarros redondeados, en su mayoría de andesitas, con contactos
tangenciales entre sí, compactados por matriz areno-limosa que no es lavada de manera fácil por la
corriente durante los episodios de crecientes importantes. La disposición de los clastos, es decir, la
estructura clasto-soportada indica unas condiciones de formación donde los clastos mayores se
depositaban uno sobre otro como carga de arrastre y posteriormente los intersticios se llenaban con
materiales más finos de textura areno - limosa, como se muestra en la Fotografía 3-4.

1-108

1-108
 Figura 1-50. Perfil de terrzas QT1 y QT2 margen izquierda del río Otún

Fuente: Elaboración propia, 2019.

Fotografía 3-4. Composición de los depósitos aluviales en el sitio del proyecto

Fuente: Elaboración propia, 2019.

Las condiciones de depositación ya descritas son evidencia de formación durante crecientes hidrológicas
normales, donde la cantidad de sólidos que llevaban las corrientes que los transportaron era relativamente
baja y no se puede calificar como formada por un flujo torrencial, que se caracteriza por cantidades
importantes de sólidos, del orden de 20 a 30% o más. En el sitio, la superficie de la terraza se encuentra

1-109

1-109
con una cubierta vegetal compuesta totalmente por guadua, por lo que su superficie está tapizada por
hojas caídas de esta especie que inhiben el crecimiento de sotobosque, pero permiten apreciar fácilmente
su superficie plana, con inclinación paralela al río (ver Fotografía 3-5).

Con el fin de hacer un estimativo preliminar de su edad de formación mediante la medida del espesor de
su horizonte edáfico A, se perforó con un auger manual desde la superficie y fue posible alcanzar más de
25 cm, lo que puede indicar una edad de formación de al menos decenas de miles de años. En el caso de
que en tiempos históricos o posteriores a la formación de dicho horizonte A, una creciente hubiese
superado este nivel, se supone que o bien lo habría erosionado, o se habría depositado algún estrato, así
fuera delgado, de limos o arenas sobre él. Si el evento hubiese sido torrencial, se encontraría sobre el A,
una capa superficial de bloques y guijarros con disposición caótica y matriz arenoarcillosa o una barra
torrencial que consiste de una serie de bloques alineados de forma paralela a la corriente que los
transportó.

Fotografía 3-5. Superficie de la Terraza T1, con auger que muestra el espesor del
horizonte A > 25 cm
Fuente: Elaboración propia, 2019.

1.10.6.3 Terraza dos (T2)

Este nivel de terraza, que es el de mayor extensión, se puede observar en el sector occidental del lote San
Cayetano y se encuentra aproximadamente a 1.7 metros de altura sobre el nivel de la Terraza T1, es decir,
a 4.9 m sobre el nivel de aguas medias del río Otún (Fotografía 3-6), con un ancho desde 15 hasta 30 m
y una longitud mayor de 800 m. En el sitio, al igual que la T1, se encuentra con una cubierta vegetal de
guadua y naranjos tecnificados, por lo que su superficie está cubierta de hojas caídas de estas especies
que limitan con su sombra el crecimiento de especies rastreras, pero permiten apreciar fácilmente su
superficie plana, con inclinación suave, paralela al río, tal como se aprecia en la Fotografía 3-6.

1-110

1-110
El espesor del horizonte A, que alcanza más de 40 cm, de acuerdo con una perforación con auger manual
realizada en su superficie, parece evidenciar una edad del orden de unas pocas decenas de miles de años
y tampoco se encuentran evidencias de haber sido alcanzada por flujos torrenciales en este lapso de
tiempo.

Fotografía 3-6. Disposición paralela al río de la Terraza T2, lote San Cayetano
Fuente: Elaboración propia, 2019.

1.10.6.4 Terraza tres (T3).

Este rasgo fluvial se encuentra únicamente en el lote Valparaíso, a unos 18 - 20 m por encima del cauce
actual como se puede observar en la Fotografía 3-1 y se aprecia no sólo con inclinación paralela al río sino
también con declive hacia él, lo que indica una edad mucho mayor que las otras terrazas más bajas ya
descritas, tal vez del orden de centenares de miles de años. Presenta unas dimensiones superficiales de
unos 100 m de longitud por 40 m de ancho, pero carece de afloramientos superficiales y su composición
solamente se podrá conocer mediante apiques que se deben realizar al momento de realizar los estudios
de suelos. Debido a su tamaño y espesor menores y a que durante la etapa de construcción va a ser
cortada en su totalidad, no es un elemento de importancia relevante para este Proyecto.

1.10.6.5 Terraza cuatro (T4).

En la parte central del lote Valparaíso, a unos 25 m de altura sobre el nivel del río (Fotografía 3-1 y Figura
1-51), se puede observar un peldaño casi plano de aproximadamente 50x30m de extensión, con bordes
completamente redondeados que indican una edad avanzada de incisión, del orden de miles de años. No
es posible conocer su composición ni textura, la que debe ser objeto de análisis durante la campaña del
estudio de suelos, pero igual que la anterior, por sus dimensiones y posición respecto a la infraestructura
a construir, es un elemento que debe desaparecer al realizar los cortes respectivos.

1-111

1-111
En la Figura 1-52 se muestra una aproximación a la geomorfología utilizando un Modelo Digital de Terreno
que ilustra parcialmente las sub-unidades, donde las terrazas T 3 y T4 se confunden con los depósitos
coluviales, debido a la depositación de cenizas y aporte coluvial incipiente en las primeras, factor que les
proporciona un declive hacia el río.

Figura 1-51. Perfil de terrza QT4 margen derecha del río Otún

Fuente: Elaboración propia, 2019.

Figura 1-52. Geoformas área de estudio

Fuente: Elaboración propia, 2019.

1-112

1-112
1.11 Avenidas Torrenciales.
Las avenidas torrenciales son crecientes excepcionales que se presentan en algunos cauces y superan
en decenas y hasta centenares de veces los caudales normales de la corriente respectiva. Estos eventos
se asocian en el área andina de Colombia a episodios extremos de lluvias que generan enjambres de
movimientos en masa; por bloqueo de cauces por movimientos en masa únicos que forman diques, donde
su posterior ruptura genera un gran caudal. En última instancia, por enjambres de movimientos en masa
cosísmicos que llegan a un cauce principal como flujos simultáneos, aumentando así su cantidad de
sólidos y, en consecuencia, su potencial destructivo.

No existe un límite preciso entre avenida torrencial y una creciente de un cauce puesto que a medida que
crece el caudal, la corriente gana velocidad y capacidad de arrastrar sólidos, que una vez incorporados al
fluido lo convierten en un medio más denso y destructivo. En general una creciente de un cauce con
contenido menor de tal vez un 20% de sólidos se puede considerar como “normal”, en tanto que con el
aumento de sólidos pasa a recibir el calificativo de torrencial. El término avenida torrencial es correcto en
el idioma español, pero tiene otros nombres dependiendo de los países, en Perú es guaico, en España
torrente y en Centroamérica crecida; en Colombia se les llama coloquialmente avalancha (que es otro
fenómeno geológico diferente, donde el aire es el medio de transporte como la nieve o un caso particular
de erupción volcánica).

Bajo un marco conceptual de la nomenclatura internacional, en el idioma inglés, el término “debris flow” se
utiliza principalmente para describir materiales transportados por un movimiento en masa que se depositan
cerca de su zona de arranque, lo que en nuestro medio recibe el nombre de depósito de vertiente, pero
también se usa para referirse a flujos que se mueven confinados a un canal donde abundan bloques de
roca y gravas, con proporciones menores de arenas y finos. El término “mudflow” se usa para referirse a
un material viscoso compuesto por gravas, arenas y finos que se mueve por el fondo de un cauce al cual
se le sumarían los bloques de rocas y sería más o menos equivalente a lo que significa el término avenida
torrencial en Colombia. En este trabajo no se pretende hacer una discusión académica de los términos y
por lo tanto se asume, la definición que encabeza este numeral.

Cada evento torrencial que transita por un cauce, deja en los sectores de menor pendiente un depósito
que es proporcional al tamaño del flujo transportado, el cual muestra unas características típicas (Tabla 1-
45), que permiten diferenciarlo de otros depósitos originados por las crecientes “normales” asociadas al
cauce (Costa, 1988). De este modo, cauces con grandes y frecuentes avenidas torrenciales forman
abanicos o terrazas aluviales de tamaño mayor que otros de menor actividad, todos ellos fácilmente
identificables mediante interpretación de imágenes remotas, dependiendo de su escala de estudio. En el
caso de cauces pequeños, con cuencas de menos de 10 km², difícilmente transitan avenidas torrenciales,
en razón a que su bajo caudal no es suficiente para transportar un volumen grande de detritos y, en
consecuencia, se forman preferentemente depósitos semejantes a los conocidos en la literatura
geomorfológica como debris flows. Estos últimos son igualmente peligrosos, pero si alcanzan a llegar a su
confluencia con otro cauce, persiste su carácter violento a lo largo de centenares de metros y aún varios
kilómetros, aunque su contenido de sólidos se va depositando aguas abajo hasta perder energía y, en
consecuencia, capacidad de destrucción.
Tabla 1-45 Diferencias entre depósitos aluviales y torrenciales

1-113

1-113
 Depósito torrencial Depósito aluvial
Estructura caótica Estructura ordenada
Clastos flotantes en matriz Clastos con contactos tangenciales entre sí
Gradación inversa Gradación normal
Clastos angulares abundantes Clastos angulares ocasionales
Estratificación irregular Estratificación plana a lenticular
Imbricación de clastos elongados
Disposición irregular de clastos elongados
Textura homogénea
Textura heterogénea

Fuente: Elaboración propia, 2019. Modificado de Costa, 1988.

El POMCA del río Otún del año 2017 indica que en su cuenca existen una serie de depósitos torrenciales
antiguos los cuales requieren una identificación detallada incluyendo dataciones adicionales a las
efectuadas por la Carder, que permitan establecer las zonas realmente involucradas con eventos antiguos
y diferenciables de los depósitos fluviovolcanicos existentes en la zona y concluye que la mayoría de los
eventos torrenciales son ocasionados por precipitaciones fuertes hacia la parte alta de la cuenca y en
específico hacia el nacimiento del río Otún. En combinación por las altas pendientes y las condiciones del
suelo, siendo estos muy meteorizados y por tanto, deleznables, generan caídas de material rocoso sobre
el cauce del rio y sobre afluentes directos tales como el rio Barbo que por eventos históricos han generado
grandes avenidas torrenciales que caen directo al rio Otún generando taponamiento y desabastecimiento
de servicios públicos para el municipio de Pereira 9. Como resultado del análisis presenta el mapa de
amenaza por avenida torrencial, el cual en el área de influencia se clasifican niveles de amenaza alta,
media y baja (Figura 1-53).

9 CONSORCIO ORDENAMIENTO CUENCA RÍO OTÚN. POMCA Ajuste plan de ordenación y manejo de la cuenca
del río Otún. Fase de diagnóstico. Pereira, 2017

1-114

1-114
Fuente: Elaboración propia, 2019. Modificado de (CONSORCIO ORDENAMIENTO CUENCA RÍO OTÚN,
2017).
Figura 1-53. Amenaza por avenida torrencial

El estado de estabilidad actual que presentan las márgenes del río, indican que históricamente sólo las
terrazas más bajas serían afectables, sin efectos mayores de socavación ni depositación como se muestra
en la caracterización del componente Geomorfológico, pues no ha habido afectaciones significativas de
este tipo en el sitio en cientos de años.

Con el fin de determinar la amenaza potencial de este fenómeno en el área de estudio y su posible
afectación para las obras del proyecto, se presenta un estudio detallado de la torrencialidad del río Otún
en el área de interés. Básicamente se trata de responder a las siguientes preguntas: ¿Por qué ocurren?,
qué tamaño tienen?, ¿dónde se originan?, ¿qué tan rápidas son? con qué frecuencia se presentan? y,
por último, ¿qué tanto podrían afectar al lote donde se encuentra el proyecto PTAR Pereira?

La metodología de evaluación de la torrencialidad utilizada en este trabajo es directa, es decir, se basa en

la interpretación de los depósitos asociados a cada cuenca, que fueron estudiados en sitios escogidos
previamente por interpretación de imágenes para ser analizados en detalle durante el trabajo de campo,
lo que permite la caracterización de los eventos de caudales importantes que han transitado por el cauce
respectivo a lo largo de miles e incluso en hasta millones de años. Para inferir las crecientes futuras y de

1-115

1-115
esta manera establecer la Amenaza por este tipo de eventos, se determinan las edades de los depósitos
más recientes de acuerdo con los conceptos de Schneuwly-Bollschweiler (2013), para lo cual se utilizan
métodos directos de correlación como son la existencia de poblaciones de líquenes, espesores de los
horizontes A y B del suelo y la longevidad de las cicatrices de movimientos en masa identificadas en las
vertientes que no sobrepasen miles de años, de acuerdo con los criterios de asignación de edades de
Parra (2016)10.

1.11.2 Torrencialidad histórica.

Esta disciplina se apoya también en los registros históricos, los cuales pueden ser valiosos en la medida
que ellos se encuentren bien compilados o sean parte de una base de datos que corrobore la ocurrencia
de los hechos en el pasado y proporcionen información acerca de un período de retorno aproximado, de
manera semejante a lo que se estudia para evaluar la amenaza sísmica. Su importancia también radica
en que un acontecimiento histórico documentado, permite convencer a la población afectable de la
posibilidad real de que se presente el suceso y de esta manera aceptar con facilidad las medidas de
protección que se tomen al respecto. La documentación fue posible con ciertas limitaciones, en razón a
que no hay un inventario o catálogo local de los desastres naturales de Pereira con suficiente detalle que
permita evaluar la magnitud de los eventos. En última instancia, se tomó como base la información del
catálogo desinventar.org, que parece ser el más completo con la información de Colombia. Cabe anotar
que la colonización de la zona de Pereira sólo se dio en la mitad del siglo XIX y, en consecuencia,
solamente existen un poco más de 150 años de historia, en contraste con la frecuencia de las avenidas
torrenciales que a veces alcanza centenares y aún miles de años, lo que significa que es muy probable
que ninguna de ellas se encuentra en los relatos históricos.

1.11.3 Estudios anteriores.

Existen sólo unos pocos trabajos que abordan específicamente el tema de las crecientes del Otún con
excepción de los informes disponibles en redhidro.org; El POMCA Otún, a pesar de ser un trabajo reciente,
aborda el tema con una metodología cualitativa que no conduce a resultados que permitan tomar
decisiones de prevención o manejo de contingencias, del cual se toman las apreciaciones de los dos
siguientes párrafos.

“La cuenca del Río Otún hasta el sitio de aforo El Cedral, según el parámetro Coeficiente de Compacidad
igual a 1.40, se encuentra en el rango de oval redonda u oval oblonga, que indica que la probabilidad de
que se presenten crecientes o avenidas torrenciales en este punto es media, ya que la concentración del
agua que llega a la cuenca va a tener relativamente los mismos tiempos de concentración, es decir que el
agua que llega a la corriente principal, producto de las lluvias, lo va hacer casi en el mismo intervalo de
tiempo.”

“El parámetro morfométricos Factor de Forma igual a 0.25, indica que la cuenca, hasta éste punto,
presenta relativamente una baja concentración de agua, lo que se traduce en una probabilidad media de
presentar crecientes. De otra parte, su forma es alargada hecho que la puede hacer menos susceptible a

10
PARRA, Eduardo, 2016. Determinacion de edades en depósitos de vertiente y avenidas torrenciales en Colombia.
Mem. XVI Congreso Colombiano de Geología, Santa Marta, 23 p.

1-116

1-116
crecidas por efecto diferencial entre los puntos de la divisoria y un punto central. Esta forma también hace
pensar que tiene baja capacidad de captación de las aguas lluvias.”

En el estudio MUNICIPIO DE PEREIRA –CARDER (2004)11, se cita “...en la cuenca del Río Otún han
ocurrido y seguirán ocurriendo avalanchas y por esto hay que tomar algunas medidas de control o
seguridad sobre la población... lo anterior se evidenció más recientemente con la reciente del Río Otún de
marzo del año 2002”. En el mismo estudio, se presenta el mapa de isoyetas, donde se aprecia un máximo
pluviométrico de 2500 mm anuales de lluvia entre La Florida y El Cedral que aparentemente se prolonga
en jurisdicción del municipio de Dosquebradas. Esta distribución se relaciona en parte con el origen de
algunas crecientes.

La Dirección de Gestión del Riesgo de Pereira –DIGER- (2019)12 hace un diagnóstico donde concluye:
“Si bien las avenidas torrenciales son uno de los eventos que mayores afectaciones ha causado en la
población del municipio, en la presente investigación se encontraron pocos estudios, que si bien han
aportado al conocimiento de la amenaza, no se han consolidado en medidas concretas de reducción y
manejo de desastres.”

Como se puede notar, no hay un estudio específico acerca de la torrencialidad potencial del río Otún, por
lo que en este trabajo se aborda esta temática sobre el total de la cuenca, con el fin de evaluar la amenaza
que pueda tener el sitio del proyecto PTAR. Es preciso advertir que este estudio no incluye la evaluación
de la amenaza por flujos de escombros producidos por grandes movimientos en masa que se puedan
producir en la cuenca, los que no parecen haberse producido al menos en los últimos 10.000 años, de
acuerdo con la fotointerpretación general del área y verificación parcial en campo.

1.11.4 Antecedentes históricos

El análisis de los eventos históricos que han ocurrido en la cuenca es básico para conocer en parte los
intervalos de recurrencia de los eventos y su magnitud o severidad con la que han afectado al entorno
físico y social. En este caso se consultaron los archivos de Desinventar.org con el fin de conocer las
inundaciones y avenidas torrenciales que se hayan presentado en la cuenca del Otún. Esta base de datos,
aunque carece de muchos detalles de los eventos, tales como duración de lluvias y alturas de inundación,
contiene casi con toda seguridad los eventos más relevantes que han ocurrido desde 1920 por su
afectación. Se pudo verificar durante su consulta que algunos de los sucesos no aparecen en el catálogo,
pero no se puede esperar que esas bases de datos tengan la totalidad de los hechos importantes en la
cuenca. Se encontraron documentos sobre casos de eventos que no aparecen en la base y varios de ellos
se describen brevemente citando la fuente respectiva, dada la importancia de los datos contenidos en
ellos.

Del catálogo Desinventar.org, se filtraron los eventos de inundación y avenidas torrenciales más
importantes, eliminando aquellos que tuvieran relación con cuencas vecinas, lluvias locales que afectaron
viviendas o sectores aislados, daños u obstrucción de alcantarillados o canales. Solamente quedaron las
siguientes fichas para el río Otún y para la Quebrada Dosquebradas, las que se podrían considerar como

MUNICIPIO DE PEREIRA – CARDER. Diagnóstico de riesgos ambientales municipio de Pereira. 65 P. 2004.

11

MUNICIPIO DE PEREIRA – DIRECCIÓN MUNICIPAL DE GESTIÓN DEL RIESGO. Estado del arte de estudios de
12

Gestión de Riesgos en el municipio de Pereira. Pag 22 Pereira, 2019.

1-117

1-117
crecientes importantes antes que eventos torrenciales (Tabla 1-46), si se tiene en cuenta que no tenían
contenidos importantes de sólidos y sus efectos fueron marginales en el sentido que no se presentaron
pérdidas humanas.

La revisión de los datos contenidas en Desinventar.org (129 fichas), indica una relación directa lluvias-
inundación, que en buena parte dice que se trata de crecientes súbitas de caños o quebradas urbanas o
semiurbanas que han tenido intervención de cauces, tal vez disminuyendo su área libre, lo que genera
gran parte de los desbordes reseñados. Numerosas fichas reportan obstrucción de caños o desagües, sin
embargo, existen vacíos importantes de información entre 2009 y 2012.
Tabla 1-46 Eventos de crecientes más relevantes

Fecha Evento Observaciones

1944/09/20 Inundación 2 viviendas afectadas, 12 damnificados, creciente del Otún
1974/05/15 Inundación 50 damnificados, desborda el río
1981/05/15 Inundación 1 víctima por creciente del río
1987/12/01 Inundación 7 viviendas destruidas, 35 damnificados, crece el Otún
1997/01/21 Inundación Daños en acueducto por creciente del río
1985/08/09 Avenida torrencial 4 heridos por bloqueo del canal de hidroeléctrica
2011/10/14 Creciente Área rural, 10 viviendas afectadas, 50 damnificados
2011/12/19 Creciente 6 viviendas afectadas en La Suiza, 45 afectados
2013/07/22 Creciente súbita Q. Dosquebradas sin datos
2014/04/27 Creciente súbita Q. Dosquebradas sin datos

Fuente: Elaboración propia, 2019. Modificado de la base de datos de Desinventar.org.

Un ejemplo de lo anterior es la creciente súbita de la quebrada El Manzano, ocurrida en noviembre de

2018 que afectó parte del corregimiento de La Florida, depositando pantano en algunas viviendas y
socavando su propio cauce. Dada la importancia de este evento, más adelante se hace un análisis en
detalle dada su valor como indicador de alerta temprana.

1.11.5 Otros eventos recientes.

Algunas crecientes del río Otún por circunstancias diversas, en especial daños, tienen reportes que
consignan algunos detalles importantes acerca de los eventos y por esta razón se incluyen los más
relevantes de ellos como apoyo para algunos conceptos utilizados en este trabajo.

Creciente del 2 de noviembre de 1926: En el Anuario Estadístico de Pereira, Entrega No1, 1924 – 1925
– 1926 y 1927 (citado por Méndez y Ojeda, 2002)13, se reporta una creciente del Otún con 28 víctimas, 10
desaparecidos y 13 heridos. La mención de desaparecidos y heridos parece indicar más bien un fenómeno
de movimiento en masa que cae parcialmente al río, en razón a que en los eventos torrenciales no son
comunes estos datos, pues en el Otún, dado su carácter de rio de montaña, las personas que hubiesen
sido arrastradas se encontrarían con relativa facilidad aguas abajo, al tiempo que personas heridas son

13MÉNDEZ, R., OJEDA, J. Evaluación de los efectos de la creciente del río Otún el día 28 de marzo de 2002. 8 p.
2002.

1-118

1-118
pocas en eventos torrenciales. El sustento de la interpretación anterior se encuentra en el siguiente párrafo
de la narrativa sobre la construcción del ferrocarril de Caldas: “El mismo paso forzado diseñado por el
Mayor Morris para el Otún, que obligó a un continuado uso de dinamita en las vertientes rocosas, fue causa
de la mayor tragedia sucedida durante la construcción del ferrocarril, cuando el 2 de noviembre de 1926,
se combinaron una creciente del río más el derrumbe de un gran talud, que cayó sobre el agua represada,
ocasionando un “tsunami” que arrastró un pequeño poblado de casuchas de trabajadores, cafés y
prostíbulos, apelado “La huella del tigre”, con más de medio centenar de víctimas”14.

Creciente del 22 de octubre 1979: Este episodio aparece en las bases de datos de Desinventar y se
reportan 7 víctimas y 40 viviendas destruidas, además de daños en acueducto y vías. En el estudio de
Parra y James (1984), se determinó un área mojada máxima en el sector de La Bananera de 133.25 m 2.
En el mismo no se hace referencia a una avenida torrencial sino a un episodio de creciente y se determina
que el desastre tuvo un agravante social consistente en la invasión de viviendas construidas en las
márgenes del río, algunas de ellas sobre muros con alturas insuficientes. Desinventar reporta que hubo un
movimiento en masa en la vereda El Cedral; según testimonio de personas de Pereira, no se conocían
eventos de este tamaño en el siglo XX. De acuerdo con los registros encontrados, este es uno de los 3
eventos de mayor magnitud ocurrido en los últimos 50 años en Pereira, con registro de 200 m3/seg en la
estación La Bananera, en operación desde 1971.

1.11.6 Análisis de torrencialidad de la cuenca.

Para el conocimiento de la amenaza por avenidas torrenciales del lote, se recorrió la cuenca del río Otún
desde el sitio La Pastora hasta su desembocadura al Cauca con el fin de determinar la tipología de los
depósitos de los afluentes que mostraron geoformas de conos o terrazas en la fase previa de
fotointerpretación y de esta manera determinar los que hayan tenido eventos de clase torrencial que
pudiesen haber alcanzado a afectar el lote. De manera semejante, se investigaron las confluencias y
sectores aguas arriba de los afluentes principales que tuviesen posibilidad de transportar flujos de grandes
volúmenes. Los sitios estudiados en detalle en cada cuenca se ilustran en la Figura 1-54. Para efectos de
este trabajo, los afluentes se describen en orden de confluencia con el Otún desde arriba hacia abajo y los
diferentes niveles de las terrazas de cada cauce se numeran desde T 0, el más bajo y reciente, hasta Tn el
más alto y antiguo.

Solamente para efectos de este trabajo, la cuenca del río Otún se divide en tres tramos que se diferencian
en la morfología fluvial del cauce y sus laderas. Otún Alto, desde la laguna del Otún hasta la
desembocadura del río San Juan, Otún Medio desde este sitio hasta el puente de la vía a Marsella y Otún
bajo, desde este puente hasta la desembocadura al Cauca. Esta división facilita el orden para la
descripción de los cauces y los depósitos encontrados en las márgenes o confluencia con el Otún y su
relación con la dinámica fluvial.

14
DE LOS RÍOS, Ricardo. Rieles y montañas: Cuando el eje cafetero empezó a marchar sobre ruedas. Rev. Ciencia
Nueva, V.1, N.1. pag 104-135. 2017

1-119

1-119
Fuente: IGAC.
Figura 1-54. Delimitación de la cuenca y sitios estudiados sobre base IGAC 1:200.000

1.11.6.2 Otún alto.

Desde su nacimiento en la Laguna del Otún hasta un poco más abajo de los 3.000 m, tal como se aprecia
en la Figura 1-55, tanto para las cabeceras del Otún como para el Barbo, se encuentra vegetación de
páramo la que se asocia directamente con baja precipitación, incapaz de generar flujos torrenciales. La
textura lisa de la superficie de este paisaje, carente de cicatrices de movimientos en masa o socavación
de cauces, evidencia una mínima incisión de los depósitos no consolidados, lo que da lugar a interpretarlo
como un espacio dentro del cual no transitan avenidas torrenciales, al menos desde el pleniglaciar, hace
más de 11. 000 años.

1-120

1-120
Fuente: Modificado de Bing, 2019.
Figura 1-55. Morrenas y cicatrices de deslizamientos de la parte alta del Otún – Barbo

Continuando desde el anterior paisaje hasta El Cedral, a 2200 m.s.n.m., se ven morrenas laterales y
depósitos fluvioglaciares de edad tal vez más antigua de 50.000 años que hasta ahora no tienen una
cartografía sistemática en la temática glaciar (Thouret, 1985). La conservación de estas geoformas hasta
esta cota es evidencia de que no se han presentado ni movimientos en masa de importancia o enjambres
de ellos, ni socavación severa que hubiese erosionado las formas glaciares existentes como las morrenas.

En la parte media de este último segmento, en la cuenca del río Barbo, ya sobre la zona de bosque
húmedo, se pueden apreciar algunas cicatrices de movimientos en masa revegetalizadas a estado de
rastrojo, a las cuales se hace referencia en el siguiente numeral.

En el depósito que conforma la zona plana de La Pastora y sus alrededores (Figura 1-56), se puede
visualizar la dinámica holocénica del cauce del Otún Alto, donde en las vertientes se nota el modelado
denudacional por movimientos en masa, de los cuales quedan cicatrices en distintos grados de
conservación de sus coronas y flancos. Estos rasgos denotan que los movimientos en masa no fueron
simultáneos ni por eventos de lluvia que hubiesen producido enjambres ni por sismos, en ambos casos
generadores de avenidas torrenciales. En la base de una de las cicatrices señaladas en la misma figura
se aprecia el cono respectivo que permite dimensionar el tamaño del movimiento en masa que lo generó.
El redondeamiento avanzado de los bordes de los diferentes depósitos que conforman esta zona de La
Pastora, incluyendo las terrazas más bajas, implica una longevidad superior a los varios miles de años
para su depositación y al tiempo evidencia la ausencia de avenidas torrenciales que hayan transitado por
el río en esta zona durante este lapso de tiempo. El puente de El Cedral, con más de 60 años de existencia
y con un área libre de 9x6=54 m2, no muestra ninguna evidencia de haber sido afectado por el tránsito de
una avenida torrencial desde su construcción hasta la fecha; los líquenes verdes que cubren algunas de
las rocas presentes en mitad del cauce parecen indicar también más de 60 años sin que haya pasado una
corriente con suficiente contenido de sólidos como para borrarlos.

1-121

1-121
Fuente: Modificado de Google earth, 2018.
Figura 1-56. Rasgos morfodinámicos y de depositación, Otún, sitio La Pastora

En conclusión, el río Otún desde su nacimiento hasta El Cedral no ha transitado avenidas torrenciales por
lo menos en las últimas decenas y tal vez centenares de años con el tamaño suficiente como para superar
el puente actual en servicio y su amenaza por torrencialidad se puede considerar como baja.

1.11.6.3 Río Bárbaro.

Este afluente nace en el Paramillo de Santa Rosa y en su parte superior contiene casi todo su paisaje
periglaciar y de páramo aferente a la cuenca del Otún, con excepción del sector NW que alimenta al río
San Eugenio. Dentro de esta zona presenta afloramientos rocosos y unos pocos arenales producto de la
gelifracción que se ocurre en la actualidad por la acumulación ocasional de nieve. Sus productos quedan
retenidos por la vegetación de gramíneas que crece más abajo en la zona de páramo, sin mostrar ninguna
área de acumulación de importancia. Cabe señalar que en este sector las precipitaciones son muy bajas,
donde rara vez superan los 40 mm/día, de acuerdo con los registros pluviométricos de la estación cercana,
La Laguna, en operación desde 1991, en la que como caso excepcional se midieron 95 mm en diciembre
23 de 2003, sin que se generaran situaciones anormales a la fecha en las instalaciones de Pez Fresco
que es el sitio con ocupación más alto y vulnerable de este cauce.

En su parte media – alta, cerca de la cota 2800 que es el máximo pluviométrico de acuerdo con la
vegetación existente, se aprecian varias cicatrices de movimientos en masa que se prolongan hasta el
cauce indicando así la formación de flujos de escombros por exceso de lluvia. Su estado de
revegetalización a rastrojo medio a alto permite asignarles una edad tal vez mayor de 30 años, de acuerdo
con el comportamiento de las vertientes en Colombia a estas alturas (Parra, 2016) y porque las cicatrices

1-122

1-122
correspondientes a la lluvia de 2007 apenas se empiezan a colonizar con rastrojo bajo en la imagen BING
2018 (Figura 1-57), aunque parte de ellas se encuentran cubiertas parcialmente por nubes. Las cicatrices
citadas evidencian un volumen relativamente pequeño de masa que cayó al cauce para alimentar un flujo
torrencial de tamaño reducido que no alcanzó el volumen de sólidos necesario para conformar una avenida
torrencial de magnitud importante. Pese a lo anterior, existen en la imagen BING huellas de arrasamiento
de la vegetación en el fondo del cauce con el mismo estado de crecimiento de la vegetación de rastrojo
que se diluyen progresivamente hasta llegar a Pezfresco, unos 3 Km aguas abajo.

Fuente: Modificado de Bing, 2019.

Figura 1-57. Huella del flujo torrencial o creciente de mayo de 2007

En el cauce y los bordes inmediatos del río en el puente de Pezfresco, Fotografía 10-1, se observan
características típicas de un río sin depósitos torrenciales como lo es la terraza T1 que tiene alrededor de
2 m de altura sobre el nivel de aguas medias, compuesta de con bloques de roca decimétricos,
redondeados, con contactos tangenciales entre sí, indicativos de un régimen aluvial. La superficie de los
bloques se encuentra limpia, indicando así que en un lapso no mayor de 20 años ha transitado una
creciente con los suficientes sólidos como para esmerilar su superficie y borrar las huellas de líquenes que
hubiesen podido crecer sobre ella. No obstante, lo anterior, en la terraza respectiva se encuentran árboles
maduros, de varias decenas de años de edad, sin huellas de abrasión por flujos de tipo torrencial que
indiquen que este nivel haya sido alcanzado en por lo menos unos 30 años. En la Figura 10 se muestra el
efecto del evento de creciente del 14 de octubre de 2017, donde se aprecia la composición de la terraza
T0, inundable con una frecuencia de 3 a 5 años, formada por limos y arenas que también demuestran su
carácter aluvial en épocas recientes.

Un criterio adicional acerca de la recurrencia de grandes caudales que pudieron bajar por el cauce es la
estratigrafía de los depósitos, en especial de los más recientes que son las terrazas más bajas. Los
horizontes edáficos A y B se desarrollan a largo plazo, del orden de centímetros por miles de años y en el
evento de que ocurra una avenida torrencial que pase sobre ellos, los cubre con sedimentos para
comenzar un nuevo desarrollo. La terraza más baja ya enunciada tiene un perfil edáfico que supera los 30

1-123

1-123
cm de horizonte A, indicativo de que hace al menos varios miles de años no ha ocurrido un evento torrencial
ni aluvial que supere este nivel y deposite materiales sobre la geoforma.

Fuente: Modificado RHDR, 2011.

Fotografía 10-1.Cauce del río Barbo en Pezfresco luego de la creciente del 14 de octubre de 2011

En conclusión, en el río Barbo se ha formado recientemente al menos una creciente de magnitud que no
alcanza a considerarse como un flujo torrencial (2007), que alcanzó a llegar relativamente diluido a la
bocatoma de Pezfresco causando daños significativos, pero su caudal no superó los 2.0 m de altura en el
sitio del puente de acceso a esas instalaciones.

En conclusión, este cauce se puede clasificar su torrencialidad en este tramo como intermedia en razón al
tamaño reducido de las crecientes históricas.

1.11.6.4 Río San Juan.

El San Juan nace al W del Paramillo de Santa Rosa sin tener influencia directa de su zona periglaciar, sino
que sus cabeceras se encuentran en zona de bosque húmedo. Este cauce no parece tener un historial
muy destacado de crecientes destructivas de acuerdo con la base de datos, pero los días 6 y 14 de octubre
de 2011 presentó crecientes importantes que produjeron afectaciones leves a moderadas en sus
márgenes15. En ambas fechas el río subió un poco más de 1 m del nivel de su cauce normal en el puente
de La Florida y en el área urbana de Pereira, tal como se muestra en la Fotografía 10-2 A y B. A pesar de
que para esa fecha ya se encontraban operando algunas estaciones de la -RHDR- Red Hidroclimatológica
del departamento de Risaralda, las mediciones correspondientes a estos eventos no son muy concluyentes

15RED HIDROMETEOROLÓGICA DEL DEPARTAMENTO DE RISARALDA. Informe ejecutivo del

evento del 14 de abril de 2011. Pag 15. Pereira, 2011.

1-124

1-124
porque la estación de aforo del San Juan fue destruida por el evento del día 6 y para el día 14 y la posición
del pluviómetro estuvo un poco por fuera del núcleo del aguacero principal, pues sólo se registraron 65.1
mm en Pezfresco, valor que no es significativo como para generar un evento destructivo en ninguno de los
cauces de la parte alta del Otún. Unas pocas cicatrices de movimientos en masa en proceso de
revegetalización a rastrojo bajo se pueden observar en la imagen BING de la parte alta de la cuenca, los
cuales pueden corresponder a este episodio de lluvia.
Una anotación adicional tiene que ver con los grandes bloques de roca presentes en el cauce de la parte
baja del río, donde ellos muestran restos de colonias bien desarrolladas de líquenes verdes que sufrieron
abrasión casi total en un episodio que bien puede ser el de 2014. Esta observación permite deducir que la
creciente anterior, de características similares, había ocurrido tal vez más de 70 años atrás.

Fuente: Modificado RHDR, 2011.

Fotografía 10-2.Huellas de creciente San Juan el 14/10/2011. B) el Otún urbano el mismo evento

Por último, los depósitos que conforman las terrazas del San Juan muestran bloque decimétricos y cantos
centimétricos de roca redondeados, clasto –soportados, con contactos tangenciales entre sí y matriz
areno-limosa suelta, cubiertos por al menos 40 cm de horizonte A mezclado con arenas y limos derivados
de crecientes no violentas del río.

En síntesis, al igual que el río Barbo, el San Juan ha transitado solamente crecientes importantes en los
últimos centenares a miles de años y en consecuencia su amenaza por avenidas torrenciales se puede
clasificar como intermedia hasta su confluencia con el Otún donde el contenido de sólidos se diluye y
pierde poder destructivo.

1.11.6.5 Otún medio.

En este tramo, el Otún no presenta depósitos torrenciales en las terrazas más bajas que son las que
reflejan la dinámica fluvial reciente; todos los depósitos, incluidos los pequeños abanicos de los afluentes
principales presentan características típicas de cauces no torrenciales como son el redondeamiento de
ellos, sus contactos tangenciales, matriz arenolimosa suelta y ausencia de bloques mayores subparalelos
al cauce o barras torrenciales. Un afloramiento típico de este tramo, cerca de La Bananera se muestra en
la Fotografía 10-3, donde en la parte inferior del depósito que forma la terraza T2, se aprecian los bloques

1-125

1-125
con alto redondeamiento y contactos tangenciales entre sí, con una cubierta no apreciable en la imagen,
compuesta por una capa de arenas y limos de desborde de inundaciones, mezclada con cenizas
volcánicas que alcanza más de 40 cm de espesor, que implica una edad mayor a varios miles de años
para el depósito de bloques.

Fuente: Elaboración propia, 2019.

Fotografía 10-3.a) Terraza T2 típica del cauce medio (urbano) del Otún. b) Marca de erosión
reciente en la margen del Otún, sector La Bananera

Al interior del casco urbano donde las geoformas naturales han sido intervenidas por completo con la
construcción de muros y explotación de material de playa, en uno los pocos sitios no intervenidos del sector
de La Bananera, se puede observar el conjunto T 0 la de inundación frecuente (3 a 5 años) y T1, la que se
inunda cada 10 a 15 años, donde esta última muestra un escalón ´que evidencia erosión reciente y ninguna
de las dos posee barra torrencial indicadora de eventos de flujo con alto contenido de sólidos.

Vale la pena advertir que los grandes bloques contenidos en el cauce son producto de la incisión
(profundización) sobre el Abanico de Pereira que es la unidad geológica que los contiene y es
predominante en las vertientes inmediatas, más que por arrastre directo en eventos de crecientes. Además
de lo anterior. El Otún a lo largo de su recorrido es un río de montaña, donde gracias a su pendiente, las
crecientes normales presentan una velocidad relativamente rápida, del orden de los 5 m/seg. o más, lo
que convierte su caudal en altamente peligroso para personas y estructuras.

De acuerdo con las observaciones anteriores, se puede concluir que el río Otún en este tramo se puede
clasificar en amenaza baja por torrencialidad. Pese a lo anterior, en este tramo se presenta la mayor
vulnerabilidad por dos causas. La primera consiste en la ocupación del cauce mayor por efectos de
invasión para vivienda no planificada que en partes ha construido muros, disminuyendo el área libre de
crecientes del río y de esta manera se aumenta el riesgo de inundación. La segunda se relaciona con el
alto porcentaje de suelo urbanizado que no permite infiltración de la lluvia, propiciando así lo que se conoce
como inundaciones rápidas (flash floods), aún en las cañadas de caudal efímero que cruzan los diferentes
barrios y que son y serán en el futuro de alto peligro si no tienen los drenajes suficientes para evacuar los
episodios de lluvias fuertes a muy fuertes características de la zona urbana de Pereira y Dosquebradas.

1-126

1-126
1.11.6.6 Río San José.

Esta corriente descarga sus aguas al Otún en un sector relativamente estrecho del fondo de su valle,
propiciado por el afloramiento de un macizo rocoso compuesto por basaltos masivos. No existe un abanico
de depositación generado por el río San José, sino un islote bajo, estrecho y alargado en dirección aguas
abajo que alberga una vivienda campestre de vieja data, tal vez de más de 50 años de haber sido
construida. En la confluencia del San José al Otún (Fotografía 10-4), se puede observar en primer plano
un bloque de roca cubierto por una colonia de líquenes verdes, lo que demuestra la carencia de eventos
torrenciales que en más de 50 a 60 años no los han destruido, conclusión que es válida para ambos
cauces.

Fuente: Elaboración propia, 2019.

Fotografía 10-4.Aspecto de la confluencia del río San José al Otún

A este río se le atribuyó el evento de inundación del 28 de marzo de 2002 (Méndez y Ojeda, 2002), que
causó daños en la vía al corregimiento de La Florida, en especial del sector de Colpapel hacia abajo. En
esta ocasión, la instrumentación pluviométrica de la estación de Pezfresco en el río Barbo registró una
lluvia de 82.5 mm, pero en Nuevo Libaré, cerca de la desembocadura del San José al Otún, sólo se
registraron 6.4 mm, lo que indica que el núcleo de la lluvia se encontraba en la parte media alta del San
José, con un área subordinada en la cuenca del río Barbo en cercanías de Pezfresco. El caudal máximo
medido en la estación La Bananera el mismo día fue de 189.2 m 3/seg, cifra que ha sido superada varias
veces así: 22 de octubre de 1979 (200 m3/seg), 6 de diciembre de 1990 (196.3 m 3/seg), 30 de mayo de
1996 (205.5 m3/seg), 26 de enero de 1999 (201,3 m 3/seg) y 16 de octubre de 2012 (210.0 m 3/seg), donde
sólo el evento resaltado en negrilla ha producido víctimas.

1-127

1-127
1.11.6.7 Quebrada Dosquebradas.

Este afluente converge con el Otún en zona suburbana de Dosquebradas, formando un pequeño abanico
que muestra una pendiente importante en dirección al Otún, mostrando así que la dinámica fluvial en el
sitio ha estado dominada por el río. El depósito formado inmediatamente aguas abajo (Fotografía 10-5),
muestra características netas de comportamiento aluvial como son: bloques con alta redondez con
contactos tangenciales entre sí, matriz subordinada que rellena intersticios e imbricación en dirección de
la corriente, lo que descarta el carácter torrencial de este cauce y además una cubierta de horizonte edáfico
A de más de 30 cm de espesor que denota su longevidad en por lo menos varios miles de años.
La planicie sobre la cual está construido el casco urbano de Dosquebradas, muestra en sus depósitos un
carácter netamente aluvial, descartando la torrencialidad de sus afluentes. Una datación C 14 obtenida de
un tronco de madera contenido en este depósito, reportó una edad superior a 37.000 años (Parra y James,
1983), indicador que corrobora lo observado en su desembocadura al Otún.

Fuente: Elaboración propia, 2019.

Fotografía 10-5.Depósito que conforma la Terraza-abanico de la Q. Dosquebradas

Esta quebrada tiene tres connotaciones que la hacen susceptible de inundaciones rápidas como son el
alto porcentaje de la cuenca urbanizado que inhibe la infiltración de la lluvia en el suelo y aumenta los
caudales pico durante lluvias moderadas a fuertes, también influye su forma relativamente circular,
permitiendo así que los caudales pico de los afluentes lleguen casi al tiempo a la corriente principal y por
último, la alta intervención del cauce con canalizaciones rectas y largas que tienen como consecuencia el
aumento severo de la velocidad de la corriente. Este aspecto deberá ser tenido en cuenta como prioritario
para los planes de desarrollo municipales con el fin de evitar pérdidas cuantiosas en infraestructura e
impactos sobre las personas que habitan zonas ribereñas. Para este efecto, la densificación y seguimiento
estricto de la red hidroclimatológica, en especial lluvias y caudales es clave para el diseño de obras de
drenaje futuras.

1-128

1-128
En conclusión, la quebrada Dosquebradas tiene un historial negativo para avenidas torrenciales y por lo
tanto su amenaza por este concepto es baja. Es preciso hacer una salvedad de la calificación anterior que
está dada por su alta susceptibilidad a inundaciones rápidas facilitadas de manera especial por las
canalizaciones que han rectificado su cauce y alto grado de urbanización.

1.11.6.8 Otún Bajo.

Desde el Puente Marsella y en dirección aguas abajo, el Otún muestra un fondo más estrecho que los
tramos anteriores y a la vez su cañón es más profundo con respecto al abanico hasta su desembocadura
al Cauca, mostrando un valle en V de unos 800 m de ancho por 100 de profundidad. Solamente tiene un
afluente de importancia que es la quebrada La Arenosa, la que, por su baja pendiente por tener su cuenca
al interior del Abanico, no tiene un cauce que sea potencialmente generador de avenidas torrenciales y
únicamente se deben presentar inundaciones súbitas.

El río Otún en su desembocadura al Cauca muestra una serie de depósitos aluviales que forman un cono
deflectado hacia aguas abajo por la corriente del Cauca, lo que indica que no se han presentado al menos
en los últimos miles de años aportes tan significativos como los de una avenida torrencial que provocaran
la formación de un abanico normal y obligar al Cauca a desviarse hacia la margen opuesta. Los depósitos
forman un complejo de terrazas aluviales con más de 5 o 6 niveles (Fotografía 10-6), los más altos de ellos
muy antiguos, tal vez con una edad de cientos de miles de años, al punto que se encuentran incisados a
colinas y se pueden apreciar en el corte de la vía antes de llegar al caserío de Estación Pereira. Todos los
niveles de terrazas, incluyendo los más antiguos muestran depósitos de características netamente
aluviales, semejantes a los otros tramos del Otún.

Fuente: Elaboración propia, 2019.

Fotografía 10-6.Sistema de terrazas aluviales del Otún en Estación Pereira

1-129

1-129
Un indicador importante acerca de las crecientes del Otún en este sitio es el puente del ferrocarril, que fue
construido tal vez al final de la década de 1920 (De los Ríos, 2017) y destruido por una creciente en 2012
como ya se mencionó en numeral anterior.

1.11.6.9 Quebrada la arenosa.

Desemboca aguas abajo del Proyecto PTAR aproximadamente en la cota 950 m y en conjunto con el río
Otún forma un abanico de unos 600 m de longitud por 250 de ancho, con ápice en el Otún, aguas arriba
de su confluencia con la Arenosa. El anterior rasgo indica un aporte menor de sedimentos por parte de la
quebrada y que al igual que los afluentes situados aguas arriba, sus episodios de crecientes son poco
significativas con respecto al caudal del Otún, en razón a que éste último traslada la carga de arrastre
hacia aguas abajo sin mayores dificultades. Una evidencia adicional se observa en la pendiente del
depósito de la desembocadura de la quebrada que muestra una pendiente subordinada a la del depósito
formado por el río.

Se concluye que esta corriente no tiene en su dinámica fluvial rasgos que indiquen la presencia de
avenidas torrenciales en los últimos miles de años y, en consecuencia, su amenaza por torrencialidad es
baja.

1-130

1-130
1.12 Conclusiones Ronda Hídrica
Dentro del proceso de definición de la ronda hídrica, se define que la misma depende de tres componentes
principales:

 Geomorfológico
 Hidrológico (dinámica de los cuerpos de agua y en ellos de los sedimentos y de los nutrientes)
 Ecosistémico.

Dentro de la metodología desarrollada por la CARDER para la definición de la ronda hídricas en los
aspectos hidrológicos e hidráulicos, se aclara que:

1. En cauces con llanura inundable ocupada, el componente hidrológico de la ronda hídrica se

puede definir por la zona que queda debajo del nivel de aguas máximas de una crecida con un
período de retorno de 15 años.
2. En cauces con llanura inundable ocupada, el componente hidrológico de la ronda hídrica se debe
definir por el espacio que requiere el flujo en un evento de mayor importancia y al menos con 100
años de período de retorno.
3. Realización de crónicas con los pobladores de la zona en donde se pueda datar las crecientes
históricas de los cauces analizados

En este sentido es de aclarar que las manchas de inundación presentadas en este estudio cumplen con
estos requisitos desde el punto de vista de hidrología e hidráulica, sin tener en cuenta componentes de
sedimentos o nutrientes.

Para el componente geomorfológico de la ronda hídrica es necesario considerar la faja comprendida entre
el cauce del río y el borde de la terraza T2 en dirección al río, pues ella no es inundable y además no existe
un borde del río que esté afectado por amenaza por movimientos en masa, que en este caso habría que
hacer otro retiro, en conclusión, las manchas de inundación abarca el componente geomorfológico.

1-131

1-131
BIBLIOGRAFÍA
ARBOLEDA, J., ZULUAGA, J. El concepto del riesgo ambiental y su evaluación. Revista Empresas
Públicas de Medellín, volumen 15, No. 3. Enero – Abril de 2005

CONSORCIO ORDENAMIENTO CUENCA RÍO OTÚN. (2017). POMCA Ajuste plan de ordenación y
manejo de la cuenca del río Otún. Fase de diagnóstico. Pereira.

Lalinde, C. (2004). Evidencias paleosísmicas en la región de Pereira - Armenia, Colombia.

SERVICIO GEOLÓGICO COLOMBIANO. (1984). Geología de la plancha 224 Pereira. Medellín.

Servicio Geológico Colombiano. (2011). Propuesta de estandarización de la cartografía geomorfológica en

Colombia.
Aldridge, B. N., & Garrett, J. M. (1973). Roughness coefficients for stream channels in Arizona.
US Geological Survey.73, 3.
Barón, A. (2003). Modelos geoespaciales de la distribución de las variables climatológicas en el
territorio colombiano. Meteorol. Colomb. ISSN 0124-6984, 81-89.
Benson, M., & Dalrymple, T. (1967). General field and office procedures for indirect discharge
measurements. U.S. Geological Survey Techniques of water resources investigations,
Book No. 3.
Bolaños, C. (1995). Regionalización de caudales mínimos en el departamento de Antioquia.
Brunner, G. W. (2016). HEC-RAS River Analysis System. Hydraulic Reference Manual. Version 5.0.
Davis, CA: USACE.
Budyco, M. J. (1974). Climate and life. Academic Press, 508.
C. C. (2017). Actualización POMCA río Otún. Risaralda.
CARDER. (2014). Condiciones hidrológicas y climáticas en el área de jurisdicción de la
Corporación Autónoma Regional de Risaralda. Risaralda.
CHOW, V. M. (1998). Applied hydrology. New York.: McGraw-Hill.
Chow, V. T. (1983). Hidráulica de los Canales Abiertos. Mexico: Editorial Diana.
Chow, V. T. (1994). National Engineering Handbook. Part 630-Hidrology. Hidrology, Cap 4.
COWAN, W. L. (1956). Estimating hydraulic roughness coefficients. Agricultural Engineerin, Vol.
37.
de ALMEIDA, I. K. (2014). Estimation on time of concentration of overland flow in watersheds: a
review. Geosciences= Geociências, 33(4), 661-67.
Dixon, W. M. (1983). Introduction to Statistical Analysis (Fourt edition). New York: McGraw-Hill
Book Company.
Feldman, A. D. (2000). Hydrologic modeling system HEC-HMS: technical reference manual. US
Army Corps of Engineers, Hydrologic Engineering Center.

1-132

1-132
GIANDOTTI, M. (1940). Previsione empirica delle piene in base alle precipitazioni meteoriche,
alle caratteristiche fisiche e morfologiche dei bacini; Applicazione del metodo ad alcuni
bacini dell’Appennino Ligure. Memorie e Studi Idrografici, v. 10, 5-13.
Grubbs, .. (1969). Procedures for Detecting Outlying Observations in. Technometrics, 1-21.
Hershfield, D. (1961). Rainfall Frequency Atlas of the United States for Durations from 30
Minutes to 24 Hours. Washington, D.C.: National Weather Bureau.
Hydrologic Engineering System . (1991). HEC 2, Water Surface Profiles, User Manual. David, CA:
US Army Corps of Engineering.
IDEAM, I. C. (2007). Metodología Corine Land Cover Adaptada para Colombia.
INVIAS. (2009). Manuel de drenaje para carreteras. Bogota.
Jaramillo, A. (1989). Relación entre la evapotranspiración y los elementos climáticos. Cenicafé,
40.
Kang, H. M., & Yusof, F. (2012). Homogeneity tests on daily rainfall series. Contemp. Math.
Scienc, 9-22.
Kendall, M. G. (1955). Rank correlation methods. London: Charles Griffin & Co. Ltd.
KIRPICH, Z. P. (1940). Time of concentration of small agricultural watersheds. Civ. Eng., 362.
Leguizamón, B. A. (2002). Modelos geoespaciales de la distribución de variables climatológicas
en el territorio colombiano. Bogotá: Universidad Nacional de Colombia. Facultad de
Ciencias.
Mann, H. B. (1945). Nonparametric tests against trend. Econometrica: Journal of the
Econometric Society, 245-259.
MATA-LIMA, H., VARGAS, H., CARVALHO, J., GONÇALVES, M., CAETANO, H., MARQUES, A., &
RAMINHOS, C. (2007). MATA-LIMA, H.; VARGAS, H.; CARVALHO, J.; GONÇALVES, M.;
CAETANO, H.;MARQUES, A.; RAMINHOS, C. Comportamento hidrológico de bacias
hidrográficas:integração de métodos e aplicação a um estudo de caso. Escola de Minas,
60.
MOTH, M. o. (1998). Hydraulics manual. Columbia: Engineering Branch. Province of British.
Pettitt, A. (1979). A non-parametric approach to the change-point problem. Applied statistics,
28(2), 126-135.
Pierson, T. (2005). Hyperconcentrated flo: transitional process between water flow and debris
flow. In Debris-flow hazards and related phenomena (pp. 159-202). Springer, Berlin,
Heidelberg.
Pulgarin, E. (2009). Fórmulas regionales para la estimación de curvas intensidad-frecuencia-
duración basadas en las propiedades de escala de la lluvia (Región Andina Colombiana).
Medellin.

1-133

1-133
Rodríguez Espinosa, D. M., Córdoba Guerrero, G., & Delgado-Granados, H. (2017). Evaluación
probabilística del peligro por lahares en el flanco NE del Volcán Popocatépetl. Boletín de
la Sociedad Geológica Mexicana. Boletín de la Sociedad Geológica Mexicana.
Rojo, J. D. (2000). Algunas metodologias para cuantificar la evapotranspiración media
multianual en Colombia. Bogota.
Samuels, P. (September de 1990). Cross section location in one-dimensional models. In
International Conference on river flood hydraulics, 339-350.
Thornthwaite, C. W. (1948). An approach toward a rational classification of climate. Vol. 66, No.
1, p. 77.
TURC, L. (1954). Le bilan d’eau des sols: relation entre les precipitations. En: l’évaporation et
l’écoulement. Ann. Agron, 491-596.
Yen, B. C. (1992). Dimensionally homogeneous Manning's formula. Journal of hydraulic
engineering, 1326-1332.

1-134

1-134