Artículos

 

Evolución post-lahárica de un canal proglaciar: garganta de Huiloac (México)

 

Post-laharic evolution of a proglacial cannel: the Huiloac creek (Mexico)

 

Nuria Andrés de Pablo^1,*, José Juan Zamorano Orozco², José Juan de Sanjosé Blasco³, Luis Miguel Tanarro García¹, David Palacios Estremera¹

 

¹ Departamento de Análisis Geográfico Regional y Geografía Física. Facultad de Geografía e Historia. Universidad Complutense de Madrid. 28040 Madrid, España.

² Instituto de Geografía, Universidad Nacional Autónoma de México, Coyoacán, 04510 México D.F., México.

³ Departamente de Expresión Gráfica, Escuela Politécnica, Universidad de Extremadura, Cáceres, España. ^*nuriand@ghis.ucm.es.

 

Manuscrito recibido: Septiembre 22, 2013.
Manuscrito corregido recibido: Octubre 17, 2013.
Manuscrito aceptado: Octubre 22, 2013.

 

Resumen

En el último periodo eruptivo del volcán Popocatépetl (19º 02' N, 98º 37' W, 5424 m) se produjeron varios lahares sin-eruptivos que modificaron la red fluvial de su ladera nororiental, donde se localiza el curso del río Huiloac. Los lahares responsables de los cambios morfológicos más importantes fueron los ocurridos en 1997 y 2001. Aunque desde entonces se han producido algunos eventos explosivos, como las recientes erupciones de mayo y junio de 2013, en la garganta de Huiloac no se ha vuelto a registrar ningún lahar relacionado con la actividad volcánica.

El presente artículo propone una metodología para detectar y cuantificar los cambios morfológicos, así como determinar la dinámica de procesos erosivos y sedimentarios en un tramo de la garganta de Huiloac, tras el paso del lahar sin-eruptivo de 2001, durante un primer periodo de siete años. Este método combina la interpretación de una serie temporal de mapas geomorfológicos y de perfiles topográficos transversales del canal, con ayuda de herramientas de CAD (Computer-aided design) y de SIG (Sistemas de Información Geográfica). El análisis de la evolución geomorfológica y topográfica se cruza finalmente con la información meteorológica disponible sobre precipitaciones.

Los resultados obtenidos muestran que la fase inicial (hasta octubre de 2002) estuvo caracterizada por la incisión y evacuación del material que colmató el cauce durante el lahar de 2001. Posteriormente, la acción geomorfológica de las aguas corrientes y la dinámica de laderas ensancharon y profundizaron el canal, aunque también se produjo sedimentación, como lo atestiguan bancos y terrazas fluvio-laháricos. La frecuencia y capacidad de los lahares secundarios, alimentados por las precipitaciones y el agua del deshielo glaciar, fueron los factores que determinaron el dominio de procesos erosivos o de sedimentación en Huiloac. De esta manera, se ha observado que en épocas menos lluviosas y con precipitaciones regulares, el canal experimentó menos variaciones morfológicas y las cantidades erosionadas son también menores, como en el periodo 2004 ‒ 2006. En cambio, las precipitaciones excepcionales en la estación seca (enero de 2002 y enero de 2004) y los aguaceros continuados al final de la estación húmeda (septiembre y octubre de 2007) activaron los procesos de erosión por incisión y por zapa lateral, lo cual produjo el vaciado de los depósitos laháricos que tapizan el canal.

Una vez concluido el primer periodo de observación, se concluye que el cauce no ha alcanzado su estabilidad, ya que sus laderas fueron muy inestables debido a la fuerte pendiente que presentaron, a pesar del considerable vaciado de depósitos que ha experimentado en este lapso temporal.

Palabras clave: cartografía geomorfológica, lahares, SIG, CAD, Popocatépetl, México.

 

Abstract

In the last eruptive period of Popocatépetl volcano (19º 02'N, 98º 37' W, 5424 m) several syn-eruptive lahars modified the fluvial network of the northeastern slope, where the Huiloac River is located. The lahars that produced the most important morphological changes took place in 1997 and 2001. Although some explosive events have taken place since then, such as the recent eruptions of May and June 2013, lahars related to volcanic activity have not been recorded in Huiloac gorge.

This article proposes a methodology to detect and quantify morphological changes and to determine the dynamics of erosion and sedimentary processes in a section of Huiloac gorgef or the first seven years after the syn-eruptive lahar of 2001. This method combines the interpretation of a temporal series of geomorphological maps and topographic profiles across the channel using CAD (Computer-aided design) and GIS (Geographic Information Systems). Finally, the analysis of the geomorphic and topographic evolution is cross-referenced with available information on rainfall.

The results show that the initial phase (up until October 2002) is characterized by incision and removal of the material that filled the river bed during the 2001 lahar. Afterwards, the geomorphological action of running water and the dynamics of slopes widened and deepened the channel, although sedimentation occurred as well, as evidenced by river-banks and lahar terraces. The frequency and capacity of secondary lahars, fed by rainfall and glacier melt water, are the factors that determine the domain of erosion or sedimentation in Huiloac. Thus, it has been observed that in seasons with less rain and regular precipitation, fewer morphological variations occurred and a smaller amount of material was eroded, such as in the period 2004 ‒ 2006. On the other hand, exceptional rainfall in the dry season (January 2002 and January 2004) or continuous rains at the end of the wet season (September and October 2007) triggered processes of erosion by incision, lateral erosion, and removal of the laharic deposits. It can be concluded, at the end of this first period of observation, that despite considerable erosion, the riverbed is still unstable due to steep slopes on the sides of the river.

Keywords: Geomorphological mapping, lahars, GIS, CAD, Popocatépetl, Mexico.

 

1. Introducción

Las erupciones volcánicas generan perturbaciones en el paisaje que afectan al sistema hidrológico circundante y en especial, al aporte de sedimentos. Major (2003) expone una lista de las principales consecuencias hidrológicas e hidráulicas que conlleva la actividad volcánica, resaltando los cambios que sufren los almacenes hídricos (en especial la nieve y la cubierta vegetal) y las alteraciones en las características, magnitud y duración de la escorrentía. Concluye que, en general, los cambios hidrológicos incrementan la superficie drenada y reducen el tiempo que tardan las aguas en escurrir por la ladera hasta los canales, de forma que la escorrentía post-eruptiva alcanza los canales más rápidamente y en mayor cantidad que la escorrentía pre-eruptiva.

Gran y Montgomery (2005) hacen especial hincapié en que la principal modificación volcánica consiste en el aporte de una mayor cantidad de sedimentos disponibles en las cabeceras de las redes de drenaje de los volcanes, que resultan muy superiores a las cantidades localizadas en sistemas fluviales no volcánicos. La existencia en las cabeceras fluviales de una capa de materiales fácilmente erosionables y la destrucción de la vegetación hacen que estas superficies experimenten una respuesta rápida ante aportaciones copiosas de agua (Segerstrom, 1950; Waldron, 1967; Kadomura et al., 1983; Janda et al., 1984; Yamamoto, 1984; Leavesley et al., 1989; Shimokawa et al., 1989).

En la mayoría de los volcanes el transporte extraordinario de sedimentos se realiza mediante flujos rápidos que se generan por lluvias intensas, fusión repentina de nieve y hielo, o liberación de aguas embalsadas. Estos procesos reciben el nombre de lahares (Smith y Fritz, 1989; Rodolfo, 1989). Se puede distinguir entre lahares primarios o sin-eruptivos y secundarios o post-eruptivos, sin relación con una erupción concreta (Lavigne, 1998; Vallance, 2000; Thouret y Lavigne, 2000). Los lahares también pueden presentar distintos tipos de flujo (avalancha de escombros, también denominado debris flow, flujo hiperconcentrado y flujo de arroyada) en relación con la proporción de agua y materiales sólidos que contengan (Smith y Lowe, 1991).

La gran mayoría de las investigaciones realizadas sobre flujos hidrovolcánicos centran su interés en el lapso de tiempo inmediatamente posterior a la erupción, cuando se produce la mayor perturbación en la red de drenaje y se desarrollan los lahares primarios. Existe una gran profusión de publicaciones que analizan las características de estos flujos en distintos volcanes como en los montes Usu, Unzen, y Sakurajima en Japón (Kadomura et al., 1983; Shimokawa y Taniguchi, 1983; Chinen y Kadomura, 1986; Mizuyama y Kobashi, 1996); en el Galunggung de Indonesia (Hamidi, 1989; Hirao y Yoshida, 1989); en el Paricutín y el Chichón de México (Inbar et al., 1994, Inbar et al., 2001); en el Ruapehu de Nueva Zelanda (Cronin et al., 1999; Hodgson y Manville, 1999); en el Mayon y el Pinatubo de Filipinas (Rodolfo, 1989; Rodolfo y Arguden, 1991; Pierson et al., 1992; Pierson et al., 1996; Daag, 1994; Daag, 2003; Major et al., 1996; Scott et al., 1996; Newhall y Punongbayan, 1996; Umbal, 1997; Tuñgol, 2002).

Otro grupo de investigaciones concentra la atención sobre las precipitaciones que desencadenan los lahares, especialmente en los parámetros de cantidad, intensidad y recurrencia (Rodolfo y Arguden, 1991; Iwamoto, 1996; Tuñgol y Regalado, 1997; Lavigne, 1998; Thouret y Lavigne, 2000; Lavigne y Thouret, 2002; Daag, 2003; Scott et al., 2005; van Westen y Daag, 2005; Barclay et al., 2007; entre otros estudios).

Los datos obtenidos en todos estos trabajos citados se han empleado como valores de entrada para realizar simulaciones de lahares suponiendo diferentes escenarios, a partir de las cuales se valoran los riesgos asociados y se determinan las medidas para reducir sus efectos (Daag y van Westen, 1996; van Westen, 1997; Renschler, 2005; Ackerman et al., 2000; Long, 2000; Muñoz-Salinas et al., 2007, 2009; Haddad et al., 2011). Pero si la precisión en los valores de entrada referentes a la cantidad de agua y material del flujo, su viscosidad o velocidad son condiciones imprescindibles para que el modelo o la simulación resulten aceptables, también lo es el empleo de una base topográfica detallada de los cauces. De hecho, los manuales de hidráulica recogen como principio básico que todo flujo está influido en gran medida por las características morfológicas del canal por el que discurre (p.e. Chow, 1959). De esta manera, la profundidad, la anchura, la pendiente longitudinal, la sinuosidad y la rugosidad (materiales que revisten lecho y paredes) son los principales parámetros del cauce que condicionan variables como la velocidad del flujo o la capacidad o competencia de la carga.

Por otra parte, los escasos trabajos dedicados a evaluar las modificaciones sufridas por la red de drenaje después de un lahar de grandes dimensiones (Janda et al., 1981; Pierson, 1985; Pierson y Scott, 1985; Rodolfo, 1989; van Westen, 1997; Garcin et al., 2005; Gran y Montgomery, 2005; Barclay et al., 2007), concluyen que el impacto producido en las características hidrológicas y geomorfológicas de la red de drenaje en los volcanes se mantiene incluso varias décadas después del evento eruptivo. En este aspecto, son relevantes las conclusiones de los estudios llevados a cabo en la red de drenaje del volcán St. Helens (Cascades Range, EEUU), la cual sufre complejos ciclos de incisión, sedimentación y ensanchamiento en una tendencia claramente no lineal (Meyer y Martinson, 1989; Simon, 1999), de forma que la inestabilidad generada por los lahares primarios en el sistema fluvial (aporte de material y cambios en la topografía) se mantiene a largo plazo y la migración efectiva de mayor cantidad de sedimentos se mantiene funcional durante décadas (Major et al., 2000; Major, 2003; Major, 2004). Atendiendo a estas condiciones, para una adecuada modelización de los flujos y para conseguir cálculos correctos de los mismos, con el fin de mejorar el sistema de alerta a la población, es necesario un seguimiento detallado de los parámetros morfológicos del canal, con especial atención a las variaciones topográficas (Muñoz-Salinas, 2007). Surge entonces, la necesidad de estudios de seguimiento detallado en los cauces de estos sistemas, que pongan de manifiesto cuáles son los procesos geomorfológicos que se suceden en ellos, no solamente como respuesta inmediata a las erupciones, sino durante toda la dinámica de adaptación a la nueva situación (lahares secundarios) que puede durar décadas. Por otra parte, la monitorización detallada de los canales proporcionará además información valiosa para la modelización de futuros lahares, ya que en todo momento ofrece datos topográficos actualizados de una morfología que presenta grandes variaciones en periodos cortos de tiempo.

En esta línea, el objetivo principal de esta investigación se centra en diseñar una metodología válida para llevar a cabo un seguimiento de las modificaciones geomorfológicas experimentadas en cauces, posteriores al paso de un lahar primario de grandes dimensiones. Así mismo, se plantea la forma de detectar los cambios morfológicos del cauce, en términos de erosión y sedimentación, y su relación con las precipitaciones. La metodología propuesta se aplica en detalle a un tramo de la garganta Huiloac, perteneciente a la red de drenaje del volcán Popocatépetl (México). Para agilizar el estudio en el tiempo, ya que este debería prolongarse mientras duren las modificaciones y hasta que el canal muestre señales claras de estabilidad, se han planificado periodos cortos consecutivos de observación de siete años (inferiores a la década).

 

2. Área de estudio

El volcán Popocatépetl (19º 01'N; 98º 37'W; 5424 m. s.n.m.) se encuentra en el sector central de la Faja Volcánica Transmexicana (FVTM) y es uno de los volcanes más activos de México, ya que se localiza en la parte frontal del arco volcánico, en el sur, hacia donde ha ido migrando la actividad volcánica en los últimos 2 Ma (Macías, 2005). En concreto, el Popocatépetl es el volcán más meridional de la Sierra Nevada, cordillera que discurre de N a S y constituye el límite oriental de la Cuenca de México. Se puede considerar también como uno de los más peligrosos, debido a que se sitúa a 65 km al sudeste de México D.F. —con 8.9 millones de habitantes en 2010 (INEGI, 2011a)—, y a 45 km de la ciudad de Puebla —1.5 millones de habitantes en 2010 (INEGI, 2011b)—.

Durante el Holoceno, los tres periodos eruptivos principales del Popocatépetl (3195 – 2830 a.C., 800 – 215 a.C., y 675 – 1095 d.C.) estuvieron caracterizados por un mismo esquema de actividad, que comenzaba con emisiones de pequeñas cantidades de ceniza (flujos y caída) y alcanzaban su máxima intensidad con la deposición de grandes cantidades de pumitas, flujos de cenizas calientes y finalmente se producían extensos lahares (Siebe et al., 1996). Macías (2005) indica que el colapso de las columnas provocó el bloqueo de la red hidrográfica del Popocatépetl. Con la adición de agua, procedente de la fusión de hielo y nieve (González et al., 1997), o de lluvias torrenciales (Macías, 2005), se generaron lahares que esparcieron la carga sedimentaria por la cuenca de Puebla. Este mecanismo fue probablemente el que generara el lahar denominado San Nicolás (González et al., 1997; González, 2000), que se produjo en la garganta Huiloac hace 1100 – 1300 años. Este flujo, con un volumen de 5 x 107 m³, se inició a media ladera del volcán y viajó unos 60 km hasta alcanzar la ciudad de Izúcar de Matamoros.

No se tiene referencia de lahares producidos entre 1095 y el inicio del actual periodo eruptivo en 1994 (Cruz-Reyna et al., 1995), aunque parece lógico pensar lo contrario (Muñoz-Salinas, 2007), ya que cualquier emisión de material supone una perturbación en el sistema fluvial circundante, que se manifiesta con un aumento del caudal y de la carga sedimentaria en los ríos.

En el presente periodo eruptivo iniciado en diciembre de 1994, los principales eventos laháricos en Huiloac han tenido lugar en 1995, 1997 y 2001. En abril de 1995 una emisión cubrió de piroclastos aproximadamente el 50 % del glaciar, lo que ocasionó una fusión repentina. El agua de fusión se canalizó por las gargantas proglaciares (Tepeteloncocone, La Espinera y Tenenepanco en la Figura 2.A) y saturó los materiales que se encontraban en ellas, una mezcla de sedimentos sin consolidar de anteriores avalanchas y ceniza volcánica (Palacios, 1995; Palacios, 1996). La masa se inestabilizó y resbaló sobre la capa congelada, que actuó como un estrato impermeable, hasta los 4020 m.s.n.m. A partir de aquí, el agua se empezó a infiltrar en los depósitos y el flujo se detuvo (Palacios et al., 1998).

Dos años después, entre el 30 de junio y el 1 de julio se produce la erupción mayor registrada en el periodo eruptivo reciente. El Centro Nacional de Prevención de Desastres de México (CENAPRED, 2001) describe que el material incandescente al caer sobre el glaciar causó abrasiones y acanaladuras. El agua de fusión, junto con el procedente de las fuertes precipitaciones que se produjeron en las laderas del volcán (Sheridan et al., 2001), movilizó los sedimentos que tapizaban el fondo de las tres gargantas proglaciares y se generó un flujo hiperconcentrado. Al ir incorporando material, se transformó en un flujo de derrubios (debris flow), para después progresivamente cambiar a un flujo hiperconcentrado en su parte distal, a la vez que depositaba su carga (Capra et al., 2004). Este lahar transportó 1.85 x 105 m³ de material sólido y agua (Muñoz-Salinas et al., 2009), alcanzó una velocidad entre 1.4 y 7.7 m/s (Muñoz-Salinas et al., 2007) y tras 21 km de recorrido, se detuvo justo a las puertas de la población de San Nicolás de los Ranchos.

El lahar de 2001 tuvo un origen diferente. El 22 de enero de 2001 el volcán emitió un flujo piroclástico que recorrió 6 km por la ladera norte del volcán. Al atravesar el glaciar, causó una fuerte abrasión, lo que provocó el deshielo parcial del mismo. El agua de fusión saturó los materiales e inició un lahar en la garganta Tenenepanco, que se continuó por Huiloac y finalizó 2 km antes de alcanzar la población de Santiago Xalitzintla. En los 14.3 km de recorrido, este lahar se comportó como un debris flow (Capra et al., 2004). El volumen lahárico transportado fue de 1.6 x 105 m³ (Muñoz-Salinas et al., 2009), y desarrolló velocidades comprendidas 1.3 y 13.8 m/s (Muñoz-Salinas et al., 2007).

A pesar de que el volcán ha registrado posteriormente eventos eruptivos destacados, como los de la primavera y el verano de 2013, que llevaron a las autoridades a elevar el nivel de alerta (CENAPRED, 2013), la garganta de Huiloac no ha vuelto a ser afectada por ningún lahar directamente relacionado con la actividad volcánica (Andrés et al., 2007).

La red hidrográfica del Popocatépetl presenta una típica distribución radial, fuertemente condicionada por la topografía cónica del volcán (Figura 1.B). La elevada permeabilidad e incoherencia de sus materiales y las fuertes pendientes favorecen la incisión de las aguas corrientes, que forman estrechas y profundas gargantas en las cotas superiores, para luego ensancharse y hacerse menos profundas en áreas más bajas. La ladera norte del volcán está drenada por tres gargantas principales (Tenenepanco, La Espinera y Tepeteloncocone) que captan las aguas de fusión del glaciar y que a partir de su confluencia reciben el nombre de Huiloac. Desde su punto más elevado en el inicio de Tenenepanco (4960 m.s.n.m.) hasta la población de Santiago de Xalitzintla (2560 m.s.n.m.), el cauce Tenenepanco-Huiloac recorre unos 18 km, primero de sur a norte y, tras salvar el obstáculo que supone la pared de la caldera volcánica de Tlamacas (3500 m.s.n.m.), de oeste a este (Figura 2.A). Tras superar la citada localidad, el curso fluvial se une a la garganta Alseseca (Muñoz-Salinas, 2007; Andrés, 2009; Tanarro et al., 2010).

Capra et al. (2004) dividen el eje Tenenepanco-Huiloac en tres sectores determinados por cambios de pendiente. El curso alto presenta pendientes superiores a 11 º y discurre entre los 4960 y los 3400 m.s.n.m., la mayor parte flanqueado por morrenas. En ocasiones el río se encaja en lavas holocenas, lo que indica un proceso de incisión reciente. A partir de los 3400 m.s.n.m. se desarrolla el curso medio, con pendientes entre 11 º y 6 º, encajado entre los paquetes de materiales piroclásticos, fluvioglaciares y laháricos. La profundidad de este tramo alcanza los 40 m y su anchura media del valle es de 20 m. Por último, al entrar en rampa acumulativa volcaniclástica que constituye el piedemonte del volcán (3150 m.s.n.m.), la pendiente se reduce por debajo de 6 º en el curso bajo.

Ante la imposibilidad de recoger los datos necesarios a lo largo de toda la garganta, debido al prolongado tiempo de exposición a los peligros volcánicos y de caída de material, se propone un tramo para realizar un análisis de detalle (Figura 2). Este se localiza en el curso medio de Huiloac, entre los 3200 y 3240 m.s.n.m., donde se pueden observar tanto procesos de erosión como de deposición, mientras que en la cabecera y en la parte distal domina uno de los dos (Chen et al., 2005). Por otra parte, el tramo elegido se encuentra en el bosque monoespecífico de oyamel (Abies religiosa), que se extiende hasta los 3400 m (Beaman, 1962) y ocupa las áreas con precipitaciones más elevadas en las laderas de las montañas del centro de México (García-Romero, 1998), por lo que existe una mayor probabilidad de que se vea afectado por pequeños lahares secundarios alimentados por el agua de precipitación, que en la Sierra Nevada supone aproximadamente 1300 ‒ 1200 mm anuales entre los 3500 y los 3300 m.s.n.m. (Beaman, 1962; Lauer y Stiehl, 1973). Además, es un área accesible a partir de la carretera que une el Paso de Cortés con Santiago de Xalitzintla (Figura 2.A).

 

3. Metodología

Las variaciones de la morfología del canal de la garganta Huiloac se han estudiado para el primer periodo de observación de siete años, desde el lahar sin-eruptivo ocurrido el 22 de enero de 2001, hasta febrero de 2008. La metodología empleada en un primer momento combina el seguimiento de las variaciones topográficas y las transformaciones geomorfológicas del sector seleccionado. Posteriormente, se contrastan estas variaciones con los datos disponibles de precipitaciones, con el fin de relacionar las características de los flujos laháricos secundarios, con los distintos ciclos erosivos o de sedimentación observados, así como con la intensidad de los procesos.

 

3.1. Obtención y tratamiento de la información topográfica

Los datos topográficos se obtuvieron en seis campañas de trabajo de campo que se realizaron en febrero y octubre de 2002; septiembre de 2003; febrero de 2004; marzo de 2006 y febrero de 2008. El primer paso consistió en el trazado de una nivelación siguiendo el recorrido del talweg, donde se marcaron los puntos en los que se localizarían posteriormente los perfiles transversales. La nivelación se realizó con la ayuda de una brújula, mira, cinta métrica y GPS para obtener las coordenadas de los puntos de intersección de la nivelación (perfil longitudinal) con los perfiles transversales. La longitud de la proyección en la horizontal del perfil longitudinal fue de 511.2 m.

En las dos primeras campañas, los perfiles se trazaron con la misma instrumentación que la nivelación (brújula, mira, cinta métrica y GPS), eligiendo un punto en el talweg y fijando los extremos en los levées del canal, marcados con varillas de hierro. Entre estos tres puntos se midieron distancias en la línea del perfil y las secciones se dibujaron en alzados sobre un eje de coordenadas x, y. (Figura 3) en papel y en un programa de CAD (MicroStation). Se tomaron 29 perfiles, denominados P1 a P29, aguas abajo desde el punto de inicio (14Q, 543339.8 m E, 2109354.0 m N: Datum, WGS84), con longitudes que varían entre 10.0 y 17.1 m. Los perfiles se realizaron en los lugares del cauce en los que se observó un cambio en la morfología del cauce, de manera que determinan 28 subtramos (de S1 a S28) con características geomorfológicas similares. Según este criterio, los perfiles no son equidistantes y los subtramos presentan longitudes entre 6.2 y 31.1 m.

A partir de la campaña de septiembre de 2003 los datos topográficos se obtuvieron con una estación total y se estimó un error máximo de +/- 2 cm. Aunque se intentó mantener la localización de todos los perfiles, marcados desde el principio con varillas, no se logró en todos los casos y se modificaron los perfiles 6, 11, 21, 23 y 24. Los datos de la estación se volcaron al programa AutoCad, donde se dibujaron los perfiles en planta y en alzado (Figura 4).

Todos los perfiles se importaron en ArcGis para llevar a cabo comparaciones entre ellos. Los datos de las campañas de febrero y octubre de 2002 se trataron de forma separada, ya que se obtuvieron con técnicas diferentes. El resto se contrastaron dos a dos: entre septiembre de 2003 y marzo de 2006; entre marzo de 2006 y febrero de 2008; entre septiembre de 2003 y febrero de 2008. Las diferencias entre dos perfiles se midieron calculando las áreas de los polígonos que resultaron de la intersección entre los perfiles de dos fechas. Se consideró que un área se correspondía con un valor de erosión, cuando la altura del perfil en el momento t (fecha más antigua) era superior a la que presentaba en el momento t + 1 (fecha más reciente), es decir, sufre una variación negativa entre los dos momentos. La sedimentación en un perfil transversal se definió como el cambio positivo en la altura entre el momento t y el t + 1 (Figura 5). La tabla de atributos de los polígonos intersección se importó a Excel, donde se sumaron las áreas de todos los polígonos resultantes de la intersección que suponían erosión (Aei) y las áreas de los polígonos que representaban deposición (Adi).

Para estimar el volumen erosionado o sedimentado entre dos perfiles se partió del supuesto de que la evacuación y deposición del material a lo largo de cada subtramo es uniforme, ya que la morfología observada presenta escasas variaciones. Se asimiló entonces el volumen del material erosionado o sedimentado al volumen de un cono truncado oblicuo y se emplearon las fórmulas:

Donde:

V_e, V_d: volumen de material erosionado o sedimentado, respectivamente (en m³).

D_i: distancia entre el perfil i y el perfil i + 1 (en m).

A_ei, A_di: Área de la variación negativa o positiva, respectivamente, experimentada entre las dos fechas consideradas en el perfil i (en m²).

A_ei+1, A_di+1: Área de la variación negativa o positiva, respectivamente, experimentada entre las dos fechas consideradas en el perfil i + 1 (en m²).

Como distancia entre perfiles (D_i) se consideró la medida siguiendo el trazado del canal, es decir, teniendo en cuenta la sinuosidad y la pendiente.

 

3.2. Obtención y tratamiento de la información de las variaciones geomorfológicas

El análisis de la evolución geomorfológica de la garganta Huiloac se basa en la comparación de una serie temporal de mapas geomorfológicos. Durante las campañas de campo utilizadas para tomar datos topográficos también se realizaron cartografías geomorfológicas del cauce. En las dos primeras campañas, febrero y octubre de 2002, se hicieron anotaciones, fotografías y dibujos con sus correspondientes descripciones, y se tomaron puntos con GPS para cartografiar posteriormente las formas. En campañas posteriores, ya se contó con una base topográfica de detalle obtenida a partir del levantamiento de 2003, para trabajar directamente en el campo. El resultado fue una colección de seis mapas geomorfológicos elaborados en febrero de 2002, octubre de 2002, septiembre de 2003, febrero de 2004, marzo de 2006 y febrero de 2008.

Se empleó una leyenda geomorfológica sencilla, que atiende al tipo de depósito en el que están labradas las formas y el proceso que las modifica (Figura 6). Así, se distingue una primera categoría (A) constituida por las formas relacionadas con el flujo reciente de lahares primarios. Los márgenes de la garganta están formados por materiales depositados por los lahares de 1997 y 2001 principalmente. La potencia de los depósitos del lahar de 1997 oscila entre los 2.5 y los 4 m y en algunos sectores se encuentran erosionados por flujos posteriores o cubiertos por depósitos del lahar de 2001, que se adaptaron "como una sábana" discontinua a las paredes del cauce. El contacto entre los depósitos laháricos y la antigua pared de la garganta, compuesta por materiales volcánicos y de antiguos lahares, lo constituye un estrecho surco (A1), que en ocasiones se ha rellenado con material erosionado de las paredes superiores, y aparece colonizado por la vegetación. Las crestas de los levées (A2) son resaltes alomados a ambos lados del canal en posiciones topográficas elevadas entre 1 y 4 metros por encima del talweg. En algunos lugares las crestas han sido erosionadas o cubiertas por depósitos de ladera, pero junto con los surcos constituyen los elementos geomorfológicos más estables.

En posiciones topográficas más bajas y enlazando con el lecho se distinguen las laderas del levée, donde se distinguen la ladera superior (A3), sin aparentes modificaciones desde la deposición del material lahárico, y la ladera inferior (A4), con huellas de erosión, en especial, las provocadas por procesos de zapa.

Existen además formas laháricas recientes (B), posteriores al lahar primario de 2001, que constituyen niveles planos, elevados entre 0.20 y 0.70 m por encima del fondo del lecho, a las que se les ha denominado terrazas laháricas (B1). En el salto vertical que presentan hacia el lecho se distingue la composición heterométrica del depósito. En otras ocasiones la forma que presentan estos depósitos laháricos recientes no supone un resalte topográfico superior a 15 cm sobre el lecho. Sus límites no están constituidos por escarpes verticales tan nítidos como en el caso de las terrazas, y su superficie tiene un aspecto ligeramente ondulado, en relación a la cantidad de materiales gruesos o finos que lo componen. Debido a la escasa altura a la que se encuentran sobre el lecho, también son retocados por la acción fluvial, por lo que se les ha denominado bancos fluvio-laháricos (B2).

Las formas anteriores destacan sobre el lecho fluvial propiamente dicho (C1), que presenta un perfil plano o ligeramente cóncavo y está surcado por los estrechos canales de las aguas bajas. El material que se observa en su superficie está compuesto por arenas gruesas y gravas. Se distinguen claramente de los bancos fluvio-laháricos, porque estos presentan perfiles convexos y están formados por depósitos heterométricos.

Mención especial merecen las formas relacionadas con los deslizamientos (D1), que afectan a los depósitos laháricos de 1997 y 2001 y, ocasionalmente, a las laderas constituidas por depósitos volcaniclásticos. Las partes bajas de la garganta acumulan material desplazado de las laderas, tanto por procesos gravitacionales, como por aquellos en los que interviene la fuerza del agua. Este material incluye vegetación arrancada, de manera que no es extraño encontrar troncos de grandes árboles en el fondo de la garganta, que actúan como trampa de sedimentos. Finalmente, se han cartografiado aquellos bloques (E1) que tienen un tamaño del eje mayor superior a 0.5 m. Estos bloques localizados en el lecho, en las superficies de las terrazas y bancos, o en las bases de las laderas son indicadores de la dinámica de la garganta, ya que su estabilidad o modificación indican la fuerza de los procesos que han tenido lugar entre las fechas de estudio.

Las cartografías se digitalizaron en ArcGis (ESRI©) y de manera individual se calculó el área ocupada por cada una de las formas de relieve y el porcentaje que representaba respecto del total. Después se fueron superponiendo por parejas de mapas, de forma que la capa resultado mostrara las áreas que habían sufrido algún cambio y las que permanecían geomorfológicamente estables. Las modificaciones posibles vienen determinadas por las combinaciones que se muestran en la matriz de la Tabla 1. De los cruces posibles, 9 no suponen ningún cambio y se han denominado:

0. Sin cambios en el surco.

1. Sin cambios en las laderas laháricas.

2. Sin cambios en las terrazas laháricas.

3. Sin cambios en los bancos fluvio-laháricos.

4. Sin cambios en el lecho fluvial.

5. Sin cambios en los depósitos de ladera.

6. Bloques estables.

Del resto de posibilidades (72) tan solo se han observado 28 cambios, excluyendo los referentes a los bloques, que se consideran de forma separada, ya que constituyen una superficie muy pequeña respecto del total del tramo de garganta analizado. Las áreas que han experimentado cambios se clasifican atendiendo al proceso inicial que origina las formas.

7. Erosión de los depósitos laháricos primarios: o bien las crestas de los levées se transforman en laderas inferiores (A2-A4), o las laderas superiores pasan a ser laderas inferiores (A3-A4);

8. Sedimentación de depósitos laháricos recientes;

8.1. Sedimentación de terraza con retroceso de laderas laháricas (A2-B1, A3-B1, A4-B1).

8.2. Sedimentación de terraza en áreas ocupadas por bancos o lecho fluvial (B2-B1, C1-B1).

8.3. Sedimentación de terraza en áreas ocupadas por depósitos de ladera (D1-B1).

9. Sedimentación fluvio-lahárica;

9.1. Sedimentación fluvio-lahárica con retroceso de laderas laháricas (A2-B2, A3-B2, A4-B2).

9.2. Sedimentación fluvio-lahárica sobre el lecho fluvial (C1-B2).

9.3. Sedimentación fluvio-lahárica a expensas de depósitos de ladera. (D1-B2).

10. Acción fluvial;

10.1. Retroceso de las laderas laháricas por acción fluvial (A2-C1, A3-C1, A4-C1).

10.2. Desaparición de nuevas formas laháricas por acción fluvial (B1-C1, B2-C1).

10.3. Desaparición de depósitos de ladera por acción fluvial a expensas (D1-C1).

11. Erosión de terrazas laháricas (B1-B2);

12. Erosión de depósitos de ladera (D1-A4);

13. Sedimentación de depósitos de ladera (A2-D1, A3- D1, A4- D1, B1-D1, B2-D1, C1-D1).

Con el fin de facilitar la interpretación de la información se han simplificado en cuatro grupos de cambios. Así, el primer grupo lo constituyen las variaciones constituidas por el retroceso o desaparición de las laderas laháricas primarias (en color naranja en la Tabla 1). En este grupo se han incluido las sustituciones de los depósitos de ladera (D1) por cualquier otra forma, ya que se considera que los D1 están constituidos por materiales de los lahares primarios.

Otro grupo (en color azul en la Tabla 1) lo constituyen los cambios que producen la desaparición o la reducción topográfica de las formas constituidas por nuevos depósitos laháricos o fluvio-laháricos. De igual manera se han considerado los cambios que resultan en una aparición o incremento en la diferencia topográfica con el fondo de las nuevas formas laháricas y fluvio-laháricas (en verde en la Tabla 1). Por último, se considera el cambio ocasionado por la adicción de un depósito de ladera sobre cualquiera de las formas definidas (en color amarillo).

La cuantificación espacial de los cambios se realiza mediante la medida de la superficie que ocupan las modificaciones y su porcentaje de presencia en el total. El resultado se presenta en matrices numéricas.

 

3.1. Obtención y tratamiento de datos sobre precipitaciones

El Servicio Meteorológico Nacional de México ha proporcionado datos de precipitaciones de las estaciones cercanas al área de estudio. Así se cuenta con las series normalizadas de 1971 – 2000 y 1981 – 2010 de las estaciones Amecameca de Juárez DGE (15007) y San Pedro Nexapa (15103); y con la serie de 1971 – 2000 para las estaciones Hueyapan E-4 (17046), Alponocan (17060), Atlixco (21012), Huejotzingo (21046), C.A.E. San Juan Tetla (21164), San Andrés Calpan (21167), San Pedro Benito Juárez E-1 (21193) y San Juan Tetla Siete Norte (21214). También ha facilitado las cantidades medias y extremas diarias de precipitación y los valores medios y extremos mensuales de precipitación referidos a las dos primeras estaciones para todo el periodo de estudio (2001 ‒ 2008); así como las cantidades máximas diarias de precipitación en el periodo 2001 – 2004 de las estaciones Hueyapan E-4, Alponocan y Atlixco; y del año 2001 de las estaciones San Andrés Calpan y San Pedro Benito Juárez E-1.

Para relacionar las precipitaciones máximas en 24 horas registradas en las estaciones cercanas con la generación de algún lahar secundario en la barranca de Huiloac se han tomado los casos estudiados por Muñoz-Salinas (2007) y Andrés (2009), donde se vinculan cantidades máximas de precipitación de entre 16.5 y 90.1 mm recogidas en estaciones cercanas, con el registro de lahares secundarios en Huiloac. Ambas autoras establecen además una relación entre el patrón de distribución espacio-temporal que presentan las precipitaciones y la formación del evento lahárico.

Por último, para las estaciones con datos disponibles se han relacionado las cantidades anuales y mensuales con las medias de la serie normal de referencia (1971 – 2000), con el fin de establecer años y meses secos y húmedos, según estén por encima o por debajo de la media de la serie. En los casos en los que no se tienen datos directos se ha recurrido a los informes del North America Drought Monitor (NADM) (2009, en línea) (Andrés, 2009).

 

4. Resultados

4.1. Dinámica erosiva/sedimentaria en Huiloac

Los resultados derivados de la aplicación de la metodología propuesta para el estudio de las variaciones topográficas en el cauce de la garganta Huiloac en términos de erosión y sedimentación se detallan en la Tabla 2.

Entre febrero de 2002 y octubre de 2002 todos los tramos en los que se dividió el área de estudio de la garganta Huiloac presentan valores superiores de erosión. El balance de material perdido entre ambas fechas es de 677.31 m³, lo que arroja una relación de 1.31 m³/mes. Los valores de sedimentación son más modestos, con una cifra absoluta de 145.45 m³ (0.28 m³/mes).

El primer subtramo (S1) presenta los valores más bajos de erosión y sedimentación (0.73 m³/mes y 0.06 m³/mes, respectivamente), influidos por el hecho de que en el P1 no se registraron cambios topográficos entre las dos fechas de estudio. Aguas abajo los valores de erosión y de sedimentación se distancian y, mientras que las cifras de deposición se muestran más o menos estables, en torno a 0.28 m³/mes; las de erosión presentan fuertes fluctuaciones. De esta forma, en los subtramos 2, 3, 8, 17 y 24 se pierden cantidades de 1.77 m³/mes, 1.74 m³/mes, 3.66 m³/mes, 2.02 m³/mes y 2.49 m³/mes, respectivamente; y los valores más bajos se registran en el S22, con 0.54 m³/mes, y en los subtramos del 11 al 14 (0.67; 0.75; 0.70 y 0.71 m³/mes).

Entre septiembre de 2003 y marzo de 2006 los valores estimados de pérdidas y aumentos de volumen en el tramo de estudio de la garganta Huiloac resultan muy similares, con 620.03 m³ y 621.22 m³, respectivamente. Esto indica tasas de erosión y sedimentación de 1.20 m³/mes. Pero estos valores totales no suponen una distribución uniforme de las variaciones, sino que esconden importantes fluctuaciones a lo largo del recorrido del cauce. Se pueden observar tres sectores donde predominan claramente las variaciones positivas de volumen (S4-S6, S12-S15 y S24-S29), y otros dos en los que son mayores las cantidades erosionadas (S2-S3 y S16-S22). En el resto del trayecto, la diferencia entre erosión y sedimentación es muy pequeña.

La erosión y la sedimentación producidas en el S1 presentan valores similares (0.66 m³/mes para la erosión y 0.65 m³/mes para la sedimentación), pero aguas abajo aumenta la erosión y disminuye la sedimentación hasta el S3 con 1.54 m³/mes erosionados, frente a 0.23 m³/mes sedimentados. Después los valores se invierten y en el S5 la sedimentación alcanza 3.40 m³/mes, frente a una cantidad erosionada de 0.45 m³/mes. Entre los perfiles 6 y 12 las cantidades erosionadas y depositadas oscilan ligeramente sin presentar picos destacados. Más adelante, ambas cantidades aumentan y en el S13 se estima que se erosionan 2.01 m³/mes y se depositan 3.24 m³/mes. A partir de aquí hasta el S23 los valores de erosión disminuyen ligeramente en los dos subtramos siguientes, para luego incrementarse y alcanzar el valor máximo registrado en el sector de estudio, en el S21 con 2.51 m³/mes. Por el contrario, las cantidades sedimentadas son cada vez menores hasta el S20, donde tan sólo se registran 0.04 m³/mes. Aguas abajo de los subtramos 20 y 21, las cantidades de material perdido y depositado tienden a igualarse, de forma que en el S23 alcanzan las cifras de 0.29 m³/mes y 0.25 m³/mes, respectivamente. Desde aquí hasta el punto más bajo del estudio (P1) el cauce se rellena con cantidades cada vez mayores de material, hasta alcanzar una tasa de 4.39 m³/mes en el S28. Por su parte, la erosión presenta valores bajos entre los perfiles 24 y 28 (0.28 m³/mes a 0.48 m³/mes) y un pequeño aumento en el S28 (0.94 m³/mes).

Entre marzo de 2006 y febrero de 2008 se registran las cantidades erosionadas más elevadas (1423.94 m³), lo que supone una tasa de 2.75 m³/mes. Suponiendo una erosión uniforme a lo largo del tiempo, el tramo considerado habría perdido 61.91 m³ de material al mes. Por otra parte, la sedimentación fue de 277.07 m³, lo que equivale a una tasa de 0.54 m³/mes y a una sedimentación mensual en todo el tramo de 12.05 m³.

En cuanto a la distribución a lo largo del espacio considerado, se pueden observar dos trayectos de mayor longitud en los que dominan valores elevados de erosión (S4-S16 y S19-S29), separados por intervalos más cortos, en los que domina ligeramente la sedimentación. También el espacio integrado por S1 y S2 presenta valores superiores de erosión, aunque no suponen un máximo.

Entre los perfiles 1 y 17 domina la erosión, excepto en el S3 donde la sedimentación presenta una tasa de 1.60 m³/mes, frente a 1.39 m³/mes de erosión. En este intervalo las cantidades erosionadas presentan dos máximos en S5 (2.09 m³/mes) y S10 (4.68 m³/mes); y un mínimo de 0.22 m³/mes en el S7. Los valores de sedimentación alcanzan sus máximos en los subtramos 3 (1.60 m³/mes) y 7 (0.59 m³/mes), para luego descender hasta la cifra de 0.06 m³/m en el S11. Desde aquí aguas abajo, la sedimentación aumenta ligeramente hasta el perfil 17.

Entre los perfiles 17 y 19 hay un cambio en la tendencia y la sedimentación supera a la erosión (2.22 m³/mes sedimentados en el S17, frente a 0.85 m³/mes erosionados). Desde el perfil 17 al 29 la sedimentación registra valores uniformes muy bajos tanto en cifras absolutas (< 5.35 m³) como relativas (< 0.32 m³/mes); mientras que las cantidades erosionadas aumentan aguas abajo (desde 2.03 m³/mes en el S19, hasta 6.21 m³/mes en el S28).

Entre septiembre de 2003 y febrero de 2008 se pierde un total de 1523.22 m³ de material en el tramo de estudio y se registra un aporte de 387.22 m³. En conjunto, la tasa de erosión fue de 2.95 m³/mes, mientras que la de sedimentación fue de 0.75 m³/mes. Las variaciones negativas de volumen son superiores a la deposición en todos los subtramos excepto en el S5, el S6 y el S17, donde la variación positiva es ligeramente superior.

Desde el punto más elevado hasta el P5, las tasas de erosión se mantienen en torno a 2 m³/mes, mientras que los valores de sedimentación experimentan un incremento desde 0.30 m³/mes en el S1 hasta 2.74 m³/mes en el S5, el valor más elevado de deposición estimado. Aguas abajo, las cantidades sedimentadas experimentan oscilaciones con dos máximos en el S13 (2.50 m³/mes) y en el S17 (1.40 m³/mes), hasta alcanzar su valor más bajo en el S21 (0.03 m³/mes). Desde aquí hasta el último perfil, la sedimentación tiende a aumentar hasta alcanzar los 2.13 m³/mes en el S28.

La amplitud que presentan las tasas de erosión es mucho mayor a lo largo de trayecto de estudio. El valor máximo se estima para el S22, con una cifra absoluta de 188.39 m³ y una tasa de 5.98 m³/mes, que disminuye ligeramente aguas arriba (en el S21 se estima una tasa de 5.97 m³/mes). Otros dos picos de erosión se localizan en el S10 (5.43 m³/mes) y en el S28 (4.76 m³/mes). Por el contrario el valor más bajo de erosión se registra en el S6 (0.64 m³/mes), con mínimos secundarios en S17 (0.73 m³/mes) y S26 (1.17 m³/mes).

 

4.2. Variaciones geomorfológicas en Huiloac

Como primera aproximación a la evolución morfológica de la garganta se analizan las variaciones de superficie que presentan las distintas formas cartografiadas en las seis fechas de estudio (Tabla 3 y Figura 7).

Según los datos, se aprecia que las formas más estables en cuanto a superficie ocupada son los surcos (A1) y las crestas de los levées (A2). Los cambios más significativos que se observan en la variación de la superficie ocupada por las distintas formas son los experimentados por las laderas (A3 y A4) y por el lecho fluvial (C1). En especial, destaca la disminución de las laderas superiores (A3), que pasan de constituir el 41.2 % de la superficie total en febrero de 2002, al 25.6 % en febrero de 2008. Las laderas inferiores (A4) también sufren pérdidas y así, disponen del 4.8 % de la superficie en el año de inicio, mientras que seis años después sólo ocupan el 1.8 %. Incluso se observa una ligera disminución en las crestas (A2), reconocidas como el elemento más estable, entre las dos últimas fechas consideradas. Se aprecia que la reducción de las laderas superiores es constante entre febrero de 2002 y septiembre de 2003, cuando se pierde una superficie de 14.81 m²/mes. Entre septiembre de 2003 y febrero de 2004, el retroceso de la ladera es mucho más rápido y pierde 61.74 m²/mes. En el siguiente intervalo considerado se frena la disminución (4.98 m²/mes), para superar los valores iniciales en el último periodo (23.36 m²/mes).

Este acortamiento de las laderas sucede paralelo a un ensanchamiento del fondo de la garganta, ocupado por el lecho fluvial (C1) y las formas constituidas por depósitos laháricos recientes ‒terrazas (B1) y bancos (B2)‒. En conjunto, estas tres formas pasan de ocupar el 31.0 % de la superficie en febrero de 2002, al 52.7 % en febrero de 2008. El lecho fluvial aumenta considerablemente su extensión desde la fecha de inicio hasta septiembre de 2003 (23.0 %, 21.0 % y 32.7 % en las tres fechas respectivas), mientras que las terrazas laháricas tienen escasa representación y se reducen con el tiempo (4.8 %, 3.6 % y 2.8 %). Sin embargo, a partir de septiembre de 2003 la tendencia se invierte y el fondo fluvial se reduce a favor de las terrazas laháricas y de los bancos fluvio-laháricos, de manera que en las dos fechas de muestreo siguientes el lecho ocupa el 25.7 % y 22.7 % de la superficie de la garganta respectivamente, y las nuevas formas se incrementan hasta alcanzar el 18.9 % y el 25.2 %, conjuntamente. Por último, en febrero de 2008 el fondo fluvial ocupa la extensión mayor en los seis años de estudio (41.8 % del total), mientras que las terrazas y los bancos experimentan una reducción hasta el 18.2 %.

Ya que los periodos de estudio no son iguales, se han calculado las tasas de variación en m²/mes que se recogen en la Tabla 4. Se observa que la morfodinámica actuante entre febrero de 2002 y febrero de 2008 en este tramo de la garganta de Huiloac ha tendido al ensanchamiento del lecho y al retroceso de las laderas, con el consiguiente aumento de la pendiente de las mismas. Dentro de esta línea general, se aprecia que las laderas retroceden a un ritmo similar en los periodos de 2002-04 y 2006-08, mientras que en los dos años intermedios se ralentizó el proceso. Por su parte, las formas que ocupan el fondo (lecho, bancos fluvio-laháricos y terrazas laháricas) aumentaron su superficie en 45.36 m²/mes en el primer periodo; en 10.26 m²/mes, en el segundo y en 16.86 m²/mes, entre 2006 y 2008. El lecho tiende a incrementar su superficie durante el primer y tercer bienio, especialmente en este último (40.83 m²/mes), y por el contrario retrocede entre 2004 y 2006 (-9.53 m²/mes). El incremento de superficie de las terrazas y bancos fue mayor en los dos primeros años (36.25 m²/mes), que en el bienio siguiente (19.80 m²/mes), para luego comenzar a disminuir a partir de 2006 a razón de -23.96 m²/mes.

El análisis espacial realizado a partir de las 6 cartografías geomorfológicas entre febrero de 2002 y febrero de 2008 (Figura 8) ofrece una combinación de modificaciones que se pueden cuantificar según la superficie planimétrica a la que afectan. Así, entre 2002 y 2008 (Tablas 5 y 6) el 37.34 % de la extensión del tramo analizado experimenta variaciones morfológicas. Las modificaciones del 62.84 % del área con cambios (23.37 % del área total) afectan a las laderas formadas sobre los depósitos laháricos de 1997 y 2001. Tan sólo el 2.01 % de la superficie de las laderas (99.69 m²) es atacada por la erosión, pero se sigue manteniendo como ladera lahárica (A4), mientras que el 36.21 % es sustituido por otras formas. Las laderas retroceden principalmente a favor del lecho fluvial (12.78 % de la superficie total) y de las formas labradas en depósitos laháricos posteriores a enero de 2001 (8.82 % de la superficie total).

Por lo tanto, en la medida en que retroceden las laderas laháricas ceden espacio a las formas que se desarrollan en el fondo del canal: lecho fluvial, bancos fluvio-laháricos y terrazas laháricas. Pero se observa que tanto bancos fluvio-laháricos como terrazas laháricas, muestran una gran variabilidad espacial, ya que sólo mantienen el 20.98 % y el 27.19 %, respectivamente, de sus superficies originales, mientras que la superficie que ocupan en 2008 duplica a la que tenían en 2002 (Tabla 2). Esto se debe a que, si bien el 70 % de su área desaparece por acción fluvial (5.46 % del total), invaden espacios cedidos por las laderas (702.21 m², 8.82 % del total) y por el propio lecho (557.42 m², 7.00 % del total).

El lecho fluvial también incrementa su presencia espacial en la garganta, del 23.0 % al 34.5 % (Tabla 2). Un 70.84 % de la superficie cartografiada como lecho en febrero de 2002 mantiene esta denominación en 2008 (Tabla 6), mientras que el resto viene a ser ocupado por las nuevas formas laháricas (terrazas y bancos). El nuevo lecho fluvial se extiende a expensas del espacio cedido por el retroceso de las laderas, de forma que el 35.78 % de la superficie ocupada por el lecho en febrero de 2008 pertenecía en la fecha inicial a los depósitos de los lahares de 1997 y 2001.

Las formas menos afectadas por los cambios son las laderas internas de la garganta constituidas por depósitos de los lahares de 1997 y 2001. El 61.54 % de la garganta que no ha experimentado cambios (38.56 % del total) (Tabla 5) está cubierto por estas formas (A2, A3 y A4). En especial, las crestas son las formas más estables, ya que el 94.85 % de la superficie que presentaban en febrero de 2002 se mantienen en febrero de 2008. Las laderas superiores (A3) conservan el 61.66 % de la superficie inicial, mientras que las más afectadas son las laderas inferiores (A4), que sólo mantienen el 12.48 % de la extensión original (Tabla 6). Por otra parte, el 26.28 % de la superficie sin alterar pertenece al lecho fluvial (Tabla 5), que tan sólo cede a otras formas el 29.16 % de su superficie inicial.

Mención especial merecen las terrazas laháricas (B1) de las que tan sólo permanece en 2008 el 27.19 % de la superficie cubierta en 2002 (Tabla 6).

 

4.3. Evolución geomorfológica de Huiloac tras los lahares primarios

Por un lado, el análisis de las variaciones topográficas ha proporcionado información sobre qué sectores de los perfiles transversales se ven afectados por la erosión, la sedimentación o permanecen estables, además de aportar cantidades estimadas de material erosionado y sedimentado en un determinado periodo de tiempo. Por otra parte, la comparación de las cartografías geomorfológicas permite cuantificar las variaciones de superficie que presentan las distintas formas cartografiadas. A partir de estos datos se puede inferir la evolución morfológica (Figura 9) que ha experimentado el tramo de estudio desde el último lahar primario.

a) Los trabajos realizados en la primera fecha (febrero de 2002) muestran que desde el último lahar primario de enero de 2001 el cauce ya ha sufrido importantes modificaciones. En las observaciones de campo se encontraron los depósitos laháricos de 2001 disectados, a la vez que se pudieron cartografiar terrazas laháricas (B1) recientes (Figura 9.A). Esta intensa actividad geomorfológica ocurrida en el año posterior al flujo primario, tuvo que estar determinada por las importantes precipitaciones de marzo, junio, septiembre del 2001. El CENAPRED (2001) informó de la existencia de lahares secundarios en marzo de 2001 y Andrés (2009) analiza minuciosamente la relación entre las cantidades recogidas en los observatorios y este flujo lahárico. Por otra parte, las copiosas nevadas de enero de 2002 proporcionaron importantes caudales a un canal rellenado de sedimentos laháricos (Andrés, 2009).

b) Entre febrero y octubre de 2002 los perfiles longitudinales del cauce indican que se produce una incisión generalizada del fondo (los puntos más bajos de los perfiles transversales descienden una media de 0.53 m), con una mayor profundización en los perfiles del 24 al 29. Tan sólo se observa que entre los perfiles 5 y 8 el punto más bajo sufre una elevación (máxima de 0.5 m). Por otra parte, en los perfiles transversales se observa que domina la erosión frente a la sedimentación (excepto en el P23). De esta manera se estima una pérdida de 677.31 m³ de material y un incremento de 145.45 m³entre los dos momentos de muestreo (Tabla 2).

Entre estas dos fechas, los cambios morfológicos (sólo el 15 % de la superficie) se centran en el fondo del cauce, donde adquieren especial desarrollo los bancos fluvio-laháricos (B2), mientras que las terrazas laháricas se ven reducidas. También se aprecia un retroceso en las laderas inferiores (A4), constituidas por los depósitos laháricos de 1997 y 2001, que dejan su lugar al lecho fluvial o nuevos depósitos de lahar.

Se produce entonces una apertura de la sección del canal, ya que domina la evacuación de los materiales, frente a la deposición. El fondo de la garganta se ensancha en una tendencia a formar perfil en "U" cada vez más abierto (Figura 9.B). La razón del predominio de la actividad erosiva hay que buscarla en los copiosos aguaceros del mes de septiembre de 2002 (con cantidades superiores a la media de la serie) y las elevadas precipitaciones del 9 de octubre (28.5 mm en 24 horas en la estación de Amecameca, por ejemplo), justo antes de la realización del trabajo de campo (Andrés, 2009).

c) Entre octubre de 2002 y septiembre de 2003, las laderas labradas en los materiales laháricos depositados en 1997 y 2001 se ven afectadas por procesos erosivos y reducen su extensión en un 5 %. El resto de las modificaciones se concentran en el fondo del cauce donde reducen su presencia depósitos fluvio-laháricos y terrazas (Figura 9.C).

La comparación de los perfiles transversales entre estas dos fechas es arriesgada, ya que se modificó la instrumentación en la obtención de datos topográficos. De todas formas, las observaciones realizadas en los trabajos de campo señalan que el fondo del cauce sufrió un relleno, lo que en parte también puede explicar la reducción de la superficie ocupada por las laderas formadas en los depósitos de los lahares primarios, de los bancos fluvio-laháricos y de las terrazas laháricas formadas por flujos recientes. Entre las dos fechas, tan solo los meses de noviembre de 2002 y junio de 2003 recibieron cantidades de precipitación ligeramente superiores a la media y el NADM considera 2003 como un año seco. Por otra parte, las cantidades registradas en 24 horas en el periodo resultan bastante modestas y existe poca probabilidad de que se haya generado un caudal elevado en la barranca o al menos con capacidad para seguir erosionando en el tramo medio de Huiloac, donde se encuentra el sector de ensayo (Andrés, 2009).

d) La situación de febrero de 2004 muestra notables cambios respecto de la de septiembre de 2003 (Figura 8). Un 22.8 % de la superficie de estudio ha sido afectada por procesos de erosión y deposición. Destaca la reducción en superficie que experimentan los depósitos de los lahares primarios de 1997 y 2001 (7.6 % del total del área analizada), en especial las laderas superiores, que empiezan a modificar su morfología al ser perturbadas por procesos erosivos. En cambio, el lecho se ensancha notablemente, a la vez que adquieren mayor desarrollo las formas asociadas a las dinámicas de flujos recientes (B1 y B2). En febrero de 2004, el lecho fluvial, los bancos fluvio-laháricos y las terrazas laháricas recientes constituyen el 48 % de la superficie de estudio, mientras que en septiembre de 2003 ocupaban el 39 % (Tabla 3). En particular, las terrazas laháricas aumentan su presencia (8.8 %).

En conjunto los cambios sucedidos entre febrero de 2002 y febrero de 2004 se pueden resumir en la erosión sufrida por las laderas constituidas por los depósitos de los lahares primarios de 1997 y 2001 y en la formación de nuevas terrazas laháricas (Figura 8). Las nuevas terrazas se generan directamente sobre el fondo del lecho (5.92 % de la superficie total) o restando espacio a las laderas inferiores (3.84 % de la superficie total). La actividad erosiva registrada en la barranca parece ser más elevada que la esperada para una temporada seca. De hecho, resultó ser una estación seca anormalmente húmeda, con una anomalía del 147 % en las precipitaciones de octubre de 2003 (NADM), máximos en 24 horas elevados (57 mm registrados en el observatorio de Amecameca el 28 de septiembre, por ejemplo) y precipitaciones copiosas en enero de 2004, en forma de nieve en los sistemas montañosos cercanos a Ciudad de México (NADM).

e) Las transformaciones geomorfológicas que tienen lugar en la garganta son menos apreciables entre febrero de 2004 y marzo de 2006, con una superficie de cambio que ocupa el 19% del área total. El retroceso de las laderas laháricas es sensiblemente inferior a los dos años anteriores (3.67 % de la superficie total) y el espacio que liberan es ocupado en un 72.21 % (210.80 m²) por nuevas terrazas y bancos fluvio-laháricos, que además restan al lecho otros 666.53 m². En marzo de 2006, los bancos fluvio-laháricos y las terrazas ocupan el 25.2 % de la garganta, lo que supone la máxima extensión para los bancos (16.2 % de la superficie total) y un aumento de 193.90 m²para las terrazas.

La comparación de los datos topográficos disponibles indica que para el total del tramo estudiado se produjeron cantidades muy similares de erosión (620.03 m³) y de sedimentación (621.22 m³) entre septiembre de 2003 y marzo de 2006 (Tabla 2). Pero erosión y sedimentación no se reparten de forma uniforme a lo largo del tramo de garganta, sino que presentan concentraciones en distintos lugares. Así, entre los perfiles 16 y 24 predomina la pérdida de material, con el máximo valor de erosión registrado en el S21 (2.51 m³/m). Se dan aquí también los valores más altos de incisión en el cauce con una variación de -1.24 m y de -0.76 m en los puntos más bajos de los perfiles 20 y 21, respectivamente. En el análisis detallado de estos perfiles se observa que entre las dos fechas se produce un retroceso de las laderas inferiores, lo que se traduce en el predominio de procesos de zapa. En especial, los perfiles 22 y 23 experimentan retrocesos en sus laderas, debido a la erosión concentrada en la base.

También son más elevados los valores de erosión que de sedimentación en los subtramos 1, 2, 3, 7, 10 y 11, donde el talweg presenta también incisiones de -1.60 m (P2), -0.59 m (P3), -0.73 m (P10). Los perfiles que enmarcan estos subtramos indican una profundización en el lecho, además de un proceso de zapa incipiente en la margen izquierda del perfil 11. El descenso en altura de todo el fondo resulta imperceptible en los mapas de variaciones geomorfológicas, ya que el lecho fluvial excavado sigue ocupando una superficie similar, por lo que no se registran cambios. Incluso la cartografía morfológica puede dar lugar a equívocos, ya que el fondo del lecho fluvial queda ocupado por nuevas terrazas laháricas en las orillas, con lo que da una falsa idea de dominio de la sedimentación.

La sedimentación predomina en el resto de los subtramos considerados y alcanza el máximo valor absoluto en el S13, aunque en relación con la longitud del canal tiene mayor incidencia entre los subtramos 24 y 28. Es en este sector donde se produce un relleno del fondo del canal de aproximadamente 1 m de espesor, que entierra las laderas inferiores y parte de las superiores. Además va acompañado de un ligero retroceso de las laderas superiores, con lo que el resultado final es una elevación del fondo y un ensanchamiento del canal, donde apenas resaltan en altura los bancos fluvio-laháricos localizados en las orillas del lecho fluvial.

El comportamiento en el patrón de precipitaciones es bastante diferente al observado hasta ahora. Aunque 2004 se considera levemente húmedo y registra máximos en 24 horas capaces de generar lahares competentes desde el punto de vista erosivo en el tramo de estudio (52.3 mm el 29 de junio de 2004 en San Pedro Nexapa), desde noviembre de 2004 a marzo de 2006 se generaliza una situación de sequía marcada, cuando se registran los meses más secos desde 1941 (NADM). Estas condiciones anormalmente secas parecen ser la causa de un balance erosión/sedimentación muy igualado al final del periodo.

f) En el último bienio de estudio, entre marzo de 2006 y febrero de 2008, los cambios geomorfológicos se dieron en el 31.12 % de la superficie de estudio. El retroceso de las laderas se incrementa respecto de la fase anterior y afecta al 7.51 % de la superficie total, de manera que la erosión llega a perturbar las crestas laháricas. Ahora son frecuentes los deslizamientos que se producen en unas laderas cada vez más verticales. También reducen su presencia las formas laháricas recientes distribuidas por el fondo de la garganta.

Se ha calculado que el tramo de estudio ha perdido 1423.94 m³ de material, mientras que sólo ha aumentado 277.07 m³ (Tabla 2). Entre los perfiles 1 y 21 el fondo se ha alimentado con el material erosionado a los depósitos de los lahares primarios, cuyas laderas han sido afectadas por deslizamientos. En cambio, entre los perfiles 22 y 29 domina la incisión en el fondo que ha vaciado este sector del relleno que sufriera en el bienio anterior, borrando así los depósitos de terrazas recientes.

Este nuevo incremento de la actividad erosiva parece estar estrechamente relacionado con un aumento de las precipitaciones en general. Por influencia del fenómeno El Niño (ENSO), las condiciones de humedad aumentan a partir de marzo de 2006 y sólo se ven frenadas en una ligera sequía en la primavera de 2007, que rápidamente es superada por un nuevo incremento en las precipitaciones hasta el final del periodo (NADM y tratamiento de datos del SMN en Andrés, 2009). Según esto se puede inferir que en los años más húmedos que la media, los flujos se muestran más competitivos, no sólo en erosionar sino en evacuar lo erosionado o depositado en el tiempo transcurrido entre las fechas de muestreo. Esto se refuerza con el hecho de que las laderas labradas en los depósitos laháricos presenten fuertes pendientes y, por lo tanto, una alta inestabilidad.

De forma independiente, los grandes bloques (E1) se han interpretado como puntos de referencia en la evolución geomorfológica de la garganta. Su posición estable sobre el fondo es indicadora de que la garganta ha sido recorrida por flujos menores, sin competencia para movilizar grandes tamaños. Por otra parte, un aumento de la presencia de bloques en el fondo se puede explicar como consecuencia de aportes desde las laderas por procesos de zapa, tal como ocurre entre febrero 2002 y febrero 2004. En otras ocasiones, su representación disminuye en relación con el dominio de procesos de sedimentación, que cubren buena parte de los bloques (periodo entre febrero 2004 y marzo 2006). Por último, entre marzo de 2006 y febrero de 2008 se observa un aumento de estos elementos en el fondo del cauce que nos da idea de la capacidad de los flujos que circularon durante esas fechas y vaciaron el cauce de los tamaños inferiores.

 

5. Discusión

Aunque tras los flujos laháricos de 1997 y 2001, el volcán Popocatépetl se ha mantenido en un periodo de calma eruptiva en el que no se han producido más lahares primarios, el seguimiento geomorfológico y topográfico llevado a cabo en este estudio muestra una importante actividad geomorfológica en la garganta Huiloac entre 2002 y 2008. Esta actividad se caracteriza fundamentalmente por procesos erosivos que evacuan los depósitos de lahares primarios, consiguiendo así un ensanchamiento y profundización del cauce, y por procesos sedimentarios que se reflejan en la formación de nuevos modelados de terrazas que colmatan el fondo del cauce. Los flujos de pequeños lahares secundarios que a menudo discurren por la garganta, del tipo dilute sediment-laden flow (Capra et al., 2010) o del tipo muddy stream flow (van Westen y Daag, 2005), son los causantes de estos cambios morfológicos.

En este sentido, Tanarro et al. (2010) establecen una relación entre las precipitaciones registradas en las proximidades del área de estudio y los momentos en los que se producen mayores cambios morfológicos en la garganta Huiloac después del lahar primario de 2001 hasta marzo de 2006. Estos autores determinan dos situaciones muy favorables para las modificaciones: a) durante la estación seca y al inicio de la estación húmeda con umbrales de precipitación inferiores a 30 mm/día, y b) al final de la estación húmeda, con aguaceros en uno o varios días consecutivos (alrededor de 70 mm/día). De acuerdo con estas premisas, encontramos que las copiosas precipitaciones de marzo, junio y septiembre del 2001, junto con las extraordinarias de enero de 2002 (Andrés, 2009), aportaron importantes caudales a un canal rellenado de sedimentos laháricos y vaciaron el fondo del cauce de la mayor parte de dichos sedimentos, como lo demuestra la morfología de febrero de 2002, donde sólo quedan algunos resaltes aislados, a modo de pináculos en el fondo, y el resto de material lahárico adosado a las laderas. De nuevo se encuentran notables cambios en octubre de 2002, ya que los flujos han vaciado todo el fondo de los depósitos laháricos primarios y han disectado las formas encontradas en el mes de febrero. En esta ocasión, las precipitaciones de la estación de lluvias fueron claramente inferiores a la media (Andrés, 2009), pero se recogieron aguaceros continuados en los meses de septiembre y octubre. Los cambios topográficos y geomorfológicos registrados en febrero de 2004 (erosión de las laderas constituidas por los depósitos de los lahares primarios de 1997 y 2001 y formación de nuevas terrazas laháricas sobre el fondo), se relacionan con la competitividad de las precipitaciones extraordinarias de enero de 2004 (Andrés, 2009; Tanarro et al., 2010), mientras que durante la sequía de 2003 los cambios fueron menores, como se observa en las cartografías de septiembre de 2003.

Entre febrero de 2004 y marzo de 2006 se ralentizan los cambios y se frena considerablemente la erosión de los depósitos de los lahares primarios que constituyen las laderas del cauce en relación con una fase seca. Los máximos de precipitación coinciden con las estaciones húmedas y, aunque alcanzan valores suficientemente altos como para alimentar lahares secundarios, estos no resultaron tan competitivos como para producir grandes cambios en el modelado del cauce. Sin embargo, entre marzo de 2006 y febrero de 2008, la capacidad erosiva de las avenidas se incrementa (1423.9 m³, Tabla 2) y se registra la mayor superficie afectada por cambios geomorfológicos (31.1 % del total). En este periodo, se supera la fase de sequía anterior y las cantidades de precipitación registradas superan ligeramente la media de la serie normal, excepto las de los meses de septiembre y octubre de 2007, que reciben importantes aguaceros ocasionados por influencia del fenómeno El Niño (ENSO) (NADM). Esta observación corrobora que las precipitaciones copiosas al final de la estación húmeda resultan especialmente competitivas en introducir cambios en el modelado de los cauces laháricos, como indican Tanarro et al. (2010).

Capra et al. (2010) advirtieron una situación similar para los flujos del volcán de Colima, donde los lahares son más frecuentes al inicio de la estación seca, con pequeñas cantidades acumuladas de precipitación (< 10 mm), que los causados por los picos más elevados de la estación húmeda (> 70 mm/h). Los autores explican esta situación aparentemente contradictoria por el "mecanismo de repelencia de agua" (Doerr et al., 2000), asociado a la presencia en los canales de acículas producidas por las coníferas que pueblan las laderas del volcán. Las acículas, cuya acumulación sobre el suelo es muy elevada en la estación seca y al inicio de la estación húmeda, facilitan la escorrentía y la consecuente formación de lahares secundarios (Capra et al., 2010). La misma causa se puede aplicar para explicar la competencia de los flujos originados por precipitaciones en la estación seca en la garganta Huiloac del volcán Popocatépetl, mientras que la capacidad de modificación de los flujos generados a finales de la estación húmeda, se debe a que el suelo se encuentra saturado y la mayor parte de las aguas precipitadas pasan a la escorrentía directamente.

Además de interesarse por las relaciones entre las precipitaciones, la formación de lahares y las modificaciones morfológicas en los cauces, los científicos también han prestado atención a determinar el tiempo que tarda la red de drenaje de los volcanes en recuperar el equilibrio geomorfológico después de ser afectada por el paso de lahares sin-eruptivos (Chow, 1959; Pierson, 1985; Pierson y Costa, 1987; Daag y van Westen, 1996; van Westen, 1997; Rodolfo, 1989; Rodolfo et al., 1989). Algunos de estos autores convienen en señalar que las mayores variaciones geomorfológicas se dan tras periodos de relativa calma eruptiva (Chow, 1959; Rodolfo, 1989; Pierson, 1985; Pierson y Scott, 1985). A partir del caso del estudio del St. Helens, se han propuesto cuatro etapas en la recuperación de los cauces tras los lahares primarios: iniciación del canal, incisión, agradación/ensanchamiento y estabilización (Meyer y Martison, 1989; Major, 2003). Además, Gran y Montgomery (2005), en su estudio de la red fluvial del Pinatubo, apuntan que la recuperación fluvial comienza con un descenso en la entrada de sedimentos a los cauces, ya que el área fuente de sedimentos está estabilizada o ha sido ya evacuada.

En el caso que nos ocupa, el volcán Popocatépetl entró en una fase de calma latente posterior a la erupción de enero de 2001, que causó el último lahar primario hasta la fecha. A partir del análisis combinado de la evolución geomorfológica y topográfica post-lahárica realizado en este trabajo se pueden distinguir distintas fases en la morfodinámica reciente de la garganta Huiloac. En el primer momento de estudio, febrero de 2002, se observan huellas de incisión en los materiales de los lahares de 1997 y 2001, por lo que ya habría comenzado la fase de incisión del canal (según la terminología empleada en Meyer y Martinson, 1989; y Major, 2003). En febrero de 2002 todavía quedan resaltes de estos depósitos aislados en el fondo del cauce, que han desaparecido en octubre de 2002. A partir de entonces, como indica la sucesión de perfiles, el cauce empieza una dinámica de ensanchamiento y profundización a expensas de los depósitos de los lahares primarios que todavía permanecen adosados en los flancos, formando las laderas del cauce. Estos procesos están asociados con corrientes de agua someras ligeramente cargadas de sedimentos (Meyer y Martinson, 1989; Hayes et al., 2002; Manville et al., 2005) y con pequeños lahares secundarios originados por precipitaciones (Barclay et al., 2007), por lo que las etapas de incisión y deposición se van sucediendo en función del tipo de flujo que discurra por el cauce. Así, durante el periodo 2004 ‒ 2006 hemos encontrado una menor área afectada por modificaciones geomorfológicas y valores de sedimentación y erosión muy similares (Tabla 2) en relación con cantidades de precipitación inferiores a la media y una distribución temporal poco favorable para la generación de lahares secundarios (sin precipitaciones en la estación seca ni picos a finales de la estación húmeda). En cambio, en el bienio 2006 ‒ 2008 una mayor superficie se ve afectada por cambios y dominan los procesos de profundización y ensanchamiento, de manera que la erosión llega a perturbar las crestas del material depositado por los lahares primarios. Estos valores de erosión más elevados se relacionan con un ligero aumento de las precipitaciones recibidas (ligeramente superiores a la media) y, sobre todo, con una prolongación de la estación húmeda debido a una situación particular (la presencia del ENSO). A esto hay que sumar la inestabilidad que han ido adquiriendo las laderas labradas en los depósitos de los lahares primarios por aumento de su pendiente, debido a los procesos de zapa.

En la situación de febrero de 2008 todavía no se puede afirmar que el cauce de la garganta Huiloac haya alcanzado una fase de estabilidad. Como se puede apreciar en los mapas geomorfológicos y en los perfiles, todavía quedan paquetes de sedimentos de los lahares primarios adosados en la parte inferior de las paredes y su disposición es muy inestable, ya que han adquirido una fuerte pendiente y están afectadas por procesos erosivos en su base. En el caso de que se produjeran fuertes precipitaciones, durante la estación seca o a finales de la estación húmeda principalmente, se daría una alta probabilidad de producción de lahares secundarios y habría que tener en cuenta la inestabilidad que presentan los depósitos de lahares primarios, más la presencia de las terrazas laháricas recientes (material suelto), a la hora de calcular los efectos de un flujo lahárico. Además, es recomendable que en futuros estudios se trataran otros aspectos que influyen en la formación y desarrollo de estos flujos, como puede ser la cantidad de materia vegetal muerta que se encuentra en el interior del cauce o las características de los depósitos del lahar primario de 2001 que todavía se mantiene disperso en la cabecera de la garganta (cantidad, disposición, permeabilidad...). Todas estas consideraciones se están teniendo en los trabajos que el mismo grupo de investigación está realizando durante el actual periodo de observación (2008 ‒ 2015).

 

6. Conclusiones

La interpretación combinada de los mapas geomorfológicos y los perfiles topográficos realizados en un tramo de la garganta Huiloac muestran que este cauce ha experimentado una constante transformación desde el suceso del lahar primario de enero de 2001. La metodología empleada no sólo ha servido para cuantificar las variaciones espacio-temporales de las distintas formas reconocidas en el cauce, sino también ha posibilitado relacionarlas con las variaciones topográficas. De este modo se establecen el tipo de procesos que han causado las modificaciones, su carácter erosivo o sedimentario, su distribución espacial y la evolución temporal; y se determina la situación actual en la que se encuentra el cauce ante un posible evento extremo.

La morfodinámica reciente del cauce de Huiloac comenzó con la incisión y evacuación del material que rellenó el cauce durante el lahar de 2001. Las fases iniciales (hasta octubre de 2002) se caracterizan por el vaciado de la parte central del cauce, mientras que las laderas inferiores quedan recubiertas por espesores variables depositados durante los lahares sin-eruptivos de 1997 y 2001. En etapas posteriores la acción geomorfológica de las aguas corrientes y la dinámica de laderas se centran fundamentalmente en ensanchar y profundizar el cauce, aunque se producen momentos intermedios de sedimentación en el fondo del canal, en los que se forman bancos fluvio-laháricos y terrazas laháricas, dependiendo del tipo de flujo.

La frecuencia y capacidad de los lahares secundarios que discurren por la garganta condicionan la evolución morfológica del cauce. A su vez, las precipitaciones son el factor determinante para la producción de estos flujos. Así, el tramo de estudio ha experimentado fases con menos variaciones morfológicas y cantidades inferiores de material erosionado (2004 ‒ 2006), coincidiendo con épocas ligeramente menos lluviosas que la media y con precipitaciones más regulares; mientras que las precipitaciones excepcionales en la estación seca (enero de 2002 y enero de 2004) y los aguaceros continuados a final de la estación húmeda (septiembre y octubre de 2007), agilizan los procesos de erosión por incisión y por zapa lateral, ocasionando así un vaciado de los depósitos de los lahares primarios.

En el periodo analizado se han observado procesos de ensanchamiento y profundización del cauce a expensas de los depósitos abandonados por los lahares primarios de 1997 y 2001, con los que alternan en segundo plano procesos de deposición y remoción que crean nuevas formas en el fondo del cauce. La intensidad de las transformaciones depende de las características de los flujos (especialmente su alimentación) y el grado de la inestabilidad de las laderas labradas en los depósitos laháricos. Por lo tanto, a pesar de haber transcurrido siete años desde el último lahar sin-eruptivo hasta la última fecha considerada en este trabajo, no se puede afirmar que el cauce haya alcanzado su estabilidad, sino que, a pesar de la considerable disminución en la cantidad de depósitos, estos se encuentran en una situación altamente inestable, debido a la fuerte pendiente que presentan. El trabajo que el mismo grupo de investigación está realizando en el actual periodo de observación (2008 ‒ 2015) proporcionará nuevos datos a la evolución morfológica de este tramo del cauce de Huiloac.

 

Agradecimientos

El presente trabajo se ha realizado con la financiación del proyecto CGL 2012-35858 (CRYOCRISIS) del Ministerio de Economía y Competitividad de España y del grupo de Trabajo BSCH/UCM: 931562 GEOGRAFÍA FÍSICA DE ALTA MONTAÑA. Los autores agradecen al Instituto de Geografía de la UNAM por su total apoyo en el trabajo de campo, así como el trabajo desinteresado y exhaustivo de numerosos de sus becarios en las tareas de monitorización.

 

Referencias

Ackerman, C.T., Evans, T.A., Brunner, G.W., 2000, HEC-GeoRAS: linking GIS to hydraulic analysis using ARC/ INFO and HEC-RAS, en Maidment, D., Djokic, D. (eds.), Hydrologic and hydraulic modeling support with Geographic Information System: Redlands, California, ESRI press, 155–176.         [ Links ]

Andrés, N., 2009, Técnicas de información geográfica aplicadas al estudio del origen de los lahares y su experimentación en estratovolcanes tropicales: Madrid, Universidad Complutense de Madrid, tesis doctoral, 476 p.         [ Links ]

Andrés, N., Zamorano, J.J., Sanjosé, J.J., Atkinson, A., Palacios, D., 2007, Glacier retreat during the recent eruptive period of Popocatépetl volcano, Mexico: Annals of Glaciology, 45, 73–82.         [ Links ]

Barclay, J., Alexander, J., Sušnik, J., 2007, Rainfall-induced lahars in the Belham Valley, Montserrat, West Indies: Journal of the Geological Society, 164, 815–827.         [ Links ]

Beaman, J.H., 1962, The timberlines of Iztaccíhuatl and Popocatépetl, México: Ecology, 43, 377–385.         [ Links ]

Capra, L., Poblete, M.A., Alvarado, R., 2004, The 1997 and 2001 lahars of Popocatépetl volcano (Central Mexico): textural and sedimentological constrains on their origin and hazards: Journal of Volcanology and Geothermal Research, 131, 351–369.         [ Links ]

Capra, L., Borselli, L., Varley, N., Gavilanes-Ruiz, J.C., Norini, G., Sarocchi, D., Caballero, L., Cortes, A., 2010, Rainfall-triggered lahars at Volcán de Colima, Mexico: surface hydro-repellency as initiation process: Journal of Volcanology and Geothermal Research, 189, 105–117.         [ Links ]

Centro Nacional de Prevención de Desastres (CENAPRED), 2001, Resumen de la actividad del volcán Popocatépetl de diciembre de 1994 a mayo de 2001 (en línea): México, disponible en <http://www.cenapred.unam.mx/es/Instrumentacion/InstVolcanica/MVolcan/Resumen/>, consultado 5 de junio de 2013.         [ Links ]

Centro Nacional de Prevención de Desastres (CENAPRED), 2013, Monitoreo volcánico, reportes diarios de la actividad volcánica desde agosto de 1997 (en línea), disponible en <http://www.cenapred.gob.mx>, consultado 9 de septiembre de 2013.         [ Links ]

Chen, J., He, Y.P., Wei, F.Q., 2005, Debris flow erosion and deposition in Jiangjia Gully, Yunnan, China: Environmental Geology, 48, 771–777.         [ Links ]

Chinen, T., Kadomura, H., 1986, Post-eruption sediment budget of a small catchment on Mt. Usu, Hokkaido: Zeitschrift für Geomorphologie N.F. Supplementband, 60, 217–232.         [ Links ]

Chow, V.T., 1959, Open-channel hydraulics: New York, McGraw-Hill Book Company, 680 p.         [ Links ]

Cronin, S.J., Neall, V.E., Lecointre, J.A., Palmer, A.S., 1999, Dynamic interactions between lahars and stream flow: A case study from Ruapehu volcano, New Zealand: GSA Bulletin, 111(1), 28–38.         [ Links ]

Cruz-Reyna, S. de la, Quezada, J.L., Peña, C., Zepeda, O., Sánchez, T., 1995, Historia de la actividad reciente del Popocatépetl (1354–1995), en Zepeda, O., Sánchez, T.A. (eds.), Volcán Popocatépetl: Estudios Realizados Durante la Crisis de 1994–1995, Secretaría de Gobernación, Sistema Nacional de Protección Civil, Centro Nacional de Prevención de Desastres (CENAPRED), México, D.F., Universidad Nacional Autónoma de México, 3–22.         [ Links ]

Daag, A.S., 1994, Geomorphic developments and erosion of the Mount Pinatubo 1991 pyroclastic flows in the Sacobia watershed, Philippines: A study using remote sensing and Geographic Information Systems (GIS): Enschede, Netherlands, International Institute for Aerospace Survey and Earth Sciences, tesis de maestría, 106 p.         [ Links ]

Daag, A.S., 2003, Modelling the Erosion of Pyroclastic Flow Deposits and the Occurrences of Lahars at Mt. Pinatubo, Philippines: Enschede, Netherlands, International Institute for Aerospace Survey and Earth Sciences, Universidad de Utrecht, tesis doctoral, 238 p.         [ Links ]

Daag, A., van Westen, C.J., 1996, Cartographic modelling of erosion in pyroclastic flow deposits of Mount Pinatubo, Philippines: ITC Journal, 1996-2, 110–124.         [ Links ]

Doerr, S.H., Shakesby, R.A., Walsh, R.P.D., 2000, Soil water repellency: its causes, characteristics and hydro-geomorphological significance: Earth-Science Review, 51, 33–65.         [ Links ]

García-Romero, A., 1998, Análisis integrado de paisajes en el occidente de la cuenca de México (la vertiente oriental de la Sierra de las Cruces, Monte Alto y Monte Bajo): Madrid, Universidad Complutense de Madrid, tesis doctoral, 600 p.         [ Links ]

Garcin, M., Poisson, B., Pouget, R., 2005, High rates of geomorphological processes in a tropical area: the Remparts River case study (Réunion Island, Indian Ocean): Geomorphology, 67, 335–350.         [ Links ]

González, A.E., 2000, Estudios de detalle estratigráfico y sedimentológico del Lahar de San Nicolás en el flanco noreste del volcán Popocatépetl: México, Universidad Nacional Autónoma de México, tesis de licenciatura, 109 p.         [ Links ]

González, A.E., Delgado, H., Urrutia, J., 1997, The San Nicolás Lahar at Popocatépetl Volcano (Mexico): a case study of a glacier-ice-melt-related debris flow, triggered by a blast at the onset of a plinian eruption (resumen), en Abstracts of International Association of Volcanology and Chemistry of the Earth's Interior, General Assembly, Puerto Vallarta, Jalisco, México, Gobierno de Jalisco, Secretaría General, Unidad Editorial, 94.         [ Links ]

Gran, K., Montgomery, D., 2005, Spatial and temporal patterns in fluvial recovery following volcanic eruptions: Channel response to basin-wide sediment loading at Mount Pinatubo, Philippines: Geological Society of America Bulletin, 117 (1–2), 195–211.         [ Links ]

Haddad, B., Pastor, M., Palacios, D., Muñoz-Salinas, E., 2011, A SPH Depth Integrated Model for Popocatépetl 2001 Lahar (Mexico): Sensitivity Analysis and Runout Simulation: Engineering Geology, 114(3–4), 312–329.         [ Links ]

Hamidi, S., 1989, Lahar of Galunggung Volcano from 1982 through 1986 (resumen), en Proceedings of International Symposium on Erosion and Volcanic Debris Flow Technology, Yogyakarta, Indonesia, Ministry of Public Works, VP1-1–VP1-23.         [ Links ]

Hayes, S.K., Montgomery, D.R., Newhall, C.G., 2002, Fluvial sediment transport and deposition following the 1991 eruption of Mt. Pinatubo: Geomorphology, 45, 211–224.         [ Links ]

Hirao, K., Yoshida, M., 1989. Sediment yield of Mt Galunggung after eruption in 1982 (resumen), en Proceedings of International Symposium on Erosion and Volcanic Debris Flow Technology, Yogyakarta, Indonesia, Ministry of Public Works, V21-1–V21-22.         [ Links ]

Hodgson, K., Manville, V., 1999, Sedimentology and flow behavior of a rain-triggered lahar, Mangatoetoenui Stream, Ruapehu volcano, New Zealand: Geological Society of America Bulletin, 111, 743–754.         [ Links ]

Inbar, M., Lugo Hubp, J., Villers Ruiz, L., 1994, The geomorphological evolution of the Paricutin cone and lava flows, Mexico, 1943– 1990: Geomorphology, 9(1), 57– 76.         [ Links ]

Inbar, M., Enriquez, A. R., Granel, J. H, 2001, Morphological changes and erosion processes following 1982 eruption of El Chichón volcano, Chiapas, México: Geomophologie, 3, 175–184.         [ Links ]

Instituto Nacional de Estadística, Geografía e Informática (INEGI), 1978, Cartografías digitales de las hojas Amecameca de Juárez E14-B41 y Huejotzingo E14-B42, escala 1:50,000 (en línea): México. D.F., Secretaría de Programación y Presupuesto, Instituto Nacional de Estadística, Geografía e Informática, disponible en <http://www3.inegi.org.mx/sistemas/descarga/>, consultado 5 de junio de 2013.         [ Links ]

Instituto Nacional de Estadística y Geografía (INEGI), 2011a, Principales resultados del Censo de Población y Vivienda 2010. Puebla (en línea), disponible en <http://www3.inegi.org.mx/sistemas/productos/>, consultado 5 de junio de 2013.         [ Links ]

Instituto Nacional de Estadística y Geografía (INEGI), 2011b, Principales resultados del Censo de Población y Vivienda 2010. Distrito Federal, México (en línea), disponible en <http://www3.inegi.org.mx/sistemas/productos/>, consultado 5 de junio de 2013.         [ Links ]

Iwamoto, M., 1996, Prevention of disasters caused by debris flows at Unzen volcano, Japan, en Slaymaker, O. (ed.), Geomorphic hazards: Chichester, John Wiley and Sons, 95–110.         [ Links ]

Janda, R.J., Scott, K.M., Martinson, H.A., 1981, Lahar movement, effects, and deposits, en Lipman, P.W., Mullineaux, D.R. (eds.), The 1980 eruptions of Mount St Helens, Washington, U.S. Geological Survey Professional Paper, 1250, 461–478.         [ Links ]

Janda, R.J., Meyer, D.F., Childers, D., 1984, Sedimentation and geomorphic changes during and following the 1980–1983 eruptions of Mount St. Helens, Washington: Shin-Sabo, 37(2), 10–21.         [ Links ]

Kadomura, H., Imagawa, T., Yamamoto, H., 1983, Eruption-induced rapid erosion and mass movements on Usu Volcano, Hokkaido: Zeitschrift für Geomorphologie N.F. Supplementband, 46, 123–142.         [ Links ]

Lauer, W., Stiehl, E., 1973, La clasificación del clima en la región de Puebla-Tlaxcala: Comunicaciones, 7, 31–36.         [ Links ]

Lavigne, F., 1998, Les lahars du volcan Merapi, Java central, Indonésie: déclenchement, budget sédimentaire, dynamique et zonage des risques associés. Clermont-Ferrand, Université Blaise-Pascal, tesis doctoral, 539 p.         [ Links ]

Lavigne, F., Thouret, J.C., 2002, Sediment transportation and deposition by rain-triggered Lahar at Merapi Volcano, Central Java, Indonesia: Geomorphology, 49, 45–69.         [ Links ]

Leavesley, G.H., Lusby, G.C., Lichty, R.W., 1989, Infiltration and erosion characteristics of selected tephra deposits from the 1980 eruption of Mount St. Helens, Washington, USA: Hydrological Sciences Journal, 34(3), 339–353.         [ Links ]

Long, W.S., 2000, Development of digital terrain representation for use in river modelling, en Maidment, D., Djokic, D. (eds.), Hydrologic and hydraulic modeling support with Geographic Information System: Redlands, California, ESRI Press, 145–154.         [ Links ]

Macías, J.L., 2005, Geología e historia eruptiva de algunos de los grandes volcanes activos de México: Boletín de la Sociedad Geológica Mexicana, Volumen Conmemorativo del Centenario Temas Selectos de la Geología Mexicana, LVII, 3, 379–424.         [ Links ]

Major, J.J., 2003, Post-eruption hydrology and sediment transport in volcanic river systems: Water Resources Impact, 5 (3), 11–15.         [ Links ]

Major, J.J., 2004, Posteruption suspended sediment transport at Mount St. Helens: Decadal-scale relationships with landscape adjustments and river discharges: Journal of Geophysical Research: Earth Surface, 109, F01002, doi:10.1029/2002JF000010.         [ Links ]

Major, J.J., Janda, R.J., Daag, A.S., 1996, Watershed disturbance and lahars on the east side of Mount Pinatubo during the mid-June 1991 eruptions, en Newhall, C.G., Punongbayan, R.S. (eds.), Fire and mud: Eruptions and lahars of Mount Pinatubo, Philippines: Quezon City and Seattle, Philippine Institute of Volcanology and Seismology and University of Washington Press, 895–919.         [ Links ]

Major, J.J., Pierson, T.C., Dinehart, R.L., Costa, J.E., 2000, Sediment yield following severe volcanic disturbance- A two-decade perspective from Mount St. Helens: Geology, 28(9), 819–822.         [ Links ]

Manville, V., Newton, E.H., White, J.D.L., 2005, Fluvial responses to volcanism: resedimentation of the 1800a Taupo ignimbrite eruption in the Rangitaiki River catchment, North Island, New Zealand: Geomorphology, 65, 49–70.         [ Links ]

Meyer, D.F., Martinson, H.A., 1989, Rates and Processes of Channel Development and Recovery Following the 1980 Eruption of Mount St. Helens, Washington: Hydrological Sciences Journal, 34, 115–127.         [ Links ]

Mizuyama, T., Kobashi, S., 1996, Sediment yield and topographic change after major volcanic activity, en Walling, D.E., Webb, B.W. (eds.), Erosion and Sediment Yield: Global and Regional Perspectives: Reino Unido, International Association of Hydrological Sciences (Exeter Symposium), 295–301.         [ Links ]

Muñoz-Salinas, E., 2007, Los lahares del Popocatépetl: obtención y tratamiento de la información para la prevención de riesgos, Madrid, Universidad Complutense de Madrid, tesis doctoral, 229 p.         [ Links ]

Muñoz-Salinas, E., Manea, V.C., Palacios, D., Castillo-Rodríguez, M., 2007, Estimation of lahar flow velocity on Popocatépetl volcano (Mexico): Geomorphology, 92, 91–99.         [ Links ]

Muñoz-Salinas, E., Castillo-Rodríguez, M., Manea, V., Manea, M., Palacios, D., 2009, Lahar flow simulations using LAHARZ program: Application for the Popocatépetl volcano, Mexico: Journal of Volcanology and Geothermal Research, 182(1-2), 13–22.         [ Links ]

Newhall, C.G., Punongbayan, R.S., 1996, Fire and mud: Eruptions and lahars of Mount Pinatubo, Philippines: Quezon City and Seattle, Philippine Institute of Volcanology and Seismology and University of Washington Press, 1126 p.         [ Links ]

North American Drought Monitor (NADM): North American Drought Monitor Maps (en línea), disponible en <http://www.ncdc.noaa.gov/temp-and-precip/drought/nadm//>, consultado 8 de agosto de 2009.         [ Links ]

Palacios, D., 1995, Rockslide processes at the North Slope of Popocatepetl Volcano: Permafrost and Periglacial Processes, 6(4), 345–359.         [ Links ]

Palacios, D., 1996, Recent Geomorphologic evolution of a glaciovolcanic active stratovolcano: Popocatepetl (Mexico): Geomorphology, 16, 319–335.         [ Links ]

Palacios, D., Zamorano, J.J., Parrilla, G., 1998, Proglacial debris flows in Popocatépetl north face and their relation to 1995 eruption: Zeitschrift Geomorphologie, 42(3), 273-295.         [ Links ]

Pierson, T.C., 1985, Initiation and flow behaviour of the 1980 Pine Creek and Muddy River lahars, Mount St. Helens, Washington: Geological Society of American Bulletin, 96, 1056–1069.         [ Links ]

Pierson, T.C., Scott, K.M., 1985, Downstream Dilution of a lahar: Transition from debris flow to hyperconcentrated streamflow: Water Resources Research, 21(10), 1511–1524.         [ Links ]

Pierson, T.C., Costa, K.M., 1987, A rheologic classification of subaerial sediment-water flows, en Costa, J.E., Wieczorek, G.E. (eds.), Debris flows/Avalanches: Process, recognition, and mitigation, Boulder, Colorado, Geological Society of America, Reviews in Engineering Geology, 8, 1–12.         [ Links ]

Pierson, T.C., Janda, R.J., Umbal, J.V., Daag, A.S., 1992, Immediate and long-term hazards from lahars and excess sedimentation in rivers draining Mt. Pinatubo, Philippines: Vancouver, Washington, U.S. Geological Survey Water-Resources Investigations Report, 92-4039, 35 p.         [ Links ]

Pierson, T.C., Daag, A.S., Reyes, P.J., Regalado, M.T.M., Solidum, R.U., Tubianosa, B.S., 1996, Flow and Deposition of Hot lahars on the East Side of Mt. Pinatubo, July-October 1991, en Newhall, C.G., Punongbayan, R.S. (eds), Fire and Mud, Eruptions and Lahars of Mt. Pinatubo, Philippines: Quezon City and Seattle, Philippine Institute of Volcanology and Seismology and University of Washington Press, 921–950.         [ Links ]

Renschler, C.S., 2005, Scales and uncertainties in using models and GIS for volcano hazard prediction: Journal of Volcanology and Geothermal Research, 139, 73–87.         [ Links ]

Rodolfo, K.S., 1989, Origin and Early Evolution of Lahar Channel at Mabinit, Mayon Volcano, Philippines: Geological Society of America Bulletin, 101, 414–426.         [ Links ]

Rodolfo, K.S., Arguden, A.T., 1991, Rain-lahar generation and sediment-delivery systems at Mayon Volcano, Philippines, en Fisher, R.V., Smith, G.A. (eds.), Sedimentation in Volcanic Settings: Tulsa, Oklahoma, SEPM Special Publication, 45,. 71–87.         [ Links ]

Rodolfo, K.S., Arguden, A.T., Solidum, R.U., Umbal, J.V., 1989, Anatomy and behaviour of a post-eruptive rain lahar triggered by a typhoon on Mayon Volcano, Philippines: Bulletin of the International Association of Engineering Geology, 40, 55–66.         [ Links ]

Scott, K.M., Janda, R.J., Cruz, E.G., Gabinete, E., Eto, I., Isada, M., Sexton, M., Hadley, K., 1996, Channel and sedimentation responses to large volumes of 1991 volcanic deposits on the east flank of Mount Pinatubo, en Newhall, C.G., Punongbayan, R.S. (eds.), Fire and mud: Eruptions and lahars of Mount Pinatubo, Philippines: Quezon City and Seattle, Philippine Institute of Volcanology and Seismology and University of Washington Press, 971–988.         [ Links ]

Scott, K. M., Vallance, J.W., Kerle, N., Macías, J.L., Strauch, W., Devoli, G., 2005, Catastrophic precipitation-triggered lahar at Casita volcano, Nicaragua: occurrence, bulking and transformation: Earth Surface Processes and Landforms, 30, 59–79.         [ Links ]

Segerstrom, K., 1950, Erosion studies at Paricutín, State of Michoacán, Mexico: U.S. Geological Survey Bulletin, 965-A, 1–64.         [ Links ]

Sheridan, M.F., Hubbard, B., Bursik, M.I., Abrams, M., Siebe, C., Macías, J.L., Delgado, H., 2001, Gauging short-term volcanic hazards at Popocatépetl: EOS, Transactions American Geophysical Union, 82, 187–188.         [ Links ]

Shimokawa, E., Taniguchi, Y., 1983, Sediment yield from hillside slope of active volcanoes: Proceedings of the Symposium on erosion control in volcanic areas, Sabo Division, Erosion Control Department, Technical Memorandum, 1908, 155–181.         [ Links ]

Shimokawa, E., Jitousono, T., Yazawa, A., Kawagoe, R., 1989, An effect of tephra cover on erosion processes of hillslopes in and around Sakurajima Volcano (resumen), en Proceedings of International Symposium on Erosion and Volcanic Debris Flow Technology, Yogyakarta, Indonesia, Ministry of Public Works, V32-1–V32-21.         [ Links ]

Siebe, C., Abrams, M., Macías, J.L., Obenholzner, J., 1996, Repeated volcanic disasters in Prehispanic time at Popocatépetl, central Mexico: Past key to the future?: Geology, 24, 399–402.         [ Links ]

Simon, A., 1999, Channel and drainage-basin response of the Toutle River system in the aftermath of the 1980 eruption of Mount St. Helens, Washington: U.S. Geological Survey Open-File Report, 96-633, 130 p.         [ Links ]

Smith, G.A., Fritz, W.J., 1989, Volcanic Influences on Terrestrial Sedimentation: Geology, 17, 375–376.         [ Links ]

Smith, G.A., Lowe, D.R., 1991, Lahars: Volcano-Hydrologic Events and Deposition in the Debris Flow-Hyperconcentrated Flow Continuum: Sedimentation in Volcanic Settings: SEPM Special Publication, 45, 59–69.         [ Links ]

Tanarro, L.M., Andrés, N., Zamorano, J.J., Palacios, D., Renschler, C.S., 2010, Geomorphological evolution of a fluvial channel after primary lahar deposition: Huiloac Gorge, Popocatépetl volcano (Mexico): Geomorphology, 122, 178–190.         [ Links ]

Thouret, J.C., Lavigne, F., 2000, Lahars: occurrence, deposits and behaviour of volcano-hydrologic flows, en Leyrit, H., Montenat, C. (eds.), Volcaniclastic rocks from magmas to sediments: Amsterdam, Gordon and Breach Science Publishers, 151–174.         [ Links ]

Tuñgol, N.M., 2002, Lahar initiation and sediment yield in the Pasig-Potrero River basin, Mount Pinatubo, Philippines: Canterbury, New Zealand, University of Canterbury, tesis doctoral, 172 p.         [ Links ]

Tuñgol, N.M., Regalado, T.S., 1997, Rainfall, Acoustic Flow Monitor records and observed lahars of the Sacobia river in 1992, en Newhall, C.G., Punongbayan, R.S. (eds.), Fire and mud: Eruptions and lahars of Mount Pinatubo, Philippines: Quezon City and Seattle, Philippine Institute of Volcanology and Seismology and University of Washington Press, 1023–1032.         [ Links ]

Umbal, J.V., 1997, Five years of lahars at Pinatubo volcano: Declining but still potentially lethal hazards: Journal of the Geological Society of the Philippines, 52(1). 1–19.         [ Links ]

Vallance, J. W., 2000, Lahars, en Sigurdsson, H. (ed.) Encyclopedia of Volcanoes: San Diego, California, Academia Press, 601–616.         [ Links ]

van Westen, C.J., 1997, Modelling erosion from pyroclastic flow deposits on Mount Pinatubo, en van Westen, C.J., Saldaña, A., Uría, P., Chávez, G., (eds.), ILWIS Applications Guide: Enschede, The Netherlands, 53–72.         [ Links ]

van Westen, C.J. y Daag, A., 2005, Analysing the relation between rainfall characteristics and lahar activity at Mount Pinatubo, Philippines: Earth Surface Processes and Landforms, 30. 1663–1674.         [ Links ]

Waldron, H.H., 1967, Debris flow and erosion control problems caused by the ash eruptions of Irazu Volcano, Costa Rica: U.S. Geological Survey Bulletin, 1241(1), 1–37.         [ Links ]

Yamamoto, H., 1984, Erosion of the 1977–1978 tephra layers on a slope of Usu Volcano, Hokkaido: EOS, Transactions American Geophysical Union, 5 (2). 111–124.         [ Links ]