インバージョンで仮定した断層面の位置を図-1に示す.断層面は,気象庁による前震の震源(北緯33.742°,東経130.808°,深さ11km)を含むように設定した.走向と傾斜は,F-netによる前震のメカニズム解(表-1)と一致するように212°と89°とした.長さと幅については,今回はすでにpublishされている既往研究に合わせることを考え,Asano and Iwata (2016)と同様,長さ13km,幅12kmとした.
インバージョンはHartzell and Heaton(1983)の方法に基づいている.13km×12kmの断層を13×12の小断層に分割し,それぞれの小断層でのモーメントレート関数は,小地震のモーメントレート関数とインパルス列との合積で表されると仮定した.インパルス列は0.25秒間隔の12のインパルスからなるものとし,このインパルスの高さをインバージョンの未知数とした.したがって破壊フロント通過後の3.0秒間だけすべると仮定したことになる.破壊フロントは,気象庁の破壊開始点から同心円状に拡大するものとした.その拡大速度については,2.1km~3.0kmの範囲で0.1km刻みで値を変えて解析を試みたところ,2.1km/sのときに最も残差が小さくなったため,最終的に2.1km/sの場合の結果を採用した.基盤のS波速度は3.55km/sとした.インバージョンには非負の最小自乗解を求めるためのサブルーチン(Lawson and Hanson,1974)を用いた.また,すべりの時空間分布を滑らかにするための拘束条件を設けた.観測波と合成波を比較する際には記録のヘッダに記載された絶対時刻の情報を用いている.
インバージョンに用いた観測点(図-1の▲)における観測波(黒または灰色)と合成波(赤)の比較(0.2-2Hzの速度波形)を図-2に示す.これらの図において観測波を黒で示した部分がインバージョンに用いた区間(10秒間)である.多くの地点で観測波と合成波は良く一致している.
|