ＹＣＵ地震学レポートNo.４４ Aug.19,95
８月１６日 ソロモン諸島地震（本震，余震） と
８月１７日 中国の核実験による地震波の解析
１。ソロモン諸島の地震
日本時間で８月１６日午後７時３０分ごろ，ソロモン諸島を震源とするＭs7.8の地
震がありました。Ｍs の値は先月のチリ北部地震と同じです。さらに半日後ぐらいに
同じ場所でＭs7.2の地震がありました。これは時系列的には余震ですが，Ｍs7.8の余
震としては異常に大きいものです。一体この地震の破壊域が本震の破壊域に含まれて
いるかどうか興味があるところです。
アメリカ地質調査所（ＵＳＧＳ）による震源情報は次の通りです。
発震時刻 震央 深さ Ｍｓ
------------------------------------------------------
本震 95/08/16 10:27:28 UT 5.7S 154.1E 33 km 7.8
余震 95/08/16 23:10:28 UT 5.7S 154.1E 71 km 7.2
ここで，ＩＲＩＳの広帯域遠地実体波を次の手順で解析しました。まず１個の点震
源を仮定して，そのメカニズムと深さを最小自乗法で決める。次いで，得られた最適
メカニズムを固定し，反復はぎとり法によって，１つの断層面上でのサブイベントの
時空分布を求める。
結果を図１，２に示します。本震，余震とも破壊パターンがよく似ています。はじ
めに小破壊があり，少し小休止のあと主破壊へと進んでいます。大地震の典型的な破
壊パターンです。問題は破壊域ですが，本震では北西方向にほぼ 120〜140kmにわたる
片方向破壊伝播が得られたのに対し，余震では逆に南東方向に60kmほどの破壊伝播モ
デルが得られました。いずれ余震分布で破壊域をチェックする必要がありますが，も
しこの結果が正しいとすると，本震と余震の破壊域はオーバーラップしていないこと
になります。
震源の深さについては，ここで得られた結果とハーバード解，ＵＳＧＳ解の間にか
なりのばらつきがあります。主な震源パラメターを表１に示します。比較のため，ハ
ーバード大学のＣＭＴ解も示します。
応力降下はプレート間大地震の平均値よりかなり低目です。この地域では１９７１
年７月に２つの巨大地震（Mw= 8.0, 8.1）が起こっていますが，これらの地震も低目
の応力降下をもっています（Lay & Kanamori,1980,PEPI)。
表１
--------------------------------------------------------------------------
Ｈ Ｍｏ Ｍｗ メカニズム Ｔ Ｓ Δσ
km 1020 Nm (strike,dip,rake) s km2 MPa
--------------------------------------------------------------------------
本震 60 4.73 7.7 (335,40,106)/(135,52,77) 60 140x60 1.5
(HRV CMT) 45 5.5 7.8 (331,57,107)/(121,37,65) 37
--------------------------------------------------------------------------
余震 40 0.52 7.1 (324,49,104)/(123,44,74) 27 60x40 1.1
(HRV CMT) 44 0.73 7.2 (310,42, 88)/(133,48,92) 16
--------------------------------------------------------------------------
Ｈ（初期破壊の深さ），Ｍｏ（モーメント），Ｔ（破壊継続時間），
Ｓ（断層面積），Δσ（応力降下）
２。中国の核実験による地震波
８月１７日昼ごろ，中国での核実験のニュースが伝わってきました。ＵＳＧＳによ
る震源は次の通りです。
発震時刻 震央 深さ ｍb
95/08/17 00:59:57.8 UT 41.6N 88.7E 0 km 6.1
”何ゆえにいま核実験を”と残念な気持ちを持ちつつ，広帯域地震観測点にダイヤル
アップしてみました。マグニチュード６相当ならはっきりと地震波がとらえられるも
のと思ったのですが，実際には，松代や中伊豆のような静かな観測点でも，そのまま
では脈動に隠れて見えません。図３は0.8Hz より低周波をカットした変位記録です。
振幅0.1〜0.2 ミクロンのＰ波とＰｃＰ波が見えます。残念ながら当方の観測点ＹＣＵ
などではノイズレベルが高いためほとんど信号が特定できません。
次に，ＩＲＩＳの広帯域遠地実体波記録のインバージョンを試みました。変位記録
を図４に示します。脈動成分を除くために， 0.3〜0.8 Hz より低周波成分をカットし
ました。震央距離は 35〜90°の範囲です。ＢＧＩＯの記録は大変きれいなのですが，
他の記録と比べて振幅が１０倍もあるので，インバージョンには使いませんでした。
震源がどのような時間関数をもっているのか全く予備知識がないのですが，図４の
波形を見るかぎり，通常の自然地震と同じような時間変化をしていますので，三角形
型の震源時間関数を用いました。また，モーメントテンソルの基底として，等方成分
を含む６成分全部を用いました。
結果を図５に示します。予想どおり，等方成分Ｉが大部分（５０％以上）のメカニ
ズム解が得られました。モーメントテンソルの全成分は次の通りです。
Ｍij ＝ 10.3 2.6 -1.5 [x1016 Nm]
2.6 1.1 -1.5
-1.5 -1.5 2.5 (x1,x2,x3)=(North,East,Down)
主値Ｍ1,2,3，スカラーモーメントＭｏ，Ｉの値は
Ｍ1,2,3 ＝ (11.2, 2.7, 0.0) x1016 Nm
Ｍｏ＝（Ｍ1＋Ｍ3）/2 ＝ 5.6 x1016 Nm
Ｉ＝（Ｍ1＋Ｍ2＋Ｍ3）/3＝ 4.6 x1016 Nm
です。また，震源時間はＴ＝0.1s（短い！）です。もちろん，これらはＴ*（T/Q）の
値にかなり依存します。ここではＴ*＝0.7sとしています。最適の深さは ｈ＝1.2 km
です。
仮に，震源をｒ=1kmの球とすると，体積増加に対応した応力変化は
Δσ＝Ｉ／Ｖ＝10 MPa
となります。
以上が主な結果ですが，これまでこの種の震源を扱ったことがないので，得られた
解が具体的にどのような物理量を表しているのかよくわかりません。
なお，等方成分無しの束縛条件を加えた場合の解を図６に示します。あきらかに波
形の一致は悪くなります。したがって，等方成分の存在はかなりはっきりしていると
言えます。
表２ 震源パラメタ
Ｈ Ｍｏ（Ｉ） Ｍw Ｔ Ｖ Δσ
----------------------------------------------------------------
1.2 km 5.6 (4.6) x1016 Nm 5.1 0.1 s 4 x109 m3 12 MPa