Could Glacial Retreat-Related Landslides Trigger Volcanic Eruptions? Insights from Mount Meager,

Gioachino Roberti, Brent Ward, Benjamin van Wyk de Vries, Nicolas Le Corvec, Swetha Venugopal, Glyn Williams-Jones, John J. Clague, Pierre Friele, Giacomo Falorni, Geidy Baldeon, Luigi Perotti, Marco Giardino, and Brian Menounos

Abstract chamber at 3–16 km depth. Based on numerical model simulations carried out to constrain the stress change, the Mount Meager, a glacier-clad volcanic complex in British failure would affect the stress field to depths of up to Columbia, , is known for its large landslides, as *6 km, with changes in effective stress of up to well as a major eruption about 2360 years ago. In 2010, *4 MPa. The change in effective stress following such after decades of glacier retreat, the south flank of Mount a landslide might destabilize the magmatic chamber and Meager collapsed, generating a huge (53 Mm3) landslide. trigger an eruption. This result also suggests that a In 2016, fumaroles formed ice caves in one of the glaciers previously documented major flank collapse may have on the complex. This glacier is bordered by a large had a role in the 2360 cal yr BP eruption. unstable slope presently moving about 3.5 cm per month. If this slope were to fail, a long-runout debris avalanche Keywords would reach the floor of the Lillooet River valley, with possible destructive effects on downstream infrastructure. Volcanic landslide Á Stress changes Á Eruption trigger Á The unloading of the volcanic edifice from an abrupt FEM Á InSAR failure would also have unknown effects on the magmatic plumbing system. From geochemical, geophysical, and petrological data, we infer the presence of a magmatic Introduction G. Roberti (&) Minerva Intelligence, 850 W Hastings St #301, Vancouver, BC, Mount Meager is a glaciated volcanic complex in the Canada e-mail: [email protected] southern Coast Mountains of British Columbia, with activity spanning the past *2 Ma (Read 1979) (Fig. 1). The most B. Ward Á S. Venugopal Á G. Williams-Jones Á J. J. Clague Earth Sciences Department, Simon Fraser University, Burnaby, recent eruption occurred 2360 years ago (Clague et al. 1995) BC, Canada and produced a tephra plume that tracked eastward into B. van Wyk de Vries Á S. Venugopal western Alberta (Hickson et al. 1999). Today, Mount Mea- CNRS, IRD, OPGC, Laboratoire Magmas et Volcans, Université ger hosts an active hydrothermal system and is the site of Clermont Auvergne, Aubiere Cedex, France frequent shallow earthquakes (Friele et al. 2008; Venugopal N. Le Corvec et al. 2017). It is known for its very large pre-historic and 111 rue Jules Ferry, 33200 Bordeaux, France historic landslides (Friele et al. 2008), of which a 53 Mm3 P. Friele debris avalanche in 2010 is the largest landslide in Canadian Cordilleran Geoscience, Squamish, BC, Canada history (Guthrie et al. 2012). Following the 2010 landslide, G. Falorni Á G. Baldeon recent studies have focused on the ongoing instability on TRE-ALTAMIRA Inc., Vancouver, BC, Canada Mount Meager and its link to recent deglaciation and climate L. Perotti Á M. Giardino warming (Roberti et al. 2018a, b). These studies have Earth Sciences Department, GeoSitLab - GIS and Geomatics identified many large unstable slopes, one of which, on the Laboratory, University of Torino, Torino, Italy west flank of Plinth Peak (Fig. 1), is of particular concern as B. Menounos it moved 3 cm during a 24-day period during the summer of Geography Program and Natural Resources and Environmental 2016 (determined by InSAR interferometry; see Roberti Studies Institute, University of Northern British Columbia, Prince George, BC, Canada

© Springer Nature Switzerland AG 2021 147 V. Vilímek et al. (eds.), Understanding and Reducing Landslide Disaster Risk, ICL Contribution to Landslide Disaster Risk Reduction, https://doi.org/10.1007/978-3-030-60319-9_15 148 G. Roberti et al. et al. 2018b). This slope is near an amphitheatre (A in collapse in the 2360 cal yr BP eruption. That collapse might Fig. 1) that hosts a retreating glacier and a fumarole field. have decompressed the volcanic plumbing system and In this study, we present a numerical model to investigate contributed to the eruption. the stress field evolution prior to and after a hypothetical collapse of the west flank of Plinth Peak. We discuss how this stress change could affect the hydrothermal and volcanic Methods plumbing systems and lead to an eruption. We also docu- ment another large edifice collapse from Job Valley, which We describe the geology and structure of the west flank of happened a few centuries before the 2360 cal yr BP erup- Plinth Peak and Job Valley from historical aerial pho- tion. This landslide originated from the above-mentioned tographs, Lidar data, literature review and field work. Field amphitheatre near Plinth Peak that is now occupied by Job work at Plinth, Job Glacier, and the mouth of Job Valley was Glacier and the fumarole field. Discussion of decompression conducted in October 2017. We collected a core sample on mechanisms after the eventual collapse of Plinth Peak and the landslide deposit at the mouth of Job Valley to constrain documentation of the collapse in Job Valley enable an the temporal relation between the landslide and the examination of the possible role of the earlier edifice 2360 cal yr BP eruption. We used a 2D elastic finite element model (FEM) with COMSOL Multiphysics® to explore the effect of the collapse of Plinth Peak on the state of stress within the Mount Meager volcanic edifice. The FEM topography was built by superimposing E-W profiles of the west flank of Plinth (B-B’ in Fig. 1) and of Job Valley 500 m to the south, assuming that the failure would remove the bulge of the west flank of Plinth Peak to reach the profile of Job Valley. Two elastic domains were defined to represent both volcanic material (Plinth Peak) and the granitic basement rocks. We used Plinth Peak’s 2360 cal yr BP eruption dacite to model the volcanic edifice and the Fall Creek Stock monzogranite to model the substratum rock (Campbell et al. 2016; Kelly Russell, personal communication 2018). Fol- lowing previous work (Le Corvec et al. 2018), the FEM defined as a half-space, is loaded with gravitational forces. Except from its surface that is set free to move, the other boundaries are defined by a roller condition restricting dis- placements normal to the boundary. COMSOL allows to compute the differential stresses and the orientation of the stress field before and after the removal of the landslide material, subsequently allowing to compute the difference between the two cases.

Results

Plinth

A prominent knob forms the west flank of Plinth Peak on the east side of Job Valley about 2 km from the Job Glacier headwall. The knob has an area of 2.7 Â 106 m2 and rises 1000 m from Job Glacier to the ridge crest (Fig. 2). DInSAR analysis of Sentinel-1 satellite images show displacement of Fig. 1 Job Valley, Plinth Peak, and the 2360 cal yr BP eruption crater. this part of the slope of 34 ± 10 mm over a 24-day period Scars, unstable slopes, landslides, eruption deposits, fumaroles, sample, fi fi from July to August 2016 (Roberti et al. 2018a) and and pro le locations are indicated. a Job Valley amphitheatre. b Pro le ± shown in Fig. 2. c 2010 landslide scar. d Job landslide deposit. Base 30 10 mm over a 12-day period in July 2017. The upper map is a hillshaded rendering of Lidar data from 2015–2016 slope extends 400 m above the lower slope (Fig. 2) and has Could Glacial Retreat-Related Landslides Trigger … 149 near-vertical cliffs and buttresses composed of horizontally the 2360 cal yr BP eruption. This debris flow unit is thus layered rhyodacitic lava flows of the Plinth assemblage interpreted as a close precursor to the eruption. (p9f, Read 1979). SqueeSAR analysis of the 1992–2000 ERS data indicate LOS (line of sight) displacements up to 17 mm yr−1. Analysis of the 2015–2018 Sentinel-1 data Discussion indicate LOS displacements of up to 10 mm yr−1. 2D stress field modeling in COMSOL shows that removal of the Plinth bulging rock mass would affect the effective stress to a depth *6kmbyupto *4 MPa (Fig. 2). Most of the decom- In view of displacement rates, geometry, geological setting, pression would occur at the toe of the slope, about 1 km and a low (1.084–1.077) factor of safety (FOS) calculated by downstream of the fumaroles. (Hetherington 2014), the west flank of Plinth Peak can be considered to be at failure and near catastrophic collapse threshold. The fumaroles on Job Glacier may accelerate Job glacier melting and bring the slope closer to catastrophic failure. The base of the slope, weakened by glacier erosion A 2-km long, >100 m thick, hummocky landslide deposit is and having been recently deglaciated, shows high defor- located at the Job Creek—Lillooet River confluence (D in mation rates. It is likely to fail first, followed quickly by Fig. 1). Considering the thickness and extent of the deposit, retrogressive failure involving the whole slope, similar to the we tentatively relate the landslide to the formation of the evolution of the 2010 collapse (Roberti et al. 2018b). The large amphitheatre at the head of Job Valley (*109 m3)(A failed rock mass would rapidly fragment, liquefy, and in Fig. 1). This amphitheatre is the only feature large enough transform into a debris flow, as observed at other large to be the source of the landslide deposit and furthermore has landslides at Mount Meager (e.g. Roberti et al. 2017) and in widespread hydrothermal alteration, consistent with the other volcanic terrain (e. g. Scott et al. 2001) with a volume altered nature of the deposit. Given the size of this deposit, of 108–109 m3 (Class 9 debris flow; Jakob 2005). The this landslide may be related to one of the distal volcanic hydropower infrastructure near the Keyhole Falls would be debris flow and hyperconcentrated flow in the valley fill. In impacted, as it is within the reach of a 107 m3 debris flow. four cores 50 km downstream from Mount Meager there is a Even populated areas in Pemberton Meadows, which are massive debris flow unit, about 4 m thick forming a 2 km within the reach of a 109 m3 debris flow, might be affected valley-wide sheet, with an age bracketed between (see Friele and Clague 2004; Friele et al. 2008; Simpson 2570 ± 40 14C yr BP and 2690 ± 50 14C yr BP, or et al. 2006; Hetherington 2014 for inundation zones of between 2540–2970 cal yr old, of unknown proximal debris flows from Mount Meager). deposit correspondence. This event predates the 2360 cal yr Decompression following the collapse of Plinth Peak BP eruption by 100–200 years. It is capped by a thin humic would influence the stress field to a depth of *6 km. This laminae, the incipient paleosol, carbon source for the sudden decompression would affect the near-surface younger bracketing age, and in turn is capped by a reversely hydrothermal system at Job Glacier, as well as the deeper graded, pumiceous hyperconcentrated flow unit related to magmatic plumbing system. The geothermal gradient in the

Fig. 2 Differential stresses (Dr) and stress field orientations before to *4 MPa and would reach a depth of *6 km. The blue and red and after a hypothetical collapse of the knob west of Plinth Peak. The cones represent the direction of the maximum and minimum compres- oblique line separates volcanic rocks to the left from basement rocks to sional stresses (r1 and r3), respectively the right. The change in effective stress after the collapse would be up 150 G. Roberti et al. area is 80–100 °Ckm-1, with recorded temperatures in impact infrastructure and communities near the volcano. The geothermal explorative drill holes of up to 270 °C (Ghom- model of stress change and a new field description of a large shei et al. 2013). Shallow earthquakes (1–6.7 km depth; rock avalanche source in upper Job Valley suggest that the Friele et al. 2008) occur below Mount Meager. This seismic past landslide may have had a role in the 2360 cal yr BP activity and the high geothermal gradient may indicate the Mount Meager eruption. The large (*109 m3) collapse that presence of magma at shallow depths (3–10 km). The stress created the amphitheatre at the head of Job Valley preceded perturbation of the system by a large flank collapse from the eruption by, at most, a few centuries. Decompression of Plinth Peak might trigger a series of processes, including the magmatic system might have allowed new magma to rise fracture and magma propagation, gas exolution, and magma and mix, leading to the 2360 cal yr BP eruption. The trigger differentiation and mixing that could lead to an eruption. for the collapse is still unknown but may be related to the Deeper magma (5–16 km, 1057–1142 ˚C, and 1.4–4.8 kbar; hydrothermal activity and continued glacial erosion. Venugopal et al. 2020) is also present below Mount Meager, as indicated by olivine melt inclusion analysis of basalts of Acknowledgements Alex Wilson and Kelly Russell provided the data the Mosaic assemblage (p10f, Read 1979). The decom- for Mount Meager rock properties used in the model. We thank Hazel pression following the landslide would not affect these Wong, Brad Woods, Erik Young, and Kai Friele for the help in the field. Financial support for the 2015 Lidar data acquisition was pro- deeper magma reservoirs unless there were a telescoping of vided by the Natural Sciences Engineering Research Council of Canada events downward through the system, and the localised through Discovery Grants to Clague and Menounos and a grant from decompression would not likely induce partial melting in the the Canada Foundation for Innovation (Menounos). Natural Resources mantle. Canada (NRCan), through Melanie Kelman, provided funds for acquisition of the 2016 Lidar dataset. Fieldwork was supported by SFU Earth Science Department Chair grants to Ward. The InSAR analyses were provided by TRE-Altamira. Financial support for Roberti's PhD Job salary was partially provided by the “End of an Arc: The Remarkable Life and Death of a Volcanic Arc” project, a French-Canadian partnership. A very large landslide such as this may have perturbed the stress field in the volcanic edifice, and may have triggered the ascent of deeper, less differentiated magma and its References mixing with a shallow, more differentiated magma, trigger- ing the 2360 cal yr BP eruption (Stasiuk et al. 1996). The Campbell ME, Porritt L, Russell JK (2016) Forensic recovery of interval between the landslide and the eruption could be transient eruption parameters for the 2360 BP fall deposit, Mount related to the time that the magmatic system took to react to Meager, British Columbia. J Volcanol Geotherm Res. 312:11–25. the decompression. The location of the vent, NE of the Job https://doi.org/10.1016/j.jvolgeores.2016.01.013 amphitheatre collapse scarp, might also have been condi- Clague JJ, Evans SG, Rampton VN, Woodsworth GJ (1995) Improved fi age estimates for the White River and Bridge River tephras, western tioned by the large collapse. The stress eld perturbation Canada. Can J Earth Sci 32:1172–1179. https://doi.org/10.1139/ after the Job Valley collapse could have conditioned the e95-096 location of the vent, lateral to the collapse scar, and triggered Delcamp A, Troll VR, van Wyk de Vries B, et al (2012) Dykes and a sequence of events that led to the eruption as observed at structures of the NE rift of Tenerife, Canary Islands: A record of stabilisation and destabilisation of ocean island rift zones. Bull other volcanoes (Delcamp et al. 2012; Zurek et al. 2015; Volcanol 74:963–980. https://doi.org/10.1007/s00445-012-0577-1 Maccaferri et al. 2017). Friele PA, Clague JJ (2004) Large holocene landslides from pylon peak. Southwestern British Columbia. Can J Earth Sci 182:165– 182. https://doi.org/10.1139/E03-089 Friele PA, Jakob M, Clague JJ (2008) Hazard and risk from large Conclusions landslides from Mount Meager volcano, British Columbia, Canada. Georisk Assess Manag Risk Eng Syst Geohazards 2:48–64. https:// The high gravitational slope deformation rates and rapid doi.org/10.1080/17499510801958711 glacier retreat near Plinth Peak are factors favouring large Ghomshei M, Sadlier-Brown T, Hassani F (2013) Current status of the fi Pebble Creek (North Meager) Geothermal Project, 11–15 August catastrophic slope failures. We applied a 2D nite element 2013, Montreal, QC. Suthwestern British Columbia. 23rd World model to assess the state of stress after a catastrophic Plinth Min Congr, 10p slope failure. Changes in effective stresses would reach a Guthrie RH, Friele P, Allstadt K et al (2012) The 6 August 2010 Mount depth of up to *6kmwithupto*4 MPa decompression. Meager rock slide-debris flow Coast Mountains, British Columbia: fi Characteristics, Dynamics, and Implications for Hazard and Risk Such changes in con ning pressure could trigger hydrother- Assessment. Nat Hazards Earth Syst Sci 12:1277–1294. https://doi. mal explosions from the shallow hydrothermal system at Job org/10.5194/nhess-12-1277-2012 Glacier. Fracture propagation, magma migration, and gas Hetherington R (2014) Slope stability analysis of mount meager, exolution following decompression could also lead to a South-Western British Columbia. MSc Thesis, Michigan Technol Univ 68 magmatic eruption. Such a landslide and eruption would Could Glacial Retreat-Related Landslides Trigger … 151

Hickson CJ, Russell JK, Stasiuk MV (1999) Volcanology of the 2350 Columbia. Landslides 15:637–647. https://doi.org/10.1007/ B.P. eruption of Mount Meager volcanic complex, British s10346-017-0901-0 Columbia, Canada: Implications for hazards from eruptions in Scott KM, Macias JL, Naranjo JA et al (2001) Catastrophic debris topographically complex terrain. Bull Volcanol 60:489–507. https:// flows transformed from landslides in volcanic terrains: mobility, doi.org/10.1007/s004450050247 hazard assessment, and mitigation strategies. Professional Paper Jakob M (2005) A size classification for debris flows. Eng Geol 1630:59. https://doi.org/10.3133/pp1630 79:151–161. https://doi.org/10.1016/j.enggeo.2005.01.006 Simpson KA, Stasiuk M, Shimamura K et al (2006) Evidence for catastrophic Le Corvec N, Muirhead JD, White JDL (2018) Shallow magma volcanic debris flows in Pemberton Valley, British Columbia. Can J Earth diversions during explosive diatreme-forming eruptions. Nat Com- Sci 43:679–689. https://doi.org/10.1139/e06-026 mun 9:1–8. https://doi.org/10.1038/s41467-018-03865-x Stasiuk M, Russell JK, Hickson C (1996) Distribution, nature, and Maccaferri F, Richter N, Walter TR (2017) The effect of giant lateral origins of the 2400 BP eruption products of Mount Meager, British collapses on magma pathways and the location of volcanism. Nat Columbia: linkages between magma chemistry and eruption. Geol Commun 8:1–11. https://doi.org/10.1038/s41467-017-01256-2 Surv Canada Bullettin 486:27. https://doi.org/10.4095/208165 Read PB (1979) Geology of Meager Creek Geothermal Area, British Venugopal S, Moune S, William-Jones G et al (2017) Gas emissions Columbia. Geol Surv Canada Open File 603, map, 1 sheet, scale and magma source of the mount meager volcanic complex, 120,000. https://doi.org/https://doi.org/10.4095/129507 Garibaldi , BC. In: IAVCEI scientific assembly. Roberti G, Friele van P, Wyk de Vries B, et al 2017 Rheological Portland, OR. 14–18 August 2017 evolution of the mount meager 2010 debris avalanche, southwestern Venugopal S, Moune S, Williams-Jones G et al (2020) Two distinct british columbia Geosphere 13(2): 369–390. 0.1130/GES01389.1 mantle sources beneath the : Insight from Roberti G, Ward B, van Wyk de Vries B, et al (2018a) Landslides and olivine-hosted melt inclusions Chem Geol 532 10.1016/j. glacier retreat at Mt. Meager volcano: Hazard and risk challenges. chemgeo.2019.119346 In: Geohazards7 Engineering Resiliency in a Changing Climate. Zurek J, Williams-Jones G, Trusdell F, Martin S (2015) The origin of Mauna June 3 - 6, 2018 Canmore, AB Loa’sNīnole Hills: Evidence of rift zone reorganization. Geophys Res Roberti G, Ward B, van Wyk de Vries B, et al (2018) Precursory slope Lett 42:8358–8366. https://doi.org/10.1002/2015GL065863 distress prior to the 2010 Mount Meager landslide, British