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

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

2019-10-21 地震の予測マップ 22日の地震列島は釧路沖でM4.Xが2件! 解説:復旧致しました、又、新規な評価プログラムを立ち上げました! 本日の予測マップは?

皆様、ご心配をお掛け致しました、無事復旧致しました、皆様の暖かいご声援に心より感謝申し上げます、又、新規な評価プログラムも立ち上げました、詳細は解説をご参照下さい、今後ともよろしくお願い申し上げます

 

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

・ 地震予測マップは、毎日、午前1時〜午前3時頃にアップしています

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

 

[こよみの計算 - 国立天文台暦計算室] 、[月の地心座標 - 国立天文台暦計算室] より:

・ 10月21日は月齢21.9日 西進Day、月距離 37.7223万kmで凹部に既に入ってます

・ 10月22日は月齢22.9日 西進Day、月距離 37.2828万km

・ 10月23日は月齢23.9日 西進Day、月距離 36.8705万kmで右肩下がりライン上です

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

月距離とは地球から月までの距離、約35〜40万kmの間を振動しています、38.7万kmとは大体の中間ラインです

 

経験上、千島海溝、日本海溝、でM6.5以上の西進破壊による巨大地震は:

・ 赤マークで覆われている領域で、西進Dayに発生します

・ 青マークで覆われている領域で、東進Dayに発生します

 

 

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

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

 10月22日14時51分、釧路沖でM4.1、深さ50km、震度2

 10月22日17時06分、釧路沖でM4.4、深さ30km、震度2

 

 22日は釧路沖でM4.Xが2件でした

 

 

 

解説: 復旧致しました、又、新規な評価プログラムを立ち上げました!

 皆様、ご心配をお掛けして申し訳ありませんでした、思ったより短期間で復旧致しました、これはバックアップしていたソースが、開発が一段落して安定化した状態であった為、月距離地震グラフや時系列地震予測グラフは入っていなかったものの、それらは記憶に新しく、短期間で復旧できた事によります

復旧に際して、大きなバグを発見し、これを修正した事と、新たに時系列地震予測の評価プログラムを作成致しましたので、結果をご報告致したく、お付き合い頂けますよう、よろしくお願い申し上げます

 

 まず、バグですが、Rのforecast機能プログラムにバグがありました、凸凹時系列地震予測グラフにおいて、入力を26期入れても25期しか現れないのは、私のバグでした、今回ソースが消えたのでゼロから作り直したら26期現れました、理由は細かい事なので述べませんが、予測は27期と28期の2期となります

本日の凸凹時系列地震予測グラフは、2019-10-18データによる予測をバグ修正版にて置き換えています

 

 続いて、予測結果の評価です、2019-10-04データから逆上ること26期分の凸凹時系列地震予測グラフを走行させ、この1年に発生したM5.0以上地震との比較を行いました、時系列地震予測を26回走行させた結果を評価しています

比較は、以下の4種類に分けています

HitP・・・来ると予測して、実際にM5.0以上が発生した

HitN・・・来ないと予測して、実際にM5.0以上は発生しなかった

Karaburi・・・来ると予測して、実際にはM5.0以上は発生しなかった

Minogashi・・来ないと予測して、実際にはM5.0以上が発生した

 

結果は長いリストと短いサマリーリストの2種を出してみました

 結果です、少し長いリストですが、東域、中域、西域に分けて:

Area-> East
HitP :
2019/08/24 to 2019/09/06: >=M5.0 ... 1 times_happen and 1.7 times_estimated
2019/08/10 to 2019/08/23: >=M5.0 ... 1 times_happen and 1.7 times_estimated
2019/07/27 to 2019/08/09: >=M5.0 ... 1 times_happen and 1.8 times_estimated
2019/06/29 to 2019/07/13: >=M5.0 ... 1 times_happen and 1.8 times_estimated
2019/06/17 to 2019/06/28: >=M5.0 ... 1 times_happen and 1.8 times_estimated
2019/06/02 to 2019/06/16: >=M5.0 ... 1 times_happen and 1.8 times_estimated
2019/05/06 to 2019/05/20: >=M5.0 ... 1 times_happen and 1.8 times_estimated
2019/04/24 to 2019/05/05: >=M5.0 ... 2 times_happen and 1.7 times_estimated
2019/04/09 to 2019/04/23: >=M5.0 ... 3 times_happen and 1.8 times_estimated
2019/02/26 to 2019/03/12: >=M5.0 ... 2 times_happen and 2.1 times_estimated
2019/02/13 to 2019/02/25: >=M5.0 ... 1 times_happen and 1.8 times_estimated
2019/01/16 to 2019/01/28: >=M5.0 ... 1 times_happen and 1.8 times_estimated
2018/12/18 to 2018/12/31: >=M5.0 ... 1 times_happen and 2.0 times_estimated
2018/11/21 to 2018/12/05: >=M5.0 ... 1 times_happen and 2.4 times_estimated
2018/10/24 to 2018/11/08: >=M5.0 ... 4 times_happen and 1.9 times_estimated
HitN : Nothing, not detected
Karaburi :
2019/10/05 to 2019/10/18: >=M5.0 Not_happen but 1.4 times_estimated
2019/09/21 to 2019/10/04: >=M5.0 Not_happen but 1.5 times_estimated
2019/09/07 to 2019/09/20: >=M5.0 Not_happen but 1.5 times_estimated
2019/07/14 to 2019/07/26: >=M5.0 Not_happen but 1.8 times_estimated
2019/05/21 to 2019/06/01: >=M5.0 Not_happen but 1.8 times_estimated
2019/03/27 to 2019/04/08: >=M5.0 Not_happen but 2.5 times_estimated
2019/03/13 to 2019/03/26: >=M5.0 Not_happen but 1.9 times_estimated
2019/01/29 to 2019/02/12: >=M5.0 Not_happen but 1.8 times_estimated
2019/01/01 to 2019/01/15: >=M5.0 Not_happen but 1.9 times_estimated
2018/12/06 to 2018/12/17: >=M5.0 Not_happen but 2.3 times_estimated
2018/11/09 to 2018/11/20: >=M5.0 Not_happen but 1.0 times_estimated
Minogashi : Nothing, not detected

Area-> Midl
HitP :
2019/08/24 to 2019/09/06: >=M5.0 ... 2 times_happen and 1.5 times_estimated
2019/08/10 to 2019/08/23: >=M5.0 ... 1 times_happen and 2.0 times_estimated
2019/07/27 to 2019/08/09: >=M5.0 ... 1 times_happen and 2.1 times_estimated
2019/07/14 to 2019/07/26: >=M5.0 ... 1 times_happen and 2.1 times_estimated
2019/06/29 to 2019/07/13: >=M5.0 ... 2 times_happen and 1.9 times_estimated
2019/06/17 to 2019/06/28: >=M5.0 ... 3 times_happen and 1.5 times_estimated
2019/05/21 to 2019/06/01: >=M5.0 ... 1 times_happen and 1.7 times_estimated
2019/04/09 to 2019/04/23: >=M5.0 ... 2 times_happen and 3.4 times_estimated
2019/02/26 to 2019/03/12: >=M5.0 ... 4 times_happen and 2.6 times_estimated
2019/02/13 to 2019/02/25: >=M5.0 ... 1 times_happen and 2.2 times_estimated
2019/01/29 to 2019/02/12: >=M5.0 ... 1 times_happen and 1.6 times_estimated
2019/01/16 to 2019/01/28: >=M5.0 ... 2 times_happen and 1.8 times_estimated
2018/12/18 to 2018/12/31: >=M5.0 ... 4 times_happen and 1.8 times_estimated
2018/12/06 to 2018/12/17: >=M5.0 ... 1 times_happen and 1.1 times_estimated
2018/11/21 to 2018/12/05: >=M5.0 ... 3 times_happen and 3.5 times_estimated
2018/10/24 to 2018/11/08: >=M5.0 ... 3 times_happen and 2.6 times_estimated
HitN :
2018/11/09 to 2018/11/20: >=M5.0 Not_happen and 0.9 times_estimated
Karaburi :
2019/10/05 to 2019/10/18: >=M5.0 Not_happen but 1.3 times_estimated
2019/09/21 to 2019/10/04: >=M5.0 Not_happen but 1.4 times_estimated
2019/09/07 to 2019/09/20: >=M5.0 Not_happen but 1.9 times_estimated
2019/06/02 to 2019/06/16: >=M5.0 Not_happen but 2.6 times_estimated
2019/05/06 to 2019/05/20: >=M5.0 Not_happen but 2.6 times_estimated
2019/04/24 to 2019/05/05: >=M5.0 Not_happen but 3.0 times_estimated
2019/03/27 to 2019/04/08: >=M5.0 Not_happen but 3.5 times_estimated
2019/03/13 to 2019/03/26: >=M5.0 Not_happen but 1.6 times_estimated
2019/01/01 to 2019/01/15: >=M5.0 Not_happen but 1.4 times_estimated
Minogashi : Nothing, not detected

Area-> West
HitP :
2019/09/21 to 2019/10/04: >=M5.0 ... 1 times_happen and 1.3 times_estimated
2019/08/24 to 2019/09/06: >=M5.0 ... 1 times_happen and 2.1 times_estimated
2019/07/27 to 2019/08/09: >=M5.0 ... 2 times_happen and 2.1 times_estimated
2019/07/14 to 2019/07/26: >=M5.0 ... 1 times_happen and 1.8 times_estimated
2019/06/29 to 2019/07/13: >=M5.0 ... 1 times_happen and 1.7 times_estimated
2019/06/17 to 2019/06/28: >=M5.0 ... 2 times_happen and 1.5 times_estimated
2019/05/06 to 2019/05/20: >=M5.0 ... 4 times_happen and 1.4 times_estimated
2019/03/27 to 2019/04/08: >=M5.0 ... 3 times_happen and 1.8 times_estimated
2019/03/13 to 2019/03/26: >=M5.0 ... 1 times_happen and 1.7 times_estimated
2019/02/26 to 2019/03/12: >=M5.0 ... 3 times_happen and 1.6 times_estimated
2019/02/13 to 2019/02/25: >=M5.0 ... 1 times_happen and 1.6 times_estimated
2019/01/29 to 2019/02/12: >=M5.0 ... 1 times_happen and 1.7 times_estimated
2019/01/16 to 2019/01/28: >=M5.0 ... 1 times_happen and 1.6 times_estimated
2019/01/01 to 2019/01/15: >=M5.0 ... 1 times_happen and 1.5 times_estimated
2018/12/18 to 2018/12/31: >=M5.0 ... 4 times_happen and 1.6 times_estimated
2018/12/06 to 2018/12/17: >=M5.0 ... 1 times_happen and 1.6 times_estimated
2018/11/21 to 2018/12/05: >=M5.0 ... 2 times_happen and 1.6 times_estimated
2018/10/24 to 2018/11/08: >=M5.0 ... 3 times_happen and 1.7 times_estimated
HitN : Nothing, not detected
Karaburi :
2019/10/05 to 2019/10/18: >=M5.0 Not_happen but 1.2 times_estimated
2019/09/07 to 2019/09/20: >=M5.0 Not_happen but 1.3 times_estimated
2019/08/10 to 2019/08/23: >=M5.0 Not_happen but 1.9 times_estimated
2019/06/02 to 2019/06/16: >=M5.0 Not_happen but 1.1 times_estimated
2019/04/24 to 2019/05/05: >=M5.0 Not_happen but 1.4 times_estimated
2019/04/09 to 2019/04/23: >=M5.0 Not_happen but 1.4 times_estimated
2018/11/09 to 2018/11/20: >=M5.0 Not_happen but 1.8 times_estimated
Minogashi :
2019/05/21 to 2019/06/01: >=M5.0 ... 1 times_happen but No 0.0 times_estimated

 

 次が抽象化したサマリーリストです:

2018/10/24 to 2018/11/08 ... 2019/10/05 to 2019/10/18, in total 26 periods evaluated
Area-> East
           HitP : 15 periods Matched with HitP condition
           HitN : Nothing, not detected
   Karaburi : 11 periods Matched with Karaburi condition
Minogashi : Nothing, not detected

Area-> Midl
           HitP : 16 periods Matched with HitP condition
           HitN : 1 periods Matched with HitN condition
   Karaburi : 9 periods Matched with Karaburi condition
Minogashi : Nothing, not detected

Area-> West
           HitP : 18 periods Matched with HitP condition
           HitN : Nothing, not detected
   Karaburi : 7 periods Matched with Karaburi condition
Minogashi : 1 periods Matched with Minogashi condition

抽象化とは、予測件数や発生地震件数は表示せず、期間単位で予測があった、地震があった、というON/OFFレベルでまとめるものです

上記サマリーリストでの件数は期間の数です、ある期間で地震の予測が1.0以上であれば「予測されていた」で、M5.0以上地震の発生件数も期間内で1件以上あれば「実際に発生した」となり、Matchedとしています、抽象化サマリーリストでは詳細な情報は欠落させています

各領域の Matched 期間数を足せば、26期間となります(1年間の予測実績となります、分かりやすいと思います)

ここで、HitP vs Karaburi が:

東域で 15 vs 11

中域で 16 vs 9

西域で 18 vs 7

というのはナカナカ良い成績であるように思われます、私は Karaburi だらけになるのでは?と思っていましたので

中域でHitNが1件あるのは立派と思います、が、西域で Minogashi が1件あるので、これでチャラでしょうか?

あと、サマリーリストに%表示を出すようにします、結局、時系列地震予測の実績は%で表示する事になると思います、ここで%表示の分子は明確ですが分母を何とするか?は少し落ち着いて考えないといけません

 

以上、読みづらい長文にお付き合い頂き、誠にありがとう御座いました

 

 

 

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

 

= 以下、2019-10-21 迄データ1年分による本日の地震の予測マップ =

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

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

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

= 地震の予測マップ・ピッチ36kmマップです =

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

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

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

 

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

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

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

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

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

根室沖M7.8〜8.5の確率80%、青森東方沖及び岩手沖北部M7.0〜7.5の確率90% 、宮城沖M7.0〜M7.5の確率90%、福島沖M7.0〜M7.5の確率50% [海溝で起こる地震 | 地震本部] 2019年1月1日算定基準日の30年間確率、以下同様

月距離地震グラフ東域は、

38.7万km以上と未満: East ["41.0%", "59.0%", 39] Days ["47.9%", "52.1%", 365]

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

暫定:凸凹時系列地震予測グラフ東域は、2019-10-18 の予測結果で、

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

ミニ解説・東域

38.7万kmを境目とし1年の凸部凹部を26期で示し、続く27期と28期を予測、予測値実数は整数化し判断

Forecast East M>=4.8: next period 27th-> 1.4, next of the next 28th-> 1.4
27th(btm)-> 2019/10/19 to 2019/11/01, 28th(top)-> 2019/11/02 to 2019/11/14
Selected ARIMA Best model index-> (0,0,0)

東域では10月19日〜11月1日迄の間にM4.8以上が1回、11月2日から11月14日の間でも1回発生するだろう、と予測しています

次の更新は 2019-11-01 です

 

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

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

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

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

青森東方沖及び岩手沖北部M7.0〜7.5の確率90% 、宮城沖M7.0〜M7.5の確率90%、福島沖M7.0〜M7.5の確率50%、茨城沖M7.0〜M7.5の確率80%、相模トラフ南関東直下地震M6.7〜7.3の確率70% [海溝で起こる地震 | 地震本部]

月距離地震グラフ中域は、

38.7万km以上と未満: Midl ["40.8%", "59.2%", 49] Days ["47.9%", "52.1%", 365]

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

暫定:凸凹時系列地震予測グラフ中域は、2019-10-18 の予測結果で、

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

ミニ解説・中域

38.7万kmを境目とし1年の凸部凹部を26期で示し、続く27期と28期を予測、予測値実数は整数化し判断

Forecast Midl M>=4.8: next period 27th-> 1.2, next of the next 28th-> 1.1
27th(btm)-> 2019/10/19 to 2019/11/01, 28th(top)-> 2019/11/02 to 2019/11/14
Selected ARIMA Best model index-> (0,1,2)

中域では10月19日〜11月1日迄の間にM4.8以上が1回、11月2日〜11月14日迄の間でも1回発生するだろう、と予測しています

次の更新は 2019-11-01 です

 

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

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

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

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

南海トラフM8〜M9クラスの確率70%〜80%、日向灘M7.1前後の確率70〜80% [海溝で起こる地震 | 地震本部]

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

月距離地震グラフ西域は、

38.7万km以上と未満: West ["29.4%", "70.6%", 34] Days ["47.9%", "52.1%", 365]

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

暫定:凸凹時系列地震予測グラフ西域は、2019-10-18 の予測結果で、

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

ミニ解説・西域

38.7万kmを境目とし1年の凸部凹部を26期で示し、続く27期と28期を予測、予測値実数は整数化し判断

Forecast West M>=4.3: next period 27th-> 1.3, next of the next 28th-> 1.3
27th(btm)-> 2019/10/19 to 2019/11/01, 28th(top)-> 2019/11/02 to 2019/11/14
Selected ARIMA Best model index-> (0,0,0)

西域では10月19日から11月1日の間でM4.3以上が1回、11月2日から11月14日の間でも1回発生するだろう、と予測しています

次の更新は 2019-11-01 です 

 

= まとめです =

 2017年の放出エネルギーは過去最低 、2018年は増加、2019年も8月現在増加しています [2019-09 地震の予測マップ ここ26年間の北西太平洋地域・地震放出エネルギー推移を観測する! - 地震の予測マップ]

 2019年がどうなるか? もうしばらく見守る必要があります

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

 

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

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

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

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

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

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

= 以上です =

 

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

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

引用:

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

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

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

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

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

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

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