2018-01-02 地震の予測マップ (北海道中央・東北太平洋側・関東・鳥取・岡山・香川を除く日本全国が注意対象)

東進マーク西進マークの圧力ポイントを識別する地震予測です。 東進西進どちらも危ないです。 東進西進識別方法は一番下に記載(モデルは現在模索中です)。

マグネチュード M6.5以上の地震を予測します。

2018年1月2日まで一年間の統計解析結果・45kmピッチ概要マップを示します。 [1,2]

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

海溝型向け45kmピッチ概要マップです。 南西諸島・小笠原諸島など陸地直下型の起きない諸島部予測におかれましては上図をご参照下さい。

マークは東に向かう圧力、マークは西に向かう圧力がある地点です。

白い飛行機マークは30日以内に発生したM5.5以上地震です。

オレンジ色飛行機マークは31日以前かつ一年以内にM5.5以上地震があった地点です。

飛行機マーク周囲にマークが出ている所は続いて起きる地震に注意です。 飛行機マーク周辺にマークが出ていなければ、エネルギーは解放され小康状態に入っています。

次は9kmピッチ詳細マップ・東域です。 陸地直下型はこちらをご参照下さい。 凡例は概要と同じです。

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

以下、9kmピッチ詳細マップ・中域です。

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

以下、9kmピッチ詳細マップ・西域です。

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

マークはM6.5以上クラスの地震を起こす可能性のある注意ポイントを示します。 あくまでも可能性で注意喚起が目的です。 M6.5以上クラスの地震は年に一度起こすか起こさないか、の頻度です。 しかし、巨大地震・過去17事例において、M6.5以上クラスはすべてマークの在る場所またはその境界で起きています。

 

45km概要マップで昨日との比較ですが、大きな変化はありませんが、伊豆半島沖に一連の青マークラインが出現しました。 以下、9km詳細マップをご参照下さい。

まず北海道周辺部には赤青マークが出ており要注意です。

東北地方では東北地方太平洋沖地震の大きな余震は収まったようで赤青マークが消えました。 但し、M5.5クラスの弱い余震はマークなしでも起きます。 また青森・秋田と日本海側には注意ポイントが出ています。

関東(一部)・鳥取・岡山・香川では赤青マークが消えていますが、その他、マークが出ている部分は日本全国すべて要注意です。 1月2日19時48分頃、東京湾でM4.2の地震がありました。 東京湾には全く赤青マークがなく、M6.5に繋がるような前震型地震とは思えません。 過去事例では全く赤青マークの出ていない所にいきなりM6.5クラスが来た事はありません。

赤青マークで覆われた南海トラフと連結する陸地部分と日向灘は長期に渡り要注意です。 九州もまだ大半が赤マークで覆われています。

相模トラフも赤青マークが出ていますが、以前は青マークに覆われて見えなかった伊豆七島がよく見えるようになってきました。 関東大震災震源領域が赤青マークで覆われたら要注意です(まだですが)。

各トラフ・海溝と関東大震災(大正・元禄)震源領域を示します(Wikipediaより引用)。

f:id:yoshihide-sugiura:20171223144936p:plain f:id:yoshihide-sugiura:20171207001138p:plain f:id:yoshihide-sugiura:20171207001059p:plain

琉球海溝                   南海トラフ                             相模トラフ

f:id:yoshihide-sugiura:20171223095633p:plain f:id:yoshihide-sugiura:20171221230842p:plain f:id:yoshihide-sugiura:20171209000948p:plain

日本海溝               千島海溝                                  大正・元禄 関東大震災震源

過去の巨大地震マップ(M7.0以上)は大変役立ちます(これは貴重な財産です)

日本の主な地震の震央。赤:M 7以上青:死者有り紫:最大震度6以上。 [3]

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

2017年・地震の放出エネルギーはこの24年間で過去最低

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

一番右は2017年1月1日〜2017年12月31日までの一年で、この24年間で最低となりました。 問題は、この低レベルがいつ上昇に転ずるか、です。 太陽サイクルリストは: [4] 

2017年11月・太陽黒点数はこの48ヶ月で過去最低

三鷹太陽地上観測さん測定の月平均黒点数・直近48ヶ月を表示(©国立天文台)。 [5]

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

一番右が11月で黒点数は48ヶ月で最低でした。 放出エネルギーとの関係を追っています。

次のM6.5以上クラスはいつ頃来るのだろうか?

ラス前のM6.5以上は2016年10月21日の鳥取県中部地震(M6.6)で、ラストは2016年11月22日の福島県沖(M7.4、東北地方太平洋沖地震の余震)でした。 この一年M6.5以上は起きていません。 放出エネルギーはこの1年間最低のレベルでしたので、今後1〜2年の内には上昇に転じM6.5以上クラスが来るものと思われます。

 

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

明日も報告致します。 (よろしくお願い致します)

東進・西進の識別には月齢を使います。 [6] 月齢は新月から満月まで(上弦月)を3分割、満月から新月まで(下弦月)を3分割、合計6分割します。 地震があれば必ずどこかの分割スロットに入ります。

あるグリッド・ポイントについて上記6スロットを準備、グリッド中心+ー230km正方の中にある地震の月齢をスロット上で”1”加算します。 地震はM3.0以上を採用しM3.0未満は無視。 M3.0でもM6.0でも同じ”1”加算効果です。 対象となる地震数が96以上となった場合のみ測定に入ります。 結果、スロットのピーク値と最小値が出ます。 最小値はゼロの場合があり、これは”1”とします。 この比率が2.5以上の場合に色を付けます。 ピーク値が新月から満月の3スロット内にあれば東進(赤)、満月から新月の3スロット内にあれば西進(青)。  これを45kmピッチまたは9kmピッチで行います。 地震は一年以内が対象で、これを過ぎると(どんなに大きくても)忘れ去られます。 即ち、ひとつのマークはマーク中心+ー230km正方かつ一年以内に起きたM3.0以上の地震96件以上の6分割月齢加算結果比率に2.5以上のアンバランスがある事を赤or青で示しています。 これに深さ120kmより浅い地震という条件が加わります。

モデルに関しましては、現在模索中です。  潮汐力を考慮しないモデル [9] はそのままにしておきます(が、間違っています)。 

 

以下、潮汐力モデル中間報告その2です。(中間報告その1はこちら: [13]

潮汐力とは完全に対象でした。 [12] これは月の在る側と反対側で潮汐力は完全に同じ大きさ、という意味です。 地球と太陽も同じでした。

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

By Krishnavedala - Own work

これですと上弦月の場合と下弦月の場合で差が出て来ません。 潮汐加速を考慮しても同じ、月と地球の共通重心を考えても同じで、差が出てきません。 従って、潮汐力ではない別の非対称な要素を探す必要があり、現在模索中です。

 

測定は、N20°E120°からN47°E150°まで深さ120kmより浅い領域を45kmピッチ(概要)及び9kmピッチ(詳細)グリッド中心+ー230km正方にてスキャン測定し表示。 但し、地震ヶ所は緯度経度で入力されQGISも緯度経度表示なので変換に次ぐ変換です(実際にはかなり複雑です)。 日本を取り巻くプレートの動きは: [7] プログラムはRuby、地球儀表示はQGIS関数、日本地図は国土地理院さん提供shapeデータをQGIS表示、グラフはPython。 Ruby&QGISRuby&Pythonインターフェイスはすべてcsv。 プログラム概略仕様は: [8]

屋久島・種子島周辺は過去7300年間大地震を起こした事のない領域で、見やすくする為マークを外しています(外す処理をしているのはココだけです)。

以上です。

 

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

免責: 本予測は個人の推論によるもので、結果について何ら責任を負うものではありません。

引用: [1] 防災科学技術研究所 Hi-net 高感度地震観測網  [2] 気象庁|震源データ

[3] 地震の年表 (日本) - Wikipedia [4] List of solar cycles - Wikipedia

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

[6] 国立天文台 天文情報センター 暦計算室 [7] 気象庁 | 地震発生のしくみ

[8] プログラム仕様 (地震の評価マップ: 全体概要です) - 地震の予測マップ

[9] 2017-12-27 地震の予測マップ (北海道中央・東北太平洋側・関東・鳥取・岡山を除く日本全国が注意対象) - 地震の予測マップ

[10] Tidal acceleration - Wikipedia [11] Tidal locking - Wikipedia [12] Tidal force - Wikipedia

[13] 2017-12-31 地震の予測マップ (北海道中央・東北太平洋側・関東・鳥取・岡山を除く日本全国が注意対象) - 地震の予測マップ