Global optimization of a numerical two-layer model using observed data: a case study of the 2018 Sunda Strait tsunami
Tóm tắt
Following the eruption of Mount Anak Krakatau, a considerable landslide occurred on the southwestern part of the volcano and, upon entering the sea, generated a large tsunami within the Sunda Strait, Indonesia, on December 22, 2018. This tsunami traveled ~ 5 km across the strait basin and inundated the shorelines of Sumatra and Java with a vertical runup reaching 13 m. Following the event, observed field data, GPS measurements of the inundation, and multibeam echo soundings of the bathymetry within the strait were collected and publicly provided. Using this dataset, numerical modeling of the tsunami was conducted using the two-layer (soil and water) TUNAMI-N2 model based on a combination of landslide sources and bathymetry data. The two-layer model was implemented to nest the grid system using the finest grid size of 20 m. To constrain the unknown landslide parameters, the differential evolution (DE) global optimization algorithm was applied, which resulted in a parameter set that minimized the deviation from the measured bathymetry after the event. The DE global optimization procedure was effective at determining the landslide parameters for the model with the minimum deviation from the measured seafloor. The lowest deviation from the measured bathymetry was obtained for the best-fitting parameters: a maximum landslide thickness of 301.2 m and a landslide time of 10.8 min. The landslide volume of 0.182 km3 estimated by the best-fitting parameters shows that the tsunami flow depth could have reached 3–10 m along the shore with a K value of 0.89, although the simulated flow depths were underestimated in comparison with the observation data. According to the waveforms, the general wave pattern was well reproduced at tide gauges during the event. A large number of objective function evaluations were necessary to locate the minimum with the DE procedure to fix the grid cell size to 20 m; this limited the accuracy of the obtained parameter values for the two-layer model. Moreover, considering the generalizations in the modeling of landslide movements, the impact landslide time and thickness must be carefully calculated to obtain a suitable accuracy.
Tài liệu tham khảo
Aburaya T, Imamura F (2002) The proposal of a tsunami run-up simulation using combined equivalent roughness. Annu J Coast Ocean Eng 49:276–280 (in Japanese)
Aida I (1978) Reliability of a tsunami source model derived from fault parameters. J Phys Earth 26:57–73
Baba T, Allgeyer S, Hossen J, Cummins PR, Tsushima H, Imai K, Yamashita K, Kato T (2017) Accurate numerical simulation of the far-field tsunami caused by the 2011 Tohoku earthquake, including the effects of Boussinesq dispersion, seawater density stratification, elastic loading, and gravitational potential change. Ocean Model 111:46–54
Choi B, Pelinovsky E, Kim K, Lee J (2003) Simulation of the trans-oceanic tsunami propagation due to the 1883 Krakatau volcanic eruption, Copernicus Publications on behalf of the European Geosciences Union. Natural Hazards and Earth System Science 3(5):321–332
Copeland R (2000) Determinate of flow resistance coefficient due to shrubs and woody vegetation, ERDC/CHL HETN-II-3, US Army Corps of Engineers, p 1–7
Dutta D, Alam J, Umeda K, Hayashi M (2007) A two-dimensional hydrodynamic model for flood inundation simulation: a case study in the lower Mekong River basin. Hydrol Process 21:1223–1237
Fine IV, Rabinovich AB, Bornhold BD, Thomson RE, Kulikov EA (2005) The Grand Banks landslide-generated tsunami of November 18, 1929: Preliminary analysis and numerical modeling. Marine Geology 215(1D2):45–57
Fukui N, Prasetyo A, Mori N (2019) Numerical modeling of tsunami inundation using upscaled urban roughness parameterization. Coast Eng 152:103534
Giachetti T, Paris R, Kelfoun K, Ontowirjo B (2012) Tsunami hazard related to a flank collapse of Anak Krakatau Volcano, Sunda Strait, Indonesia. Geol Soc Lond Spec Publ 361:79–90
Grilli ST, Harris JC, Tayebeh S, Bakhsh T, Masterlark TL, Kyriakopoulos C (2012) Numerical simulation of the 2011 Tohoku tsunami based on a new transient FEM co-seismic source: Comparison to far- and near-field observations. Pure Appl Geophys 170(6–8):1333–1359
Grilli S, Tappin D, Carey S, Watt S, Ward S, Grilli A (2019) Modelling of the tsunami from the December 22, 2018 lateral collapse of Anak Krakatau volcano in the Sunda Straits. Indonesia Sci Rep 9:11946
Gylfadottir SS, Kim J, Helgason JK, Brynjolfsson S, Hoskuldsson A, Johannesson T, Harbitz CB, Løvholt F (2017) The 2014 Lake Askja rockslide-induced tsunami: Optimization of numerical tsunami model using observed data. J Geophys Res Oceans 122:4110–4122
Harbitz CB, Parker G, Elverhøi A, Marr JG, Mohrig D, Harff PA (2003) Hydroplaning of subaqueous debris flows and glide blocks: analytical solutions and discussion. J Geophys Res Solid Earth 108(B7):2349
Heidarzadeh M, Krastel S, Yalciner AC (2014) The state-of-the-art numerical tools for modeling landslide tsunamis: a short review. In: Krastel S. et al. (eds), Submarine mass movements and their consequences. Advances in natural and technological hazards research, 37, Chapter 43 (pp. 483–959). USA: Springer International Publishing
Heidarzadeh M, Ishibe T, Sandanbata O, Muhari A, Wijanarto A (2020) Numerical modeling of the subaerial landslide source of the 22 December 2018 Anak Krakatoa volcanic tsunami Indonesia. Ocean Eng 20:195
Heinrich P, Boudon G, Komorowski JC, Sparks RSJ, Herd R, Voight B (2001) Numerical simulation of the December 1997 Debris Avalanche in Montserrat. Lesser Antilles Geophys Res Lett 28(13):2529–2532
Iglesias O, Lastras G, Canals M, Olabarrieta M, González M, Aniel-Quiroga I (2012) The BIG’95 submarine landslide-generated tsunami; a numerical simulation. J Geol 120(1):31–48
Imamura F (1995) Review of tsunami with a finite difference method, In: Long-wave runup models. World Scientific Pub Co Inc: Singapore, pp 25–42
Imamura F, Imteaz MA (1995) Long waves in two-layers: governing equations and numerical model. Sci Tsunami Hazards 13(1):3–24
Ioki K, Tanioka Y, Yanagisawa H, Kawakami G (2019) Numerical simulation of the landslide and tsunami due to the 1741 Oshima-oshima eruption in Hokkaido, Japan. J Geophys Res Solid Earth 124:1991–2002
Iverson RM, Denlinger RP (2001) Flow of variably fluidized granular masses across three-dimension terrain 1. Coulomb mixture theory. J Geophys Res 106(1):537–552
Japan Society of Civil Engineering (2002) Tsunami assessment method for nuclear power plants in Japan, 72. https://www.jsce.or.jp/committee/ceofnp/Tsunami/eng/JSCE_Tsunami_060519.pdf.
Kirby J, Shi F, Nicolsky D, Misra S (2016) The 27 April 1975 Kitimat, British Columbia, submarine landslide tsunami: A comparison of modeling approaches. Landslides 13:1421–1434
Kowalik Z, Murty TS (1993) Numerical modeling of ocean dynamics. World Scientific Pub, Singapore, p 481
Lastras G, De Blasio FV, Canals M, Elverhøi A (2005) Conceptual and numerical modeling of the Big’95 debris flow, western Mediterranean Sea. J Sediment Res 75:784–797
Liu PLF, Woo SB, Cho YS (1998) Computer program for tsunami propagation and inundation. Technical report; Cornell University: Ithaca, NY, USA
Ma G, Shi F, Kirby J (2012) Shock-capturing non-hydrostatic model for fully dispersive surface wave processes. Ocean Model 43–44:22–35
Macías J, Vázquez JT, Fernández-Salas LM, González-Vida JM, Bárcenas P, Castro MJ (2015) The Al-Borani submarine landslide and associated tsunami. A modelling approach. Mar Geol 361:79–95
Maeda T, Tsushima H, Furumara T (2016) An effective absorbing boundary condition for linear long-wave and linear dispersive-wave tsunami simulations. Earth Planets Space 68(63):1–14
Maeno F, Imamura F (2011) Tsunami generation by a rapid entrance of pyroclastic flow into the sea during the 1883 Krakatau eruption. Indonesia J Geophys Res 116:B09205
Muhari A, Heidarzadeh M, Susmoro H, Nugrohh HD, Kriwati E, Suparto WAB, Imamura F, Arikawa T (2019) The December 2018 Anak Krakatau volcano tsunami as inferred from post-tsunami field surveys and spectral analysis. Pure Appl Geophys 176(12):5219–5233
Nakata K, Katsumata A, Muhari A (2020) Submarine landslide source models consistent with multiple records of the 2018 Palu tsunami, Sulawesi, Indonesia. Earth Planets Space 72:1–16
Pakoksung K, Suppasri A, Imamura F (2018) Systematic evaluation of different infrastructure system for tsunami defense in Sendai City. Geoscience 8:173
Pakoksung P, Suppasri A, Imamura F, Athanasius C, Omang A, Muhari A (2019) Simulation of the submarine landslide tsunami on 20 September 2018 in Palu Bay, Sulawesi Island, Indonesia, using a two-layer model. Pure Appl Geophys 176(8):3323
Paris A, Heinrich P, Paris R, Abadie S (2020) The December 22, 2018 Anak Krakatau, Indonesia, landslide and tsunami: preliminary modeling results. Pure Appl Geophys 177:571–590
Patra AK, Bauer AC, Nichita CC, Pitman EB, Sheridan MF, Bursik M, Rupp B, Webber A, Stinton AJ, Namikawa LM, Renschler CS (2005) Parallel adaptive numerical simulation of dry avalanches over natural terrain. J Volcanol Geotherm Res 139(1–2):1–21
Pitman EB, Nichita CC, Patra A, Bauer A, Sheridan M, Bursik M (2003) Computing granular avalanches and landslides. Phys Fluids 15:3638
PUSHIDROSAL-Naval Hydrographic and Oceanographic Center (2019) Pushidrosal assesses bathymetric change and geomorphology of Anak Krakatau Volcano in Bahasa. https://www.pushidrosal.id/berita/6344/PUSHIDROSAL-TELITI-PERUBAHAN-KONTUR-KEDALAMAN-DAN-GEOMORFOLOGI-GUNUNG-ANAK-KRAKATAU/. Accessed 20 Mar 2019
Satake K (1995) Linear and nonlinear computations of the 1992 Nicaragua earthquake tsunami. Pure Appl Geophys 144:455–470
Shi F, Kirby J, Harris J, Geiman J, Grilli S (2012) A high-order adaptive time-stepping TVD solver for Boussinesq modeling of breaking waves and coastal inundation. Ocean Model 43–44:36–51
Shigihara Y, Goto D, Imamura F, Kitamura Y, Matsubara T, Takaoka K (2006) Hydraulic and numerical study on the generation of a subaqueous landslide-induced tsunami along the coast. Nat Hazards 39:159–177
Skvortsov A (2002) Numerical simulation of landslide generated tsunamis with application to the 1975 failure in Kitimat Arm, British Columbia, Canada. PhD Thesis, School of Earth and Ocean Sciences, University of Victoria, British Columbia.
Storn R, Price K (1997) Differential evolution—a simple and efficient heuristic for global optimization over continuous spaces. J Glob Optim 11:341–359
Suppasri A, Koshimura S, Imamura F (2011) Developing tsunami fragility curves based on the satellite remote sensing and the numerical modeling of the 2004 Indian Ocean tsunami in Thailand. Nat Hazards Earth Syst Sci 11:173–189
Syamsidik B, Luthfi M, Suppasri A, Comfort LK (2020) The 22 December 2018 Mount Anak Krakatau Volcanogenic Tsunami on Sunda Strait Coasts, Indonesia: tsunami and damage characteristics. Nat Hazards Earth Syst Sci 20:549–565
Tinti S, Pagnoni G, Zaniboni F (2006) The landslides and tsunamis of the 30th of December 2002 in Stromboli analysed through numerical simulations. Bull Volcanol 68(5):462–479
Tinti S, Chiocci FL, Zaniboni F, Pagnoni G, de Alteriis G (2011) Numerical simulation of the tsunami generated by a past catastrophic landslide on the volcanic island of Ischia, Italy. Mar Geophys Res 32:287–297
Titan2D (2016) Titan2D mass-flow simulation tool. https ://github.com/TITAN2D/titan 2d. Accessed 6 June 2019
Titov VV, González FI (1997) Implementation and testing of the method of splitting tsunami (MOST) model, NOAA Technical Memorandum ERL PMEL-112, p. 11
Wang X, Liu PF (2006) An analysis of 2004 Sumatra earthquake fault plane mechanisms and Indian Ocean tsunami. J Hydraul Res 44:147–154
Watts P, Grilli ST, Kirby JT, Fryer GJ, Tappin R (2003) Landslide tsunami case studies using a Boussinesq model and a fully nonlinear tsunami generation model. Nat Hazards Earth Syst Sci 3:391–402
Ye L, Kanamori H, Rivera L, Lay T, Zhou Y, Sianipar D, Satake K (2020) The 22 December 2018 tsunami from flank collapse of Anak Krakatau volcano during eruption. Sci Adv 6(3):1377