Combining slope stability and groundwater flow models to assess stratovolcano collapse hazard

Flank collapses are a well-documented hazard at volcanoes. Elevated pore-fluid pressures and hydrothermal alteration are invoked as potential causes for the instability in many of these collapses. Because pore pressure is linked to water saturation and permeability of volcanic deposits, hydrothermal alteration is often suggested as a means of creating low-permeability zones in volcanoes. Here, we seek to address the question: What alteration geometries will produce elevated pore pressures in a stratovolcano, and what are the effects of these elevated pressures on slope stability? We initially use a finite element groundwater flow model (a modified version of OpenGeoSys) to simulate `generic' stratovolcano geometries that produce elevated pore pressures. We then input these results into the USGS slope-stability code Scoops3D to investigate the effects of alteration and magmatic intrusion on potential flank failure. This approach integrates geophysical data about subsurface alteration, water saturation and rock mechanical properties with data about precipitation and heat influx at Cascade stratovolcanoes. Our simulations show that it is possible to maintain high-elevation water tables in stratovolcanoes given specific ranges of edifice permeability (ideally between 10-15 and 10-16 m2). Low-permeability layers (10-17 m2, representing altered pyroclastic deposits or altered breccias) in the volcanoes can localize saturated regions close to the surface, but they may actually reduce saturation, pore pressures, and water table levels in the core of the volcano. These conditions produce universally lower factor-of-safety (F) values than at an equivalent dry edifice with the same material properties (lower values of F indicate a higher likelihood of collapse). When magmatic intrusions into the base of the cone are added, near-surface pore pressures increase and F decreases exponentially with time ( 7-8% in the first year). However, while near-surface impermeable layers create elevated water tables and pore pressures, they do not necessarily produce the largest or deepest collapses. This suggests that mechanical properties of both the edifice and layers still exert a significant control, and collapse volumes depend on a complex interplay of mechanical factors and layering.