６－４

東海・紀伊半島・四国における短期的スロースリップイベント（2015 年 5 月～2015 年 10 月）
Short-term slow slip events in the Tokai area, the Kii Peninsula and the Shikoku
District, Japan (from May 2015 to October 2015)
産業技術総合研究所
Geological Survey of Japan, AIST
防災科学技術研究所
National Research Institute for Earth Science and Disaster Prevention

歪・傾斜・地下水の変化から短期的 SSE の断層モデルを推定したイベントについて，その解析結果を報告する．
2015 年 5 月 7 日頃から 10 日頃にかけて，三重県中部〜伊勢湾にかけて深部低周波微動活動が観測された（第 1 図）．第 2 図は同時期の微動活動周辺の産総研・気象庁・防災科研の観測点における多成分歪・水位・体積歪・傾斜の観測結果である．これらの結果は BAYTAP-G により気圧応答成分，潮汐成分およびホワイトノイズ成分を取り除き，さらに 2015 年 5 月 1 日から 5 月 6 日までのデータを用いて 1 次トレンドを除去したものである．
第 2 図の網掛部分に示した変化を用いて短期的 SSE の断層面を推定した結果を第 3 図に示す．
モーメントマグニチュードは Mw 5.7 と推定された．
今回の活動域付近における最近の短期的 SSE には，和歌山県・奈良県県境付近〜愛知県東部の領域で 2015 年 1 月 1 日頃から 21 日頃にかけて発生したイベント（発生順に Mw 5.6, 5.7, 5.8, 5.6,
5.2, 5.8; 第 3 図の灰色矩形 1A-1F），愛知県東部〜愛知県・長野県県境付近で 2015 年 4 月 4 日頃から 28 日頃にかけて発生したイベント（発生順に Mw 5.8, 5.7, 5.7, 5.3; 第 3 図の灰色矩形 2 および
3A-3C）がある．
2015 年 5 月 14 日頃から 16 日頃，および 5 月 26 日頃から 6 月 4 日頃にかけて，豊後水道から四国西部の領域で深部低周波微動活動が観測された（第 4 図）．第 5 図は同時期の微動活動周辺の産総研・防災科研の観測点における多成分歪・傾斜・地下水の観測結果である．これらの結果はBAYTAP-G により気圧応答成分，潮汐成分およびホワイトノイズ成分を取り除き，さらに 2015 年
5 月 1 日から 5 月 12 日までのデータを用いて 1 次トレンドを除去したものである．
第 5 図に示した変化の状況と微動活動の推移から，期間を[A] 2015 年 5 月 14 日〜5 月 16 日，[B1] 5 月 26 日 〜5 月 29 日午前，[B2] 5 月 29 日午後〜5 月 31 日午前および[B3] 5 月 31 日午後〜6
月 2 日に分け，網掛部分に示した変化を用いて短期的 SSE の断層面を推定した結果を第 6-9 図に示す．モーメントマグニチュードは順に Mw 5.5, 5.6, 5.9, 5.7 と推定された．なお前述の通り微動活
動は 6 月 4 日頃まで見られたが，観測値に対応する明瞭な変化が見られなかったことから，解析は6 月 2 日までとしている．
今回の活動域付近における最近の短期的 SSE には，2014 年 12 月 26 日頃から 2015 年 1 月 12 日頃にかけて発生したイベント（順に Mw 6.1, 5.8, 5.9; 第 6-9 図の灰色矩形 1A-1C），2015 年 3 月 24

－ 255 －

日頃から 27 日頃にかけて発生したイベント（順に Mw 5.8, 5.3; 第 6-9 図の灰色矩形 2A, 2B）がある．
2015 年 6 月 4 日頃から 7 日頃にかけて，奈良県南部で深部低周波微動活動が観測された（第 10図）．第 11 図は同時期の微動活動周辺の産総研の観測点における多成分歪の観測結果である．これらの結果は BAYTAP-G により気圧応答成分，潮汐成分およびホワイトノイズ成分を取り除き，さらに 2015 年 5 月 20 日から 6 月 3 日までのデータを用いて 1 次トレンドを除去したものである．
第 11 図の網掛部分に示した変化を用いて短期的 SSE の断層面を推定した結果を第 12 図に示す．
モーメントマグニチュードは Mw 5.7 と推定された．
今回の活動域付近における最近の短期的 SSE には，同じ領域で 2015 年 1 月 7 日頃から 10 日頃にかけて発生したイベント（Mw 5.8; 第 12 図の灰色矩形 1），西側の領域で 2015 年 3 月 14 日頃から 16 日頃にかけて発生したイベント（Mw 5.3; 第 12 図の灰色矩形 2）がある．
2015 年 7 月 18 日頃から 25 日頃にかけて，三重県中部で深部低周波微動活動が観測された（第13 図）．第 14 図は同時期の微動活動周辺の産総研・気象庁・防災科研の観測点における多成分
歪・傾斜の観測結果である．これらの結果は BAYTAP-G により気圧応答成分，潮汐成分およびホ
ワイトノイズ成分を取り除き，さらに 2015 年 7 月 1 日から 7 月 15 日までのデータを用いて 1 次ト
レンドを除去したものである．
第 14 図の網掛部分に示した変化を用いて短期的 SSE の断層面を推定した結果を第 15, 16 図に示す．このイベントは微動活動の推移に対応して 2 つに分けられ，推定されたモーメントマグニチュ
ードは順に Mw 5.5, 5.6 であった．
今回の活動域付近における最近の短期的 SSE には，2014 年 7 月 3 日頃から 14 日頃にかけて発生したイベント（順に Mw 5.4, 5.6, 5.5, 5.8, 5.7, 5.7; 第 15, 16 図の灰色矩形 1A-1F），2015 年 1 月 1
日頃から 5 日頃にかけて発生したイベント（順に Mw 5.6, 5.7; 第 15, 16 図の灰色矩形 2A, 2B）がある．
2015 年 9 月 1 日頃から 5 日頃にかけて，豊後水道から四国西部の領域で深部低周波微動活動が観測された（第 17 図）．第 18 図は同時期の微動活動周辺の産総研・防災科研の観測点における多成分歪・傾斜の観測結果である．これらの結果は BAYTAP-G により気圧応答成分，潮汐成分およびホワイトノイズ成分を取り除き，さらに 2015 年 8 月 16 日から 31 日までのデータを用いて 1 次
トレンドを除去したものである．
第 18 図に示した変化の状況と微動活動の推移から，期間を[A] 2015 年 9 月 1 日〜3 日，[B] 9 月
4 日〜5 日に分け，網掛部分に示した変化を用いて短期的 SSE の断層面を推定した結果を第 19, 20
図に示す．モーメントマグニチュードは順に Mw 6.0, 5.5 と推定された．
今回の活動域付近における最近の短期的 SSE には，2015 年 3 月 24 日頃から 27 日頃にかけて発生したイベント（順に Mw 5.8, 5.3; 第 19, 20 図の灰色矩形 1A, 1B），2015 年 5 月 14 日頃から 6 月
2 日頃にかけて発生したイベント（順に Mw 5.5, 5.6, 5.9, 5.7; 第 19, 20 図の灰色矩形 2A-2D）がある．
2015 年 10 月 25 日頃から 27 日頃にかけて，三重県中部〜奈良県東部で深部低周波微動活動が観－ 256 －

測された（第 21 図）．第 22 図は同時期の微動活動周辺の産総研・気象庁・防災科研の観測点における多成分歪・傾斜・体積歪の観測結果である．これらの結果は BAYTAP-G により気圧応答成分，
潮汐成分およびホワイトノイズ成分を取り除き，さらに 2015 年 10 月 10 日から 24 日までのデータを用いて 1 次トレンドを除去したものである．
第 22 図の網掛部分に示した変化を用いて短期的 SSE の断層面を推定した結果を第 23 図に示す．
モーメントマグニチュードは Mw 5.8 と推定された．
今回の活動域付近における最近の短期的 SSE には，2015 年 5 月 7 日頃から 10 日頃にかけて発生したイベント（Mw 5.7; 第 23 図の灰色矩形 1），2015 年 7 月 18 日頃から 25 日頃にかけて発生したイベント（順に Mw 5.5, 5.6; 第 23 図の灰色矩形 2A, 2B）がある．
解析方法
短期的 SSE の断層面推定には，それぞれの観測点の水平歪 4 成分・体積歪・傾斜 2 成分・地下水圧の記録を用いる．地下水圧は，O1 および M2 分潮の振幅を BAYTAP-G [Tamura et al., 1991]により計算し，GOTIC2 [Matsumoto et al., 2001]により推定した地球個体潮汐および海洋荷重潮汐（O1および M2 分潮）との振幅比を用いて，体積歪に変換する．歪・傾斜・地下水ともに，観測波形から BAYTAP-G により，気圧応答成分，潮汐成分およびホワイトノイズ成分を取り除く．また，イベント直前の期間を用いて 1 次トレンドも取り除く．微動活動も参考にして，数時間～半日単位で活動開始・終了時期を判断し，その期間の変化量を短期的 SSE による変化量とする．その際，歪については Matsumoto et al. [2010]の手法で理論潮汐歪を用いてキャリブレーションを行っている．
断層面の推定は，計算時間の短縮と，推定された結果の一意性を確認するために 2 段階で行う．
断層面推定は板場ほか[2012]の手法を用いた．フィリピン海プレート境界面上[弘瀬ほか， 2007]に多数の断層面を仮定してグリッドサーチにより推定する．仮定した断層面上のすべりによって各観測点で期待される歪変化の計算には Okada [1992]のプログラムを用いる．1 段階目には，断層面のサイズは固定（幅・長さ共に 20 km），断層面の位置（0.1°間隔）およびすべり量（1～100 mm の間で 1 mm 間隔）のみ可変として広範囲で計算を行う．1 段階目の結果を示す図では，それぞれの断層面において最適なすべり量を与えたときの，観測値と計算値（期待値）との残差分布を示している．これにより，短期的 SSE が生じている可能性が高い領域を絞り込むとともに，推定された結果の任意性を確認することが出来る．2 段階目には，１段階目で絞り込んだ領域（＝残差が小さい領域）付近で，位置及びすべり量に加えて，断層面の長さを 10～80 km，幅を 10～50 km，それぞれ 1 km 間隔で可変として計算を行なう．その結果，観測値との残差が最小となる断層面が 1 つ計算されるが，計算に使用している観測点数が少ない場合や，断層面と観測点配置の関係によっては任意性が高くなるので注意が必要である．なお，異種観測値を統合して解析するため，観測点ごとに残差をノイズレベルによって規格化している．ノイズレベルは，気圧応答，潮汐成分およびホ
ワイトノイズ成分を取り除いた後（微動活動が活発な期間および周辺の日雨量 50 mm を超える時
期を除く）の 24 時間階差の 2σとした．
深部低周波微動の検出・震源決定には，エンベロープ相関法を用いている．
（落唯史・板場智史・小泉尚嗣・松本則夫・北川有一・武田直人・木村尚紀・木村武志・松澤孝
紀・汐見勝彦）－ 257 －

謝辞
短期的 SSE の断層モデル推定には，気象庁，静岡県の多成分歪計および体積歪計の記録を使用しました．気象庁の歪計データを解析する際には，気象庁によるキャリブレーション係数を使用しました．微動の解析には，気象庁，東京大学，京都大学，名古屋大学，高知大学，九州大学の地震波形記録を使用しました．低周波地震の震央位置表示には，気象庁の一元化カタログを使用しました．ここに記して感謝します．
参考文献
弘瀬冬樹, 中島淳一, 長谷川昭 (2007), Double-Difference Tomography 法による西南日本の 3 次元地震波速度構造およびフィリピン海プレートの形状の推定, 地震 2, 60, 1-20.
板場智史, 松本則夫, 北川有一, 小泉尚嗣, 松澤孝紀, 歪・傾斜・地下水統合解析による短期的スロースリップイベントのモニタリング, 日本地球惑星連合 2012 年大会, 千葉, 5 月, 2012.
Matsumoto, K., T. Sato, T. Takanezawa, and M. Ooe, GOTIC2: A Program for Computation of Oceanic
Tidal Loading Effect, J. Geod. Soc. Japan, 47, 243-248, 2001.
Matsumoto, N., O. Kamigaichi, Y. Kitagawa, S. Itaba, and N. Koizumi (2010), In-situ Calibration of
Borehole Strainmeter Using Green’s Functions for Surface Point Load at a Depth of Deployment, Eos,Trans. AGU, Abstract G11A-0626.
Okada, Y. (1992), Internal deformation due to shear and tensile faults in a half-space, Bull. Seismol. Soc.
Am., 82, 1018-1040.
Tamura, Y., T. Sato, M. Ooe and M. Ishiguro (1991), A procedure for tidal analysis with a Bayesian
information criterion, Geophys. J. Int., 104, 507-516.

－ 258 －

136ņ

137ņ

27

ņ

04

11

RK3

35

NSZ

RK4ANO

03E)
1, 2, 3, 4 1(N168E, N213E, N303E)

TYE

1E)
1, 2, 3, 4 (N356E, N86E, N131E, N221E)

NSZ1

5E)
1, 2, 3, 4 (N290E, N20E, N65E, N155E)

2.0x
10ŧ
[rad]

6E)
1, 2, 3, 4 (N141E, N231E, 1(, N6E)

N, E (+: down)

30[mm/h]
0

2.0x
10ŧ
[rad]

3.0x
10ŧ

2.0x
10ŧ
[rad]
2.0x
10ŧ
[rad]

0E, N30E)
1, 2, 3, 4 (N165E, N255E, N300E,

50
[/h]
0

3.0x
10ŧ

N, E (+: down)

N, E (+: down)

N, E (+: down)$68+ +LŧQHW

2.0x
10ŧ
[rad]

1, 2, 3, 4 (N65E, N155E, N200E,
E, N290E)

3.0x
10ŧ

27

N, E (+: down)

ITA1$6++ +LŧQHW

2.0x
10ŧ
[rad]

5, 6, , 8 (N65E, N155E, N200E,
E, N290E)

50 km

ANO2

7+1+ +LŧQHW

2.0x
10ŧ
[rad]

3E)
1, 2, 3, 4 (N38E, N128E, 1(, N263E)

3.0x
10ŧ

34

ANO2

2.0x
10ŧ
[rad]

ANO1

ICU

HGM

ANO1

1.0x
10ŧ

3.0x
10ŧ

ICU1

ņ

0.2
[m]

3.0x
10ŧ

MYM

ITA

30
1010
[mm/h] [hPa]
0
960

3.0x
10ŧ

ITA

RK2

TYS1

3.0x
10ŧ

TYS

RT0

ņ

ITA

1010
[hPa]
960

RT1

*12+ +LŧQHW

N, E (+: down):$7+ +LŧQHW

N, E (+: down).57+ +LŧQHW

N, E (+: down)

RK3, RK4, RT0, RT1
25

26

2

28

29

30

01

02

03

0E, N30E)
5, 6, , 8 (N165E, N255E, N300E,

04

05

0 08
06 0
015/05
2015/05

RK3, RK4, RT0, RT1
25

26

2

28

29

30

01

02

03

04

05

04

137

138

05/14

11

12

09

10

11

12

13

50
[/h]
0

第 2 図 紀伊半島〜東海地方における歪・傾斜・地下水観測結果(2015/4/25 00:00 - 2015/5/14 00:00 (JST))
Fig. 2
Observed strain, tilt and groundwater changes at the observationsites and number of deep low frequency tremors in Kii
Peninsula and the Tokai district from 0:00 April 25 to 0:00
May 14, 2015.

2015/5/7-10
11
136
04/25

06 0 08015/05
2015/0510

2015/5/7-10

3.0x
10ŧ
50
[/h]
0

09

Date

第 1 図 2015/4/25 00:00-5/14 00:00 の紀伊半島〜東海地方における深部低周波微動の時空間分布図．エンベロープ相関法
により決定．
Fig. 1
Epicentral and space-time distributions of deep low frequency
tremors in Kii Peninsula and the Tokai district from 0:00 April
25 to 0:00 May 14, 2015.
2015/5/7-10

132ņ

(a)

ASHH
ASUH

35.0°

GNOH

NSZ

ANO
34.5°

ITA

RS4

ņ

2.5

TYS

TYE

2.0

WATH

KOC

33

1.5

50 km

Slip: 25 mm, Mw: 5.7
136.0°

136.5°

137.0°

137.5°

138.0°(b1)

35.5°
tremor (AIST)
LFE (JMA)

1D

2

35.0°

1A

GNOH
ANO

34.0°

(b3)

Obs.
ASHH
ASUH

1B

ITA

NSZ

3C

Calc.

NSZ

1E
TYS
3B

TYE

11

25

137.0°

137.5°

ICU

138.0°

[B1]
[B2]
[B3]

01

ANO

Tilt
1x10ŧ rad

[A]

ŧ

TYE

WATH

136.5°

0

TYS

3A
THNH

KRTH

136.0°x10ŧ
ANO1ANO2

18

132
05/01

ŧ

Obs.
Calc.
ICU
Lat. 34.97° Lon. 136.93° Dep. 26 km Len. 39 km Wid. 34 km
Strike 203° Dip 13° Rake 78° Slip 9 mm Mw 5.7

135.5°

04

ANO1,2

(b2)

1F

34

SSK

UWA

ICU

1C

01

TSS

34.0°

34.5°

25

ņ

KRTH

135.5°

18

NHK

RS2
RS1

THNH

11

MAT

RS3

residual

tremor (AIST)
LFE (JMA)

04

ņ

35.5°

ŧ

2x10ŧ strain
Expansion

133
Date

Obs.
Calc.

Contraction

第 3 図 2015/5/7-10 の歪・傾斜・地下水変化（第 2 図網掛）を説明する断層モデル．地図中の太丸印は解析に使用した観測点（黒：産総研，緑：気象庁，青：防災科研 Hi-net）を示す．小さい青色丸印
は産総研による深部低周波微動，小さい緑色丸印は気象庁一元化カタログによる低周波イベントの震央位置を示す．(a) プレート境界面に沿って分布させた 20x20 km（大きさ固定）の矩形断層面に対して，おのおのの断層面で最適なすべり量を選んだときの残差の総和をその断層面の位置に示した分布図．赤色矩形が残差の総和が最小となる断層面．(b1) (a)の断層面付近でのグリッド
サーチにより推定された断層面（赤色矩形）と断層パラメータ．矢印は傾斜の観測値と計算値の比較．灰色矩形は最近周辺で発生した短期的 SSE の推定断層面．
（1A: 2015/1/1-3 (Mw 5.6), 1B:
2015/1/3-5 (Mw 5.7), 1C: 2015/ 1/7-10 (Mw 5.8), 1D: 2015/1/12-14 (Mw 5.6), 1E: 2015/1/15−18 (Mw
5.2), 1F: 2015/1/19−21 (Mw 5.8), 2: 2015/04/04-08 (Mw 5.8), 3A: 2015/4/14-15 (Mw 5.7), 3B:
2015/4/16-17 (Mw 5.7), 3C: 2015/4/18 (Mw 5.3)）．(b2) 主歪の観測値と(b1)に示した断層モデルから求めた計算値との比較．(b3) 体積歪の観測値と(b1)に示した断層モデルから求めた計算値との比
較．
Fig. 3 The fault model for the strain, tilt and groundwater changes from May 7 to May 10, 2015. The black,green and blue circles show the observation sites of AIST, JMA and NIED Hi-net. The blue and green
small circles show hypocenters of deep low frequency tremors estimated by AIST and low frequency
earthquakes estimated by JMA, respectively. (a) The distribution of the sum of residuals between
observed and calculated values that are caused by the uniform slip on a rectangular fault on the plateinterface. The length and width of the rectangular fault are fixed as 20 km and 20 km, and only the slip
amount that minimizes the sum of residuals is estimated. The red rectangle shows the fault patch withminimum residual. (b1) The red rectangle shows the estimated fault model. The gray rectangles show
the fault models of the recent events (1A: 2015/1/1-3 (Mw 5.6), 1B: 2015/1/3-5 (Mw 5.7), 1C: 2015/
1/7-10 (Mw 5.8), 1D: 2015/1/12-14 (Mw 5.6), 1E: 2015/1/15−18 (Mw 5.2), 1F: 2015/1/19−21 (Mw 5.8), 2:
2015/04/04-08 (Mw 5.8), 3A: 2015/4/14-15 (Mw 5.7), 3B: 2015/4/16-17 (Mw 5.7), 3C: 2015/4/18 (Mw5.3)). The observed and calculated tilt changes are also shown on the map by the red and the blue arrows.
(b2) The observed and calculated principal strain changes. (b3) The observed and calculated volumetricstrain changes.

06/07[A] 2015/5/14-16[B1] 2015/5/26-29AM[B2] 2015/5/29PM-31AM[B3] 2015/5/31PM-6/2

第 4 図 2015/5/1 00:00-6/7 00:00 の四国地方における深部低周波
微動の時空間分布図．エンベロープ相関法により決定．
Fig. 4
Epicentral and space-time distributions of deep low frequency
tremors in the Shikoku district from 0:00 May 1 to 0:00 June 7,2015.

－ 259 －

13

50
[/h]
0

UWA

1010
[hPa]
960

UWA1

[A] 2015/5/14-1650[mm/h]
0

1, 2, 3, 4 (N222E,
N312E, 1(,, 1(
2
(N2

(a)

tremor (AIST)
LFE (JMA)

1.10

34.0°

2.0x
10ŧ

1.05

OOZH
IKTH

33.5°

1, 2, 3 (N355E,(N355E N85E, 1(

UWAH
UWA

2.0x
10ŧ

YNDH
HIYH

residual

SSK1

SSK

1.00

33.0°
0.95

N, E

22=+ +LŧQHW

3.0x
10ŧ
[rad]

Slip: 14 mm, Mw: 5.5

32.5°

132.5°

133.0°(b2)
tremor (AIST)
LFE (JMA)

3.0x
10ŧ
[rad]
3.0x
10ŧ
[rad]

1B

2B

33.5°

N,, E

8:$+ +LŧQHW

1A

34.0°

N, E

,.7+ +LŧQHW

0.90

133.5°(b1)

N, E+,<+ +LŧQHW

3.0x
10ŧ
[rad]

132.0°

1C

UWAH
UWA

Calc.

x10ŧ
UWA1

ŧ

2x10ŧ strain

YNDH
HIYH

(b3)

0

UWA

SSK

OOZH
IKTH

Obs.

Expansion

SSK

Contraction

N, E<1'+ +LŧQHW

33.0°

UWA1

32.5°

0.2
[m]
RS0, RS1, RS2, RS3
3

50
[/h]
0

04

11

18

[A] 2015/05

25
5

01

50
[/h]
0

[B1] [B2] [B3]

第 5 図 四国地方における歪・傾斜・地下水観測結果 (2015/5/1 00:00 - 2015/6/7 00:00 (JST))
Fig. 5 Observed strain, tilt and groundwater changes at the observation sites and number of deep low frequency
tremors in the Shikoku district from 0:00 May 1 to 0:00 June 7, 2015.

Tilt
1x10ŧ rad/DWr/RQr'HSNPLen. 17 km Wid. 17 km
6WULNHr'LSr5DNHrSlip 19 mm Mw 5.5
132.0°

132.5°

133.0°

Obs.
Calc.

133.5°[B2] 2015/5/29PM-31AM

(a)

(a)

tremor (AIST)
LFE (JMA)

OOZH

UWA

34.0°

HIYH

SSK

2.5

OOZH

2.5

YNDH

IKTH

33.5°

residual

UWAH

3.0

tremor (AIST)
LFE (JMA)

3.0

34.0°

IKTH

Obs.
Calc.

第 6 図 2015/5/14-16 の歪・傾斜・地下水変化（第 5 図網掛[A]）を説明する断層モデル．地図中の太丸印
は解析に使用した観測点（黒：産総研，青：防災科研 Hi-net）を示す．小さい青色丸印は産総研による深部低周波微動，小さい緑色丸印は気象庁一元化カタログによる低周波イベントの震央位置を示す．(a) プレート境界面に沿って分布させた 20x20 km（大きさ固定）の矩形断層面に対して，おのおのの断層面で最適なすべり量を選んだときの残差の総和をその断層面の位置に示した
分布図．赤色矩形が残差の総和が最小となる断層面．(b1) (a)の断層面付近でのグリッドサーチにより推定された断層面（赤色矩形）と断層パラメータ．矢印は傾斜の観測値と計算値の比較．灰
色矩形は最近周辺で発生した短期的 SSE の推定断層面．（1A: 2014/12/26-31 (Mw 6.1), 1B:
2014/12/31-2015/1/3 (Mw 5.8), 1C: 2015/1/10-12 (Mw 5.9), 2A: 2015/3/24PM-26AM (Mw 5.8), 2B:
2015/3/26PM-27 (Mw 5.3)）．(b2)主歪の観測値と(b1)に示した断層モデルから求めた計算値との比較．(b3)体積歪の観測値と(b1)に示した断層モデルから求めた計算値との比較．
Fig. 6 The fault model for the strain, tilt and groundwater changes from May 7 to May 10, 2015. The black andblue circles show the observation sites of AIST and NIED Hi-net. The blue and green small circles show
hypocenters of deep low frequency tremors estimated by AIST and low frequency earthquakes estimated by
JMA, respectively. (a) The distribution of the sum of residuals between observed and calculated values
that are caused by the uniform slip on a rectangular fault on the plate interface. The length and width of
the rectangular fault are fixed as 20 km and 20 km, and only the slip amount that minimizes the sum of
residuals is estimated. The red rectangle shows the fault patch with minimum residual. (b1) The redrectangle shows the estimated fault model. The gray rectangles show the fault models of the recent events(1A: 2014/12/26-31 (Mw 6.1), 1B: 2014/12/31-2015/1/3 (Mw 5.8), 1C: 2015/1/10-12 (Mw 5.9), 2A:
2015/3/24PM-26AM (Mw 5.8), 2B: 2015/3/26PM-27 (Mw 5.3)). The observed and calculated tiltchanges are also shown on the map by the red and the blue arrows. (b2) The observed and calculated
principal strain changes. (b3) The observed and calculated volumetric strain changes.

[B1] 2015/5/26-29AM

33.5°

ŧ

2A

UWAH
UWA

YNDH
HIYH

residual

3.0x
10ŧ
[rad]

SSK

2.0

2.0

33.0°

33.0°
1.5

Slip: 12 mm, Mw: 5.5

32.5°

132.0°

132.5°

133.0°(b1)(b2)
tremor (AIST)
LFE (JMA)

1A

34.0°

2B
1C

UWAH
UWA

YNDH
HIYH

SSK

32.5°

Obs.
Calc.

132.5°

133.0°

133.5°

(b2)
1A

34.0°

1B

2B

33.5°

1C

Expansion

B1
IKTH

UWAH
UWA

33.0°

Tilt
1x10ŧ rad

32.5°

133.5°

第 7 図 2015/5/26-29AM の歪・傾斜・地下水変化（第 5 図網掛[B1]）を説明する断層モデル．地図中の太丸印は解析に使用した観測点（黒：産総研，青：防災科研 Hi-net）を示す．小さい青色丸印は産総研による深部低周波微動，小さい緑色丸印は気象庁一元化カタログによる低周波イベントの震央位置を示す．(a)プレート境界面に沿って分布させた 20x20 km（大きさ固定）の矩形断層面に対して，おのおのの断層面で最適なすべり量を選んだときの残差の総和をその断層面の位置に示した分布図．赤色矩形が残差の総和が最小となる断層面．(b1)(a)の断層面付近でのグリッドサーチ
により推定された断層面（赤色矩形）と断層パラメータ．矢印は傾斜の観測値と計算値の比較．
灰色矩形は最近周辺で発生した短期的 SSE の推定断層面．（1A: 2014/12/26-31 (Mw 6.1), 1B:
2014/12/31-2015/1/3 (Mw 5.8), 1C: 2015/1/10-12 (Mw 5.9), 2A: 2015/3/24PM-26AM (Mw 5.8), 2B:
2015/3/26PM-27 (Mw 5.3), A: 2015/5/14-16 (Mw 5.5)）．(b2)主歪の観測値と(b1)に示した断層モデルから求めた計算値との比較．(b3)体積歪の観測値と(b1)に示した断層モデルから求めた計算値との比較．
Fig. 7 The fault model for the strain, tilt and groundwater changes from May 26 to May 29 AM, 2015. The blackand blue circles show the observation sites of AIST and NIED Hi-net. The blue and green small circles
show hypocenters of deep low frequency tremors estimated by AIST and low frequency earthquakes
estimated by JMA, respectively. (a) The distribution of the sum of residuals between observed andcalculated values that are caused by the uniform slip on a rectangular fault on the plate interface. The
length and width of the rectangular fault are fixed as 20 km and 20 km, and only the slip amount thatminimizes the sum of residuals is estimated. The red rectangle shows the fault patch with minimum
residual. (b1) The red rectangle shows the estimated fault model. The gray rectangles show the fault
models of the recent events (1A: 2014/12/26-31 (Mw 6.1), 1B: 2014/12/31-2015/1/3 (Mw 5.8), 1C:
2015/1/10-12 (Mw 5.9), 2A: 2015/3/24PM-26AM (Mw 5.8), 2B: 2015/3/26PM-27 (Mw 5.3), A:
2015/5/14-16 (Mw 5.5)). The observed and calculated tilt changes are also shown on the map by the red
and the blue arrows. (b2) The observed and calculated principal strain changes. (b3) The observed andcalculated volumetric strain changes.

Obs.

Calc.

(b3)

UWA

SSK

OOZH

2x10ŧ strain

YNDH
HIYH

SSK

Expansion
Contraction

A/DWr/RQr'HSNPLen. 39 km Wid. 39 km
6WULNHr'LSr5DNHrSlip 6 mm Mw 5.6
132.0°

133.0°

tremor (AIST)
LFE (JMA)

Contraction

2A

132.5°(b1)

2x10ŧ strain

A
33.0°

(b3)

132.0°

SSK

OOZH
IKTH

Calc.

UWA

1B

33.5°

Obs.

1.5

Slip: 20 mm, Mw: 5.6

32.5°

133.5°

2A

Obs.
Calc.

Tilt
1x10ŧ rad/DWr/RQr'HSNPLen. 59 km Wid. 31 km
6WULNHr'LSr5DNHrSlip 11 mm Mw 5.9
132.0°

132.5°

133.0°

133.5°

第 8 図 2015/5/29PM-31AM の歪・傾斜・地下水変化（第 5 図網掛[B2]）を説明する断層モデル．地図中の太丸印は解析に使用した観測点（黒：産総研，青：防災科研 Hi-net）を示す．小さい青色丸印は産総研による深部低周波微動，小さい緑色丸印は気象庁一元化カタログによる低周波イベントの震央位置を示す．(a)プレート境界面に沿って分布させた 20x20 km（大きさ固定）の矩形断層面に対して，おのおのの断層面で最適なすべり量を選んだときの残差の総和をその断層面の位置に示した分
布図．赤色矩形が残差の総和が最小となる断層面．(b1)(a)の断層面付近でのグリッドサーチにより推定された断層面（赤色矩形）と断層パラメータ．矢印は傾斜の観測値と計算値の比較．灰色矩形は最近周辺で発生した短期的 SSE の推定断層面．
（1A: 2014/12/26-31 (Mw 6.1), 1B: 2014/12/31-2015/1/3 (Mw
5.8), 1C: 2015/1/10-12 (Mw 5.9), 2A: 2015/3/24PM-26AM (Mw 5.8), 2B: 2015/3/26PM-27 (Mw 5.3), A:
2015/5/14-16 (Mw 5.5), B1: 2015/5/26-29AM (Mw 5.6)）．(b2)主歪の観測値と(b1)に示した断層モデルから求めた計算値との比較．(b3)体積歪の観測値と(b1)に示した断層モデルから求めた計算値との比較．
Fig. 8 The fault model for the strain and tilt changes from May 29 PM to May 31 AM, 2015. The black and blue
circles show the observation sites of AIST and NIED Hi-net. The blue and green small circles show
hypocenters of deep low frequency tremors estimated by AIST and low frequency earthquakes estimated by
JMA, respectively. (a) The distribution of the sum of residuals between observed and calculated values
that are caused by the uniform slip on a rectangular fault on the plate interface. The length and width of
the rectangular fault are fixed as 20 km and 20 km, and only the slip amount that minimizes the sum of
residuals is estimated. The red rectangle shows the fault patch with minimum residual. (b1) The redrectangle shows the estimated fault model. The gray rectangles show the fault models of the recent events(1A: 2014/12/26-31 (Mw 6.1), 1B: 2014/12/31-2015/1/3 (Mw 5.8), 1C: 2015/1/10-12 (Mw 5.9), 2A:
2015/3/24PM-26AM (Mw 5.8), 2B: 2015/3/26PM-27 (Mw 5.3), A: 2015/5/14-16 (Mw 5.5), B1:
2015/5/26-29AM (Mw 5.6)). The observed and calculated tilt changes are also shown on the map by the
red and the blue arrows. (b2) The observed and calculated principal strain changes. (b3) The observedand calculated volumetric strain changes.

－ 260 －

135ņ[B3] 2015/5/31PM-6/2
(a)

136ņ

25

ņ

01

08

RK3
tremor (AIST)
LFE (JMA)

RK1

1.8

OOZH
UWAH
UWA

YNDH
HIYH

1.6

SSK

ņ

residual

IKTH

33.5°

ITA

RK2

2.0

34.0°

MYM

RK0

HGM

1.4

34

ICU

KST

33.0°
1.2

50 km

Slip: 9 mm, Mw: 5.4

32.5°

132.0°

132.5°

133.0°

133.5°(b1)(b2)
tremor (AIST)
LFE (JMA)

1A

34.0°

Obs.

Calc.

(b3)

x10ŧ
UWA1

25

0

UWA

1B
B2
2B

33.5°

B1
IKTH

1C

01

SSK

OOZH
UWAH
UWA

2015/6/4-7

2x10ŧ strain

YNDH
HIYH

ŧ

08

Expansion

SSK

Contraction

A

ŧ
33.0°

32.5°

2A

Obs.
Calc.

Tilt
1x10ŧ rad135

Lat. 33.60° Lon. 133.15° Dep. 27 km Len. 80 km Wid. 26 km
Strike 228° Dip 12° Rake 103° Slip 5 mm Mw 5.7
132.0°

132.5°

133.0°

133.5°

ŧ

1030
980

ICU
ICU1

第 10 図 2015/5/20 00:00-6/15 00:00 の紀伊半島における深部低周波微動の時空間分布図．エンベロープ相関法により決定．
Fig. 10
Epicentral and space-time distributions of deep low
frequency tremors in Kii Peninsula from 0:00 May 20 to 0:00
June 15, 2015.

2015/6/4-7

40

(a)

[mm/h]

1, 2, 3, 4 (N141E, N231E,
1()
31E, 1(, 1

137

06/15

Date

Obs.
Calc.

第 9 図 2015/5/31PM-6/2 の歪・傾斜・地下水変化（第 5 図網掛[B3]）を説明する断層モデル．地図中の太丸印は解析に使用した観測点（黒：産総研，青：防災科研 Hi-net）を示す．小さい青色丸印は産総研による深部低周波微動，小さい緑色丸印は気象庁一元化カタログによる低周波イベントの震央位置を示す．(a)プレート境界面に沿って分布させた 20x20 km（大きさ固定）の矩形断層面に対して，おのおのの断層面で最適なすべり量を選んだときの残差の総和をその断層面の位置に示した分布図．赤色矩形が残差の総和が最小となる断層面．(b1)(a)の断層面付近でのグリッドサーチ
により推定された断層面（赤色矩形）と断層パラメータ．矢印は傾斜の観測値と計算値の比較．
灰色矩形は最近周辺で発生した短期的 SSE の推定断層面（1A: 2014/12/26-31 (Mw 6.1), 1B:
2014/12/31-2015/1/3 (Mw 5.8), 1C: 2015/1/10-12 (Mw 5.9), 2A: 2015/3/24PM-26AM (Mw 5.8), 2B:
2015/3/26PM-27 (Mw 5.3), A: 2015/5/14-16 (Mw 5.5), B1: 2015/5/26-29AM (Mw 5.6) , B2:
2015/5/29PM-31AM (Mw 5.9)）．(b2)主歪の観測値と(b1)に示した断層モデルから求めた計算値との比較．(b3)体積歪の観測値と(b1)に示した断層モデルから求めた計算値との比較．
Fig. 9 The fault model for the strain, tilt and groundwater changes from May 31 PM to June 2, 2015. The blackand blue circles show the observation sites of AIST and NIED Hi-net. The blue and green small circles
show hypocenters of deep low frequency tremors estimated by AIST and low frequency earthquakes
estimated by JMA, respectively. (a) The distribution of the sum of residuals between observed andcalculated values that are caused by the uniform slip on a rectangular fault on the plate interface. The
length and width of the rectangular fault are fixed as 20 km and 20 km, and only the slip amount thatminimizes the sum of residuals is estimated. The red rectangle shows the fault patch with minimum
residual. (b1) The red rectangle shows the estimated fault model. The gray rectangles show the fault
models of the recent events (1A: 2014/12/26-31 (Mw 6.1), 1B: 2014/12/31-2015/1/3 (Mw 5.8), 1C:
2015/1/10-12 (Mw 5.9), 2A: 2015/3/24PM-26AM (Mw 5.8), 2B: 2015/3/26PM-27 (Mw 5.3), A:
2015/5/14-16 (Mw 5.5), B1: 2015/5/26-29AM (Mw 5.6) , B2: 2015/5/29PM-31AM (Mw 5.9)). The
observed and calculated tilt changes are also shown on the map by the red and the blue arrows. (b2) The
observed and calculated principal strain changes. (b3) The observed and calculated volumetric strain changes.

[hPa]

136

05/20

0

tremor (AIST)
LFE (JMA)

3.5

34.5°

2.5
34.0°

HGM

HGM1

1, 2, 3, 4 1(1(,(, N112E, N202E)

ICU

2.0

1.5

33.5°

KST
Slip: 48 mm, Mw: 5.9
135.0°

5.0x
10ŧ

135.5°

136.0°

136.5°

137.0°

(b1)

34.0°

Calc.

HGM

1

2.0x
10ŧ

Obs.

ICU

34.5°

1, 2, 3, 4 (N310E, N40E,
1()
0E, N85E, 1

1.0

(b2)
tremor (AIST)
LFE (JMA)

KST1

residual

3.0

2.0x
10ŧ

2
HGM

KST

ICU

2x10ŧ strain
Expansion

33.5°

40

[/h]

0

RK0, RK1, RK2, RK340

Lat. 34.13° Lon. 135.97° Dep. 35 km Len. 40 km Wid. 10 km
Strike 233° Dip 23° Rake 108° Slip 27 mm Mw 5.7[/h]

20 21 22 23 24 25 2 2 28 29 30 31 01 02 03 04 05 0 0 08 09 10 11 12 13 14
2015/05
2015/0

02015/6/4-7

第 11 図 紀伊半島における歪観測結果 (2015/5/20 00:00 2015/6/15 00:00 (JST))
Fig. 11 Observed strain changes at the observation sites and numberof deep low frequency tremors in Kii Peninsula from 0:00
May 20 to 0:00 June 15, 2015.

Contraction

KST

135.0°

135.5°

136.0°

136.5°

137.0°

第 12 図 2015/6/4-7 の歪変化（第 11 図網掛）を説明する断層モデル．地図中の太丸印は解析に使用した
産総研の観測点を示す．小さい青色丸印は産総研による深部低周波微動，小さい緑色丸印は気象庁一元化カタログによる低周波イベントの震央位置を示す．(a)プレート境界面に沿って分布させた 20x20 km（大きさ固定）の矩形断層面に対して，おのおのの断層面で最適なすべり量を選んだときの残差の総和をその断層面の位置に示した分布図．赤色矩形が残差の総和が最小となる断層面．(b1)(a)の断層面付近でのグリッドサーチにより推定された断層面（赤色矩形）と断層パ
ラメータ．灰色矩形は最近周辺で発生した短期的 SSE の推定断層面（1: 2015/1/7-10 (Mw 5.8), 2:
2015/3/14-16 (Mw 5.3)）．(b2)主歪の観測値と(b1)に示した断層モデルから求めた計算値との比較．
Fig. 12 The fault model for the strain changes from June 4 to June 7, 2015. The black circles show theobservation sites of AIST. The blue and green small circles show hypocenters of deep low frequency
tremors estimated by AIST and low frequency earthquakes estimated by JMA, respectively. (a) The
distribution of the sum of residuals between observed and calculated values that are caused by the uniform
slip on a rectangular fault on the plate interface. The length and width of the rectangular fault are fixed as
20 km and 20 km, and only the slip amount that minimizes the sum of residuals is estimated. The redrectangle shows the fault patch with minimum residual. The black circles show the observation sites of
AIST. The blue and green small circles show hypocenters of deep low frequency tremors estimated by
AIST and low frequency earthquakes estimated by JMA, respectively. (b1) The red rectangle shows theestimated fault model. The gray rectangles show the fault models of the recent events (1: 2015/1/7-10(Mw 5.8), 2: 2015/3/14-16 (Mw 5.3)). The observed and calculated tilt changes are also shown on themap by the red and the blue arrows. (b2) The observed and calculated principal strain changes.

－ 261 －

136ņ

06

ņ

13

20

27

TYS

RT0

ņ

35

NSZ

RK4 ANO

ITA

1000
[hPa]
950

RT1

ANO1

1, 2, 3, 4 (N38E, N128E, 1(,
N26
63E)

(, N263E)
6

NSZ1

1, 2, 3, 4 (N290E, N20E, N65E,
N15
5
5
55E)
5E, N155E)

2.0x
10ŧ

RK3

ITA

RK2

2.0x
10ŧ

MYM
ņ

1, 2, 3, 4 (N65E, N155E,
N200E, N290E)
E,
E N200E,

34

ICU

HGM

40
[mm/h]
0

2.0x
10ŧ

50 km

5, 6, , 8 (N65E, N155E,
N200E, N290E)
E,
E N200E,

06

2.0x
10ŧ

13
2.0x
10ŧ
[rad]

[A1]

20

[A2]

27

136
07/01

137
07/29

N, E

ITA1

N, E

2.0x
10ŧ
[rad]
2.0x
10ŧ
[rad]

[A1] 2015/7/18-19[A2] 2015/7/20-25

Date

ANO2

2.0x
10ŧ
[rad]

第 13 図 2015/7/1 00:00-7/29 00:00 の紀伊半島〜東海地方における深部低周波微動の時空間分布図．エンベロープ相関法
により決定．
Fig. 13
Epicentral and space-time distributions of deep low
frequency tremors in Kii Peninsula and the Tokai district
from 0:00 July 1 to 0:00 July 29, 2015.

2.0x
10ŧ
[rad]
50
[/h]
0

856+ +LŧQHW

N, E*12+ +LŧQHW

N, E:$7+ +LŧQHW

N, E

RK2, RK3, RK4, RT0
01 02 03 04 05 06 0 08 09 10 11 12 13 14 15 16 1 1
18
8 19 220
0 21 2
22
22
23
3 24
24 25
25 26 2 28
[A1]
[A2]

50
[/h]
0

第 14 図 紀伊半島〜東海地方における歪・傾斜観測結果 (2015/7/1 00:00 - 2015/7/29 00:00 (JST))
Fig. 14 Observed strain and tilt changes at the observation sites and number of deep low frequency tremors in Kii
Peninsula and the Tokai district from 0:00 July 1 to 0:00 July 29, 2015.

[A2] 2015/7/20-25[A1] 2015/7/18-19
(a)

(a)

tremor (AIST)
LFE (JMA)

tremor (AIST)
LFE (JMA)

0.9

35.0°

1.8

35.0°

NSZ

GNOH

ANO

URSH
WATH

1.6

residual

0.7

URSH

34.5°

ITA

WATH

residual

ANO

34.5°

NSZ

GNOH

0.8

0.6
1.4
34.0°

34.0°

0.5

Slip: 10 mm, Mw: 5.4
135.5°

136.0°(b1)

136.5°

137.0°

Slip: 11 mm, Mw: 5.4
137.5°

tremor (AIST)
LFE (JMA)

GNOH

1A

2A

1F

Obs.

Calc.

35.0°

1D

2x10ŧ strain

1E

Expansion

A1

136.5°

137.0°

137.5°

(b2)

1C

NSZ

2B

1B
URSH
ITA

WATH

1D
1F

Contraction

1E

34.0°

Lat. 34.63° Lon. 136.76° Dep. 27 km Len. 60 km Wid. 41 km
Strike 225° Dip 14° Rake 100° Slip 2 mm Mw 5.5
136.0°

2A

ANO

34.5°

Tilt
1x10ŧ7 rad

137.0°

1A

GNOH
NSZ

WATH

Obs.
Calc.

136.5°

tremor (AIST)
LFE (JMA)

1B
URSH

34.0°

135.5°

136.0°

NSZ

2B

ANO

1C

135.5°

(b1)

ANO

35.0°

34.5°

(b2)

Obs.
Calc.

Tilt
1x10ŧ radLat. 34.86° Lon. 136.89° Dep. 28 km Len. 66 km Wid. 28 km
Strike 211° Dip 11° Rake 86° Slip 4 mm Mw 5.6

137.5°

135.5°

第 15 図 2015/7/18-19 の歪・傾斜変化（第 14 図網掛[A1]）を説明する断層モデル．地図中の太丸印は解析に使用した観測点（黒：産総研，緑：気象庁，青：防災科研 Hi-net）を示す．小さい青色丸印は産総研による深部低周波微動，小さい緑色丸印は気象庁一元化カタログによる低周波イベントの震央位置を示す．(a)プレート境界面に沿って分布させた 20x20 km（大きさ固定）の矩形断層面に対して，おのおのの断層面で最適なすべり量を選んだときの残差の総和をその断層面の位置に示した分布図．赤色矩形が残差の総和が最小となる断層面．(b1)(a)の断層面付近でのグリッドサーチにより推定された断層面（赤色矩形）と断層パラメータ．灰色矩形は最近周辺で発生した短期的 SSE の推定断層面（1A: 2014/7/3-4 (Mw 5.4), 1B: 2015/7/5-6 (Mw 5.6), 1C: 2015/7/7-8 (Mw 5.5),
1D: 2015/7/9-10 (Mw 5.8), 1E: 2014/7/11-12 (Mw 5.7), 1F: 2015/7/13-14 (Mw 5.7), 2A: 2015/1/1-3 (Mw 5.6),
2B: 2015/1/3-5 (Mw 5.7)）．(b2)主歪の観測値と(b1)に示した断層モデルから求めた計算値との比較．
Fig. 15 The fault model for the strain and tilt changes from July 18 to July 19, 2015. The black, green and blue
bold circles show the observation sites of AIST, JMA and NIED Hi-net, respectively. The blue and green
small circles show hypocenters of deep low frequency tremors estimated by AIST and low frequency
earthquakes estimated by JMA, respectively. (a) The distribution of the sum of residuals between
observed and calculated values that are caused by the uniform slip on a rectangular fault on the plateinterface. The length and width of the rectangular fault are fixed as 20 km and 20 km, and only the slip
amount that minimizes the sum of residuals is estimated. The red rectangle shows the fault patch withminimum residual. (b1) The red rectangle shows the estimated fault model. The gray rectangles show thefault models of the recent events (1A: 2014/7/3-4 (Mw 5.4), 1B: 2015/7/5-6 (Mw 5.6), 1C: 2015/7/7-8(Mw 5.5), 1D: 2015/7/9-10 (Mw 5.8), 1E: 2014/7/11-12 (Mw 5.7), 1F: 2015/7/13-14 (Mw 5.7), 2A:
2015/1/1-3 (Mw 5.6), 2B: 2015/1/3-5 (Mw 5.7)). The observed and calculated tilt changes are also
shown on the map by the red and blue arrows. (b2) The observed and calculated principal strain changes.

136.0°

136.5°

137.0°

137.5°

第 16 図 2015/7/20-25 の歪・傾斜変化（第 14 図網掛[A2]）を説明する断層モデル．地図中の太丸印は解析に使用した観測点（黒：産総研，緑：気象庁，青：防災科研 Hi-net）を示す．小さい青色丸印は産総研による深部低周波微動，小さい緑色丸印は気象庁一元化カタログによる低周波イベントの震央位置を示す．(a)プレート境界面に沿って分布させた 20x20 km（大きさ固定）の矩形断層面に対して，おのおのの断層面で最適なすべり量を選んだときの残差の総和をその断層面の位置に示した分布図．赤色矩形が残差の総和が最小となる断層面．(b1)(a)の断層面付近でのグリッドサーチにより推定された断層面（赤色矩形）と断層パラメータ．灰色矩形は最近周辺で発生した短期的 SSE の推定断層面（1A: 2014/7/3-4 (Mw 5.4), 1B: 2015/7/5-6 (Mw 5.6), 1C: 2015/7/7-8 (Mw 5.5),
1D: 2015/7/9-10 (Mw 5.8), 1E: 2014/7/11-12 (Mw 5.7), 1F: 2015/7/13-14 (Mw 5.7), 2A: 2015/1/1-3 (Mw
5.6), 2B: 2015/1/3-5 (Mw 5.7), A1: 2015/7/18-19 (Mw 5.5)）．(b2)主歪の観測値と(b1)に示した断層モデルから求めた計算値との比較．
Fig. 16 The fault model for the strain and tilt changes from July 20 to July 25, 2015. The black, green and blue
bold circles show the observation sites of AIST, JMA and NIED Hi-net, respectively. The blue and green
small circles show hypocenters of deep low frequency tremors estimated by AIST and low frequency
earthquakes estimated by JMA, respectively. (a) The distribution of the sum of residuals between
observed and calculated values that are caused by the uniform slip on a rectangular fault on the plateinterface. The length and width of the rectangular fault are fixed as 20 km and 20 km, and only the slip
amount that minimizes the sum of residuals is estimated. The red rectangle shows the fault patch withminimum residual. (b1) The red rectangle shows the estimated fault model. The gray rectangles show thefault models of the recent events (1A: 2014/7/3-4 (Mw 5.4), 1B: 2015/7/5-6 (Mw 5.6), 1C: 2015/7/7-8(Mw 5.5), 1D: 2015/7/9-10 (Mw 5.8), 1E: 2014/7/11-12 (Mw 5.7), 1F: 2015/7/13-14 (Mw 5.7), 2A:
2015/1/1-3 (Mw 5.6), 2B: 2015/1/3-5 (Mw 5.7), A1: 2015/7/18-19 (Mw 5.5)). The observed andcalculated tilt changes are also shown on the map by the red and blue arrows. (b2) The observed andcalculated principal strain changes.

－ 262 －

132ņ

24

ņ
RS4

ņ

07

NHK

RS2

14
34

MAT

RS3
RS1

31

UWA

1010
[hPa]
960


(,
312E, 13
1, 2, 3, 4 (N222E, N312E,
1(,
1(

UWA1

2.0x
10ŧ

KOC
SSK

UWA

100
[mm/h]
0

TSS1

33

ņ

E, N21E
E 1(1, 2, 3, 4 (N246E, N336E
N336E,
N21E,

2.0x
10ŧ

TSS
0

0E, 1
0E
55E,, N130
N131, 2, 3, 4 (N355E, N85E,
N130E,
1(

SSK1

50 km
2.0x
10ŧ

24

2.0x
10ŧ
[rad]

31

[A]
[B]

07

132
08/16

133

2.0x
10ŧ
[rad]

09/16

Date

50
[/h]
0

[A] 2015/9/1 12:00 - 9/4 0:00
[B] 2015/9/4 0:00 - 9/6 0:00

第 17 図 2015/8/16 00:00-9/16 00:00 の四国地方における深部低周波微動の時空間分布図．エンベロープ相関法により決定．
Fig. 17 Epicentral and space-time distributions of deep low
frequency tremors in the Shikoku district from 0:00 August
16 to 0:00 September 16, 2015.

N, E

0,6+ +LŧQHW

N, E

22=+ +LŧQHW

N, E

RS0, RS1, RS2, RS304 05
0 06 0 08 09 10 11 12 13 14 15
1 02 03
03 04
16 1 18 19 20 21 22 23 24 25 26 2 28 29 30 31 01
2015/08[A] [B] 2015/09

50
[/h]
0

[A] 2015/9/1 12:00 - 9/4 0:00
[B] 2015/9/4 0:00 - 9/6 0:00

第 18 図 四国地方における歪・傾斜観測結果 (2015/8/16 00:00 2015/9/16 00:00 (JST))
Fig. 18 Observed strain and tilt changes at the observation sites andnumber of deep low frequency tremors in the Shikoku district
from 0:00 August 16 to 0:00 September 16, 2015.

[B] 2015/9/4-5[A] 2015/9/1PM-3

(a)

(a)

1.8

tremor (AIST)
LFE (JMA)

tremor (AIST)
LFE (JMA)

34.0°

34.0°
2.5

SSK

UWA

2.0

33.5°

residual

33.5°

1.6

UWA

33.0°

1.4

residual

14

2.0x
10ŧ
[rad]

8:$+ +LŧQHW

33.0°
1.2

TSS

1.5

Slip: 33 mm, Mw: 5.8

32.5°
132.0°(b1)

132.5°

133.0°

133.5°

tremor (AIST)
LFE (JMA)

Obs.

132.0°

(b1)

Calc.

132.5°

2B

34.0°

2A

33.0°

1A
Obs.
TSSCalc.

2A

SSK

Tilt
1x10ŧ rad

UWA

Calc.

SSK
2x10ŧ strain

1A

Expansion
Contraction

Obs.
Calc.

Lat. 33.28° Lon. 132.55° Dep. 28 km Len. 70 km Wid. 50 km
133.0°

Obs.

TSS

2D

33.0° A

Tilt
1x10ŧ rad

Expansion
Contraction

Lat. 33.35° Lon. 132.59° Dep. 31 km Len. 35 km Wid. 35 km

32.5° Strike 215° Dip 17° Rake 90° Slip 9 mm Mw 6.0
132.5°

2B

33.5°

SSK

UWA

1B

2x10ŧ strain

132.0°

(b2)

2C
TSS

2D

33.5°

133.5°

UWA

2C
1B

133.0°

tremor (AIST)
LFE (JMA)

UWA

34.0°

Slip: 11 mm, Mw: 5.4

32.5°

(b2)

32.5° Strike 224° Dip 13° Rake 99° Slip 5 mm Mw 5.5

133.5°

132.0°

第 19 図 2015/9/1PM-3 の歪・傾斜変化（第 18 図赤色[A]）を説明する断層モデル．地図中の太丸印は解析に使用した観測点（黒：産総研，青：防災科研 Hi-net）を示す．小さい青色丸印は産総研による
深部低周波微動，小さい緑色丸印は気象庁一元化カタログによる低周波イベントの震央位置を示す．(a)プレート境界面に沿って分布させた 20x20 km（大きさ固定）の矩形断層面に対して，おのおのの断層面で最適なすべり量を選んだときの残差の総和をその断層面の位置に示した分布
図．赤色矩形が残差の総和が最小となる断層面．(b1)(a)の断層面付近でのグリッドサーチにより推定された断層面（赤色矩形）と断層パラメータ．灰色矩形は最近周辺で発生した短期的 SSEの推定断層面（1A: 2015/3/24PM-26AM (Mw 5.8), 1B: 2015/3/26PM-27 (Mw 5.3), 2A: 2015/5/14-16(Mw 5.5), 2B: 2015/5/26-29AM (Mw 5.6), 2C: 2015/5/29PM-31AM (Mw 5.9), 2D: 2015/5/31PM-6/2(Mw 5.7)）．(b2)主歪の観測値と(b1)に示した断層モデルから求めた計算値との比較．
Fig. 19 The fault model for the strain and tilt changes from September 1 PM to September 3, 2015. The blackand blue bold circles show the observation sites of AIST and NIED Hi-net, respectively. The blue andgreen small circles show hypocenters of deep low frequency tremors estimated by AIST and low
frequency earthquakes estimated by JMA, respectively. (a) The distribution of the sum of residuals
between observed and calculated values that are caused by the uniform slip on a rectangular fault on theplate interface. The length and width of the rectangular fault are fixed as 20 km and 20 km, and only theslip amount that minimizes the sum of residuals is estimated. The red rectangle shows the fault patchwith minimum residual. (b1) The red rectangle shows the estimated fault model. The gray rectangles
show the fault models of the recent events (1A: 2015/3/24PM-26AM (Mw 5.8), 1B: 2015/3/26PM-27(Mw 5.3), 2A: 2015/5/14-16 (Mw 5.5), 2B: 2015/5/26-29AM (Mw 5.6), 2C: 2015/5/29PM-31AM (Mw
5.9), 2D: 2015/5/31PM-6/2 (Mw 5.7)). The observed and calculated tilt changes are also shown on themap by the red and blue arrows. (b2) The observed and calculated principal strain changes.

132.5°

133.0°

133.5°

第 20 図 2015/9/4-5 の歪・傾斜変化（第 18 図緑色[B]）を説明する断層モデル．地図中の太丸印は解析に使用した観測点（黒：産総研，青：防災科研 Hi-net）を示す．小さい青色丸印は産総研による深部低周波微動，小さい緑色丸印は気象庁一元化カタログによる低周波イベントの震央位置を示す．
(a)プレート境界面に沿って分布させた 20x20 km（大きさ固定）の矩形断層面に対して，おのおのの断層面で最適なすべり量を選んだときの残差の総和をその断層面の位置に示した分布図．赤色矩形が残差の総和が最小となる断層面．(b1)(a)の断層面付近でのグリッドサーチにより推定された断層面（赤色矩形）と断層パラメータ．灰色矩形は最近周辺で発生した短期的 SSE の推定
断層面（1A: 2015/3/24PM-26AM (Mw 5.8), 1B: 2015/3/26PM-27 (Mw 5.3), 2A: 2015/5/14-16 (Mw 5.5),
2B: 2015/5/26-29AM (Mw 5.6), 2C: 2015/5/29PM-31AM (Mw 5.9), 2D: 2015/5/31PM-6/2 (Mw 5.7), A:
2015/9/1PM-3 (Mw 6.0)）．(b2)主歪の観測値と(b1)に示した断層モデルから求めた計算値との比較．
Fig. 20 The fault model for the strain and tilt changes from September 4 to September 5, 2015. The black andblue bold circles show the observation sites of AIST and NIED Hi-net, respectively. The blue and green
small circles show hypocenters of deep low frequency tremors estimated by AIST and low frequency
earthquakes estimated by JMA, respectively. (a) The distribution of the sum of residuals between
observed and calculated values that are caused by the uniform slip on a rectangular fault on the plateinterface. The length and width of the rectangular fault are fixed as 20 km and 20 km, and only the slip
amount that minimizes the sum of residuals is estimated. The red rectangle shows the fault patch withminimum residual. (b1) The red rectangle shows the estimated fault model. The gray rectangles show thefault models of the recent events (1A: 2015/3/24PM-26AM (Mw 5.8), 1B: 2015/3/26PM-27 (Mw 5.3), 2A:
2015/5/14-16 (Mw 5.5), 2B: 2015/5/26-29AM (Mw 5.6), 2C: 2015/5/29PM-31AM (Mw 5.9), 2D:
2015/5/31PM-6/2 (Mw 5.7), A: 2015/9/1PM-3 (Mw 6.0)). The observed and calculated tilt changes arealso shown on the map by the red and blue arrows. (b2) The observed and calculated principal strainchanges.

－ 263 －

135°

136°

RK3

26

1030
[hPa]
980
1030
[hPa]
980

35

NSZ

RK4ANO

RK1

MYM
KST1

1, 2, 3, 4 (N310E, N40E, N85E, 1(HGM1

1, 2, 3, 4 1(1(, N112E, 1(ICU1

1, 2, 3, 4 (N141E, N231E, 1(, 1(ANO1

1, 2, 3, 4 (N38E, N128E, 1(, 1(1030
50[mm/h] [hPa]
980
0
1030
50
[mm/h] [hPa]
980
0
2.0 x
10ŧ
[rad]

2.0 x
10ŧ

2.0 x
10ŧ
[rad]
2.0 x
10ŧ
[rad]
2.0 x
10ŧ
[rad]

2.0 x
10ŧ

MYM

RK0

[A]

HGM

2.0 x
10ŧ

ITA

RK2

34°

19

TY

RT0

35°

12

137°

RT1

34

HGM ICU

NSZ1

KST

2.0 x
10ŧ
[rad]

2.0 x
10ŧ

2.0 x
10ŧ
[rad]

1, 2, 3, 4 (N290E, N20E, N65E, 1(

2.0 x
10ŧ

2.0 x
10ŧ
[rad]


(
5, 6, , 8 (N65E, N155E, N200E, 1(

50 km

2.0 x
10ŧ
[rad]

2.0 x
10ŧ

2.0 x
10ŧ
[rad]

ANO1, ANO2
0.2
[m]

12

50
[/h]
0

2.0 x
10ŧ
50
[/h]
0

19

135

136

10/10

10 11 12 13 14 15 16 1 18 19 20 21 22 23 24 25 26 2 28 29 30 31
2015/10137
11/01[A] 2015/10/25 00:00 - 10/28 00:00

Date

第 21 図 2015/10/10 00:00-11/1 00:00 の紀伊半島〜東海地方における深部低周波微動の時空間分布図．エンベロープ相関法により決定．
Fig. 21 Epicentral and space-time distributions of deep low
frequency tremors in Kii Peninsula and the Tokai district
from 0:00 October 10 to 0:00 November 1, 2015.

[A] 2015/10/25-27
(a)

GNOH

NSZ

ANO
URSH

34.5°

MGWH

KRTH

34.0°

2.0

WATH

1.5

ICU

residual

tremor (AIST)
LFE (JMA)

HGM
KHOH
NKTH
KST

33.5°

1.0

Slip: 27 mm, Mw: 5.7(b1)

135.5°

136.0°

136.5°

GNOH

tremor (AIST)
LFE (JMA)

137.5°

(b2)

NSZ

2B

ANO1

ANO1,2

Obs.

(b3)
Calc. x10ŧKST

0

ANO1ANO2

URSH

2A

34.5°

137.0°

KRTH

MGWH

HGM

WATH

ŧ

ICU
34.0°

ICU

ANO

HGM
KHOH
NKTH

Obs.
Calc.

KST

33.5°

Tilt
1x10ŧ rad

ŧ

NSZ

Lat. 34.33° Lon. 136.43° Dep. 26 km Len. 37 km Wid. 37 km
Strike 208° Dip 18° Rake 83° Slip 10 mm Mw 5.8
135.5°

136.0°

136.5°

137.0°

137.5°

ŧ

2x10ŧ strain

50
[/h]
0

MYM
ANO2

N, E

ITA1

N, E

N, E

.57+ +LŧQHW

N, E

0*:+ +LŧQHW

Obs.
Calc.

Expansion
Contraction

第 23 図 2015/10/25-27 の歪・傾斜・地下水変化（第 22 図赤色[A]）を説明する断層モデル．地図中の太丸
印は解析に使用した観測点（黒：産総研，緑：気象庁，青：防災科研 Hi-net）を示す．小さい青色丸印は産総研による深部低周波微動，小さい緑色丸印は気象庁一元化カタログによる低周波イ
ベントの震央位置を示す．(a)プレート境界面に沿って分布させた 20x20 km（大きさ固定）の矩形断層面に対して，おのおのの断層面で最適なすべり量を選んだときの残差の総和をその断層面の位置に示した分布図．赤色矩形が残差の総和が最小となる断層面．(b1)(a)の断層面付近でのグリッドサーチにより推定された断層面（赤色矩形）と断層パラメータ．灰色矩形は最近周辺で発生した短期的 SSE の推定断層面（1: 2015/5/7-9 (Mw 5.7), 2A: 2015/7/18-19 (Mw 5.5), 2B:
2015/7/20-25 (Mw 5.6)）．(b2)主歪の観測値と(b1)に示した断層モデルから求めた計算値との比較．
Fig. 23 The fault model for the strain, tilt and groundwater changes from October 25 to October 27, 2015. The
black, green and blue bold circles show the observation sites of AIST, JMA and NIED Hi-net, respectively.
The blue and green small circles show hypocenters of deep low frequency tremors estimated by AIST and
low frequency earthquakes estimated by JMA, respectively. (a) The distribution of the sum of residuals
between observed and calculated values that are caused by the uniform slip on a rectangular fault on theplate interface. The length and width of the rectangular fault are fixed as 20 km and 20 km, and only theslip amount that minimizes the sum of residuals is estimated. The red rectangle shows the fault patchwith minimum residual. (b1) The red rectangle shows the estimated fault model. The gray rectangles
show the fault models of the recent events (1: 2015/5/7-9 (Mw 5.7), 2A: 2015/7/18-19 (Mw 5.5), 2B:
2015/7/20-25 (Mw 5.6)). The observed and calculated tilt changes are also shown on the map by the redand blue arrows. (b2) The observed and calculated principal strain changes.

－ 264 －

50[mm/h]
050[mm/h]
0N, E

1.7+ +LŧQHW

.+2+ +LŧQHW

N, E

856+ +LŧQHW

N, E:$7+ +LŧQHW

N, E*12+ +LŧQHW

N, E

RK0, RK1, RK2, RK3
10 11 12 13 14 15 16 1 18 19 20 21 22 23 24 25 26 2 28 29 30 31
2015/1050[/h]0[A] 2015/10/25 00:00 - 10/28 00:00

第 22 図 紀伊半島・東海地方における歪・傾斜・地下水観測結果(2015/10/10 00:00 - 2015/11/1 00:00 (JST))
Fig. 22 Observed strain, tilt and groundwater changes at theobservation sites and number of deep low frequency tremorsin Kii Peninsula and the Tokai district from 0:00 October 10
to 0:00 November 1, 2015.

[A]

26

RK0, RK1, RK2, RK3[A]

HGM