Thursday, February 19, 2026

#SP12-12 日向灘の特異震源の解析(ミクロ的解析)

 [目的] 

12の特異震源の12番目、日向灘にある#12-12について、該当グリッド近傍の19震源を気象庁有感地震データベースから抽出、解析することにより幾何的特徴を見出し、その物理的意味を考察する。

[概要]

本解析では、同じく日向灘にある12の特異震源#12-2と同様に特異震源近傍の狭い空間領域に分布する複数の震源系列を対象とした。当該領域設定の下では、各震源系列を近似する回帰直線群において、傾きおよび切片の相関が1となり、その傾きおよび切片から導出される幾何的中心と特異震源の中点との位置関係、およびを回帰直線群の幾何的特徴有無を確認した。

更には特異震源#12-2との対比、およびマクロ的な日向灘におけるM7以上の分布から導きうる可能性について検討した。

[近傍震源リスト]

12の特異震源の1つ#SP12-12のグリッド(東経131.9度・北緯31.6度・震源深さ37km)と同一緯度・経度(0.1度単位)の微小範囲で深さ違いの18震源のデータリストは以下の通り。データが多すぎても少なすぎても解析が困難なので、グリッド要素センター値の緯度経度に対して±0.075度の範囲で抽出した。


オレンジ色は12の特異震源(東経131.7度北緯31.7度深さ31km)の2震源#12-13と#12-18で、後者は数2段目のM7.1でM7以上の12の特異震源の条件を満たす。黄色はいずれかの高相関線に属する震源。

[震源推移図]

上記リストの19の震源推移を下図に示す。

推移図各点のラベルはリストの連番に対応、オレンジ色の四角の2震源は特定震源#SP12-12のグリッドの1段目#12-13(M4)と2段目#12-18(M7.1)。その中点よりやや左側の赤いプロットが幾何的中心。

 [震源推移の相関群・非相関群]

下図に高相関ライン10本を赤線で追記した(R² = 0.9984以上)。近傍を赤線の通っていない非相関群は#12-1/#12-14/#12-17の3震源。

相関線群については右上の領域に集中している以外、特徴は見当たらず。各相関線の傾きと切片の相関図を以下に示す。相関式の負の傾きと切片が相関線群の幾何的中心となる。


高相関線群に含まれない非相関群の震源推移は#SP12-2同様に逆時計周り。一方で#SP12-2とは異なり、対角のマグネチュードの大きさはバラバラ(M4.8/M3.4)で、対角線の中点通過は3点のみなので不明。
 


[特異震源とその中点、および幾何的中心の位置関係]

下図に特異震源とその中点、および幾何的中心の位置関係を示す。


この座標系での傾きはm= 0.4571で角度にすると約24.6°となる。もう一つの日向灘の特異震源でのM7/M5.2中点と幾何的中心では90°(中点の真南)で異なるため、二つの特異点の役割は異なる可能性がある。

[他の12の特異震源群と日向灘M7以上の震源群との角度比較]
上記の幾何的中心と特異震源#SP12-12の線分方向(m= 0.4571/約24.6°)と比較するため、以下の角度を確認した。

 HG7-3(#SP12-2)HG7-4の線分方向
   傾きm= 1.4606/約55.6°
 HG7-5HG7-6(#SP12-12)の線分方向
   傾きm=1.4981/約56.3°
 #SP12-2 ↔ #SP12-12の線分方向
   傾きm=−0.5555/約ー29.1°
 #SP12-2↔#SP12-3の線分方向
   傾きm=0.4239/約23.0°
 #SP12-8~#SP12-10の線分方向(相関係数0.9999の震源群)
   傾きm= 0.4239/約23.0°
 #SP12-11↔#SP12-12の線分方向
   傾きm=0.6020/約31.4° 

以上から角度をグループ分けすると次のようになる。
 負の傾き群 約-29°(#SP12-2 ↔ #SP12-12)
 浅傾斜群 約23°(#SP12-2↔#SP12-3、#SP12-8~#SP12-10)、
      約25°(#SP12-12中点↔同幾何的中心)、
      約31°(#SP12-11↔#SP12-12)
 急傾斜群 約56°(HG7-3(#SP12-2)↔HG7-4、及びHG7-5↔HG7-6(#SP12-12))
 直角   90°(#SP12-2中点↔同幾何的中心)南北方向
ここで注目されるのは急な傾き群(マクロ)と負の傾き(超マクロ構造線)が成す角で、
      56°-(-29°)=85°
日向灘エリア(マクロ)でほぼ直角に近い角度となっている。
一方、#SP12-12近傍での主軸方向約25°は超マクロの浅傾斜群約23°とほぼ一致。

[結論]
日向灘M7以上のHG7系列の解析により、傾き約1.48(約56°)の急傾斜軸が安定して確認された。この方向は#SP12-2と#SP12-12を結ぶ(超)マクロ構造線と約85°の角度差を示し、ほぼ直交関係にある。一方、幾何的値中心軸(傾き0.4571, 約24.6°)は浅傾斜群(傾き0.4239, 約23.0°)と実質的に一致した。

Saturday, February 14, 2026

Analysis of Singularity Around Hyuga Nada(日向灘) - the 12 Singular Epicenter (#SP12-2) (micro-analysis)

  [Object] 

Especially for the one of the 12 Singular Epicenter  (#SP12-2), I tried the micro analysis in the narrow area around the #SP12-2. From all epicnters of the felt earthquake database of JMA (Japan Meteorological Agency) since 1919-Jan2026 , 18 epicenters are picked up as the common grids as the same grid latitude and longitude. By this analysis, the physical meaning will be lead from the geometric features.

[Conceptual framework] 

In this analysis, we examine multiple earthquake source sequences distributed within a narrowly defined spatial region in the vicinity of a singular source.

Under such a constrained regional setting, it is mathematically inevitable that, for the set of regression lines approximating each source sequence, the correlations between their slopes and intercepts become unity.

What is essential, however, is not this constraint itself, but the mode of intersection exhibited by the ensemble of regression lines under this condition.

In practice, the regression lines do not intersect at a single exact point. Rather, a small number of representative lines intersect in the vicinity of a central region, while the remaining lines are distributed around this region in an envelope-like manner.

In this study, we define the reference point formed near this central region as a fixed point, and interpret it as a central structure by which the earthquake source distribution is geometrically constrained.

This fixed point does not appear as a strictly defined mathematical point, but instead as a stable region with a finite spatial width. Such a feature is interpreted as reflecting a state in which the transmission of stress and strain is not entirely free, but instead the system as a whole is undergoing partial synchronization while maintaining a preferred directional tendency.

[The list of the 18 epicenters neighborhood of #SP12-2 in Hyuga Nada(日向灘) ]

18 epicenters are picked up as the common grids as the same grid latitude and longitude with #SP12-2 ( 31.6° N, 131.9° E, 37km depth). See the list as below,

These 18 epicenters have the common grid of N Lat/E Lon as the rounded values. #SP12-2 data are showed as orange colored. #2-1 is the 1st stage of M of #SP12-2 as the initial occurence of the grid)(The number of the stage of M (段数) is "1"). #2-3 is the 2nd stage of M of #SP12-2 as the bigger than the 1st stage in the same grid(段数 is "2") . #2-10 is the same stage of M of #SP12-2 as the less than the initial value of the 2nd stage(M5.2 <M7) in the same grid(段数 is "2". But, not the initial value). Initial values of each stage of M as 段数 are with red bold font.   

[Time-series plot of the 18 earthquake epicenters (N Lat, E Lon)]

Time-series plot of the 18 earthquake epicenters (N Lat, E Lon) is as below,


Each label value shows the Magnitude value. The start point was M5.2(#2-1) with orange colored scqure as the one of #SP12-2 in 1930. The route went with CCW until M3.9(#2-4).  The route went straight back the same way on M3.9(#2-4), and extened until M5.1 (#2-5). The route continued to M5.2(#2-10) as the 2nd step of #SP12-2. Since M5.2(#2-10) until M5.3(#2-16), the route went CCW. Then, it changed to CW  until the goal "M3.8(#2-18)" in 2024.



[Time-series plot of the 18 earthquake epicenters (N Lat, E Lon) "with correlation lines"]
Red correlation lines are added on the above graph as below,


There are 7 red lines as below,

3 red lines cross with #2-6(M3.9)<L1-L3>

l1:#2-5(M4.3)・#2-6(M3.9)・#2-7(M5.1)
    y = -0.2971x + 70.773 R² = 0.9999

l2:#2-6(M3.9)・#2-11(M5.2)・#2-16(M5.3)
    y = 0.2x + 5.2033 R² = 1

l3:#2-6(M3.9)・#2-12(M4.5)・#2-13(M3.7)
    y = 1.2342x - 131.21 R² = 0.9986

2 red lines cross #2-3(M7)<L4/L5>

l4:#2-3(M7)・#2-5(M4.3)・#2-9(M4.6)
    y = -1.5329x + 233.81 R² = 1

l5:#2-3(M7)・#2-15(M3.5)・#2-18(M3.8)
    y = 3.0822x - 374.86 R² = 0.9996

2 red lines cross #2-10(M5.2)<L6/L7>

l6:#2-8(M5.1)・#2-10(M4.3)・#2-12(M4.5)・#2-18(M3.8)

    y = 0.2765x - 4.9011 R² = 0.9995

l7:#2-9(M4.6)・#2-10(M5.2)・#2-17(M4.5)
    y = 31.5633333333333 R² = na

※This line isn't a correlation line because of horizontal one.

[Meta correlation between slopes and intercepts among L1-L6]

Meta correlation between slopes and intercepts among L1-L6 always has the correlation coefficient as the gragh below,


 It's a truism for any lines in narrow area. Even though the natural resilt, we can find out the center of the epicenter from both the slope and intercept of meta correlation formula as the geometric center among them.(It's at 131.89E, 31.596N).  In the time-series plot with correlation lines, the lines shows "the single geometric center" of them as the rotation cencer or convergence point of the lines. It suggests the center shouldn't only the result of the calculation, but the fixed point as "great tough and huge rock mass" (extra high intensity asperity) phisically. In fact, there has been no felt earthquake more than M4.2 within a radius of 4km since 1977.

[Location plots among 3 Singularity points/ "the single geometric center"]

Location plots among 3 Singularity points/ "the single geometric center" are shown as below,


The single geometric center among L1-L6 is located on 560m south of the midpoint between #2-3(M7) and #2-10(M5.2) . They are not located directly below the top of the triangle formed L4/L5. The structure doesn't consist of a single line, but also of some main lines and envelopes. This distance means the tilt between them in depth. #2-2(M4.7) is the geometric boundary condition of 0km, and controls the degree of freedom of this structure. The single geometric center acts not as "the Load oapplication point", but as the fixed point. It was led out as the center of correlationship of the structure. It doesn't suggest the instant break of the structure, but the metastable statement of it as the being fixed the transmission path of the stress. 

[L7 (Lat N = const.)  treatment]
 Because of enable of the difinition as a correlation for contant lines of Lat, there should be no discussion of the intrerpretation about L7 in this blog.

[Time-series plot of the epicenters except ones with high correlations on L1-L7 ]
There are only 4 epicenters with no colored in the 18 epicenter's list. These epicenters,  #2-1 (M5.2)/ #2-2(M4.7) /#2-4(M5.3)/#2-14 (M4.6), are located as below, 


There are two pairs of the silmilar M on the diagonal epicenters, and one line segment between M4.7/M4.6 goes through the midpont between M5.2/M5.3.  The other one goes through near the midpoint of the other one.

It suggest the bipolar destructions between each diagonal epicenter as the sililar magnitudes. It also suggest that the center of them should not be the destruction point, but the pivot point of the stess point.

It is interpreted that the appromaxinate location of midpoint of these diagonal lines should show the ballanced system spatially. It shows the stress transmission should be not random ones. It is interpreted phisically that the center should not accuate as the maximum destructive point, but also the "pivot" of the applied torques.

[Conclusion]
There have been some M7+ earthquakes in Hyuga Nada area (日向灘). The fact suggest this area around the structure around the Singularity should have the ability of the transmission/redistributuin of the similar huge stress.
Though the Singularity itself has the less possibility of a origin of the next M7+, it should have the greater possibility of becoming the stress bypass along the north-south direction  relatively because of the weaker struction in some future.   

Monday, February 9, 2026

#SP12-2 日向灘の特異震源の解析(ミクロ的解析)

 [目的] 

12の特異震源の2番目、日向灘にある#12-2について、経度・緯度が共通のグリッドで深さ違いの18震源を気象庁有感地震データベースから抽出、解析することにより幾何的特徴を見出し、その物理的意味を考察する。

[概要] 

本解析では、特異震源近傍の狭い空間領域に分布する複数の震源系列を対象とした。当該領域設定の下では、各震源系列を近似する回帰直線群において、傾きおよび切片の相関が1となること自体は数学的に必然である。

重要なのは、この制約条件のもとで、回帰直線群が示す交会様式である。

実際には、すべての直線が厳密に一点で交差するのではなく、少数の代表的ラインが中心近傍で交差し、その他のラインはその周囲に包絡的に分布する構造を示す。

本稿では、この中心近傍に形成される参照的な点を不動点と定義し、震源分布が幾何学的に拘束される中心構造として解釈する。

当該不動点は厳密な点解ではなく、微小な空間幅をもった安定領域として現れており、これは応力・ひずみの伝達方向が完全に自由ではなく、一定の方向性を保ったまま系全体で同期化しつつある状態を反映したものと考えられる。

[震源推移図]

12の特異震源の1つ#SP12-2のグリッド(東経131.9度・北緯31.6度・震源深さ37km)と同一緯度・経度(0.1度単位)の微小範囲で深さ違いの18震源のデータリストと推移図は以下の通り。推移図各点のラベルはマグネチュードで、オレンジ色の四角の3震源は#SP12-2のグリッドの1段目初期値(M5.2)・2段目初期値(M7)・2段目(M5.2)。画面中央上側付近の特異震源M5.2(#2-1)が始点、画面左下角のM3.8(#2-18)が終点と連番を付与。従って、ここでは特異震源#12-2の3点は#2-1/#2-3/#2-10(真の特異震源M7以上は#2-3のみ)。



 [震源推移群の相関群・非相関群]

下図にこの推移図で直線上にある震源群を赤線で図示した。


赤色の線群は全部で7本あり、一覧表で赤字で示した#2-6(M3.9)を通る線が3本、特異震源#2-3(M7)を通る2本、特異震源#2-10(M5.2)を通る2本でここではl1~L7と呼ぶ。各震源群と相関式・相関係数は以下の通り。

#2-6(M3.9)を通る3本 l1~L3

l1:#2-5(M4.3)・#2-6(M3.9)・#2-7(M5.1)
    y = -0.2971x + 70.773 R² = 0.9999

※相関直線の中間点#2-5(M4.3)から#2-6(M3.9)に推移、来た経路を折り返した延長線上に#2-7(M5.1)がある連続発生の3点。

l2:#2-6(M3.9)・#2-11(M5.2)・#2-16(M5.3)
    y = 0.2x + 5.2033 R² = 1

※特異震源#2-10(M5.2)から#2-11(M5.2)に進み、#2-6(M3.9)上を素通りして#2-16(M5.3)までは逆時計回りの推移上の2点#2-11(M5.2)・#2-16(M5.3)と、#2-6(M3.9)を結ぶ線。

l3:#2-6(M3.9)・#2-12(M4.5)・#2-13(M3.7)
    y = 1.2342x - 131.21 R² = 0.9986

※上記の特異震源#2-10(M5.2)~#2-16(M5.3)までの逆時計回りの推移上の2点#2-12(M4.5)・#2-13(M3.7)と、#2-6(M3.9)の近似直線。

#2-3(M7)通る2本 l4・L5

l4:#2-3(M7)#2-5(M4.3)・#2-9(M4.6)
    y = -1.5329x + 233.81 R² = 1

※特異震源#2-3(M7)からL1の中間点#2-5(M4.3)を通り、L5と共通の#2-9(M4.6)を通る線

l5:#2-3(M7)#2-15(M3.5)・#2-18(M3.8)
    y = 3.0822x - 374.86 R² = 0.9996

※特異震源#2-3(M7)から2点#2-15(M3.5)・#2-18(M3.8)を通る直線上。これら2点間#2-15~#2-18の経路は時計回りで#2-15は逆時計回りと時計回りの反転ポイントとなっている。

#2-10(M5.2)通る2本 l6~L7

l6:#2-8(M5.1)#2-10(M4.3)・#2-12(M4.5)・#2-18(M3.8)

    y = 0.2765x - 4.9011 R² = 0.9995

※#2-8(M5.1)から#2-10(M4.3)を通り、L3と共通の#2-12(M4.5)、L5と共通の#2-18(M3.8)を通る線。

l7:#2-9(M4.6)・#2-10(M5.2)・#2-17(M4.5)
    y = 31.5633333333333 R² = 定義不可

L4と共通の#2-9(M4.6)から特異震源#2-10(M5.2)まで推移線を辿り、その延長線上の#2-17(M4.5)を通る北緯一定の線。さらに延長した近傍にはl6の中間点がある。

l1~L6の傾きと切片の相関図

 l1~L6の傾きと切片は完全に相関係数1の直線上に乗る。狭い範囲での任意の直線で成り立つ関係であり、当然の結果だが、その相関式の傾きと切片から震源群の幾何的中心を見出すことができる。その座標は東経131.89度北緯31.596度であり、前述の推移図では#2-6(M3.9)付近のこの幾何的中心を唯一の「幾何学的な要(かなめ)」として、この点を中心に回転、あるいは収束するよう震源が推移しているように見えている。事実、2024年に付近でM7.1が発生したにも関わらず、幾何的中心の周囲 4km 以内で1977年以来 M4.2以上の地震が起きていない。これは、この地点が単なる計算上の中心ではなく、物理的に「極めて強固で巨大な岩塊(超高強度アスペリティ)」である不動点である事を示している。


特異震源3震源と幾何的中心の位置関係は下図の通り。特異震源3点(橙)のうちの#2-3(M7)・#2-10(M5.2)の中点からのl1~L6の幾何的中心との距離は真南に約560m。完全な垂直ではなく深部から浅部へと傾斜した震源配列が形成されている。この配列は
単一の直線ではなく、少数の代表ラインと包絡構造によって特徴づけられる。
#2-2 は深さ0km幾何学的境界条件として、この構造の自由度を制限している。不動点は荷重集中点ではなく、回帰構造の中心として導かれる参照点であり、これらの関係は即時的破壊を示すものではなく、応力伝達経路が固定化しつつある準安定状態を示唆する


L7(北緯一定)の震源ライン

…北緯一定(31.5633333333333度)のため、相関係数などの定量評価が定義できないため、
本稿では構造的解釈の対象からは除外する。

線群l1~l6に含まれない震源群(非相関群)の震源推移

…冒頭の一覧表で線群l1~l6の震源群(黄色)に含まれない無色の#2-1(M5.2)/#2-2(M4.7)/#2-4(M5.3)/#2-14(M4.6)4震源の時計回りの推移を下図に示す。対角での規模がほぼ同様、かつ#2-1(M5.2)/#2-4(M5.3)の対角線を結ぶ線分がもう一組の対角を結ぶ線分をちょうど2等分にしている。もう一組も2等分ではないが、それに近い比率を示している。

これらの事象からみると対角に配置された震源がほぼ同じ規模で発生する「双極破壊」が起きている。この事は中心点は破壊点ではなく、応力の支点(ピボット)である事を示唆している。

対角線が中間点を通る意味としては、系が空間的にバランスしている。このことは応力伝達が非ランダムであることを示唆している。物理的には中心点は最大破壊点ではなく、トルクがかかる“軸”として機能している。

[結論]

日向灘周辺では過去にM7規模の地震が実際に発生しており、この事実は当該領域が同程度のエネルギーを伝達・再配分し得る構造を有していることを示している。

特異震源自体が次のM7級発震点となる可能性は高くない一方で、本震源域においては、現状の相関構造を踏まえると、将来的にひずみのバイパスは構造的に弱い南北方向に生じる可能性が相対的に高いと考えられる。

Saturday, January 31, 2026

Analysis of M6.6+ in Hyuga Nada(日向灘) including both #SP12-2&#SP12-12 of (macro-analysis)

 [Summary] 

There are 2 special earthquakes of the 12 Singularities (#SP12-2&#SP12-12) in Hyuga Nada. For the macro analysis of them, I have syudied about M6.6+ in Hyuga Nada based on downloaded data (1919~Jan.2026) from JMA Seismic Intensity Database. The analysis showed not only  the counterclockwise loop of earthquakes similar with Case 1 (壱岐・隠岐) or Case 2 (the 12 Singularities), but also clockwise loop after the turning point. The study showed the possiility of the next M7+ as the closing point of the M7 loop around #SP12-2 on the exteded line from #SP12-12.  The next M7+ might become the trigger of Nankai Trough earthquake as M8+.

[macro-analysis] 
M7+ in Hyuga Nada area is shown as the circles inside of blue rectangular in the figure (downloaded from JMA) as below, 

The listed earthquakes of M7+ are listed with the data of the "date/time/seismic intensity/rounded longitude(E)/rounded latitude(N)/depth/ Number of Stage of M/M/elongitude(E)/rounded latitude(N)"  as the below list. The series numbers (#) are given to them as HG7-1 ~ HG7-6 (old~new). Orange colored one is M7 of #SP12-2 (HG7-3), and yellow colored one is M7.1 of #SP12-12(HG7-6). 


The 6 of M7+ in Hyuga Nada are plotted on the chart of Lng ve Lat of epicenters with the labels of each Magnitude as below,
M7 of #SP12-2 (HG7-3) is plotted as orange dot, and M7.1 of #SP12-12(HG7-6) is plotted as yellow dot. The loop of M7+ is counterclockwise(CCW) as similar with Case 1 (壱岐・隠岐) or Case 2 (the 12 Singularities),  but the loop hasn't closed. 

Characteristic of M7+ on the loop
The loop started with M7.1(HG7-1) as counterclockwise(CCW), and moved to M7.2 towards to NNE. Then, it went back with the sharp angle to M7 (#SP12-2) with CW. The next line from it to M7.5(HG7-4) formed the closed a smal closed triangle area. The red line (Southeast side) is almost parallel with the another red line (Northwest side) betweem M7.1 (HG7-5)/M7.1(#SP12-12/HG7-6). In fact, the SE side red line was drawn as the parallel line with the NW red line as the copy.   Each red line includes one of the 12 Singularities (#SP12-2/#SP12-12). 

The correlation of SE red line with HG7-1/HG7-3/HG7-4
 Please see the chart as below,
The correlation coefficient(R2) beteen the Lng(E) and Lat(N) is 0.9991 (very high value). Furthermore, there are higher R2 beteen Lng(E)/Lat(N) and M as 0.9999 /0.9996. It means the  size of these M7+ on the line are almost specified perfectly on the location of the Lng(E)/Lat(N).  




As shown with the descripions so far,  earthquakes occurrences of these should be not random ones, but  be with the predestined harmony.

Characteristic of M7+ on the loop
For the additional data, I picked up M6.6+ in Hyuga Nada (日向灘). 
The new list is as below, 
(The serial numers are updated as ##1~##13.)




The loop started with M7.1(##1) with blue lable of M6.9 as counterclockwise(CCW). The lines until M7.5(##5) with blue label are formed the small closed triangle area with sharp edges. The CCW loop continued until M6.7(##6) with yellow label. Then, the loop chnge the turning direction from CCW to CW. The new loop with CW has contined until to the end point of M6.6(##13).  

For the seeking the root cause of the turning on M6.7(##6), I added the correlation red lines with high R2 (0.9942~1) on the above figure as below,


3 red correlation lines are focuced on M6.7(yellow label) with red label. It means M6.7 is the focuced point of stress from different directions, and also the tough node common with differrent lines of big earthquake occurrences. So, I guessed the turning should have occurred on this tough point from CCW to CW.
Next, I would like to focus on the 2 red solid lines.
The red solid line with negative slope includes 4 points as correlation line with high R2 (Other lines have the high R2 among only 3 points). As explained, this line includes 2 points of the 12 Singularities (M7#SP12-2: orange dot and M7.1#SP12-12: yellow dot). Thease singularities are also the start point and the end point of the M7+loop in Hyuga Nada(日向灘). If the M7+ loop extend from M7.1#SP12-12 to(or over) M7#SP12-2, the huge closed loop will be formed, and will might mean the separation from the base plate. In result, a huge earthquake should occurr similar with "茂木ドーナツ(Mogi's Donut)" effect. 
The high R2 correlation line includng 2 singularities makes the possibility higher greatly. 


Another red solid line, with positive slope, has the perfect R2 among ##1/##5/##8 as 1. Similar with the high R2 correlation line among the M7+ plots, another high correlations there are. But, the high R2 are not between Lng/Lat vs M, but between Lng/Lat vs depth among  M7.1(HG7-1)/M7(HG7-3)/M7.5(HG7-4) in M7+ loop. 
It means Ms of the grids on the perfect correlation line shoud be fixred not only by the Lng(E)/Lat(N), but also the depth fatefully.




As dicussed above, there are multi high R2 correlation lines in Hyuga Nada. They restrict the epicenters position/M/the turning directions(CCW or CW), and seems to be never random ones.

Estimation of next M7+
In the case of the M7+ loop closing by the next M7+,  the estimated M might could be calcutated as below,

Mmax(N=7)=real max value (M7.5)+log10(chain reaction lines number:7)/1.5
               =7.5+log10(7)/10≒8.1

(Based on several hypothesis as belows,
  * The released energy might be the same as all charged including the closed area, and be the total of all high R2 correlation lines which had the same energy with M7.5 (real max value in Hyuga Nada). The idea is similar with "茂木ドーナツ(Mogi's Donut)" theory.  )  

Caution
If N=6/N=5, the estimated value becomes around M8.0 anyway. But,I coudn't fix the occurrence date. So, I'm not sure if it could be tomorrow or 100 years later, or so on.

#SP12-2&12 日向灘の震源の解析(マクロ的解析)

[概要] 
2つの特異震源を含む日向灘周辺のM6.6以上の震源推移がこれまでの壱岐・隠岐と北伊豆地震#SP12-1の解析同様の逆時計周りの推移と途中からの時計周りを示しており、日向灘での次のM7以上の震源が2つの特異震源を結ぶ線上で発生、南海トラフ地震のトリガーとなる可能性を示しています。

[マクロ的解析] 
12の特異点の震源のうち、#SP12-2及び#SP12-12の2つの震源を含む日向灘の震源域を内包する矩形エリアについて解析します。
下図は気象庁有感地震データベースでの日向灘を内包する矩形エリアでのM7以上の震源(1919~2026/1月)とそのリスト、震源推移図です。リストのオレンジの震源が#SP12-2の特異震源M7、黄色の震源が#SP12-12で一昨年2024年に発生した特異震源M7.1です。リストでは推移順を明確にするため、HG7-1~HG7-6の連番を付与しました。



M7以上の震源推移の特徴
震源推移図では壱岐・隠岐や#SP12-1 静岡県伊豆地方の震源の解析で見られた逆時計回りの推移がここでも確認されました。北東端のM7.5(HG7-4)への軌跡で最初に小さな三角形の領域を逆時計周りで閉じた後、今度は大きな四角形の領域を囲むように逆時計周りで推移しています。赤い線はM7.1(HG7-5)から黄色の特異震源M7.1(HG7-6)を通る北西側直線、及び、その南東側に平行線をもう一つのオレンジの特異震源(HG7-3)と北東端のM7.5(HG7-4)近傍を通るように引いたもので、南東側のM7.1(HG7-1)/M7(HG7-3)/M7.5(HG7-4)の近傍を通る直線がほぼ北西側の赤線と平行であることが分かります。
北東端M7.5を除き、規模がM7~M7.1と近い規模の地震が揃っているのも特徴の一つです。

南東側の赤の直線上に位置するM7.1(HG7-1)/M7(HG7-3)/M7.5(HG7-4)では東経vs北緯の相関が高い(0.9991)だけでなく、東経・北緯のそれぞれがMに対しての相関もそれ以上の線形性(0.9999/0.9996)があります。Mの時間や深さに対する依存は見られず、Mはほぼ完全に震源の東経・北緯位置に依存しています。




※ グラフタイトルに誤りがあり、投稿後に訂正しました。

ここまで示したように少なくとも一部の地震の発生は決して偶然ではなく、厳密なルールに基づき予定調和的に発生してると見る方が妥当と考えられます。

M6.6以上の震源推移の特徴
もう少し解析データを増やすためにM6.6以上の震源のリストとその震源推移図を以下に示します。




リストのオレンジと黄色の地震は12の特異震源ですが、震源を増やしたため、##1~13の連番を振りなおしています。そのため、ここでは##4と##12が特異震源となります。
震源推移図のラベルは各マグネチュードで、中央やや右下のM6.9の青いラベルが##1で##2はM7.1、右斜め上の##3のM7.2、折り返す形で##4のM7、再び急角度の折り返しで##5のM7.5、更に急角度の折り返しで##6のM6.7(黄色ラベル)。ここまでの推移はこれまで見られたように逆時計回りの推移が続いています。しかし、以降の推移は反転して時計回りが最後の##13のM6.6まで続いています。つまり、黄色ラベルのM6.7が反転ポイントとなっています。

この謎を解くためにこれらの震源の東経vs北緯の相関の高いラインを確認した結果を下図に示します。赤の実線及び点線で示した相関ラインは相関係数0.9942~1です。


逆時計周りの反転ポイントM6.7(黄色ラベル)では3本の相関ラインが交差しており、子の震源の場所が各方向からの応力の集中点であると同時に強固な結節点であることを示唆しており、そのために周回方向の反転がここで生じたものと推測されます。
一方、赤で示した負の傾きの赤実線は二つの特異震源を含む相関ラインは4点からなります(他のラインは3点の相関)。冒頭で示したように12の特異震源のうちの二つの震源M7(#SP12-2:オレンジ)とM7.1(#SP12-12:黄色)は日向灘のM7以上のループにおける始点(HG7-1)と終点(HG-6)であり、これらを結ぶM6.6以上の高相関ラインは次のM7以上の地震がM7ループを閉じる形で始点(HG7-1)付近で発生する可能性が高い事を示唆していると推測されます。


もう一つの正の傾きの赤実線は相関係数1で先ほどのM7震源群中の高相関線M7.1(HG7-1)/M7(HG7-3)/M7.5(HG7-4)で東経北緯vsM で別の高い相関があったのと同様に、こちらでは東経北緯vs深さにも高相関係数があります。
<訂正・補足>
M7震源群で一部に確認された高相関は東経 vs 北緯、及び東経/北緯 vs M、
M6.6震源群で一部に確認された高相関は東経 vs 北緯、及び東経/北緯 vs 深さです。
グラフタイトル見直ししています。確認不足で申し訳ありません。 

これはライン上の震源の規模Mが東経北緯だけでなく深さに対しても厳しく制約されている事を示します。





ここまで示したように震源位置の東経北緯には(高相関係数を持つ)震源ラインが複数存在し、そこで発生する地震の規模と震源推移を制約(逆時計周りor時計周り)しており、ランダムの発生とは言えないと考えます。

冒頭で示したM7の推移ループが始点のM7(12の特異点の2番目#SP12-2)付近で閉じた場合、「茂木ドーナツ(Mogi's Donut)」と同様にその囲まれた領域の蓄積エネルギーが放出される事、この領域での各グリッドの最大蓄積エネルギーをこの領域の最大M7.5である事、異なる方向の領域内の7本の震源ラインが各最大M7.5の蓄積エネルギーを持つ事、を仮定するとこれらの震源ラインが連動して全放出された場合の想定規模Mmax(N=7)は以下の試算となります。

Mmax(N=7)=現状の最大規模(実績値M7.5)+log10(連動する震源ライン数)/1.5
               =7.5+log10(7)/10≒8.1

なお、N=6/N=5の場合でもMmax(N=6/5)=8.0となり、M8クラスの可能性を示していますが、発生時期については来年か100年後か不明です。

Wednesday, January 28, 2026

Case2 /Analysis of Singularity Around Izu Region - the 12 Singular Epicenters (#SP12-1 to 12)

Interpretation of the 12 Singular Epicenters (#SP12-1 to 12)

Before explaining the screening method, we must define "Stages."

  • 1st Stage: The first earthquake in a grid (0.1° Lon/Lat, 1 km depth).

  • 2nd Stage: The first earthquake to exceed the magnitude of the 1st stage.

  • 3rd Stage: The first earthquake to exceed the magnitude of the 2nd stage.

In the JMA database, 94% of M7.0+ earthquakes are 1st stage. Cases of 2nd or 3rd stages are extremely rare (only 12 cases total). We define these 12 cases of M7.0+ occurring as 2nd or 3rd stages as the "12 Singular Epicenters."




Among these 12, two groups of three consecutive earthquakes (Blue: #SP12-4 to 6; Yellow: #SP12-8 to 10) lie almost perfectly on straight lines with correlation coefficients of 1.0 and 0.9999.



  • Blue Group: Located between Urakawa-oki and Offshore of Ibaraki Pref., moving nearly North-South on the landward side of the Japan Trench (1978–1982).


  • Yellow Group: Follows an interesting line including the 1995 Kobe Earthquake, the 2016 Kumamoto Earthquake, and the 2016 Fukushima M7.4.


The remaining six epicenters (Orange) form a counter-clockwise loop enclosing a triangular region, similar to the Iki-Oki(壱岐・隠岐) case. This loop started with the 1930 North Izu Earthquake (near a plate triple junction) and traveled along the Nankai Trough, Izu-Ogasawara Trench, and Japan Trench, eventually reaching the 2024 Hyuga Nada M7.1.


E-Lng vs M chart / N-Lat vs M : They shows counter-clockwise direction moving similarly.




The similarity in the singularities shown here suggests that earthquake occurrence is not random, but follows specific rules or groups. The following analysis will further deepen this conviction.

#SP12-12 日向灘の特異震源の解析(ミクロ的解析)

 [目的]  12の特異震源 の12番目、日向灘にある#12-12について、該当グリッド近傍の19震源を気象庁有感地震データベースから抽出、解析することにより幾何的特徴を見出し、その物理的意味を考察する。 [概要] 本解析では、同じく日向灘にある 12の特異震源#12-2 と同様...