Avenidas Torrenciales Metodo PDF

Revista Ingenierías Universidad de Medellín
EVALUACIÓN DE ZONAS DE AMENAZA POR AVENIDAS

TORRENCIALES UTILIZANDO METODOLOGÍAS CUALITATIVAS.
CASO DE APLICACIÓN A LA QUEBRADA DOÑA MARÍA
Luis Javier Montoya Jaramillo*

Sandra Milena Silva Arroyave**
Juan Esteban González***
Recibido: 26/09/2008
Aceptado: 07/05/2009
RESUMEN
En este artículo se presenta la aplicación de una metodología para la evalua-
ción de zonas de amenaza por avenidas torrenciales, implementada en la cuenca de
la quebrada Doña María. Para esta evaluación de amenazas se utilizaron sistemas
de información geográfica y la técnica de los residuos ponderados, para evaluar
el peso relativo que tenían diferentes factores como la topografía, la pendiente, la
geología y geomorfología, la red de drenaje acumulada, la distancia a los cauces y
el índice de Melton para cada una de las microcuencas. Los resultados obtenidos
permitieron identificar que las microcuencas que mayor amenaza presentaron fueron
la Arenala, la Tablaza, la Limona y la Manguala. Se verificó la amenaza obtenida
mediante reconocimientos de campo y mediante comparación con los resultados
de modelaciones hidráulicas para un caudal de período de retorno de 100 años
incrementado en un 40%.
Palabras clave: avenidas torrenciales, zonas de amenaza, quebrada Doña María.
*
Ingeniero civil, magíster en aprovechamiento de recursos hidráulicos, profesor asistente, Grupo de Investigación y mediciones ambientales
GEMA, Universidad de Medellín, Cra 87 No 30-65, ljmontoya@udem.edu.co, 3405216
**
Joven investigadora, Grupo GEMA, Universidad de Medellín, Cra 87 No 30-65, smsilva@udem.edu.co, 3405216
***
Profesor catedrático, Grupo GEMA, Universidad de Medellín, Cra 87 No 30-65, jegonzalez@udem.edu.co, 3405216
Revista Ingenierías Universidad de Medellín, volumen 8, No. 15, pp. 11-29 - ISSN 1692-3324 - julio-diciembre de 2009/146 p. Medellín, Colombia
12 Luis Javier Montoya Jaramillo, Sandra Milena Silva Arroyave, Juan Esteban González
EVALUATION OF HAZARD ZONES BY TORRENT FLOW

USING QUALITATIVE METHODOLOGIES. APPLICATION
TO THE CASE OF DOÑA MARIA BASIN
ABSTRACT
This paper shows the application of a methodology for evaluation of hazard
zones by torrent flow, it was done at Doña Maria basin. Hazard evaluation was done
using geographical information system, and the residual pondered weight method,
for the evaluation of different factors as topography, slope, geology and geomor-
phology, cumulative flow, distance to the rivers and Melton index for each ones to
the affluent basin. The results obtained allow to identify that the little basin that
more hazard presented were la Arenala, la Tablaza, la Limona y la Manguala. These
results also were verified through field recognitions y comparison with the results
of hydraulic modeling for a return period flow of 100 years increased in 40%.
Key words: torrent flow, hazard zones, Doña Maria basin
Universidad de Medellín
Evaluación de zonas de amenaza por avenidas torrenciales utilizando metodologías cualitativas. Caso de aplicación ... 13
1. INTRODUCCIÓN (QWpUPLQRVJHQHUDOHVVHFXHQWDFRQSRFRPR-
nitoreo y registro sobre eventos de este tipo, lo
Las avenidas torrenciales son fenómenos
que dificulta datar adecuadamente los mismos
naturales que aunque tienen una ocurrencia rela-
y contar con una escala temporal que permita
tivamente baja, por su gran potencial destructivo,
asociar un período de retorno a los mismos o
son de particular interés para evaluar los riesgos
una probabilidad de ocurrencia.
naturales a los que se encuentra sometida una
determinada cuenca. Entre las razones que existen (OFRQFHSWRGHDYHQLGDVWRUUHQFLDOHVWLHQHJUDQ
para que este tipo de evento sea tan dañino están variedad de significados, esto debido a que
su naturaleza casi impredecible, la rapidez a la cual diferentes disciplinas han trabajado el tema
ocurre, su corta duración y su largo período de y, por consiguiente, cada uno hace su aporte
retorno así como su distribución poco uniforme distinto al concepto.
en el espacio y el tiempo [1]. ([LVWHQGLIHUHQWHVFDXVDVTXHGDQRULJHQDXQD
En Colombia por las características topográfi- avenida torrencial; entre ellas se pueden men-
cas de la zona andina, este fenómeno se presenta cionar una alta precipitación que genere una
con una relativa mayor frecuencia. Según INGEO- serie de deslizamientos, un deslizamiento de
MINAS “las avenidas torrenciales son una amenaza gran magnitud que origine represamiento en
natural muy común y posiblemente la menos el flujo, la ocurrencia de sismos que produzcan
estudiada en Colombia; sin embargo causan gran- deslizamientos, o la ocurrencia de erupciones
des pérdidas en vidas humanas e infraestructura. volcánicas.
Solamente entre 1985 y 1995, mas de 40 avenidas (VGLItFLODVRFLDUXQDSUREDELOLGDGGHRFXUUHQ-
torrenciales altamente destructivas ocurrieron en cia a los fenómenos detonantes de las avenidas.
el país, con más de 200 pérdidas de vidas humanas 1RUPDOPHQWHQRVHFXHQWDFRQHVWLPDFLRQHVGH
y cuantiosos daños a viviendas, sistemas de gene- umbrales entre intensidad de la precipitación
ración de energía, carreteras y tierras cultivables” y generación de deslizamientos. Los umbrales
[1]. Polanco y Bedoya [2], utilizando la base de para la generación de avenidas torrenciales de
datos de DesInventar, refieren que en Antioquia orígenes diferentes a altas precipitaciones son
(exceptuando los municipios del valle de Aburrá), de difícil definición.
un 15% de los desastres naturales reportados entre /RV IOXMRV GH ORGRV JHQHUDGRV VRQ GH GLItFLO
los años de 1920 y 1999 corresponde a avenidas modelación. Los parámetros como la densidad
torrenciales. de la mezcla, el volumen de sedimentos, el
Para planificar adecuadamente la ocupación caudal de la avenida son de difícil estimación.
del territorio y en particular una cuenca hidrográfi- (VWLPDUHOYROXPHQGHVHGLPHQWRVDVRFLDGRVD
ca, es necesario tener en cuenta la amenaza asociada un evento con cierta probabilidad de ocurren-
a avenidas torrenciales y especialmente en cuencas cia puede ser difícil.
urbanizadas en las que sus posibles efectos adversos
/DPRGHODFLyQPDWHPiWLFDGHOIHQyPHQRQR
sobre la población, la economía y la infraestructura
se encuentra del todo establecida.
urbana son muy altos.
Uno de los principales inconvenientes para de- La mayoría de los trabajos realizados en ave-
limitar la amenaza por este tipo de eventos consiste nidas torrenciales se concentran en estudiar las
en la dificultad de predecir su ocurrencia. Entre características de zonas en las que se ha presentado
estos se pueden mencionar: un evento torrencial. González et al. [1] estudian
la avenida torrencial y los movimientos en masa Se concluyó que las zonas urbana y sub-urba-
relacionados que ocurrieron el 31 de enero de na, al pie de los cerros noroccidentales de Ibagué,
1994 en la cuenca del río Fraile, fenómeno que pueden ser severamente afectadas por flujos de
generó severas consecuencias para los habitantes tierra y avalanchas originados en la zona monta-
del poblado de Florida (Valle del Cauca). En ñosa, así como por inundaciones. A largo plazo,
dicho estudio se analizaron diversas causas del la posibilidad de daños es mucho mayor ya que
fenómeno, y se encontró que muchos aspectos aumenta la probabilidad de que sucedan eventos
geomorfológicos (como el tamaño de la cuenca; como sismos y lluvias de gran magnitud capaces
el intervalo entre la precipitación y la descarga de detonar deslizamientos más numerosos y de
pico; la forma de la cuenca; la localización en altas mayor tamaño que los ocurridos en años recientes.
montañas con relieves abruptos y clima húmedo; Piedrahíta y Hermelín [4] analizan la avenida
la variabilidad en la descarga; y la red de drenajes) torrencial del río Tapartó (municipio de Andes,
influyeron en la generación del evento. Antioquia), ocurrida en la noche del 25 al 26 de
Se concluyó que la avenida torrencial en el río abril de 1993, la cual dejó depósitos de los flujos
Fraile en 1994 fue un evento de gran magnitud de lodos transportados por el río con un espesor
con un período de retorno posible medido en promedio de 30 cm y causó socavación del lecho
cientos de años. Algunas características geomor- del río. Dichos autores asociaron la ocurrencia de
fológicas intrínsecas hacen de esta cuenca un área dicho evento a fuertes lluvias de gran intensidad,
propensa a la generación de avenidas torrenciales; altas pendientes y desprendimientos de material
sin embargo, la ocurrencia de una inundación de sólido en la parte alta de la cuenca. Mediante el
grandes proporciones, que deje grandes depósitos análisis de 7 depósitos previos en la cuenca, en-
torrenciales, requiere una combinación particular cuentran una recurrencia de estos eventos de 280
de precipitación y de saprolito removido. años para este tipo de avenidas en el río Tapartó.
Barreto y Ruiz [3] realizan un estudio Hermelín et al. [5] presentan una descripción
sobre la amenaza por avenidas torrenciales de la avenida torrencial de San Carlos (Antioquia)
en el sector noroccidental de la ciudad de en el año de 1990.
Ibagué (departamento de Tolima). Se elaboró Flores et al. [6] presentan un estudio sobre
un inventario de desastres y se analizó la la cuantificación e identificación de los tipos de
cronología de los movimientos de masa y las movimientos de masa ocurridos en la cuenca expe-
avenidas torrenciales. Con la base cartográfica rimental de San José de Galipán y su contribución
(ortofotografía y modelo digital del terreno) y al desastre de Vargas (Venezuela) de 1999.
trabajo de campo, se levantó información de En una primera etapa se caracterizó el sitio de
capas de geología y cobertura de la zona. De igual estudio, se estudiaron las características físicas y hu-
modo se tomaron muestras de material superficial manas de predisposición al riesgo de movimientos
y del lecho de las quebradas, se corroboró el uso de masas y flujos torrenciales, se realizó un trabajo
y cobertura del suelo, se realizó un levantamiento de fotointerpretación, el cual, conjuntamente
topográfico de secciones en los cauces y, al con visitas de campo y la información geológica
inventariar los procesos denudativos, se obtuvo recolectada, permitió identificar y cartografiar los
un mapa de susceptibilidad a movimientos de diferentes procesos encontrados. Se encontró que
masa y avenidas torrenciales, el cual fue la base la causa principal de estos flujos torrenciales estuvo
para establecer criterios de calificación de cada asociada a los siguientes factores: presencia de altas
uno de los factores dentro de la evaluación. pendientes, ocurrencia de lluvias de baja intensi-
dad por un período de dos semanas produciendo la ción máxima, sismo, aumento de caudales de
saturación de los suelos, un evento extraordinario la corriente, entre otros. El porcentaje de denu-
donde se registró una precipitación acumulada de dación depende de la magnitud de la anomalía,
911 mm en tres días y la presencia de abundantes el uso del suelo, la pendiente del terreno y la
sedimentos, rocas meteorizadas, espesores variables susceptibilidad a los deslizamientos y/o a la
de suelos, y caídas de bloques. erosión, la cual depende de las características
Diferentes enfoques se pueden adoptar para geológicas. La profundidad de denudación
enfrentar el problema de estimar zonas de amenaza depende de la profundidad del perfil de me-
ante avenidas torrenciales. Una primera alternativa teorización y de la morfología del terreno.
cuando se cuenta con eventos recientes, posibles &iOFXOR GH ORV FDXGDOHV GH DJXD VH HVWLPDQ
de reconstruir, es asociar dicho evento a un de- los caudales de agua disponible la forma y
terminado período de ocurrencia y delimitar las características del probable hidrograma para
zonas afectadas por el evento. El problema de esto la precipitación máxima esperada, suponiendo
es que no en todas las zonas se cuenta con eventos solamente aportes de agua sin tener en cuenta
recientes y en este caso la delimitación sería con el sólidos. Este análisis requiere de información
fin de impedir que se ubique población en zonas sobre precipitación, respuesta de las corrientes,
afectadas, pero no permitiría prevenir a los habitan- uso del suelo y morfología de la cuenca.
tes ante la ocurrencia de un determinado evento.
&iOFXORGHOYROXPHQGHVHGLPHQWRVVHFDOFXOD
Siguiendo la misma idea, una posibilidad es
el volumen de sedimentos que se incorporaría
buscar depósitos aluvio-torrenciales en la cuenca
a la avalancha al momento de la precipitación
y delimitarlos como zonas de riesgo; sin embargo,
máxima, especificando el volumen de sedimen-
el hecho de que en el pasado se hayan presentado
tos finos y gruesos.
eventos en un lugar determinado no es garantía
de que en el futuro se presenten también (aunque &iOFXORGHORVFDXGDOHVGHDYHQLGDWRUUHQFLDOVH
puede ser un buen indicador), pues en la cuenca debe realizar para cada sección o punto de con-
se podrían haber presentado cambios que faciliten trol incluyendo la suma de líquidos y sólidos.
que nuevas zonas se encuentren amenazadas, o que 'HILQLFLyQGHOWLSRGHIOXLGRHOWLSRGHIOXLGR
zonas que antes no lo estaban hayan disminuido depende de la concentración y tipo de sedi-
su amenaza. Además, algunos depósitos podrían mentos. El tipo de fluido determina cuál de
haber sido lavados (algunos resultantes de flujos los modelos matemáticos es el recomendado
de pequeña magnitud), por lo que se podrían para la modelación.
disfrazar algunas zonas. Esta metodología podría
&DOLEUDFLyQGHOPRGHORHVQHFHVDULRFRPSDUDU
ser completada considerando la geomorfología
los resultados del modelo con información his-
de la cuenca y dando alguna conectividad a los
tórica de la zona de estudio y con los resultados
diferentes depósitos. Otra alternativa consiste en
obtenidos mediante otras metodologías.
la modelación. Suárez [7] considera que dentro de
las variables de importancia para la modelación de Aunque actualmente se encuentra disponible
las avenidas torrenciales se encuentran: gran cantidad de modelos, con diferentes grados
$QiOLVLV GH ODV FDUDFWHUtVWLFDV JHROyJLFDV \ de complejidad, predecir los datos necesarios para
geotécnicas: se requiere definir el porcentaje y alimentarlos es difícil. Sin embargo, la modela-
profundidad de denudación de cada subcuenca ción aporta valiosos elementos al monitoreo del
en el momento de ocurrencia de la precipita- comportamiento de las avenidas torrenciales, del
mismo modo que permite una mejor comprensión rocas de todos los tamaños y lodo húmedo que
de los procesos involucrados y apoya la definición aparece súbitamente en el fondo de un cañón pre-
de zonas de riesgos. cedido por un ruido sordo. Mientras la pared pasa
La modelación matemática hace una contri- el canal permanece lleno, con un torrente de es-
bución importante al mejoramiento en la gene- combros compuestos de lodo y bloques, rodando,
ración de advertencias por fenómenos peligrosos chocándose y moliéndose entre sí. Cuando llega a
[8]. Un tercer enfoque, muy cualitativo, es el que un abanico aluvial arrasa estructuras, vehículos y
se presenta en este artículo, en el cual a partir de cubre carreteras, cultivos y campos como un tapete
sistemas de información geográfica se realiza una de escombros que lentamente se va deteniendo
ponderación de diferentes factores con el fin de dejando un montículo frontal”.
identificar zonas sujetas a mayor probabilidad de Suárez [7] considera que una avenida torren-
amenaza, entendiendo estas como aquellas zonas cial es “el fenómeno de erosión en masa que inclu-
que recogen mayor número de factores de riesgo. ye el desprendimiento, transporte y depositación
Con el fin de esclarecer un poco el proceso de las de grandes masas de suelo, por acción del agua.
avenidas torrenciales se parte de algunas defini- Las masas en movimiento se comportan como un
ciones previas, que permitan delimitar un poco fluido pero su comportamiento es diferente al de
el fenómeno que se piensa estudiar. Luego se los flujos de sedimentos, especialmente debido a la
describe la metodología empleada, se presentan alta concentración de sólidos”. Una terminología
los resultados de la aplicación y la verificación de común para clasificar las avenidas torrenciales se
estos mediante reconocimiento en campo y de basa en criterios físicos que tienen en cuenta la
modelaciones hidráulicas. viscosidad y la turbulencia del flujo. De acuerdo
con este criterio los flujos se dividen en tres cate-
2. DEFINICIONES PREVIAS gorías: flujos de lodo, flujos hiperconcentrados y
flujos de detritos.
El concepto de avenidas torrenciales tiene Según la base de datos de Desinventar de la
gran variedad de significados debido a que di- Corporación Osso [11] “una avenida torrencial es
ferentes disciplinas han trabajado el tema y por un flujo violento de agua en una cuenca, a veces
consiguiente cada uno hace su aporte distinto al reportado como creciente (súbito, rápido), o como
concepto. En términos sencillos se puede definir torrente”. Se aplica cuando “la avenida transporta
una avenida torrencial como “el aumento del troncos de árboles y/o abundantes sedimentos
caudal en un cauce con volúmenes excepcionales, desde finos hasta bloques de roca”. “Pueden ser
en el cual, el fluido además de agua contiene una generados por lluvias, por ruptura de represamien-
mezcla de escombros compuesta por suelo, roca y tos o por abundantes deslizamientos sobre una
material vegetal. Este fenómeno está restringido a cuenca. Excluye los aludes porque estos implican
cauces relativamente pequeños de ríos de montaña desprendimiento de hielo o nieve”.
y no se produce en ríos con cuencas grandes, de- Piedrahíta y Hermelín [4] presentan la
bido a que los agentes naturales que las provocan siguiente definición: “los eventos hidrológicos
afectan sólo áreas pequeñas” [9]. torrenciales consisten en la ocurrencia súbita de
De acuerdo con Johnson citado en Castillo caudales pico de corta duración y gran velocidad,
[10] se define un flujo de escombros (avenida que transitan a lo largo de canales en cuencas hi-
torrencial con poca agua y alta cantidad de rocas, drográficas de montaña, normalmente de tamaño
suelo y escombros) como “una pared de bloques, reducido”.
Origen y causas de las avenidas torrenciales Aburrá (también llamado río Medellín). Tiene
su nacimiento en la vereda Astilleros (sur occi-
Las avenidas torrenciales se generan por causas
dente del municipio de Medellín, corregimiento
hidrometeorológicas (lluvias concentradas), sísmi-
de San Antonio de Prado) y su desembocadura
cas (enjambres de deslizamientos cosísmicos), de
al río Medellín está ubicada en el municipio de
inestabilidad de vertientes (bloqueo de un cauce
Itagüí. En la parte alta y media de la quebrada se
por un deslizamiento y posterior ruptura del di-
encuentra el corregimiento de San Antonio de
que), o por erupciones volcánicas y deshielo o por
Prado (municipio de Medellín) y en la parte baja
acumulación de capas gruesas de cenizas sueltas.
se encuentra el municipio de Itagüí. La figura 1
Las avenidas torrenciales también se pueden pre-
presenta la localización geográfica de la cuenca de
sentar por otras causas como son deshielo al final
la quebrada Doña María.
del invierno o lluvias concentradas asociadas a
Los principales parámetros geomorfológicos
ciclones tropicales [9].
para la cuenca de la quebrada Doña María se
3. DESCRIPCIÓN DE LA ZONA presentan en la tabla 1, en la cual se destacan sus
DE ESTUDIO altas pendientes.
La cuenca de la quebrada Doña María está Tabla 1. Parámetros morfométricos de la

localizada en el departamento de Antioquia en quebrada Doña María
la zona centro, al sur occidente de la cuenca del Parámetro Valor
río Aburrá, en la ladera occidental del valle, entre Área de la cuenca [km ] 2
75,82
los municipios de Medellín, Itagüí y La Estrella. Orden de Horton – Strahler 5
La cuenca tiene un área aproximada de 75,82 Magnitud 1371
km2, y es uno de los principales tributarios al río
Longitud total de canales [km] 229,14
Longitud del canal principal [km] 25,84
Cota máxima de la cuenca [m] 3105
Cota en la salida de la cuenca [m] 1522
Distancia al centroide [km] 13,43
Perímetro [km] 49,075
Pendiente media [%] 57,69
Pendiente promedio del canal
5,57
principal [%]
Fuente: AMVA y UNAL [13]
Se tienen registros de eventos con caracterís-

ticas de avenidas torrenciales recientes en algunas
microcuencas afluentes como la Limona, la Man-
guala, la Justa, la Ospina y la Cabuyala, las cuales
presentan fajas de material torrencial, circunscrito
)LJXUD/RFDOL]DFLyQJHRJUiÀFDGHODPLFURFXHQFDGH
Figura 1. Localización geográfica de la a los cauces activos de los afluentes mayores [14],
microcuenca de la quebrada Doña María. con eventos recientes de tipo torrencial reportados
Fuente: AMVA y CTA [12] para el caso de la Cabuyala (16 de junio de 1990).
El estudio realizado por AMVA y UNAL [13], rentes condiciones, además, permite calcular pesos
hace mención a la presencia de depósitos de tipo relativos a cada una de las variables a considerar en
aluvio-torrencial en la quebrada la Candela, la el desencadenamiento de las avenidas torrenciales.
Limona, y la Manguala, y de depósitos intercalados La propuesta parte de la combinación de
(intercalación de depósitos de flujos de lodos, de diversos parámetros previamente normalizados
flujos de escombros, aluviotorrenciales, y aluvia- y se basa en la forma como dichos parámetros
les, en zonas deprimidas), en quebradas como la afectan, de forma individual, la intensidad de un
Macana, y la Manguala. evento torrencial. Dichos mapas normalizados, se
combinaron mediante un sistema de pesos ponde-
rados, donde los valores de dichos pesos relativos
4. DATOS Y MÉTODOS
se obtuvieron a partir de un análisis jerárquico,
Para el estudio de las avenidas torrenciales y entre las diferentes parejas de parámetros consi-
la definición de los niveles de amenaza ante dicho derados. Se realizó un análisis de sensibilidad de
evento, se empleó una metodología basada en el los parámetros considerados dentro del proceso
análisis multicriterio de diversos parámetros que de evaluación multicriterio, del cual se obtuvo
dan razón de la posible intensidad de una avenida un grupo último de pesos relativos a considerar.
torrencial. Hay que aclarar que la amenaza por Los límites de los niveles de amenaza empleados
avenidas torrenciales integra gran cantidad de se definieron teniendo en cuenta un análisis de
fenómenos hidrológicos, físicos, hidráulicos, entre zonas críticas evidenciadas a partir del ejercicio
otros, y en general una descripción completa de de sensibilidad antes mencionado. Los resultados
todos estos fenómenos es compleja, por lo que se finales se cotejaron con la información histórica
presenta una metodología que intenta combinar disponible, y con información de campo, para
diferentes factores para, de forma general, iden- verificar la pertinencia de los resultados. A con-
tificar posibles zonas de amenaza por avenidas tinuación se describe el manejo de cada variable:
torrenciales en una cuenca, que serían donde Para el caso de la cuenca de la quebrada Doña
varios factores se combinan. No se pretende dar María, se utilizó el modelo digital de elevación
una descripción completa del fenómeno o realizar (MDE), el mapa de geología y geomorfología, la
un análisis absoluto del mismo, así que las zonas red de drenaje acumulada, el mapa de pendientes,
obtenidas mediante esta metodología deben ser el mapa de distancia a los cauces y el mapa con el
interpretadas simplemente como referencias de índice de Melton para las microcuencas presentes
análisis para un manejo de las mismas. en la zona de estudio [13].
La metodología utilizada se basa en el método
presentado por Ordóñez y Martínez [15], quienes 4.1 Factores considerados
proponen analizar riesgos naturales (identificación en la ponderación
de zonas de amenaza por inundaciones y por
El modelo digital de elevación (MDE) es un
deslizamientos) mediante los sistemas de infor-
mapa que se obtiene a partir de la cartografía
mación geográfica, utilizando técnicas de análisis
de la zona de estudio, y consiste básicamente en
multivariado. Para la aplicación de este modelo se
construir un mapa digital en el cual, por medio de
empleó el software IDRISI Kilimanjaro, el cual tiene
una extrapolación o interpolación de las curvas de
la capacidad de superponer varias capas temáticas
nivel presentes en el mapa cartográfico, se permite
con el fin de obtener un resultado a partir de dife-
a cada píxel adoptar un valor sobre la altura del
terreno; de esta forma se obtiene un mapa con las Se le dio entonces mayor importancia a aque-
diferentes elevaciones presentes para este caso en la llas microcuencas que presentan valores altos según
cuenca de la quebrada Doña María. Para el estudio el índice de Melton; de esta forma, dichas micro-
se empleó un tamaño de píxel de 10 x 10 metros. cuencas tendrán mayor relevancia al momento de
La red de drenaje acumulada se obtiene a partir ejecutar el modelo.
del MDE mediante la opción Flow la cual permite,
a medida que el cauce se dirige a su desembocadura, 4.2 Funciones de normalización
acumular el número de celdas que descarguen en y transformación
dicho punto. Se quiso con este temático asignar a
cada píxel en la microcuenca un peso de acuerdo Los factores considerados para la ponderación
con la cantidad de puntos que descargasen en ella tienen diferentes unidades y representan otros
y luego, mediante el proceso de normalización para efectos, por lo que no se puede hacer una compa-
la red de drenaje, se asignó una calificación mayor ración directa entre ellos. Para poder combinarlos,
a los puntos en que drenarán un mayor número de es necesario llevarlos todos a una escala común. En
celdas. El índice de Melton se calculó para las mi- este caso se llevarán todas las variables a una escala
crocuencas que conforman la quebrada Doña María entre 0 y 1, donde 0 representa un valor en el cual
clasificando la torrencialidad de cada una de ellas; la variable no favorece la generación de amenazas
este cálculo hace parte del mapa de microcuencas. por avenidas torrenciales y 1 un valor en el cual el
Hay que aclarar que este índice no fue diseñado favorecimiento de amenazas por avenidas torren-
para cuencas en zonas ecuatoriales; sin embargo, ciales sería muy alto.
se utilizó como índice de torrencialidad. Para cada Además de asignar valores entre 0 y 1, se quiere
microcuenca se empleó la siguiente ecuación: representar el comportamiento de la variable, de
tal forma que sea consistente y si al aumentar una
índice de Melton = (H max – H min) × S–0.5 (1)
determinada variable se aumenta el favorecimiento
Donde Hmax y Hmin (en km) son la cota de de la amenaza, entonces el índice se debe aumentar
altura máxima y mínima de la microcuenca y S, y si, por contrario, al aumentarlo, se disminuye el
es el área de la microcuenca expresada en km2. índice, debe disminuir. Para facilitar este proceso de
Una microcuenca presenta torrencialidad cuando normalización, se utilizaron funciones sigmoidales
el índice de Melton presenta valores superiores a crecientes o decrecientes, como las presentadas
0.5 según [16]. Para esta capa temática se realizó en la figura 2. Dichas funciones permiten asignar
un análisis para comparar la torrencialidad entre umbrales inferiores y superiores y presentan un
las microcuencas. Dicho análisis se hizo aplicando comportamiento más o menos lineal en zonas
la siguiente ecuación: intermedias, permitiendo ajustar lo más aproxima-
(IMi – M min) damente posible el comportamiento de las variables
Comparación Torrencialidad = (2)
(M max – M min) a la ocurrencia del evento. Por ejemplo, en el caso
Donde: IMi = índice de Melton de la de la pendiente (para este estudio en específico)
microcuenca cuando ésta tiene valores superiores al 70%, la
IM min = Valor Mínimo del índice probabilidad de ocurrencia es casi nula (se califica
de Melton de todas las con 0); cuando la pendiente varía entre 5% y 70% la
microcuencas amenaza aumenta casi de forma lineal; pero cuando
IM max = Valor Máximo del índice se presentan valores inferiores a 5% la amenaza es
de Melton de todas las alta y, por consiguiente, se adopta un valor de 1.
microcuencas
1
Para normalizar la amenaza de las capas te-
0.8
máticas de geología, geomorfología e índice de
Melton de la microcuencas, no se aplicó ninguna
0.6
función matemática, pues se tenían pocos valores
Índice
y se procedió a calificar cada una de sus caracte-

0.4 rísticas de acuerdo con el rango establecido para
todas las capas temáticas (0-1) según el criterio de
0.2
expertos en el tema.
0
0 10 20 30 40 50 60
4.3 Evaluación multicriterio
(a) sigmoidal creciente
Para evaluar las zonas de amenaza se realiza
1 una ponderación entre los diferentes factores, de
0.9 acuerdo con la siguiente expresión:
0.8
0.7 M = ¦ pi xi (3)
0.6
Donde: M = índice de amenaza
Índice
0.5
0.4
p = pesos asignados a cada factor
0.3
0.2
x = factores
0.1
El análisis se realiza mediante la aplicación
0
0 50 100 150 200 250 300 350 400 450 500 MCE (evaluación multicriterio) la cual permite
Distancia a los ríos (m)
introducir las capas temáticas y sus pesos relativos y
(b) sigmoidal decreciente el programa de forma automática realiza el álgebra
de mapas. El resultado que arroja la evaluación
Figura 2. Función sigmoidal decreciente multicriterio es un mapa con valoraciones entre
Fuente: elaboración propia adaptada de [15] 0 – 1, donde 0 significa la mínima amenaza por
avenidas torrenciales mientras que 1 es la máxima
La tabla 2 presenta las capas temáticas que amenaza del mencionado evento.
se normalizaron, el rango en el que se realizó la
normalización y la función aplicada en el proceso. 4.4 Asignación de pesos relativos
Tabla 2. Funciones matemáticas empleadas en la Para la ocurrencia de avenidas torrenciales

normalización intervienen diferentes variables, sin embargo, el
aporte de cada una de ellas en el desencadenamien-
Rango to del evento no es igual, es posible que algunas
Capa Temática Función matemática
de Valor
variables sean más relevantes que otras. Por tal
Pendiente 0–1 Sigmoidal decreciente
motivo, es necesario determinar cuáles son las
Red de drenaje variables con mayor y menor importancia en la
0–1 Sigmoidal creciente
acumulada
ocurrencia del fenómeno en estudio. Como herra-
Distancia a los mienta de ayuda para la toma de esta decisión, el
0–1 Sigmoidal decreciente
cauces
software IDRISI Kilimanjaro posee una función de
Fuente: elaboración propia ponderación denominada comparación por pares,
Tabla 3. Escala de Valoración de pesos relativos
Valoración Significado
1/4 El factor fila es moderadamente menos importante que el factor de la respectiva columna.
1/6 El factor fila es moderadamente menos importante que el factor de la respectiva columna.
1/9 El factor fila es extremadamente menos importante que el factor de la respectiva columna
1 Igual importancia del factor de la fila y de la columna
4 El factor fila es moderadamente más importante que el factor de la respectiva columna
6 El factor fila es fuertemente más importante que el factor de la respectiva columna
9 El factor fila es extremadamente más importante que el factor de la respectiva columna
Fuente: Valores establecidos en el software IDRISI Kilimanjaro, [15]
la cual permite, de forma cuantitativa, identificar Se realizó un análisis de sensibilidad con di-
cuál o cuáles son las variables con mayor peso en ferentes ponderaciones, con el fin de comparar la
la ocurrencia. Este modelo consiste en una matriz sensibilidad frente a diferentes ponderaciones; para
de doble entrada donde se listan, tanto en filas ello se generaron seis escenarios diferentes. En la
como en columnas, las capas temáticas. El análisis tabla 4, se presentan los valores de ponderación
de valoración se hace confrontando la relevancia empleados para obtener dichos escenarios.
del factor fila respecto al factor columna. La tabla
3 presenta dicha escala de valoración.
Tabla 4. Valores de ponderación empleados

Escenario Escenario Escenario Escenario Escenario Escenario
Capa temática
1 2 3 4 5 6
Pendiente normalizado 0.1667 0.3093 0.3581 0.2 0.5058 0.1683
Red de drenaje normalizado 0.1667 * * * * *
Microcuencas normalizado 0.1667 0.1048 0.1145 0.2 0.1640 0.4944
Geología normalizado 0.1667 0.4129 0.3581 0.2 0.0829 0.1662
Geomorfología normalizado 0.1667 0.0487 0.0490 0.2 0.0462 0.0394
Distancias a los cauces 0.1667 0.1243 0.1202 0.2 0.2011 0.1317
Σ 1 1 1 1 1 1
* La capa temática no fue implementada en el modelo

Fuente: elaboración propia.
5. RESULTADOS Y DISCUSIÓN resultados entre dichos mapas, se seleccionaron

4 puntos de control en diferentes zonas de la
5.1 Análisis cualitativo de las zonas de amenaza cuenca. Un resumen de los escenarios ensayados
Con las diferentes ponderaciones ensayadas se presenta en la tabla 5, donde se enuncian las
se obtuvieron 6 mapas de amenaza; para comparar características de los mapas de amenaza generados.
Tabla 5. Características de las ponderaciones utilizadas

Escenario Características
Se emplearon las seis variables de análisis (modelo digital de elevación MDE, geología, geomorfología, red
Escenario 1
de drenaje, microcuencas, pendiente, distancia a los cauces) con un valor de ponderación igual para todos.
No se tuvo en cuenta la capa de la red de drenaje, ya que ésta incidía fuertemente en la identificación de
Escenario 2 las zonas de avenidas torrenciales y las asocia a la red hídrica. Para este caso se dio mayor importancia a
la capa de geología.
No se tuvo en cuenta la capa de la red de drenaje y se dio mayor importancia tanto a la geología como a
Escenario 3
la pendiente, por consiguiente sus pesos ponderados fueron los más altos e iguales.
No se tuvo en cuenta la capa de la red de drenaje y se le asignan pesos ponderados iguales a todas las
Escenario 4
demás capas temáticas, con el fin de identificar la influencia de la red de drenaje en el análisis de amenaza.
Se dio mayor importancia a la capa de pendientes, mientras que a la capa de geología se le dio menor
Escenario 5
importancia en la ponderación.
Se asignó mayor peso de ponderación a la capa microcuencas (índice de Melton) con el objetivo de conocer
Escenario 6
la importancia de dicha capa en la ocurrencia de avenidas torrenciales.
Fuente: elaboración propia
En la figura 3 se presentan los mapas de amenaza obtenidos con los escenarios de ponderación 1 y 6.
821000 824000 827000 83000 833000 821000 824000 827000 83000 833000
N N
W E W E
1184000
1184000
1184000
1184000
S S
1180000
1180000
1180000
1180000
1176000
1176000
1176000
1176000
1172000
1172000
821000 824000 827000 83000 833000 821000 824000 827000 83000 833000
amenaza 1 amenaza 6
0-0.2 0-0.2
0.2-0.4 0.2-0.4
0.4-0.6 0.4-0.6
0.6-0.8 0 385 770 1,540 2,310 3,080 0.6-0.8 0 385 770 1,540 2,310 3,080
Meters 0.8-1 Meters
0.8-1
(a) Escenario 1 (b) Escenario 6
Figura 3. Mapas de amenaza con los escenarios de ponderación 1 y 6

Fuente: AMVA y CTA [12]
Tabla 6. Valores de amenaza para diferentes puntos de control

en la cuenca de la quebrada Doña María
Punto de control Escenario 1 Escenario 2 Escenario 3 Escenario 4 Escenario 5 Escenario 6
Punto 1 confluencia
quebrada Astilleros
0.6384 0.6639 0.6890 0.7370 0.7450 0.4046
con la quebrada Doña
María
Punto 2 confluencia
quebrada La Manguala
0.9291 0.9192 0.9239 0.9040 0.9045 0.6549
María
Punto 3 confluencia
quebrada La Tablaza
0.7691 0.9317 0.9272 0.8510 0.8514 0.5285
María
Punto 4 parte baja
de la quebrada Doña
M a r í a , c e rc a n a s u 0.5787 0.8082 0.7961 0.5854 0.6146 0.3901
desembocadura al río
Medellín
Los resultados de la evaluación y comparación En la tabla 6 se puede observar que para el

de los valores del índice de la amenaza se presentan escenario seis, los valores de amenaza son más
en la tabla 6; es posible observar que no existe una bajos y muy diferentes en relación con los otros;
variación importante en los valores obtenidos en esto indica que el Índice de Melton (variable con
los puntos 1, 2 y 3, salvo en el caso del mapa de ame- mayor ponderación para este mapa) está disminu-
naza 6 (mayor importancia el índice de Melton), yendo considerablemente la amenaza en toda la
donde se presenta una disminución considerable cuenca.
de la cuantificación de la amenaza, indicando que Para los puntos de control 1 y 2 los valores de
el proceso de evaluación es sensible a la manipula- amenaza para cada mapa son muy similares; esto
ción de dicho parámetro. Para el caso del punto 4, indica que las zonas altas y parte de la media no
se observó una gran sensibilidad a la variación del son muy sensibles a los cambios en las variables de
nivel de importancia de los parámetros considera- análisis y a los pesos de ponderación. En los puntos
dos, y se obtuvieron los valores de amenaza más de control 3 y 4 se logra observar la influencia de
elevados, cuando se asumió la geología como el las variables de análisis empleadas en el análisis y
parámetro de mayor importancia en el análisis. El los pesos de ponderación. Los mapas de amenaza
menor valor se obtiene cuando el índice de Melton 2 y 3 presentan los valores máximos de amenaza lo
posee la mayor importancia. Se observó también, que indica que los mapas temáticos de geología y
que la red de drenaje no influye de una manera pendiente están incidiendo considerablemente en
considerable en la valoración de la amenaza, y que la ocurrencia de las avenidas torrenciales, mientras
en los puntos de control, es indiferente su inclusión que otras variables como la red de drenaje y la
o no en el proceso de análisis. geomorfología no son tan relevantes.
Como resultados finales de este cuadro, se lo que está de acuerdo con lo expuesto por Mi Río
observa que la amenaza por avenidas torrenciales et al. [14] en cuanto que las corrientes ubicadas en
está muy localizada en la parte media-baja de la la margen derecha, son las que presentan historial
cuenca de la quebrada Doña María, coincidiendo de eventos torrenciales en la microcuenca.
con las zonas donde se encuentran los depósitos Porcentaje de la Cuenca en Amenaza Alta
y donde las pendientes favorecen la depositación Mapa de Amenaza 1
de materiales. En la parte media-alta de la cuenca, 70.000
61.073
60.000
aunque pueden presentarse zonas muy pequeñas 50.000
con valores de amenaza altos, no son muy repre-
Porcentaje
40.000
sentativos y de forma general se puede decir que 30.000
la ocurrencia de avenidas torrenciales es muy baja 20.000

14,789
en comparación con la zona baja de la cuenca. 10.000 8,273
10,658
2,645 1,197 1,127 0,071
Aunado a lo anterior y tomando como base 0,079 0,089 0,000
0
0-0,09 0,1-0,19 0,2-0,29 0,3-0,39 0,4-0,49 0,5-0,59 0,6-0,69 0,7-0,79 0,8-0,89 0,9-0,99 1
las evidencias históricas evaluadas en la literatura Rangos de valores de amenaza
disponible para la zona en estudio, cotejando las

Figura 4. Histograma de distribución de pixeles
mismas con los resultados obtenidos del mapa de
de la microcuenca Doña María, para diferentes
amenaza 1, en cuanto a la distribución espacial
rangos de valoración de la amenaza, en el mapa
de las zonas con mayor valor de la amenaza, en la
de amenaza 1
escala de 0 a 1, y su posición en los rangos estudia- Fuente: AMVA y CTA [12]
dos en el histograma de la figura 6 (abscisas rangos 821000 824000 827000 830000 833000
de valores de amenaza; ordenadas, porcentaje de
pixeles pertenecientes a un rango específico), se
seleccionaron los rangos de valores de 0 a 0.5 N
1184000
1184000
como de amenaza baja (94.87% del área de la mi- W E
crocuenca), de 0.5 a 0.7 como de amenaza media
S
o moderada (3.84 % del área de la microcuenca),
y de 0.7 a 1 como de amenaza alta por avenidas
1180000
1180000
torrenciales (1.29% del área de la microcuenca),
obteniéndose así el mapa de las zonas con amenaza
por avenida torrencial, para la microcuenca de la
1176000
1176000
quebrada Doña María, resultante de la metodolo-

gía de análisis multivariado; ahí se observa que la
zona de amenaza alta la comprende la quebrada
Doña María en el tramo que fluye por la llanura
1172000
aluvial de la misma; las quebradas la Arenala, y la 821000 824000 827000 83000 833000
Tablaza, en su tramo medio-bajo; las quebradas la Amenaza por avenidas torrenciales

Isabela, la Cascada, y la Montañita (sobre la margen baja
izquierda), y las quebradas la Limona, la Cabuyala, media

alta
0 480 960 1.920 2.880 3.840
Meters
la Manguala, y el tramo final de la Justa, y de la
Muñoz, y una proporción importante de cauces Figura 5. Mapa de amenaza por avenidas
con amenaza media o moderada, en la vertiente torrenciales, metodología de análisis multicriterio
sur de la quebrada Doña María (margen derecha), Fuente: AMVA y CTA [12]
5.2 9HULÀFDFLyQGHORVUHVXOWDGRV las condiciones de fronteras determinadas para el

Con el fin de comparar la validez de los resultados modelo hidráulico de la microcuenca, empleando
se empleó un modelo de tránsito hidráulico, en estado un régimen de flujo subcrítico, el cual se evaluó
uniforme y permanente (con el modelo HEC-RAS), teniendo como referente, el no superar velocidades
con el fin de simular bajo ciertas consideraciones, del orden de los 8 m/s (velocidad media para un
el flujo de una posible avenida torrencial, en la flujo de detritos [17]).
red de drenaje que conforma la cuenca. Para dicha Para asignar el nivel de amenaza sobre la red
modelación se realizó un incremento de un 40% en de drenaje modelada, se empleó el criterio expues-
el caudal con un período de retorno de 100 años, to por López et al. [18], el cual fue empleado en la
se estimó un coeficiente de rugosidad de Manning región de Cerro Grande en Venezuela, teniendo
de acuerdo con los valores de este coeficiente [7] como base, los datos colectados a partir de la tra-
y empleados en China para modelar avalanchas. gedia provocada por los flujos de detritos ocurridos
Dichos valores están en función del tipo de flujo, en el Estado de Vargas en 1999. Se implementó la
de las características del cauce (material de lecho propuesta para un período de retorno del orden
y alineamiento del mismo), de la pendiente media de los 100 años, donde para una intensidad del
del cauce, y de la profundidad del flujo. Para el caso fenómeno alta, se obtiene un nivel de amenaza
de la microcuenca de la quebrada Doña María, se alto, para una intensidad del fenómeno media, se
emplearon los datos correspondientes a un flujo obtiene un nivel de amenaza medio, y para una
con predominio de partículas gruesas, para un intensidad del fenómeno baja, se obtiene un nivel
canal con curvas y gradas, con un material del lecho de amenaza bajo. La intensidad del fenómeno
conformado por piedras de diámetro medio entre en la zona catalogada como la llanura aluvial de
0.3 y 0.5 m, con pendientes medias del cauce entre el la quebrada Doña María se estimó a partir de la
8% y el 15%. El valor de la profundidad de flujo, se tabla 7, empleando los ráster de profundidad y
obtuvo de una corrida preliminar del modelo, para velocidad para dicha zona. Para el resto de la red
las condiciones propuestas para estudiar la creciente de drenaje, se determinó como zona de amenaza
con período de retorno de 100 años, en lo referente alta, la correspondiente al dominio de la mancha
a la asignación por la amenaza por inundación. Así, de inundación, para las condiciones especiales de
se asumieron valores del coeficiente de rugosidad, modelación, siguiendo así, el criterio expuesto en
del orden de los 0.1, 0.125, 0.167, y 0.25, para Mi Río et al. [14] para la determinación del nivel
profundidades cercanas a los 0.5 m, 1 m, 2 m, y 3 de amenaza alta ante inundación, para un período
m, respectivamente. Adicionalmente, se emplearon de retorno de 100 años.
Tabla 7. Valoración de la intensidad de un flujo de lodo o de un flujo de escombros
Profundidad máxima Producto de la velocidad máxima (v),

Grado de intensidad Condicional
(h en metros) y la profundidad máxima (h) (m2/s)
Alta >1 ó vh > 1
Media 0.2 < h < 1 y 0.2 < vh < 1
Baja 0.2 < h < 1 y 0.2 < vh < 1
Se evaluó el nivel de amenaza, considerando que en la llanura aluvial los datos coincidieron más
el dominio de la mancha de inundación, como la o menos bien, pero que el modelo hidráulico no
zona de amenaza alta por avenidas torrenciales, representaba las zonas presentadas por el modelo
excepto en el tramo correspondiente a la llanura multicriterio para las zonas de montaña, pues el
de inundación de la quebrada Doña María, donde primero consideró solamente la inundación gene-
se utilizó el criterio de comparación de la profundi- rada y no los deslizamientos asociados. La figura
dad máxima del flujo versus el producto de dicha 7 presenta una combinación de los resultados,
profundidad, con la velocidad máxima del flujo, donde en la zona de montaña, y de piedemonte,
propuesto por López et al. [18] y se obtuvo, así, se acogió la evaluación obtenida a partir de la eva-
el mapa de amenaza por avenidas torrenciales, a luación multicriterio, y en la zona correspondiente
partir de la modelación. En la figura 6 se presentan a la llanura aluvial de la quebrada Doña María, se
los resultados obtenidos. agregaron ambos resultados, teniendo como crite-
rio el dominio de la amenaza alta, sobre la amenaza
820000 824000 828000 832000
media y baja; así, se obtuvo el mapa (a asignación)
final de la amenaza por avenida torrencial.
1186000
1186000
1183000
1183000
1180000
1180000
1177000
1177000
Riesgo Avenidas torrenciales

Nivel de amenaza
Alto
1174000
1174000
Bajo
Medio
1171000
1171000
820000 824000 828000 832000
Amenaza por
avenidas torrenciales
0 950 1.900 3.800 5.700 7.600
alta Meters N
W E
Figura 6. Zonas de amenaza por avenidas
S Mapa de amenaza
torrenciales, a partir del modelo de tránsito 0 0.5 1 2 3 4 por avenida torrencial
hidráulico HEC-RAS Kilómetros
Fuente: AMVA y CTA [12]
5.3 Combinación de los resultados Figura 7. Asignación final de las zonas de

En la figura 7 se presenta una combinación amenaza por avenidas torrenciales, en la
de los resultados de la evaluación multicriterio y microcuenca de la quebrada Doña María.
del modelamiento hidráulico. En esta se puede ver Fuente: AMVA y CTA [12]
De la figura 7 hay que mencionar que la timbres de Itagüí. En la parte baja de este afluente,
metodología asignó como amenaza alta algunos se aprecia una cantidad de rocas de diámetros del
sitios en donde no se poseen registros históricos orden de los 0.5 m (que no corresponden a residuos
de avenidas torrenciales (especialmente en las que- del lavado de la matriz, figura 8a), que por su forma,
bradas la Isabela, La Montañita, y la Cascada), sin tamaño, textura, estructura, etc., son indicativas de
embargo, dichos resultados pueden ser evidencia la gran competencia de la quebrada, y dan pie a
de zonas que muestran un potencial importante pensar que dichas partículas pudieron haber sido
a desencadenar y sufrir los efectos de eventos de depositadas por eventos del tipo torrencial. En la
tipo torrencial, que deben ser corroborados por quebrada la Limona, la cual se encuentra en la
estudios de geomorfología fluvial y estratigráficos parte media de la microcuenca, sobre la margen
antes de planificar el uso de dichas zonas, pues en derecha, de la quebrada Doña María, se observó la
ellas se presenta la unión de factores que propician presencia de rocas de tamaño considerable (figura
la aparición de eventos torrenciales. 8b) lo que posiblemente indica que estas han sido
transportadas por eventos extraordinarios, tal vez
5.4 Recolección de evidencias en campo del tipo torrencial, apoyado esto en las evidencias
Con el fin de constatar en campo las tenden- reportadas de depósitos aluvio-torrenciales en las
cias mostradas por los modelos, se realizó un reco- proximidades del cauce activo de la quebrada.
nocimiento a las zonas evaluadas como de amenaza Del trabajo de campo antes descrito, se evi-
alta por el análisis multicriterio, estas fueron las denció de forma general, que dichas corrientes
quebradas la Arenala, la Tablaza, la Limona, y la presentan indicios sobre el cauce activo, evidencias
Manguala. de un comportamiento torrencial, y son candi-
La quebrada la Tablaza se encuentra en la parte datas a ubicarse en el rango de la amenaza alta a
baja de la microcuenca de la quebrada Doña María, moderada, comportamiento similar al obtenido
y desemboca a ésta en cercanías de la empresa Cur- previamente.
Figura 8. (a) Quebrada la Tablaza, detalle material depositado en el lecho del cauce.
(b) Quebrada la Limona, detalle material del lecho
Fuente: Elaboración propia
6. CONCLUSIONES [3] G. Barreto, and D. Ruiz, Estudio de la amenaza por

avenidas torrenciales en el sector noroccidental de la ciudad
Se ha presentado la aplicación de una me- de Ibagué (departamento de Tolima), Ibagué, 2001.
todología cualitativa de valoración de zonas de [4] I. Piedrahíta, and M. Hermelín, «La avenida torrencial
amenaza por avenidas torrenciales, en la cual me- del río Tapartó (Antioquia) 1993,» Desastres de origen
diante ponderación de diferentes factores como la natural en Colombia 1979-2004, M. Hermelín, ed.,
topografía, la pendiente, la distancia a los cauces, pp. 109-120, Medellín: Universidad Eafit, 2005.
el flujo acumulado, la geología y la geomorfología [5] M. Hermelín, O. Mejía, and E. Velásquez, «El evento
y el índice de Melton para la microcuenca, se iden- torrencial de septiembre de 1990 en San Carlos (Co-
tifican zonas de amenaza a avenidas torrenciales. lombia) y sus consecuencias catastróficas,» Desastres
naturales en América Latina, J. Lugo and M. Inbar, eds.,
A pesar que se contaba con poca información
pp. 335-348, México: Fondo de Cultura Económica,
para el estudio de la quebrada Doña María, fue
2002.
posible realizar un análisis en el que se obtuvieron [6] H. Flores, R. García, K. Bukor et al., Estudio geomor-
resultados acordes a las evidencias históricas, y a las fológico del evento de aludes torrenciales de 1999 en la
condiciones actuales de los cauces estudiados en cuenca experimental de San José de Galipán, Caracas:
corrientes como la Limona, la Manguala, la Justa, Universidad Central de Venezuela, 2004.
y la Cabuyala, en lo que a asignación del nivel de [7] J. Suárez, Control de Erosión en Zonas Tropicales, Buca-
amenaza se refiere. ramanga: Universidad Industrial de Santander, 2001.
[8] M. Mironova, and V. Yablonsky, «A mathematical
model of shear debris flow,» Erosion, debris flows and en-
7. AGRADECIMIENTOS
vironment in mountain regions (Proceedings of the Chengdu
Al Centro de Ciencia y Tecnología de Antio- Symposium, July 1992): IAHS, 1992.
quia (CTA) a través del convenio interinstitucio- [9] USGS, Geological Survey Water-Supply Paper 2275, 1984.
nal Cátedra del Agua, quienes fueron los líderes [10] L. F. Castillo Navarro, “Aplicación de un modelo nu-
del proyecto “Plan de ordenación y manejo de la mérico de flujos de escombros y lodo en una quebrada
en el Perú,” Facultad de Ingeniería Civil, Universidad
microcuenca de la quebrada Doña María” que
Nacional de Ingeniería, Lima, 2006.
permitió la realización de este trabajo. A las otras
[11] C. OSSO. «Escudriñando en los desastres a todas las
entidades involucradas en el proyecto, particular-
escalas,» agosto, 2008; http://www.osso.org.co/docu/
mente al Área Metropolitana del Valle de Aburrá. publicac/1999/escudrinhando/definiciones.pdf.
[12] AMVA, and CTA, «Plan de Ordenación y Manejo de
REFERENCIAS la microcuenca de la quebrada Doña María, municipio
de Itagüí. (en proceso de edición),» 2008.
[1] J. González, O. Chávez, and M. Hermelín, «Aspectos [13] AMVA, and UNAL, Estudio para la reglamentación del
geomorfológicos de la avenida torrencial del 31 de aprovechamiento, uso y administración de las aguas de la
enero de 1994 en la cuenca del río Fraile y sus fenóme- quebrada Doña Maria, incluidas todas las subcuencas o
nos asociados,» Desastres de origen natural en Colombia microcuencas de sus quebradas afluentes, hasta su desem-
1979-2004, M. Hermelín, ed., Medellin: Universidad bocadura en el río aburra, localizada entre los Municipios
Eafit, 2005. de la Estrella, Itagüí y Medellín., 2007.
[2] C. Polanco, and G. Bedoya, “Compilación y análisis de [14] MIRIO, CORANTIOQUIA, and UNAL, «Metodo-
los desastres naturales reportados en el departamento logía para la formulación de planes de ordenación y
de Antioquia exceptuando los municipios del Valle manejo de microcuencas (PIOM),» 2002, p. 430.
de Aburrá - Colombia, entre 1920-1999,” Ingenieria y [15] C. Ordóñez, and R. Martínez, Sistemas de Información
ciencia, vol. 1, no. 1, pp. 45-65, 2005. Geográfica, Aplicaciones prácticas con IDRISI KILI-
MANJARO al análisis de riesgos naturales y problemáticas [17] D. Vandine, Debris flow control structures for forest Engi-
medioambientales, Madrid: RA-MA, 2003. neering.: Ministry of Forest, British Columbia, 1996.
[16] A. Iroume, “Transporte de sedimentos en una cuenca [18] J. López, R. Garcia, M. Bello et al., “Analysis and Map-
de montaña en la Cordillera de los Andes de la Nove- ping of Debris Flow in Alluvial Fans (Part I: Methodo-
na Región de Chile,” Bosque (Valdivia), vol. 24, no. 1, logy Description).” in XX Congreso latinoamericano
pp. 125-135, 2003. de Hidráulica, La Habana, Cuba, 2002.

Avenidas Torrenciales Metodo PDF

Cargado por

Copyright:

Formatos disponibles

Avenidas Torrenciales Metodo PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Avenidas Torrenciales Metodo PDF

Cargado por

Copyright:

Formatos disponibles

Revista Ingenierías Universidad de Medellín

EVALUACIÓN DE ZONAS DE AMENAZA POR AVENIDAS

Luis Javier Montoya Jaramillo*

EVALUATION OF HAZARD ZONES BY TORRENT FLOW

La cuenca de la quebrada Doña María está Tabla 1. Parámetros morfométricos de la

Se tienen registros de eventos con caracterís-

y se procedió a calificar cada una de sus caracte-

Tabla 2. Funciones matemáticas empleadas en la Para la ocurrencia de avenidas torrenciales

Tabla 3. Escala de Valoración de pesos relativos

Fuente: Valores establecidos en el software IDRISI Kilimanjaro, [15]

Tabla 4. Valores de ponderación empleados

* La capa temática no fue implementada en el modelo

5. RESULTADOS Y DISCUSIÓN resultados entre dichos mapas, se seleccionaron

Tabla 5. Características de las ponderaciones utilizadas

(a) Escenario 1 (b) Escenario 6

Figura 3. Mapas de amenaza con los escenarios de ponderación 1 y 6

Tabla 6. Valores de amenaza para diferentes puntos de control

Fuente: elaboración propia

Los resultados de la evaluación y comparación En la tabla 6 se puede observar que para el

la ocurrencia de avenidas torrenciales es muy baja 20.000

disponible para la zona en estudio, cotejando las

quebrada Doña María, resultante de la metodolo-

Tablaza, en su tramo medio-bajo; las quebradas la Amenaza por avenidas torrenciales

izquierda), y las quebradas la Limona, la Cabuyala, media

5.2 9HULÀFDFLyQGHORVUHVXOWDGRV las condiciones de fronteras determinadas para el

Tabla 7. Valoración de la intensidad de un flujo de lodo o de un flujo de escombros

Profundidad máxima Producto de la velocidad máxima (v),

Fuente: elaboración propia

Riesgo Avenidas torrenciales

820000 824000 828000 832000

Fuente: AMVA y CTA [12]

5.3 Combinación de los resultados Figura 7. Asignación final de las zonas de

6. CONCLUSIONES [3] G. Barreto, and D. Ruiz, Estudio de la amenaza por

También podría gustarte