火山 第 61 巻 ( 2016)
第 2 号 449-457 頁

解説・紹介

楕円体体積震源モーメントテンソル評価のための Web 計算ツール
─ 3 つのモデルの順計算と逆計算：（1）膨張，（2）流体移動，（3）流体移動後の圧力回復─
水 野 尚 人＊・亀伸 樹＊＊・市 原 美 恵＊＊（2015 年 10 月 19 日受付，2016 年 1 月 15 日受理）

Ellipsoidal Volume Source Calculation Webtool for Moment Tensor Evaluation ofThree Models: (1) Expansion, (2) Simple Movement of Magma and(3) Pressure Recovery after Simple Movement
Naoto MIZUNO＊, Nobuki KAME＊＊ and Mie ICHIHARA＊＊
1．はじめに : 楕円体体積震源のモーメントテンソルと MTCR ダイアグラム
Mizuno et al. (2015) は，2 つの楕円体空隙間をマグマが裂 P : Penny Shaped Crack）の MT の線形結合の取り得る範囲である三角形領域 SPC と比べると，1 つの EX 震源の MT の取り得る範囲 SXPXBS はこれより小さい範囲

移動する火山震源モデルを考えた．圧力平衡と質量保存になり，2 つ以上の EX 震源のトータル MT の取り得るの流体力学的拘束条件を考慮して，その震源過程がモー範囲は三角形領域 SPC となり慣習的範囲と一致する．

メントテンソル（Moment tensor, MT）でどのように表現

これにより楕円体空隙が膨張する震源過程の MT の取りされるかを理論的に考察した．そこでは，弾性体中にあ得る範囲が明確に限定された．

る 1 つの楕円体空隙が単純に膨張した際の体積震源のMizuno et al. (2015) では，次に，2 つの空隙間をマグマ

MT 表 現（Eshelby, 1957 ; Aki and Richards, 1980, 2002,

が単純に（マグマの弾性変形のみ考慮するという意味）

2004 ; Davis, 1986）の理論研究成果に基づき，空隙間の移動して圧力平衡に達するモデルを考えた（Simple Move-

流体力学相互作用の理論を展開している．

ment : SM）．高圧の楕円体から低圧の楕円体にマグマがMizuno et al. (2015) では，まず，1 つの楕円体の単純な移動した結果の 2 つの空隙のトータル MT を，マグマの膨張（Expansion : EX）による MT の取り得る範囲を解析解

体積弾性率 km が周辺岩石の体積弾性率 k より小さい場（Eshelby, 1957 ; Davis, 1986）を用いて MTCR ダイアグラ合（3-2 節）に，楕円体形状の関数として理論式を導いム 図（Moment-tensor component-ratio diagram : Trasatti etた．まず，特別な場合として高圧楕円体を球形状に限定

al., 2009, 2011 ; Amoruso and Crescentini, 2013）上で示し

して，移動先の楕円体形状を変化させた場合の SM のた上で，2 つの楕円体が膨張する場合の MT の取り得るトータル MT の取り得る範囲を MTCR 図上に整理した

範囲を MTCR ダイアグラム上で整理した（Fig. 1）．これ

ところ，EX の MT と全く重複しない別の領域となる大を既往研究で慣習的に用いられる 3 つの単純幾何形状の変興味深い結果を得た（Fig. 1 の領域 PX'LDZB'X'P）．こ空隙（Fig. 1 の球 S : Sphere，無限長の円柱 C : Cylinder, 亀の高圧球体 SM により楕円体体積震源の MT の取り得る＊〒113-0033 東京都文京区本郷 7-3-1
東京大学大学院理学系研究科地球惑星科学専攻
Department of Earth and Planetary Science, Graduate
School of Science, The University of Tokyo, 7-3-1
Hongo, Bunkyo-ku, Tokyo 113-0033, Japan.
＊＊〒113-0032 東京都文京区弥生 1-1-1

東京大学地震研究所
Earthquake Research Institute, The University of Tokyo,
1-1-1 Yayoi, Bunkyo-ku, Tokyo 113-0032, Japan.
Corresponding author: Mie Ichiharae-mail : ichihara@eri.u-tokyo.ac.jp

450

水野尚人・亀伸樹・市原美恵
範囲が大幅に拡張された．任意の高圧楕円体形状の場合にもトータル MT の取り得る範囲を整理したところ，この 高 圧 球 体 SM の 範 囲 か ら の ず れ は 僅 か で あ っ た（Mizuno et al., 2015 の Fig. 3a 参照）．
しかしながら，観測データの MT は EX と SM から外れる範囲（Fig. 1 の領域 BXPX'B'CB）に入るものも報告されている．Fig. 1 に示されたキラウエア火山の地殻変動（Davis, 1986），浅間山の長周期振動（Maeda, 2009 ;Maeda and Takeo, 2011）の MT がその例である（Fig. 1 のK と A）．キラウアエア火山の MT は複数 EX の範囲（3-1 節の ΔSPC）には入っているが，浅間山の MT はいずれにも該当しない．Mizuno et al. (2015) では，浅間山の MT を説明するために高圧球体 SM を考えた後に更に圧力回復が起きるモデル（Pressure recovery : PR）を考えた．SM において高圧球体に起きる圧力低下に対してその何 % が回復するかを圧力回復率 Δp'% として与えた
場合，移動先楕円体の幾何形状と圧力回復率を用いてトータル MT が解析的に表現できる．PR は浅間山のMT の解釈に適用され，もっともらしい結果を得た．PRの取り得る範囲は，Fig. 1 の三角形領域 ΔSLZ となり EXと SM を完全に含む．また，この領域は MT の対角成分の和が非負の場合に取り得る MT 範囲でもある．
我々は，Mizuno et al. (2015) において MT の解析的表現が得られたこれらの 3 つのモデルに対して，誰でも簡便に MT を評価できるように Web 計算ツールを開発しFig. 1. MTCR diagram of the possible ranges of the
MTCR for (i) an ellipsoidal volume source, EX :SXPXBS, (ii) simple magma movement, SM : PX 'LDZB ' X ' P and (iii) pressure recovery, PR : Δ SLZ,
where S, B, and P indicate the MTCR values for the
expansion of a sphere, a prolate ellipsoid, and a pennyshaped crack, respectively, and C, D, and L are for theconventional cylinder model, the double couple,
CLVD, respectively. X' and B' are the total MTs in SMfrom S to X and B, respectively, and Z indicates thecorner position (1, −2). The aspect ratio of a1/a3 of anellipsoid is varied along the solid contour lines for aconstant a2/a3 as indicated by bold labels, while a2/a3
is varied along the dashed contour lines for a constant
a1/a3 as indicated by normal labels. The representativevalues of observed moment tensor for the 1970
deformation at Kilauea, Hawaii (Davis, 1986), and thatfor VLP observed at Asama, Japan (Maeda, 2009 ;Maeda and Takeo, 2011), are shown by K with a solid
triangle and by A with an open triangle, respectively.
This figure is modified from figure 3 of Mizuno et al.
(2015).

た．Mizuno et al. (2015) の解析的表現はいずれも，1 つの楕円体震源の単純膨張の解（Eshelby, 1957 ; Davis, 1986）に依拠し，楕円積分が含まれる．我々の Web 計算ツー
ルは，この楕円積分を意識することなく楕円体の形状を画面上にて入力するだけで MT 出力が得られる順計算，
MT 入力に対して形状が得られる逆計算の両方に対応している．なお，Web 計算ツールは，SM では高圧球体のみを対象とし（3-2 節），PR では圧力回復率を追加で入力する（3-3 節）．
この Web 計算ツールにより，地震波形データや地殻
変動データの逆解析で得られた MT を楕円体震源モデルへ対応させることができる．火山震源過程の推定に役立
てられることを願い，3 つの楕円体体積震源モデルに対応する Web 計算ツールの使い方をそれぞれ練習問題を交えながら説明する．
2．楕円体体積震源 Web 計算ツール : 共通設計
我々の Web 計算ツールは，http://www.eri.u-tokyo.ac.jp/
ichihara/mizuno/index.html にて公開されている．まず最初に Fig. 2（a）の画面が現れるので最上部のメニュー（Expansion（EX）, Simple Movement（SM）, Pressure Recovery楕円体体積震源モーメントテンソル評価のための Web 計算ツール

451

で（1, M22/M11, M33/M11）の形で出力される．すなわち，
MT の独立情報は M22/M11, M33/M11 の 2 つである．楕円体形状（a3≥a2≥a1）に対して一般的に M11≥M22≥M33となるが，MT 成分の大きさが逆転する部分がわずかに存在し（領域 XP と PX'）この場合，MTCR ダイアグラム≡（M33/M11, M22/M11）に取上の点は（M '22/M11, M '33/M11）
り直して，常に M '22/M11≥M '33/M11（ダイアグラム右下）に点をプロットする仕様になっている．
逆計算では（M11, M22, M33）の対角 3 成分を大きい順（M11 ≥ M22 ≥ M33）に 入 力 す る と，楕 円 体 の 幾 何 形 状（a2/a3, a1/a3）が出力される．上記の MT 成分の逆転領域
では，逆転が考慮され常に a2 /a3 ≥a1 /a3 で出力される．
順計算と逆計算ともに 2 つの数字の入力に対して 2 つの数字が出力される．なお，PR では圧力回復率 Δp'%（3-3節）を追加で指定する必要がある．
Web 計算ツールでは楕円体形状に依存する特徴的な
2 つの量 Ψ =ΔV C/ΔV Tと K C=k ellip/k が同時に表示される．ここで，ΔV Tは応力ゼロ歪みによる楕円体の体積変化，ΔV C は実際の体積変化，k は媒質の体積弾性率，（ΔV C/V）は楕円体空隙の実効的な体積弾性率，
k ellip≡ΔP/
ΔP は圧力変化，V は楕円体体積である．ΔV T は体積震源のモーメントテンソルに現れる仮想的体積変化であり，K C は楕円体震源の体積変化に対する実効的圧力変化の指標である．最近，体積震源の MT 表現において
ΔV C ではなくΔV T が現れる必然性について表現定理の観点から説明された（Ichihara et al., 2016 ; 日下部・他，
2016）．なお，EX ツールのみの付録機能として，計算結果の MT の MTCR ダイアグラム（M22/M11, M33/M11）上Fig. 2. Ellipsoidal volume source calculation webtool : (a)Display of the main menu, and (b) Display of ‘Simple
Movement (SM)’ webtool.

での位置情報と楕円体の断面形状が表示される．
3．各 Web 計算ツールの使用法と練習問題
3-1

（PR）
）の 3 つから計算したいモデルを 1 つ選ぶ．
全てのモデルで 1 つの楕円体空隙の幾何形状が入出力

Expansion（EX） : 楕円体の膨張モデル

1 つの楕円体空隙の膨張モデルによる MT 計算にはFig. （2 a）のメニュー画面から Expansion（EX）を選ぶ．EXに関係する．デカルト座標系（x1, x2, x3）を考え，楕円体の Web 計算ツール画面において楕円体形状（a2/a3, a1/a3）

は座標軸と平行な 3 軸（a1, a2, a3）を持つ．楕円体形状を 入 力 し て 下 向 き ↓ を ク リ ッ ク す る と EX の MTを入力して MT を出力することを順計算，MT を入力し（1, M22 /M11, M33 /M11）が計算される（順計算）．逆に，て楕円体形状を出力することを逆計算と呼ぶ．計算対象（M11, M22, M33）を入力して上向き↑をクリックすると楕とする MT は，EX では 1 つの楕円体の MT であり，SM

円体形状（a2/a3, a1/a3）が計算される（逆計算）．楕円体と PR では球体と楕円体のトータル MT である．

の幾何形状に対応する Ψ と K Cの値が表示される．

例として Fig. 2（b）に SM Web 計算ツールの画面を示

逆計算における入力成分の大小順序は M11≥M22≥M33

す．楕円体形状に関する独立情報は（a2 /a3, a1 /a3）のとする．これを満たさない場合は，↑ボタンが赤色とな

2 つであり，この形状比の値を順計算の際に入力する

りエラーであることを示す．また，MTCR ダイアグラム．すると，楕円体の置かれ（1≥a2 /a3 ≥a1 /a3 >0 とする）

上で SXPXBS 領域に入っていない MTCR を入力してた座標系での MT の対角成分（M11, M22, M33）が M11でも，↑ボタンが赤色となりエラーであることを示す．

規格化された比（moment tensor component ratio, MTCR）

Web 計算ツールを使わなくても，楕円体形状と MTCR

452

水野尚人・亀伸樹・市原美恵の関係は MTCR ダイアグラム上の等値線の値により示

成分の大小が入れ替わる領域 XP では，MT 成分を共通されており（Fig. 1），概数はダイアグラムのみを用いて設計に従い並べ替える．

読み取ることが出来る．領域 SXPXBS において，実線はMTCR ダイアグラム上での，トータル MT の位置は，

a2/a3=const，点線は a1/a3=const の等値線を表している．

それぞれの要素の MTCR ダイアグラム上の点（M22/M11,

例えば，（a2/a3, a1/a3）の値を決めたなら，それぞれの値に

M33 /M11）と（M '22 /M '11, M '33 /M '11）と を 結 ぶ 線 分 を，

近い等値線の交点を探す．そこが対象とする楕円体膨張

tot
tot
/M 11
につ
M '11 : M11 で内分した点になる．これは，M 22

震源の MTCR である．楕円体形状を示す a2 /a3 =const.といて，

a1 /a3 =const. の等値線が集まる領域 XP 付近は，モーメ
ントテンソルを楕円体震源と対応させる際に特に注意が必要となる領域である．ここでは，形状が大きく変化しても MTCR はほとんど変化しない．すなわち，MTCR

M 22+M ' 22
M tot
22
=
M 11+M ' 11
M tot
11
M 22
M ' 22
M 11
M ' 11=+,
M 11+M ' 11 M 11 M 11+M ' 11 M ' 11

（3）から形状を求めることが非常に困難な領域である．
3-1-1

tot
tot
であり，また M 33
/M 11についても明らかに同じことがExercise EX1
17

モ ー メ ン ト テ ン ソ ル 解 析 か ら M11 = 4. 00 × 10 ,言えるからである．さて，この分割比 M '11 : M11 が何に（単位 N m）の値が決定
M22=2.80×1017, M33=2.40×1017

依存しているかを Mizuno et al. (2015) に沿って考えてみされたとする．このとき，適切な Web 計算ツールを選よう．

び，楕円体形状と体積変化を求めてみよう．

楕円体 1 と楕円体 2 の MT の値に対して，MTCR ダイ

MTCR が（2. 80/4. 00, 2. 40/4. 00）=（0. 700, 0. 600）とアグラム上の位置が決まり，

MTCR ダイアグラム（Fig. 1）における領域 SXPXBS に
あ る の で EX Web 計 算 ツ ー ル を 選 び，M11 = 4. 00,

α=1+

M22=2.80, M33=2.40 を入力する．逆計算器のボタン↑

M ' 22 M ' 33
M 22 M 33+, α'=1++,
M 11 M 11
M ' 11 M ' 11（4）

Ψ =0.617, K =0.621 となる．楕円体の形状は EX Webの値がそれぞれ確定する．一方， ΔV T，ΔP，および，ΔV T 'と ΔP' の関係（Eq.（1），および，Ψ と K C の定義式）

計算ツールの数字が答となる．

からを ク リ ッ ク す る と a2 / a3 = 0. 459, a1 / a3 = 0. 294,
C

実際の体積変化 ΔV C を求めるには，仮想的体積変化
T

ΔV と MT の対角成分の総和の関係式
M ii=3kΔV T,（1）を用いる（Mizuno et al., 2015）．ここで k は媒質の体積弾性率である．これと Ψ（=ΔV C/ΔV T）の値を組み合わせ
ると，実際の体積変化（2）となる．k=20 GPa とすると Eq.（2）よりΔVC =（4.00+
97

3（3×20×10 ）〜0.946×10 m と2.80+2.40）
×10 ×0.617/
73

なる．これは東京ドームの体積 1.2×10 m のおよそ 3/4の関係式が得られる．
次に，上の議論を共通の過剰圧 ΔP=ΔP' を持つ 2 つM ' 11 V ' Ψ K C α,=
M 11
V Ψ 'K C' α'

（6）

あるいは，
V ' M ' 11 Ψ 'K C ' α',=
M 11 Ψ K C α'V

（7）

2 つの楕円体膨張と MTCR ダイアグラム上のとなる．これは，内分点の位置は，楕円体の形状だけで決定される量の比 Ψ '/Ψ，K C '/K C，α '/α に加えて，体積内分点の関係比 V '/V が重要になることが分かる（Mizuno et al., 2015）．

に相当する．
3-1-2（5）の楕円体のそれぞれの MT に適用すると，(M 11+M 22+M 33)
M ii
Ψ,
Ψ=
ΔV C=
3k
3k
17

3V ΔP
3kΔV T=,M 11
Ψ K C M 11
3kΔV T '3V ' ΔP'α'==,
M ' 11
Ψ 'K C' M ' 11

α=

主軸の向きの揃った 2 つの楕円体（楕円体 1 と楕円体 2）

V ' ≫ V と す る と，小 さ い 方 の 楕 円 体 が 薄 い ク ラ ッ ク

が，共通の過剰圧 ΔP を受けて膨張する場合を考えよう． （K C '≫K C）の場合を除き M11'≥M11となるので，内分点
楕円体 1 に関する量を上添え字無しで，楕円体 2 に関すは（M '22/M '11, M '33/M '11）の方に引っ張られる．

る量を上添え字 ' 付きで，両者のトータルを上添え字 tot

以上のことから，1 つの楕円体の膨張の MT は MTCR付きで表わす．それぞれの MT の成分（M11, M22, M33）

ダイアグラムにおいて領域 SXPXBS しか取り得ないが，と（M '11, M '22, M '33）は，同じ方向の成分とし，また，

2 つの楕円体の膨張を考えるとそれらの主軸の向きがM11 ≥M22 ≥M33，および M '11 ≥M '22 ≥M '33 とする．MT

揃っている場合には，三角形領域 ΔSPB の任意の値を取楕円体体積震源モーメントテンソル評価のための Web 計算ツール

453

り得る．また，軸の回転も含めると，三角形領域 ΔSPCの任意の値を取り得る（Mizuno et al., 2015）．ただし，個々の楕円体の形状については任意性があり，MTCR ダ
イアグラム上の 1 点の位置情報から形状を一意に決めることはできない．
3-1-3

Exercise EX2

Davis (1986) によって決められた Kilauea 火山の体積
震源の MTCR は K（0.691, 0.535）であった．これは，単一の楕円体の膨張では説明できない領域にある．これを，MTCR 値（1, 2/3）を持つ細長い楕円体 B と，領域
XP にある薄い楕円体の組み合わせで考えてみよう．
例えば，B（1, 2/3）と K（0.691, 0.535）を通る直線を延長すると（0.420, 0.4195）を通る（この付近は等値線が集まってくる領域であり，0.4195 の数値だけ有効数字が 4 桁に選んだ．例えば，これを 0.420 に丸めると全く異なる楕
円体形状になる）．この点を P1とする（Fig. 3a）．P1の楕
円 体 形 状 を EX Web 計 算 ツ ー ル を 用 い て 求 め る と，
a2/a3=0.562, a1/a3=0.141, Ψ =0.776, K C=0.289 となる．
また，α=1.84 である．一方，B に対しては，Ψ '=0.625,
K C ' =0.600，α '=2.67 と求められる．K は線分 BP1 を
0.533 : 0.467 に内分する点である（1×0.467+0.420×
0.533=0.691，2/3×0.467+0.4195×0.533=0.535）．すな
わ ち，B と P1 の 楕 円 体 の MT の 最 大 成 分 を そ れ ぞ れ
M '11 と M11 と す る と，M '11 /M11 = 0. 467/0. 533 で あ る．
Eq.（7）から，
V ' 0.467 0.625×0.600 2.67=
=2.13,
V
0.533 0.776×0.289 1.84

（8）と，2 つの楕円体の体積比が決められる．
0.711）と P（0.367,
0.350）
これに対して，K を B（1.00,
2
2
を 0.488 : 0.512 に内分する点と考えてみよう（Fig. 3b）．
P2 は，軸の長さが a2/a3=0.0833, a1/a3=0.0112 の細い棒状で，Ψ =0.895，K C=0.117，α=1.72 である．一方，B2は，軸の長さの比が a2/a3=a1/a3=0.310 の卵型で，Ψ '=0.597，K C '=0.675，α '=2.71 が得られる．これらの値から，
V ' 0.512 0.597×0.675 2.71=
=6.36,
V
0.488 0.895×0.117 1.72

Fig. 3. Supplemental plots on MTCR diagram for Exercise
EX2, determining Kilauea point K by EX : (a) Case with
an ellipsoid of B (1, 2/3) as a starting MTCR point, (b)Case with an ellipsoid of B2 (1.00, 0.711) as a starting
MTCR point. The corresponding goal points P1 and P2,
and the aspect ratio of each ellipsoid are also shown.
Contour lines and symbols are shown as similar to Fig. 1.

（9）

となる．
ここで取り上げた 2 つの場合のように，1 つの MTCR

Movement（SM）を選ぶ．SM の Web 計算ツールの画面（Fig. 2 (b）) が現れるので，楕円体形状（a2/a3, a1/a3）を入力して下向き↓をクリックすると SM モデルの（1, M22に対して，形状や体積比の大きく異なる組み合わせで観/M11, M33/M11）が計算される（順計算）．一方，（M11, M22,

測値が説明できることになる．

M33）を入力して上向き↑をクリックすると楕円体形状

Simple Movement（SM）: マグマの単純移動モデル（a2 / a3, a1 / a3）が 計 算 さ れ る（逆 計 算）．逆 計 算 で は，

高圧の球体から低圧の楕円体へのマグマ移動モデルにMTCR 入力値が SM の MT の取り得る範囲である領域

よる MT 計算には，Fig. 2（a）のメニュー画面から Simple

PX'LDZB'X'P（Fig. 1）に入っていなければ↑が赤色とな

3-2

454

水野尚人・亀伸樹・市原美恵

りエラーであることを示す．計算結果は，楕円体の Ψと K C の値と共に表示される．
EX モデルの場合と同様に，Web 計算ツールを使わなくても，Fig. 1 の MTCR ダイアグラム上に引かれた楕円
体形状の等値線により概数がわかるようになっている．
領域内部の実線と点線は，移動先の膨張する楕円体形状の a2/a3=const.と a1/a3=const. の等値線をそれぞれ表している．
（a2/a3, a1/a3）の値を決めたなら，それぞれの値に近い等値線の交点を探す．そこが対象とする SM のトータル MT の MTCR である．楕円体形状を示す a2/a3=const. の等値線が集まってくる領域 PX'付近は，モーメントテンソルから楕円体形状を決定することが困難な領域である．
3-2-1

マグマ単純移動モデル（SM）の図学的理解

SM モデルではマグマ移動により高圧の球体は収縮し，低圧の楕円体は膨張し圧力が釣り合う．以下，高圧の球体に関係する量を上添え字 ' 付きで，低圧の楕円体
は上添え字無しで表す．
SM モデルで注意すべきことは，球体の収縮量の絶対値（ ΔV C ' ）と楕円体の膨張の絶対値（ ΔV C ）は，必ずしも同じ値にならないことである．まず，このことについて考えてみよう．Mizuno et al. (2015) は，移動後の圧力平衡と，マグマの圧縮性を考慮した質量保存から，以下の関係を導いた．



kK C
+1
km
Ψ
ΔV
=A.
T =
C'
Ψ'
kK
ΔV
+1
km
T'



（10）

真の体積変化が  ΔV C '/ΔV C = ΔV T '/ΔV T  ×（Ψ '/Ψ）であることを思い起こすと，  ΔV C ' = ΔV C  となるのは，Eq.

Fig. 4. A supplemental plot on MTCR diagram for
Exercise PR, determining Kilauea point K by PR. M
and B3 indicate the MTCR values for an ellipsoidconsidered in SM and EX, respectively. Contour lines
and symbols are shown as similar to Fig. 1.

（10）において，km≫k のとき，又は，KC '=KC の場合であることが分かる．マグマと周辺岩石の体積弾性率を考えとなり，それぞれの体積とは関係なく，それぞれの幾何

ると，km≪k の方が適当であり，本 Web 計算ツールもこ形状のみで決まる．この性質を用いて，Fig. 1 においてれを仮定している．また，球体と同じ値の KC を持つ楕SM の MT を決定することができる．

円体形状は球体以外になく，球体から球体への流体の移

Mizuno et al. (2015) は，マグマの単純移動モデル（SM）

動は，MT 自体が符号を反転して等しくなり，重ね合わにおける球体と楕円体のトータル MT の MTCR 値を，せると 0 となる．すなわち点震源として捉えた場合にMTCR ダイアグラムから図学的に求める方法を提案しは，意味のないものになってしまう．

ている．球体膨張と楕円体膨張の MTCR ダイアグラム

次に，体積震源を特徴づける量である MT の対角成分上での位置をそれぞれ S，B3 であるとする．マグマの移の和（Mii）について考えてみよう．MT の各要素 Mij i=jを決めるのは，ΔV T 'とΔV T である．本 Web ツールでは，

動により球体の収縮と楕円体の膨張のトータル MT のT'

MTCR 上の位置（M）は，SM : MB3=1−xmov : xmov の外km≪k を仮定したが，この時 Eq.（10）からΔV とΔV の分点になる（Fig. 4）．xmov は負の値で以下の式で与えら比はれる（Mizuno et al., 2015）．

ΔV T'

KC Ψ

 ΔV = K ' Ψ ' =A,
T

C

T

（11）

楕円体体積震源モーメントテンソル評価のための Web 計算ツール

x mov=α3A movF, F=− =− ,movαα1+A F

A mov=Ψ KC
9
Ψ K C.
 =4
Ψ K CA=Ψ KC
9= K CΨ ,Ψ 'K C ' 4

ことに注意する．これを Web 計算ツールのの値を組み合わせると，

一般的な高圧楕円体形状に対しても Eq.（12）の 3 つの値

A=

9
9K CΨ
= ×0.563×0.777=0.984, (1−A)=0.016,
4
4
（16）

は全て形状だけで決まり，それぞれの容器の大きさによらない．また，移動によるトータル MT の Mii が与えられると，高圧側の収縮量  ΔV C '  と，低圧側の膨張量  Δ
V C  も形状のみから決まることは，Exercise SM1 に示すとおりである．

M ii=(4.00+0.620−2.80)×10 17=1.82×10 17,
ΔV T=

れぞれの容器の大きさの比が効いてくる．Mizuno et al.
(2015) は，質量保存と圧力平衡から関係式





,

（13）

（18）
ΔV T '=−AΔV T=−1.87×10 8 m3,（19）

ΔV C=ΨΔV T=0.563×1.90×10 8=1.07×10 8 m3,（20）

5
ΔV C '=ΨΔV T '= ×(−1.87)×10 8=−1.04×10 8 m3,
9（21）を得ている（Mizuno et al. (2015) の Eq.（37）を変形）．
km≪k の場合，薄い亀裂（KC が非常に小さい）を含まない

（17）

M ii
1.82×10 17=1.90×10 8 m3,=
3k(1−A) 60×10 9×0.016

一方で，高圧側と低圧側の圧力変化比については，そ1
1+k m kK C
1
1+
k m kK C '（15）となり（Mizuno et al., 2015），膨張と収縮で打ち消し合う（12）

ここで高圧球体の値 Ψ ′=5/9 と K C '=4/5 を用いたが，

 ΔP'
V
=
 ΔP
V'

455

となる．

限り，上の式の（ ）内は 1 で近似でき， ΔP' / ΔP 〜V/V 'Exercise EX1 で は，Δ V C = 0. 946 × 107 m3 で あ っ た．

となる．楕円体形状の変化に対して，こうして決定された外分点の位置を結んだものが Fig. 1 において MTCR

Exercise SM1 で は マ グ マ 移 動 に よ る 体 積 変 化（ΔV C =1.07×108 m3 と ΔV C '=−1.04×108 m3）が 10 倍ダイアグラムの SM 範囲に示された等値線である．

程度生じているにも関わらず，両者は打ち消し合い，トー

なお，CLVD（Fig. 1 の L）のように，Mii〜0 であるよ
うな MT は，流体移動を表わすと解釈されることが多い

タル MT における Mii =1.82×1017 N m は Exercise EX1の Mii=9.20×1017 N m より小さいことに注意する．

が，実際には同じ形状間，あるいは特定の形状のセット

3-2-3の間の流体移動でなければ，Mii〜0 を満たすことができ以上を踏まえて，注意すべき問題として，ダブルカッ

ないことをここに強調しておく．
3-2-2

プル（Fig. 1 の D）のモーメントテンソルを入力して逆Exercise SM1

解析を行ってみる．モーメントテンソル逆解析から M11

モ ー メ ン ト テ ン ソ ル 逆 解 析 か ら M11 = 4. 00 × 1017,
17

Exercise SM2 : Double couple

17=1.00×1017，M22=0.00×1017，M33=−1.00×1017（単位M22=0.620×10 , M33=−2.80×10 （単位 N m）の値がN m）の値が決定されたとする．このとき，適切な Web

決定されたとする．このとき，適切な Web 計算ツール

計算ツールを選び，楕円体形状と体積変化を求めてみよを選び，楕円体の形状と体積変化を求めてみよう．

う．

MTCR が（0.620/4.00, −2.80/4.00）=（0.155, −0.700）とMTCR が（0.00, −1.00）と MTCR ダイアグラム（Fig.

MTCR ダイアグラム（Fig. 1）における SM の範囲にある

1）における SM の範囲にあるので SM Web 計算ツール

ので SM Web 計算ツールを選ぶ．Web 計算ツールにを選び，M11=1.00, M22=0.00, M33=−1.00 を入力する．

M11=4.00, M22=0.620, M33=−2.80 を入力する．↑をク↑ボタンをクリックすると赤色となりエラーであること

リ ッ ク す る と a2 / a3 = 0. 796, a1 / a3 = 0. 694, Ψ = 0. 563,を示す．これは，球体収縮側と楕円体膨張側の MT を足K C=0.777 を得る．マグマ移動先の楕円体形状についてし合わせると 0 となる場合は，楕円体が球形状以外になはこの数字が答となる．

く，Mii〜0 となる場合は直線 LDZ 上で実現し得ないか

T'
T
ΔV （>0）とすると，仮想的体積変化を ΔV （<0）とM ii=3k(ΔV T '+ΔV T)=3kΔV T(1−A),

らである。SM ツールにおいては CLVD の MT を含み
Mii〜0 となる場合は全てエラーとなる．

トータル MT の対角成分の総和の関係式は，（14）

Web 計算ツールの対象外であるが，コイン状の薄いクラックから，それに直交する同じ形状のクラックに流体456

水野尚人・亀伸樹・市原美恵が移動する場合には，意味のあるダブルカップル解を得

0.535）を，PR モデルで説明することを試みる．この場ることができる（Mizuno et al., 2015）．例えば，高圧側（収合にも，不確定性は多分にあるが，Web 計算ツールの適縮側）の MT が（M '11, M '22, M '33）=（−1, −1/3, −1/3），

用範囲内で考える．即ち，高圧側の容器は球体で，その=（1/3, 1, 1/3）である
低圧側（膨張側）が（M11, M22, M33）

体積は低圧側の楕円体容器よりも十分に大きく，マグマ

ような場合を考えれば，自明であるだろう．

の体積弾性率 km は周囲弾性体の k よりも十分に小さい．

Pressure recovery（PR）: マグマ移動後に圧力回復

3-3

するモデル

Fig. 2（a）のメニュー画面において Pressure Recovery を選び，PR Web 計算ツールに M11 =1.00，M22 =0.691，

SM の後に圧力回復する PR モデルの MT 計算には，

M33=0.535 を入力する．圧力回復を Δp '=50.0 % と入力

Fig. 2（a）のメニュー画面から Pressure Recovery（PR）

すると，楕円体の形状が a2/a3=0.686，a1/a3=0.511 とな

ツールを選ぶ．PR の Web 計算ツールにおいては順計算

る．この形状の楕円体が膨張する場合の MTCR は，と逆計算の両方の場合に圧力回復率 Δp '% の値（3-3-1

B3（0.814, 0.721）であることが，EX Web 計算ツールに節）を追加で入力する必要がある．（a2/a3, a1/a3）と Δp '%よって計算される．さらに，この形状を SM Web 計算

を 入 力 し て 下 向 き ↓ を ク リ ッ ク す る と PR の MTCR

ツ ー ル に 入 れ て 計 算 す る と，行 き 先 の MTCR は，（1, M22 /M11, M33 /M11）が計算される（順計算）．逆に，M（0.0712, −0.396）となる．SMB3 の位置より外分比（M11, M22, M33）と Δp '% を入力して上向き↑をクリック

xmov =−3.99 であり，Eq.（12）の第 1 式と第 2 式より，

すると（a2/a3, a1/a3）が計算される（逆計算）．計算結果

Amov =0.9466 となる．また，B3 の楕円体に対して計算は，楕円体の Ψ と K C の値と共に表示される．PR ツー

さ れ た，Ψ = 0. 579，K C = 0. 727 を Eq.（12）の 第 3 式
9Amov= 4 Ψ K C に代入しても Amov が同じ値が得られ，計算が正しく行われていることが確かめられる．

ル の 対 象 と な る の は MTCR ダ イ ア グ ラ ム に お け る
Mii〜0 を除く三角形領域 SLZ 全体である．
3-3-1

K を，SB3 の 外 分 点 と し て 見 た と き の 分 割 比 は

マグマ移動後に圧力回復するモデル（PR）の

obs

x

図学的理解
SM の後（Fig. 4 の M），両方の容器（Fig. 4 の S と B3）が圧力平衡を保ちながら増圧する場合を考える．PR の=−0.661 である．対応する Aobs =0.471〜Amov /2 となっており，Eq.（22）より 50 % の圧力回復を与えたことと整合的な結果が得られている．

最終的なトータル MTCR 値を R とすると，R は，0 % の圧力回復の場合には M と同じ位置となる．100 % の圧

4．まとめ

力回復（高圧球体が元の圧力に戻る）の場合には，実質

我々は Mizuno et al. (2015) において考察された楕円体

的には低圧楕円体の膨張の寄与しか見られないので，B3空隙に関係する 3 つの体積震源モデル（膨張 : EX，マグの位置になる．より一般的に，0 から 100 % の圧力回復

マの単純移動 : SM，単純移動後の圧力回復 : PR）の MTに対して R は直線 MB3 上を M から B3 に移動することを評価する Web 計算ツールを作成した．楕円体の形状が示される（Mizuno et al., 2015）．100 % を超えると，直を入力すると MTCR 出力を得られる順計算，MTCR を線 B3S 上において S に近づく．

入力すると楕円体形状が出力される逆計算がそれぞれ可Mizuno et al. (2015) は，移動の場合と同じように，外分能である（PR では圧力回復率を追加で入力する）．結果

比 SR : RB3 = 1 − xobs : xobs と，Aobs の 関 係 に つ い て，は Ψ と K C と共に画面に表示される．Web 計算ツール

Eq.
（12）の，xmov と Amov の関係が成り立つことを利用し，は誰にでも簡単に使ってもらえることを目指した．地震波形データや地殻変動データの逆解析から得られた MT

ΔP rec
+1
A obs
ΔP'
,
mov =
rec
A
ΔP
+1
ΔP

（22）

に対して，この Web 計算ツールを適用して火山震源過程の解釈に役立てて欲しい．
Web 計算ツールのプログラムは著作権フリーである．

の関係を導いた．ここで，ΔP rec は圧力回復値であり，Web 計算ツールトップページ（Fig. 2）を参照されたい．

Δp '=ΔP rec/ ΔP '  を百分率で表したものを Δp 'と定義す組 み 込 ま れ て い る 楕 円 積 分 に 関 し て は Mizuno et al.

る．高圧の球体の体積が，低圧の楕円体の体積よりも十分(2015) の Appendix に詳しく説明されている．興味があに大きい時（V/V '→ 0），Eq.（13）より  ΔP

rec

< ΔP ' ≪ΔP

る方はそちらも参照してほしい．

であり，Eq.（22）の右辺分母は，1 で近似できる．PR
Web 計算ツールでは，この近似を使うことにする．
3-3-2

Exercise PR

Exercise EX2 と 同 じ Kilauea の MTCR の K（0. 691,

謝辞

日下部哲也氏との議論は大変有用でした．本稿の改訂にあたり匿名の査読者 1 名，および編集委員の橋本武志

楕円体体積震源モーメントテンソル評価のための Web 計算ツール
氏のコメントは非常に有意義でした．記して感謝いたします．N. K. は JSPS 科研費 #21107007，および文部科学省による「災害の軽減に貢献するための地震火山観測研究計画」の助成を受けました．
引用文献
Aki, K. and Richards, P. G. (1980) Quantitative Seismology
1st ed. Freeman, New York, 57-60 p.
Aki, K. and Richards, P. G. (2002) Quantitative Seismology
2nd ed. University Science Books, Sausalito, 53-58 p.
安芸敬一・リチャーズ，P. G.（上西幸司・亀 伸樹・青地
秀雄 訳）(2004) 地震学 定量的アプローチ．古今書院，
東京，69-76 p.
Amoruso, A. and Crescentini, L. (2013) Analytical models of
volcanic ellipsoidal expansion sources. Ann. Geophys., 56,
doi : 10.4401/ag-6441.
Davis, P. M. (1986) Surface deformation due to inflation of anoriented triaxial ellipsoidal cavity in an elastic half-space,with reference to Kilauea volcano, Hawaii. J. Geophys.
Res., 91, 7429-7438.
Eshelby, J. D. (1957) The determination of the elastic field of
an ellipsoidal inclusion and related problems. Proc. R. Soc.
London, Ser. A, 241, 376-396.
Ichihara, M., Kusakabe, T., Kame, N. and Kumagai, H. (2016)
On volume source representations based on the representa-

457

tion theorem. Earth Planets Space, doi 10. 1186/s40623016-0387-3.
日下部哲也・亀 伸樹・市原美恵・熊谷博之 (2016) 表現定
理とグリーン関数 (2) ─体積震源のモーメントテンソ
ル表現─．地震，68，169-176，doi : 10.4294/zisin.68.169.
Maeda, Y. (2009) Very-long-period pulses at Asama volcano
inferred from dense seismic observation. Ph.D. thesis, The
University of Tokyo.
Maeda, Y. and Takeo，M. (2011) Very-long-period pulses at
Asama volcano, central Japan, inferred from dense seismic
observations. Geophys. J. Int., 185, 265-282.
Mizuno, N., Ichihara, M. and Kame, N. (2015) Moment
tensors associated with the expansion and movement of
fluid in ellipsoidal cavities. J. Geophys. Res., doi : 10.1002/
2015JB012084.
Trasatti, E., Cianetti, S., Giunchi, C., Bonafede, M., Agostinetti,
N. P., Casu, F. and Manzo, M. (2009) Bayesian source
inference of the 1993-1997 deformation at mount Etna
(Italy) by numerical solutions. Geophys. J. Int., 177, 806814.
Trasatti, E., Bonafede, M., Ferrari, C., Giunchi, C. and Berrino,
G. (2011) On deformation sources in volcanic areas :modeling the Campi Flegrei (Italy) 1982-84 unrest. Earth
Planet. Sci. Lett., 306, 175-185.
（編集担当

橋本武志）