No abstract

On Dec. 22, 2018, at approximately 20:55–57 local time, Anak Krakatau volcano, located in the Sunda Straits of Indonesia, experienced a major lateral collapse during a period of eruptive activity that began in June. The collapse discharged volcaniclastic material into the 250 m deep caldera southwest of the volcano, which generated a tsunami with runups of up to 13 m on the adjacent coasts of Sumatra and Java. The tsunami caused at least 437 fatalities, the greatest number from a volcanically-induced tsunami since the catastrophic explosive eruption of Krakatau in 1883 and the sector collapse of Ritter Island in 1888. For the first time in over 100 years, the 2018 Anak Krakatau event provides an opportunity to study a major volcanically-generated tsunami that caused widespread loss of life and significant damage. Here, we present numerical simulations of the tsunami, with state-of the-art numerical models, based on a combined landslide-source and bathymetric dataset. We constrain the geometry and magnitude of the landslide source through analyses of pre- and post-event satellite images and aerial photography, which demonstrate that the primary landslide scar bisected the Anak Krakatau volcano, cutting behind the central vent and removing 50% of its subaerial extent. Estimated submarine collapse geometries result in a primary landslide volume range of 0.22–0.30 km 3 , which is used to initialize a tsunami generation and propagation model with two different landslide rheologies (granular and fluid). Observations of a single tsunami, with no subsequent waves, are consistent with our interpretation of landslide failure in a rapid, single phase of movement rather than a more piecemeal process, generating a tsunami which reached nearby coastlines within ~30 minutes. Both modelled rheologies successfully reproduce observed tsunami characteristics from post-event field survey results, tide gauge records, and eyewitness reports, suggesting our estimated landslide volume range is appropriate. This event highlights the significant hazard posed by relatively small-scale lateral volcanic collapses, which can occur en-masse , without any precursory signals, and are an efficient and unpredictable tsunami source. Our successful simulations demonstrate that current numerical models can accurately forecast tsunami hazards from these events. In cases such as Anak Krakatau’s, the absence of precursory warning signals together with the short travel time following tsunami initiation present a major challenge for mitigating tsunami coastal impact.

The uncertainties associated with tephra thickness measurements are calculated and implications for volume estimates are presented. Statistical methods are used to analyse the large dataset of Walker and Croasdale J Geol Soc 127:17-55, 1971 of the Fogo A plinian deposit, São Miguel, Azores. Dirichlet tessellation demonstrates that Walker and Croasdale's measurements are highly clustered spatially and the area represented by a single measurement ranges between 0.5 and 10 km 2 . K-means cluster analysis shows that lower thickness uncertainties are associated with closely spaced measurements. Re-examination and analysis of Fogo A fall deposits show thickness uncertainties are about 9% for measured thickness while uncertainty associated with natural variance ranges, between 10 and 40%, with an average error of 30%. Correlations between measurement uncertainties and natural variance are complex and depend on a unit's thickness, position within a succession and distance from source. Normative error increases as tephra thickness decreases. The degree to which thickness measurement error impacts on volume uncertainty depends on the number of measurements within a given dataset and their associated uncertainty. The uncertainty in volume associated with thickness uncertainty calculated herein for Fogo A is 1.3%, equivalent to a volume of 0.02 km 3 . However uncertainties associated with smaller datasets can be much larger; for example typically exceeding 10% for less than 20 data points.

The 39 ka Campanian Ignimbrite (CI) super-eruption was the largest volcanic eruption of the past 200 ka in Europe. Tephra deposits indicate two distinct plume forming phases, Plinian and co-ignimbrite, characteristic of many caldera-forming eruptions. Previous numerical studies have characterized the eruption as a single-phase event, potentially leading to inaccurate assessment of eruption dynamics. To reconstruct the volume, intensity, and duration of the tephra dispersal, we applied a computational inversion method that explicitly accounts for the Plinian and co-ignimbrite phases and for gravitational spreading of the umbrella cloud. To verify the consistency of our results, we performed an additional single-phase inversion using an independent thickness dataset. Our better-fitting two-phase model suggests a higher mass eruption rate than previous studies, and estimates that 3/4 of the total fallout volume is co-ignimbrite in origin. Gravitational spreading of the umbrella cloud dominates tephra transport only within the first hundred kilometres due to strong stratospheric winds in our best-fit wind model. Finally, tephra fallout impacts would have interrupted the westward migration of modern hominid groups in Europe, possibly supporting the hypothesis of prolonged Neanderthal survival in South-Western Europe during the Middle to Upper Palaeolithic transition.

Tephra deposits in the deep sea can survive undisturbed for long periods of time and, on regional scales, tend to be much better preserved than their subaerial counterparts. In this study, grain size distributions and thicknesses of tephra deposits from the Campanian Ignimbrite (CI) eruption (39 000 yr BP; magnitude c. 7.7) preserved in thirty-three deep sea cores are analysed to infer key eruption parameters. Distal deep sea tephra thickness data show an exponential decrease with distance from source. Such trends are difficult to identify in distal subaerial data owing to reworking and limited exposure. We find that tephra grain size distributions are much less affected by depositional environment than thickness, with trends that are consistent across distal subaerial, lacustrine and deep sea environments. The CI layer exhibits bimodal grain size distributions to distances of c. 1000 km, after which it becomes unimodal. Such trends can be related to different mechanisms of tephra transport in the atmosphere, whereby at proximal to medial distances the Plinian and co-ignimbrite phases produce distinct plumes. Within 150 and 900 km from source, Plinian tephra constitutes 40 + 5% of the deposit volume. Beyond this region where coarse particles are deposited, the plumes merge and fines derived largely from co-ignimbrite elutriation are spread in the atmosphere at velocities greater than the settling velocities of the particles.