Cascading events during the 1650 tsunamigenic eruption of Kolumbo volcano

Volcanic eruptions can trigger tsunamis, which may cause significant damage to coastal communities and infrastructure. Tsunami generation during volcanic eruptions is complex and often due to a combination of processes. The 1650 eruption of the Kolumbo submarine volcano triggered a tsunami causing major destruction on surrounding islands in the Aegean Sea. However, the source mechanisms behind the tsunami have been disputed due to difficulties in sampling and imaging submarine volcanoes. Here we show, based on three-dimensional seismic data, that ~1.2 km³ of Kolumbo's northwestern flank moved 500-1000 m downslope along a basal detachment surface. This movement is consistent with depressurization of the magma feeding system, causing a catastrophic explosion. Numerical tsunami simulations indicate that only the combination of flank movement followed by an explosive eruption can explain historical eyewitness accounts. This cascading sequence of natural hazards suggests that assessing submarine flank movements is critical for early warning of volcanogenic tsunamis.

A scenario-based volcanic hazard assessment for the Mount Meager Volcanic Complex, British Columbia

The Mount Meager Volcanic Complex (Mount Meager) is a glacier-clad stratovolcanic system in southwestern British Columbia which last erupted over 2400 years ago (VEI 4). While this is Canada's most recent major explosive eruption, most past research on Mount Meager has focused on its numerous and large volume landslides and thus the volcanic hazard characteristics remain understudied. Here we present a suite of scenario-based hazard maps and an assessment addressing a range of potential future explosive eruptions and associated hazards. In order to overcome limited knowledge of the eruptive history, numerical models have been used to simulate the primary syneruptive hazards of concern (dome-collapse pyroclastic density currents, lahars and tephra fallout) largely utilizing eruption parameters from analogous volcanoes, i.e., glacier-clad stratovolcanoes in a subduction zone setting. This study provides a framework for similar volcanic hazard studies where geologic data is limited, funds are minimal, and access is difficult. Furthermore, this sets the stage for recognizing volcanic hazards in the Canadian landscape, providing a resource to prepare for and mitigate potential impacts well in advance of a crisis situation.

Effect of Boulder Size on Debris Flow Impact Pressure Using a CFD-DEM Numerical Model

Debris flows (DFs) are dangerous events that can cause the complete destruction of buildings and infrastructure, such as bridges; DFs therefore represent a high risk to public safety in exposed areas. The impact pressures due to these flows are essentially determined by the flow height, velocity and density, but other parameters that are less often considered are also involved. We developed a numerical model to evaluate the impact pressure of mass flows, focusing on a better description of the influence of the blocks transported in these flows: the block size strongly influences the impact pressure, which has a strong effect on structural damage. The numerical model proposed considers a staggered, loosely one-way granular–fluid coupling based on a distinct-element-method code, using the separate simulation results of a computing fluid dynamics code used to model the fluid phase. This model estimates the impact pressure distribution due to blocks at the local scale of the obstacle; the pressure due to the fluid phase can be added afterwards. The pressure applied by the DF increased with the maximum block size for a given set of DF characteristics: velocity, height and apparent density. The vulnerability of a given structure depends on the intensity of DFs: the pressure applied on the structure is one of considerable intensity. The existing vulnerability functions are interpreted in the light of the results obtained with the numerical model. This interpretation highlights the need to integrate new parameters in the intensity to better evaluate structures’ vulnerability to debris flows.

Pyroclastic Density Current Hazard Assessment and Modeling Uncertainties for Fuego Volcano, Guatemala

On 3 June 2018, Fuego volcano experienced a VEI = 3 eruption, which produced a pyroclastic density current (PDC) that devastated the La Réunion resort and the community of Los Lotes, resulting in over 100 deaths. To evaluate the potential hazard to the population centers surrounding Fuego associated with future PDC emplacement, we used an integrated remote sensing and flow modeling-based approach. The predominate PDC travel direction over the past 15 years was investigated using thermal infrared (TIR) data from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) instrument validated with ground reports from the National Institute of Seismology, Volcanology, Meteorology, and Hydrology (INSIVUMEH), the government agency responsible for monitoring. Two different ASTER-derived digital elevation model (DEM) products with varying levels of noise were also used to assess the uncertainty in the VolcFlow model results. Our findings indicate that the recent historical PDC travel direction is dominantly toward the south and southwest. Population centers in this region of Fuego that are within ~2 km of one of the volcano’s radial barrancas are at the highest risk during future large eruptions that produce PDCs. The ASTER global DEM (GDEM) product has the least random noise and where used with the VolcFlow model, had a significant improvement on its accuracy. Results produced longer flow runout distances and therefore better conveys a more accurate perception of risk. Different PDC volumes were then modeled using the GDEM and VolcFlow to determine potential inundation areas in relation to local communities.

Scenario-Based Pyroclastic Density Current Invasion Maps at Poorly Known Volcanoes: A Case Study from Changbaishan (China/North Korea)

Changbaishan volcano (China/North Korea; last eruption in 1903 AD) was responsible for a Volcanic Explosivity Index (VEI) 7 eruption in 946 AD. Approximately 186,000 people live around Changbaishan and 2,000,000 tourists/year visit the volcano. An unrest occurred between 2002 and 2006. Despite the relevant hazard, the eruptive history is poorly known, a condition common to many volcanoes worldwide. Here, we investigate the extension of the areas potentially affected by pyroclastic density currents (PDCs) in case of future eruptions following a scenario-based approach. We perform energy cone runs referred to four scenarios from columns of height 3, 10, 20 and 30 km at different vents. By using global datasets on PDCs, we produce spatial probability maps of PDCs invasion. Empirical laws between covered areas, PDC travelled distances, and heights of collapse are provided. In scenarios 3 and 4, PDCs expand at distances up to 42 km and 85 km, respectively. In scenarios 1 and 2, PDCs invade the touristic area and few main roads. Severe effects emerge from scenarios 3 and 4 with the interruption of the China–North Korea land and aerial connections and PDC. Our approach may serve as guide for the rapid evaluation of the PDC-related hazard at poorly known volcanoes.

Modeling of partial dome collapse of La Soufrière of Guadeloupe volcano: implications for hazard assessment and monitoring

Over the past 9,150 years, at least 9 flank collapses have been identified in the history of La Soufrière of Guadeloupe volcano. On account of the volcano’s current unrest, the possibility of such a flank collapse should not be dismissed in assessing hazards for future eruptive magmatic as well as non-magmatic scenarios. We combine morphological and geophysical data to identify seven unstable structures (volumes ranging from 1 × 10⁶ m³ to 100 × 10⁶ m³), including one that has a volume compatible with the last recorded flank collapse in 1530 CE. We model their dynamics and emplacement with the SHALTOP numerical model and a simple Coulomb friction law. The best-fit friction coefficient to reproduce the 1530 CE event is tan(7°) = 0.13, suggesting the transformation of the debris avalanche into a debris flow, which is confirmed by the texture of mapped deposits. Various friction angles are tested to investigate less water-rich and less mobile avalanches. The most densely populated areas of Saint-Claude and Basse-Terre, and an area of Gourbeyre south of the Palmiste ridge, are primarily exposed in the case of the more voluminous and mobile flank collapse scenarios considered. However, topography has a prominent role in controlling flow dynamics, with barrier effects and multiple channels. Classical mobility indicators, such as the Heim’s ratio, are thus not adequate for a comprehensive hazard analysis.

Statistical emulation of landslide-induced tsunamis at the Rockall Bank, NE Atlantic

Statistical methods constitute a useful approach to understand and quantify the uncertainty that governs complex tsunami mechanisms. Numerical experiments may often have a high computational cost. This forms a limiting factor for performing uncertainty and sensitivity analyses, where numerous simulations are required. Statistical emulators, as surrogates of these simulators, can provide predictions of the physical process in a much faster and computationally inexpensive way. They can form a prominent solution to explore thousands of scenarios that would be otherwise numerically expensive and difficult to achieve. In this work, we build a statistical emulator of the deterministic codes used to simulate submarine sliding and tsunami generation at the Rockall Bank, NE Atlantic Ocean, in two stages. First we calibrate, against observations of the landslide deposits, the parameters used in the landslide simulations. This calibration is performed under a Bayesian framework using Gaussian Process (GP) emulators to approximate the landslide model, and the discrepancy function between model and observations. Distributions of the calibrated input parameters are obtained as a result of the calibration. In a second step, a GP emulator is built to mimic the coupled landslide-tsunami numerical process. The emulator propagates the uncertainties in the distributions of the calibrated input parameters inferred from the first step to the outputs. As a result, a quantification of the uncertainty of the maximum free surface elevation at specified locations is obtained.

Post-eruptive flooding of Santorini caldera and implications for tsunami generation

Caldera-forming eruptions of island volcanoes generate tsunamis by the interaction of different eruptive phenomena with the sea. Such tsunamis are a major hazard, but forward models of their impacts are limited by poor understanding of source mechanisms. The caldera-forming eruption of Santorini in the Late Bronze Age is known to have been tsunamigenic, and caldera collapse has been proposed as a mechanism. Here, we present bathymetric and seismic evidence showing that the caldera was not open to the sea during the main phase of the eruption, but was flooded once the eruption had finished. Inflow of water and associated landsliding cut a deep, 2.0-2.5 km³, submarine channel, thus filling the caldera in less than a couple of days. If, as at most such volcanoes, caldera collapse occurred syn-eruptively, then it cannot have generated tsunamis. Entry of pyroclastic flows into the sea, combined with slumping of submarine pyroclastic accumulations, were the main mechanisms of tsunami production.