DETERMINATION OF THE INITIAL SEA SURFACE DEFORMATION CAUSED BY A TSUNAMI USING GENETIC ALGORITHM

Iyan Eka Mulia,Toshiyuki Asano

doi:10.9753/icce.v34.currents.10

Abstract

The characteristics of a tsunami is largely determined by its initial source, which can be defined as a perturbation to the equilibrium state of the ocean hydrodynamic. For tsunamis generated by earthquakes, the initial source (as a form of sea surface deformation) is driven by a coseismic vertical displacement of the seabed. Assuming that the rupture process occurs instantaneously, the sea surface deformation can be inferred from recorded tsunami waveforms. Tsunami waveform inversion using the Green’s function technique with least squares optimization is the most widely used method. However, this method can sometimes suffer from ill-posed problems, because of the non-uniqueness of the inversion results. This study proposes a genetic algorithm that helps the traditional least squares method find an optimum solution. We applied the method to an artificial tsunami source with a complex profile. We achieved significant improvements by determining the optimum spatial distribution of the unit sources (the Green’s function) prior to the inversion. In contrast with the regular Green’s function that has equidistant unit sources, our method generates a random spatial distribution. This random characteristic leads to a better approximation of the initial sea surface deformation.

Highlights

A conservative approach for simulating tsunami sources is to use a simplified elastic dislocation model based on an analytical formula (Okada 1985)
We propose a new approach that helps the least squares method (LSM) find stable and optimal solutions, by determining the optimal position or spatial distribution for the unit sources located around the tsunami source or epicenter
We conducted numerical experiments to determine the initial sea surface deformation caused by a tsunami

Summary

Introduction

A conservative approach for simulating tsunami sources is to use a simplified elastic dislocation model based on an analytical formula (Okada 1985). This method requires an assumption of a uniform slip distribution throughout the rupture area, so the rupture complexity cannot be well represented. A more realistic model of a tsunami source uses tsunami waveform inversion. It considers that the recorded waveforms exhibit the signature of the initial profile of the tsunami, over a specific time period (Satake 1987). We can estimate the heterogeneous slip distribution observed in nature

Methods

Results

Conclusion