図2aの地震の震源分布は、1983年1月6日から2012年1月30日までに発生したマグニチュードが3.0以上の地震で、それらは気象庁の全国震源カタログと一元化震源カタログから集めたものである。水色の点で示された各地震の震央は、マグニチュード(MAG)が3.0以上3.5未満で、黒色の丸印の震央は、MAGが3.5以上の地震である。MAGが5以上で6未満の地震には、橙色の三角印、MAGが6以上の8個の地震の地震は、薄橙色で塗りつぶした丸印と発生年代のラベルを付けてある。同じ年代に発生した地震には、アルファベットのaとbとで区別している。1983年に発生した三朝地震は、1983bとラベルが付けられている。2つの群発地震は、発生年代の前に群発(Swarm)のSが付けられている。深発地震に関しては、MAGが3.0以上3.5未満の地震は黒色の円内を、MAGが5以上で6未満の地震は、黒色の三角内を、MAGが6以上の地震は黒色の大きな円内を、それぞれ青色で塗りつぶして表示してある。直前にピークを迎える地震発生(歪エネルギー密度)のサイクルである。左側の縦軸目盛は、NCI(m,70)、NCD(m,70)グラフの目盛である。右側の縦軸の6以上の目盛は、MAGの時系列[d(MAG,m)]をMAGとラベル付けした矢印の高さのグラフの目盛である。従って、MAGのグラフは、マグニチュードである矢印の高さが6以上の大地震になると表示される。右側の縦軸目盛にLONとラベル付けされている範囲は、震源パラメータLONの時系列d(LON,m)を点グラフにしたLONの目盛で、経度が131.5度〜136.5度に相当する。大地震の余震や群発地震は、震源が略同じなので、それら余震や群発地震の震源の経度は殆んど変化しない。従って、d(LON,m)の点グラフの変動は無くなる。図2a−図2cに使用されたマグニチュードが6以上の大地震の発生年代ラベルは、そのまま図4の時系列d(MAG,m)に表記されている。2つの群発地震S1984とS2001は、時系列d(LON,m)の点グラフに矢印と共にラベルが付けられている。図2aの地震の震源分布は、1983年1月6日から2012年1月30日までに発生したマグニチュードが3.0以上の地震で、それらは気象庁の全国震源カタログと一元化震源カタログから集めたものである。水色の点で示された各地震の震央は、マグニチュード(MAG)が3.0以上3.5未満で、黒色の丸印の震央は、MAGが3.5以上の地震である。MAGが5以上で6未満の地震には、橙色の三角印、MAGが6以上の8個の地震の地震は、薄橙色で塗りつぶした丸印と発生年代のラベルを付けてある。同じ年代に発生した地震には、アルファベットのaとbとで区別している。1983年に発生した三朝地震は、1983bとラベルが付けられている。2つの群発地震は、発生年代の前に群発(Swarm)のSが付けられている。深発地震に関しては、MAGが3.0以上3.5未満の地震は黒色の円内を、MAGが5以上で6未満の地震は、黒色の三角内を、MAGが6以上の地震は黒色の大きな円内を、それぞれ青色で塗りつぶして表示してある。カタログと気象庁カタログから、図2の中国・四国・近畿を含む緯度32度-36度、経度131.5度-136.5度の小さな領域内に発生した地震で、深発地震を除き、マグニチュード(M)が、M≧Mc(Mc=3.5)で選択される地震のみ抽出する。抽出した地震の震源要素を、地震の発生した順番に並べ、次の様に、時系列データを作成する。最初の地震(0番目)から2番目まで、次に与えられているとする。これら3つの地震の各震源要素に震源時の差である間隔時間INTを加え、震源要素を(LAT、LON、DEP、INT、MAG)とした時系列データの作成は、それぞれの単位を、LATとLONが度、DEPがkm、INTが時間(60分)とすると、次の様になる。0番目=(LAT=34.453、LON=132.607、DEP=27.0、INT=X、MAG=3.5)、1番目=(LAT=33.817、LON=131.843、DEP=83.0、INT=49.0186、MAG=4.0)、2番目=(LAT=33.468、LON=132.448、DEP=48.0、INT=217.7151、MAG=4.1)0番目のINT=Xは、その前の地震情報がないので震源時(発生時刻)の差である間隔時間INTが算出できない事を意味する。従って、1番目が時系列データの始まりとなる。この時系列の表記方法を次に述べる。時系列にする各震源要素を任意パラメータcで表記し、時系列の時間tを地震の発生順序j(j=1、2、3、・・、m、・・)とし、各震源要素の時系列を、[c]=[c(1)、c(2)、c(3)、・・、c(m)、・・]とする。各時系列[c]は、震源要素cを、LATとすると、[LAT]=[LAT(1)=33.817、LAT(2)=33.468、・・・]、LONとすると、[LON]=[LON(1)=131.843、LON(2)=132.448、・・・]、DEPとすると、[DEP]=[DEP(1)=83.0、DEP(2)=48.0、・・・]、INTとすると、[INT]=[INT(1)=49.0186、INT(2)=217.7151、・・・]、MAGとすると、[MAG]=[MAG(1)=4.0、MAG(2)=4.1、・・・]、と表記される。地震の発生順序を示すインデックスjの1、2、・・・が各時系列の時間tとなる。これら震源要素cの時系列[c]を、次の数式(1)で表記する。d(MAG,1)となる。説明に戻る図3aの地震発生の変化のランダム変化を平滑し、最大の変動値が1となるように規格化する。例として、震源要素cが、INTとDEPの時系列d(INT,m)とd(DEP,m)を取り上げる。平滑は、移動平均でも累積加算でも良い。ランダム変動のみならず、予兆と予知の図3c-図3eで観測したA(c,m)の略70イベントの周期変動も平滑する場合は、移動平均個数を70イベントとしたとした時系列<d(INT,m)>の過去の最大値がそれぞれ1となるように規格化した時系列を、NCI(m,70)とする。同様に、<d(DEP,m)>を規格化した時系列を、NCD(m,70)とする。これらNCI(m,70)とNCD(m,70)は、図4に示すように、大地震発生と同期している。サイクルのピークは、大地震発生前に出現する。規格化したNCI(m,70)は、1994年4月6日に、又、NCD(m,70)は、1994年1月31日にピークに到達後、共に減少し始め、更に急減少する過程で、1995年1月17日にM7.2の兵庫県南部地震(1995)が発生した。1997年6月25日に発生したM6.6の山口県東部地震(1997)、そして、2000年10月6日に発生したM7.2の鳥取県西部地震(2000)に観察されている。2001年1月12日に発生した兵庫県北部の群発地震(S2001)直後の2001年3月23日に発生したM6.7の安芸灘地震(2001)の場合、観察されない。5.9地震の発生と、M5.4以上の地震の多発とに関連している。km)が発生している。発生と同期していて、サイクルのピークは、必ず、発生の前に起きる。(Accelerated-Moment-Release)現象[P11]である。従って、臨界サイクルをモニターし、これら大地震、巨大地震の発生時期を予測できる。サイクルのピークに到達する様子をモニターする事は、任意に選択した対象地域の地殻の応力が、臨界状態に到達しているかどうかを、監視している事になる。何時発生しても不思議でない危険な状態を意味する。比例する事が証明されている[P1]。地震発生変化のランダム現象を平滑化した移動平均個数の70をより少なくし、例えば、40とすると、AMR現象を、大地震や巨大地震が発生する時刻をより早く検出できる。国土地理院がオンライン公開している気象庁一元化処理の震源要素と電子基準点のF3座標値です。