[背景]
気象庁の有感地震全データを東経北緯を0.1度単位深さを１ｋｍ単位のグリッドに分けた場合、Ｍ７以上の震源、かつ段数２段目あるいは３段目以上のグリッド震源は１２震源しか存在せず、１２の特異震源と呼んでいる。その中の発生順で７番目の三陸沖Ｍ７．１の震源#SP12-7に関して、緯度経度は共通グリッドの深さ違いの１５震源について解析を行い、特異震源との関連性を見出そうとするものである。
[データ範囲と特異点近傍震源データ（#SP12-7との深さ違いグリッド）]
気象庁データ（1919/01-2026/05）から#SP12-7との深さ違いグリッドの１５震源のデータを得て、これらを対象に解析を行う。このエリア内に１２の特異震源#SP12-7 の１段目/２段目の初期値Ｍ４．８(##7)/Ｍ７．１(##11)も含まれている。なお、##11までは２０１１年の東日本大震災より以前、#12～#15が以降となっている。
[特異点近傍の震源推移]
下の図１に示すように##1～##9までは時計周り、##-9で反転、##9～##14では反時計周り、##14で変曲し、##15に至る。##7が特異震源#SP12-7の一段目、##11が二段目。
[特異点近傍の震源間相関直線とその傾き、およぴ幾何的中心]
１０本の震源間の高相関式とその傾き・切片・（東経に対する）角度・相関係数は以下の通りで図１の震源間推移に赤線で重ね書きしたものを図２に示す。Ｌ１はほぼ完全相関でR² = 1。Ｌ１/Ｌ５はほぼ平行、各々特異震源の一段目と二段目を含むＬ１０/Ｌ８はほぼ直交。なお、共に図示されたＧＣとＧＣ(all)については後述する。
これら10本の相関式の傾きと切片の相関式から得られるメタ相関式 y = -143.13x + 39.88の傾きと切片から震源群の幾何的中心の座標ＧＣが得られる。一方、時系列に連続した２つの震源を結ぶ１４本の線分群から得られる幾何的中心をGC(all)とすると各GCは以下の通り。それらの座標も図２で追加した
１０相関式から算出 GC （143.13, 39.88）
１４本の線分群から算出 GC(all) (143.13, 39.895)
ただし、##6/##7は同一東経であり傾きの分母が0となるため、ＧＣ(all)の計算から除外しており、実際は１３本からの算出。
その結果、幾何的中心GCはＬ１/Ｌ４の交点の北東付近、ＧＣ(all)はL２～Ｌ４およびＬ８に囲まれた領域に位置しており、最多５本の相関直線群が集中する##5に対して、それぞれ東側と南側に位置している。
[震源間相関直線Ｌ１～Ｌ１０におけるパラメーターとの相関]
下表のようにパラメーターとして、一つ前のM7の地震との線分の傾き（北緯差ΔN/東経差ΔE）、線分切線分傾きΔN/ΔEを追加した。
特異震源近傍の震源間相関直線Ｌ１～Ｌ１０におけるパラメーター間の相関高めの図（東経 vs 北緯 以外）を以下に抜粋した。
Ｌ１(##1/##4/##5)…東経 vs 北緯のみ高相関で、相関係数は１。
Ｌ２(##1/##3/##15)…東経 vs 北緯のみ高相関で、相関係数は0.9963。
Ｌ３(##2/##3/##5)…東経 vs 北緯の他、東経/北緯 vs 深さでも高相関。
Ｌ４(##2/##9/##10)…東経 vs 北緯の他、日時 vs 深さでも高相関。
Ｌ５(##2/##9/##12)…東経 vs 北緯の他、日時 vs 東経/北緯でも高相関。
Ｌ６(##3/##4/##6/##11)…東経 vs 北緯では４震源間の高相関（R² = 0.9999 ）。３震源ではその他の高相関もあり。
Ｌ６(##4/##6/##11)…東経 vs 北緯は完全相関(R² = 1)の他、日時 vs Ｍ、東経/北緯 vs 北緯差でも高相関。
Ｌ６(##3/##6/##11)…東経 vs 北緯 の高相関以外に、日時 vs Ｍ、東経/北緯 vs 北緯差でも高相関。
Ｌ６(##3/##4/##11)…東経 vs 北緯のみ高相関（R² = 0.9999 ）。
Ｌ６(##3/##4/##6)…東経 vs 北緯の高相関以外に、東経/北緯 vs Ｍでも高相関。
Ｌ７(##5/##7/##9)…東経 vs 北緯のみ高相関で完全相関（R² =1 ）。
Ｌ８(##5/##10/##11)…東経 vs 北緯の高相関以外に、東経/北緯 vs Ｍ、東経/北緯 vs 東経差ΔEでも高相関。
Ｌ９(##6/##12/##14)…東経 vs 北緯のみ高相関で完全相関（R² =0.9998 ）。
Ｌ１０(##7/##12/##13)…東経 vs 北緯のみ高相関で完全相関（R² =0.9980 ）。
[震源間相関とその物理的挙動の推測]
２震源間線分でもＬ１・Ｌ５に水平・垂直な直交系の下表の線分群があり、もう１つの直交系Ｌ８/Ｌ１０、４震源を含むＬ６が震源推移の主要骨格を成しており、その他の震源相関直線を省略したものを図３に示す。
■前半サイクル：時計回りのエネルギー分散と深度上昇（##01 〜 ##08）
【Step 1】##1 → ##2 → ##3 → ##4 （1960年3月〜1963年10月）
震源の推移: 最も深い領域（深さ30km・M5.7）の##1から開始し、##2（19km）→ ##3（6km）→ ##4（2km・M4.3）へと、時系列に沿って狭い範囲で浅部へと突き上げる急激な深度上昇（上昇リレー）を伴って、東経北緯平面の三角形の領域を囲う時計回りの移動 。
拘束メカニズム: ##1と##4は、最多の相関直線が集中するハブ的な震源##5とともに、完全相関ライン L1（相関係数 = １、角度76.3°）を形成 。 ##02と##03は、##05とともに L3（角度-19.7°）を形成し、このＬ３上では東経北緯平面だけでなく「東経/北緯 vs 深さ」でも相関係数 0.9978の高い相関で3次元的に線形拘束されており、Ｌ１とＬ３は96.0°と直角に近い角度で交差、とくにＬ３上の線分##2/##5(角度-14.0°)であり、Ｌ１との角度は90.3とほぼ直交。さらに線分##3/##4はその南西側の##6[Step2]/##11[Step3]へ相関係数0.9999で伸びるＬ６に拘束されている。
物理的解釈: 深部（30km）での極大値M5.7（##1）によって生じたストレスが、Ｌ１とＬ３の特定の物理的レールに沿って急速に上方移行し、最浅部（2km）の##4に達して一旦このエリアのストレスが開放された 。
【Step 2】##4 → ##5 → ##6 → ##7 → ##8 （1963年10月〜1970年9月）
震源の推移: 最浅部##04に達した後、東経北緯平面上で時計回りの軌道を維持しつつ南下、幾何学的中心GCの北側にあるハブ震源##5：17km、更に##06：20kmへと再沈降 。その後、##6から急激に最浅部の特異震源（0km・M4.8）の##7へと上昇、再び東側の##8（14km）へと沈降・推移 。[Step 1]の三角形を囲う時計回りの推移に対して、約３倍の面積の四角形の領域を囲う時計回りの推移。
拘束メカニズム: やや浅い(17km)##05は全体推移の軸的な役割を果たしており、Ｌ７（##5/##7/##9、角度-46.0°）はごく浅い(0km)特異震源##7を幾何学的に拘束、さらに深部領域##09（25km）をも拘束している（##9は[Step3]）。
物理的解釈: 浅部での狭い範囲での弱い部分での破壊が進んで表層でのストレス遷移が一段落後、ストレスは系の構造的中心のやや浅い(17km)のハブ震源（##5）に回帰、深部領域20kmの##6で浅部～深部のストレス蓄積が限界に達し、深部から一気に地表付近の開放端でもある特異震源１段目の##7（0km深さ）でストレスを開放（M4.8） 。##8への遷移は、[step2]の時計回りサイクルの始点##4近くの東端へ戻っていくプロセス 。
■ 中盤サイクル：反転・収束と特異震源２段目M7.1の発現（##8 〜 ##11）
【Step 3】##8 → ##9 → ##10 （1970年9月〜1974年2月）
震源の推移: ##8から南側のハブ#5(17km)近くの深部領域##9（25km）へ移行、##5～##9でほぼ四角の領域が閉じる「反時計回り」への軌道反転が発生 。その後、Ｌ４(##2/##9/##10)に沿って、更に深い##10（27km）へと推移 。
拘束メカニズム: Ｌ４（角度53.6°)は##2/##9/##10の各震源の東経北緯を拘束するだけでなく、深さも日時と共に線形で深くなる線形で拘束（相関係数0.9938）。
物理的解釈: 1960年以来続いたごく浅い～やや浅い震源深さでの時計回りのクローズドな応力巡回が限界に達し、系全体の幾何学的底面（25〜27kmの深部プレート境界）へストレスがが集中・蓄積したプロセス 。時間軸と深度が同期（Ｌ４）しながら反転・収束へ向かう##11のM7.1への準備期間 。
【Step 4】##10 → ##11 （1974年2月〜1989年11月(特異震源２段目M7.1）)
震源の推移: 約15年間の長い静穏期（エネルギー蓄積期間）を経て、系の南西端かつ最浅部（0km）において、本グリッド最大の特異震源 ##11で（M7.1、段数2）発生 。
拘束メカニズム: ##11は、最多４震源を高相関0.9999で結ぶ Ｌ６（##3/##4/##6/##11、角度40.3°）および Ｌ８（##5/##10/##11、角度23.4°）の交点。Ｌ６が東経北緯だけでなく、##4以降では「日時 vs M」が相関係数１の完全相関で拘束されており、##11の発生日時は予定調和的に決定されており、場所についてもＬ８との交点近くと１つ手前の##10発生時点でほぼ拘束されていた。
物理的解釈: [Step2]の##7 同様、約15年かけて時計回りで囲われた東部領域の浅部～深部に蓄積されたエネルギーが限界に達し、行き場を失ったその前の深部##10（27km）から一気に南東端の地表付近の特異震源２段目の##11（0km深さ）でストレスを開放（M7.1）。中心軸的な##05を支点として、東部の浅部～深部に蓄積されたエネルギーが南西端のアスペリティ（##11）へ一気に集中し、M7.1の巨大なアスペリティ破壊（特定震源グリッド２段目）に至った 。
■ 後半サイクル：M7.1後の反時計回りによる再調整と変曲（##11 〜 ##15）
【Step 5】##11 → ##12 → ##13 （1989年11月〜2026年4月）
震源の推移: ##11のM7.1後、2011年東日本大震災の超巨大応力変化を跨ぎ、約35年以上の時間をかけてＬ１０に沿って、##12（19km）、さらに北西の##13（16km）へ推移 。##8以降で##14まで反時計回りで推移。##15で変曲。
拘束メカニズム: ##12は、過去の反転起点である##2、##9とともに Ｌ５（角度76.4°）を形成 。このＬ５は、前半の急上昇ラインＬ１と「ほぼ平行」であり、図３に示したようにＬ１/Ｌ５と垂直な線分が多数あることからこの領域では構造的にこれらの直交直線群の格子状断層群が潜在的に存在、震源推移と共に主要な南北のストレス伝達の震源ラインがＬ１からＬ５に推移したと推測される 。Ｌ５上では「日時 vs 東経/北緯」が 高相関（相関係数0.9916）で時間とともに南東へ等速推移 。##13は、特異震源１段目（##7）、ポストM7.1震源（##12）とともに Ｌ１０（角度-64.5°）を形成し、これはＬ８（##10/##11を含む）と「ほぼ直交」
物理的解釈: ##11および東日本大震災によって系の応力バランスが大きく更新された後、大震災以降（##12〜）のストレスは、かつてのＬ１と平行な「新・並行ルート（Ｌ５）」を等速で移動しながら、特異点１段目の応力軸（Ｌ１０）を逆行・刺激する形で、反時計回りにエネルギーを再配分 。
【Step 6】##13 → ##14 → ##15 （2026年4月〜2026年5月：直近の急激な変動）
震源の推移: 2026年に入り、##13から2週間で西側の##14（17km・M5.4・段数2）へシフトし、さらに12日間で最西端の##15（16km・M4.9）へと変曲・横移動 。
拘束メカニズム: ##14は、##6、##12とともに Ｌ９（角度-41.7°）によって線形拘束（相関係数 0.9999）。最近の到達点##15は、初期の##1、##3とともにＬ２（角度11.0°）のライン上に拘束され、系の西端に至る。
物理的解釈: 反時計回りの過程で ##11のM7.1が発生、余震として#11を含むＬ８にほぼ垂直な##7を含むＬ１０上の#12/#13に再配分。#12以降では中心軸#5(17km)とほぼ同じ深さの19～16kmで推移、20km台を中心に再チャージに推移。M7.1の準安定化の過程で#14まで続いた反時計回りが中断し、#15に至る。
■ 結論（メカニズムの総括）
この15震源の推移メカニズムは、単なるランダムな群発地震ではなく、「Ｌ１/Ｌ５の平行断裂スリット」と「Ｌ８/Ｌ１０の直交応力軸」という、グリッド内にあらかじめ存在する幾何学的な構造骨組み（スキャフォールド）に沿って、応力が時計回り・反時計回りに流動する循環システムとみなされる 。
