skip to main content

Title: Modelling of the tsunami from the December 22, 2018 lateral collapse of Anak Krakatau volcano in the Sunda Straits, Indonesia

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 more » result in a primary landslide volume range of 0.22–0.30 km3, 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 occuren-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.

« less
; ; ; ; ; ; ; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Scientific Reports
Nature Publishing Group
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract As demonstrated at Anak Krakatau on December 22 nd , 2018, tsunamis generated by volcanic flank collapse are incompletely understood and can be devastating. Here, we present the first high-resolution characterisation of both subaerial and submarine components of the collapse. Combined Synthetic Aperture Radar data and aerial photographs reveal an extensive subaerial failure that bounds pre-event deformation and volcanic products. To the southwest of the volcano, bathymetric and seismic reflection data reveal a blocky landslide deposit (0.214 ± 0.036 km 3 ) emplaced over 1.5 km into the adjacent basin. Our findings are consistent with en-masse lateral collapse with a volume ≥0.175 km 3 , resolving several ambiguities in previous reconstructions. Post-collapse eruptions produced an additional ~0.3 km 3 of tephra, burying the scar and landslide deposit. The event provides a model for lateral collapse scenarios at other arc-volcanic islands showing that rapid island growth can lead to large-scale failure and that even faster rebuilding can obscure pre-existing collapse.
  2. Abstract

    Improvement of our understanding of the role of ground deformation due to viscoelastic relaxation following eruptions is important, as the generated signals can resemble renewed magma inflow. We study post‐eruptive unrest at the subglacial Bárðarbunga volcano, Iceland, after a caldera collapse and major magma drainage in 2014–2015. Elevated seismicity began about 6 months after the eruption ended, including nine Mlw> 4.5 earthquakes. Global Navigation Satellite System and Sentinel‐1 Interferometric Synthetic Aperture Radar geodesy are applied to evaluate post‐eruptive ground deformation from 2015 to 2018. Horizontal velocities locally exceed 10 cm/year and rapidly decay with distance away from the caldera. We explore two end‐member models and their combination to explain the post‐eruptive deformation field: 1) viscoelastic relaxation caused by the co‐eruptive caldera collapse and magma withdrawal, and 2) renewed magma inflow. We find parameter combinations for each model that explain the observed ground deformation. The purely viscoelastic relaxation model, consisting of a half‐space composed of a 7‐km thick elastic layer on top of a viscoelastic layer with a viscosity of 3.0 × 1018 Pa s reproduces broadly the observations. A simple magma inflow model consisting of a single point source with an inflow rate of 1 × 107 m3/year at 0.7 km depth broadly fits the observations, butmore »may be unrealistic. A more elaborate model of magma inflow into a 10‐km deep sill combined with slip on the caldera ring fault explains the observations well. Our results suggest that the co‐eruptive deformation field is likely influenced by viscoelastic relaxation, renewed magma inflow, or a combination of both processes.

    « less
  3. On 22 December 2018, a devastating tsunami struck Sunda Strait, Indonesia without warning, leaving 437 dead and thousands injured along the western Java and southern Sumatra coastlines. Synthetic aperture radar and broadband seismic observations demonstrate that a small, <~0.2 km 3 landslide on the southwestern flank of the actively erupting volcano Anak Krakatau generated the tsunami. The landslide did not produce strong short-period seismic waves; thus, precursory ground shaking did not provide a tsunami warning. The source of long-period ground motions during the landslide can be represented as a 12° upward-dipping single-force directed northeastward, with peak magnitude of ~6.1 × 10 11 N and quasi-sinusoidal time duration of ~70 s. Rapid quantification of a landslide source process by long-period seismic wave inversions for moment-tensor and single-force parameterizations using regional seismic data available within ~8 min can provide a basis for future fast tsunami warnings, as is also the case for tsunami earthquakes.
  4. Fault friction is central to understanding earthquakes, yet laboratory rock mechanics experiments are restricted to, at most, meter scale. Questions thus remain as to the applicability of measured frictional properties to faulting in situ. In particular, the slip-weakening distance d c strongly influences precursory slip during earthquake nucleation, but scales with fault roughness and is challenging to extrapolate to nature. The 2018 eruption of Kīlauea volcano, Hawaii, caused 62 repeatable collapse events in which the summit caldera dropped several meters, accompanied by M W 4.7 to 5.4 very long period (VLP) earthquakes. Collapses were exceptionally well recorded by global positioning system (GPS) and tilt instruments and represent unique natural kilometer-scale friction experiments. We model a piston collapsing into a magma reservoir. Pressure at the piston base and shear stress on its margin, governed by rate and state friction, balance its weight. Downward motion of the piston compresses the underlying magma, driving flow to the eruption. Monte Carlo estimation of unknowns validates laboratory friction parameters at the kilometer scale, including the magnitude of steady-state velocity weakening. The absence of accelerating precollapse deformation constrains d c to be ≤ 10 mm, potentially much less. These results support the use of laboratory frictionmore »laws and parameters for modeling earthquakes. We identify initial conditions and material and magma-system parameters that lead to episodic caldera collapse, revealing that small differences in eruptive vent elevation can lead to major differences in eruption volume and duration. Most historical basaltic caldera collapses were, at least partly, episodic, implying that the conditions for stick–slip derived here are commonly met in nature.« less
  5. Abstract

    The Baiyun slide complex contains geological evidence for some of the largest landslide ever discovered in the continental slopes of the South China Sea. High‐resolution seismic data suggest that a variety of landslides with varied scales have occurred repeatedly in this area. The largest landslide reconstructed from bathymetric and seismic data has an estimated spatial coverage of ~5,500 km2and a conservative volume of ~1,035 km3. Here, using geomorphological and geotechnical data, we construct a series of probable landslide scenarios and assess their tsunamigenic capacity. By treating the slides as deformable mudflows, we simulate the dynamics of landslide movements. The simulated landslide motions match the geophysical observations interpreted in previous studies. Particularly, we are able to reproduce the spatial distribution of observed runout, including the distance, shape, and deposit thickness, for the most credible slide scenario. We investigate tsunami impacts generated by different slide scenarios and highlight the importance of initial water depth, sliding direction, and nearshore bathymetry. The worst‐case scenario is capable of producing basin‐wide tsunami, with maximum wave amplitudes reaching ~5 m near Hong Kong and Macau, 1–3 m in western Philippines, and at least 1 m along central Vietnam, southeast Hainan, and southern Taiwan. The most noticeablemore »phenomenon we observed is that the southern Chinese coast is the hardest‐hit region in all the simulated scenarios regardless of the diverse slide features. We conclude that the persistence of high tsunami impact is caused by the unique bathymetric feature of the wide continental shelf in front of southern China.

    « less