Modeling the Tsunami of 25 October 2010 in Sumatra
Disaster Control Research Center, Tohoku University

*

*

Useful Information

General

CNN
Reuters

Emergency Response

ReliefWeb
OCHA Situation Reports

Seismology

USGS

Tsunami Observation

WCATWC web site

*

Modeling the Tsunami of 25 October 2010 in Sumatra
DCRC, Tohoku University

Tsunami Model Description

Geverning Equation : Linear Shallow Water Equations (Far-field model)
Numerical Scheme : Leap-frog Finite Difference Method (TUNAMI-CODE of Tohoku University)
Spatial Grid Size : 1 arc-min. (Far-field model)
Bathymetry Data : GEBCO
Population Data : Landscan

*

Historical events

Historical Earthquake and tsunami in the west coast of Sumatra. White dotted lines are the historical faults, indicated by year in italics (Natawidjaja et al., 2006); transparent rupture zones are the result of recent seismic activity, indicated by year in bold (Sieh et al., 2008), and the seismic gap indicated by the yellow to red area (Natawidjaja, 2009)

*

Tsunami Source Model

Case 1, Case 2 and Case 3

The differences among cases are the location of the fault origin.

Mw = 7.8
Fault Length / Width : 95 km / 48 km
Dislocation : 4.9 m
(Wells D.L. & Coppersmith K.J. 1994.)
Source Mechanism (Strike, Dip, Slip) = (319, 7, 98) Reference :Global CMT
Depth : 12 km

*

Seismic Deformation Model Result

For original resolution, click the image below.

Case 1


*

Case 2


*

Case 3


*

Tsunami Model Result

Modeled tsunami height , click the image below

Case 1


*

Case 2


*

Case 3


*

Tsunami Propagation time (min)

Tsunami Propagation time (case1 and case 2), click the image below

Case 1


*

Case 2


*

Case 3


*

Tsunami Model Result and exposed population

Modeled tsunami height (case 3) and Global (Large scale) Population data, click the image below


Modeled tsunami height (case 3) and Village (Detail) Population data, click the image below

The maximum wave height is calculated using a 90m grid size, based on SRTM data for topography. The model using a fully non-linear shallow water equation solved by a finite difference scheme in TUNAMI-N2. The black bold text showing the name of village, total population and number of casualties on each village. The maximum wave height distribution is shown by the color blue to red.


Potential Tsunami Exposure Estimation

Accumulated number of population at risk (PTE) and estimated number of death against the maximum tsunami height

Number of population at risk along the coastal grid (30 sec or 928 m) was estimated over the simulated maximum tsunami height in order to evaluate a number of potential tsunami exposure (PTE), as exposed population. Results are shown in the Figure 1 below. In brief, there were about 14,000 people at risk of the maximum tsunami height of almost 6 m. Number of death can be estimated by two proposed functions. Tsunami death ratio was first proposed by Kawata (1997) using data from fives historical tsunami events in Japan. Data in Banda Aceh, Indonesia was then added by Koshimura et al. (2009) as shown in Figure 2 (Red dot). Using the curve of Banda Aceh (upper regression curve), number of death is estimated at 534. Moreover, number of death can also be estimated using the proposed fragility curve (Koshimura et al., 2009) as shown in Figure 3. Applying Figure 3, 205 deaths is estimated. Main different of these two curves is that, for example, the death ratio is about 1 % at 1 m tsunami height in Figure 2 where the death ratio starts if tsunami height larger than 4 m in Figure 3.

Figure 1


Death rate as a function of tsunami height (Kawata, 1997 and Koshimura et al., 2009)

Figure 2


Figure 3

Fragility curve for fatality (Koshimura et al., 2009)


References

1. Kawata Y. [1997] “Prediction of loss of human lives due to catastrophic earthquake disaster,” Japan society for natural disaster science, 16(1), 313. (In Japanese)

2. Koshimura, S., Namegaya, Y. and Yanagisawa, H.: Tsunami Fragility - A new measure to assess tsunami damage, Journal of Disaster Research, 4, 479−488, 2009

*

Tsunami Waveforms

For original figure (.png), click the image below.

Tidal station off the western coast of Sumatra


Waveforms at Padang station (Case 2)


Waveforms at Padang station (Case 3)


Mid-ocean Tsunami Propagation Model Animation

Click to play(1.8MB)

Case 2


The verification of the numerical model of 2010 Mentawai tsunami with the field survey data (Case 3)


Tsunami run-up Animation

Click to play(26.37MB)

Case 3 a (Non - Linear Model )

Click to play(17.75MB)

Case 3 b (Non - Linear Model )

Click to play(28.02MB)

Case 3 c (Non - Linear Model )