Field Survey and Numerical Modelling of the December 22, 2018 Anak Krakatau Tsunami

On December 22, 2018, the eruption and flank collapse of the Anak Krakatau volcano generated a tsunami in the Sunda Strait causing catastrophic damage to uninhabited coastlines proximal to the source. Along the heavily populated shores of Banten and Lampung provinces in Java and Sumatra, tsunami waves caused severe damage, extensive inundation and more than 430 deaths. An international tsunami survey team (ITST) deployed 6 weeks after the event documented the tsunami effects including runup heights, flow depths and inundation distances, as well as sediment deposition patterns and impacts on infrastructure and the natural environment. The team also interviewed numerous eyewitnesses and educated residents about tsunami hazards. This ITST was the first to visit and document the extreme tsunami effects on the small islands in the Sunda Strait closest to Anak Krakatau (Rakata, Panjang, Sertung, Sebesi and Panaitan). Along the steep slopes of Rakata and Sertung islands, located less than 5 km from and facing directly towards the southwest flank of Anak Krakatau, all of the dense coastal vegetation was stripped to bare earth up to elevations of more than 80 m, while on the northeast tip of Sertung Island, facing away from the source, a single tree remained standing after flow depths of > 11 m above ground struck there. The runup distributions on the islands encircling Anak Krakatau highlight the directivity of the tsunami source suggesting that the collapse occurred towards the southwest. This manifested as tsunami runup of 10 m that stripped away coastal forests to bare rock for up to 400 m inland in the Ujung Kulon National Park, located 50 km to the south-southwest. Inundation and damage were mostly limited to within 400 m of the shoreline, likely the result of the relatively short wavelengths caused by the landslide generated tsunami. A significant variation in tsunami impact was observed along the shorelines of the Sunda Strait, with runup heights rapidly decreasing with distance from the inferred tsunami source. To model the event we applied a hot-start initial condition that roughly reproduced the measured tsunami runup heights along Rakata and Sertung. The waveforms were then propagated through the Sunda Straight using a Boussinesq-type wave model. The results showed a good fit to the observed heights along the Java and Sumatra coastlines, the northern