火山 第 48 巻 (2003)
第 1 号 1ῌ10 頁

論 説

ダイクの貫入活動と断層運動の相互作用῏
有限要素法による数値シミュレῐション
和 賀 俊 治῍῎橋 本学῍῍ῒ2000 年 5 月 11 日受付ῌ 2002 年 11 月 20 日受理ΐInteraction between Dyke Intrusion and Movements of Surrounding Faults:
A Numerical Simulation with Finite Element Method
Shunji W6<6῍ and Manabu H6H=>BDID῍῍
We investigate the e#ect of movement of nearby faults caused by an intrusion of dyke on volcanic activity.
We examine two factors as indicators of the e#ect of faulting. One is a change in stress concentration at the dyke
tip and the other is a change in opening of dyke. In this study, a ﬁnite element method is applied under a twodimensional plane stress state to incorporate geometrical relationship between the dyke and the fault. Dyke andfault are modeled as an open crack with internal pressure and a shear crack with Coulomb friction, respectively.
Geometrical parameters considered are an angle, a distance between the dyke and the fault, and the length of fault.
One of the most important results is as follows; when a fault with high angle against the dyke strike extends beyond
the dyke tip along the propagation direction, and when the fault with low angle extends to it, the dyke propagate
and open more easily than the case without fault. We apply our result to recent volcanic activities in Japan
associated with the dyke intrusion, the 1986 Izu-Oshima and 1989 o# Ito eruptions, and we suggest that changes
in volcanic activities may have been attributed to the occurrence of moderate-sized earthquakes near the dyke tip.
1. は じ め にた῍ このケῐスでは M 5.5 の震源域から断層は海底噴火

ダイク貫入を伴う火山活動ではῌ やや規模の大きい地を引き起こしたとされるダイク北西端付近に推定された震の発生に伴い火山活動に変化がみられるケῐスがあ(Okada and Yamamoto, 1991)῍ これらの事例はῌ ダイク

る῍ 1986 年伊豆大島火山噴火では外輪山北西斜面からのの貫入活動が地震の発生と何らかの関係があることを示

割れ目噴火に先立って M 4.0ῑ5.1 の地震が島北西沿岸

唆する῍

付近で多発した ῒ和賀 ῎ 山岡ῌ 1996ΐ῍ Hashimoto and

地震の発生に関係するダイクの影響は次の 2 つが考え

Tada (1990) が地殻変動解析より推定した開口割れ目北られる῍ 一つはダイクの開口量の変化である (Fig. 1

西端はこの群発域付近にありῌ 特に最大地震の M 5.1 の(a))῍ ダイクの開口量が変化するとῌ ダイク中のマグマ

地震は開口割れ目北西端のすぐ西側で発生した῍圧に変化が生じる῍ これによりマグマの運動が活性化さ1989 年伊東沖海底噴火ではῌ マグマが地下から上昇をれるという考えである ῒ例えば大倉῎安藤ῌ 1994ΐ῍

続けῌ それに伴い激しい群発地震が発生した῍ 一連の活も う 一 つ は ダ イ ク 先 端 の 応 力 変 化 で あ る (Fig. 1

動中 M 5.5 の地震が発生しῌ 数日後海底噴火が発生し(b))῍ 破壊力学によればῌ クラックはクラック先端の応力集中が媒質の物性である破壊じん性を越えると伸展し῍ ῔253ῌ0006 茅ヶ崎市堤 1ῌ1ῌ8ῌ406
Tsutsumi 1ῌ1ῌ8ῌ406, Chigasaki, Kanagawa 253ῌ
0006, Japan.
῍῍ ῔611ῌ0011 宇治市五ヶ庄
京都大学防災研究所῎地震予知研究センタῐ
Disaster Prevention Research Institute, Kyoto University, Gokasho, Uji, Kyoto 611ῌ0011, Japan.
Corresponding author: Manabu Hashimotoe-mail: hasimoto@rcep.dpri.kyoto-u.ac.jp

はじめる῍ ダイク付近に存在する断層がダイク先端付近の応力集中を変化させῌ 結果としてダイクが拡大もしくは停止することが考えられる῍
ダイクは断層との相互作用によりダイク先端の応力集中および開口量に変化が生じῌ 結果としてダイクの貫入
活動が変化することが予想される῍ そこでῌ 本論はこれらに注目して定量的な評価を行う῍
ここではῌ 断層とダイクを簡単な 2 次元クラックで表和賀俊治῎橋本 学

2Fig. 1. Factors which are likely to change activity
of dyke due to an interaction between the dyke
and fault, (a) change of a dyke shape, (b)change in stress concentration at a dyke tip.

現しῌ ダイク近傍に存在する断層を想定する῍ 有限要素法によりクラックの応力集中および形状を計算しῌ 得られた結果から 1986 年伊豆大島噴火および 1989 年伊東沖Fig. 2. Description of adaptive gap elements in
MSC/NASTRAN for the judgment of contactof both planes of the fault.

海底噴火のケ῏スでみられた地震とダイクの貫入活動の関係を考察する῍

大きい場合に接触状態とみなす῍ このときῌ アダプティ
ブ GAP 要素は周囲より大きい軸剛性 (Ka) をもつ状態

2. 計 算 方 法になる῍ 一方ῌ 差が Uo と比べて小さい場合ῌ 開口状態とダイクの形状およびダイク先端の応力集中を詳細に調みなしῌ アダプティブ GAP 要素は小さな軸剛性 (Kb) をべるためῌ 本論では 2 次元平面応力を仮定した非線形静解析によりこれらを評価する῍ 数値計算はῌ 京都大学大型計算機センタ῏の有限要素法プログラム MSC/NASTRAN を使用した῍持つ状態になる῍
2ῌ3 せん断クラックの滑り基準
MSC/NASTRAN はῌ 滑り基準として Coulomb の摩擦法則を採用している῍ 2 次元平面のときῌ アダプティ

2ῌ1 ダイクのモデル化

ブ GAP 要素における力の成分はῌ 要素座標系におけるダイクは開口クラックでモデル化する῍ 開口クラック

変位で表現される軸圧縮力 (Fx) およびせん断方向の摩面上の節点は同座標に異なる節点をもつ二重節点とし擦力 (Fy) でῌ これらはアダプティブ GAP 要素の状態にてῌ これらに等分布荷重を与えることで内圧を表現す

より次のように計算される῍ 接触がなく GAP が開いてる῍ Aki et al. (1977) はῌ キラウエア火山下のマグマ圧いるときῌを 2 MPa と見積もった῍ Wadge (1977) はῌ エトナ火山

FxΐKb uxῌ0 および Fyΐ0(1)における側噴火の圧力として 15 MPa と推定した῍ しかでせん断剛性がない῍ ここでῌ Kb は開いている状態でのしマグマ圧はよく分っておらずῌ 観測値にもばらつきが軸剛性である῍ GAP が閉じて固着しているときῌ

ある῍ ここではῌ これらの事例を参考にしてマグマ圧を10 MPa と仮定した῍

Fyῌms Fx(2)

ここでῌ FxΐKa ux, FyΐKt uy でῌ Ka は閉じている状態で

2ῌ2 断層のモデル化の軸剛性ῌ ms は静摩擦係数ῌ Kt はせん断剛性である῍ 一

断層はせん断強度をもつせん断クラックでモデル化す

方ῌ GAP が閉じて滑りがあるときῌ

る῍ せん断クラックの内面は X, Y 軸方向ともに変位拘束されない自由表面として表現されῌ 固着῎すべり῎開口の変形が生じる῍ MSC/NASTRAN ではῌ 2 節点間の断層面の接触状態をシミュレ῏トするのにアダプティブ
GAP 要素を用いる ῑ日本エムエスシ῏ῌ 1994ῒ῍

せん断クラック面の接触はῌ 2 節点間に配置されたアFy῔mk Fx(3)

ここでῌ FxΐKa ux でῌ mk は動摩擦係数である῍
2ῌ4 モデルῌパラメ῍タ

計算に用いた各パラメ῏タを Table 1 に示す῍ ダイク
貫入に伴う地震の震源は数キロと浅いものが多い ῑ山
岡῎他ῌ 1988; Ukawa and Tsukahara, 1996ῒ῍ この深さでダプティブ GAP 要素の状態で判別される῍ アダプティ

は断層面に作用する垂直応力は 100ῐ200 MPa と推定さ

ブ GAP 要素の概念を Fig. 2 に示す῍ アダプティブ GAPれる ῑ山口῎西松ῌ 1991ῒ῍ ここではῌ Byerlee (1978) に要素はῌ 節点 Ga, Gb を通る方向を X 軸としてῌ X 軸方

よる法線応力とせん断応力の関係を考慮して静摩擦係数

向の変位差 (Ua῍Ub) が予め設定した初期開口量 Uo よりを 0.6 とする῍

ダイクの貫入活動と断層運動の相互作用 有限要素法による数値シミュレション

3

Table 1. Model parameters adopted inthis analysis.

Fig. 3. Idealized model and coordinate systems.

一般に 滑りが生じている時の動摩擦係数は静摩擦係数よりも小さい そこで 滑りが生じる状態をより現実

Table 2. Geometric parameters used in this analysis.

に近づけるため動摩擦係数を採用し その値を 0.4 とした
マグマは地下数キロにあるマグマ溜まりから上昇すると考えられる ここでは この深さを考慮して地殻の剛性率として 20 GPa を採用した また ポアソン比は0.25 を仮定したがない場合の開口クラックのモド I 型応力拡大係数
KῑI により正規化した値 KI/KῑI について評価する これ

2ῌ5 応力拡大係数の評価法

により KI/KῑI 1 のとき 開口クラックはせん断クラッ

二次元平面においてクラック先端付近の変位分布は

クがない場合と比較して伸展しやすい状態 K I /K ῑI

クラックの代表寸法と比較して十分小さい範囲では円筒座標により近似的に次式で表される

1のとき 伸展しにくい状態にあると解釈できる
2ῌ6 座標系と幾何学パラメῌタui(Kn/2G) (r/2p)1/2 fi(q) ix, y nI, II, III (4)計算に用いた開口クラックおよびせん断クラックの座

ここで Kn は応力拡大係数で n は変形モド G およ

標系を Fig. 3 に モデルの幾何学パラメタを Table 2

び fi(q) はそれぞれ剛性率変位関数 q は X 軸から右に示す 開口クラックおよびせん断クラックの長さをそ周りに測った角度 i は X 軸又は Y 軸方向 r はクラッ

れぞれ Lo, Ls と表現し 開口クラックとせん断クラッ

ク先端からの距離を示す 変位関数 fi(q) は開口モド

クの位置を表現する XῐY 座標と 開口クラック せん断

型 (nI) について

クラック上の節点の位置を表現する XoῐYo 座標 Xsῐ
Ys 座標をそれぞれ独立に与える f は X 軸に対する角q῎qῌ
fx(q)cosk12 sin22῏
2῍fy(q)sinqῌq῎k12 cos22῍
2῏

度を表し 開口クラックの走向から反時計周りに計った(5)k(3v)/(1v)と表現される v はポアソン比である (4) および (5) 式

値に相当する
数値計算を行うにあたって 4 点四角要素を使用した
節点数および要素数はモデルにより若干異なり それぞれ 1400 2000 および 1200 1800 である 数値計算の精よりクラック先端付近の変位は座標により一義的に決ま度を確保するため 開口クラック周辺の要素は正方形にる そこで 計算により求められた節点の変位を上式になるよう配慮した モデル外部境界は 境界面に直交す

代入して KI を評価する 変位は開口クラック端の手前る方向の変位を拘束した 外部境界は開口クラック端のに設定された節点を選んだ この節点とクラック端まで応力変化に影響を及ぼす これらを避けるため 開口クの距離と 開口クラックの長さ比は 1/40 である

ラック端から外部境界までの距離を開口クラックの長さ

地殻の破壊じん性はクラックの伸展 または停止を決の 5 倍とした

定する重要な物性値と考えられる 過去の研究で火山岩

数値計算の効率を高めるため 開口クラックとせん断の破壊じん性を見積もった事例があるが 例えば Atkin-

クラック端との間に Y 軸方向へ 0.15 Lo の距離をとっson, 1984  実際にはよく分かっておらず KI の絶対値た なお 開口クラックとせん断クラックの距離 せん

で評価することは困難である そこで せん断クラック

断クラックの長さはそれぞれ Lo で正規化をおこない和賀俊治 橋本 学

4計算された開口量についても せん断クラックがない場

は殆ど変化がなかった このように 開口クラックとせ合の開口クラックの Xo/Lo0.5 に配置された節点のん断クラックのなす角が 90に近いほど開口クラックは

開口量 Yo () で正規化した

伸展しやすい状態になり これに対応して開口量も増加

計算は まず X/Lo0.0, Ls/Lo1.0, に固定して fする

を 0150 の間で 15 刻みに変えて行った 次に X/

3ῌ2

Lo を 1.00.5 の範囲で変えて せん断クラック長にmodel-B

KI/KῌI の計算結果を Fig. 5(a) に示す X/Lo

0 のとよる効果も考慮して Ls/Lo0.25, 1.0, 4.0 の 3 ケスをき f15と f45, 75, 90 の X/Lo1.0 では

それぞれについて計算した f は 15, 45, 75, 90 の 4

KI/KῌI にほとんど変化がなく 1 に近い値になる それ以ケスを想定した

外は KI/KῌI

3. 解 析 結 果

大きくなる

1 で 開口クラックが伸展しやすい状態になり X/Lo がゼロに近づくにしたがい KI/KῌI は急激に
3ῌ1model-A

KI/KῌI の計算結果を Fig. 4(a) に示す 計算した f の範囲内では KI/KῌI 1 であり 開口クラックは伸展しや

例えば Ls/Lo1.0 として f45, 75, 90 のケス
をみると f は低角よりも高角のほうがせん断クラックの効果が表れて KI/KῌI が大きくなる

すい状態にある KI/KῌI は f90で最大となり 開口ク

また 例えば f45の Ls/Lo0.25, 1.0, 4.0 と f90

ラックとせん断クラックが直交するように配置されるのそれを比較すると 後者のほうが Ls/Lo が大きいほど時 せん断クラックの効果が最も大きい

KI/KῌI はより大きくなる傾向にあり せん断クラックの開口量の分布を Fig. 4(b) に示す 開口量は f90の長さによる影響が表れる なお Ls が Lo より短い場合

とき最大で KI/KῌI の最大に対応する 一方 f が低角にを想定して Ls/Lo0.25 を計算したが この場合 f やなるほど開口量は小さくなり f15, 45 のモデルでX/Lo の違いに関係なく KI/KῌI は 1 に近い値となった

Fig. 4. Results in model A, (a) stress intensity factors, and (b) deformations of both planes of the open crack.
They are normalized by the stress intensity factor and maximum opening value (Xo/Lo0.5) without a shear
crack respectively.

ダイクの貫入活動と断層運動の相互作用 有限要素法による数値シミュレション

5Fig. 5. Results in model-B, (a) stress intensity factors, and (b) deformations of both planes of the open crackwith Ls/Lo1.0, respectively. They are normalized by the stress intensity factor and maximum opening value(Xo/Lo0.5) without a shear crack, respectively.

X/Lo0 では せん断クラックの伸長方向は開口ク

Xo/Lo1 の範囲で 開口クラック端付近に限られたラックと交差する配置になる この配置では f によりf が低角だと X/Lo0.5 で大きな開口を伴うが その

K I /K ῌI が異なる f15 のケスでは K I /K ῌ1となI

分布はせん断クラックに近い一部の範囲に限られる 例り このとき X/Lo が大きいほど KI/KῌI は急激に大きくえば f15のケスでは 0.4 Xo/Lo1.0 の範囲でなる 例えば f15, X/Lo0.5 では KI/KῌI が約 40%

開口が生じ せん断クラックがない場合と比べて開口量増加した これに対して f45, 75, 90 になると せは約 2 倍になった 一方 f が高角では X/Lo0 で開口

ん断クラックの伸長方向が開口クラック端近くで交差す

量は最大になる この場合 f が低角のケスに比べてる配置となる X/Lo0.125, 0.25 で K I /K ῌI は減少の傾

開口クラック全体が開口する傾向がみられた

向となり このうち f75, 90 のケスでは X/Lo
0.25 で KI/KῌI 1 となった 一方 X/Lo0.25 になると

4. 議

K I/K ῌI は増加の傾向となり f45, 75 のモデルでは

本論は 簡単な二次元クラック モデルを用いてせん

K I /K ῌ1となった Ls/Lo についてみると Ls/Lo
I

断クラックの影響による開口クラックの開口量および先論

1.0, 4.0 では KI/KῌI がほとんど同じ値になり せん断ク

端の応力集中の変化を調べた 本論は開口クラックに伴

ラックの長さにあまり影響されない

う静的な変形について計算したものであり せん断ク

開口量の分布を Fig. 5(b) に示す ここには X/Lo

ラックの動的な破壊プロセスによる効果について考慮し0, 0.25, 0.5 の計算結果について示した これらのうち

ていない マグマの移動速度 (Hauksson, 1983), 割れ目

開口量が減少するのは f75, 90 の X/Lo0.25 の噴火の拡大速度 荒牧 早川 1984; 阿部 1988  震源となるケスに対応する
ケスで これらは KI/Kῌ1Iの移動速度 山岡

開口量の減少は例えば f90, X/Lo0.25 では 0.7

son, 1979

他 1988; Brandsdottir and Einars-などから ダイクの拡大速度はおよそ 0.2

和賀俊治῎橋本 学

6Fig. 5.
0.7 m/s のオ῏ダ῏と推定されるのに対してῌ 断層の平均的な破壊伝播速度は 1ῐ4 km/s であり ῑ例えばῌ 阿Continuedの距離とῌ 双方のなす角に依存する῍
結果を考慮するとῌ 例えばῌ ダイク周辺にダイク長と部ῌ 1988ῒῌ ダイクの拡大速度は地震の破壊伝播速度に比同程度以上の長さをもつ断層がありῌ 双方のなす角が小

べてきわめて小さいと考えられる῍ よって断層運動の動さくῌ 断層の伸長方向がダイクと交わる配置でῌ またῌ

的効果は無視できῌ ここで示した計算結果は第一近似的双方のなす角が大きくῌ 断層の伸長方向がダイクと交わに有効であると考える῍

らない配置であるとῌ ダイクの開口量が増大する῍ 前者

結果を Fig. 6 に模式的に示す῍ せん断クラックによるのケ῏スではῌ ダイクの一部でῌ 後者のケ῏スではダイ

開口クラック開口量 (Fig. 6(a)) および先端の応力集中

クの広範囲で開口を伴うことが考えられる῍ 一方ῌ 双方(Fig. 6(b)) はῌ 主として開口クラックとせん断クラックのなす角が大きくῌ 断層の伸長方向がダイクの先端付近

ダイクの貫入活動と断層運動の相互作用῏ 有限要素法による数値シミュレῐション

7

Fig. 6. A schematic idea of the results in this paper, (a) deformation of the dyke and (b) a change in stressconcentration when a fault exists around the dyke.

で交わる配置だとῌ ダイク端の応力集中が弱められて拡大が停止することが考えられる῍

活動の説明を試みる῍
4ῌ1

1986 年伊豆大島噴火

以上よりῌ ダイクと地震の発生場所からダイクが受けこの事例ではῌ 大島を北西ῌ南東方向に横切るダイク

た断層の影響を検討してみる῍ ここではῌ ダイクの貫入の貫入が推定された ῒ山岡 ῎ 他ῌ 1988; Hashimoto and活動に伴い規模の大きな地震が発生した 1986 年伊豆大Tada, 1990ΐ῍ またῌ 割れ目噴火開始後ῌ 大島北西部にお島噴火と 1989 年伊東沖海底火山噴火の事例を取り上げ

い て M 5.1 を 最 大 と す る や や 規 模 の 大 き な 地 震る῍ 特にダイクと断層の位置に注目しῌ その配置から期(M 4ῑ) が多発した (Fig. 7(a))῍ これらのほとんどはῌ待されるダイクの貫入活動の変化を検討して一連の火山推定されたダイク北西端付近から大島北西岸付近の狭い和賀俊治橋本 学

8Fig. 7. (a) A distribution of epicenters and a dyke location, (b)the geometrical location between the fault anddyke associated with the 1986 eruption of Izu-Oshima volcano. Gray and black lines denote the dyke location
inferred from Hashimoto and Tada(1990). The epicenter of the M 5.1 earthquake and a region (dashed line)
west o# the shore, where earthquakes (M 4.0) swarmed after the beginning of the ﬁssure eruptions, wasinferred from Waga and Yamaoka (1996). Others were determined by E. R. I. Slip direction of the fault in (b)was referred to Yamaoka et al. (1988).

範囲で発生しており 和賀山岡 1996 この領域で断議論はできない Fig. 7(a) では ダイク北西端は元町直層が形成されたことを示す M と断層長の経験則から推下付近に求められ (Hashimoto and Tada, 1990), ダイク

定すると M 5.1 の地震では約 4 km の断層が形成された端付近に断層があるようにみえる これは ダイクが拡と考えられる

大を続けて断層と交差する配置に至り 結果 ダイク先

この地震が発生した時 大島北西部に貫入したダイクは割れ目噴火開始直前の地震活動域 山岡他 1988

端の応力集中が弱められ 拡大が停止したものと解釈できる

噴火割れ目の位置長さ 阿部 1988 から長さ約 234ῌ2

km 程度と推定され 大島北西部の断層はダイクの 2 倍

この事例では 地震活動およびダイクの貫入の様子が程度の長さだったと考えられる 山岡他 (1988) によ

時間変化とともに報告されている (Ukawa and Tsuka-

1989 年伊東沖海底噴火

れば 地震が群発した領域は北西ῌ南東を主圧縮軸とす

hara, 1996; Okada and Yamamoto, 1991) Fig. 8(a) はるストライクスリップ型が支配的なメカニズム解であOkada and Yamamoto (1991) が推定したダイクと最大る これらを考慮すると断層とダイクのなす角は約 40

地震 (M 5.5) の配置である これによれば 約 3 km のとなる

ダイクと約 6 km の断層が約 20 という低角でダイクの以上のように解釈すると ダイクと断層の配置は Fig.

北西端付近で交差する (Fig. 8(b)) この配置は Fig. 57(b) になる この図の配置に従うと Fig. 5(a) の f(b) の f15で X/Lo0.25 に近いと考えられ ダイク

スに近く 断層の伸長方

北西端付近に局部的な開口が推定される 地震発生の数45で X/Lo

0.125 でのケ

向はダイクと交わらない この配置だと ダイクは拡大

日後に海底噴火が発生したが その位置はダイク北西端

しやすい状態になり大島北西部への貫入が推定される付近だった 開口量の変化に伴うマグマ圧の変化を考慮

カルデラ北部で始まった割れ目噴火は M 5.1 の地震発する必要があるが Okada and Yamamoto (1991) が指生 30 分後に外輪山北側斜面でも割れ目噴火が生じた

摘したように 断層が噴火に寄与した可能性が考えられ

これは 地震発生後もダイクが拡大を続けたことを示しるており ダイクが断層の影響を受けた可能性がある た
だし 今回のモデルでは位置に不確定性があり決定的なダイクの貫入活動と断層運動の相互作用῏ 有限要素法による数値シミュレῐション

9Fig. 8. (a) A spatial distribution of fault (M 5.5) and dyke locations; (b) the geometrical locations (twodimensional model) between the fault and dyke associated with the 1989 eruption o# Ito (Okada andYamamoto, 1991).

5. 今後の課題(1)

断層の伸長方向がダイクと交わらない配置でῌ

本論はῌ 簡単な 2 次元クラック῎モデルにより断層が双方のなす角度が高角かῌ もしくはῌ 断層の伸長方

ダイクに与える影響を評価しῌ ダイクの活動を伴った火向がダイクと交わる配置で双方のなす角度が低角の

山活動の変化を検討した῍ 本論で示したようにῌ ダイク

場合ῌ ダイクは拡大しやすい状態になりῌ このときと断層の幾何学的な配置はῌ ダイクの活動に影響を及ぼす要素と考えられる῍ ダイクの受ける影響はῌ 特にῌ 断層がダイク端付近に位置する場合にあらわれῌ 互いの配置のわずかな違いで受ける影響は異なることからῌ 両者の関係をさらに調べるには断層とダイクの位置を詳細に把握する必要がある῍

開口量も増大する῍(2)断層の延長がダイク端で交わる配置でῌ ダイク

は拡大しにくい状態になるがῌ 双方のなす角度が高角ほどその効果が大きい῍
(3)断層とダイクのなす角が高角ほど断層の長さによる影響が大きくなる῍

またῌ 本論はダイクの形状と応力集中の変化に着目したがῌ ダイクと断層の相互作用により断層自身も影響を謝受けてῌ これが地震の発生と関わっている可能性があ

本研究は和賀が京都大学大学院在学中および京都大学

辞る῍ さらにῌ 今回の計算ではῌ マグマ圧の変化῎地殻の防災研究所における研修として行ったものである῍ 本研

破壊強度などを考慮していない῍ 加えてῌ 火山活動は広究を行うにあたってῌ 石原和弘教授ῌ 岡田篤正教授はじ域応力場に強く支配されると考えられることからῌ 今

め多くの方の助言を頂いた῍ 特にῌ 京都大学大型計算機

後ῌ これらの要素も含めて詳細な関係を考察する必要がセンタῐ῎伊藤彰朗技官には汎用プログラムの扱い方である῍

ご指導を頂いた῍ 名古屋大学大学院環境学研究科地震῎
火山観測研究センタῐ῎山岡耕春助教授には議論で有益

6. ま と めな助言を頂いた῍ 東海大学海洋研究所῎楠本成寿博士にダイク周辺の断層がダイクに与える影響に着目しῌ ダは有限要素法の助言や動機付けを頂いた῍ この場を借りイクと断層間の距離および双方のなす角度について有限要素法により影響を評価した῍ 得られた結果を以下にまとめる῍て感謝致します῍和賀俊治῎橋本 学

10

引用文献
阿部勝征 (1988) 伊豆大島 1986 年割れ目噴火の拡大速度῍ 火山ῌ 33, S16῎S19.
Aki, K., Fehler, M. and Das, S. (1977) Source mechanismsof volcanic tremor : Fluid-driven crack models and theirapplication to the 1963 Kilauea eruption. J. Volcanol.
Geotherm. Res., 2, 259῎287.
Atkinson, B. K. (1984) Subcritical crack growth in geologic materials. J. Geophys. Res., 89, 4077῎4114.
Ando, M. (1974) Seismo-tectonics of the 1923 Kanto
Earthquake. J. Phys. Earth, 22, 263῎277.
荒牧重雄 ῎ 早川由紀夫 (1984) 1983 年 10 月 3ῌ4 日三宅島噴火の経過と噴火様式῍ 火山ῌ 29, S24῎S35.
Brandsdsdottir, B. and Einarsson, P. (1979) Seismic activity associated with the September 1977 deﬂation of
Kraﬂa volcano in north-eastern Iceland. J. Volcanol.
Geotherm. Res., 6, 197῎212.
Byerlee, J. D. (1978) Friction of rock. Pageoph, 116: 615῎
626.
Hashimoto, M. and Tada, T. (1990) Crustal deformations
associated with the 1986 ﬁssure eruption of Izu-Oshima
volcano, Japan, and their tectonic signiﬁcance. Phys.
Earth Planet. Inter., 60, 324῎338.

Hauksson, E. (1983) Episodic rifting and volcanism at
Kraﬂa in North Iceland: Growth of large ground ﬁssures
along the plate boundary. J. Geophys. Res., 88, 625῎636.
日本エムエスシῐ (1994) MSC/NASTRAN 非線形解析
ハンドブック῍ 4 章ῌ 24῎37.
大倉敬宏῎安藤雅孝 (1994) 1991 年ピナツボ噴火は 1990
年フィリピン地震に引き起こされたか῏ 地学雑ῌ
103, 464῎470.
Okada, Y. and Yamamoto, E. (1991) Dyke intrusion modelfor the 1989 seismovolcanic activity o# Ito, central
Japan. J. Geophys. Res., 96, 10361῎10376.
Ukawa, M. and Tsukahara, H. (1996) Earthquake swarms
and dike intrusions o# the east coast of Izu Peninsula,central Japan. Tectonophysics, 253, 285῎303.
Wadge, G. (1977) The storage and release of magma on
Mount Etna. J. Volcanol. Geotherm. Res., 2, 361῎384.
和賀俊治῎山岡耕春 (1996) 割れ目噴火に伴う地震の特質῍1986 年伊豆大島火山噴火の場合῍῍ 火山学会秋季
大会講演予稿集ῌ A4.
山岡耕春῎渡辺秀文῎坂下至功 (1988) 1986 年伊豆大島噴火前後の地震活動῍ 火山ῌ 33, S91῎S101.
山口梅太郎῎西松裕一 (1991) 岩石力学入門῍ 東京大学
出版会ῌ 331 pp.
ῑ編集担当 松本 聡ῒ