Propagation of Tsunamis Generated by Seabed Motion with Time-History and Spatial-Distribution: An Analytical Approach

Taehwa Jung; Sangyoung Son

doi:10.9765/KSCOE.2018.30.6.263

Journal of Korean Society of Coastal and Ocean Engineers > Volume 30(6); 2018 > Article

시간이력 및 공간분포를 지닌 지반운동에 의한 지진해일 발생 및 전파: 해석적 접근

Article

J Korean Soc Coast Ocean Eng 2018; 30(6): 263-269.

Published online: December 31, 2018

DOI: https://doi.org/10.9765/KSCOE.2018.30.6.263

시간이력 및 공간분포를 지닌 지반운동에 의한 지진해일 발생 및 전파: 해석적 접근

정태화^*, 손상영^**

^*한밭대학교 건설환경공학과

^**고려대학교 건축사회환경공학부

Propagation of Tsunamis Generated by Seabed Motion with Time-History and Spatial-Distribution: An Analytical Approach

Taehwa Jung^*, Sangyoung Son^**

^*Department of Civil and Environmental Engineering, Hanbat National University

^**School of Civil, Environmental, and Architectural Engineering, Korea University 145 Anam-ro, Seongbuk-gu, Seoul 02841, Korea

Corresponding author: Sangyoung Son, School of Civil, Environmental, and Architectural Engineering, Korea University 145 Anam-ro, Seongbuk-gu, Seoul 02841, Korea, Tel: +82-2-3290-4865, Fax: + 82-2-3290-5999, sson@korea.ac.kr

Received November 5, 2018 Revised December 6, 2018 Accepted December 12, 2018

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

요약

해저에서 지진 및 산사태 등에 기인하는 수심변화는 자유수면의 변동을 유발한다. 이러한 자유수면의 변동은 해안으로 전파하고 천수현상에 의해 파고가 상승하면서 큰 피해를 야기한다. 지진해일의 생성 및 전파과정을 모의하기 위해 다양한 수치해석모형이 개발된 바 있다. 대부분의 지진해일모형은 해저 지반의 움직임이 수면에 즉시적이고 동일하게 전달된다는 가정하에 초기수면을 결정한다. 하지만, 이 접근 방법은 시간적 및 공간적 차이를 지니고 발생하는 해저지진의 특성을 고려하지 못하며, 이는 그대로 해석오차에 반영된다. 본 연구에서는 Hammack(1973)의 해석해를 응용하여 수심이 시간차를 가지고 변할 때 이로 인한 초기 자유수면의 변동과 이렇게 변화된 수면의 전파 특성을 모의하고 그 결과를 검토하였다. 이는 다양한 유형의 지반운동과 자유수면 형성과의 인과관계 및 상승된 자유수면의 이동속도 등의 관계규명에 도움이 될 것이다.

핵심용어: 시간이력 수심변화, 지진해일, 해석해, 지진해일모형, 초기 파형, Laplace-Fourier 변환

Abstract

Changes in water depth caused by underwater earthquakes and landslides cause sea surface undulations, which in turn propagate to the coast and result in significant damage as wave heights normally increase due to the wave shoaling process. Various types of numerical models have been developed to simulate the generation and propagation of tsunami waves. Most of tsunami models determine the initial surface of the water based on the assumption that the movement of the seabed is immediately and identically transmitted to the sea surface. However, this approach does not take into account the characteristics of underwater earthquakes that occur with time history and spatial variation. Thus, such an incomplete description on the initial generation of tsunami waves is totally reflected in the error during the simulation. In this study, the analytical solution proposed by Hammack (1973) was applied in the tsunami model in order to simulate the generation of initial water surface elevation by the change of water depth with time history and its propagation. The developed solution is expected to identify the relationship among various type of seabed motions, initial surface undulations, and wave speeds of elevated water surfaces.

Keywords: time-delayed bottom variation, tsunami, analytical solution, tsunami model, initial surface elevation, Laplace-Fourier transformation

1. 서 론

해저에서의 지각변동 또는 산사태는 수심의 변동을 야기하며 이러한 수심의 변화는 자유수면으로 전달되어 지진해일이라고 일컫는 수면변화를 일으킨다. 상승된 수면은 천수 및 굴절 현상 등을 겪으면서 해안으로 전파하고 이에 의해 증가된 파고는 연안 지역에 막대한 피해를 야기한다.

지진해일의 발생과 그로 인한 피해는 역사적으로 오랜기간 동안 관측 및 기록되어 왔다(Braynt, 2014; Gusiakov, 2009). 한 연구를 통해 보고된 바에 따르면 지구상 20여개국의 100여개 도시가 지진해일의 위험성에 심각하게 노출되어 있다고 한다(Wilson et al., 2013). 2000년대에 발생한 지진해일의 대표적 사례를 소개하면 다음과 같다. 우선 2004년 12월 26일 인도네시아 수마트라 섬 서부 해안의 40 km 지점에서 발생한 규모 9.3의 인도양 수마트라 지진에 의해 약 33 m에 이르는 지진해일고가 해안에 전파되었으며, 이로 인해 약 23만명이 목숨을 잃고 5만 명 이상이 실종, 169만 명 이상의 난민이 발생하는 등의 인명피해를 입었다. 이는 역사상 3번째로 큰 규모의 대지진으로 기록되기도 하였다(Grilli et al., 2007). 한편, 2011년 3월 11일에 일본 동북부 태평양 연안에서 발생한 규모 9.0 도호쿠-오키 지진해일은 38.9 m의 최대 지진해일고를 기록함에 따라 막대한 인명 및 재산피해를 야기했다(Bryant, 2014). 약 2만 명의 사망 및 실종자와 6천 여명 부상자가 발생하였으며 총 피해액이 330조 원에 이르는 등 그 피해가 천문학적이었다. 가장 최근인 2018년 9월 28일에 인도네시아 팔루시에서 약 77 km 떨어진 곳에서 일어난 규모 7.5의 술라웨시 지진은 약 2100명의 사망피해를 낳기도 하였다. 한 가지 주목할 만한 점은 술라웨시 지진의 경우 섭입대(subduction zone)에서 주로 발생하는 지진과 달리 지진해일이 발생하기 힘들다고 알려진 주향이동단층(strike-slip, 수평 이동하는 단층)이었음에도 불구하고 약 4~7 m의 지진해일고를 기록하였다는 점이다. 그 원인에 대해서는 만의 형태를 지닌 피해지역의 지형, 해저산사태의 연쇄작용 등 현재 과학계에서 다양한 주장이 제기되고 있는 실정이며 정확한 원인규명을 위한 후속연구가 진행 중이다.

이처럼 지진해일의 피해는 불의 고리라고 일컫는 환태평양 조산대를 중심으로 하여 여러 나라에서 지속적으로 발생하고 있으며, 그 피해 역시 엄청난 규모를 나타내기에 정확한 예측모형개발과 방재시스템 구축에 많은 노력을 기울이고 있다. 지진해일을 모의 할 때, 통상적으로 장파이론에 기반을 둔선형 혹은 비선형의 비분산성 천수방정식을 사용하게 되고, 이에 따라 다양한 모형들이 구축되어 왔다. 몇 가지 예를 들면, Titov and Synolakis(1998)가 개발한 MOST(Method Of Splitting Tsunami), Wang(2009)이 개발한 COMCOT(COrnell Multi-grid COupled Tsunami model), Kowalik et al.(2005)가 개발한 GTM(Global Tsunami Model), George and LeVeque(2006)이 만든 TsunamiClaw(Conservation Laws) 등이 있다. 한편 지진해일 등의 장파해석에 Boussinesq 방정식 혹은 Navier-Stokes 방정식을 이용하는 등 파랑의 분산성과 동수압을 고려함으로써 정확도 높은 해석결과를 얻으려는 노력도 병행되고 있다(Son et al., 2011; Grilli et al., 2007; Saito and Furumura, 2009; Pringle et al., 2018).

이상에서 소개한 지진해일과 관련된 지금까지의 연구는 해저 수심의 변화가 자유 수면에 즉시적이고 동일하게 전달된다는 가정하에 초기파형을 결정하고 파랑의 전파양상을 수치모의하는 방식이 대부분이다. 하지만 최근 관측장비의 발전과 분석기법의 개발을 통해 지진발생의 세밀한 과정이 규명됨에 따라 지진시 해저지반의 운동이 시간이력과 공간분포 특성을 지니고 있음이 알려졌고, 이에 따라 지진해일 또한 이를 반영하여 생성된다고 밝혀졌다. Suppasri et al.(2010)과 Ohmachi et al.(2001)은 파열 속도(rupture velocity)가 느리고 단층에 가까운 지역일수록 시간이력을 고려함에 따른 지진해일의 차이가 증가함을 확인하였다. 또한 Imai et al.(2010)은 수분의 시간차를 두고 발생하는 대규모 지진의 경우 단층운동의 시간이력에 따라 지진해일이 해안에서 증폭될 수 있음을 보였다. Ishii et al.(2005) 및 Suzuki et al.(2011)은 각각 2004년 수마트라지진 및 2011년 동일본 대지진에 대하여 단층파괴과정의 시간이력을 실제적으로 확인하였으며 이에 따라 시간이력의 중요성을 강조하기도 하였다. 더군다나 기존의 접근법은 상대수심(파장에 대한 수심비)이 매우 작은 경우에만 성립하며 지진해일의 파장이 짧거나 수심변화의 지속시간이 긴 경우에는 적용할 수 없다(Lee et al., 2015; Saito and Furumura, 2009).

수심변화에 의한 초기 수면변위를 예측하는 연구는 주로 해석적인 방법을 통해 이루어졌다. Kajiura(1963)는 Green 함수를 이용하여 바닥의 변화가 근거리 및 원거리에서의 수면변위에 끼치는 영향을 검토하였으며, Hammack(1973)은 속도 포텐셜을 이용하여 수심변화가 자유수면에 끼치는 영향을 해석적으로 연구하였다. 이후 Dutykh et al.(2006), Dutykh and Dias(2007), Kervella et al.(2007) 및 Dutykh et al.(2013) 등은 Hammack(1973)의 접근법을 응용하여 다양한 수심변위 조건에서 수면변위를 계산하였다.

이상의 논문들을 통해서 다양한 수심변화에 의해서 발생하는 수면변화 및 기타 수반되는 여러가지 물리량들에 대해서는 연구가 수행되었으나 바닥수심의 변화가 시간차를 가지고 발생하는 경우에 대해서는 아직까지 연구된 바가 없다. 따라서 본 연구에서는 Hammack(1973)의 해석해를 확장하여 영역별로 시간차를 가지고 바닥수심이 변하는 경우에 대해 해석해를 제시하여 수면변위를 계산하였다. 본 연구에서는 2개의 영역에서 대해서만 다루었으나 이를 다수의 영역으로 확장할 수 있으며 또한 수심변화가 발생하는 시간 또한 다수로 나누어 계산할 수 있다. 이와 같은 접근법은 공간 및 시간 영역을 세밀하게 나누어 계산할 경우, 실제에서 발생하는 바닥 수심의 변화를 좀 더 정확하게 모의할 수 있다는 장점이 있다.

2. 해석해 유도

비압축성 유체의 비회전 흐름을 가정하면, 속도포텐셜을 이용하여 지배방정식을 다음과 같이 표현할 수 있다.

(1)

∇2ϕ=0 in D

여기서 φ는 속도포텐셜을 의미하며 D는 해석 영역을 의미한다. 선형화된 자유수면 경계조건은 다음과 같이 표현이 가능하며,

(2)

∂ϕ∂z=∂η∂t at z=0

(3)

∂ϕ∂t+gh=0 at z=0

이를 합쳐 다음과 같은 단일화된 자유수면 경계조건식을 구할 수 있다. 여기서 η는 자유수면, g는 중력가속도, h는 정수면으로부터 측정된 수심을 의미한다.

(4)

∂2ϕ∂t2+g∂ϕ∂z=0 at z=0

Fig. 1과 같이 시간에 따른 수심의 변화를 고려하면 선형화된 바닥경계조건은 다음과 같이 된다.

(5)

∂ϕ∂z=∂ζ∂t at z=−h

식(1), (4) 및 (5)에 대하여, 시간에 관해서는 Laplace 변환을 적용하고 공간에 대해서 Fourier 변환을 적용하면, 각각 다음과 같이 나타낼 수 있다.

(6)

∂ϕ¯(k,z,s)2∂2z−k2ϕ¯(k,z,s)=0

(7)

∂ϕ¯(k,−h,s)∂z=sζ¯(k,s)

(8)

∂ϕ¯(k,0,s)∂z+s2gϕ¯(k,0,s)=0

여기서 f는 다음과 같이 정의된다.

(9)

f¯(k,s)=∫−∞∞dx∫0∞eikxe−stf(x,t)dt

식(6)~(8)을 이용하여 ϕ를 구한후, Fourier 및 Laplace 역변환을 취하면 수심변화에 의한 자유수면을 다음과 같이 계산할 수 있다(Hammack, 1973; Dutykh and Dias, 2007).

(10)

η(x,t)=12∫−∞∞{limΓ→∞12πi∫μ−iΓμ+iΓs2e−ikxestζ¯(k,s)(s2+ω2)coshkhds}dk

Hammack(1973)은 원점을 중심으로 대칭인 직사각형 형태의 지형이 바닥에서 상승하거나 하강하는 경우에 대한 선형파의 해석해를 제시한 바가 있다. 이를 좀 더 일반화하여, Fig. 1과 같이 특정시간(t₁ →t₂) 동안 특정범위(x₁ ≤ x ≤ x₂)에서 최대변위(ζ₀)가 발생한다고 하면 ζ는 다음과 같이 나타낼 수 있다.

(11)

ζ(x,t)=ζ0×{0,t<t1t−t1t2−t1t1≤t≤t2in x1≤x≤x21,t>t2

또는

(12)

ζ(x,t)=ζ0{H(t−t2)(1−t−t1t2−t1)+H(t−t1)t−t1t2−t1}[H(x−x1)−H(x−x2)]for t≥t1

여기서 H는 단위 계단함수를 의미한다. 식(12)를 식(9)에 대입하여 ζ를 구하면 다음과 같이 되며,

(13)

ζ¯(k,s)=ζ0(sinkx2−sinkx1k−icoskx2−coskx1k)1−e−s(t2−t1)s2(t2−t1)

이를 다시 식(10)에 대입하면 최종적으로 Fig. 1과 같이 주어진 조건에서 바닥지형변화에 따른 수면변위를 구할 수 있다.

(14)

η(x,t)=ζ0π(t2−t1)∫0∞ABdk

여기서 A 및 B는 식(13)에서 공간과 시간에 해당하는 부분으로 실수값으로 각각 다음과 같이 주어진다.

식(14)의 적분값은 해석적으로 구하는 것이 힘들기 때문에 수치적분을 통해 구할 수 있다. 이 경우 Filon(1930)이 제안한 구적법(quadrature formula)를 사용하면 효율적으로 계산할 수 있다.

(15)

A=Re[e−ikxcoshkh{sinkx2−sinkx1k−icoskx2−coskx1k}]=1kcoshkh{coskx(sinkx2−sinkx1)−sinkx(coskx2−coskx1)}

(16)

B=Re[limΓ→∞12πi∫μ−iΓμ+iΓs2ests2+ω2⋅1−e−s(t2−t1)s2ds]=H(t−t1)sinω(t−t1)−H(t−t2)sinω(t−t2)ω

2.1 Hammack(1973)의 해와의 연관성

Hammack(1973)은 |x| ≤ b에서 수심변화가 발생하는 경우에 대하여 다음과 같은 해석해를 제안하였다.

(17)

ηHammack(x,t)=ζ0π∫0∞[2coskxsinkbkcoshkh]{(α2α2+ω2)(−e−at+cosωt+ωαsinωt)}dk

여기서 [ ]에 해당하는 부분이 식(14)의 A와 대응되는 부분이다. x₂ = b 및 x₁ = − b로 놓으면 식(15)는 Hammack(1973)의 해와 동일하게 된다. 반대로 Hammack(1973)의 해석해로부터 식(15)를 유도할 수 있다. Hammack(1973)은 수심변화 구간의 중심을 원점으로 설정하였다. 따라서, 본 연구에서의 좌표체계를 Fig. 2와 같이 x'로 가정하면 다음과 같은 관계식이 성립한다.

(18)

x = x' − (x2' + x1')/2 및 b =(x2' + x1')/2

식(18)을 식(17)의 [ ]에 대입하여 정리하면,

(19)

2coskxsinkbkcoshkh=2kcoskhcos(kx′−k(x2'+x1')2)sin(k(x2'−x1')2)

삼각함수에 대한 합과 곱의 공식을 사용하면

(20)

cos(kx′−k(x2'+x1')2)sin(k(x2'−x1')2)=coskx′(sinkx2'−sinkx1')2−sinkx′(coskx2'−coskx1')2

이 되어 식(19)는 식(15)와 동일하게 된다. 식(14)의 시간에 해당하는 B 부분은 Dutykh and Dias(2007)의 결과와 비교하여 검증할 수 있다. Dutykh and Dias는 초기부터 해저 지형이 움직이는 경우인 경우에 대한 시간항의 해석해를 제시한 바가 있다. 이는 B의 t₁이 0이 되는 경우로서, t₁ =0으로 놓으면 본 연구의 시간항은 그들의 결과와 동일하게 된다.

2.2 두 영역에서 수심변화가 있는 경우

Fig. 3과 같이 두 영역에서 시간차를 가지고 수심변화가 발생하는 경우 수심은 다음과 같이 표현할 수 있다.

(21)

ζ(x,t)=ζ1{H(t−t1,e)(1−t−t1,st1,e−t1,s) +H(t−t1,s)t−t1,st1,e−t1,s}[H(x−x1,s)−H(x−x1,e)] +ζ2{H(t−t2)(1−t−t2,st2,e−t2,s) +H(t−t2)t−t2,st2,e−t2,s}[H(x−x2,s)−H(x−x2,e)]

앞선 과정을 반복하며 주어진 수심변화로 인한 수면변화를 구해보면 다음과 같다.

(22)

η(x,t)=ζ1π(t1,e−t1,s)∫0∞A1B1dk+ζ2π(t2,e−t2,s)∫0∞A2B2dk

여기서 A₁, B₁, A₂ 및 B₂는 각각 다음과 같다.

(23)

A1=1kcoshkh{coskx(sinkx1,e−sinkx1,s)−sinkx(coskx1,e−coskx1,s)}

(24)

B1=H(t−t1,s)sinω(t−t1,s)−H(t−t1,e)sinω(t−t1,e)ω

(25)

A2=1kcoshkh{coskx(sinkx2,e−sinkx2,s)−sinkx(coskx2,e−coskx2,s)}

(26)

B2=H(t−t2,s)sinω(t−t2,s)−H(t−t2,e)sinω(t−t2,e)ω

식(22)는 동일한 방법을 사용하여 n개 영역에서 수심변화가 발생하는 경우로 확장할 수도 있다.

3. 결과 분석

본 연구의 접근법에 대한 검증은 Hammack(1973)이 실험값과 비교함으로써 증명하였기 때문에 생략하였다. 먼저, Fig. 1과 같이 하나의 영역에서 바닥수심이 변하는 경우에 대하여 자유수면 변화를 모의하였다. 계산 조건은 Hammack(1973)과 동일하게 수심(h) 0.5 m, 원점에서 영역의 우측경계까지의 거리(b) 0.61 m 및 바닥수심의 최대 변동폭(ζ₀)은 0.05 m로 설정하였다.

Fig. 4는 주어진 조건에서 바닥이 상승하는 시간(t_rise)에 변화를 주면서 특정지점에서의 수면변위와 특정 시간대에서의 수면분포를 도시한 그림이다. 무차원화된 바닥상승시간(t'=trise·gh/h)의 값은 0.1, 1.0 및 10.0을 사용하였으며 Fig. 4(a)는 원점(x = 0)에서 시간에 따른 수면변위를 계산한 값이며 Fig. 4(b)는 t = t_rise에서 위치에 따른 수면분포를 계산한 결과이다. 수심변화가 비교적 짧은 시간에 발생하는 경우, 수위는 수심이 증가하는 동안 꾸준히 상승하다가 수심 변화가 멈춘 뒤 점진적으로 감소하는 모습을 보였다. 그러나 수심 변화가 느리게 발생하는 경우에는 수심이 상승하는 동안 상승된 수면이 좌우로 이동하면서 0.2~1.0 의 시간대에서는 거의 수면변화가 발생하지 않았다. 이는 Fig. 4(b)를 통해서도 확인할 수 있다. t'의 값이 증가할수록 수면의 진행이 크게 발생하여 10.0인 경우에는 이미 많이 진행을 하여 거의 전 영역에서 비슷한 수면 변위를 보이는 것을 알 수 있다.

다음으로 Fig. 3과 같이 두 개의 영역에서 시간차를 가지고 수심 변화가 발생하는 경우에 대하여 수면변위를 계산해보았다. 계산 조건은 Fig. 4의 경우와 유사하게 설정하였다. 수심은 0.5 m로 설정하였으며, ζ₁ = ζ₂ = ζ₀ = 0.05 m 그리고 x_{1, e} − x_{1, s} = x_{2, e} − x_{2, s} = b = 0.61m로 설정하였다.

Fig. 5는 두 영역이 이웃한 상황 (즉, x_{2, s} = x_{1, e})에서 시간차를 가지고 상승하는 경우에 수면분포를 모의한 결과이다. 수심의 상승시간은 두 영역 모두 t' = 1.0으로 하였으며, 파가 어느 정도 진행을 한 t = 3t_rise인 순간에 위치에 따른 수면변위를 계산하였다. 시간차가 없는 경우(실선)와 비교해보면, t_delay = 2.0t_rise인 경우에는 최대수면변위는 상승하였으며 t_delay = 4.0t_rise인 경우에는 최대수면변위가 오히려 감소하였다. 아울러 최대값이 발생하는 지점 역시 모두 달랐다.

Fig. 6은 다른 조건 등을 동일하게 두고 두 영역의 거리에 변화를 주면서 수면변위를 계산한 결과이다. t' = 1.0 및 t_delay = 2.0t_rise인 상황에서 두 영역 사이의 거리(x_{2, s} − x_{1, e})를 0, 1b 및 2b로 변화를 주면서 값을 계산하였다. Fig. 5의 결과와 유사하게 거리에 따라 최대수면변위의 값과 최대변위가 발생하는 지점이 달라지는 것을 확인할 수 있다.

4. 결 론

본 연구에서는 Hammack(1973)의 해석해를 응용하여 수심 변화가 발생하는 경우에 초기 수면변위의 생성 및 전파 특성을 검토하였다. 수심이 상승하는 시간이 짧으면 수심변화에 의해 상승된 수면이 이동하기 전에 수심 상승이 발생하여 수심의 상승과 수면의 상승이 비례하면서 거의 동시에 발생하였다. 그러나 수심이 상승하는 시간이 긴 경우에는 수심이 상승하기 전에 이미 상승된 수면이 좌우로 이동하여 수심이 상승하는 동안에 더 이상의 수면상승이 발생하지 않았다.

두 개의 영역에서 시간차를 가지며 수심변화가 발생하는 경우에는 두 영역간의 거리차 및 수심이 변하는 시간차에 따라서 최대 수면변위의 값이랑 위치가 달라지는 것을 확인할 수 있다. 이는 먼저 상승한 수심으로 인해 발생한 수면이 이동하면서 두 번째로 상승한 수심으로 인한 수면변위와 중첩하여 발생하는 것으로 경우에 따라 증가하거나 감소하였다.

실제 해역에서 지진해일이 발생하는 경우, 대부분의 지진해일 수치모형은 바닥의 변화를 자유수면의 초기파형으로 가정하기 때문에 이러한 시간차에 의해 발생하는 효과를 고려하지 못하고 있다. 본 연구를 확장하여 바닥변화와 자유수면과의 인과관계 및 상승된 자유수면의 이동속도 등을 규명한다면 향후 지진해일 수치모형연구에 도움이 될 것으로 기대한다.

감사의 글

본 연구는 행정안전부 재난안전기술개발사업단 연구사업(MOIS-재난-2015-03)과 기상청 기상 · 지진 See-At 기술개발연구사업(KMI 2018-09510)의 지원으로 수행되었습니다.

Fig. 1.

Definition sketch of bottom variation with time.

Fig. 2.

Relation between Hammack’s paper and present study.

Fig. 3.

Definition sketch of bottom variaon with time in two regions.

Fig. 4.

Normalized free surface elevation for different up-rising time; (a) x =0; (b) t = t_rise.

Fig. 5.

Normalized free surface elevation for different time-delayed bottom variation.

Fig. 6.

Normalized free surface elevation for different gaps between two bottom variation regions.

References

Bryant, E. (2014). Tsunami. Springer.

Dutykh, D, Dias, F. (2007). Water waves generated by a moving bottom. Tsunami and Nonlinear Waves, In: Kundu A, ed. Geo Sciences. New York: Springer, p 63-94.

Dutykh, D., Dias, F., Kervella, Y. (2006). Linear theory of wave generation by a moving bottom, Comptes Rendus de l’Academie des Sciences Paris, Series I, 343, 499-504.

Duthkh, D., Mitsotakis, X., Gardeil, X., Dias, F. (2013). On the use of the finite fault solution for tsunami generation problem, Theoretical and Computational Fluid Dynamics, 27(1–2):177-199.

Filon, L.N.G. (1930). On a quadrature formula for trigonometric integrals, Proceedings of the Royal Society of London A, 49, 38-47.

George, D.L., LeVeque, R.J. (2006). Finite volume methods and adaptive refinement for global tsunami propagation and local inundation, Science of Tsunami Hazards, 24, 319-328.

Grilli, S.T., Ioualalen, M., Asavanant, J., Shi, F., Kirby, J.T., Watts, P. (2007). Source constraints and model simulation of the December 26, 2004, Indian Ocean Tsunami, Journal of Waterway, Port, Coastal, and Ocean Engineering, 133(6):414-428.

Gusiakov, V.K. (2009). Tsunami history: recorded, The Sea, 15, 23-53.

Hammack, J.L. (1973). A note on tsunamis: their generation and propagation in an ocean of uniform depth, Journal of Fluid Mechanics, 60(4):769-799.

Imai, K., Satake, K., Furumura, T. (2010). Amplification of tsunami heights by delayed rupture of great earthquakes along the Nankai trough, Earth, Planets and Space, 62(4):427-432.

Ishii, M., Shearer, P.M., Houston, H., Vidale, J.E. (2005). Extent, duration and speed of the 2004 Sumatra-Andaman earthquake imaged by the Hi-Net array, Nature, 435(7044):933.

Kajiura, K. (1963). The leading wave of a tsunami. Bulletin Earthquake Research Institute, Tokyo University, 41, 535-571.

Kervella, Y., Dutykh, D., Dias, F. (2007). Comparison between three-dimensional linear and nonlinear tsunami generation models, Theoretical and Computational Fluid Dynamics, 21, 245-269.

Kowalik, , et al. (2005). Numerical modeling of the global tsunami: Indonesian Tsunami of 26 December 2004, Science of Tsunami Hazards, 23(1):40-56.

Lee, J.-W., Park, E.H., Park, S.-C., Woo, S.-B. (2015). Development of the global tsunami prediction system using the finite fault model and the cyclic boundary condition, Journal of Korean Society of Coastal and Ocean Engineers, 27(6):391-405.

Ohmachi, T., Tsukiyama, H., Matsumoto, H. (2001). Simulation of tsunami induced by dynamic displacement of seabed due to seismic faulting, Bulletin of the Seismological Society of America, 91(6):1898-1909.

Pringle, W.J., Yoneyama, N., Mori, N. (2018). Multiscale coupled three-dimensional model analysis of the tsunami flow characteristics around the Kamaishi Bay offshore breakwater and comparisons to a shallow water model, Coastal Engineering Journal, 60(2):200-224.

Saito, T., Furumura, T. (2009). Three-dimensional tsunami generation simulation due to sea-bottom deformation and its interpretation based on the linear theory, Geophysical Journal International, 178(2):877-888.

Son, S., Lynett., P., Kim, D.-H. (2011). Nested and multi physics modeling of tsunami evolution from generation to inundation, Ocean Modelling, 38, 96-113.

Suppasri, A., Imamura, F., Koshimura, S. (2010). Effects of the rupture velocity of fault motion, ocean current and initial sea level on the transoceanic propagation of tsunami, Coastal Engineering Journal, 52(2):107-132.

Suzuki, W, Aoi, S, Sekiguchi, H, Kunugi, T. (2011). Rupture process of the 2011 Tohoku-Oki mega-thrust earthquake (M9. 0) inverted from strong-motion data, Geophysical Research Letters, 38(7.

Titov, V.V., Synolakis, C.E. (1998). Numerical modeling of tidal wave runup, Journal of Waterway, Port, Coastal, and Ocean Engineering, 124(4):157-171.

Wang, X. (2009). User manual for COMCOT version 1.7 (first draft). Cornel University, p 65.

Wilson, R.I., Admire, A.R., Borrero, J.C., Dengler, L.A., Legg, M.R., Lynett, P., McCrink, T.P., Miller, K.M., Ritchie, A., Sterling, K., Whitmore, P.M. (2013). Observations and impacts from the 2010 Chilean and 2011 Japanese tsunamis in California (USA), Pure and Applied Geophysics, 170(6–8):1127-1147.