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

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

2019-11-02 地震の予測マップと時系列予測グラフ 3日の地震列島は福島沖でM5.1, 宮城沖でM4.3, 千葉,宮城沖でM3.X! 解説:東域と西域の範囲を調整する! 本日の予測マップ&グラフは?

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

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

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

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

 

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

・ 11月02日は月齢4.5日 東進Day、月距離 38.9047万km 38.7万kmを超えた日

・ 11月03日は月齢5.5日 東進Day、月距離 39.4327万km

・ 11月04日は月齢6.5日 東進Day、月距離 39.8741万km 右肩上がり上昇ライン

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

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

月の地球に対する公転周期は27.3日で、新月から次の新月前日までの月齢一月は29.5日と、月齢の方が2日少々長いのは、地球の太陽に対する公転運動を月の公転運動が追いかけている為です

 

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

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

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

 

 

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

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

 11月3日00時33分、宮城沖でM4.3、深さ30km、震度2

 11月3日02時39分、千葉北西部でM3.8、深さ110km、震度1

 11月3日12時30分、宮城沖でM3.3、深さ70km、震度1 本日二回目

 11月3日14時02分、福島沖M5.1、深さ30km、震度3

 

  3日、宮城沖と千葉内陸部での地震が収まりません

早速ですが、福島沖のM5.1は、空間予測(救急マーク)、時系列予測(11月2日〜11月14日中域でHIt_Pを予測、的中率の過去実績は68.0% )に当てはまっています

 

 

解説: 東域と西域の範囲を調整する!

 東域の範囲を東へシフトさせ国後島近海を含むよう変更、及び、西域の範囲を南へ多少シフトさせました、ので、ご説明致したく、お付き合い頂けますよう、よろしくお願い致します

 

 結果:

東域については、M4.8以上の地震数は34 ➡ 46に増加し、Hit_P的中率は57.7% ➡ 61.5%に改善されました、ARIMA指数(p,d,q)は同じ(0,0,0)でした

西域については、M4.3以上の地震数は37 ➡ 41に増えましたが、Hit_P的中率は変化なしの72.0%でした、ARIMA指数も同じ(0,0,0)です

多少作為的な感じもありますが、まぁ、良しとしておきましょう、本文の東域と西域をアップデートしてあります

 

 やはりARIMA指数(p,d,q)についてもう少し理解を深める必要がありそうです、dは定常状態に至る次数で定常状態に至った時の、pが連続性の次数、qが周期性の次数、という事のようですが、ここで周期性は良いのですが、問題は連続性です、二つ問題があります

・ 昨日はゼロで、今日もゼロだった、で、明日もゼロだろうか?

地震の場合、明日は起きる可能性が高くなります

・ 昨日は起きて、今日も起きた、で、明日は?

これは、地震の場合、本当に分からない!

これは内部エネルギー状態を持たないARIMAモデルの限界を示しており、Rのforecast機能がこれらの連続性にどう対処しているのか?を見極める必要がありそうです

前者のゼロ入力の場合、感度を上げてゼロ連続入力を少なくしてしまう方法を現在は採用し、問題を回避しようとしていますが、西域では二ヶ所ほど連続ゼロ入力が出てしまっています、ARIMAのグラフを見て頂けると分かります、もう少し感度を上げる必要があるのでしょうか?

ここで後者については特に手を打っていません、ARIMAにお任せ、です

この辺りは、調べます、なかなか南海トラフ深層横断地震帯へは行けそうにありませんが、ここはキチッと把握しておく必要があります

 

 尚、当面次の更新2019-11-14まで、範囲の微調整と感度の微調整が続く可能性があります

 

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

 

 

 

 

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

 

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

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

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

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

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

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

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

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

 

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

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

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

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

f:id:yoshihide-sugiura:20191104001550p: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 ["45.7%", "54.3%", 46] Days ["47.7%", "52.3%", 365]

東域でこの1年に発生したM4.8以上は46個で、45.7%と54.3%の比率で38.7万以上(凸部)とそれ未満(凹部)に別れる、この1年365日の凸部凹部の比率は47.7%と52.3%で、凹部での地震発生比率が少し高い(偏っている事が分かる

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

時系列地震予測グラフ(Rのforecast機能)東域は、2019-11-01 の予測結果で、

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

ミニ解説・東域

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

Forecast East M>=4.8: next period 27th-> 1.7, next of the next 28th-> 1.7
27th(top)-> 2019/11/02 to 2019/11/14, 28th(btm)-> 2019/11/15 to 2019/11/29
Selected ARIMA Best model index-> (0,0,0) out of 6 trials

東域では11月2日〜11月14日迄の間にM4.8以上が1回、11月15日から11月29日の間でも1回発生するだろう、と予測しています、(0,0,0)は白色雑音に近くwith non-zero meanが出た場合は単純平均となります、下記の過去直近1年の予測実績より、M4.8以上が発生すると予測された期でM5.0以上が発生する確率は61.5%です、次の更新は 2019-11-14 です

過去1年26期のM5.0発生実績と予測結果を比較する抽象化サマリーは以下の通りで、

抽象化とは「来ると予測してM5.0以上が来た(これがHIt_P)」場合その期をONとする数え方で個数を問わない、Hit_P Rateとは、Hit_P数を、Hit_P数と空振りKaraburi数を足した数で割った%、即ちスウィングして当たった%です

過去直近1年の予測実績を示す抽象化サマリー東域は、M5.0以上の発生事例に対し

Hit_P : 16 periods Matched with Hit_P condition
Hit_N : Nothing, not detected
Karaburi : 10 periods Matched with Karaburi condition
Minogashi : Nothing, not detected
Hit_P Rate of the East in the year 2018/11/09 ... 2019/11/01-> 16/(16+10) = 61.5%

東域の過去的中率、M5.0以上について過去1年間2018/11/09 ... 2019/11/01までの26期に渡るHit_P Rate、は61.5%です

 

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

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

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

f:id:yoshihide-sugiura:20191104001653p: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 ["42.2%", "57.8%", 45] Days ["47.7%", "52.3%", 365]

中域でこの1年に発生したM4.8以上は45個で、42.2%と57.8%の比率で38.7万以上(凸部)とそれ未満(凹部)に別れる、この1年365日の凸部凹部の比率は47.7%と52.3%で、凹部での地震発生比率が少し高い(偏っている)事が分かる

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

時系列地震予測グラフ(Rのforecast機能)中域は、2019-11-01 の予測結果で、

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

ミニ解説・中域

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

Forecast Midl M>=4.8: next period 27th-> 1.0, next of the next 28th-> 1.0
27th(top)-> 2019/11/02 to 2019/11/14, 28th(btm)-> 2019/11/15 to 2019/11/29
Selected ARIMA Best model index-> (0,1,1) out of 18 trials

中域では11月2日〜11月14日迄の間にM4.8以上が1回、11月15日〜11月29日迄の間でも1回発生するだろう、と予測しています、(0,1,1)は18回のトライアルの結果選択され、1,1は、1次で定常状態化しその際1次の周期性が検出された予測を意味しています、下記の過去直近1年の予測実績より、M4.8以上が発生すると予測された期でM5.0以上が発生する確率は68.0%です、次の更新は 2019-11-14 です

過去1年26期のM5.0発生実績と予測結果を比較する抽象化サマリーは以下の通りで、

抽象化とは「来ると予測してM5.0以上が来た(これがHIt_P)」場合その期をONとする数え方で個数を問わない、Hit_P Rateとは、Hit_P数を、Hit_P数と空振りKaraburi数を足した数で割った%、即ちスウィングして当たった%です

過去直近1年の予測実績を示す抽象化サマリー中域は、M5.0以上の発生事例に対し

Hit_P : 17 periods Matched with Hit_P condition
Hit_N : 1 periods Matched with Hit_N condition
Karaburi : 8 periods Matched with Karaburi condition
Minogashi : Nothing, not detected
Hit_P Rate of the Midl in the year 2018/11/09 ... 2019/11/01-> 17/(17+8) = 68.0%

中域の過去的中率、M5.0以上について過去1年間2018/11/09 ... 2019/11/01までの26期に渡るHit_P Rate、は68.0%です

 

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

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

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

f:id:yoshihide-sugiura:20191104001758p: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 ["41.5%", "58.5%", 41] Days ["47.7%", "52.3%", 365]

西域でこの1年に発生したM4.3以上は41個で、41.5%と58.5%の比率で38.7万以上(凸部)とそれ未満(凹部)に別れる、この1年365日の凸部凹部の比率は47.7%と52.3%で、凹部での地震発生比率は少し高い(偏っている)事が分かる

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

時系列地震予測グラフ(Rのforecast機能)西域は、2019-11-01 の予測結果で、

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

ミニ解説・西域

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

Forecast West M>=4.3: next period 27th-> 1.5, next of the next 28th-> 1.5
27th(top)-> 2019/11/02 to 2019/11/14, 28th(btm)-> 2019/11/15 to 2019/11/29
Selected ARIMA Best model index-> (0,0,0) out of 6 trials

西域では11月2日から11月14日の間でM4.3以上が1回、11月15日から11月29日の間でも1回発生するだろう、と予測しています、(0,0,0)は白色雑音に近くwith non-zero meanが出た場合は単純平均となります、下記の過去直近1年の予測実績より、M4.8以上が発生すると予測された期でM5.0以上が発生する確率は72.0%です、次の更新は 2019-11-14 です

過去1年26期のM5.0発生実績と予測結果を比較する抽象化サマリーは以下の通りで、

抽象化とは「来ると予測してM5.0以上が来た(これがHIt_P)」場合その期をONとする数え方で個数を問わない、Hit_P Rateとは、Hit_P数を、Hit_P数と空振りKaraburi数を足した数で割った%、即ちスウィングして当たった%です

過去直近1年の予測実績を示す抽象化サマリー西域は、M5.0以上の発生事例に対し

Hit_P : 18 periods Matched with Hit_P condition
Hit_N : Nothing, not detected
Karaburi : 7 periods Matched with Karaburi condition
Minogashi : 1 periods Matched with Minogashi condition
Hit_P Rate of the West in the year 2018/11/09 ... 2019/11/01-> 18/(18+7) = 72.0%

西域の過去的中率、M5.0以上について過去1年間2018/11/09 ... 2019/11/01までの26期に渡るHit_P Rate、は72.0%です

 

= まとめです =

 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] こよみの計算 - 国立天文台暦計算室