Tsunami-HySEA model background and validation

HySEA codes have been developed by EDANYA Group (https://edanya.uma.es) from the Universidad de Málaga (UMA) for more than a decade and they are in continuous evolution and upgrading. Initially, the software was developed and the numerical algorithms implemented published under no particular name. Several developments have been published in peer-reviewed international journals since 2005 where different analytical and experimental test cases have been presented. Ordered by year of publication, some of these model developments can be found in [1]-[7]. In September 2013, at ITS 2013, held at Göcek (Turkey) the ensemble of all the geophysical codes developed by the EDANYA group were named under the common name of HySEA for the very first time, and Landslide-HySEA (for a case study of an aerial landslide) and Tsunami-HySEA were presented in the two separate contributions [8] and [9], respectively. The work recently published in Marine Geology in [10], in this case for a submarine landslide, is the first peer-reviewed paper where a HySEA code is named as such.

Tsunami-HySEA combines robustness, reliability and good accuracy in a model based on a GPU faster than real time (FTRT) implementation. It has been severely tested, and in particular has passed all tests in Synolakis et al. (2008), but also other laboratory tests and proposed benchmark problems. Some of them can be found in references [1]-[7] cited above. Synolakis et al. (2008) test cases for Tsunami-HySEA can be found at [11] in Spanish. Benchmark problems 1, 4, 6, 7, and 9 can be found in the report submitted to the MMS/NTHMP [12]. Among these benchmarks the Monai Valley test case, one of the most emblematic test cases for tsunami codes, was presented at EGU 2013 in [13]. A comparison of Tsunami-HySEA numerical results with the MOST model for LANTEX 2013 scenario and tsunami impact on Puerto Rico coasts can be found at [14]. Tsunami-HySEA has also participated in the benchmarking exercise “NOAA/NTHMP MMS Benchmarking Workshop: Tsunami Currents”, see [15] in NTHMP (2016). Besides all this, much effort has been made to develop a specific Tsunami-HySEA code suitable for tsunami computations in the framework of Tsunami Early Warning Systems (TEWS). The work done to develop a tsunami code for computations in times much faster than real time was presented at EGU 2014 in [16] and at Perspectives of GPU Computing in Physics and Astrophysics in [17]. The TEWS version of Tsunami-HySEA is currently the core numerical code at INGV (Instituto Nazionale di Geofisica e Vulcanologia) of the Italian TEWS. It has also been adopted by the JRC (Joint Research Centre of the European Commission) as one of the numerical codes of their Tsunami Alert System and it will be used as computing code by the Spanish IGN (Instituto Geográfico Nacional) in its TEWS. In order to be approved by the INGV as the numerical code for their system, the Tsunami-HySEA model has passed through an exhaustive one-year validation process at INGV. In fact, sixteen different TVD (Total Variation Diminishing) and non-TVD numerical schemes were implemented and tested in search for robustness, computational speed and suitable numerical results, in order to make a final suitable choice for the code to be adopted by the Italian TEWS. A comparative study containing most of these numerical schemes can be found at [18]. An important application of the Tsunami-HySEA model to the emblematic real case of Tohoku 2011, can be found in [19]. In December 2014, Lisbon 1775 and Puerto Rico 1918 earthquake scenarios were simulated with Tsunami-HySEA and presented at the “Experts Meeting Workshop on Tsunami Modeling and Mitigation” held at Cartagena de Indias (Colombia). The corresponding presentation can be downloaded at [20].

See References

 

 Model Equations

 Tsunami HySEA