Probabilistic Assessment of Tsunami Recurrence in the Indian Ocean

The Indian Ocean is one of the most tsunamigenic regions of the world and recently experienced a mega-tsunami in the Sumatra region on 26 December 2004 (MW 9.2 earthquake) with tsunami intensity I (Soloviev-Imamura intensity scale) equal to 4.5, causing heavy destruction of lives and property in the Indian Ocean rim countries. In this study, probabilities of occurrences of large tsunamis with tsunami intensities I ≥ 2.0 and I ≥ 3.0 (average wave heights H ≥ 2.83 m and H ≥ 5.66 m, respectively) during a specified time interval were calculated using three stochastic models, namely, Weibull, gamma and lognormal. Tsunami recurrence was calculated for the whole Indian Ocean and the special case of the Andaman-Sumatra-Java (ASJ) region, excluding the 1945 Makran event from the main data set. For this purpose, a reliable, homogeneous and complete tsunami catalogue with I ≥ 2.0 during the period 1797–2006 was used. The tsunami hazard parameters were estimated using the method of maximum likelihood. The logarithm of likelihood function (ln L) was estimated and used to test the suitability of models in the examined region. The Weibull model was observed to be the most suitable model to estimate tsunami recurrence in the region. The sample mean intervals of occurrences of tsunamis with intensity I ≥ 2.0 and I ≥ 3.0 were calculated for the observed data as well as for the Weibull, gamma and lognormal models. The estimated cumulative and conditional probabilities in the whole Indian Ocean region show recurrence periods of about 27–30 years (2033–2036) and 35–36 years (2039–2040) for tsunami intensities I ≥ 2.0 and I ≥ 3.0, respectively, while it is about 31–35 years (2037–2041) and 41–42 years (2045–2046) for a tsunami of intensity I ≥ 2.0 and I ≥ 3.0, respectively, in the ASJ region. A high probability (>0.9) of occurrence of large tsunamis with I ≥ 2.0 in the next 30–40 years in the Indian Ocean region was revealed.