10－３

平成28年（2016年）熊本地震
The 2016 Kumamoto Earthquake
気象庁
Japan Meteorological Agency【地震活動の概要】
2016年4月14日21時26分からの「平成28年（2016年）熊本地震」の地震活動は，2016年5月以降，全体として減衰傾向がみられる．熊本県熊本地方及び阿蘇地方においては，減衰しつつも活動は継続している．大分県中部の活動は低下した（第1図～第2図）．布田川断層帯・日奈久断層帯周辺や大分県中部における詳細な地震活動経過をみても，特段の活動域の広がりや活発化はみられない（第3図～第6図）．熊本地方及び阿蘇地方を小領域に分割すると，地震活動の減衰の傾向は領域によって若干の差が認められる（第7図）．
2016年5月以降で最大規模の地震は，8月31日19時46分に熊本県熊本地方（布田川断層帯（宇土区間）周辺）の深さ13kmで発生したM5.2の地震（最大震度5弱）であり，翌日9月1日には，この地震とほぼ同じ場所で，M4.8の地震（最大震度4）が発生した（第8図）．これらの地震の発震機構は，南北方向に張力軸を持つ正断層型であり，4月14日以降に周辺でみられていた発震機構と同様であった（第9図）．これらの地震発生後1週間程度の地震活動について，波形相関を用いたDouble-Difference法による震源再決定を行ったところ，北方向に深くなる震源分布がみられた（第10図）．また，国土地理院による熊本地震の余効変動モデル1)を用いて8月31日（M5.2）の地震に対する影響（ΔCFF）を求めたところ，この地震は，潮汐応答を大きく超えない程度の誘発域に位置していることが分かった（第11図）．
【地震活動の推移】
熊本県熊本地方・阿蘇地方における地震活動の推移から，今後の一定規模以上の地震の発生頻度を推定するため，現状の活動レベルを，活動初期（4月24日～5月13日）の活発な時期と比較した．
その結果，5月14日から約3週間の活動レベルはM6以上の地震回数でみると0.041回/日となり，活動初期の0.17回/日に比べ，およそ4分の1に減少した（第12図）．また，6月14日から約3週間の活動レ
ベルはM5以上の地震回数でみると0.31回/日となり，活動初期の1.77回/日に比べ，およそ6分の1に減少した（第13図）．熊本地震発生前後での活動レベル（1日あたりのM5以上の発生回数）を比較す
ると，熊本地震発生直後の活動初期では約150～200倍であったが，その後，約30～50倍にまで低下した（第14～15図）．
【過去の地震活動の事例】
熊本地震の活動の推移予測のために参照した，日本国内で発生した内陸地震での活動経過の過去事例を示す．2004年（平成16年）新潟県中越地震や2011年4月11日の福島県浜通りの地震では，最大規模の地震から1～2ヶ月後に最大震度5弱を観測する地震が発生した（第16図）．九州地方ではM6クラスの地震が続発した事例が複数あり（第17図），1975年の熊本県阿蘇地方から大分県西部にかけての地震活動や1997年の鹿児島県薩摩地方の地震活動では2～3ヶ月の間をおいて同程度の規模の－ 431 －

地震が発生した（第18図）．
【余震確率】
4月14日のM6.5の地震について，翌15日に余震確率の情報として，「M6（所によって震度6弱程
度）以上の余震発生確率が3日間で20％」と発表したが，16日にM7.3の地震が発生し，さらに地震活動が熊本県から大分県にかけての広域に広がったことから，通常の本震－余震型ではなくなったと判断し，余震確率の情報の発表を取りやめた．ここでは，その後，地震調査研究推進本部で確認された「大森・宇津公式を拡張した本震の余震活動とその二次余震系列の活動を重ね合わせたモデル2) 3)」等に基づく余震確率算出手法4)によって，7月11日以降の3日間の余震確率を計算した結果を示している（第19図～第21図）．
参考文献
1)

国土地理院：第212回地震予知連絡会国土地理院資料，2016．

2)

Ogata, Y.：Estimation of the Parameters in the Modified Omori Formula for Aftershock Frequencies by

the Maximum Likelihood Procedure, J. Phys. Earth, 31, 115-124, 1983.
3)

宇津徳治：地震活動総説，876, 1999．

4)

地震調査研究推進本部：大地震後の地震活動の見通しに関する情報のあり方，2016．

－ 432 －

「平成 28 年（2016 年）熊本地震」
今期間（2016 年５月～10 月）、「平成 28 年（2016 年）熊本地震」の活動で最大規模の地震は、
８月 31 日 19 時 46 分に熊本県熊本地方の深さ 13km で発生した M5.2 の地震（最大震度５弱）であった。発震機構は、南北方向に張力軸を持つ正断層型であった。この地震により軽傷２人などの被害が発生した（９月７日 12 時 00 分現在、総務省消防庁による）。また、今期間における阿蘇地
方での最大規模の地震は５月５日 10 時 40 分に発生した M4.9 の地震（最大震度４）であった。
熊本県熊本地方（領域ａ３）及び阿蘇地方（ａ２）における「平成 28 年（2016 年）熊本地震」の一連の地震活動は、全体として引き続き減衰しつつも、活動は継続している。大分県中部（ａ１）の活動は低下した。
今期間に震度１以上を観測した地震は 1,099 回（最大震度５弱：２回、最大震度４：17 回、最大震度３：79 回、最大震度２：258 回、最大震度１：743 回）発生した。
今回の一連の地震活動により、死者 140 人、負傷者 2,586 人、住家全壊 8,305 棟などの被害が発生した（11 月２日 10 時 00 分現在、総務省消防庁による）。

震央分布図（2016 年４月 14 日 21 時～10 月 31 日、深さ０～20km、M≧2.0）
2016 年５月以降の地震を濃く表示

ａ

Ａ

CMT
布田川断層帯別府－万年山断層帯

ａ１
ａ２
CMT

ａ３

Ｂ

日奈久断層帯CMT

CMT

図中の細線は地震調査研究推進本部の長期評価による活断層を示す
M6.0 以上の地震と領域 a１、a２における４月 14 日以降の最大規模の地震に吹き出し（黒色）をつけている。また、a２、a３における５月以降の最大規模の地震に吹き出し（赤色）をつけている。
※１ M7.3 の地震の発生直後に発生したものであり、M の値は参考値。

第1図 地震活動
Fig.1 Seismic activity.

－ 433 －

領域ａ内の時空間分布図（Ａ－Ｂ投影）（2016 年４月 14 日 21 時～10 月 31 日）

領域 a1,a2,a3 内のＭ－Ｔ図及び回数積算図（2016 年４月 14 日 21 時～10 月 31 日）

Ａ

ａ１

ａ２
Ｂ
領域ａ内のＭ－Ｔ図及び回数積算図ａ３「平成 28 年（2016 年）熊本地震」の地震回数
最大震度別回数期間
4/14-4/30
5/1-5/31
6/1-6/30
7/1-7/31
8/1-8/31
9/1-9/30
10/1-10/31
合計

１

２

３

４

５弱 ５強 ６弱 ６強

７

1722
344
147
85
77
49
41
2465

859
134
51
19
28
16
10
1117

323
43
14
8
3
7
4
402

98
8
4
1
2
2

10

115

12

5

3

2

2

5

3

2

2

1
1

合計

累計

3024
529
217
113
111
74
55
4123

3024
3553
3770
3883
3994
4068
4123

1223

震度１以上の日別地震回数グラフ（2016 年４月 14 日 21 時～10 月 31 日）

第2図 つづき
Fig.2 Continued.

－ 434 －

布田川断層帯周辺の活動
2016年4月14日～10月31日、M≧2.0、深さ0～20km
緑：4月～7月、赤：8月～10月
M5.0以上の地震、または最大震度５弱以上の地震を濃く表示。
【領域a】布田川断層帯全域【領域ｂ】宇土区間・宇土半島北岸区間

震央分布図【領域c】布田川区間

領域a内のMT図・回数積算図
B

領域a
C
A

領域c

領域b

領域a内の時空間分布図（B-A投影）
B

A

領域a内の断面図（A-B投影）

B

A

宇土区間・宇土半島北岸区間
領域b内の断面図（B-C投影）
B
C

B

領域b内の時空間分布図（B-C投影）C
領域b内のMT図・回数積算図
領域b内の深さの時間変化

第3図 地震活動（布田川断層帯周辺）
Fig.3 Seismic activity (around Futagawa fault).

－ 435 －

布田川断層帯周辺の活動
2016年4月14日～10月31日、M≧2.0、深さ0～20km
緑：4月～7月、赤：8月～10月
M5.0以上の地震、または最大震度５弱以上の地震を濃く表示。
【領域a】布田川断層帯全域【領域ｂ】宇土区間・宇土半島北岸区間【領域c】布田川区間

布田川区間
領域c内の断面図（ B-C投影）
B
C

B

領域c内の時空間分布図（B-C投影）

C
領域c内のMT図・回数積算図
領域c内の深さの時間変化

第4図 つづき
Fig.4 Continued.

－ 436 －

日奈久断層帯周辺の活動
2016年4月14日～10月31日、M≧2.0、深さ0～20km
緑：4月～7月、赤：8月～10月
M5.0以上の地震、または最大震度５弱以上の地震を濃く表示。
【領域a】日奈久断層帯全域【領域ｂ】高野－白旗区間

震央分布図【領域c】日奈久区間

領域a内のMT図・回数積算図

領域b
B

C

領域c

領域a
A

領域a内の時空間分布図（B-A投影）
B

領域a内の断面図（A-B投影）
A
B

A

高野－白旗区間

日奈久区間

領域b内の断面図（B-C投影）
B
C

領域c内の断面図（ B-C投影）
B
C

領域b内のMT図・回数積算図

領域c内のMT図・回数積算図

第5図 地震活動（日奈久断層帯周辺）
Fig.5 Seismic activity (around Hinagu fault).

－ 437 －

大分県中部から伊予灘にかけての活動

領域ａ内の時空間分布図（A-B投影）
A

震央分布図（2016年4月14日21時～10月31日24時、マグニチュード全て、深さ20km以浅）
７月24日までに発生した地震を赤く、それ以降を青く表示。

Aａ
別府－万年山断層帯

B

B（※）熊本地方で発生したM7.3の地震の発生直後に発生したものであり、Mの値は参考値。
4月14日～7月24日のM2.0未満の震源には未精査の自動処理結果が含まれており、今後の精査により変更される場合がある。

領域ａ内のM-T図
及び回数積算図（Ｍ≧1.7）

震央分布図（1997年10月～2016年10月、Ｍ≧1.7、深さ20km以浅）
2016年4月14日21時より前を灰色、それ以降は上図と同様に表示。

Ｃｂ

領域ｂ内の時空間分布図（Ｃ-Ｄ投影）
Ｃ

Ｄ

参考

領域ｂ内のM-T図
及び回数積算図Ｄ（※）熊本地方で発生したM7.3の地震の発生直後に発生したものであり、Mの値は参考値。

第6図 地震活動（大分県中部～伊予灘）Fig.6 Seismic activity (from the central part of Oita Prefecture to Iyo-nada).

－ 438 －

ａ

ｂ

ｃ

ｄ

ｅ

ｆ

ｇ

ｈ

「平成28年（2016年）熊本地震」（領域別の活動）
震央分布図（2016年４月14日～2016年10月31日24時00分、深さ０～20km、M2.0以上）

ｅ
ｄ
ｃ

ｂ
ａ

ミャンマー

ｆ
ｇ
ｈ
丸の大きさはマグニチュードの大きさを表す。‐‐‐‐‐‐は地震調査研究推進本部の長期評価による活断層の位置を表す

上図の震央分布図内の地震活動経過図

横軸は時間、縦軸はマグニチュード、縦棒のついた丸は地震発生時刻とマグニチュードの大きさを表す。

第7図 地震活動（領域別）
Fig.7 Seismic activity (each subdivided region).

－ 439 －

布田川断層帯・宇土区間周辺の活動
2016年7月25日～9月1日08時、M≧1.0、深さ0～20km
8月15日～8月30日の地震を緑で表示
8月31日以降の地震を赤で表示

布田川断層帯宇土半島北岸区間

震央分布図

青色矩形内の時空間分布図
（A-B投影）布田川断層帯宇土区間

A

青色矩形内のMT図・回数積算図（上段）
深さの時間変化（下段）
MT図・回数積算図

B

今回の地震
深さの時間変化

青色矩形内の断面図（A-B投影）

今回の地震

7月

第8図 2016年8月31日 熊本県熊本地方の地震（M5.2）
Fig.8 The earthquake (M5.2) in Kumamoto region of Kumamoto Prefecture on August 31, 2016.

－ 440 －

8月

布田川断層帯・日奈久断層帯付近の発震機構解
震央分布図（2016年4月14日21時～9月1日、M≧3.0、20km以浅）ｂ
初動発震機構解の断層の型別（分類はFrohlich（1992）による）に以下の色分けで表示。
赤：逆断層型
青：正断層型
緑：横ずれ断層型
灰：その他

ａ
ｃ

なお、図の期間中には未処理のデータがある。

ｄ

領域別の三角ダイアグラム・Ｔ軸方位
（領域ｂ）

（領域ａ）

合計
横ずれ断層型
正断層型
逆断層型
その他

37個
3個（ 8.1％）
28個（75.7％）0個（ 0.0％）6個（16.2％）

合計 26個横ずれ断層型 2個（ 7.7％）
正断層型
9個（34.6％）
逆断層型
1個（ 3.8％）
その他
14個（53.8％）

8月31日M5.2の地震を強調して表示

（領域ｄ）

（領域ｃ）

合計
横ずれ断層型
正断層型
逆断層型
その他

36個
10個（27.8％）8個（22.2％）
1個（ 2.8％）
17個（47.2％）

第9図 つづき（発震機構解の分布）
Fig.9 Continued (distribution of focal mechanism).

－ 441 －

合計
横ずれ断層型
正断層型
逆断層型
その他

56個
24個（42.9％）
7個（12.5％）
4個（ 7.1％）
21個（37.5％）

平成28年８月31日

熊本県熊本地方の地震（宇土区間の活動）波形相関ＤＤ法による再決定震源

一元化震源（ルーチン）
震央分布図（2016年８月31日０時～９月７日12時、
深さ０～20km、M0.5以上）

ａ

震央分布図（2016年８月31日０時～９月７日12時、
深さ０～20km、M0.5以上）

震源フラグ：
K、k、Aａ

9/7 Ｍ3.9の
自動初動解

M5.2の地震より前の震源を
青、同地震以降の震源を赤で表示。
赤・茶線は地震調査研究推進本部の長期評価による活断層の位置を表す。

上図の領域ａ内のＭＴ図

上図の領域ａ内の断面図（南北投影）

上図の領域ａ内の断面図（南北投影）

宇土区間
地表トレー
スの分布
範囲

第10図 つづき（波形相関を用いたDouble-Difference法による震源）
Fig.10 Continued (hypocenter distribution by the Double-Difference method based on waveform
correlation).

－ 442 －

震源フラグ：
K、A熊本地震余効変動モデルの2016年8月31日熊本地方の地震（M5.2）に対するΔCFF
熊本地震余効変動モデルは国土地理院の暫定モデル（第212回地震予知連絡会資料）による。
2016年8月31日熊本地方の地震のメカニズムは気象庁初動解の２面（パラメータ下記）を使用。

断層面解１（strike277，dip31，rake‐77）に対する計算断層面解２（strike81，dip60，rake‐98）に対する計算深さ１３km（M5.2の一元化
震源の深さ）で計算

1E‐7

1E‐7

震央

1E‐6

震央

1E‐6

1E‐6
1E‐6

今回の震源付近のΔCFFは5E‐7 （剛性率30GPaとすると0.15bar）で促進域

細線は新編日本の活断層による断層トレース

今回の震源付近のΔCFFは9E‐7 （剛性率30GPaとすると0.27bar）で促進域

第11図 つづき（静的応力変化（ΔCFF））
Fig.11 Continued (changes in static stress (ΔCFF)).

－ 443 －

熊本地方・阿蘇地方の地震活動状況の比較（2016年4月24日～6月3日）
震央分布図（Ｍ≧2.0、深さ20km以浅） ・ＭＴ図（領域ａ内）① 2016年4月24日～5月13日② 2016年5月14日～6月3日

ａ

ａ

4/25

4/30

5/5

5/15

5/10

5/20

5/25

5/30

地震発生間隔の累積曲線（領域ａ内、Ｍ≧3.0）

規模別度数分布（領域ａ内、Ｍ≧2.0）

回

回
1000

● Ｎ（Ｍ以上積算回数）○ ｎ（Ｍ別回数）

①

①

129回
最大間隔0.77日

100

Ｍ≧3.0の地震発生状況から推定されるM6.0以上の回数＝0.17回

10

1

回
1000
0.1

日

②

回

100
0.01

0

1

2

3

4

5

6

Ｍ

43回

② 最大間隔3.09日

10

1

Ｍ≧3.0の地震発生状況から推定されるM6.0以上の回数＝0.041回

0.1

0.01

0

1

2

3

4

5

第12図 地震活動の推移
Fig.12 Change of seismic activeness.

－ 444 －

6Ｍ日

熊本地方・阿蘇地方の地震活動の比較（2016年4月24日～5月13日と6月14日～7月3日）
震央分布図（Ｍ≧2.0、深さ20km以浅） ・ＭＴ図（領域ａ内）① 2016年4月24日～5月13日② 2016年6月14日～7月3日

ａａ

規模別度数分布（領域ａ内、Ｍ≧2.0）
回Ｍ≧3.0の地震発生状況から推定される
M5.0以上の回数＝1.8回

回

Ｍ≧3.0の地震発生状況から推定される
M5.0以上の回数＝0.31回

1000

1000

ｂ値0.90（σ＝0.20）

ｂ値0.93（σ＝0.08）

100

100

6/14～7/3の発生状況から推定されるM5.0以上の回数は4/24～
5/13に比べると
約６分の１

10

10

1

1

● Ｎ（Ｍ以上積算回数）○ ｎ（Ｍ別回数）

0.1

0

1

2

3

4

5

6

Ｍ

0.1

0

第13図 つづきFig.13 Continued.

－ 445 －

1

2

3

4

5

6

Ｍ

熊本地方・阿蘇地方の地震活動の比較（熊本地震以前と2016年6月14日～7月3日）
震央分布図（Ｍ≧4.0、深さ30km以浅） ・ＭＴ図（領域ａ内）③ 1923年1月1日～2016年4月13日

震央分布図（Ｍ≧2.0、深さ20km以浅） ・ＭＴ図（領域ａ内）② 2016年6月14日～7月3日

ａａ

規模別度数分布（領域ａ内、Ｍ≧2.0）

規模別度数分布（領域ａ内、Ｍ≧4.0）
回

約93年間のＭ≧4.0の地震発生状況から推定される
M5.0以上の回数＝約15回（20日間では0.0087回）

20日間のＭ≧3.0の地震発生状況から推定されるM5.0以上の回数＝0.31回

回
1000

1000

ｂ値0.99（σ＝0.08）

100

20日間で推定されるM5.0以上の回数の比（熊本地震発生前を1とする）
2016/4/24～5/13
約197倍

10

ｂ値0.90（σ＝0.20）

100

10

2016/5/14～6/3
約49倍
1

2016/6/14～7/3
約35倍

0.1

0

1

2

3

4

5

6

Ｍ

● Ｎ（Ｍ以上積算回数）○ ｎ（Ｍ別回数）

第14図 つづきFig.14 Continued.

－ 446 －

1

0.1

0

1

2

3

4

5

6

Ｍ

熊本地方・阿蘇地方の地震活動の比較（熊本地震以前と2016年6月14日～7月3日）
震央分布図（Ｍ≧5.0、深さ30km以浅） ・ＭＴ図（領域ａ内）③ 1923年1月1日～2016年4月13日

震央分布図（Ｍ≧2.0、深さ20km以浅） ・ＭＴ図（領域ａ内）② 2016年6月14日～7月3日

ａａ

規模別度数分布（領域ａ内、Ｍ≧2.0）
約93年間のＭ≧5.0の地震発生回数
20回（20日間では0.012回）

回

20日間のＭ≧3.0の地震発生状況から推定されるM5.0以上の回数＝0.31回

1000

ｂ値0.90（σ＝0.20）

20日間で推定されるM5.0以上の回数の比（熊本地震発生前を１とする）
2016/4/24～5/13 ：約151倍

2016/5/14～6/3

：

約37倍

2016/6/14～7/3

：

約26倍

● Ｎ（Ｍ以上積算回数）○ ｎ（Ｍ別回数）

100

10

1

0.1

0

第15図 つづきFig.15 Continued.

－ 447 －

1

2

3

4

5

6

Ｍ

本震の１～２ヶ月後に最大震度５弱を観測する余震が発生した事例●2004年（平成16年）新潟県中越地震
上：時空間分布図（ＡＢ投影、左図矩形内）
下：ＭＴ図（左図矩形内）

震央分布図（2004年10月23日～2005年1月22日、M≧2.0、深さ25km
以浅）ＡＡ

長岡平野西縁
断層帯

Ｂ
１１月
２００４年

六日町断層帯・北部

１月
２００５年

2004年10月23日
M6.8

Ｂ
十日町断層帯

１２月

六日町断層帯・南部

2004年12月28日
M5.0

本震の66日後に東北東約10kmでM5.0の地震が発生し、最大震度５弱を観測した。
１１月
２００４年

１２月

１月
２００５年

●2011年4月11日の福島県浜通りの地震
震央分布図（2011年4月11日～7月10日、M≧2.0、深さ25km以浅）

上：時空間分布図（ＣＤ投影、左図矩形内）
下：ＭＴ図（左図矩形内）
Ｃ

Ｃ

Ｄ
４月

Ｄ

2011年4月11日
M7.0

５月

６月

７月

2011年5月6日
M5.2
2011年5月25日
M5.0

本震震央の北北東約20kmで、本震の25日後にM5.2の地震、本震の44日後にM5.0の地震が発生し、それぞれ最大震度５弱を観測した。

４月

上図の発震機構解は全てCMT解。活断層トレースは地震本部の長期評価による。

第16図 過去の地震活動における続発事例
Fig.16 Sequencial occurring earthquakes in the past seismic activity.

－ 448 －

５月

６月

７月

九州地方における内陸・浅い地震の震央分布（1889年以降）
②
③
①
④

⑤

④
①

③

②

⑤

（地震予知連：大竹名誉委員の指摘）①1889年の熊本の地震の際にはその後６年でM６クラスが４回発生②1916年にはM６クラスが10ヶ月の間をおいて発生③1975年にM６クラスの地震が３ヶ月の間隔で発生※④⑤はその他の規模のやや大きい活動
第17図 つづきFig.17 Continued.

－ 449 －

九州地方において近接地域で２～３ヶ月の間をおいて同程度の地震が発生した事例●1975年の熊本県阿蘇地方から大分県西部にかけての地震活動
上：時空間分布図（ＡＢ投影、左図矩形内）
下：ＭＴ図（左図矩形内）

震央分布図（1975年1月1日～5月31日、M≧3.0、深さ20km以浅）ＡＡ
別府－万年山
断層帯

Ｂ
Ｂ

1月

２月

３月

４月

1975年4月21日
M6.4

1975年1月23日
M6.1

熊本県阿蘇地方のM6.1の地震の88日後に北東へ約20km離れた大分県西部でＭ6.4が発生。熊本県阿蘇地方の地震では負傷者10人、住家全壊16棟、大分県西部の地震では負傷者22人、住家全壊58棟の被害がそれぞれ発生した（被害は「九州地域の活断層の長期評価」より）。

５月

1月

２月

３月

４月

５月

３月

４月

５月

６月

７月

４月

５月

６月

７月

●1997年の鹿児島県薩摩地方の地震活動
震央分布図（1997年3月1日～7月31日、M≧2.0、深さ20km以浅）

北

出水断層帯

南

東

西
右上：時空間分布図（南北投影、左図矩形内）
右中：時空間分布図（東西投影、左図矩形内）
右下：ＭＴ・回数積算図（左図矩形内）

３月

1997年3月26日
M6.6

1997年5月13日
M6.4

M6.6の地震の48日後に西南西約6kmでM6.4の地震が発生。
M6.6の地震では負傷者31人、住家全壊4棟など、M6.4の地震では負傷者43人、住家全壊4棟などの被害がそれぞれ発生した（被害は「九州地域の活断層の長期評価」より）。

上図の発震機構解は全てCMT解。活断層トレースは地震本部の長期評価による。 ３月

第18図 つづきFig.18 Continued.

－ 450 －

４月

５月

６月

７月

熊本県熊本地方と熊本県阿蘇地方を合わせた領域の余震発生確率
2016 年 4 月 14 日 21 時 26 分（M6.5 発生）～7 月 8 日 11 時 00 分までの M3.6 以上のデータをもとに、7 月 11 日 17 時 00 分での確率を算出。
M5.0 以上の地震が発生する３日間確率は 4％
M4.5 以上の地震が発生する３日間確率は 12％

大分県中部●：M 大→M 小へ
地震回数を積算

阿蘇地方

したもの○：M ごとの地
熊本地方

震回数

地震積算回数（回）

地震発生確率（％）

地震回数積算の予測

M4.5 以上の地震発生の３日間確率
12％

7/11
17 時

実際の活動（－）
・ETAS モデル（－）
・改良大森モデル（－）の比較

地震積算回数（回）

モンテカル
ロ法による
予測の幅

・改良大森の AIC は-1745.74 ETAS の AIC は-1664.43
・ETAS による M5.0 以上の地震が発生する３日間確率は 4％・ETAS による M4.5 以上の地震が発生する３日間確率は 15％→予測期間

第19図 余震確率
Fig.19 Probability of the aftershock occurring.

－ 451 －

熊本県熊本地方の余震発生確率
2016 年 4 月 14 日 21 時 26 分（M6.5 発生）～7 月 8 日 11 時 00 分までの M3.6 以上のデータをもとに、7 月 11 日 17 時 00 分での確率を算出。
M5.0 以上の地震が発生する３日間確率は 4％
M4.5 以上の地震が発生する３日間確率は 10％

大分県中部●：M 大→M 小へ
地震回数を積算

阿蘇地方

したもの○：M ごとの地

熊本地方

震回数

地震積算回数（回）

地震発生確率（％）

地震回数積算の予測

M4.5 以上の地震発生の３日間確率
10％

7/11
17 時

実際の活動（－）
・ETAS モデル（－）
・改良大森モデル（－）の比較

地震積算回数（回）・改良大森の AIC は-1197.61 ETAS の AIC は-1102.95
・ETAS による M5.0 以上の地震が発生する３日間確率は 4％・ETAS による M4.5 以上の地震が発生する３日間確率は 12％

モンテカル
ロ法による
予測の幅

→予測期間

第20図 つづきFig.20 Continued.

－ 452 －

大分県中部の余震発生確率
2016 年 4 月 16 日 01 時 25 分（M7.3 発生）～7 月 8 日 11 時 00 分までの M3.0 以上のデータをもとに、7 月 11 日 17 時 00 分での確率を算出。
M4.0 以上の地震が発生する３日間確率は 0.1％

大分県中部●：M 大→M 小へ
地震回数を積算

阿蘇地方

したもの○：M ごとの地

熊本地方

震回数

地震積算回数（回）

地震発生確率（％）

地震回数積算の予測

M4.0 以上の地震発生の３日間確率
0.1％
7/11
17 時

実際の活動（－）
・ETAS モデル（－）
・改良大森モデル（－）の比較

地震積算回数（回）・改良大森の AIC は-379.33 ETAS の AIC は-292.60
・ETAS による M4.0 以上の地震が発生する３日間確率は 0％

モンテカル
ロ法による
予測の幅

→予測期間

第21図 つづきFig.21 Continued.

－ 453 －