Slope stability analysis of the Pacaya Volcano, Guatemala, using limit equilibrium and finite element method

The Pacaya volcanic complex is part of the Central American volcanic arc, which is associated with the subduction of the Cocos tectonic plate under the Caribbean plate. Located 30 km south of Guatemala City, Pacaya is situated on the southern rim of the Amatitlan Caldera. It is the largest post-caldera volcano, and has been one of Central America’s most active volcanoes over the last 500 years. Between 400 and 2000 years B.P, the Pacaya volcano had experienced a huge collapse, which resulted in the formation of horseshoe-shaped scarp that is still visible. In the recent years, several smaller collapses have been associated with the activity of the volcano (in 1961 and 2010) affecting its northwestern flanks, which are likely to be induced by the local and regional stress changes. The similar orientation of dry and volcanic fissures and the distribution of new vents would likely explain the reactivation of the pre-existing stress configuration responsible for the old-collapse. This paper presents the first stability analysis of the Pacaya volcanic flank. The inputs for the geological and geotechnical models were defined based on the stratigraphical, lithological, structural data, and material properties obtained from field survey and lab tests. According to the mechanical characteristics, three lithotechnical units were defined: Lava, Lava-Breccia and Breccia-Lava. The Hoek and Brown’s failure criterion was applied for each lithotechnical unit and the rock mass friction angle, apparent cohesion, and strength and deformation characteristics were computed in a specified stress range. Further, the stability of the volcano was evaluated by two-dimensional analysis performed by Limit Equilibrium (LEM, ROCSCIENCE) and Finite Element Method (FEM, PHASE 2 7.0). The stability analysis mainly focused on the modern Pacaya volcano built inside the collapse amphitheatre of “Old Pacaya”. The volcanic instability was assessed based on the variability of safety factor using deterministic, sensitivity, and probabilistic analysis considering the gravitational instability and the effects of external forces such as magma pressure and seismicity as potential triggering mechanisms of lateral collapse. The preliminary results from the analysis provide two insights: first, the least stable sector is on the south-western flank of the volcano; second, the lowest safety factor value suggests that the edifice is stable under gravity alone, and the external triggering mechanism can represent a likely destabilizing factor.

SLOPE STABILITY ANALYSIS OF MOUNT MEAGER, SOUTH-WESTERN BRITISH COLUMBIA, CANADA

The Mount Meager Volcanic Complex (MMVC) in south-western British Columbia is a potentially active, hydrothermally altered massif comprising a series of steep, glaciated peaks. Climatic conditions and glacial retreat has led to the further weathering, exposure and de-buttressing of steep slopes composed of weak, unconsolidated material. This has resulted in an increased frequency of landslide events over the past few decades, many of which have dammed the rivers bordering the Complex. The breach of these debris dams presents a risk of flooding to the downstream communities. Preliminary mapping showed there are numerous sites around the Complex where future failure could occur. Some of these areas are currently undergoing progressive slope movement and display features to support this such as anti-scarps and tension cracks. The effect of water infiltration on stability was modelled using the Rocscience program Slide 6.0. The main site of focus was Mount Meager in the south- east of the Complex where the most recent landslide took place. Two profiles through Mount Meager were analysed along with one other location in the northern section of the MMVC, where instability had been detected. The lowest Factor of Safety (FOS) for each profile was displayed and an estimate of the volume which could be generated was deduced. A hazard map showing the inundation zones for various volumes of debris flows was created from simulations using LAHARZ. Results showed the massif is unstable, even before infiltration. Varying the amount of infiltration appears to have no significant impact on the FOS annually implying that small changes of any kind could also trigger failure. Further modelling could be done to assess the impact of infiltration over shorter time scales. The Slide models show the volume of material that could be delivered to the Lillooet River Valley to be of the order of 109 m3 which, based on the LAHARZ simulations, would completely inundate the valley and communities downstream. A major hazard of this is that the removal of such a large amount of material has the potential to trigger an explosive eruption of the geothermal system and renew volcanic activity. Although events of this size are infrequent, there is a significant risk to the communities downstream of the complex.