[背景]
気象庁の有感地震全データを東経北緯を0.1度単位深さを１ｋｍ単位のグリッドに分けた場合、Ｍ７以上の震源、かつ段数２段目あるいは３段目以上のグリッド震源は１２震源しか存在せず、１２の特異震源と呼んでいる。その中の発生順で４番目の宮城県沖の震源#SP12-4について、その解析を行い、その束縛された自由度を下げた状態から地震の特性を見出そうとするものである。
[データ範囲と震源データ（同一グリッドおよび深さ違い）]
気象庁データ（1919/01-2026/01）における宮城県沖の全震源を含む矩形エリアから特異震源Ｍ７以上#SP12-4とグリッド座標中の同一経度緯度、深さ違いを抽出。#SP12-4-##1～##16の16組の震源データを得て、これらを対象に解析を行う。このエリア内に１２の特異震源#SP12-4の１段目・２段目が含まれている。
[震源推移]
下の図１に示すように##2～##9までは左手周り方向に推移、##9～##10で一旦、反転、##10～##15は再び左手周り方向に推移。##3と##4は１２の特異震源#SP12-4の１段目と２段目。
［震源間相関直線と幾何的中心]
図２として、図１の震源群に高相関係数の相関直線群Ｌ１～Ｌ７、および幾何的中心gc(all)、gc、特異震源##3/##4の中点を追加した。
各相関式と相関係数は以下の通り。
L1 : y = 1.0887x - 116.61 R² = 0.9998 (##1/##3/##12)
L2 : y = -0.3473x + 87.555 R² = 0.9968 (##2/##6/##7)
L3 : y = 1.9868x - 244.31 R² = 0.9967 (##3/##4/##10)
L4 : y = 0.8784x - 86.735 R² = 0.9983 (##4/#5/##7)
L5 : y = 1.2461x - 138.95 R² =1 (##6/##11/##16)
L6 :y = 2.0773x - 257.2 R² = 0.9998 (##7/##9/##14)
L7 :y = 0.5582x - 41.201 R² =0.9999 (##7/##12/##13)
これら７本の相関式の傾きと切片の相関式から得られるメタ相関式y = -142.18x + 38.177の傾きと切片から震源群の幾何的中心の座標 gc=（142.18, 38.177）は特異震源##3/##4を結ぶ線と##7/##8との交点付近、時系列に連続した２つの震源を結ぶ15本の線分群による幾何的中心をgc(all) = (142.21, 38.209)はＬ１とｌ６の交点付近にある。 ##3/##4の中点とGCを結ぶ線分の角度６４．２度と ##3/##4を結ぶ線分の角度６３．４度でほぼ同等。
Ｌ１～Ｌ７での震源群の東経・北緯データから追加で東経差ΔE・北緯差ΔN (例 L2-L1)、差分傾きΔN/ΔEのパラメータを算出。日時・深さ・Ｍを加えて、パラメーター間で１に近い以下の相関係数一覧は以下の通り。太字は相関係数0.999以上、斜字は３震源のうち、２震源が同一あるいは近接しているために疑似的に相関の高い可能性があり。
東経vs北緯および３点中２点が同一か近接の場合を除く高相関係数のグラフ(0.99以上)は以下の通り。ただし、相関係数があまり高くなくても相関傾きが互いに近い場合にグラフを併記の場合あり。
Ｌ１～Ｌ７を含む高相関組数は以下の通り。太字は相関係数が0.9998以上。
Ｌ１…０組（ただし、東経差ΔE・北緯差ΔN、差分傾きΔN/ΔEのパラメータの計算値無。）
Ｌ２…２組（深さ vs 東経差 ΔE）
Ｌ３…１組（東経 E vs 北緯差 ΔN）
Ｌ４…０組
Ｌ５…３組（日時 vs 東経 E、日時 vs 北緯 N、Ｍ vs 北緯差ΔＮ）
Ｌ６…５組（東経 E vs 北緯差 ΔN、北緯 N vs 北緯差 ΔN、東経 E vs 差分傾き ΔN/ΔE、北緯 N vs 差分傾き ΔN/ΔE、北緯差 ΔN vs 差分傾き ΔN/ΔE）
Ｌ７…２組（東経 E vs 北緯差 ΔN、北緯 N vs 北緯差 ΔN）
東経・北緯の空間的な幾何的中心のgcやgc(all)から最も離れているＬ５で日時との高い相関(R² =1 /0.9999)が見られる事は偶然ではない可能性あり。
No comments:
Post a Comment