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.

Revised Minoan eruption volume as benchmark for large volcanic eruptions.

Despite their global societal importance, the volumes of large-scale volcanic eruptions remain poorly constrained. Here, we integrate seismic reflection and P-wave tomography datasets with computed tomography-derived sedimentological analyses to estimate the volume of the iconic Minoan eruption. Our results reveal a total dense-rock equivalent eruption volume of 34.5 ± 6.8 km³, which encompasses 21.4 ± 3.6 km³ of tephra fall deposits, 6.9 ± 2 km³ of ignimbrites, and 6.1 ± 1.2 km³ of intra-caldera deposits. 2.8 ± 1.5 km³ of the total material consists of lithics. These volume estimates are in agreement with an independent caldera collapse reconstruction (33.1 ± 1.2 km³). Our results show that the Plinian phase contributed most to the distal tephra fall, and that the pyroclastic flow volume is significantly smaller than previously assumed. This benchmark reconstruction demonstrates that complementary geophysical and sedimentological datasets are required for reliable eruption volume estimates, which are necessary for regional and global volcanic hazard assessments.

Slow slip source characterized by lithological and geometric heterogeneity.

Slow slip events (SSEs) accommodate a significant proportion of tectonic plate motion at subduction zones, yet little is known about the faults that actually host them. The shallow depth (<2 km) of well-documented SSEs at the Hikurangi subduction zone offshore New Zealand offers a unique opportunity to link geophysical imaging of the subduction zone with direct access to incoming material that represents the megathrust fault rocks hosting slow slip. Two recent International Ocean Discovery Program Expeditions sampled this incoming material before it is entrained immediately down-dip along the shallow plate interface. Drilling results, tied to regional seismic reflection images, reveal heterogeneous lithologies with highly variable physical properties entering the SSE source region. These observations suggest that SSEs and associated slow earthquake phenomena are promoted by lithological, mechanical, and frictional heterogeneity within the fault zone, enhanced by geometric complexity associated with subduction of rough crust.

Submarine slope failures due to pipe structure formation.

There is a strong spatial correlation between submarine slope failures and the occurrence of gas hydrates. This has been attributed to the dynamic nature of gas hydrate systems and the potential reduction of slope stability due to bottom water warming or sea level drop. However, 30 years of research into this process found no solid supporting evidence. Here we present new reflection seismic data from the Arctic Ocean and numerical modelling results supporting a different link between hydrates and slope stability. Hydrates reduce sediment permeability and cause build-up of overpressure at the base of the gas hydrate stability zone. Resulting hydro-fracturing forms pipe structures as pathways for overpressured fluids to migrate upward. Where these pipe structures reach shallow permeable beds, this overpressure transfers laterally and destabilises the slope. This process reconciles the spatial correlation of submarine landslides and gas hydrate, and it is independent of environmental change and water depth.

Arctic megaslide at presumed rest.

Slope failure like in the Hinlopen/Yermak Megaslide is one of the major geohazards in a changing Arctic environment. We analysed hydroacoustic and 2D high-resolution seismic data from the apparently intact continental slope immediately north of the Hinlopen/Yermak Megaslide for signs of past and future instabilities. Our new bathymetry and seismic data show clear evidence for incipient slope instability. Minor slide deposits and an internally-deformed sedimentary layer near the base of the gas hydrate stability zone imply an incomplete failure event, most probably about 30000 years ago, contemporaneous to or shortly after the Hinlopen/Yermak Megaslide. An active gas reservoir at the base of the gas hydrate stability zone demonstrate that over-pressured fluids might have played a key role in the initiation of slope failure at the studied slope, but more importantly also for the giant HYM slope failure. To date, it is not clear, if the studied slope is fully preconditioned to fail completely in future or if it might be slowly deforming and creeping at present. We detected widespread methane seepage on the adjacent shallow shelf areas not sealed by gas hydrates.