Numerical modeling of tsunami waves generated by the flank collapse of the Cumbre Vieja Volcano (La Palma, Canary Islands): Tsunami source and near field effects

[1] In this work, we study waves generated by the potential collapse of the west flank of the Cumbre Vieja Volcano (CVV; La Palma, Canary Island, Spain) through numerical simulations performed in two stages: (i) the initial slide motion and resulting free surface elevation are first calculated using a 3D Navier-Stokes model; (ii) generated waves are then input into a 2D (horizontal) Boussinesq model to further simulate propagation to the nearby islands. Unlike in earlier work on CVV, besides a similar extreme slide volume scenario of 450 km3, in our simulations: (i) we consider several slide scenarios featuring different volumes (i.e., 20, 40, 80 km3), which partly result from a geotechnical slope stability analysis; (ii) we use a more accurate bathymetry; and (iii) an incompressible version of a multiple-fluid/material Navier-Stokes model. We find wave trains for each scenario share common features in terms of wave directivity, frequency, and time evolution, but maximum elevations near CVV significantly differ, ranging from 600 to 1200 m (for increasing slide volume). Additionally, our computations show that significant energy transfer from slide to waves only lasts for a short duration (order 200 s), which justifies concentrating our best modeling efforts on the early slide motion phase. The anticipated consequences of such wave trains on La Palma and other Canary Islands are assessed in detail in the paper.

[1]  Stephan T. Grilli,et al.  Numerical Simulation of the 2011 Tohoku Tsunami Based on a New Transient FEM Co-seismic Source: Comparison to Far- and Near-Field Observations , 2013, Pure and Applied Geophysics.

[2]  V. Titov,et al.  A nested-grid Boussinesq-type approach to modelling dispersive propagation and runup of landslide-generated tsunamis , 2011 .

[3]  K. Kelfoun,et al.  Numerical modelling of the tsunami triggered by the Güìmar debris avalanche, Tenerife (Canary Islands): Comparison with field-based data , 2011 .

[4]  F. Montagna,et al.  3D numerical modeling of landslide-generated tsunamis around a conical island , 2011 .

[5]  K. Kelfoun,et al.  Landslide-generated tsunamis at Réunion Island , 2010 .

[6]  G. Bellotti,et al.  Inclusion of landslide tsunamis generation into a depth integrated wave model , 2010 .

[7]  J. Kirby,et al.  Numerical simulation and first-order hazard analysis of large co-seismic tsunamis generated in the Puerto Rico trench: near-field impact on the North shore of Puerto Rico and far-field impact on the US East Coast , 2010 .

[8]  Stephan T. Grilli,et al.  Numerical simulation of waves generated by landslides using a multiple-fluid Navier–Stokes model , 2010 .

[9]  Hongqiang Zhou,et al.  Extended Fourth-Order Depth-Integrated Model for Water Waves and Currents Generated by Submarine Landslides , 2010 .

[10]  Stéphane Abadie,et al.  Vague générée par un glissement de terrain influence de la forme initiale et de la déformabilité du glissement , 2010 .

[11]  M. Ioualalen,et al.  Reconstructions of the coastal impact of the 2004 Indian Ocean tsunami in the Khao Lak area, Thailand , 2009 .

[12]  Per A. Madsen,et al.  Tsunami generation, propagation, and run-up with a high-order Boussinesq model , 2009 .

[13]  E. Tric,et al.  3D NUMERICAL SIMULATIONS OF WAVES GENERATED BY SUBAERIAL MASS FAILURES: APPLICATION TO LA PALMA CASE , 2009 .

[14]  Hermann M. Fritz,et al.  Hybrid modeling of the mega‐tsunami runup in Lituya Bay after half a century , 2009 .

[15]  W. Ryan,et al.  Global Multi‐Resolution Topography synthesis , 2009 .

[16]  Finn Løvholt,et al.  Oceanic propagation of a potential tsunami from the La Palma Island , 2008 .

[17]  S. Grilli,et al.  The Papua New Guinea tsunami of 17 July 1998: anatomy of a catastrophic event , 2008 .

[18]  Stephan T. Grilli,et al.  Experimental Study of Tsunami Generation by Three-Dimensional Rigid Underwater Landslides , 2007 .

[19]  Stephan T. Grilli,et al.  Source Constraints and Model Simulation of the December 26 , 2007 .

[20]  J. Kirby,et al.  Modeling the 26 December 2004 Indian ocean tsunami : Case study of impact in Thailand - art. no. C07024 , 2007 .

[21]  P. Sedwick,et al.  Elevated marine deposits in Bermuda record a late Quaternary megatsunami , 2006 .

[22]  Pierre Lubin,et al.  Three-dimensional Large Eddy Simulation of air entrainment under plunging breaking waves , 2006 .

[23]  J. Carracedo,et al.  Tsunami deposits related to flank collapse in oceanic volcanoes : The Agaete Valley evidence, Gran Canaria, Canary Islands , 2006 .

[24]  Qin Chen,et al.  Fully Nonlinear Boussinesq-Type Equations for Waves and Currents over Porous Beds , 2006 .

[25]  K. Kelfoun,et al.  Numerical modeling of the emplacement of Socompa rock avalanche, Chile , 2005 .

[26]  Stephan T. Grilli,et al.  Tsunami Generation by Submarine Mass Failure. I: Modeling, Experimental Validation, and Sensitivity Analyses , 2005 .

[27]  G. Fryer,et al.  Tsunami Generation by Submarine Mass Failure. II: Predictive Equations and Case Studies , 2005 .

[28]  Stefano Tinti,et al.  The 30 December 2002 landslide-induced tsunamis in Stromboli: sequence of the events reconstructed from the eyewitness accounts , 2005 .

[29]  Costas E. Synolakis,et al.  Runup and rundown generated by three-dimensional sliding masses , 2005, Journal of Fluid Mechanics.

[30]  Stephan T. Grilli,et al.  Mechanical models of the 1975 Kalapana, Hawaii earthquake and tsunami , 2005 .

[31]  Willi H. Hager,et al.  Near Field Characteristics of Landslide Generated Impulse Waves , 2004 .

[32]  J. Lénat,et al.  Recurrence of major flank landslides during the last 2-Ma-history of Reunion Island , 2004 .

[33]  M. Quecedo,et al.  Numerical modelling of impulse wave generated by fast landslides , 2004 .

[34]  B. Eakins,et al.  Calculated volumes of individual shield volcanoes at the young end of the Hawaiian Ridge , 2003 .

[35]  Robert A. Dalrymple,et al.  Boussinesq modeling of longshore currents , 2003 .

[36]  Stephan T. Grilli,et al.  Landslide tsunami case studies using a Boussinesq model and a fully nonlinear tsunami generation model , 2003 .

[37]  P. Liu,et al.  A numerical study of submarine–landslide–generated waves and run–up , 2002, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[38]  Stephan T. Grilli,et al.  Development of a 3D numerical wave tank for modeling tsunami generation by underwater landslides , 2002 .

[39]  F. Legros.  The mobility of long-runout landslides , 2002 .

[40]  F. Imamura,et al.  Giant landslides, mega-tsunamis, and paleo-sea level in the Hawaiian Islands , 2001 .

[41]  S. Day,et al.  Cumbre Vieja Volcano—Potential collapse and tsunami at La Palma, Canary Islands , 2001 .

[42]  I. Kimio.  Shimabara-Shigatusaku earthquake and topographic changes by Shimabara catastrophe in 1792 , 2000 .

[43]  Philip Watts,et al.  Tsunami features of solid block underwater landslides , 2000 .

[44]  Simon Day,et al.  Later stages of volcanic evolution of La Palma, Canary Islands: Rift evolution, giant landslides, and the genesis of the Caldera de Taburiente , 1999 .

[45]  J. Caltagirone,et al.  Mécanisme de génération du jet secondaire dans un déferlement plongeant , 1998 .

[46]  D. Masson,et al.  A giant landslide on the north flank of Tenerife, Canary Islands , 1995 .

[47]  Oldrich Hungr,et al.  A model for the runout analysis of rapid flow slides, debris flows, and avalanches , 1995 .

[48]  G. Wei,et al.  A fully nonlinear Boussinesq model for surface waves. Part 1. Highly nonlinear unsteady waves , 1995, Journal of Fluid Mechanics.

[49]  Lin Jiang,et al.  Numerical modeling of an underwater Bingham plastic mudslide and the waves which it generates , 1993 .

[50]  P. Leblond,et al.  The coupling of a submarine slide and the surface waves which it generates , 1992 .

[51]  Carl B. Harbitz,et al.  Model simulations of tsunamis generated by the Storegga Slides , 1992 .

[52]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[53]  Stephan T. Grilli,et al.  A high-order adaptive time-stepping TVD solver for Boussinesq modeling of breaking waves and coastal inundation , 2012 .

[54]  Jeffrey C. Harris,et al.  Numerical Simulation of Tsunami Generation by the Potential Flank Collapse of the Cumbre Vieja Volcano , 2011 .

[55]  Stephan T. Grilli,et al.  tsunami: Case study of impact in Thailand , 2007 .

[56]  Michael L. Gittings,et al.  SAGE CALCULATIONS OF THE TSUNAMI THREAT FROM LA PALMA , 2006 .

[57]  J. Carracedo,et al.  Massive flank failures and tsunamis in the Canary Islands: past, present, future , 2003 .

[58]  R. Wynn,et al.  Canary Islands Landslides and Tsunami Generation: Can We Use Turbidite Deposits to Interpret Landslide Processes? , 2003 .

[59]  James T. Kirby,et al.  Chapter 1 Boussinesq models and applications to nearshore wave propagation, surf zone processes and wave-induced currents , 2003 .

[60]  Richard E. Thomson,et al.  Numerical Modeling of Tsunami Generation by Submarine and Subaerial Landslides , 2003 .

[61]  Neil C. Mitchell,et al.  Slope failures on the flanks of the western Canary Islands , 2002 .

[62]  C. Mader.  MODELING THE LA PALMA LANDSLIDE TSUNAMI , 2001 .

[63]  J. Kirby,et al.  BOUSSINESQ MODELING OF WAVE TRANSFORMATION, BREAKING, AND RUNUP. II: 2D , 2000 .

[64]  Stephan T. Grilli,et al.  PII: S0955-7997(99)00021-1 , 1999 .

[65]  J. Kirby,et al.  BOUSSINESQ MODELING OF WAVE TRANSFORMATION , BREAKING , AND RUNUP . II : 2 D By , 1999 .

[66]  Robin T. Holcomb,et al.  Large landslides from oceanic volcanoes , 1991 .

[67]  Pierre Cochonat,et al.  Importance des dépo gravitaires dans la mise en place d'un système volcano-sédimentaire sous-marin (Volcan de la Fournaise, Ile de la Réunion) , 1990 .

[68]  R. LeVeque,et al.  Numerical methods for conservation laws , 1990 .

[69]  William R. Normark,et al.  Prodigious submarine landslides on the Hawaiian Ridge , 1989 .

[70]  L YoungsD,et al.  Time-dependent multi-material flow with large fluid distortion. , 1982 .

[71]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[72]  K. Goda,et al.  A multistep technique with implicit difference schemes for calculating two- or three-dimensional cavity flows , 1979 .

[73]  J. Caltagirone,et al.  International Journal for Numerical Methods in Fluids Efficient Solving Method for Unsteady Incompressible Interfacial Flow Problems , 2022 .