地震データの統計的性質とは、地震の発生時刻、発生間隔、マグニチュード、空間分布などに対して統計学的手法を適用し、その特徴や変動傾向を分析する研究分野である。地震学では、確率分布、探索的データ解析（EDA）、時系列解析、Q–Qプロットなどの統計的手法が、地震活動の記述や地域差の把握、地震発生過程の理解を目的として広く利用されている。これらの手法は、観測データの整理や特徴抽出、異常の検出などに用いられ、地震学の補助的な分析手段として位置づけられている。
歴史・背景 地震データに対する統計的な研究は、20世紀初頭から徐々に発展してきた。初期の研究では、地震の発生頻度や規模の分布を記述するために、経験則や単純な統計量が主に用いられていた。その後、観測網の整備と地震カタログの充実に伴い、地震活動を確率論的に扱う試みが進み、発生時刻やマグニチュードの分布特性を解析する研究が広がった。
20世紀後半には、地震学における統計解析はより体系化され、時系列解析、空間統計、確率過程モデルなどが導入されるようになった。また、地震活動の地域差や時間変化を把握するために、分布の適合度検定や探索的データ解析（EDA）などの手法が応用され、地震データの特徴を多面的に評価する枠組みが整えられた。
21世紀に入ると、計算機環境の発展と統計解析ソフトウェア（Rなど）の普及により、大規模な地震カタログを対象とした詳細な解析が容易になった。これにより、地震発生過程の理解や地震活動の変動性を評価する研究が進展し、統計学的手法は地震学の補助的な分析手段として定着している。
概要 地震カタログに統計的手法を適用する研究が進展し、地震の再発間隔、マグニチュード、震源分布に関する特徴的な統計的特性が報告されている。これらの研究は、従来の地震活動モデルや経験則とは異なる分布特性や時間的変動を明らかにすることがあり、地震発生過程の理解や前兆パターンの評価に関する議論に影響を与えている。
データソース
気象庁地震カタログ 公開地震データベース
統計的性質
地震発生間隔の分布 サムネイル|代替文=QQplots|A. 指数分布をｘ軸に、ソートしたデータをy軸にしたQ-Qplot. これが直線になるのは、データがその分布に従うことを示す。 B. 正規分布をｘ軸に、ソートしたデータをy軸にしたQ-Qplot. これが直線になるのは、データがその分布に従うことを示す。 C. 正規分布する乱数のヒストグラムを片対数表示した。データのなかで特に大きいところだけがリニアに減衰しているようにみえるが偶然であって、数学的な意味はない。
日本の気象庁地震カタログを用いた解析では、地震の発生時刻を時系列に並べ、隣接する地震の発生間隔を調べると、発生間隔が指数分布に従うことが報告されている。 指数分布は、事象が独立かつ記憶を持たないランダム過程で生じるときに現れる分布であり、この結果は「個々の地震は本質的にランダムな現象として扱える」と解釈されている。
マグニチュードの分布 地震カタログの解析では、有感地震のデータを選別せずに扱うと、マグニチュードの分布が正規分布に従うと報告されている。 これは、マグニチュード頻度が片対数グラフ上で直線的に減少するとする グーテンベルク–リヒター則とは対照的である。
一部の研究によれば、片対数プロットで見られる直線性は グラフ上の見かけの効果であり、特に小規模地震を大幅に除外した場合にのみ現れるとされる。こうした選択的なデータ利用は チェリー・ピッキングとみなされ、統計的に誤解を招く可能性が指摘されている。
グーテンベルク–リヒター則は地震のエネルギー推定の基礎として広く用いられてきたため、もし実際の分布が正規分布的であるなら、従来の解釈の再検討が必要になる可能性がある。正規分布はランダム数の和の分布であり、マグニチュードは地震のエネルギーの対数値なので、エネルギーは複数の因子の相乗作用で決定されていることになるだろう。
空間分布（震源の三次元構造） サムネイル|A.3次元表示した震源の位置。 B. 瀬戸界面と、下に潜り込んでいる太平洋プレートの関係。 これらの動かせるバージョンのhtmlファイルはここにある。C. 界面の位置。方程式で表せるので、地図上でもまぎれなく位置が表示できる。 D. 界面からみえるプレートの位置。Barnes 2003の推定とはかなり異なる結果となっている。 日本周辺の地震の多くはプレート境界で発生する。震源の深さ・緯度・経度を三次元的に可視化する解析では、日本周辺の地震活動が 東日本と西日本の二つの主要な境界領域に分類されると報告されている。これらの境界は概ね平面として扱える。
東日本の境界 a. 北海道地域 — 千島海溝沿い。太平洋プレートが北米プレートの下に約 42° で沈み込む。 b. 三陸地域 — 日本海溝沿い。太平洋プレートがユーラシアプレートの下に約 31° で沈み込む。 c. 小笠原地域 — 伊豆・小笠原海溝沿い。太平洋プレートがフィリピン海プレートの下に約 46° で沈み込む。
西日本の境界 d. 九州地域 — フィリピン海プレートがユーラシアプレートの下に約 80° で沈み込む。 e. 琉球地域 — 南西諸島海溝沿い。沈み込み角は約 64°。 f. 南西諸島地域 — 先島諸島付近で約 63° の沈み込み。
西日本の境界は東日本よりも不規則で、複数のセグメントに分かれるとされる。
東西日本をつなぐ浅い境界 両地域は浅い境界でつながっているとされ、この部分では太平洋プレートがフィリピン海プレートの下に沈み込む一方、フィリピン海プレートはユーラシアプレートの下に沈み込めないと解釈されている。そのため、西日本のプレート境界は直線的に沈み込まず、傾斜した構造を示す。
g.瀬戸地域 — フィリピン海プレートが沈み込めない浅い境界であり、圧縮によりエネルギーが蓄積する可能性が指摘されている。
幾何学的特徴: 境界がほぼ平面であるため、主成分分析（PCA）により重心や法線ベクトルを推定でき、境界面を方程式として表現できる。これにより、日本周辺のプレート配置を明確に示すことが可能となる。
南海トラフ仮説への含意 一部の研究では、これらの結果が従来の 南海トラフ巨大地震の想定位置と整合しない可能性があると指摘されている。 従来の南海トラフは太平洋側に位置するとされるが、この領域はフィリピン海プレートがユーラシアプレートの下に沈み込まないとされる。最も近い境界は 瀬戸境界であり、神戸・大阪・京都などの大都市圏に近い。
1995 年の阪神・淡路大震災はこの瀬戸境界でおきている。しかしこの地震でも、統計的異常が観測されたと報告されている。
時間変化と異常
大地震前の統計的異常 代替文=QQplotの曲がり|サムネイル|200x200ピクセル|A. 2011年東北地方太平洋沖地震の前に観測された地震発生間隔のQ–Qプロット。原点付近のデータ密度が高い部分で曲がりがみられ、この特徴は両対数グラフ（挿入図）でより明瞭になる。 B マグニチュードのQQplot。東北地方で性質が変化しているために曲がっている。東北地震はきわめておおきなマグニチュードであったが、さほどにおおきな外れ値ではないことがわかる。 地震データの分布特性を把握することは、地域間の比較や観測値の評価を一貫して行ううえで重要である。通常時には、地震の発生間隔やマグニチュードの分布は比較的安定しているとされる。しかし、大地震の前にはこれらの分布に顕著な変化がみられると報告されている。
一部の研究では、大地震の前に特定地域で地震回数が増加し、発生間隔が短くなること、また別の地域ではマグニチュードの平均値や標準偏差が上昇することが指摘されている。これらの変化は Q–Qプロットの曲がりとして現れ、より大きな地震の発生可能性が高まっている状態を示す統計的異常として解釈されることがある。
一方、火山活動に伴う地震では、類似した規模の地震が多数発生する傾向があり、その結果として発生間隔は短くなるものの、マグニチュードの標準偏差は小さくなるKonishi, Tomokazu. Earthquake Swarm Activity in the Tokara Islands (2025): Statistical Analysis Indicates Low Probability of Major Seismic Event. GeoHazards 2025, 6, 52. https://doi.org/10.3390/geohazards6030052。この特徴は、大地震前にみられる分布変化とは異なる振る舞いであるKonishi, Tomokazu. "Seismicity of the Noto Peninsula: Spatial Patterns, Shallow Seismic Zones, and Possible Volcano-Related Signals." Trends in Ecological and Indoor Environmental Engineering, vol. 4, no. 2, 2026, pp. 22–34. doi:10.62622/TEIEE.026.4.2.22-34.。
余震の減衰 代替文=地震の収束|サムネイル|200x200ピクセル|地震の収束。A.時間ごとの回数。B.マグニチュードの平均。 大地震後には余震が発生し、その時間減衰は伝統的に 大森–宇津式で表され、余震発生率は時間の逆数に比例して減少するとされる。
しかし、一部の解析では、観測される減衰が半対数グラフ上で直線的に見えると報告されており、これは半減期的な減衰に類似すると解釈されている。この挙動は、ばねの減衰のような緩和過程に近い物理機構を示唆する可能性がある。
また、単位時間あたりの平均マグニチュードも半対数的に直線的に減衰するとされるが、その半減期は頻度のものとは異なる。これは、地震の頻度とエネルギーが異なるメカニズムで制御されている可能性を示すと解釈されている。
減衰が弱い場合、十分なエネルギーが解放されていない可能性があり、比較的短期間で大きな地震が発生する可能性が指摘されている 。
批判と限界 これまで考えられてきたものは、Exploratory Data Analysis (EDA)Tukey, J.W. Exploratory Data Analysis; Addison-Wesley Pub. Co.: Reading, MA, USA; London, UK, のような近代的な統計学を使ってこなかった。そのためにデータの分布や変化について、あまり厳密に調べることができなかったのだと思われる。Rが普及することで、こうした方法への障壁は下がっている。
現時点では予知の経験が不足していて、どの程度に偽陽性や偽陰性が起きるのかがまだ判断できない。これがもっともよくある批判であろう。まだ不明なこともおおい。
関連項目
地震学 地震予知 プレートテクトニクス
脚注