伊豆衝突帯北縁部における震源及びメカニズム解分布の特徴
行竹洋平＊・武田哲也＊ 2・吉田明夫＊Detailed spatial distribution of hypocenters and focal mechanisms in the collision zone of the Izu-Bonin arc
byYohei YUKUTAKE , Tetsuya TAKEDA ＊ 2 and Akio YOSHIDA ＊
＊

1.

はじめに

138.5˚

139.0˚

139.5˚

山梨・神奈川県境付近の丹沢地域は、伊豆衝突帯北縁
36.0˚

部に位置し、フィリピン海プレートが複雑な形状で沈み込む常時地震活動が高い領域である。この領域におけるフィリピン海プレートの形状については、震源分布、3次元速度トモグラフィー法、地震波構造探査法の結果に基づき推定されてきた（例えば、Ishida, 1992; Matsubara

35.5˚

et al., 2005; Sato et al., 2005）．しかしながら、この領域におけるプレート構造が現在発生している地震活動とどのような関係にあるかについて十分な理解に至っていな
い．

35.0˚

この地域のテクトニクスを明らかにするために、震源およびメカニズム解を高精度に決定することが重要となる。本研究では、この地域に発生する地震の震源および50 km

メカニズム解を定常地震観測点の波形データを用い高精度に決定することを試み、詳細な分布の特徴について報図１使用した観測点分布。▼は温泉地学研究所、■は告する。この地域のメカニズム解については、行竹ほか

防災科研 Hi-net、▽は気象庁の定常観測点の位置を（2006）において推定されているが、温泉地学研究所の表す。灰色線は、活断層の地表トレースを示す。点線

ルーチン処理に利用されている観測点しか用いていなか

四角領域は、図２で示された領域に対応する。

った。この研究では、他研究機関も含むより広範囲の観測点データを使用し、さらに地震波の P 波初動極性だ
けでなく、振幅値も取り入れたメカニズム解決定を行い、までに発生した 4500 イベントの地震について解析を行

決定精度の向上を行った。なお、本論では丹沢地域と伊った。これらの地震のマグニチュード範囲は、0.0 から豆衝突帯北縁部は同じ領域を指していると見なし、以下5.0 である。

議論を進めていく。
2.1 Double Difference 法を用いた高精度相対震源決定2.

データおよび手法

筆 者 ら は、 上 記 の 地 震 の 震 源 位 置 を、Double

筆者らは、震源およびメカニズム解を決定するにあた

Difference 法（ 以 下、DD 法 ）（Waldhauser and Ellsworth,

り、温地研観測点に加え、防災科学技術研究所高感度地

2000）を用いて決定した。なお、DD 法の詳細について

震観測網（Hi-net）、気象庁地震観測点による、107 点のは、行竹・棚田（2008）に記載されているため、本論で定常地震観測点を使用した（図 1）。図 1 の点線で示さの説明は省略する。DD 法を行うための初期震源位置は、れた矩形領域内において、2001 年 1 月から 2008 年 5 月

hypomh 法（Hirata and Matsu'ura, 1987）を用いて決定した。

＊神奈川県温泉地学研究所＊ 2 防災科学技術研究所〒 250-0031〒 305-0006

神奈川県小田原市入生田 586
茨城県つくば市天王台 3 − 1

報告 , 神奈川県温泉地学研究所報告 , 第 42 巻 ,9-18, 2010



(a)

(b)
A
C
E

20

35.4˚

D
F

0

Y, km

B

35.2˚

−20

138.8˚0

139.4˚

139.2˚

139.0˚1020depth, km

0Depth, km 2030

図２ （a）震央分布。（b）南北深さ断面。プロットの色の濃淡は震源の深さを表す。点線矩形領域は、図３において示される、A-B、C-D および E-F 断面の位置を表す。

その際、平賀（1987）による丹沢地域の一次元速度構造を使用した。

2.2 メカニズム解決定メカニズム解の決定にあたっては精度を高めるため、

本 研 究 で は、2 つ の イ ベ ン ト の 走 時 差（Double

Ide et al.,（2003）の手法に従い P 波初動極性データに加Difference） を 相 互 相 関 処 理（ 波 形 相 関 デ ー タ ） な ら

えて、P 波ならびに SH 波の変位スペクトルにおける低びに検測時刻（カタログデータ）から求め、初期震源

周波側の振幅情報も用いた。P 波極性が 12 観測点以上

位置の再決定に用いた。手動検測値から求めた Double

で読み取られている地震について、メカニズム解の決定

difference データは、P 波走時に対して 55 万ペア S 波走を試みた。上記の手法で、決定されたメカニズム解の各時に対して 50 万ペアである。加えて、波形相互相関処

観測点での理論振幅値と観測振幅値との比から観測点振理によって得られたより高精度な Double difference デー

幅補正値を求めた。そして、決定された観測点振幅補正

タも震源決定に用いた。相互相関処理には、P 波および値を用いて、
再び上記の手順でメカニズム解を決定した。

S 波の手動で読み取られた到達時刻を含む 0.75 秒間のこれにより、地表地盤による振幅増幅あるいは減衰特性

速度波形を用い、3-20Hz の帯域のバンドパスフィルタの影響が補正され、より精度よくメカニズム解を決定すーを使用し、相関係数が 0.8 以上ある Double difference

ることが可能になる。その結果、822 イベントのメカニ

データのみを使用した。相互相関処理によって得られたズム解を決定することができた。ただし、地震の規模がDouble difference データは、P 波走時に対して 18 万ペア

大きくなると（断層サイズが大きくなると）、P 波初動S 波走時に対して 9 万ペアである。DD 法を適応するこに関係する断層形状と P 波初動後に記録される最大振とにより、Double difference データの RMS 残差は手動幅に関係する断層形状との違いが大きくなる可能性が考

検測値に対して 192ms から 130ms へ、相互相関データえられる。その場合、上記の P 波初動と P 波 S 波振幅に対して 87 ms から 6ms へそれぞれ減少した。

とを併せた解析では、解が不安定になる可能性が出てくる。そこで、マグニチュードが 3 を超える中規模地震については、F-net によるモーメントテンソル解カタログを使用した。

B0

結果

3.1 震源分布
DD 法により決定された震央分布および南北深さ断面を図 2 に示す。丹沢地域直下の深さ 10km から 30km のDepth, km

3.

(a) A

Region A

20

範囲に、震源分布のかたまり（震源クラスター）が明瞭に見え、この地域の地震活動が高いことを反映している。
この震源クラスターと丹沢地域におけるプレート構造との議論をするため、図 2 の矩形にて示された A-B、C-D

Region B
(b) C

D

0

方向は、この地域のフィリピン海プレートの走向（例えば、Sato et al., 2005）と概ね平行になるように設定した。
図 3 から、丹沢地域直下では西側と東側で震源分布の特徴が大きくことなることが分かった。東側の（図 3 にDepth, km

および E-F 領域における深さ断面を図 3 に示す。断面のRegion A

20

おいて Region A と示された）領域では、震源はほぼ水平な面上に集中して分布する。一方で、西側の（図 3 において Region B と示された）領域では、震源は Region

Region B
(c) E

F

0

Depth, km

A において発生した地震と比較して深い領域でかつ厚みをもって分布する。

3.2 メカニズム解空間分布の特徴
次に、地震活動とフィリピン海プレート構造との関係

20をさらに厳密に議論するため、メカニズム解の空間分

Region B

布の特徴を詳細に調べた。最初に Region A で発生した地震のメカニズム解の特徴を最も合理的に表す、メカ−20

ニズム解（Reference model A）を次の手法により推定し

0
km

20

た。仮定したリファレンス解とそれぞれの地震のメカニ
ズム解との類似性を評価するために、Kagan 角（Kagan,01993）を用いた。Kagan 角は、2 つのメカニズム解が完全に一致する場合 0 度となり、類似性がなくなるほど大

図３

10
20
depth, km

30

A-B、C-D および E-F 断面における震源深さ分布。

きくなる。Region A において発生した地震のメカニズ

プロットの色の濃淡は震源の深さを表す。点線で囲まム解との Kagan 角の 2 乗和が最小となるようなメカニれた領域は、それぞれ RegionA および B の位置を表す。

ズム解（Reference Model A）をグリットサーチにて推定した。
Reference Model A の断層パラメータおよびメカニズ

表１

Reference Model A および B の断層パラメータ

ム解を表 1 および図 4 にそれぞれ示す。グリットサーチの結果、東西走向に近く、北北東方向に低角に傾斜する、
低角逆断層型のメカニズム解が Reference Model A として推定された。Reference Model A の断層走向および傾

Reference Model A

strike287

dip rake
32 161

Reference Model B

195

30

40

斜は、この地域におけるフィリピン海プレート形状（例
えば、Sato et al., 2005）および 1923 年大正関東地震の断るフィリピン海プレートの相対的な運動方向である北西層モデル（Matsu'ura et al., 1980）と概ね一致する。さらに、−南東方向（Seno, 1993）と調和的である。

そのすべり方向（rake 角）はユーラシアプレートに対す

それぞれの地震のメカニズム解の Reference Model A

Reference Model A

(a) A

(d)

A

B

D

(e) C

D

F

(f ) EF

Depth, km

0Reference Model B

B

20

(b) C
Depth, km

0

20

(c) E

Depth, km

0

20

-20

0
x, km

20

-20

0
x, km
0

図４

20

40
80
kagan angle, degree

120

Reference Model A および B に対する、それぞれのメカニズム解の Kagan 角の大きさ。
（a）-（c）Reference

Model A に対する Kagan 角の大きさ。
（d）-（f）Reference Model B に対する Kagan 角の大きさ。プロット色の濃淡は、Kagan 角の大きさを表す。図の下端に下半球投影で表示した Reference Model A および B を示す。に 対 す る Kagan 角 を、A-B、C-D お よ び E-F 断 面 上 にRegion A のものとはことなる事を示唆している。

プ ロ ッ ト し 図 4（a）-（c） に 示 す。 さ ら に、Region A

Region B において発生した地震のメカニズム解の特お よ び B に お け る Kagan 角 の 頻 度 分 布 を 図 5（a） に徴を最も合理的に表すリファレンスとなるメカニズム

示 す。Region A に お い て 面 上 に 分 布 す る 地 震 の 多 く

解（Reference Model B）を、同様の手段で推定した。リが、Reference Model A に対して小さな Kagan 角を持つ。

ファレンスモデル B の断層パラメータおよびメカニズ

一 方 で Region B で 発 生 し た 地 震 の 多 く は、Reference

ム解を表 1 および図 4 にそれぞれ示す。Reference Model

Model A に対して大きな Kagan 角をもち、Kagan 角の頻

B は、南北走向に近い節面をもつ低角逆断層型であり、度分布のピークは 70 度付近にある（図 5（a））
。このこ

Reference Model A とは明瞭にことなる。A-B、C-D およRegion B で発生した地震のメカニズム解の特徴が、とは、

び E-F 断面で発生した地震のメカニズム解の Reference



(a)

(b)

Reference Model B

Normalized Number of events

Reference Model A

RegionA
Region B0
図52040 60 80 100
Kagan angle, degree

120

0

20

40 60 80 100
Kagan angle, degree

120

Region A および RegionB において発生した地震のメカニズム解に対する Kagan 角の頻度分布。Y 軸は最大頻度で正規化されている。
（a）Reference Model A（b）Reference Model B に対する Kagan 角の頻度分布。
Model B に対する Kagan 角を、図 4（d）-（f）に示す。

より詳細に見るため、Kagan 角が小さい 30 度以内のメRegion B に お い て Reference Model B に 対 し て 小 さ な

カニズム解のみを図 4 と同様な深さ断面にプロットした

Kagan 角を持つ地震が多く見られる。一方で、Region A（図 6）。図 6（a）（b）から、Region A において震源分布で発生した地震のメカニズム解は、Reference Model B

上端付近の灰色点線周辺に Reference Model A に対してに対して Kagan 角が大きくなる。

Kagan 角の小さなメカニズム解が集中して分布している

Region A および B において決定されたすべてのメカ

ことが分かる。以上の結果から、Region A で発生した地

ニズム解をそれぞれ、図 7 および 8 に下半球投影を用い

震の多くは、沈み込むフィリピン海プレート境界あるいて示す。図 7 および 8 では、それぞれ Reference Modelはそのごく近傍で発生していることが結論付けられる。

A および B に対する Kagan 角の小さいメカニズム解か

プレート境界面は、図 6（a）（b）において示した灰色

ら順にプロットした。Region A で発生した地震のメカ点線付近に存在しているのかもしれない。

ニズム解は、東西走向に近い節面を持つ低角逆断層型が一方で、Region B で発生した地震は、Region A のよ

多く見られる、一方、Region B において南北走向に近うな明瞭な面上にならんだ震源分布が見られず、厚み

い西側に傾斜する節面を持つ低角逆断層型の地震が多くを持って分布する。さらに、リファレンス解は Region

含まれることが分かる。また、どちらの領域においても、A のものとはことなり、南北走向に近くかつ西傾斜の低一部正断層型および横ずれ断層型の地震が含まれる。

角逆断層型であり Region A と明らかにメカニズム解の特徴がことなる（図 4 および図 7）。ただし、Reference
Model B に対して Kagan 角の小さなメカニズム解のみプ

4. 考察および結論
Region A において発生した地震は、フィリピン海プ

ロットした図 6 から、C-D 断面（図 6（e））の西側の灰

レートの走向に沿ってほぼ水平な面上に集中して分布し色点線付近においては Kagan 角の小さなメカニズム解ていることが明らかになった（図 3）。Reference Modelが集中して分布している場所が存在していることが分

A は、過去の研究から推定された丹沢地域東部のフィリかった。一方で、その南側の E-F 断面では（図 6（f））、

ピン海プレート境界面の形状と一致する節面をもち、かKagan 角の小さなメカニズム解が分散して分布していつそのすべり方向が関東地域のフィリピン海プレートのる。これらの結果に基づいて、Region B で発生した地運動方向と著和的である。また、永井・棚田（2010）は、

震活動の発生原因については、
以下の解釈が考えられる。

Region A 付近においてフィリピン海プレート境界を示

フィリピン海プレート形状が Region B において局所的

唆する、低速度と高速度域の境界が存在することを明らに変化し、南北走向で西側に傾斜する形状をもっているかにした。プレート境界とメカニズム解分布との関係を可能性が考えられる。図 6（e）の C-D 断面における灰Reference Model A (Kagan angle ≤ 30°) Reference Model B (Kagan angle ≤ 30°)
(a) A
(d)
B
B
0

A

Depth, km

Depth, km

0

20

20

(b) C

(e)

D

0

C

D

E

F

Depth, km

Depth, km

0

20

(c)

20

E

F

(f )

0
Depth, km

Depth, km

0

20

20

-20

0
x, km

20

-20

0
x, km
0

図6

40
80
kagan angle, degree

20

120

Reference Model A および B に対して、Kagan 角が 30 度以内のメカニズム解のみをプロット。詳細は図４と同じ。

色点線付近で発生した Kagan 角の小さなメカニズム解

域のテクトニクスについてより定量的に議論することがはプレート境界で発生した地震を反映しているのかもし必要となる。

れない。一方、震源が分散して分布していることから、
Region B で発生した地震の多くはフィリピン海プレー

5. まとめト内部で発生しると考えられる。図 6（f）E-F 断面でみ

本研究では、伊豆衝突帯北縁部にあたる丹沢地域およられた Kagan 角の小さなメカニズム解のかたまりはプ

びその周辺域で発生した地震の震源位置およびメカニズ

レート内部に形成されたフラクチャー構造を反映していム解について、定常地震観測点の地震波形データを用るのかもしれない。

今後は、応力インバージョン法にいて高精度に決定した。震源決定には DD 法を適用し、よりこの地域に作用する応力場を詳細に検討し、この地

Double difference データの推定には手動検測値および波20041226

20040602

20040913

20030511

20060817

20070329

20040216

20040921

20080420

20050614

M=1.8 K=2.9

M=2.1 K=4.9

M=2.5 K=5.0

M=1.9 K=7.2

M=2.6 K=9.8

M=1.9 K=9.8

M=1.9 K=9.9

M=1.9 K=10.2

M=1.9 K=10.6

M=1.7 K=10.8

20040227

20050625

20060104

20050515

20070722

20070209

20030927

20051116

20050821

20070615

M=2.6 K=10.8

M=1.6 K=11.1

M=2.4 K=11.4

M=1.9 K=11.9

M=1.6 K=12.0

M=1.5 K=12.1

M=1.8 K=12.6

M=1.6 K=13.3

M=1.9 K=13.4

M=1.9 K=13.8

20070529

20040603

20070417

20070921

20030509

20060328

20050410

20080401

20070405

20060209

M=1.6 K=13.8

M=2.4 K=13.9

M=1.9 K=14.3

M=1.9 K=14.4

M=2.7 K=15.2

M=1.6 K=15.4

M=1.9 K=15.5

M=1.7 K=15.5

M=1.6 K=15.8

M=2.1 K=15.8

20060311

20060219

20010202

20060428

20050718

20060225

20031202

20070210

20080330

20041003

M=1.8 K=15.9

M=1.8 K=16.3

M=4.7 K=16.5

M=2.2 K=16.6

M=2.3 K=17.1

M=2.1 K=17.2

M=2.1 K=17.8

M=2.0 K=18.1

M=2.3 K=18.5

M=1.7 K=18.7

20040923

20040602

20041029

20050513

20070420

20060414

20020328

20050307

20060722

20030713

M=1.9 K=18.8

M=2.3 K=18.9

M=1.7 K=18.9

M=2.0 K=19.1

M=2.0 K=19.1

M=1.6 K=19.9

M=4.2 K=20.1

M=2.1 K=20.6

M=1.5 K=22.7

M=2.0 K=22.9

20070129

20040603

20040602

20050804

20070117

20060325

20050820

20060710

20070512

20080209

M=1.7 K=23.2

M=1.6 K=23.4

M=2.2 K=24.9

M=1.9 K=25.2

M=2.0 K=26.8

M=1.8 K=27.1

M=1.6 K=28.6

M=2.2 K=30.5

M=1.8 K=30.5

M=1.5 K=30.6

20061229

20030813

20061201

20070712

20041127

20060122

20070720

20060705

20050528

20030331

M=2.5 K=31.4

M=1.8 K=31.7

M=2.3 K=31.9

M=3.8 K=33.4

M=1.9 K=33.9

M=1.8 K=34.6

M=1.6 K=35.0

M=1.9 K=35.2

M=1.5 K=36.4

M=2.2 K=37.1

20080420

20060904

20071003

20071120

20070712

20080218

20051001

20070625

20070702

20070106

M=1.8 K=37.2

M=1.9 K=37.5

M=1.8 K=38.3

M=1.9 K=38.5

M=1.9 K=38.6

M=1.8 K=38.7

M=1.7 K=39.8

M=1.6 K=41.0

M=1.6 K=42.0

M=1.7 K=42.3

20050727

20080106

20080129

20060224

20080216

20050728

20080324

20060101

20040913

20080518

M=1.9 K=43.5

M=2.0 K=43.6

M=1.5 K=46.8

M=1.9 K=46.8

M=1.6 K=47.7

M=1.8 K=47.9

M=1.8 K=49.6

M=1.4 K=51.1

M=1.9 K=52.1

M=1.6 K=52.9

20040725

20050315

20051123

20050629

20040811

20051123

20050412

20060817

20050401

20050114

M=1.6 K=57.6

M=2.0 K=60.0

M=2.3 K=65.3

M=1.9 K=69.0

M=2.0 K=69.4

M=1.7 K=73.1

M=1.5 K=73.1

M=1.8 K=74.9

M=1.7 K=75.2

M=1.8 K=75.8

20030711

20051101

20061021

20051116

20070518

20070206

M=1.9 K=76.4

M=2.0 K=76.7

M=2.0 K=77.2

M=1.8 K=80.8

M=1.9 K=85.8

M=1.7 K=86.1

図7

Region A において発生した地震のメカニズム解を下半球投影にて表示。各メカニズム解プロットの上に記載されている値は、上段が年月日、下段がマグニチュードと Kagan 角の大きさを示す。Reference ModelA に対して、Kagan
角が小さいメカニズム解から順にプロットされている。

20060722

20070429

20050925

20060911

20031229

20050729

20060623

20050926

20050418

20050613

M=2.5 K=6.2

M=1.5 K=6.9

M=2.0 K=7.0

M=1.8 K=8.2

M=1.9 K=8.4

M=1.9 K=9.3

M=1.6 K=9.6

M=2.4 K=10.2

M=2.0 K=11.1

M=1.8 K=11.7

20060127

20050907

20050203

20050412

20050428

20060809

20050716

20080312

20050331

20031226

M=1.8 K=12.5

M=2.4 K=12.7

M=1.8 K=12.9

M=2.0 K=13.4

M=2.4 K=13.4

M=1.8 K=13.5

M=2.4 K=13.9

M=2.5 K=14.5

M=2.1 K=14.6

M=2.1 K=15.8

20031210

20070525

20070410

20040807

20041031

20070903

20030616

20070427

20080513

20070803

M=2.2 K=16.0

M=1.6 K=16.5

M=2.3 K=17.3

M=1.5 K=17.5

M=3.6 K=18.6

M=1.7 K=18.6

M=2.2 K=18.9

M=1.8 K=19.4

M=1.8 K=19.9

M=1.9 K=20.0

20060527

20041213

20061227

20050315

20030624

20070514

20070304

20050413

20071103

20051010

M=1.7 K=20.2

M=2.0 K=20.3

M=1.9 K=20.9

M=2.0 K=21.3

M=1.8 K=21.5

M=1.6 K=21.8

M=3.7 K=22.0

M=2.1 K=22.3

M=1.7 K=22.3

M=1.7 K=22.7

20051229

20070626

20050412

20070909

20070721

20080310

20040811

20071220

20051201

20060407

M=2.2 K=23.0

M=2.0 K=23.0

M=2.0 K=23.3

M=1.9 K=23.7

M=1.7 K=23.9

M=1.8 K=24.0

M=2.0 K=24.1

M=1.7 K=24.2

M=2.2 K=24.4

M=2.1 K=24.8

20051123

20060225

20070208

20070304

20050401

20051018

20050920

20050612

20080229

20070923

M=2.3 K=24.8

M=1.8 K=24.8

M=1.6 K=25.0

M=2.0 K=26.5

M=1.7 K=26.8

M=1.7 K=26.9

M=1.7 K=27.0

M=1.6 K=27.3

M=2.0 K=28.5

M=2.2 K=28.9

20070312

20050317

20071109

20060209

20031215

20061203

20040509

20070925

20071108

20011113

M=1.8 K=29.0

M=1.7 K=29.2

M=2.4 K=29.3

M=1.9 K=29.3

M=1.8 K=29.4

M=2.5 K=29.8

M=2.4 K=30.2

M=1.8 K=30.4

M=1.6 K=30.9

M=4.0 K=31.6

20040507

20070803

20060919

20051101

20071224

20050218

20050731

20070128

20070116

20050529

M=2.4 K=32.1

M=1.6 K=32.3

M=1.9 K=32.8

M=2.0 K=32.8

M=1.9 K=33.8

M=2.1 K=34.2

M=4.5 K=35.2

M=1.7 K=35.4

M=2.5 K=35.5

M=1.9 K=36.0

20070917

20051120

20030516

20050602

20050704

20070914

20060224

20050414

20050416

20050412

M=1.6 K=37.0

M=1.7 K=37.3

M=2.0 K=37.7

M=2.4 K=38.4

M=1.9 K=39.8

M=2.2 K=41.4

M=1.6 K=41.5

M=2.5 K=41.5

M=1.8 K=42.2

M=4.0 K=42.7

20030806

20070808

20060614

20060127

20031203

20060415

20080401

20020325

20030824

20040322

M=2.8 K=43.2

M=2.1 K=44.3

M=1.5 K=44.6

M=1.8 K=50.1

M=2.2 K=51.4

M=2.0 K=51.7

M=2.1 K=51.8

M=3.6 K=52.2

M=2.4 K=54.7

M=1.8 K=58.0

20050331

20051120

20070518

20080517

2004110820070216

20050621

20051003

20070218

20070312

M=2.0 K=58.8

M=1.9 K=58.9

M=1.9 K=59.4

M=1.8 K=60.1

M=1.7 K=61.3

M=1.8 K=61.3

M=1.7 K=61.3

M=1.9 K=64.1

M=1.9 K=64.7

M=1.8 K=64.9

図8

Region B において発生した地震のメカニズム解を下半球投影にて表示。Reference Model B に対して、Kagan 角が小さいメカニズム解から順にプロットされている。



形相互相関処理を用いた。メカニズム解決定には、P 波

K（2005）Low-velocity oceanic crust at the top of

初動値のほかに、P 波振幅値および SH 波振幅値を用い

the Philippine Sea and Pacific plates beneath the

た。

Kanto region, central Japan, imaged by seismic

その結果、丹沢山地域の東側では震源はフィリピン海
プレートの走向に沿った面上に集中して分布し、メカニ

t o m o g r a p h y, J . G e o p h y s . R e s . , 11 0 , B 1 2 3 0 4 ,doi:10.1029/2005JB003673.

ズム解はフィリピン海プレート形状およびその運動方向

Matsu’ura, R. S. and I. Karakama（2005）, A point-processと調和的であることが分かった。一方、丹沢山地の西側

analysis of the Matsushiro Earthquake swarm sequence :では震源は東側より深い場所に、厚みを持って分布する

The effect of water on earthquake occurrence, Pure Appl.

ことが分かった。さらに、そこで発生した地震のメカニ

Geophys., 162, 1319-1345.

ズム解は、東側におけるメカニズム解の特徴とはことな

永井悟・棚田俊收（2010）伊豆衝突帯北東部における 3

り、南北走向に近い西傾斜の低角逆断層が多く発生して次元地震波速度構造と地震活動、日本地球惑星科学

いることが明らかになった。東側で発生する地震の震源

連合 2010 年大会予稿集、SSS024-04．

分布およびメカニズム解の特徴は、この領域で地震はフSato, H., Hirata, N., Koketsu, K., Okaya, D., Abe, S.,

ィリピン海プレートの境界面あるいはそのごく近傍で発Kobayashi, R., Matsubara, M., Iwasaki, T., Ito, T., Ikawa,

生していることを示している。西側で発生する地震の特T., Kawanaka, T., Kasahara, K. and Harder, S.（2005）

徴については、この領域では地震はフィリピン海プレー

Earthquake Source Fault Beneath Tokyo, Science, 309,

ト内部で発生していること、プレート形状が局所的に変462-464.
Seno, T.（1993）Relative plate motions in the Kanto-Tokai

化していることを示唆していると考えられる。

District, J. Geogr., 102, 374-380.
謝辞
気象庁観測点における地震波形データを使用した。F.

行竹洋平・棚田俊收（2008）箱根火山で発生する群発地

Waldhauser 博士には、hypoDD プログラムのソースコー

震の精密震源決定、神奈川県温泉地学研究所報告書、ドを提供していただいた。メカニズム解の決定には、東40、35-50.

京大学の井出哲博士から提供していただいたプログラム

行竹洋平・棚田俊收・本多亮・伊東博・原田昌武（2006）

を使用した．また F-net のモーメントテンソル解カタロ

神奈川県西部地域におけるメカニズム解の決定、神グを使用した．本報告中の図は、Generic Mapping Tools

奈川県温泉地学研究所報告書、38、69-76.

（GMT）（Wessel and Smith, 1995）により作成した。

Waldhauser, F. and W. L. Ellsworth（2000）A doubledifference earthquake location algorithm: Method and

参考文献

application to the Northern Hayward fault, Bull. Seismol.

平賀士郎（1987）箱根火山と箱根周辺海域の地震活動．

Soc. Am., 90, 1352-1368.
Wessel, P. and W. H. F. Smith（1995）New version of the

神奈川県温泉地学研究所報告，18，149-273．
Hirata, N. and M.Matsu’ura（1987）Maximum-likelihood
estimation of hypocenter with original time eliminatedusing nonlinear inversion technique, Phys. Earth Planet.
Interior, 47, 50-61.
Ide, S., G.C. Beroze, S. G. Prejean and W. L. Ellsworth（2003）, Apparent break in earthquake scaling due to path andsite effects on deep borehole recordings, J. Geophys.
Res., 108（B5）, 2271, doi:10.1029/2001JB001617.
Ishida, M.（1992）Geometry and relative motion of thePhilippine Sea plate and Pacific plant beneath the KantoTokai district, Japan, J. Geophys. Res., 97, 489-513.
Kagan, Y. Y.（1993）3-D rotation of double-couple earthquake
sources, Geophys. J. Int., 106, 709-716.
Matsubara, M., Hayashi, H., Obara, K. and Kasahara,generic mapping tools released, Eos Trans. AGU, 76,329.

