地震の予測マップと発震日予測

地震発生ヶ所をポイント予測し、度数分布から発震日を確率予測する

2020-02-05 地震の予測マップと時系列予測グラフ 6日の地震列島は, 茨城沖でM5.6, 宮城沖でM3.6! 解説:南関東_監視域の度数分布_発震予測確率を改良しました! 

常に最新記事なら [こちら最新!] をアクセスし、ブックマーク!

 防災科研さんから午前0時に2日前の詳細データが公開され、もって1年分のデータ解析を行なっています、題名先頭にある日付が解析データ1年分の最終日です、記事は毎日、午前1時〜3時頃にアップされます

 

 

[こよみの計算 - 国立天文台暦計算室] さんより:

⚫ 2月5日は月齢10.7日 東進Day

⚫ 2月6日は月齢11.7日 東進Day

⚫ 2月7日は月齢12.7日 東進Day

東進Dayとは新月から満月前日まで、西進Dayとは満月から新月前日までの日々

 

各領域の時系列データは排他的にして重複を排除しており、西域が南海トラフ監視領域を100%包含するので最も優先度が高く、順に以下の如くです

 西域_時系列:西域そのものであり、フィリピン海プレートの影響直接

 中域_時系列:中域から西域を除き、日本海溝から太平洋プレートの影響直接

 東域_時系列:東域から中域を除き、千島海溝から太平洋プレートの影響直接

= 最新地震情報2月6日(M3.0以上かつ震度1以上)です =

Yahooさん [4] より掲載(元データは [気象庁] さん)、マップ上震源位置

 6日14時30分、宮城沖でM3.6、深さ50km、震度1 中域_時系列

 6日20時20分、茨城沖M5.6、深さ40km、震度2 中域_時系列

 

 6日、茨城沖でM5.6が発生しました、茨城沖とその周辺は警戒です

 

 M6.0以上地震注意情報です

 この4年間、M6.0以上は各域合計で31件発生しており、平均発震間隔は41.5日、最後の発震 2019/08/29 M6.1 青森東方沖 から2月7日現在で162日が経過しています、既に平均の4倍近くになっており、全域において警戒が必要です

 この4年間、西域_時系列では木星の衝合期間中にM6.0以上が11件中10件発生しており(発生確率90.9%)、2月7日現在は合の期間中で、それは2月15日まで続きますので注意が必要です

 

 

解説: 南関東_監視域の度数分布_発震予測確率を改良しました!

 ARIMA時系列予測に代わる方式として、度数分布_発震確率予測を構築し南関東_監視域に適用(暫定版)した訳ですが、これをポリッシュし、一気に東/中/西_時系列の予測を置き換える予定です、そこで、まずは改良した南関東_監視域をご説明致したく

 

 宇宙動画と画像と話題の名ハンターであると共に本予測マップ&グラフの名コメンテータでもある「まさき りお(id:ballooon)」さんから、先日の南関東_監視域の度数分布_発震予測確率(暫定版)についてブクマコメを頂戴しており:

南関東_監視域は、、、30日までは50%、45日からは75%、60日以降100%の確率とする、ですね?難しい・・

というご指摘で、実は私も少し難しい(複雑)かな?と思っていた次第でして、ここは正規分布における標準偏差に立脚した確率%表示に変更致しました

 

 即ち、Wiki [標準偏差 - Wikipedia] より正規分布標準偏差との関係を示すと、

https://upload.wikimedia.org/wikipedia/commons/thumb/8/8c/Standard_deviation_diagram.svg/1024px-Standard_deviation_diagram.svg.png

平均 0, 標準偏差 σ 正規分布確率密度関数。この分布に従う確率変数が 0 ± σ の間に値をとる確率はおよそ 68% であることが読み取れる。

になる訳ですが、ここで右肩下り分布では上記確率密度分布を中央から左右に折りたたんだ密度分布を適用し、

➡ 標準偏差σの68%面積 これが基本!

➡ 68%の半分の34%面積

➡ 80%面積

➡ 100%面積

を含むbinsの右端の面積%4点を、そこに至る発震確率と致しました

 

 こうして南関東_監視域の度数分布_発震予測確率に適用しますと:

f:id:yoshihide-sugiura:20200207024401p:plain

➡ 標準偏差σの68%面積  そのまま68%で出現

➡ 半分の34%面積 ➡ 46%となって出現

➡ 80%面積   ➡ 89%となって出現

➡ 100%面積 ➡ これは100%となります

驚くべき事に、68%がbinsの右エッジに出ました、このbinsの幅20は私が指定している訳ではなく、Rが自動で設定してくるのです、明らかに、Rは正規分布標準偏差σを意識してヒストグラムを出しています! 面白いです、偶然でこのようになる訳がありません

尚、グリーンのラインは発震間隔(インターバル)の平均値(日数)で、参考までに出しました、が、発震確率%はすべて面積から出しています

現在、南関東_監視域のM5.0以上の発震確率は46%と予測される訳ですが、降水確率46%で雨は降る、という感覚でして、まぁ受け入れられるレベルではないでしょうか?

 

 地震の発生はポアソン分布に従うというのが一般に言われているのですが、南関東のように狭い領域に着眼した場合のポアソン分布は、右肩下り分布で正規分布のY軸対称に漸近する、と言えそうです

 

 これにて問題が無いようでしたら、ARIMA予測をすべて置き換えます、2〜3日位かかるか?と思っています

尚、時系列予測を度数分布_発震予測確率に置き換える為、月距離と1年26期は必要なくなりますので月距離の説明は削除しました、しかし、東進Dayと西進Dayは予測マップで使っているので、月距離地震グラフは満月新月地震グラフに名称変更ました、又、すべてのマップ&グラフは今後毎日更新となります

 

 

 

 

そして被災地は今... [happy-ok3の日記] 地震・豪雨・台風と、被災地の現状をレポートするhappy-ok3 さんの考えさせられるブログです、関心を持ち続けて欲しい

 

= 以下、2020-02-05 迄データ1年分による本日の地震の予測マップ&時系列予測グラフ =

マークは東進圧力、マークは西進圧力を示す圧力方向の解析表示です

救急マークがM5.0以上の地震予測する注意ポイントで、6kmマップにあります

[防災科学技術研究所 Hi-net 高感度地震観測網]、[気象庁|震源データ] を参照しています

= 地震の予測マップピッチ36kmマップとM6.0以上地震分析です =

東進西進圧力表示・ピッチ36km予測マップです

f:id:yoshihide-sugiura:20200207023908p:plain

ピンクの小さな★マークは、南海トラフ巨大地震発生ヶ所で、西から、1854安政南海M8.4、1946昭和南海M8.4、1707宝永M8.6、1944昭和東南海M8.2、1854安政東海M8.4

 

M6.0以上地震発生の分析:

 過去直近4年間の東中西_時系列におけるM6.0以上深さ130km未満の発震状況分析ですが、
東域_時系列 -> 4 回、中域_時系列 -> 16 回、西域_時系列 -> 11 回
であり全体で、合計31回、発震間隔は、平均で41.5日、最短で1日、最長で267日、です
最後の発震 2019/08/29 Mag_6.1, Dph_21km, 青森東方沖 から、2020/02/05現在、160日が経過しています

 

 以下のグラフは南海トラフ_監視域について、M4.8以上・500km未満・各日365日の総和を取ったグラフです、総和を取ると周期性が現れますが、周期は総和期間を変えると変化しますので、これはあくまでも参考グラフです、そこに上記31回のM6.0以上発震を東/中/西でE/M/Wマークしてみると、こうなりますというグラフです(E/M/Wマークは監視域内の発震では無い事に御注意下さい!)

f:id:yoshihide-sugiura:20200207023938p:plain

参考情報ですが、カウント4でM6.0以上の発震数は2個と少なく、全体比率は 2/31 = 6.5%
それは、以下の中域_時系列における2件で、日本海溝における発震:
・ 発生日:2019/08/04 中域_時系列、M6.4、福島沖
・ 発生日:2019/08/29 中域_時系列、M6.1、青森東方沖
です

 

 木星の衝合日に±48日を加えた97日間は、西域_時系列でM6.0以上の地震が多発します(強い相関があります)
過去4年間に発生した、M6.0以上地震数は東中西_各時系列総和で 31 回
東域と中域_時系列で、±48日以内のM6.0以上発震数は10、地震数は20、確率は、10/20 = 50.0%
西域_時系列における、±48日以内のM6.0以上発震数は10、地震数は11、確率は、10/11 = 90.9%
現在は 合の期間中で、それは2020/02/15まで続きます

以下、相関グラフで、木星の衝合については本文最後に説明があります

f:id:yoshihide-sugiura:20200207024020p:plain

±48日以内に存在した西域_時系列のM6.0以上地震は:
1 : 2016/04/01 West-> [ 33.32, 136.37, Dph_28km, Mag_6.5, : 南海トラフ_紀伊半島南東沖 ]
2 : 2016/04/14 West-> [ 32.73, 130.80, Dph_11km, Mag_6.5, : 熊本地震 ]
3 : 2016/04/15 West-> [ 32.70, 130.77, Dph_6km, Mag_6.4, : 熊本地震 ]
4 : 2016/04/16 West-> [ 32.75, 130.75, Dph_12km, Mag_7.3, : 熊本地震 ]
5 : 2016/10/21 West-> [ 35.37, 133.85, Dph_10km, Mag_6.6, : 鳥取中部地震 ]
6 : 2017/11/16 West-> [ 32.35, 140.75, Dph_46km, Mag_6.0, : 八丈島東方沖 ]
7 : 2018/04/09 West-> [ 35.18, 132.58, Dph_12km, Mag_6.1, : 島根広島県境 ]
8 : 2018/06/18 West-> [ 34.85, 135.62, Dph_13km, Mag_6.1, : 大阪北部地震 ]
9 : 2019/01/08 West-> [ 30.57, 131.17, Dph_30km, Mag_6.4, : 種子島付近 ]
10 : 2019/05/10 West-> [ 31.80, 131.98, Dph_25km, Mag_6.3, : 日向灘 ]

±48日以内から外れた西域_時系列のM6.0以上地震は:
1 : 2018/07/07 West-> [ 35.17, 140.60, Dph_56km, Mag_6.0, : 千葉東方沖 ]
です

 

 

= 地震の予測マップ・ピッチ6kmマップとポイント予測マップと月距離地震グラフと時系列地震予測グラフです =

東進西進圧力表示・ピッチ6km予測マップ東域です  凡例は36kmマップと同じ

f:id:yoshihide-sugiura:20200207024057p:plain

次がポイント予測・東域 救急マーク がM5.0以上の発生予測注意ポイントです

f:id:yoshihide-sugiura:20200207024112p:plain

根室沖M7.8〜8.5の確率80%程度 [北海道地方の地震活動の特徴 | 地震本部] 2020年1月1日算定基準日の30年間確率、以下同様

満月新月地震グラフ東域:この1年の東域_時系列データM5.0以上です

f:id:yoshihide-sugiura:20200207024125p:plain

26期予測グラフ(M5.0以上数は3リミット)東域2020-02-04 の予測結果で、

f:id:yoshihide-sugiura:20200206032024p:plain

26期予測グラフとは、凸部凹部の期に発生したM5.0以上(西域ではM4.8以上)地震数を最大3で制限し、各期に至る1年26期分の総和を26期分表示し、続く27/28期を予測するグラフです、以下同様

東域M5.0以上の発生予測 : 続く27/28期の予測は、 27th-> Hit_P, 28th-> Hit_P
27th-> 2020/02/05 to 2020/02/18, 28th-> 2020/02/19 to 2020/03/03

東域の過去1年Hit_P予測実績は93.3%、予測実績詳細は:

Hit_P : 14 期分が、Hit_P 状況にマッチしました
Hit_N : 11 期分が、Hit_N 状況にマッチしました
空振り : 1 期分が、空振り状況にマッチしました
見逃し : 発生せず
Hit_P 予測の実績は、14/(14+1) = 93.3%

Hit_P状況とはHit_P予測(M5.0以上が発生すると予測)して実際にM5.0以上が発生した状況、空振り状況とはHit_P予測してM5.0以上が発生しなかった状況、Hit_N状況とはHit_N予測(M5.0以上は発生しないと予測)して実際にM5.0以上は発生しなかった状況、見逃し状況とはHit_N予測して実際にはM5.0以上が発生した状況、以下同様

 

東進西進圧力表示・ピッチ6km予測マップ中域です 凡例は36kmマップと同じ

f:id:yoshihide-sugiura:20200207024147p:plain

次がポイント予測・中域 救急マーク がM5.0以上の発生予測注意ポイントです

f:id:yoshihide-sugiura:20200207024204p:plain

青森東方沖及び岩手沖北部M7.0〜7.5の確率90%程度以上 、宮城沖M7.0〜M7.5の確率90%程度、福島沖M7.0〜M7.5の確率50%程度、茨城沖M7.0〜M7.5の確率80%程度 [東北地方の地震活動の特徴 | 地震本部]

満月新月地震グラフ中域、この1年の域_時系列データM5.0以上です

f:id:yoshihide-sugiura:20200207024220p:plain

26期予測グラフ(M5.0以上数は3リミット)中域2020-02-04 の予測結果で、

f:id:yoshihide-sugiura:20200206032155p:plain

中域M5.0以上の発生予測 : 続く27/28期の予測は、 27th-> Hit_P, 28th-> Hit_P
27th-> 2020/02/05 to 2020/02/18, 28th-> 2020/02/19 to 2020/03/03

中域の過去1年Hit_P予測実績は75.0%、予測実績詳細は:

Hit_P : 12 期分が、Hit_P 状況にマッチしました
Hit_N : 10 期分が、Hit_N 状況にマッチしました
空振り : 4 期分が、空振り状況にマッチしました
見逃し : 発生せず
Hit_P 予測の実績は、12/(12+4) = 75.0%

 

東進西進圧力表示・ピッチ6km予測マップ西域です  凡例は36kmマップと同じ

f:id:yoshihide-sugiura:20200207024244p:plain

次がポイント予測・西域 救急マーク がM5.0以上の発生予測注意ポイントです

f:id:yoshihide-sugiura:20200207024303p:plain

茨城沖M7.0〜M7.5の確率80%程度、相模トラフ南関東直下地震M6.7〜7.3の確率70%程度、南海トラフM8〜M9クラスの確率70%〜80%、日向灘M7.1前後の確率70〜80% [関東地方の地震活動の特徴 | 地震本部] と[九州・沖縄地方の地震活動の特徴 | 地震本部]

ピンクの小さな★マークは、南海トラフ巨大地震発生ヶ所で、西から、1854安政南海M8.4、1946昭和南海M8.4、1707宝永M8.6、1944昭和東南海M8.2、1854安政東海M8.4

満月新月地震グラフ西域、この1年の西域_時系列データM4.8以上です

f:id:yoshihide-sugiura:20200207024318p:plain

26期予測グラフ(M4.8以上数は3リミット)西域2020-02-04 の予測結果で

f:id:yoshihide-sugiura:20200206032415p:plain

西域M5.0以上の発生予測 : 続く27/28期の予測は、 27th-> Hit_P, 28th-> Hit_N
27th-> 2020/02/05 to 2020/02/18, 28th-> 2020/02/19 to 2020/03/03

西域の過去1年Hit_P予測実績は80.0%、予測実績詳細は:

Hit_P : 12 期分が、Hit_P 状況にマッチしました
Hit_N : 11 期分が、Hit_N 状況にマッチしました
空振り : 3 期分が、空振り状況にマッチしました
見逃し : 発生せず
Hit_P 予測の実績は、12/(12+3) = 80.0%

 

● 南関東_監視域におけるM5.0以上の発震状況の予測です

南関東は、地震本部さん2020/01/01資料で、相模トラフ南関東直下地震M6.7〜7.3の確率70%、と予測されている領域です

以下、度数分布_発震予測確率のグラフです

f:id:yoshihide-sugiura:20200207024401p:plain

ここでは過去4年間に29回の発震があり(間隔を見るので複数同日発震は1となります、また発震間隔数は発震回数より1少なくなります)、発震間隔は、平均で48.4日、最短で1日、最長で132日です
最後の発震 2020/02/01 M5.3 Dph62km SOUTHERN IBARAKI PREF から、2020/02/05現在、4日が経過しています
現在、発震確率の状態は46%で、それは2020/03/13まで続きます

 

 

木星の衝合の説明:

木星の衝と合とは、国立天文台さん [暦Wiki/惑星/合と衝 - 国立天文台暦計算室] より、衝は外惑星についてのみ起こる現象で、太陽・地球・外惑星がその順に直線に並ぶ状態です

http://eco.mtk.nao.ac.jp/koyomi/wiki/CFC7C0B1C5B7BEDD.png

木星は約12年で太陽を公転しており、地球は1年で公転、衝は約1年に1回出現、正反対の合も約1年に1回出現、衝と合をあわせて2回/約1年出現で、衝の地球と木星の距離は約5.0au、合では約7.0au(auは太陽と地球の距離)、衝の日付は「ほんのり光房」さんのブログ:http://kuusou.asablo.jp/blog/2018/05/09/8847763 を参照し、合の日付は衝間の1/2を取りました

ここで、衝合日を含む±48日間の97日内で何%のM6.0以上地震が発震したか?を16年で調べます、衝と合で年2878日となり、2878/365x16=49.3%ですので、50%を大幅に超えれば木星衝合日を含む91日間とM6.0以上は相関があり、50%前後で相関はない、50%を大幅に下回れば逆相関がある、となります、4年の場合は50%を少し超えますが、16年を基準として±48日に設定しています

 

= まとめです =

 2017年の放出エネルギーは過去最低 、2018年は増加、2019年も増加でした [2019-12 地震の予測マップ ここ26年間の北西太平洋地域・地震放出エネルギー推移を測定する! - 地震の予測マップと時系列予測グラフ]

 ここで赤い救急マークは、M5.0以上の地震ヶ所をピンポイントで予測しています

 時系列予測では、各域において直近将来約14日間地震発生を予測しています

 どちらも、2019山形沖地震M6.8、2018胆振地震M6.7、2018大阪北部地震M6.1、2016熊本地震の前震M6.5、2011年東北太平洋地震の前震M7.3、を正しく予測しています

 

という北海道東部沖の千島海溝沿いで、東日本大震災のようなマグニチュード(M)9級の超巨大地震が、いつ起きても不思議はないとの見方を示した。
北海道東部沖の千島海溝沿いで、東日本大震災のようなマグニチュード(M)9級の超巨大地震が、いつ起きても不思議はないとの見方を示した。
政府の地震調査委員「超巨大地震:北海道沖、M9地震予測 本州にも被害の恐れ」といった報道をするのであれば、注意喚起領域でM5.5クラスが来た時には「注意喚起情報」を出されるべきではないか、と思います。
政府の地震調査委員会

最後まで読んで頂き、ありがとう御座いました

・ 東進西進の原理、東進西進の識別方法、等の説明はこちら [テクニカル事項]

・「地震の予測マップ」のデータ更新タイミングの説明はこちら [データ更新タイミング]

・ 国土地理院さん提供の地殻変動マップはこちら [地殻変動情報] javaがインストールされている必要があります

・「太陽黒点数の推移を追う!」は別ブログへ [太陽黒点数の推移を追う]

= 以上です =

 

謝辞: 本予測は「気象庁文部科学省が協力してデータを処理した結果」の「気象庁一元化処理震源要素」データ一年分(暫定)を「防災科学技術研究所」サイトよりダウンロードして解析しています。 [2] このデータによって初めて一般にリアルタイム解析が可能となったもので、構築にご尽力頂きました各国立大学、各官庁と関連する機関、都道府県と関連する機関、等の関係各位殿に深く謝意を述べさせて頂きます。

免責: 本予測は個人の推論によるもので、プログラムバグやデータ解釈ミス等も含め、ここで表示された結果について何ら責任を負うものではありません。

引用:

[1] スロースリップ - Wikipedia プレートがゆっくりと移動し大きな破壊を伴わずにエネルギー解放する現象ですが、プレート深部には応力歪が伝搬され、これが原因で深部のスラブ内領域で通常の地震が発生します。

[2] 防災科学技術研究所 Hi-net 高感度地震観測網

[3] 気象庁|震源データ

[4] 地震情報 - Yahoo!天気・災害

[5] 過去巨大地震マップ - 地震の予測マップ

[6] 国立天文台 太陽観測科学プロジェクト 三鷹太陽地上観測

[7] こよみの計算 - 国立天文台暦計算室