Caracterización del área

El municipio de Dabeiba se localiza en el occidente del departamento de Antioquia sobre la vertiente del río Atrato, en la parte norte de la Cordillera Occidental de los Andes colombianos (Figura 1 ). El área afectada tiene una extensión aproximada de 22,6 km² y una altura promedio de 1320 m s. n. m., y se ubica al sur de la cabecera municipal del municipio de Dabeiba, en las veredas El Botón, Mohán, Dabeiba Viejo, El Espinazo y El Pital.

Figura 1 Delimitación de las cuencas hidrográficas de las quebradas El Mohán y El Botón, indicadas con contornos rojos, y la red de drenaje representada por líneas azules. El recuadro insertado en la parte superior izquierda indica la ubicación general del área de estudio dentro del contexto nacional de Colombia. 

Geológicamente, el área está conformada por rocas ígneas volcánicas, rocas sedimentarias detríticas y de origen químico de ambiente marino, así como depósitos cuaternarios (Figura 2). Las principales unidades geológicas corresponden al Complejo Cañasgordas de edad cretácica superior, donde se encuentran el Miembro Urrao y el Miembro Nutibara. El Miembro Urrao está compuesto por bloques tectónicos de diabasas y basaltos, junto con paquetes de chert y calizas. Por su parte, el Miembro Nutibara se caracteriza por la presencia de bloques de lodolitas, limolitas y arenitas líticas (^{Gómez et al., 2010}; ^{SGC, 2013}). Los suelos predominantes en el área se clasifican como limos inorgánicos (ML), arenas limosas (SM) y limos de alta plasticidad (MH), según el Sistema Unificado de Clasificación de Suelos (SUCS) (^{Corpourabá, 2018}).

Figura 2 Mapa geológico del área de estudio a escala 1:25.000. Las unidades geológicas se representan con distintos colores según la leyenda. 

La Figura 3 presenta el análisis de cobertura del suelo en ambas cuencas, con patrones distintivos. En la cuenca El Mohán, el bosque predomina, cubre el 61% del área total; seguido por zonas de pastos que abarcan un 29%; áreas de suelo desnudo que representan un 9 %, y zonas urbanas que comprenden menos del 1%. En contraste, la cuenca El Botón presenta una distribución diferente, con un 40% de cobertura boscosa, 34% de pastos, 24% de suelo desnudo y un 2% de áreas urbanas.

Figura 3 Distribución de los diferentes tipos de cobertura y uso del suelo en la zona de estudio, incluidos bosques, pastos y cultivos, suelos desnudos y áreas antrópicas. 

Caracterización morfométrica de las cuencas

La cuenca El Mohán es de orden 5, con un área de 16,6 km2 y un perímetro aproximado de 25 km. Presenta un factor de forma muy alargada y un relieve de 1,7 km. El tiempo de concentración es de 56 minutos, lo cual indica la posibilidad de ocurrencia de flujos hiperconcentrados (^{Arango et al., 2021}). El coeficiente de torrencialidad de la cuenca es de 27,85, lo que sugiere la posible ocurrencia de flujos de detritos (^{Arango et al., 2021}).

Por otro lado, la cuenca El Botón es de orden 4, con un área de 5,9 krtf y un perímetro de 12,757 km. También presenta un factor de forma alargada y un relieve de 1,5 km. El tiempo de concentración es de 31 minutos, lo cual indica la capacidad de generar diferentes tipos de flujos (^{Arango et al., 2021}).

Ambas cuencas presentan características torrenciales y una geometría alargada, con valores bajos de elongación, circularidad y factor de forma, así como una pendiente media de 12°. Estas características indican un fuerte control estructural en el desarrollo de ambas cuencas. La curva hipsométrica también señala un estado de rejuvenecimiento relativo, con valores de 0,5 para la cuenca El Mohán y 0,4 para El Botón. Una síntesis de los índices obtenidos para ambas cuencas se presenta en la Tabla 1.

Evento hidrometeorológico detonante

Dentro de la zona de estudio no existen estaciones pluviográficas con resolución temporal horaria (o mayor), por esta razón, se caracterizó el evento de precipitación utilizando los datos de una estación pluviográfica del Instituto de Hidrología, Meteorología y Estudios Ambientales (IDEAM), denominada Dabeiba 2 AUT - 11117050, con registros cada 10 minutos. Esta estación se encuentra a 7,5 km del área de estudio. Además, se utilizaron estimaciones satelitales de lluvia del producto IMERG (Integrated Multi-satellite Retrievals for GPM) en su versión V06, que presenta una resolución temporal de media hora (^{Huffman, 2021}). Los datos del IMERG corresponden a estimaciones de la precipitación obtenidas a partir de sensores de la constelación de satélites GPM (Global Precipitation Measurement Missions) y están calibrados con estaciones pluviográficas en tierra. La resolución espacial de estos datos es de aproximadamente 11 km x 11 km para el trópico.

En la Figura 4 se presenta la evolución temporal de la precipitación entre las 00:00 horas del 13 de noviembre de 2020 y las 12:00 horas del 14 de noviembre de 2020. En la Figura 5 se muestra la evolución temporalde la precipitación promedio en la zona de las cuencas de las quebradas El Mohán y El Botón, utilizando el producto IMERG.

Figura 4 Precipitación medida en intervalos de 10 minutos (barras grises) y precipitación acumulada (línea roja) en la estación Dabeiba 2 (AUT-11117050) durante el evento lluvioso entre las 00:00 horas del 13 de noviembre y las 12:00 horas del 14 de noviembre de 2020. 

Figura 5 Precipitación media estimada cada media hora (barras blancas) y precipitación acumulada (línea roja discontinua) con datos IMERG en la región de las cuencas de las quebradas El Mohán y El Botón, entre las 00:00 horas del 13 de noviembre y las 12:00 horas del 14 de noviembre de 2020. 

En las figuras 4 y 5 se observan dos eventos de precipitación antes de las 00:00 del 14 de noviembre. Según la información satelital de IMERG, el primer evento se inició alrededor de las 13:00 horas, con su intensidad máxima en torno a las 16:30 horas, y se redujo posteriormente hasta las 18:30 horas. Los datos de la estación pluviográfica indican que este evento se intensificó rápidamente a partir de las 16:00 horas, y alcanzó su máxima intensidad dentro de la siguiente hora. Luego de este primer evento, se desarrolló un segundo evento de precipitación, considerablemente más fuerte, con un pico de intensidad registrado entre las 19:00 (según IMERG) y las 20:00 horas (según la estación pluviográfica), probablemente asociado a un sistema de convección profunda. Tras el pico de intensidad, la precipitación disminuyó hasta las 23:00 horas. A partir de ese momento, se desarrolló un nuevo evento de precipitación hasta las 4:00 a. m. del día siguiente. Este último evento fue clave y probablemente detonó las avenidas torrenciales en las cuencas de las quebradas El Mohán y El Botón. La máxima intensidad de la precipitación cerca de la medianoche se observa en los datos de la estación pluviográfica, pero no se ve reflejada en los datos de IMERG, probablemente debido a su carácter promedio del área.

Desde el punto de vista hidrológico, no solo fue relevante el evento de lluvia cerca de la medianoche, sino también los eventos de lluvia precedentes, que incrementaron la humedad del suelo, disminuyeron su capacidad de infiltración y aumentaron la escorrentía superficial.

Una ventaja de la información satelital, como la de IMERG, es su amplia cobertura espacial, lo que permite examinar los eventos de precipitación en un contexto regional y nacional. En la Figura 6A y Figura 6B se presentan mapas de la precipitación acumulada entre las 12:00 del 13 de noviembre y las 4:00 del 14 de noviembre de 2020, para todo Colombia y un enfoque en el noroccidente del país. Se observa que las secciones norte de la Cordillera Occidental y de la costa Pacífica presentaron acumulaciones significativas de precipitación en comparación con el resto del país. Ninguna otra región de Colombia tuvo acumulados de lluvia similares a los registrados en la Cordillera Occidental. La avenida torrencial del 13-14 de noviembre de 2020 en las cuencas de las quebradas El Mohán y El Botón se enmarca en un contexto meteorológico favorable para la generación de acumulados de precipitación significativos en la Cordillera Occidental.

Figura 6 Distribución espacial de la precipitación acumulada obtenida mediante datos satelitales IMERG. A. Totalidad del territorio colombiano. B. Enfoque ampliado de la región del noroccidente de Colombia. La caja roja destaca el área de estudio que incluye las cuencas de las quebradas El Mohán y El Botón. Los colores representan los valores de precipitación acumulada; varían desde áreas con bajas precipitaciones (naranja) hasta valores máximos superiores a 120 mm (azul oscuro). 

Durante el día del evento, se observó abundante nubosidad en el país, especialmente en la región Andina, la región Pacífica y en la costa del Pacífico y Caribe. Un aspecto relevante fue la gran cantidad de nubosidad asociada a la depresión tropical Iota, localizada al norte de Colombia en ese momento. La proximidad de la depresión tropical Iota generó cambios en los patrones de circulación atmosférica, e impactó el desarrollo de los sistemas de precipitación en el norte del país (^{Meza-Núñez, 2022}). Este ciclón tropical intensificó el Chorro del Chocó y la circulación de bajo nivel desde el Pacífico en dirección oeste-este, y favoreció la advección de humedad hacia la Cordillera Occidental. En la tropósfera media (700 hPa), también se observó una perturbación en el campo de circulación debido a la depresión tropical Iota, lo cual facilitó el ingreso de humedad a niveles medios hacia las cordilleras.

La Figura 7, basada en el reanálisis ERA5 del European Centre for Medium-Range Weather Forecasts (ECMWF) (^{Hersbach et al, 2023}), muestra la humedad relativa promedio y el campo de viento a 500 hPa. Se observan valores altos de humedad relativa sobre la cordillera, indicadores de un ambiente atmosférico favorable para el desarrollo de nubes con grandes cantidades de agua precipitable, lo cual se refleja en las lluvias abundantes que ocurrieron ese día y que detonaron las avenidas torrenciales en Dabeiba. Además, los vectores de viento indican el transporte de humedad en niveles medio-altos desde el Pacífico hacia el noroccidente de Colombia, facilitado por la depresión tropical Iota.

Figura 7 Distribución espacial de la humedad relativa (en %) y el campo de viento a 500 hPa sobre Colombia el 13 de noviembre de 2020, según datos del modelo ERA5. Las áreas en tonos azulados representan zonas con alta humedad relativa, mientras que los tonos anaranjados indican menores porcentajes de humedad. Los vectores negros ilustran la dirección e intensidad del viento a nivel de 500 hPa. 

En regiones tropicales, el desarrollo de estos sistemas de precipitación se ve favorecido por la alta disponibilidad de humedad y el calentamiento atmosférico rápido debido a la radiación solar. En una zona de orografía compleja como el extremo norte de la Cordillera de los Andes, donde se encuentra el municipio de Dabeiba, la convección se intensifica debido al ascenso orográfico de masas de aire húmedo propiciado por la topografía montañosa. Cuando las condiciones atmosféricas son favorables, se generan sistemas convectivos a mesoescala, con una extensión de cientos de kilómetros, que son capaces de detonar eventos de alta intensidad, como fue el caso de Salgar en 2015 (^{Hoyos et al, 2019}).

Descripción del enjambre de movimientos en masa y del evento torrencial

Los movimientos en masa que se presentaron en las cuencas El Mohán y El Botón se describen como estrechos y de poca profundidad respecto a la longitud de la superficie de falla y la propagación de dichos eventos. La profundidad de la superficie de falla varía entre 1 y 3 metros, lo cual corresponde al contacto entre el suelo residual y el saprolito, y presenta una geometría plana y paralela a la ladera.

La ubicación de la superficie de falla, en el contacto entre el horizonte del suelo residual y el horizonte saprolítico, señala la formación de un nivel freático colgado. Este nivel freático se genera por el cambio de permeabilidad y conductividad hidráulica entre ambos horizontes de suelo. El suelo residual está influenciado generalmente por raíces de la vegetación, la intervención de animales y procesos químicos que generan macroestructuras y conductos preferenciales, que incrementan en varios órdenes de magnitud la conductividad hidráulica. Debido a esto, durante el proceso de infiltración del agua en el suelo, se genera un flujo subsuperficial con altura Zw. Al aumentar Zw, también aumentan las presiones de poros positivas, lo cual reduce los esfuerzos efectivos y, por lo tanto, disminuye la resistencia al cortante.

La falla de la ladera bajo presiones de poros positivas se confirma por el comportamiento tipo flujo del material fallado, generado por un proceso denominado licuefacción estática, en el cual el material del suelo, bajo condiciones de saturación o cercanas a la saturación, falla rápidamente, causa una pérdida de resistencia al cortante y deforma completamente el material.

Del evento ocurrido entre los días 13 y 14 de noviembre de 2020 en las cuencas de las quebradas El Mohán y El Botón, se cartografiaron un total de 175 movimientos en masa, como se observa en la Figura 8. Estos movimientos fueron inventariados a través de sensores remotos, como imágenes satelitales de Google Earth con resolución submétrica posterior a la fecha del evento.

Figura 8 Distribución espacial de los deslizamientos desencadenados por el evento de precipitación en la región de Dabeiba, Antioquia. Los elementos del inventario incluyen las coronas y cuerpos de los movimientos en masa. 

Las unidades geológicas donde se desencadenaron los procesos son basaltos El Botón (57%), sedimentitas de Dabeibaviejo (26%) y el Complejo Cañasgordas -Miembro Urrao (17%). Estas unidades se caracterizan por estar compuestas por suelos residuales de baja cohesión. En términos de pendiente, las laderas afectadas presentan inclinaciones en un rango entre 28° y 43°, con direcciones predominantes hacia el este y el oeste.

Como se observa en la Figura 8, los movimientos en masa se localizaron principalmente en la parte media de la cuenca, lo cual puede estar relacionado con la distribución espacial reducida del evento de lluvia o las mayores intensidades del evento. El área afectada por los procesos varía desde 0,06 hasta 4,33 hectáreas, con una sumatoria total de 75,4 hectáreas. Las longitudes de los movimientos oscilan entre los 30 y los 700 metros. De los procesos registrados, el 58% corresponde a deslizamientos traslacionales (algunos con un componente de flujo), mientras que el 42% corresponde a flujos de lodos y detritos. Espacialmente, el inventario se encuentra distribuido principalmente en las partes centrales de ambas cuencas. Aproximadamente el 54% de los movimientos en masa ocurrieron en la cuenca El Botón, donde predominan los deslizamientos traslacionales con cerca del 62%. En la cuenca El Mohán se registró poco más del 45% de los movimientos, con predominancia de deslizamientos traslacionales, que representan el 54% del total.

Descripción del evento torrencial

No existen registros visuales del flujo torrencial que permitan describir directamente el evento. Por lo tanto, se realizó un trabajo de campo detallado para identificar las evidencias asociadas a fenómenos torrenciales dentro de la cuenca. Adicionalmente, se modeló el evento torrencial mediante el programa Transient Rainfall Infiltration and Grid-based Regional Slope-Stability analysis (TRIGRS) (^{Baum et al, 2008}), lo cual permitió comparar los resultados con las evidencias obtenidas en campo.

En las cuencas El Mohán y El Botón se registraron 44 indicios de la ocurrencia de avenidas torrenciales a lo largo del canal, incluidos 9 sitios con cobertura de arena limosa sobre llanuras o terrazas, 12 puntos con marcas de erosión en corteza de árboles, 18 sitios con remanentes de represión y 5 sitios con sedimentos adosados en paredes del canal o en cortezas de árboles. En la Figura 9A se presenta una superficie de arena limosa sobre una llanura aluviotorrencial, donde se observa material fino depositado sobre los primeros niveles como llanuras o terrazas bajas, que indican avenidas muy recientes. La Figura 9B muestra remanentes de represión, una de las evidencias más frecuentes en estas cuencas. Estos remanentes son áreas donde persisten materiales que obstruyeron o continúan obstruyendo el flujo de la cuenca, típicamente compuestos por grandes bloques de roca y troncos de árboles. Las marcas de erosión en la vegetación son comunes debido a la gran cantidad de bosque en el lecho de la quebrada, donde el paso de un evento torrencial deja evidencia del nivel alcanzado por el flujo (Figura 10A). Por último, otra evidencia de torrencialidad son los sedimentos adosados contra la vegetación o las paredes del canal, ilustrados en la Figura 10B, donde se observa un canto incrustado contra el tronco de un árbol.

Figura 9 A. Superficie areno-limosa sobre llanura fluviotorrencial. B. Remanentes de represamientos en la quebrada El Mohán. 

Figura 10.A Marcas de erosión en la vegetación. B. Sedimentos adosados en la vegetación. 

Un elemento fundamental del flujo torrencial es el tiempo de concentración de cada cuenca. El tiempo de concentración representa el tiempo que tarda una gota de agua en recorrer desde el punto más alejado de la cuenca hasta el punto de salida (^{Smith y Vélez-Otálvaro, 1997}). Este tiempo se da en función de las características morfométricas de cada cuenca, presentadas previamente. Numerosos autores han desarrollado formulaciones empíricas para este cálculo. En este trabajo se calcularon los tiempos de concentración bajo diferentes formulaciones empíricas, incluidas Kirpich, Williams, Johstone y Témez, tal como se presenta en ^{Vélez-Upegui y Botero-Gutiérrez (2011)}. Para la cuenca El Mohán se obtuvieron tiempos de concentración de entre 37 y 136 minutos, con un valor medio de 55 minutos. Para la cuenca El Botón se encontraron valores de entre 16 y 70 minutos, con un valor medio de 31 minutos.

Para la estimación de caudales máximos se calcularon las curvas Intensidad-Duración-Frecuencia (IDF) y se utilizó el método de hidrógrafa unitaria del Servicio de Conservación de Suelos de Estados Unidos (Soil Conservation Service - SCS). En la cuenca El Mohán, se estimaron precipitaciones de 74 mm/hora, con caudales de hasta 100 m³/s; mientras que en la cuenca El Botón se estimaron precipitaciones de 52 mm/hora, con caudales de 65 m³/s, ambos para periodos de retorno de 100 años. En el escenario más crítico evaluado, con un periodo de retorno de 500 años, se estimaron caudales máximos de 146 m³/s para la quebrada El Mohán y de 73 m³/s para la quebrada El Botón.

En cuanto a los tiempos de tránsito de la creciente, en la cuenca El Mohán, el tiempo total es de aproximadamente 240 minutos (4 horas), con el pico de la creciente en unos 80 minutos desde el inicio de la precipitación efectiva. La fase de recesión tras el pico máximo de caudal tarda alrededor de 160 minutos (2 horas y media). Para la quebrada El Botón, el tiempo total de tránsito de la creciente es de aproximadamente 130 minutos (2 horas y 10 minutos), con el pico de la creciente alrededor de los 40 minutos desde el inicio de la precipitación efectiva, y una recesión de unos 90 minutos (una hora y media). Estos valores incluyen únicamente los aportes de escorrentía superficial derivados de la precipitación efectiva (caudal líquido), sin considerar los aportes sólidos generados en la cuenca.

Para la estimación del volumen de sedimentos transportados en el evento se utilizó el modelo TRIGRS (^{Baum et al, 2008}), con el objetivo de estimar el volumen de sedimentos generado por movimientos en masa detonados por lluvias. Además, se empleó un modelo de erosión que utiliza la ecuación universal para la pérdida de suelo. El modelo TRIGRS, escrito en Fortran, evalúa la distribución temporal y espacial de movimientos en masa superficiales detonados por lluvia, calcula los cambios transitorios de la presión de poros y su impacto en la variación del factor de seguridad. Los parámetros usados para la implementación del modelo fueron las condiciones topográficas de la zona, el mapa de espesor del suelo, las propiedades mecánicas e hidrológicas del suelo, las condiciones de lluvia como factor detonante y las condiciones iniciales del agua en el suelo.

El mapa de espesor del suelo se obtiene a partir del modelo matemático propuesto por ^{Catani et al. (2010)}, donde el espesor se estima relacionando el ángulo de la pendiente. Los parámetros geotécnicos se obtuvieron a partir de la parametrización definida para las zonas geotécnicas homogéneas (Tabla 2). Con respecto a los parámetros hidráulicos, se considera el proceso de infiltración bajo condiciones saturadas, por ende, solo se emplea la conductividad hidráulica saturada (Ks), la cual se estima mediante aproximaciones dadas por funciones de pedotransferencia (PTF) con base en la textura del suelo, consignada en la Tabla 2.

En cuanto al modelo de erosión, se utilizó el modelo clásico de RUSLE (por sus siglas en inglés, Revised Universal Soil Loss Equation) de acuerdo con la implementación recientemente propuesta por ^{Tefera et al.} (manuscrito en revisión), que considera 6 factores: factor de erosividad de la lluvia (R), factor de erodabilidad del suelo (K), factor de longitud y pendiente de la ladera (LS), factor de cobertura (C) y factor de prácticas de conservación (P). Los resultados del modelo muestran que para la quebrada El Mohán, los caudales sólidos estimados fueron de 130 m3/s tanto para periodos de retorno de 100 como de 500 años, sin variaciones significativas entre ambos escenarios. En el caso de la quebrada El Botón, se estimaron caudales sólidos de 71 m3/s para un periodo de retorno de 100 años, mientras que, para un periodo de retorno de 500 años, el caudal sólido aumentó a 110 m3/s, lo cual evidencia una mayor sensibilidad a eventos extremos en esta cuenca.

Los resultados de ambos modelos fueron integrados para obtener los volúmenes totales de sedimentos transportados durante el evento. La validación de los resultados se realizó mediante la comparación del inventario de 175 movimientos en masa documentados con las zonas de falla identificadas por el modelo. Los resultados muestran una correspondencia espacial significativa, particularmente en las áreas donde se concentraron los mayores daños durante el evento del 13-14 de noviembre de 2020. Adicionalmente, se contrastaron las alturas de flujo modeladas con las evidencias de campo documentadas, incluidas las 44 evidencias de avenidas torrenciales registradas, que comprenden superficies areno-limosas sobre llanuras o terrazas, marcas de erosión en cortezas de árboles, remanentes de represamiento y sedimentos adosados en paredes de canal.

La simulación hidrodinámica se llevó a cabo mediante la herramienta r.avaflow desarrollada por ^{Mergili et al. (2017)}, una herramienta computacional para la simulación de flujos, avalanchas y otros procesos concatenados, que emplea un modelo bifásico (sólido y fluido) o de única fase con fricción tipo Voellmy.

La funcionalidad central de r.avaflow consiste en la redistribución de la masa y el momento empleando un modelo de flujo dinámico y un esquema numérico a través de un sistema de seis ecuaciones diferenciales para un modelo bifásico. Este modelo emplea la plasticidad de Mohr-Coulomb para el esfuerzo del sólido. El esfuerzo del fluido es modelado como un esfuerzo viscoso no newtoniano de un volumen sólido con gradiente de fracción mejorado. La transferencia de impulso interfacial generalizada incluye resistencia viscosa, flotabilidad y masa virtual inducida por la aceleración relativa entre las fases.

Para la simulación de una avenida torrencial en la zona de estudio se evaluaron los periodos de retorno de 2; 2,33; 5; 10; 50; 100; 300 y 500 años a través de las hidrógrafas del flujo para el caudal líquido y los volúmenes de sedimentos estimados de movimientos en masa y erosión de laderas. Estos sedimentos son distribuidos en el tiempo a través de la hidrógrafa de creciente de cada periodo de retorno.

En las figuras 11A y 11B se presentan los resultados obtenidos para las cuencas El Botón y El Mohán bajo un escenario de periodo de retorno de 500 años. En ambas cuencas, los resultados muestran que la mayoría de las acumulaciones significativas (superiores a 1 m) se encuentran concentradas en el cauce principal, lo cual indica la trayectoria preferencial del flujo hacia la parte baja del valle. Se puede observar que, aunque las áreas inundadas sean similares en extensión, la cuenca de la quebrada El Botón presenta áreas donde el flujo excede los 3 m con mayor frecuencia en comparación con la quebrada El Mohán. Lo que señala una mayor disponibilidad de sedimentos o mayores pendientes que facilitan el movimiento de flujo torrencial en la cuenca El Botón.

Comparado con la Figura 12, los resultados de modelación se ajustan a las características cartografiadas del evento del 13-14 de noviembre de 2020, el cual estuvo limitado al valle estrecho de la cuenca, con pequeñas áreas donde el valle se ensancha levemente, y donde se registraron espesores de flujo superiores a 3 m.

Figura 11 Resultados de la simulación de los flujos para las cuencas de las quebradas A. El Botón y B. El Mohán bajo un escenario de periodo de retorno de 500 años. Los colores indican la altura del flujo (en metros), desde menos de 0,1 m (en blanco) hasta alturas superiores a 3 m (en azul). 

Figura 12 Movimientos en masa, zonas de depósito y zonas de tránsito del material movilizado en el evento de avenida torrencial ocurrido en las cuencas de las quebradas El Mohán y El Botón durante el 13 y 14 de noviembre de 2020. Las zonas de depósito y tránsito están resaltadas en rosa y amarillo, respectivamente. 

La simulación fluidodinámica reveló que los flujos de escombros generados por estos movimientos en masa se concentraron sobre todo en los cauces principales, y alcanzaron alturas de flujo de hasta 3 m en algunos tramos. Sin embargo, es importante reconocer las limitaciones de estos modelos. La principal radica en la homogeneidad asumida en las propiedades del suelo dentro de cada unidad geológica, cuando en realidad existe una variabilidad espacial significativa. Otras limitaciones incluyen la incertidumbre en la estimación del espesor del suelo, la resolución temporal de los datos de lluvia disponibles, la simplificación del proceso de infiltración y la no consideración de efectos tridimensionales en la estabilidad.

La modelación integrada de los procesos de movimientos en masa y flujo torrencial mediante los modelos TRIGRS y r.avaflow ha permitido una comprensión más completa de la dinámica del evento del 13-14 de noviembre de 2020 en las cuencas El Mohán y El Botón. Los resultados del modelo TRIGRS muestran que las zonas de mayor susceptibilidad a movimientos en masa coinciden con las áreas de mayor pendiente y menor cobertura vegetal, particularmente en las partes medias y altas de ambas cuencas. Esto se correlaciona con el inventario de movimientos en masa observados durante el evento, y valida la capacidad predictiva del modelo. Por su parte, la simulación fluidodinámica con r.avaflow reveló que los flujos de escombros se concentraron sobre todo en los cauces principales, y alcanzaron alturas de flujo de hasta 3 m en algunos tramos. La interacción entre estos dos procesos es evidente en la distribución espacial de los depósitos fluviotorrenciales recientes, que se concentran en las zonas de transición entre las laderas y los valles.