Simposio Latinoamericano de Control de Erosión 2.002

MODELO PARA EL PRONOSTICO DE LA AMENAZA POR DESLIZAMIENTOS EN TIEMPO REAL

Agrólogo REINALDO SÁNCHEZ LÓPEZ

Ing. Civil LUIS FERNANDO URREGO LARA

Instituto de Hidrología, Meteorología y Estudios Ambientales, IDEAM, Colombia

Contratista, IDEAM, Colombia Geól. GERMÁN VARGAS CUERVO

Geól. RUTH MAYORGA MÁRQUEZ Ex funcionario, IDEAM, Colombia Instituto de Hidrología, Meteorología y Estudios Ambientales, IDEAM, Colombia

RESUMEN

En este documento se presenta una síntesis del modelo de susceptibilidad general del terreno a los deslizamientos y el modelo de pronóstico de la amenaza por deslizamientos detonados por lluvias, en tiempo real, que el IDEAM ha implementado para escalas generales. De igual manera se describen algunos aspectos metodológicos para el tratamiento de la información fuente y la espacialización del modelo de susceptibilidad del terreno a deslizamientos así como la lluvia crítica detonante. Finalmente se describe el proceso de automatización, captura y tratamiento de información generada por el modelo de susceptibilidad del terreno, la lluvia detonante, la salida de productos y sus aplicaciones para los diferentes sectores del país.

INTRODUCCIÓN El Instituto de Hidrología, Meteorología y Estudios Ambientales, Ideam, en el marco de sus funciones y como nodo central del Sistema de Información Nacional Ambiental y miembro del comité técnico nacional de apoyo al Sistema Nacional para la Atención y Prevención de Desastres, está implementado un modelo para el pronóstico de la amenaza por "deslizamientos" detonados por lluvias en tiempo real.

Los términos más utilizados para denominar globalmente estos fenómenos son: movimientos en masa, movimientos de terreno, deslizamientos, deslizamientos de terreno, deslizamientos de vertiente, movimientos de ladera, derrumbes, volcanes, etc. En general se definen nueve grupos principales: Caídas, flujos, deslizamientos, volcamientos, propagación lateral, hundimientos, reptación, movimientos complejos y avalanchas (Ideam 2001). ANTECEDENTES Y JUSTIFICACIÓN

El término "deslizamientos" en este trabajo se utiliza como genérico para denominar los movimientos en masa y hace referencia a desplazamientos de terreno y no incluye los procesos erosivos denudacionales. Un movimiento en masa puede ser definido como todo desplazamiento hacia abajo (vertical o inclinado en dirección del pie de una ladera) de un volumen de material litológico (roca, formación superficial o suelo) importante, en el cual el principal agente es la gravedad y que puede o no incluir el efecto del agua.

Colombia es un país que por sus características geográficas, geológicas, geomorfológicas, climáticas y de usos del suelo hacen que los deslizamientos, sean de gran impacto y se presenten con mucha frecuencia. Según estudios realizados (Castellanos, 1996) las causas de los eventos desastrosos (inundaciones, deslizamientos, avalanchas, etc.) son en un 96% meteorológicas. El mencionado autor concluye que el 56% de las causas son por lluvias prolongadas o persistentes, el 37 % por lluvias intensas de corta duración y en menor porcentaje los huracanes.

Simposio Latinoamericano de Control de Erosión 2.002

Estudios recientes realizados por el Ideam, de eventos extremos de variabilidad climática (El Niño – La Niña), se observa como se incrementa la frecuencia y la intensidad del daño de los deslizamientos, especialmente en eventos húmedos haciéndose evidente la relación lluvia deslizamiento (Ideam, 2001). Dado el problema que representan estos fenómenos en Colombia, se han adelantado investigaciones sobre la ocurrencia de los deslizamientos, estabilidad de taludes, realización de inventarios, sistemas de monitoreo y recientemente modelos de pronóstico que buscan hacer explícita tanto espacial como temporalmente la amenaza relativa que representan los deslizamientos. METODOLOGÍA El modelo esta basado en cuatro fases principales: a) caracterización y zonificación de la susceptibilidad de los terrenos a los deslizamientos, b) análisis de los eventos históricos para la estimación de las lluvias detonantes c) caracterización y zonificación de umbrales de lluvias criticas detonantes d) automatización del modelo para el pronóstico y alertas basado en mapas dinámicos de amenaza. Para la fase de zonificación de la susceptibilidad del terreno a deslizamientos se implementó un modelo heurístico semicuantitativo con doce variables o componentes de terreno: litología, densidad de fracturamiento, inclinación de la pendiente, índice de relieve relativo, Sistemas Morfogénicos, morfodinámica (procesos erosivos), suelos y su relación con la dinámica del agua (infiltración, capacidad de almacenamiento y regulación de humedad), suelos y su constitución (textura, materia orgánica, naturaleza y tipo de las arcillas como coloides del suelo) y coberturas vegetales. Para el análisis de la lluvia detonante se evaluó dos efectos: la lluvia precedente a la ocurrencia del deslizamiento, relacionados a los eventos históricos tipo y la lluvia de corta duración o del día del evento (con análisis de Intensidad, duración y frecuencia IDF). La caracterización y espacialización de los umbrales se esta implementando mediante el modelo de la lluvia precedente (numero de días y cantidad) y mapa de subregiones climáticas, el modelo y mapa de intensidades en 24 horas. La metodología para la caracterización y zonificación de umbrales de lluvias críticas detonantes está basada en el cálculo de lluvia acumulada en distintos periodos, clasificados en intervalos, según su nivel de

probabilidad estadística. Con el procesamiento de la información pluviométrica diaria de la estación más cercana al sitio de ocurrencia del evento, se sistematizan los datos diarios de tal forma que permitan el análisis de probabilidades por medio de la distribución Gumbel (distribución que más se ajusta, porque fue desarrollada para su aplicación en problemas de valores extremos). La distribución de Gumbel permite encontrar la probabilidad de excedencia o el periodo de retorno para cualquier valor a partir de los datos de la serie histórica, calculando para cada valor de la serie histórica la probabilidad de excedencia utilizando la ecuación de Weibull (Castellanos, 1996). FASES DEL MODELO Zonificación de la susceptibilidad del terreno a deslizamientos La estabilidad de un suelo es una cualidad que expresa su mayor o menor tendencia a permanecer in situ. La cualidad opuesta es la inestabilidad y se define como la tendencia al desplazamiento pendiente abajo o como grado de susceptibilidad a procesos morfodinámicos como los movimientos en masa y la erosión (WAY, 1973). Los parámetros o variables ambientales que se tuvieron en cuenta para el análisis y zonificación de la susceptibilidad a los deslizamientos son de orden Geológico, geomorfológico, edafológico y de coberturas vegetales. Geología (Litología) La geología es una de las principales variables ambientales que intervienen en la estabilidad de las laderas o taludes. Esta puede ser valorada a partir de varias propiedades como: composición mineralógica o tipo de material (roca o formación superficial), textura, estructura, grado de meteorización y grado de fracturamiento. Para este alcance del modelo el factor litología se evaluó a partir de la caracterización de la información cronoestratigráfica (roca - tiempo) presentada en el en el mapa geológico de Colombia (Ingeominas 1988). Fracturamiento y Meteorización Teniendo en cuenta que este análisis es regional, el grado de fracturamiento de los materiales litológicos se estimó a partir de la cartografía de fallas

Simposio Latinoamericano de Control de Erosión 2.002

geológicas del mapa Ingeominas (1988).

geológico

nacional

de

Como criterio de valoración de este parámetro se tomó la susceptibilidad directamente proporcional a la densidad de fracturamiento.

estimación de la estabilidad y susceptibilidad de los suelos a procesos morfodinámicos (movimientos en masa) se realizó con base en los siguientes dos criterios: Estructura y dinámica del agua.

Morfología

Estructura y composición

La génesis y forma del relieve constituye en algunos casos un factor de estabilidad de las laderas o taludes. Para la valoración de la susceptibilidad de este parámetro se utilizó la información del estudio de Sistemas Morfogénicos de Colombia (Ideam Unal, 1996).

Un suelo es más estable y menos susceptible a los movimientos en masa mientras su estructura esté garantizada por una buena aglutinación o cementación, propiciada por los coloides del suelo llámense materia orgánica, arcillas, carbonatos, óxidos e hidróxidos de hierro y aluminio entre otros. La buena estructuración del suelo permite una buena aireación, una buena permeabilidad y una mayor resistencia al corte.

Indice de relieve relativo Este indicador representa la rugosidad del terreno involucrando el aspecto de pendiente de laderas. El índice de relieve relativo se define como la mayor diferencia de altura en un área determinada. (IRr = dHmax/A en m/km2). Densidad de drenaje Un gran número de movimientos en masa ocurre en áreas de mayor densidad de drenaje. La densidad de drenaje define el grado de disección de una ladera o un terreno cualquiera por unidad de área. Este grado de disección representa áreas de mayor o menor susceptibilidad a los movimientos en masa. El mapa de densidad de drenaje se elaboró a partir del mapa de red de drenaje digital a escala 1:500.000 tomando una malla de 100 km2. Los valores de susceptibilidad son directamente proporcionales con el incremento de la densidad. Suelos Para obtener la susceptibilidad de la variable suelos se tuvo en cuenta de manera cualitativa variables como materia orgánica, tipo de arcilla, capacidad de regulación y de almacenamiento de agua en los suelos. La estabilidad del suelo está controlada por condiciones intrínsecas de orden químico, físico y biológico como: la textura, la estructura, los coloides (calidad mineralógica y física de las arcillas y de la materia orgánica), la profundidad del perfil, la densidad aparente, la porosidad, y la plasticidad entre otras, condiciones que determinan la resistencia al corte y dinámica de factores externos como el agua, la pendiente, el uso y la cobertura. La

La información de los coloides y aglutinantes del suelo se tomó de la información de los mapas de arcillas y materia orgánica (IGAC, 1987). Las arcillas como aglutinantes se evaluaron a partir de su mineralogía y para cada uno de los tipos sus características químicas como cuerpos que poseen una capacidad de intercambio y atracción entre partículas mediante: fuerzas electrostáticas Coulombicas con interacciones arcilla - arcilla, arcilla - óxidos y arcillas - humus y las fuerzas de atracción y dispersión (fuerzas de Van der Waal’s). Desde el punto de vista físico las arcillas se interpretan desde su cualidad de consistencia y plasticidad, considerada como los umbrales mediante el cual los minerales arcillosos pasan de estado sólido - plástico y liquido. Dinámica del agua Un segundo criterio para la evaluación de la estabilidad o susceptibilidad de los suelos, es la dinámica del agua en el perfil, la cual es controlada por las propiedades fisicoquímicas y biológicas y que a diferencia de la erosión, no actúa un agente de transporte (agua, hielo, viento); sin embargo el agua se encuentra íntimamente asociado a la estabilidad inicialmente incrementando el peso de la masa a desplazarse y al ayudar al flujo descendente, una vez iniciado el movimiento en masa, ya que estrecha el límite plástico y en consecuencia el proceso de fluidez. Igualmente ocasiona un aumento de la presión del aire contenido en los poros y ejercen una fuerza que se manifiesta en desplazamiento de los materiales. El suelo propiamente dicho es uno de los receptores y reguladores del agua, factor muy dinámico que altera los estados de equilibrio relativo

Simposio Latinoamericano de Control de Erosión 2.002

o de estabilidad. La dinámica del agua en el suelo se puede estimar mediante la capacidad de almacenamiento y regulación. Dentro del ciclo hidrológico, el agua que pasa por los suelos puede: infiltrase, escurrirse, ascender y almacenarse, de acuerdo con las características fisicoquímicas y biológicas. Para el análisis de estabilidad y en lo relacionado con los movimientos en masa, estas dinámicas se interpretan a partir de los indicadores de capacidad de almacenamiento y retención de humedad, considerando el agua como una masa que incrementa el peso y por lo tanto facilita los movimientos en masa. Existen mecanismos por los cuales el agua es retenida y almacenada en el suelo y están relacionadas a fuerzas derivadas de las interfases líquido - aire o sólido - líquido del sistema suelo agua. Algunos minerales del suelo como las arcillas pueden almacenar el agua dentro de su estructura cristalina y provoca la expansión y contracción volviéndolas muy dinámicas ante los cambios de humedad. Para obtener el mapa de susceptibilidad del parámetro suelo se calificaron de 1 a 5 cuatro variables: materia orgánica, tipo de arcilla, capacidad de regulación y capacidad de almacenamiento de agua en los suelos.

Cobertura vegetal y uso El tipo, la densidad, la capacidad de interceptación, el área de protección de cobertura vegetal constituye un factor de resistencia o favorecimiento de procesos morfodinámicos como la erosión y los movimientos en masa. En áreas de alta pendiente y de coberturas vegetales de ciclos biológicos muy lentos y frágiles, toda intervención de estas, hace que se aumente la susceptibilidad. Las coberturas vegetales fueron evaluadas cualitativamente a partir del estudio de coberturas vegetales de Colombia (Ideam 2000), teniendo en cuenta las variables de cobertura foliar, capacidad de almacenamiento hídrico y profundidad radicular la cual se relaciona con la estructura (arbórea, arbustiva, rastrojos, pastos, etc.). Determinación del peso de los factores Teniendo en cuenta que cada parámetro o factor de terreno influye diferentemente sobre la estabilidad o susceptibilidad de los terrenos a los movimientos en masa, se analizaron espacial y digitalmente los diferentes mapas de susceptibilidad temáticos y se establecieron los siguientes pesos en porcentaje. Tabla 1 Parámetro o factor

Peso en %

morfodinámica Este parámetro muestra la dinámica de las laderas con procesos de erosión y o remoción en masa activos y/o observados o que presentan señales que permitieron identificar su génesis y actividad histórica. Se establecieron cinco categorías de erosión de muy baja o nula a muy alta y se establecen valores de susceptibilidad directamente proporcionales con el grado de intensidad (Ideam, 2000). Grado de Inclinación y longitud de la pendiente Las pendientes es uno de los principales factores dinámicos y particularmente de los movimientos en masa, ya que determinan la cantidad de energía cinética y potencial de una masa inestable, por lo cual se calificó la susceptibilidad con los intervalos de pendiente; aumentando la susceptibilidad a mayor grado de inclinación de la pendiente en función de su energía cinética y potencial.

L

Litología

15

Df

Densidad de fracturamiento

15

M

Morfología

10

Dd

Densidad drenaje

S

Suelos

10

P

Pendientes

15

Ie

Morfodinámica (intensidad de erosión)

15

Cv

Cobertura vegetal

10

Total

de

10

100

Simposio Latinoamericano de Control de Erosión 2.002

Tabla 1. Valores de ponderación de los componentes o variables del terreno (Ideam 2001). ZONIFICACIÓN DE LA SUSCEPTIBILIDAD GENERAL DEL TERRENO Una vez establecidos semicuantitativamente los diferentes valores de susceptibilidad y los pesos de los parámetros o variables ambientales se elaboró el mapa de susceptibilidad general del terreno a los movimientos en masa mediante una análisis estadístico multivariado, en el cual son sumados digitalmente los ocho mapas de susceptibilidad temática. Este procedimiento se realizó mediante la utilización del SIG, utilizando el siguiente algoritmo: ST= StL (0.15)+ StDf (0.15)+ Stm (0.10)+ StDd (0.10)+ Np Sts (0.10) + Stp (0.15)+StIe (0.15)+ StCv (0.10) Np

De donde: ST =

Susceptibilidad total del terreno

StL= Susceptibilidad del parámetro litológico StDf = Susceptibilidad del parámetro densidad de fracturamiento

la metodología para determinar la relación lluvia – deslizamiento, consideran que existen dos tipos de efectos principales de la lluvia: “El efecto acumulativo de la lluvia de largo plazo, varios días, y el efecto inmediato de la lluvia de corto plazo, que se presenta justo antes de la ocurrencia del deslizamiento”. Para el efecto de este estudio se asume la lluvia de largo plazo de 1 a 180 días antes del deslizamiento dependiendo del sitio, y la de corto plazo del orden de 1 a 24 horas. Análisis de la lluvia antecedente – larga duración Este método el cual aporta un elemento en la predicción de movimientos en masa basado en la precipitación acumulada o antecedente de semanas previas a la ocurrencia del evento y difiere según las condiciones locales de los suelos y del régimen climático, por lo que el método debe ser aplicado a cada sitio en estudio. El estudio consiste en establecer relaciones estadísticas de la lluvia con el deslizamiento, cuantificando los componentes de intensidad, duración y frecuencia con el fin de sintetizar un modelo de comportamiento de la zona por estudiar y así contribuir al pronóstico de los mismos; para esto se seleccionan los sitios con eventos de deslizamientos causados por lluvia, que cuenten con buena información pluviométrica.

Stm = Susceptibilidad del parámetro morfología StDd = Susceptibilidad del parámetro densidad de drenaje Sts = Susceptibilidad del parámetro suelos Stp =

Susceptibilidad del parámetro pendientes

StIe = Susceptibilidad del parámetro morfodinámica (intensidad de erosión) StCv = Susceptibilidad cobertura vegetal.

del

parámetro

NP = Número de parámetros CARACTERIZACIÓN UMBRALES DE DETONANTES

Y

ZONIFICACIÓN DE LLUVIAS CRÍTICAS

Evaluación de la lluvia como evento detonante de deslizamientos De acuerdo con los estudios de La The Japanese Geotechnical Society (1997) y D´Orsi et al. (1987), en

Para cuantificar las relaciones de intensidad duración o volúmenes de agua relacionados con la ocurrencia de deslizamientos, se debe contar con series históricas de lluvia completas, con un buen período de registro a nivel diario (en lo posible mayores a 15 años). Estas series se analizan estadísticamente utilizando métodos de inferencia estadística y aplicando la teoría de las probabilidades. La metodología está basada en el cálculo de lluvia acumulada en distintos períodos, clasificados en intervalos, según su nivel de probabilidad estadística. Con el procesamiento de la información pluviométrica diaria de la estación más cercana al sitio de ocurrencia del evento, se sistematizan los datos diarios de tal forma que permitan el análisis de probabilidades por medio de la distribución Gumbel (distribución que más se ajusta, porque fue desarrollada para su aplicación en problemas de valores extremos). La distribución de Gumbel permite encontrar la probabilidad de excedencia o el periodo de retorno para cualquier valor a partir de los datos de la serie histórica, calculando para cada valor de la

Simposio Latinoamericano de Control de Erosión 2.002

serie histórica la probabilidad de excedencia utilizando la ecuación de Weibull (Castellanos, 1996). El método comprende la elaboración de la curva regresiva de lluvias acumuladas cuyo procedimiento general a seguir es según Castellanos (1996); Castellanos y González (1996, 1997) y Wieczorek et al. (1987), el siguiente: ? ? A partir de la serie histórica de precipitación diaria regresiva para el evento se elabora una serie de sumas acumuladas regresivas iniciando desde el día en que ocurrió el deslizamiento hasta los 3 a 6 meses anteriores (dependiendo del análisis de cada caso).

volúmenes de agua que han caído como precipitación en determinadas épocas y su relación con la ocurrencia de deslizamientos. Entonces se puede establecer que para los valores superiores a la lluvia crítica se aumenta la posibilidad de ocurrencia de deslizamientos en la región de estudio. AUTOMATIZACIÓN DEL MODELO PARA EL PRONÓSTICO DE DESLIZAMIENTOS BASADO EN MAPAS DINÁMICOS DE AMENAZA. Mediante programación en Arc/Info macro leguaje (aml), se procesa la componente de susceptibilidad en formato vector asociada a la componente detonante (lluvia).

? ? Se gráfica una curva en forma regresiva para el evento a partir de lluvias críticas, tomando como origen el día en que ocurrió el deslizamiento. Desde el día en que ocurrió el movimiento con un valor de lluvia hasta el día en que no hay lluvia, para este momento se lee la precipitación acumulada en mm, y la duración en días (umbral de precipitación). De la curva se pueden obtener varios umbrales de precipitación (es decir, después de períodos cortos de no lluvia se puede reactivar otro periodo de lluvia) los cuales dependen de la intensidad de lluvia de la zona.

La componente lluvia se toma automáticamente de la base datos de la red hidrometeorológica y se procesa espacialmente mediante procesos de interpolación "Kriging y Quintic" de Arc/Info, con la que se obtiene la cobertura de lluvia antecedente, detonante y pronosticada.

? ? Construcción de series parciales de la información diaria (con el periodo completo de registro) con los valores de la duración de los umbrales: es una serie de sumas acumuladas para cada duración (D), es decir que se suma la serie de lluvia diaria en intervalos según cada duración.

Modelo de susceptibilidad deslizamientos

? ? De las anteriores series parciales se construye series de máximos multianuales. ? ? Obtención de umbrales críticos con sus respectivos periodos de retorno y cálculo de la relación lluvia - deslizamiento para cada fenómeno de remoción en masa: con los umbrales de lluvia se determinan los periodos de retorno de cada uno de ellos a partir de distribuciones de densidad de probabilidad para valores extremos y luego se toma el umbral más crítico. Se toma como la lluvia acumulada crítica representativa de cada estación la de mayor periodo de retorno, pues es el evento menos frecuente. En resumen se tiene que con una buena información pluviométrica en la zona, es posible cuantificar la relación intensidad – duración de lluvia o los

El mapa de amenaza se obtiene mediante un algoritmo aplicado a las dos componentes anteriormente descritas. RESULTADOS del

terreno

a

los

Como resultado del modelamiento se obtuvo el mapa nacional de la susceptibilidad del terreno a los deslizamientos a escala 1:500.000 (mapa 1), y se puede establecer que el 6,7% del territorio nacional presenta susceptibilidad muy alta, 6.9% susceptibilidad alta, el 7.2% moderada, el 7.7% baja, el 8.1% muy baja y el 63.5% presenta una susceptibilidad nula. Susceptibilidad del terreno a los deslizamientos

6.70%

7.00% 7.06% 7.61%

63.54%

Muy Alta

8.10%

Alta

Moderada

Baja

Muy Baja

Nula

Mar Caribe

Providencia

San Andrés

Mapa 1. Susceptibilidad del terreno a los deslizamientos

Panamá

Venezuela

Océano Pacífico

Ecuador Brasil

LEYENDA Susceptibilidad nula Susceptibilidad muy baja Perú

Susceptibilidad baja Susceptibilidad moderada Susceptibilidad alta Susceptibilidad muy alta

Kilometros

Fuente: IDEAM, 2002

200

0

Simposio Latinoamericano de Control de Erosión 2.002

Modelo de lluvia detonante de deslizamientos

REFERENCIAS

Según los análisis estadísticos de la lluvia antecedente a los deslizamientos tipo y las estaciones hidrometeorológicas representativas y de influencia se propone para Colombia el modelo de lluvia crítica para deslizamientos, mediante la ecuación:

CASTELLANOS J., R. (1996).- Lluvias Críticas en la Evaluación de Amenazas de Eventos de Remoción de Masa - Tesis de Magister en Geotecnia Universidad Nacional de Colombia - Departamento de Ingeniería Civil. Santafé de Bogotá.

Y = 25.265 ln (x) -95.569 R2 = 0.4369 Modelo de pronóstico de deslizamientos Mediante el procesamiento automático de la información de los modelos de susceptibilidad y de lluvia crítica (antecedente, detonante y esperada) se elaboran mapas instantáneos que permiten la espacialización del pronóstico de la amenaza por deslizamientos.

CASTELLANOS J., R. Y GONZALEZ G., A. (1996): Relaciones entre la Lluvia Anual y la Lluvia Crítica que Dispara Movimientos en Masa. IX Jornadas Geotécnicas de la Ingeniería Colombiana. pp. 4.62 – 4.70. Santa Fe de Bogotá. CASTELLANOS J., R. Y GONZALEZ G., A. (1997): Algunas Relaciones de Precipitación Crítica – Duración de Lluvias que Disparan Movimientos en Masa en Colombia. 2ed. Panamerican Symposium on Landslides – ABMS, Vol. 2, pp. 863 – 878. Río de Janeiro.

CONCLUSIONES El modelo y mapas de amenaza para el pronóstico de deslizamientos permite identificar espacial y temporalmente en forma temprana la amenaza por deslizamientos y puede ser utilizado como herramienta en el proceso de toma de decisiones para disminuir el riesgo por estos fenómenos. El pronóstico temprano de la amenaza por deslizamientos permitirá orientar acciones de prevención y mitigación de los efectos que puedan ocasionar los deslizamientos. El modelo permitirá avanzar en el proceso de implementación de la cultura de prevención de desastres que está implementando la nación y el mundo en general. El modelo deberá ser usado como una herramienta de pronóstico de orden nacional y para mayor resolución y precisión del pronóstico se debe continuar trabajando con información fuente de mayor resolución. Agradecimientos Los autores agradecen al Dr. Carlos Castaño Uribe, Director General y a todo el personal técnicocientífico del Instituto de Hidrología, Meteorología y Estudios Ambientales, IDEAM que colaboraron en el desarrollo de este trabajo. A las directivas nuestro agradecimiento por apoyar la participación en el evento y autorizar la publicación del presente artículo.

COATES, D. (1990): Geología Ambiental. 27p. Ingeominas, 87 p., Bucaramanga. IDEAM (1996) Mapa de coberturas vegetales uso y ocupación del territorio. Escala 1:500.000. Bogotá, Colombia. IDEAM UNAL. (1996). Sistemas morfogénicos del territorio Colombiano. Mapa a escala 1:500.000. Bogotá, Colombia. IDEAM, 1998. El Medio Ambiente en Colombia. Santa Fe de Bogotá, Colombia, Julio de 1998. IDEAM, 2001. Estudio Ambiental de la Cuenca Magdalena Cauca y Elementos para su ordenamiento Territorial. Bogotá, Colombia. IDEAM, 2001. Los Fenómenos Cálido y Frío del Pacifico (El Niño - La Niña) y su incidencia en la estabilidad de laderas en Colombia. Bogotá, Colombia. INGEOMINAS (1988). Mapa geológico de Colombia. Escala 1:1.500.000 THE JAPANES GEOTECHNICAL SOCIETY (1997). Manual for zonation on areas susceptible to rain _ induced slope failure. Asian technical committee of geotechnology for natural hazard in ISSMFE, Tokyo. VARGAS CUERVO G., (1999). Guía técnica para zonificación de la susceptibilidad por movimientos en masa. GTZ PRG, Villavicencio, Meta. 197 p.