断層モデルによる地殻変動計算プログラム
DC3D0 / DC3D

（例１） 1986年伊豆大島噴火

1986年11月伊豆大島三原山の山腹で生じた割れ目噴火の際に、国土地理院によって全島をめぐる水準測量が実施され、複雑な上下変動の様子が実測されました。
発表されたばかりのOkada (1985)による地殻変動の理論式がこれに適用され、伊豆大島の直下に長さ15km、幅10km、開口量2mのマグマ貫入が推定されました［橋本・多田(1988), 火山, 33, S136-S144］。 

この時、防災科学技術研究所が御神火茶屋に設置していた傾斜計は噴火直前に大きな傾斜変化を記録し、これもマグマ上昇に伴う浅部への貫入としてモデル化されました［Yamamoto et al. (1991), J. Phys. Earth, 39, 165-176］。 

（例２） 1989年伊東沖海底噴火

伊豆半島伊東市の沖合では、1978年以来，間欠的に激しい群発地震が繰り返されてきましたが、1989年7月には海底で小規模なマグマ水蒸気爆発が起こりました。

この事件に伴って、伊豆半島東部周辺では傾斜計、体積ひずみ計、GPSなどによって大きな地殻変動が観測され、これらをOkada (1985)に基づく理論地殻変動と比較することにより、マグマ貫入の時間的な経過を含む定量的モデルが提出されました［Okada and Yamamoto (1991), J. Geophys. Res., 96, 10361-10376］。 

（例３） 1992年ランダース地震，1995年兵庫県南部地震

1992年6月、米国カリフォルニア州で発生したランダース地震(M7.3)に伴う地殻変動が、人工衛星ERS-1に搭載された合成開口レーダ(SAR)の画像を干渉させることにより初めて捉えられ、科学雑誌Natureの表紙を飾りました。
Okada (1985)から計算される理論的な干渉画像と比較することにより、地震断層のすべり量分布の推定などが行われました［Massonnet et al. (1993), Nature, 364, 138-142］。 

同様の手法は、その後さまざまな地震・火山活動に適用されるようになり、1995年兵庫県南部地震(M7.3)でも、干渉SARによる地殻変動パターンが地震断層モデルと整合することが確かめられました［国土地理院(2000)，『地震予知連絡会30年のあゆみ』, p.20］。 

（例４） 2000年伊豆諸島群発地震

2000年6月末に始まった三宅島の火山活動は、その後、新島・神津島周辺海域に及ぶ大規模な群発地震活動に発展し、大きな地殻変動を伴いましたが、これら一連の活動は、Okada (1985)による理論式を用いて、複数の開口断層とずり断層の組合せによってモデル化されました［村上ほか(2001),月刊地球, 23, 791-800］。 

この激しい群発地震活動が、懸念される「東海地震」の発生にどのような影響を与えるかが憂慮されたため、上記のモデルによって東海地震の想定断層面上に及ぼすCFF（クーロン破壊関数）の変化（：ずり応力変化、：法線応力変化、：実効摩擦係数）が、Okada (1992)による内部変形の式を用いて評価されました（気象庁資料，2000/9/1）。

その結果、これらの地震により生じる応力変化は、日々繰り返される潮汐現象による応力の数分の１程度であり、直接的な影響はないであろうとの判断がなされました。

（例５） 2000年鳥取県西部地震

2000年10月6日に発生した鳥取県西部地震（M7.3）では、震源断層に沿った場所以外に、東西の約20km離れた位置でも顕著な余震活動が発生し、注目されました。
右の図は、本震周辺の地下8kmのレベルにおいて、本震と同じタイプの地震を起こす「ずり応力」がどのように分布するかを、Okada(1992)による内部変形の式を用いて計算したものです［遠田(2002), 地学雑誌, 111, 233-247］。

暖色系で示される地域では断層がすべり易くなり、寒色系の地域ではすべりにくくなることを示しています。震源から離れた東西の余震活動は「ずり応力」が増加した地域に一致しており、本震による応力変化によってトリガーされたものと解釈されました。このように、大きな地震のあとで余震が発生しやすい場所と、しにくい場所とを識別することができるようになりつつあります。 

（例６） 2004年スマトラ島沖地震

2004年12月26日に発生したスマトラ島沖地震（M9.1）はインド洋に大津波を起こし、数10万人もの犠牲者を出しましたが、その約3ヶ月後の2005年3月28日、すぐ東隣りのニアス島付近で再びM8.7の大地震が誘発されました（右図は、八木による両地震の震源域を示す）。

このような地震が誘発される可能性は、3月17日に発行された科学雑誌NATUREの中で事前に指摘されていました［McCloskey et al. (2005), Nature, 434, 291］。このような推論は、スマトラ島沖地震の発生によって周辺に生じる応力場をOkada(1992)による内部変形の式を用いて評価した結果、得られたものです。 

（例７） 2011年東北地方太平洋沖地震

2011年3月11日に発生した東北地方太平洋沖地震（M9.0）では、巨大な津波により2万人近い死者・行方不明者を生じました。この地震により，東日本では広域にわたって非常に大きな地殻変動を生じ、周辺の広い範囲において活発な余震活動が続きました。

このように巨大な地震の発生は、近隣地域での地震や火山噴火など、様々な地学現象を誘発することが懸念されます。
下の図は、Okada (1992)による内部変形の式を用いて、周辺海域での海溝型地震や内陸の活断層型地震に及ぼすCFF（クーロン破壊関数）の変化を評価した結果です［Toda et al. (2011), Earth Planets Space, 63, 725-730］。

暖色系で示される海域や活断層では地震の発生が促進され、寒色系の海域や活断層では地震の発生が抑制されることを示しています。 

（１） これまでに提出された地表変形および内部変形の理論式

断層モデルには、その断層面の向き（垂直，水平，斜め）や断層運動の型（横ずれ，縦ずれ，開口）によって色々なタイプがあり、また媒質としては、弾性定数が等しい場合（ポアソン固体）と任意の場合とがあります。さらに、計算すべき物理量としては、地表変形の場合、変位３成分・歪４成分・傾斜２成分、内部変形の場合、変位３成分・歪９成分が必要です。

従来、これらのすべてを網羅する完全に一般的な表現式はありませんでしたが、Okada (1985)およびOkada (1992)によって，断層面の向き・断層運動の型・媒質を選ばず、かつすべての物理量を計算できる一般的な理論式のセットが提出されました（詳細については，§１プログラム概要にある参考論文を参照してください）。 

（２） 断層モデルによる地表変形の理論式 （Okada, 1985）

Strike-slip, Dip-slip, Tensile fault による表現式に統一感があり、全体として非常に簡潔かつ美しい表現式となっている点が特徴です。 　 

（３） 断層モデルによる内部変形の理論式 （Okada, 1992）

本来であれば非常に複雑となるの計算式を、以下のに変換することで簡潔なものとし、かつ全体を一枚の表形式でコンパクトに表現した点が大きな特徴です（詳細については、2ダウンロードにある「Derivation of Table6」を参照してください）。

（４） 従来の表現式との比較例 （その1 ： 縦ずれ断層，開口断層による上下変位成分）

縦ずれ断層による地表上下変位については、ポアソン媒質（）の場合の式が Savage and Hastie (1966) により、また、開口断層による地表上下変位については、任意媒質（）の場合の式が Davis (1983) により提出されています。
これら従来の表現式とOkada (1985)の式が数学的に等価であることは証明済みです。また、両者は数値的にも同一の結果を与えることを確認しています。Okada (1985)の式には統一性があり、かつ、はるかにコンパクトな表現式となっていることがわかります。

（５） 従来の表現式との比較例 （その2 ： 横ずれ断層による水平歪成分）

横ずれ断層による地表水平歪については、ポアソン媒質（）の場合の式が 山崎(1975) により提出されています。こちらについても、従来の表現式とOkada (1985)の式が数学的に等価であることは証明され、また、両者は数値的にも同一の結果を与えることを確認しています。Okada (1985)の式ははるかにコンパクトな表現式となっていることがわかります。