[背景]
気象庁の有感地震全データを東経北緯を0.1度単位深さを１ｋｍ単位のグリッドに分けた場合、Ｍ７以上の震源、かつ段数２段目あるいは３段目以上のグリッド震源はこれまで１２震源しか存在しなかったが、先月２０２６年６月の岩手県沖でのＭ７．２が２段目で発生し、１３番目の特異震源となった。そのため、これまで１２の特異震源と呼称していたが、今後は段数の特異震源と呼称し、これまで通りに発生順に#SP-* (* は連番)の略号を付与する。今回の段数の特異震源は#SP-13となる。この震源を含むマクロ的な震源域である岩手県沖エリアでのＭ６．５以上の解析を行い、特異震源との関連性を見出そうとするものである。
[データ範囲と震源データ（Ｍ6.5以上）]
気象庁データ（1919/01-2026/06）における岩手県沖の全震源を含む矩形エリアからＭ6.5以上を抽出。#IW-1～#IW-17の１７震源データを得て、これらを対象に解析を行う。このエリア内には岩手県の段数の特異震源#SP-13 (#IW-17)だけでなく、三陸沖での#SP-7 の２段目Ｍ７．１(#SR-9/#IW-7)も含まれている。
[Ｍ６．５以上の震源推移]
下の図１に示すように#IW-1～#IW-5までは反時計周り、#IW-5で反転、#IW-4～#IW-7では時計周り、#IW-7で再反転、#IW-6～#IW-14まで再び反時計回り、#IW-7で反転、#IW-13～#IW-15まで時計回り、#IW15でほぼ反転折り返して延長線上の#IW-16へ推移。#IW-15～#IW-17までは時計周り。#IW-7が段数の特異震源の７番目の#SP-7、先日の段数の特異震源の１３番目#SP-13が#IW-17。なお、#IW-10～#IW-13は２０１１年の東日本大震災当日の余震群であり、#IW-14～#IW-16も２０１１～２０１２年に発生の余震群と見られる。今回の#IW-17は当エリアでの約１５年ぶりのＭ６．５以上の震源となっている。
[Ｍ6.5以上の震源間相関直線とその傾き、およぴ幾何的中心]前記１７震源を直線で結ぶ１２本の震源間の高相関式とその傾き・切片・（東経に対する）角度・相関係数は以下の通り。４震源を含む場合、より相関係数の高い組合せを採用した。
Ｌ４～６・Ｌ８はほぼ完全相関でR² = 1。他にもＬ３・Ｌ１１の相関係数も0.9998～0.9999で完全相関に近い。Ｌ１（27.8°）/Ｌ１０（26.8°）はほぼ平行。三陸沖のの段数の特異震源#IW-7を含む高相関直線はない。
これら１２本の相関式の傾きと切片の相関式から得られるメタ相関式y = -142.32x + 39.262の傾きと切片から震源群の幾何的中心の座標ＧＣが得られる。一方、時系列に連続した２つの震源を結ぶ１６本の線分群から得られる幾何的中心をGC(all)とすると各GCは以下の通り。
１２相関式から算出 GC （143.32, 39.262）
１６本の線分群から算出 GC(all) (142.57, 39.577)
最後の震源#IW-17を含む唯一の相関直線Ｌ１２を除く、１１相関式による幾何的中心をＧＣ0、最後の線分#IW-16/#IW-17を除く線分間の幾何的中心をＧＣ0(all)とすると各々の座標は以下の通り。
１１相関式から算出 ＧＣ0 （142.02, 39.119）
１６本の線分群から算出 ＧＣ0(all) (142.58, 39.815)
高相関直線群とこれらの幾何的中心の座標群も図２に追加した。
その結果、幾何的中心GCは相関直線Ｌ５上の#IW-12の真北付近、ＧＣ(all)は相関直線Ｌ５上の線分#IW-5/#IW6と線分#IW13/#IW-14が交差する付近に位置している。
直近の段数の特異震源#IW-17の影響を除いたＧＣ0、ＧＣ0(all)からはそれぞれ大きく推移しており、ＧＣ0→ＧＣは大きく南西方向に、ＧＣ0(all)→ＧＣ(all)はほぼ真北に推移。ＧＣ0(all)は#IW-6～#IW-10で囲われた北側の四角形エリアの南側の底辺（線分#IW-6/#IW-7）付近、ＧＣ(all)は#IW-4～#IW-7で囲われた三角形エリアの底辺付近に位置している。
[M6.5以上の震源間相関直線Ｌ１～Ｌ１２におけるパラメーターとの相関]
下表のようにパラメーターとして、一つ前のM7の地震との線分の傾き（北緯差ΔN/東経差ΔE）、線分切線分傾きΔN/ΔEを追加した。
Ｍ6.5以上の震源間相関直線Ｌ１～Ｌ１２における東経vs北緯以外のパラメーター間の相関高めの図を以下に抜粋した。
Ｌ１(#IW-1/#IW-2/#IW-8)…東経 vs 北緯の相関係数は0.9982。他の高相関はなし。
Ｌ２(#IW-1/#IW-3/#IW-4)…東経 vs 北緯の相関係数は0.9989。他の高相関はなし。
Ｌ３(#IW-1/#IW-12/#IW-15)…東経 vs 北緯の相関係数は0.9998。他に日時vsMでも相関係数１の高相関だが、#IW-12/#IW-15のMが同じため。
Ｌ４(#IW-1/#IW-10/#IW-15)…東経 vs 北緯の相関係数１の完全相関。その他に東経/北緯 vs 北緯差ΔNも高相関。
Ｌ５(#IW-1/#IW-8/#IW-11)…東経 vs 北緯の相関係数１の完全相関。他の高相関はなし。
Ｌ６(#IW-2/#IW-14/#IW-15)…東経 vs 北緯の相関係数１の完全相関。他の高相関はなし。
Ｌ７(#IW-3/#IW-6/#IW-14)…東経 vs 北緯の相関係数0.9979。その他、東経 vs 東経差ΔNでも高相関。
Ｌ８(#IW-4/#IW-5/#IW-14)…東経 vs 北緯の相関係数１の完全相関。その他、東経/北緯 vs Mでも高相関。
Ｌ９(#IW-4/#IW-9/#IW-13)…東経 vs 北緯の相関係数0.9999。#IW-17を入れた４震源の相関でも相関係数は0.9994（IW-9とIW-17が近接しているため）。３震源では東経/北緯 vs 東経差ΔEでも高相関。
Ｌ１０(#IW-5/#IW-10/#IW-12)…東経 vs 北緯の他の相関係数0.9972。#IW-10/#IW-12は北緯が近接しており、日時 vs 北緯の高相関はそれが原因の可能性もあり。北緯差ΔN vs 線分傾きΔN/ΔEでも高相関。
Ｌ１１(#IW-8/#IW-13/#IW-15)…東経 vs 北緯の他の相関係数0.9999。その他、日時 vs 線分傾きΔN/ΔEでも高相関。
Ｌ１２(#IW-10/#IW-16/#IW-17)…東経 vs 北緯の他の相関係数0.9992。その他、日時 vs Mでも高相関。
[震源間相関とその物理的挙動の推測]
第1期：#IW-1 〜 #IW-5 外郭形成
【期間】1928年 〜 1968年（約40年間） / 軌跡：反時計回り
個別震源の挙動
エリア北東の#IW-1（1928年、M7.0）に続き、南西部の#IW-2（1931年、M7.2）が発生。
その後、再び北東の#IW-3（1935年、M7.1）、#IW-4（1960年、M6.7）、そして西側の#IW-5（1968年、M7.2）へと、エリアの外郭をなぞるように反時計回りのループを形成。
駆動メカニズム
この時期は、未成熟で均一化されていない断層面に対し、プレートの沈み込み帯が「面全体の歪みキャパシティ（枠組み）」を構築していたフェーズ。 数理的には、北側バウンダリであるL1障壁（角度27.8°）の起点（#IW-1, #IW-2）がこの時点で強固にロックされ、エリア内のエネルギーが外側に逃がさない構造的な「外壁」を形成。
第2期：第1の反転と臨界化フェーズ（#IW-5 〜 #IW-7）
【期間】1968年 〜 1989年（約21年間） / 軌跡：時計回りへの第1反転
個別震源の挙動
1968年の#IW-5（M7.2、深さ20km）をピボット（回転軸）として、それまでの反時計回りから時計回りへと反転。
応力は深部へと急降下し、地殻深部72kmの#IW-6（1987年、M6.6）へと推移。
そのわずか2年後、今度は一転してプレート最浅部（深さ0km）での#IW-7（1989年、M7.1）の三陸沖の段数特異点#SP-7が発生。
駆動メカニズム
#IW-5での反転は、水平方向の応力分配が限界に達し、垂直方向（Z軸）へ応力が転進した状態。深部72km（#IW-6）から極浅部0km（#IW-7）への「垂直推移」は、地殻を上下に貫く断層セグメントの伝達パスの連結形成を反映。
特に#IW-7（#SP-7）という特異点の出現は、システム全体が単なるエネルギー蓄積状態から、連鎖破壊の臨界モードへの移行に対する物理的シグナル。
第3期：広域トリガーによる強制的整列フェーズ（#IW-7 〜 #IW-14）
【期間】1989年 〜 2011年（約22年間） / 軌跡：再反転と巨大連鎖（反時計回り）
個別震源の挙動
長いM6.5クラスの静穏期間後、2011年3月11日の東日本大震災（M9.0）の外部からのメガトリガーにより、システムは強制的に再反転。
地震当日のわずか数時間の間に、#IW-10（M6.5、29km）、#IW-11（M7.4、32km）、#IW-12（M6.6、17km）、#IW-13（M6.7、24km）がドミノ倒しのように連続発生。
翌日には#IW-14（M6.6、15km）がエリア南西部で発生。
駆動メカニズム
アスペリティの自然な力学的飽和を待たず、M9.0の超巨大な歪み変化（荷重)により、システムは強制的に「既存の幾何学的レールを高速で埋める」状態に移行。ランダムに推移したように見える震源群は、実際には相関係数1の完全線形のＬ４（#IW-1/#IW-10/#IW-15）およびＬ５（#IW-1/#IW-8/#IW-11）に位置した震源群を中心に推移。事前に設計されたようにＬ４/Ｌ５の「地殻の溝（ガイドレール）」に沿って位置している震源群は整列。
第4期：秩序再編と最新特異点フェーズ（#IW-15 〜 #IW-17）
【期間】2011年後半 〜 2026年（約15年間） / 軌跡：時計回りへの再編
個別震源の挙動
大震災直後の混乱期を経て、#IW-15（2011年11月、M6.7）、#IW-16（2012年、M6.6、52km）へと推移し、軌跡は再び時計回りへシフト。
そして、前回のM6.5以上から約15年間の長い沈黙（幾何学的エネルギー充填）を経て、2026年6月に深層部（44km）で最新の#IW-17（M7.2）の段数特異点#SP-13が発生。
駆動メカニズム
#IW-17の物理的役割は、大震災（第3期）が残した広域的な歪みの「不整合」を事後調整し、システムを定常状態へ戻す復元・再ロックモードとみられる。 数理的には、#IW-17は東日本大震災当日の#IW-10および2012年の#IW-16とを結ぶ縦方向のＬ１２（角度-85.4°、R²=0.9992）によって拘束。これにより、幾何学的中心はGC0から南西のGC（143.32, 39.262）へと大きく推移。北側の浅部ブロックに残留していた応力は44kmの深部の段数の特異震源へと「荷重転嫁」され、エリアを幾何学的な安定へと収束。
その他、Ｌ８(#IW-4/#IW-5/#IW-14)における東経/北緯 vs Mの相関係数0.9998やＬ９(#IW-4/#IW-9/#IW-13)の相関係数0.9999～1の縛りもあり、地震発生はランダムではなく、予定調和的に発生していることを伺わせる。
