Abstract Volcanic-induced tsunamis have a potentially devastating impact, especially in densely populated and/or touristic coastal areas. Kolumbo submarine volcano (Greece) experienced in 1650 CE an explosive eruption with eyewitnesses’ accounts of major tsunamis along the coasts of Santorini (Thera) and other islands. We present a scenario-based tsunami hazard assessment at this volcano based on existing simulations from literature and new simulations of tsunamis triggered by a less investigated but important mechanism, i.e. submarine landslides on the volcano flanks or within its crater. Simulations results show that the remobilization of a landslide volume of 150–300 Mm3 inside the crater can produce tsunami waves larger than 10 m high along the NE coast of Thera and of the order of 5 m along the E and SE coasts. The expected tsunami arrival time ranges from 2–3 minutes along the NE coast of Thera up to 8–10 min on its SE coast. Such scenarios produce inundation areas consistent with those reconstructed for the 1650 CE event, and tsunami waves propagating inland at velocities from 2 to 12 m/s. Simulation results also suggest that, given the landslide parameters assumed, it is unlikely to mobilize a landslide with a large volume from the SW-facing Kolumbo crater slopes, given the relatively gentle topo-bathymetry of this area. The study findings are relevant based on the outcomes of the expert elicitation exercise carried out in parallel, which indicate that chances of having waves larger than 1 m high on the NE coast of Thera have median probabilities of 50–60%.
Scientists monitoring active volcanoes are increasingly required to provide decision support to civil authorities during periods of unrest. As the extent and resolution of monitoring improves, the process of jointly interpreting multiple strands of indirect evidence becomes increasingly complex. Similarities with uncertainties in medical diagnosis suggest a formal evidence-based approach, whereby monitoring data are analysed synoptically to provide probabilistic hazard forecasts. A statistical tool to formalize such inferences is the Bayesian Belief Network (BBN). By explicitly representing conditional dependencies between the volcanological model and observations, BBNs use probability theory to treat uncertainties in a rational and auditable manner, as warranted by the strength of the scientific evidence. A retrospective analysis is given for the 1976 Guadeloupe crisis, using a BBN to provide inferential assessment of the state of the evolving magmatic system and probability of incipient eruption. Conditional dependencies are characterized quantitatively by structured expert elicitation. Analysis of the available monitoring data suggests that at the height of the crisis the probability of magmatic intrusion was high, in accordance with scientific thinking at the time. The corresponding probability of magmatic eruption was elevated in July and August 1976 and signs of precursory activity were justifiably cause for concern. However, collective uncertainty about the future course of the crisis was also substantial. Of all the possible scenarios, the most likely outcome evinced by interpretation of observations on 31 August 1976 was 'no eruption' (mean probability 0.5); the chance of a magmatic eruption/blast had an estimated mean probability of ~0.4. There was therefore no evidential basis for asserting one scenario to be significantly more likely than another. Our analysis adds objective probabilistic expression to the volcanological narrative at the time of the 1976 crisis, and demonstrates that a formal evidential case could have supported the authorities' concerns about public safety and decision to evacuate. Revisiting the episode highlights many challenges for modern, contemporary decision making under conditions of considerable uncertainty, and suggests the BBN is a suitable framework for marshalling multiple, uncertain observations, model results and interpretations. The formulation presented here can be developed as a tool for ongoing use in the volcano observatory.
A non-parametric statistical approach is used to assess the global recording rate for large (M4+) stratovolcano eruptions in a modern database, LaMEVE (v3.1). This approach imposes minimal structure on the shape of the recording rate through time. We find that recording rates have declined rapidly, going backwards in time. Prior to the year 1600 they are below 50 %, and prior to 1100 they are below 20 %. Even in the recent past, e.g. the 1800s, they are likely to be less than 100 %. The assessment for very large (M5+) eruptions is more uncertain, due to the scarcity of events.