I regarded the record update of the specific depth per 10km of the specific latitude coordinate as the state changed. In fact, huge earthquakes like M9 bring out many aftershocks suddenly.
But, it seems each grid of (i, j,k)=(unit longitude,unit latitude, 10km) may be the pie of domino as the energy capacitor. So far, I have treated earthquake data as the combination of latitude coordinate and depth(per 10km). And, the max of magnitude per place is regarded as the max limitation reflecting the information of the geological structure.
For the domino analysis for some areas having only shallow epicenters, I must have changed the mesh of grids more detailed. But, it means the analysis must become more difficult and more complex. I finally have realized I should return to trace the maximum magnitude change trend per 1km in chronological order as data observation.
The result of it for 三陸沖 (Sanriku Oki) including the epicenter of M9.0 Great East Japan Earthquake in 2011 is as below,
Sorry for the tiny letters. This table includes all magnitude data of 0-50km depths of the 三陸沖 (Sanriku Oki) area between 1919-2025/12/12. Rows are dates of occurrence date, and columns are epicenter's depth per 1 km. Values of the table are the max magnitude at that point. Blank means no data or less than M6.0. Right side columns are the max values of 1km-50km or 0-50km.
This table shows that max value updates had been done in all depths between 1km-24km(the epicenter depth of M9) except 23km until M9 occurrence day of 2011/3/11. The blank cells of the upper section show it.
This table is an excerpt, with 1-50km columns and the bottom part omitted. The last 4 columns consist of 2 parts of "including M9"/"Except M9". Each column part consists of 2 columns of "Max for all"/"Max except 0km". The reason why M9 was omitted in one part is that its inclusion would make the max of all subsequent occurrence dates M9 (Refer to the bottom section of the table below). Magnitude values are shown by using a color scale. The color scale ranges from light green (M6.3) to red (M9).
The comparison between the last adjacent columns indicates that magnitude values of 0km are different from others. It may be one of the factors of outliers of G-R law. The difference between "0km" and others should be caused by "0km" being open end without a deeper layer. To emphasize it, magnitudes of 0km are shown only on the occurrence date in "0km" column.
Below is the bottom section of the table. After M9 (red) in 2011/3/11, max for all data became all M9.0 jumping up from M7.3. Max except M9 became M7.3 increased from M7.2. It is another factor of outliers of G-R law. As G-R law is valid strictly in micro areas, I believe it should be also valid in macro areas.
The reason why the law seemingly fails should be that charged energy as stress was released in M9 like a domino. The latest value of the max magnitude per 1 km should become the risk index of the next earthquake's magnitude/occurrence timing.
Below is the middle section of the table focused on before and after M9. Yellow colored cells show the 2nd longest silent term of 0km (The 1st longest one has been continuing since 2011/4/14 as shown in the table above). It indicates the 2nd one may be one of the significant signs of M9 as the charging term of stress energy.
As I wrote above, max value updates had been done in all depths between 1km-24km except 23km until M9. Especially, 8km/9km/12km/29km/36km increased the day before or 2days ago of M9 occurrence as indicated by thick frames. M6.2 of 12km was an aftershock of M7.3 in 2011/3/9. M7.3 of 8km jumped up from M6.8 as shown by the sudden color change 2days ago of M9. Others in 2011/3/10 seems they were also the aftershocks of M7.3. The last 3 updates of 1km-24km except 23km occurred between the last 2 days.
It seems the rapid updates around the shallow depth of 10km may be valid as the sign of the huge earthquake. Especially, the value M6.8 of 9km had no update as M6.2 since 1929. It may be likened to the last column of 9km falling down to support East Japan just before M9.
