地震の予測マップと発震間隔ヒストグラム

地震発生ヶ所をポイント予測し、発震間隔ヒストグラムで将来を予測する

2019-12-18 地震の予測マップと時系列予測グラフ 19日の地震列島は青森東方沖でM5.5震度5弱, 茨城沖でM4.3, 青森西方沖でM3.7! 解説:南関東_監視域の発振分析・日単位の予測本文へ組込みを考える!

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

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

 地震の予測マップは、毎日更新です

 地震の時系列予測ARIMAグラフは、2019-12-11に更新されました、次回の更新は2019-12-26で、約14日間隔で更新されます

 

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

⚫ 12月18日は月齢21.0日 西進Day、月距離 37.0387万km

⚫ 12月19日は月齢22.0日 西進Day、月距離 37.0282万km 凹部ボトム

⚫ 12月20日は月齢23.0日 西進Day、月距離 37.0631万km

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

月距離とは地球から月までの距離で、約35〜40万kmの間を振動しています、38.7万kmが大体の滞在時間(日数)中間となる境界ラインで、上部を凸部、下部を凹部と称しています、この振動の周期は月の公転周期27.3日となっており、1年で13周期あり中央38.7万kmで分けた凸部凹部の数は26個/年(これを1年26期と言う)になり、1期の日数は約14日となります

 

 

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

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

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

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

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

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

 12月19日04時15分、青森西方沖でM3.7、深さ10km、震度1 中域_時系列 

 12月19日09時23分、茨城沖でM4.2、深さ90km、震度3 西域_時系列

 12月19日15時21分、青森東方沖でM5.5、深さ50km、震度5弱 中域_時系列

 

 19日は、何と青森西方沖でM3.7でした! ブログ開始して青森西方沖というのは始めてです、深さ10kmという事で北米プレート地殻内地震です、太平洋プレートスラブ内ではありませんでした、これで納得です

青森東方沖でM5.5、震度5弱が発生しました、深さ50kmですので太平洋プレートスラブ内地震です(中域_時系列ARIMAHit_Pを予測)、南関東_監視域の分析プログラム化に続き、青森東方沖での監視域設定を急ぎたい、と考えております、ここは地震本部さん資料で、青森東方沖及び岩手沖北部M7.0〜7.5の確率90% 、が予測されている領域で引き続き警戒が必要です

 

 

解説: 南関東_監視域の発振分析・日単位の予測本文へ組込みを考える!

 いよいよ南関東_監視域の地震分析をどうのようにして予測本文に組み込むべきか?の考察に入ります、ので、お付き合い頂けますでしょうか? よろしくお願い申し上げます

 

 まず、西域の境界を東方へ移動し、南関東_監視域と重ねました、マップを以下に示します、マークは南関東_監視域で分析の参考にしている評価サンプル地震です

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

一番北の茨城北部M6.3のみが中域_北米プレート地殻内で発生した地震、残りは西域_フィリピン海プレートスラブ内で発生した地震となります

 

 ここで、南関東_監視域は深さ30km以上100km未満、西域は深さ0km以上130km未満であり、完全に西域に包含されており、以下考察です:

 南関東で時系列_地震予測をすると、ARIMA西域_時系列予測に包含され、両者の確率論理(AND)を考えざるを得ない

 例えば現在西域_時系列ARIMA予測はHit_Nであり、この1年のHit_N予測実績は85.7%であり14.3%が「見逃し」比率で、この範囲内の南関東におけるHit_Pを予測する、という事になり、かつ、南関東の予測はやるとすれば日単位なのである

 これはとてもやってられない、予測に予測を被せて%を出す(これは自然と出てき来てしまう)のは現実論ではない

 やるとすれば、ARIMA予測で西域_時系列にHit_Pが出た場合に限り、南関東における発振分析を日単位で提示する、という事であろう、これで%論は一応回避できる

 それでも2日前のデータを使って、その日の発振分析を行なう、以上は出来なくて、、、これを毎日行なう程度の事しか出来ない

 

 こんな線が見えて来たので、南関東_監視域の発振分析・日単位をプログラム化する事にします、尚、ブログがドンドン長くなるので、現在、文章をドンドン短くしています(これでも、ですが)

 

以上、分かりづらい内容にお付き合い頂き、誠にありがとう御座います、感謝です

 

 

   

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

 

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

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

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

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

= 地震の予測マップ・ピッチ36kmマップとM6.0以上地震トレンド分析グラフです =

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

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

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

 

M6.0以上地震発生トレンド分析・南海トラフ監視領域_時系列グラフ4年間です

南海トラフ監視域の発生地震から、各域M6.0以上地震との相関を分析します

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

上図のX,Vは、Wiki掲載 [地震の年表 (日本) - Wikipedia] 2016年以降地震
・ 2016/04/01 南海トラフM6.5、深さ28km
・ 2019/05/10 日向灘M6.3、深さ25km
を加えたものです
4年間各域M6.0以上と総数は、{:East=>5, :Midl=>16, :West=>11}、TOTAL-> 32
カウント4でM6.0以上の発振数は2個と少なく、全体比率-> 6.3%
それは、以下の中域_時系列における2件で3.11の余震と思われます
・ 発生日:2019/08/04 中域_時系列、M6.4、OFF FUKUSHIMA PREF
・ 発生日:2019/08/29 中域_時系列、M6.1、OFF AOMORI PREF

分析: 巨大地震は若いカウントで発振し、その後、カウント4に至る事はあっても、カウント4で巨大地震が初期発振する事は無い、及び、初期巨大地震はカウント0からカウント3へ上昇中に発振し、カウント4からカウント0に至る下降中に発生する事はない、と分析できます

又、直近7日間の南海トラフ監視領域カウントは:
2019/12/12 カウント-> 3
2019/12/13 カウント-> 3
2019/12/14 カウント-> 3
2019/12/15 カウント-> 3
2019/12/16 カウント-> 3
2019/12/17 カウント-> 3
2019/12/18 カウント-> 3
となっており、カウント4から3に下がり始めた状態で、今後しばらくの間、M6.0以上の初期巨大地震は発生しないと思われます(日本海溝において発生する3.11起因の余震を除く、従って栃木北部と茨城北部は警戒が必要)
ここで初期巨大地震とはM6.0以上の地震で、誘発地震でもなく余震でもなく、数十年かそれ以上に一度発生する巨大地震の事で、3.11・2016熊本・2016鳥取中部・2018大阪北部・2018胆振・2019山形沖、が相当し私の造語で一般用語ではなく、地震とは、初期巨大地震に始まり誘発地震と余震を繰り返しやがて次の初期巨大地震に至る一連の過程である、という事になります

 

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

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

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

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

f:id:yoshihide-sugiura:20191220013717p: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年間確率、以下同様

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

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

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

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

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

東域M5.0以上発生予測 : 続く27/28期の予測は、 27th-> Hit_N, 28th-> Hit_P
27th-> 2019/12/12 to 2019/12/26, 28th-> 2019/12/27 to 2020/01/07

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

  Hit_P : 11 期分が、Hit_P 状況にマッチしました
  Hit_N : 9 期分が、Hit_N 状況にマッチしました
空振り : 5 期分が、空振り状況にマッチしました
見逃し : 1 期分が、見逃し状況にマッチしました
Hit_N 予測の実績は、9/(9+1) = 90.0%

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:20191220013755p:plain

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

f:id:yoshihide-sugiura:20191220013812p: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% [海溝で起こる地震 | 地震本部]

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

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

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

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

中域M5.0以上発生予測 : 続く27/28期の予測は、 27th-> Hit_P, 28th-> Hit_P
27th-> 2019/12/12 to 2019/12/26, 28th-> 2019/12/27 to 2020/01/07

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

  Hit_P : 10 期分が、Hit_P 状況にマッチしました
  Hit_N : 10 期分が、Hit_N 状況にマッチしました
空振り : 3 期分が、空振り状況にマッチしました
見逃し : 3 期分が、見逃し状況にマッチしました
Hit_P 予測の実績は、10/(10+3) = 76.9%

 

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

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

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

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

南海トラフ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:20191220013921p:plain

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

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

西域M5.0以上発生予測 : 続く27/28期の予測は、 27th-> Hit_N, 28th-> Hit_P
27th-> 2019/12/12 to 2019/12/26, 28th-> 2019/12/27 to 2020/01/07

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

  Hit_P : 11 期分が、Hit_P 状況にマッチしました
  Hit_N : 12 期分が、Hit_N 状況にマッチしました
空振り : 1 期分が、空振り状況にマッチしました
見逃し : 2 期分が、見逃し状況にマッチしました
Hit_N 予測の実績は、12/(12+2) = 85.7%

 

= まとめです =

 2017年の放出エネルギーは過去最低 、2018年は増加、2019年も現在増加しています [2019-10 地震の予測マップ ここ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] こよみの計算 - 国立天文台暦計算室