[背景]
気象庁の有感地震全データを東経北緯を0.1度単位深さを１ｋｍ単位のグリッドに分けた場合、Ｍ７以上の震源、かつ段数２段目あるいは３段目以上のグリッド震源は１２震源しか存在せず、１２の特異震源と呼んでいる。その中の６番目に発生した茨城県沖Ｍ７．０の震源#SP12-6と同一グリッド、およぴ深さ違いの合わせての１７震源群について、震源推移の解析を行い、特異震源との関連性やマクロ解析との関連を見出そうとするものである。なお、１２の特異震源の中でも３段目は１番目の#SP12-1静岡県伊豆地方のＭ７．３と１１番目の#SP12-11福島県沖Ｍ７．４、そしてこの#SP12-6のＭ７．０の３震源のみである。
[データ範囲と震源データ]
気象庁データ（1919/01-2026/01）における浦河沖の全震源を含む矩形エリアから同一グリッドおよぴ深さ違いを抽出。##1～##17の１７震源データを得て、これらを対象に解析を行う。このエリア内に１２の特異震源SP12-6の１～３段目初期値（##6/##8/##9）および３段目(##10)が含まれている。最大値は##9のＭ７．０。
[特異震源近傍の震源推移]
変曲の多いマクロでの震源推移に対し、下の図１に示すように##1～##7、##7～##9、##9～##17で反時計周り（##7、##9は変曲点）。
[特異震源近傍での震源間相関直線と幾何的中心]
震源間の高相関式と相関係数は以下の通り。Ｌ１２は東経一定の震源ライン。
Ｌ１(##1/##7/##17) y = -6.7432x + 993.6 R² = 0.9973
Ｌ２(##1/##8/##14) y = -4.25x + 639.61 R² = 0.9954
Ｌ３(##2/##7/##9) y = -0.5258x + 110.83 R² = 0.9991
Ｌ４(##2/##14/##16) y = -2.0756x + 330.91 R² = 1
Ｌ５(##4/##6/##12) y = 0.4859x - 32.791 R² = 1
Ｌ６(##6/##11/##13) y = 0.2524x + 0.3634 R² = 0.9969
Ｌ７(##4/##10/##14) y = -0.4798x + 104.35 R² = 0.9974
Ｌ８(##8/##11/##15) y = -1x + 178.17 R² = 1
Ｌ９(##5/##11##15) y = -0.7368x + 140.81 R² = 0.9983
Ｌ１０(##9/##14/##17) y = 2.8408x - 367.06 R² = 1
Ｌ１１(##12/##15/##16) y = 1.2015x - 134.39 R² = 0.9971
Ｌ１２(##3/##5/##6) x=142.033333333333
これらＬ１～Ｌ１１の１１本の相関式の傾きと切片の相関式から得られるメタ相関式の傾きと切片から震源群の幾何的中心の座標ＧＣが得られる。一方、時系列に連続した２つの震源を結ぶ１６本の線分群から得られる幾何的中心をGC(all)とすると各GCは以下の通り。それらの座標も図２で追加した。ただし、GC(all)の算出に当たり、前後の震源間の東経差が０の場合は傾きが算出不能のため、２組の線分データは算出から除外している。
１１相関式から算出 GC (141.98, 36.211)
１６本の線分群から算出 GC(all) （142, 36.191）
その結果、幾何的中心GCは##1の南東近く、GC(all)はＬ４と線分##6/##7の交点近くの北東付近に位置している。
震源間相関直線が最も集中しているのは##14で、次いで##6・##8・##15の３点。これらのうちの##8と##15を含む高相関直線Ｌ８(##8/##11/##15)は相相関係数 R² =1 で完全相関の y = -1x + 178.17で傾きがちょうど－１（傾き４５°）となっており、１２の特異震源の２段目Ｍ４．９（##8）を含んでいる。
[震源間相関直線Ｌ１～Ｌ１１におけるパラメーターとの相関]
下表のようにパラメーターとして、一つ前のM7の地震との線分の傾き（北緯差ΔN/東経差ΔE）、線分切線分傾きΔN/ΔEを追加した。
震源間相関直線Ｌ１～Ｌ１１におけるパラメーター間の相関高めの相関図（東経 vs 北緯 以外）を以下に抜粋した。
Ｌ１(##1/##7/##17)…東経 vs 北緯の他、日時 vs 深さ、東経/北緯 vs 北緯差ΔNでも高相関
Ｌ２(##1/##8/##14)…東経 vs 北緯の他、東経/北緯 vs M、線分傾きΔN/ΔEでも高相関
Ｌ３(##2/##7/##9)…東経 vs 北緯のみ高相関（相関係数0.9991）
Ｌ４(##2/##14/##16)…東経 vs 北緯の相関係数１
Ｌ５(##4/##6/##12)…東経 vs 北緯の相関係数１
Ｌ６(##6/##11/##13)…東経 vs 北緯の他、日時 vs 線分傾きΔN/ΔEでも高相関
Ｌ７(##4/##10/##14)…東経 vs 北緯の他、東経/北緯 vs 深さでも高相関
Ｌ８(##8/##11/##15)…東経 vs 北緯の相関係数１で傾きが－１ちょうど。
Ｌ９(##5/##11/##15)…東経 vs 北緯のみ高相関（相関係数0.9991）
Ｌ１０(##9/##14/##17)…東経 vs 北緯の相関係数１の他、東経/北緯 vs 北緯差 ΔN、東経差ΔE vs 線分傾きΔN/ΔEでも高相関
Ｌ１１(##12/##15/##16)…東経 vs 北緯の他、日時 vs 東経/北緯、東経差 ΔE vs 北緯差 ΔNでも高相関
Ｌ１２(##3/##5/##6)…東経 一定の他、日時 vs 深さでも高相関
[震源間相関直線群と傾き]
震源間相関直線Ｌ１～Ｌ１１、およびＬ１２(東経一定)の東経に対する傾きは以下の通り。
Ｌ１(##1/##7/##17) y = -6.7432x + 993.6 -81.6°
Ｌ２(##1/##8/##14) y = -4.25x + 639.61 -76.8°
Ｌ３(##2/##7/##9) y = -0.5258x + 110.83 -27.7°
Ｌ４(##2/##14/##16) y = -2.0756x + 330.91 -64.3° （Ｌ５と直交）
Ｌ５(##4/##6/##12) y = 0.4859x - 32.791 25.9° （Ｌ４と直交）
Ｌ６(##6/##11/##13) y = 0.2524x + 0.3634 14.2°
Ｌ７(##4/##10/##14) y = -0.4798x + 104.35 -25.6°
Ｌ８(##8/##11/##15) y = -1x + 178.17 -45.0°
Ｌ９(##5/##11##15) y = -0.7368x + 140.81 -36.4°
Ｌ１０(##9/##14/##17) y = 2.8408x - 367.06 70.6°
Ｌ１１(##12/##15/##16) y = 1.2015x - 134.39 50.2°
Ｌ１２(##3/##5/##6) x=142.033333333333 90.0°
下の図２にＬ１～Ｌ１１、およびＬ１２(東経一定)、ＧＣ、ＧＣ(all)を追加。
１７の震源のうち、２つの震源の全組合せの線分群が成す角を確認した結果、線分の成す角が０°/±９０°/１８０°の線分群と±４５°/１３５°の線分群、さらにＬ１～Ｌ１２の中で唯一、直交関係にあるＬ４/Ｌ５を追加したものを図３に示す。
Ｌ４/Ｌ５/Ｌ８はいずれも相関係数が１となっている。さらに線分の成す角が０°/±９０°/１８０°となっている線分群に含まれない##4・##11・##14は±それぞれ##4がＬ５に、##11が４５°/１３５°の線分群に、#14がＬ４に含まれている。前述のようにＬ４/Ｌ５はＬ１～Ｌ１２の中で唯一、互いに直交関係を有する組合せとなっている。
[震源間相関とメカニズム推定]
図３のように多くの２組の震源が同経度、同緯度に位置しており、格子状構造を形成している。
フェーズI：マクロ応力の流入と局所フレームへの完全固着（1941年〜1979年）
マクロ領域でM7.3を含む東北東-西南西主せん断面（マクロL6）が固着状態にあるため、広域歪みが特異震源近傍のミクログリッドへ流入し、最適な破壊経路を「探る」フェーズ。
##1 (24km, M4.8) → ##2 (23km, M5.0)
推移: 1941年の約1時間で推移。
メカニズム: 局所アスペリティの浅発層（周縁部）にマクロ応力が到達し、初期の微小な亀裂が発生。
##2 → ##3 (79km, M6.1)
推移: 約20年後の1961年、深さ79kmの深部方向へ推移、規模もM6.1に拡大。
メカニズム: 浅部で受け止めた応力が、プレート境界面に沿って深部（スラブ内）への亀裂を形成。
##3 → ##4 (29km, M5.0)
推移: わずか10時間後、再び29kmの浅部へ折り返し（移動角度104°）。
メカニズム: 深部に形成された亀裂を介して元の深さ方向に折り返し。この急激な上下運動により、特定震源グリッド近傍のミクロ領域にマクロからのストレスが蓄積。
##4 → ##5 (40km, M4.5)
推移: 約18年間の沈黙の後、1972年に東端の壁（ミクロL12）と南端ラインの直交交点に##5が発生。
メカニズム: この18年間、マクロから流入したストレスが構造的な壁としてのＬ１２の、特に南端コーナーに蓄積。長期の蓄積の結果、局所の「テンション・ボルト」としての張力蓄積が限界に到達。
フェーズII：ねじれ応力の極限化と最大剪断の発生（1982年7月）
マクロ構造からの押し込みにより、局所グリッド内部が反時計回りの締め付けが続いた結果、ボルトの外郭から崩壊のフェーズに移行。
##5 → ##6 (30km, M4.7)
推移: 東端の壁（Ｌ１２）に沿って真北（0°方向）へ移動。
メカニズム: 壁に沿って応力が滑動し、北端の壁に到達。ここでストレスの逃げ場を失い、特異震源の一段目M4.7の破壊が発生。
##6 → ##7 (40km, M5.2)
推移: わずか約5時間後、135°の方向へ深さを増して突進。
メカニズム: 東と北の壁にロックされたねじれ応力が限界を超え、最も滑りやすい最大剪断応力方向（45°/135°系）へと内部破断が発生（Ｌ８の形成）。
##7 → ##8 (30km, M4.9)
推移: さらに約10時間後、真北（90°方向）へ変曲し、再び深さ30kmへ。
メカニズム: 剪断方向へ逃げた応力が、共役関係にある0°/90°系の制御を再度受け、特異震源の2段目M4.9の破壊が発生。3段目M7.0へ向けた軌道修正の段階に至る。
フェーズIII：テンション・ボルトの完全破断（1982年7月23日）
局所の限界破断であり、同時に「マクロ構造の安全弁の消失」を意味する。
##8 → ##9 (30km, M7.0)
推移: ##8の約18時間後、180°方向（真西）へ一気にスライドし、特異震源3段目となるM7.0が炸裂。
メカニズム: 局所アスペリティの最も強固なコアの破壊。この局所空間の最大歪みの解放により、マクロ構造のＬ６（マクロでの主せん断面）への荷重転移を開始。
フェーズIV：局所の残骸処理と準安定化（1982年〜2008年）
ピンが吹き飛んだ後の局所グリッドが、近傍の微細クラックを介して応力開放が進行するフェーズ。
##9 → ##10 (30km, M5.1)
推移: M7のわずか40分後、北端の壁（##6/##10）上へ27°で移動。
メカニズム: 既存の強固なフレームワーク（北端の壁）への即座な応力開放による、安定化への最初のプロセス。
##10 → ##11 → ##12 → ##13 → ##14 → ##15
推移: 反時計回り推移を維持しながら、深さを40km→47km→53km→42km→69kmと深部へシフトさせつつ、Ｌ４/Ｌ５(局所共役断層）との交点などを経由して推移。
メカニズム: マクロの主応力から解放された局所グリッドが、自らの内部に形成した「ローカルなダンパー（L4/L5の十字構造）」を用いて、残存した回転トルクを時間をかけての事後処理・吸収のプロセス。2008年の##15（69km）をもって、局所の準安定状態は一段落 。
フェーズV：巨大外部トリガーによるマクロ崩壊とミクロへの逆流入（2011年〜）
2011年3月11日のM9.0によってマクロのＬ６（マクロでの主せん断面）が臨界を超え、IB13でM7.6誘発、その破壊エネルギーが再びミクログリッドに流入のフェーズ。
##15 → ##16 (57km, M5.2)
推移: 東日本大震災の9日後（2011/3/20）に発生。Ｌ４およびＬ１１のライン上に位置。
メカニズム: 1982年の系列の単なる続きではなく、マクロ領域での最大値IB13（M7.6）の巨大破壊によって茨城県沖全体の応力場が激変したことによる「逆フィードバック」。強烈な外部衝撃により、かつてのミクロな傷跡（Ｌ４/Ｌ１１などの交点）に再度、滑りが発生。
##16 → ##17 (48km, M5.2)
推移: 約9ヶ月後（2011/12/5）に発生。Ｌ１およびＬ１０の交点に位置。
メカニズム: マクロからの逆流入応力が、ミクロ空間に残されていた別の幾何学的ネットワーク（Ｌ１/Ｌ１０）を最終的に刺激し、反時計回りの推移に沿って一連の局所的・広域的な応力の応答が収束方向。
[メカニズム推定まとめ（マクロ一部含む）]
① マクロからミクロへのエネルギー流入と「反時計回り」のねじれ
太平洋プレートが約24度の斜面（スラブ上面）に沿って西北西に沈み込むマクロな圧縮応力（最大主応力軸）が、領域に定常的に流入 。
特異点「#SP12-6」のミクロ領域（##1〜##17）では、このマクロ応力を受け止めた際、局所的な構造（東端の壁となるＬ１２など）に歪みがトラップされ、空間内で「反時計回りの局所的なねじれ（トルク）」が発生 。このねじれのプロセスが、##1〜##7などの反時計回りの一方方向の回転推移として反映。
② 蓄積期から臨界期への移行（直交フレームのバトンタッチ）
蓄積期（0°/90°系の作動）:
応力が徐々に溜まる時期は、地質境界に沿った構造的バウンダリ（エネルギー障壁）が優位に機能 。
極限ねじれ期（45°/135°系の作動）:
ねじれが限界に達すると、最後まで歪の残った「最大剪断応力方向」へと破壊がシフト。Ｌ８（傾き-1、-45.0°ちょうど）上に現れた特異震源2段目（##8、M4.9）は、この局所的な限界突破のサイン（前震）。
③ 特異震源の破断（安全弁の作動）とマクロへの荷重転移
1982年、特異震源近傍のミクロ領域の限界到達により、特異点グリッド内で##9（＝マクロ解析におけるIB11、M7.0）が発生 。
この特異震源の破断は、局所の「ストッパー」が外れたことを意味します。エネルギーは消滅したわけではなく、この特異点を端点とするマクロの相関ライン（Ｌ１０）を介して、周辺のマクロ領域へ徐々に荷重転移（応力再分配）を起こし、将来の最大規模地震の発生場所（IB13周辺の結節点・ハブ）の絞り込み（ロック）が進行 。
④ 巨大外部トリガーによるマクロ破断とミクロへの逆流入（散逸フェーズ）
2011年3月11日、系外の巨大地震（M9.0）による動的応力がシステムを直撃し、限界まで歪みを溜めていたマクロの主せん断面（Ｌ６）が滑ることで、本エリア最大規模のIB13（M7.6）が誘発 。その後のミクロでの発生はそれまでの反時計回りの推移
No comments:
Post a Comment