大地震（本震）が起きた後これに続いて発生する余震数は、本震発生からの経過時間にほぼ反比例して減少していくことが知られています。しかし実際には余震が引き起こす余震（二次余震）が存在するなど、複雑な時間変化を示します。

　こうした余震活動の推移を再現するために、地震発生を司る摩擦の特性から応力と地震発生の関係を表すモデルが導入され、様々な工夫が試みられています。

　このモデルの従来のものでは二次余震の影響が考慮されていませんでした。そこでこれを考慮する改良を施しました。一例として2004年10月に起きた新潟県中越地震の余震データに適用した結果を示します。従来モデルでは本震発生から単純に1日あたりに発生すると期待される地震数が減少していくだけですが（図1(a)、緑線）、二次余震の影響を入れられるよう改良したことにより複雑な余震発生の様子が作り出せています（赤線）。さらに改良したモデルが実際の余震発生によく合っていることも分かります（図1(b)）。

　このようにモデルを工夫することで、より現実的な余震活動を再現することが出来ました。ところでこの例では背景となる地中の応力が時間一定で上昇すると仮定しています。しかし、大地震後、その周辺域ではゆっくりとしたすべり（余効変動）が起きることが知られており、これに対応した背景応力の変化も取り入れられるような工夫を現在試みています（図2）。このようにモデルの改良を積み重ね、現実的な地震の起こり方に近づけていくことで、将来的には地震予測に繋がると考えています。