skip to main content

Title: The 22 December 2018 tsunami from flank collapse of Anak Krakatau volcano during eruption
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.
; ; ; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Science Advances
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 geometriesmore »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
  2. SUMMARY On 2020 May 2, an Mw = 6.6 earthquake struck about 63 km south of Ierapetra in Crete, Greece. The earthquake generated a small tsunami which agitated local harbours. We studied this event in the context of earthquakes with seismic records in 1908, 1910, 1923, 1952, 2009 and 2013, all of similar magnitudes located south of Crete. Based on an energy-to-moment ratio, our analysis suggests that this event was neither slow nor fast, hence appropriate for using scaling laws to infer seafloor deformations. We also performed a field survey, three days after the event and present field observations from seven locations, including the island of Chrisi, where our highest measurement of 0.95 m was located. Runup along the coast of southern Crete ranged from 0.24 to 0.87 m. One tide gauge record is available for this event, and we did image analysis to obtain accurately timed water surface elevations from eyewitness videos and images. We undertook high-resolution hydrodynamic simulations using published moment tensor solutions to identify the source of the tsunami. Simulations were performed with two models, MOST (a nonlinear shallow water model) and COULWAVE (a Boussinesq-type model), to infer how different approximations of the parent equations of motion affect predictions for tsunamismore »of this size, which are fairly common in the Eastern Mediterranean and routinely trigger Tsunami Service Providers to issue warning messages. Based on the inter-model comparison, we conclude that the shallow-water equations are adequate in modelling this event at the distances considered, suggesting that such codes can be used to infer the tsunami source and to estimate tsunami impacts. Last, our field work revealed lack of knowledge of tsunami hazards, as most eyewitnesses remained near the waterfront, filming the associated unusual water motions instead of taking shelter on high ground.« less
  3. Abstract. Surficial mass wasting events are a hazard worldwide. Seismic and acoustic signals from these often remote processes, combined with other geophysical observations, can provide key information for monitoring and rapid response efforts and enhance our understanding of event dynamics. Here, we present seismoacoustic data and analyses for two very large ice–rock avalanches occurring on Iliamna Volcano, Alaska (USA), on 22 May 2016 and 21 June 2019. Iliamna is a glacier-mantled stratovolcano located in the Cook Inlet, ∼200 km from Anchorage, Alaska. The volcano experiences massive, quasi-annual slope failures due to glacial instabilities and hydrothermal alteration of volcanic rocks near its summit. The May 2016 and June 2019 avalanches were particularly large and generated energetic seismic and infrasound signals which were recorded at numerous stations at ranges from ∼9 to over 600 km. Both avalanches initiated in the same location near the head of Iliamna's east-facing Red Glacier, and their ∼8 km long runout shapes are nearly identical. This repeatability – which is rare for large and rapid mass movements – provides an excellent opportunity for comparison and validation of seismoacoustic source characteristics. For both events, we invert long-period (15–80 s) seismic signals to obtain a force-time representation of the source. We model themore »avalanche as a sliding block which exerts a spatially static point force on the Earth. We use this force-time function to derive constraints on avalanche acceleration, velocity, and directionality, which are compatible with satellite imagery and observed terrain features. Our inversion results suggest that the avalanches reached speeds exceeding 70 m s−1, consistent with numerical modeling from previous Iliamna studies. We lack sufficient local infrasound data to test an acoustic source model for these processes. However, the acoustic data suggest that infrasound from these avalanches is produced after the mass movement regime transitions from cohesive block-type failure to granular and turbulent flow – little to no infrasound is generated by the initial failure. At Iliamna, synthesis of advanced numerical flow models and more detailed ground observations combined with increased geophysical station coverage could yield significant gains in our understanding of these events.« less
  4. Abstract

    The Kalapana, Hawaii,MW7.7 earthquake on November 29, 1975 generated a local tsunami with at least 14.3 m runup on the southeast shore of Hawaii Island adjacent to Kilauea Volcano. This was the largest locally generated tsunami since the great 1868 Ka'u earthquake located along‐shore to the southwest. Well‐recorded tide gauge and runup observations provide a key benchmark for studies of statewide tsunami hazards from actively deforming southeast Hawaii Island. However, the source process of the earthquake remains controversial, with coastal landsliding and/or offshore normal or thrust faulting mechanisms having been proposed to reconcile features of seismic, geodetic, and tsunami observations. We utilize these diverse observations for the 1975 Kalapana earthquake to deduce a compound faulting model that accounts for the overall tsunamigenesis, involving both landslide block faulting along the shore and slip on the island basal décollement. Thrust slip of 4.5–8.0 m on the offshore décollement produces moderate near‐field runup but controls the far‐field tsunami. The slip distribution implies that residual strain energy was available for the May 4, 2018MW7.2 thrust earthquake during the Kilauea‐East Rift Zone eruption. Local faulting below land contributes to geodetic and seismic observations, but is non‐tsunamigenic and not considered. Slip of 4–10 m on landslide‐like faults, whichmore »extend from the Hilina Fault Zone scarp to offshore shallowly dipping faults reaching near the seafloor, triples the near‐field tsunami runup. This compound model clarifies the roles of the faulting components in assessing tsunami hazards for the Hawaiian Islands.

    « less
  5. SUMMARY We recently found the original Omori seismograms recorded at Hongo, Tokyo, of the 1922 Atacama, Chile, earthquake (MS = 8.3) in the historical seismogram archive of the Earthquake Research Institute (ERI) of the University of Tokyo. These recordings enable a quantitative investigation of long-period seismic radiation from the 1922 earthquake. We document and provide interpretation of these seismograms together with a few other seismograms from Mizusawa, Japan, Uppsala, Sweden, Strasbourg, France, Zi-ka-wei, China and De Bilt, Netherlands. The 1922 event is of significant historical interest concerning the cause of tsunami, discovery of G wave, and study of various seismic phase and first-motion data. Also, because of its spatial proximity to the 1943, 1995 and 2015 great earthquakes in Chile, the 1922 event provides useful information on similarity and variability of great earthquakes on a subduction-zone boundary. The 1922 source region, having previously ruptured in 1796 and 1819, is considered to have significant seismic hazard. The focus of this paper is to document the 1922 seismograms so that they can be used for further seismological studies on global subduction zones. Since the instrument constants of the Omori seismographs were only incompletely documented, we estimate them using the waveforms of the observed records,more »a calibration pulse recorded on the seismogram and the waveforms of better calibrated Uppsala Wiechert seismograms. Comparison of the Hongo Omori seismograms with those of the 1995 Antofagasta, Chile, earthquake (Mw = 8.0) and the 2015 Illapel, Chile, earthquake (Mw = 8.3) suggests that the 1922 event is similar to the 1995 and 2015 events in mechanism (i.e. on the plate boundary megathrust) and rupture characteristics (i.e. not a tsunami earthquake) with Mw = 8.6 ± 0.25. However, the initial fine scale rupture process varies significantly from event to event. The G1 and G2, and R1 and R2 of the 1922 event are comparable in amplitude, suggesting a bilateral rupture, which is uncommon for large megathrust earthquakes.« less