こ う えい フ ォ ー ラ ム第 22 号 / 2014.3

東京湾内における 2011 年東北地方太平洋沖地震津波の伝播特性および計算条件の違いによる計算精度の検討
STUDY ON INFLUENCES OF NUMERICAL CONDITIONS OF TSUNAMI SIMULATION AND
CHARACTERISTICS OF TSUNAMI PROPAGATION OF THE 2011 OFF THE PACIFIC COAST OFTOHOKU EARTHQUAKE TSUNAMI IN TOKYO BAY
小園裕司 * ・ 野島和也 * ・ 桜庭雅明 *Yuji KOZONO, Kazuya NOJIMA and Masaaki SAKURABA

This paper presents a study on influences of computational conditions on numerical simulationof Tsunami generated by the 2011 from Tohoku-Paciﬁc earthquake. In Tokyo bay, the characteristicreaction of water elevation to the Tohoku tsunami was studied. The relationship between the own
natural period of Tokyo bay and the wave period of the tsunami was discussed. Two conditions aredealt with for the comparison of computational accuracy. The ﬁrst one is the difference of shape of theinitial wave, i.e. difference of the fault models. The second one is the difference of area of computational
domain. Each of computational results were compared with observed water elevation in Tokyo bay.
Inﬂuences of the differences on the accuracy of the tsunami simulation in Tokyo bay were identiﬁed
from results.

Keywords : T sunami Simulation, the 2011 off the Paciﬁc Coast of Tohoku Earthquake, TokyoBay,the Own Natural Period,Computational Conditions, The Wave Period

いにまで検証を行なっていないのが現状である。

1． はじめに数値計算による津波の予測を行う際、 適切に計算モデル ・

東京湾周辺は、 過去に 1498 年明応東海地震津波、 1605

計算条件のチューニングを行う必要がある。 その方法の一つ

年 慶 長 東 海 地 震 津 波、 1703 年 元 禄 関 東 地 震 津 波、 1707に既往津波の再現計算が考えられる。 その際、 津波計算の外年宝永東海 ・ 南海地震津波、 1854 年安政東海地震津波、力となる断層モデルの設定は重要な要因となる。 その他の計

1960 年チリ地震津波や、 2011 年東北地方太平洋沖地震津算精度に影響を及ぼす要因として、 計算格子の解像度と計算波 （以下、 東北津波） などの多くの地震津波が来襲する津波範囲の設定も重要となる。 詳細領域でのある程度の計算精度常襲地域である。 2011 年 3 月 11 日に発生した東北津波では、を維持するためには適切な計算範囲と計算格子の解像度を有

東北地方沿岸をはじめ、 広範囲にわたり甚大な被害が生じた。

する事が重要である。

波源から 500km 以上離れた東京湾外の銚子港で 5.0m 前後、以上を踏まえ本研究では、 2011 年に発生した東北津波を東京湾口に位置する第二海堡で 1m 程度、 湾央に位置する対象に、 実績潮位データから東京湾に来襲した津波特性につ横須賀港や横浜港で 2m 弱の津波が観測された。 また同いて検討した。 また、 数値計算を用いて、 複数の断層モデル湾内では 1 波目よりも 2 波目が最大波となる事が確認されておによる再現性を検討するとともに、 湾内の津波の伝播特性につり、 最大波が第 1 波となっている湾外の験潮所とは異なる水位いて検討した。 最後に、 同モデルを用いて計算範囲の違いに変動を示した。 このような湾の特性や津波の伝播過程を明らか

よる計算精度の差異について検討した。

1）, 2）にすることは、 首都機能や多大なインフラを持つ沿岸都市域の防災減災対策において有用な知見となると考えられる。
同地震津波については、 東北地方の潮位データや被害状
3）, 4）

2． 研究方法（1） 観測値による東京湾における東北津波の特性 （検討 1）。

東北地震によって東京湾に来襲した津波について各験潮所

これらの断層モデルについては東北地方を中心に再現性のにおける 1 分間隔の実績潮位データを整理した。 これらのデー検証が行われているが、 東京湾周辺での検証事例は少ない。

タを用いてスペクトル解析を行い東京湾に来襲した津波のピー

東京湾周辺における再現性や有用性については、 佐々木ら 3）

ク周期について整理するとともに、 同湾の固有周期に関して検などによって行われているが、 各断層モデルのバージョンの違討した。 なお、 スペクトル解析の対象時間については、 2011

況に基づき再現した津波断層モデルが提案されている

年 3 月 11 日午後 2 時から翌日の午前 4 時までの 14 時間と*技術本部 中央研究所 総合技術開発部

した。
1

東京湾内における 2011 年東北地方太平洋沖地震津波の伝播特性および計算条件の違いによる計算精度の検討（2） 数値解析に用いる断層モデルの精度 （検討 2）
東北津波を対象に複数の断層モデルを用いて東京湾周辺の津波の再現計算を行い、 東京湾内の伝播状況を把握するとともに、 各断層モデルの再現性を検証した。
本検討で採用する断層モデルは国内で最も用いられている
藤井 ・ 佐竹モデル 3） および東北大モデル 4） とした。 藤井 ・ 佐

3． 数値解析法（1） 津波の数値解法
津波伝播の数値計算には、 基礎方程式に非線形長波方程式を採用し leap-frog 法、 スタッガード格子に基づいた平面 2次元の浅水流モデルを用いた 6）。 海底摩擦にはマニング式を用いている。

竹モデルについては海溝軸に 40 以上の小断層を配置し、 東地震発生後の津波の初期波形については、 断層モデルか北地方沿岸域における津波実績波形のインバージョン解析によら Okada （1992） の方法 7） により算出される海底地殻変動り各断層のすべり量を推定している。 また東北大モデルについの鉛直成分を海上面に与える方法を用いた。ては、 海溝軸に 10 枚の小断層を配置し、 すべり量以外のパラメータについては藤井 ・ 佐竹モデルに基づき設定している。 す（2） 計算条件

べり量については東北地方の浸水域に再現できるように調整し

1） 数値解析に用いる断層モデルの精度について （検討 1）ている。 両断層モデルは数回のバージョンアップを図っており、断層モデルについては東北津波を対象として、 藤井 ・ 佐竹

本検討においては各断層モデルの複数のバージョンを対象としモデル Ver. 4.0、 Ver. 4.2、 Ver. 6.2、 Ver. 7.0、 Ver. 8.0、 東て比較検討した。 各種の断層モデルで、 どの程度計算結果に北大モデルの Ver. 1.0、 Ver. 1.1 および Ver. 1.2 の計 8 つの違いがあるのかを明らかにすることを目的としている。

モデルを採用した。 なお、 藤井 ・ 佐竹モデルの Ver. 6.2、 Ver.
7.0、 Ver. 8.0 については、 従来の断層モデルとは異なり、 プ（3） 数値解析による計算範囲の影響 （検討 3）

レートの破壊伝播速度とすべりの時間変化を考慮することが提東京湾周辺の詳細領域を変化させた場合、 東京湾内の計算案されている。 本研究でもこれに倣い、破壊伝播速度 2.0km/s、結果に及ぼす影響を比較した。 東京湾外は伊豆諸島や小笠原各小断層におけるすべりの時間変化を設定した。 その他の断諸島など東京都島嶼部が点在している。 また東から来襲する津層モデルについては、 破壊伝播速度と各小断層におけるすべ波に対して壁となるように伊豆半島が位置しており、 反射 ・ 屈折りの時間変化は従来通り考慮しないものとした。 計算範囲につなど非常に複雑な波の挙動を示すことが予測される。 本研究でいては、 図－ 1 中の黒枠の領域 （領域 1） および図－ 2 中はこれらの影響を加味しつつ 4 通りの計算範囲を設定して、 湾の赤い枠の領域 （領域 2） とオレンジ色の領域 （領域 3） に内の計算結果を比較検討した。

示す通り東北地方から東海地方を広域のメッシュとして、 東京湾内を詳細の検討範囲として、 1,350m、 450m、 150m 格子を採用し、 各メッシュ間は領域接続を行った。 X 方向、 Y 方

領域1領域3ケース1

領域3ケース2

領域3ケース4<RNRKDPD<RNRVXND

領域2領域3ケース3

領域3領域2

図－ 1 検討 2 における計算領域
2

図－ 2 検討 3 における計算領域と験潮所位置

こ う えい フ ォ ー ラ ム第 22 号 / 2014.3

向の分割数については表－ 1 に示す通りである。 水深値、 標
高値については広域に ETOPO1 8） を、 東京湾周辺などの詳4． 検討結果

細域に数値地図 50m メッシュ （標高） および海底地形デー（1） 東北津波の東京湾内の津波来襲状況 （検討 1）

タを用いた。 なお、 本検討においては最小メッシュサイズが図－ 3 に東京湾内に位置する横須賀港、 横浜港における150m であるため、 横断方向に十分なメッシュ数を確保できな潮位変動のスペクトル解析を行った結果を示す。 これらの結果

いことを踏まえ、 河川について考慮は行なっていない。 また、より、 各験潮所でのピーク周波数は 1.5 × 10 － 2 （× 60Hz）

横須賀港の水深については海底地形デジタルデータが沖合の付近であり、 東日本大震災時に東京湾に来襲した津波周期はデータまでしか存在していなかったため、 港湾内の水深データ67 分程度が最も卓越していた事を確認した。 既往の研究 9） よについては沖合の水深値を延長する形で水深データを補間作

り東京湾には 60 分～ 90 分程度の固有振動周期が存在する成した。 海底摩擦についてはマニングの粗度係数 0.025 とし、事が報告されており、 今回の東北津波においては、 湾内で湾潮位条件については津波来襲時の東京湾における天文潮位水振動を引き起こしていた可能性があったと考えられる。

T.P.+0.125m とした。 境界条件については、 陸域と海域の境界は完全反射とし、 領域 1 の沖側境界は自由透過とした。 再現計算時間については地震発生から 12 時間とした。

（2） 数値解析に用いる断層モデルの精度 （検討 2）
図－ 4 に藤井 ・ 佐竹モデル、 東北大モデルの鉛直地殻変

2） 数値解析による計算範囲の影響について動量 （津波波源） の一覧を、 図－ 5 に本解析より得られた藤断層条件については、 検討 2 の結果を踏まえ藤井 ・ 佐竹

井 ・ 佐竹 Ver. 4.2 における東京湾周辺の津波伝播状況を示

モ デ ル Ver. 4.2 と し た。 計 算 範 囲 に つ い て は、 領 域 1、 2す。 これらの結果より、 東北沖で発生した津波は地震発生後と も に 検 討 2 と 同 様 の 範 囲 と し て、 領 域 3 の み を 変 化 さ せ、35 分程度で銚子に到達し、 60 分程度で東京湾入り口に到達

図－ 2 に示す通り以下に示す範囲について検討した。

し湾の奥に向かって湾内を伝播する。 地震発生後 120 分程●

ケース 1: 関東全域を包括●

ケース 2: 東京湾内と周辺の島嶼部を包括●

ケース 3: 東京湾と浦賀水道を包括 （検討 2 と同範囲）央で合成する。 また、 東京湾周辺には島嶼部がありそれらが●ケース 4: 東京湾のみ

湾外で複雑な反射 ・ 屈折を引き起こす事を確認した。 図－ 6度経過すると第 2 波が東京湾入り口に到達する。 その後東京湾内で、 湾奥で反射した第 1 波と湾奥に進行する第 2 波が湾

なお、 検討 2、 3 における計算結果は図－ 2 中に示していに各験潮所における津波波形についての実績値と計算値の比る横須賀港、 横浜港、 京浜港および第二海堡の実測津波波

較を、 表－ 2 に各験潮所における最大津波水位について実形

1） ,2）と比較した。
1.0e+2

表－ 1 その他の計算条件
計算格子間隔 領域 1 ： 1,350m （1300 × 1500）（分割数）

領域 2 ： 450m （800 × 788）
領域 3 ： 150m （1337 × 1481）

時間増分間隔 すべての領域でΔt = 0.1 秒
基礎方程式
計算手法

非線形長波方程式

1.0e-21.0e-31.0e-41.0e-51.0e-70.01

Frequency *60 f(Hz)1.0e-80.00

0.10

Yokohama

1.0e+01.0e-11.0e-21.0e-31.0e-41.0e-51.0e-6
1.0e-71.0e-80.00

1.0e-11.0e-6

1.0e+1

Power Spectrum(m2sec)

対象地震津波 2011 年東北地震津波断層モデル
藤井 ・ 佐竹モデルおよび東北大モデルの各バージョン
津波初期波源 Okada（1992） の 方 法 に よ り 断 層 モ デ ルから地殻変動量・初期津波波形を算出
計算再現時間 地震発生から 12 時間

1.0e-41.0e+2

陸側境界条件 領域 1 ～ 3 は完全反射境界

潮位条件1.0e-31.0e-70.00

沖合境界条件 領域 1 で自由透過
領域間
水位 ・ 流量を接続

水深データ

1.0e-21.0e-6

時間 ： リープフロッグ法

Manning の粗度係数 0.025
ETOPO1、海底地形データ、
数値地図50m メッシュ標高
T.P.+ 0.125m

1.0e-11.0e-5空間 ： スタッガード格子

海底摩擦

1.0e+0Power Spectrum(m2sec)

内容

1.0e+11.0e+0Power Spectrum(m2sec)

項目

1.0e+2

Yokosuka

1.0e+10.01

Frequency *60 f(Hz)

0.10



図－ 3 スペクトル解析結果
3

東京湾内における 2011 年東北地方太平洋沖地震津波の伝播特性および計算条件の違いによる計算精度の検討
績値と計算値の誤差率を示す。 これらの結果から、 いずれのVer. 6.2、 Ver. 7.0、 Ver. 8.0 は Ver. 4.0、 Ver. 4.2 と比較し断層モデルにおいても実績波形を概ね再現可能であることをて過大に評価していること、 横須賀港や京浜港などの第一波確認した。

については藤井 ・ 佐竹モデルの Ver. 4.2 が最も再現性が高い事を確認した。

最大波の津波到達時間については、 破壊伝播速度と各小断層におけるすべりの時間変化を考慮した藤井 ・ 佐竹モデル

いずれの断層モデルにおいても、 横須賀港の第 1、 第 2 波の Ver. 6.2、 Ver. 7.0、 Ver. 8.0 がその他のモデルと比較し後の引き波が過大評価されている点については、 当該地点のて若干ではあるが実績値の再現性が高い。 最大波について水深値を沖合の水深値を延長する形で補間して作成していたは、 藤井 ・ 佐竹モデルはいずれも平均 13 ～ 14% 前後、 東め、 水深値が過大評価の可能性がある。 そのため、 流体の底北大モデルは平均 25% 前後の誤差率となることを確認した。

面摩擦の影響が小さくなることで水面形の変化が大きくなった

これは津波波形からも分かる通り、 東北大モデルは藤井 ・ 佐と考えられる。 この点については今後、 港湾計画図等の詳細竹モデルと比較して若干過少評価であることからも確認される。

なデータを入手して水深値と計算結果の改善を図っていく予定

なお、 藤井 ・ 佐竹モデル間では最大ピーク値こそ大きな差である。

は見られないが、 短周期な波が卓越する 20 時以降になると、藤井･佐竹v4.0

藤井･佐竹v4.2

藤井･佐竹v6.2

藤井･佐竹v7.0

0.0
0.0
0.0

-2.0

0.0
4.0 10.0-2.0 0.0

0.0

4.0
6.0
8.0

8.0 10.0
6.0-2.0 8.0
6.0
4.0
2.0

0.0

2.0

0.0

-2.0

4.0

0.0

0.0

2.0

0.0

東北大v1.1

東北大v1.0

0.0

0.0

2.0 0.0

0.0

0.0

藤井･佐竹v8.0

6.0

2.0

4.0

東北大v1.2

-2.0

-2.0

2.0
4.0

0.0

-2.0

2.0
0.0

0.0

6.0

-2.0
-4.0

8.0 8.0
6.0
10.0
4.0

0.0

-2.0
-4.0

0.0

8.0
6.0
4.0

-2.0
-4.00.0

2.0

2.0

0.0

0.0
8.0
6.0
4.00.0

8.0
6.0
4.0

2.0
0.0

0.0

0.0

図－ 4 波源モデル一覧 

Choshi

Yokohama

Yokohama

KeihinYokosuka

Choshi

Yokosuka

KeihinDainikaihou

Yokosuka

第１，２波の合成

Choshi

Yokohama

Yokohama

KeihinDainikaihou

Choshi

第１波先端
Dainikaihou

KeihinYokosuka

Dainikaihou

第１波先端

第２波先端

第１波先端

35min (15:21)

60min (15:46)

120min (16:46)

図－ 5 津波の伝播状況
4180min (17:46)

こ う えい フ ォ ー ラ ム第 22 号 / 2014.3
2.0<RNRVXND計算(藤井･佐竹v8.0)

計算(藤井･佐竹v7.0)

計算(藤井･佐竹v6.2)

計算(藤井･佐竹v4.2)

計算(藤井･佐竹v4.0)

計算(東北大v1.2)

計算(東北大v1.1)

計算(東北大v1.0)

実績(横須賀港)

津波高 P

1.0

0.0

-1.0

-2.0
14:30

15:30

16:30

17:30

18:30

19:30

20:30

21:30

22:30

23:30

0:30

1:30

2:30

時刻 （時：分）2.0<RNRKDPD

津波高 P

1.0

計算(藤井･佐竹v8.0)

計算(藤井･佐竹v7.0)

計算(藤井･佐竹v6.2)

計算(藤井･佐竹v4.2)

計算(藤井･佐竹v4.0)

計算(東北大v1.2)

計算(東北大v1.1)

計算(東北大v1.0)

実績(横浜港)

0.0

-1.0

-2.0
14:30

15:30

16:30

17:30

18:3019:30

20:30

21:30

22:30

23:30

0:30

1:30

2:30

時刻 （時：分）2.0.HLKLQ

津波高 P

1.0

計算(藤井･佐竹v8.0)

計算(藤井･佐竹v7.0)

計算(藤井･佐竹v6.2)

計算(藤井･佐竹v4.2)

計算(藤井･佐竹v4.0)

計算(東北大v1.2)

計算(東北大v1.1)

計算(東北大v1.0)

実績（京浜港）

0.0

-1.0

-2.0
14:30

15:30

16:30

17:30

18:30

19:30

20:30

21:30

22:30

23:30

0:30

1:30

2:30

時刻 （時：分）

2.0'DLQLNDLKR

津波高 P

1.0

計算(藤井･佐竹v8.0)

計算(藤井･佐竹v7.0)

計算(藤井･佐竹v6.2)

計算(藤井･佐竹v4.2)

計算(藤井･佐竹v4.0)

計算(東北大v1.2)

計算(東北大v1.1)

計算(東北大v1.0)

実績（第二海堡）

0.0

-1.0

-2.0
14:30

15:30

16:30

17:30

18:30

19:30

20:30

21:30

22:30

23:30

0:30

1:30

2:30

時刻 （時：分）図－ 6 断層モデルの違いによる水位経時変化の比較（最上段 ： 横須賀港、 中段上 ： 横浜港、 中段下 ： 京浜港、 最下段 ： 第二海堡）（3） 数値解析の計算範囲の影響 （検討 3）

ケース 2 については、 ピーク値や波形などケース 1 とほぼ

図 － 7 に ケ ー ス 1 ～ 4 ま で の 津 波 波 形 の 計 算 値 を 示 す。

同様の計算結果を示した。 ケース 3 についてもケース 1、 ケー

図中には実測値もあわせて示す。 本検討では、 最も計算領域

ス 2 とほぼ一致する結果となったが横須賀港での第 4 波以降が大きいケース 1 の計算結果をより真値に近い値として、 そのに若干ではあるがケース 1、2 とは異なる挙動を確認した。 ケー

他の計算ケースの結果を比較検討した。

ス 4 については、 第 1 波～第 3 波まではケース 1 と同様の津5

東京湾内における 2011 年東北地方太平洋沖地震津波の伝播特性および計算条件の違いによる計算精度の検討
表 –2 津波最大値の誤差

断層モデル

津波最大値の誤差率 （%）
横浜港

藤井･佐竹Ver.4.0
藤井･佐竹Ver.4.2
藤井･佐竹Ver.6.2
藤井･佐竹Ver.7.0
藤井･佐竹Ver.8.0
東北大Ver.1.0
東北大Ver.1.1
東北大Ver.1.2

京浜港㻝㻚㻣㻡㻚㻣
㻢㻚㻤
㻞㻚㻥
㻠㻚㻜
㻞㻜㻚㻡
㻞㻝㻚㻠
㻝㻢㻚㻝

横須賀港
㻟㻠㻚㻞
㻞㻤㻚㻝
㻞㻣㻚㻜
㻟㻡㻚㻞
㻟㻢㻚㻠
㻟㻥㻚㻠
㻟㻤㻚㻡
㻠㻜㻚㻟

㻝㻜㻚㻠
㻝㻠㻚㻡㻝㻡㻚㻣㻡㻚㻝㻟㻚㻥㻝㻟㻚㻥
㻝㻠㻚㻞
㻤㻚㻥

第二海堡

2.0

<RNRVXND㻝㻟㻚㻟㻝㻟㻚㻡㻝㻟㻚㻟㻝㻟㻚㻥
㻝㻠㻚㻠
㻞㻡㻚㻡
㻞㻢㻚㻝㻞㻞㻚㻣計算(CASE1)

計算(CASE2)

計算(CASE3)

計算(CASE4)

実績(横須賀港)

1.0

津波高 P

平均値

㻣㻚㻜㻡㻚㻣㻟㻚㻢㻝㻞㻚㻠㻝㻟㻚㻟㻞㻤㻚㻝㻟㻜㻚㻝㻞㻡㻚㻢

0.0

-1.0

-2.0
14:30

15:30

16:30

17:30

18:30

19:30

20:30

21:30

22:30

23:30

0:30

1:30

2:30

時刻 （時：分）2.0<RNRKDPD

計算(CASE2)

計算(CASE3)

計算(CASE4)

実績(横浜港)

1.0

津波高 P

計算(CASE1)

0.0

-1.0

-2.0
14:30

15:30

16:30

17:30

18:30

19:30

20:30

21:30

22:30

23:30

0:30

1:30

2:30

時刻 （時：分）図－ 7 計算範囲の違いによる水位経時変化の比較 （上段 ： 横須賀港、 下段 ： 横浜港）

波挙動を示すが、 それ以降の短周期の波が他のケースと比較して異なる挙動を示した。 これらの計算結果は東京湾外で知見が得られた。
●

スペクトル解析によって東京湾内における東北津波固の影響が東京湾内の計算結果に寄与したのではないのかと考有周期が既往の報告の固有振動周期と同等の値を示しえられる。 特に湾外に比較的近い横須賀の第 4 波以降では、ていることが分かった。 これにより湾水振動を起こしてい東京湾外の反射の影響を受けるためケース 4 よりケース 1 のた可能性を示した。

計算結果がより実測値を再現していることを確認した。

●

東北津波は地震発生から 60 分程度で東京湾に到達し、およそ 60 分周期で来襲を繰り返す。 また複数の津波が東京湾に浸入し湾内で波が合成される。 これらのプ5． おわりにロセスより東京湾内では 2 波以降が最大波となることを明らかにした。

本研究では、 東北津波を対象に、 東京湾内における津波の伝播状況と特性を把握するとともに、 複数の断層モデルによる●測波形を概ね再現可能である。

再現性を検証した。 また、 東京湾周辺の計算領域を変化させてその計算精度の差異について検討した。 その結果、 以下の6

いずれの断層モデルにおいても、 東京湾においては実●

破壊伝播速度と断層におけるすべりの時間変化を考慮こ う えい フ ォ ー ラ ム第 22 号 / 2014.3

した断層モデルは、 考慮していないモデルと比較して4）津波の津波到達時間の再現性が高い。
●

東北地方太平洋沖地震を対象とした津波シミュレーションの実施
東北大学モデル Ver.1.0

東京湾内の挙動を適切に捉えるためには、 湾外の挙動も捉えるための計算範囲の設定を行う必要がある。

今村文彦、 越村俊一、 大家隆行、 馬渕幸雄、 村嶋陽一 （2011） :5）

佐々木敦、 伊藤一教、 鈴木祟之、 Retno Utami Agung Wiyono、
織田幸伸、 高山百合子、 羽角華奈子、 古田敦史、 高木泰士 :2011

今後は、 その他の断層モデルについても再現性について検年東北地方太平洋沖地震津波の東京湾内での伝播特性と被害状
況、 海岸工学論文集、 Vol.68、 No.2、 pp.261-265.

討するとともに、 また、 適切な計算範囲の設定方法について検討していく予定である。

6）

後藤智明、 小川由信 :Leap-Frog 法を用いた津波の数値計算法、東北大学工学部土木工学科、 52p、1982

参考文献
1）

国土交通省 : 海上保安庁、 http://www.kaiho.mlit.go.jp/

2）

国土交通省港湾局他 : ナウファス、

7）

AND TENSILE FAULTS IN A HALF-SPACE, Bulletin of the
Seismological Society of America, Vol.82, No.2, pp.1018-1040, 19928）

http://www.mlit.go.jp/kowan/nowphas/index.html
3）

Earthquake, Earth Planets Space , Vol.63, pp.815-820

NATIONAL GEOPHSICAL DATA CENTER :ETOPO1
Global Relief Model,Fujii, Y, K Satake, S.Sakai, M.Shinohara and T.Kanazawa. :Tsunami source of the 2011 off the Pacific coast of Tohoku

Okada.Y.:INTERNAL DEFORMATION DUE TO SHEAR

http://www.ngdc.noaa.gov/mgg/global/9）

日野幹雄、 日野和子 : 東京湾の進入長波に対する振動応答特性、
第 11 回海岸工学講演会講演集、 pp.98-107、 1964

7