経験的グリーン関数を用いた波形インバージョンにより2009年8月11日駿河湾の地震（Mj6.5）の破壊過程を推定した．グリーン関数としては，本震波形と余震波形の群遅延時間の類似性を考慮し，表-1に示す8月13日12:42の余震の記録を用いた．震源を取り囲むように存在するK-NETとKiK-netの15地点を選び（図-1），そこでの0.2-2Hzの速度波形をインバージョンのターゲットとした．表層地盤の非線形挙動の影響を避けるため，K-NET観測点では表層20mの平均S波速度が300m/sを下回る観測点を避けるようにした．またKiK-net観測点では地中の記録を用いた．

図-1　インバージョンで仮定した断層面と観測点の位置．

■は本震の破壊開始点（気象庁），
□は解析に用いた余震の震央（気象庁），
ドットは本震発生後24時間の余震分布（気象庁），
▲はインバージョンに用いた観測点をそれぞれ示す．
×は東名高速盛土のり面崩落箇所

気象庁の震源位置にF-netのメカニズム解を置いてNS成分とEW成分の理論的なラディエーション係数を計算してみると，図-2に示すようになり，EW成分はラディエーションの節に近い観測点が多いことがわかる．ラディエーションの節に近い観測点では，本震および余震のメカニズムのわずかな変化に応じてラディエーション係数が鋭敏に変化し，解析が不安定になりやすい．そこで，ここではNS成分の速度波形を解析に用いることにした．解析には本震波形のS波を含む10秒間を用いた．

図-2　気象庁の震源位置にF-netのCMT解を置いた場合の理論的なラディエーション係数．
（左）EW成分．（右）NS成分．

インバージョンで仮定した断層面の位置を図-1に示す．この断層面は，気象庁の震源（北緯34.785°，東経138.498°，深さ23km）を含むように設定し，走向と傾斜は，F-netのCMT解の2つの節面のうち，気象庁の余震分布（図-1のドット）とより調和的な走向307°，傾斜47°を採用した．断層面の長さと幅は気象庁の余震分布を参考に長さ20km，幅15kmとした．
インバージョンはHartzell and Heaton（1983）の方法に基づいている．20km×15kmの断層を20×15の小断層に分割し，それぞれの小断層でのモーメントレート関数は，余震のモーメントレート関数とインパルス列との合積で表されると仮定した．インパルス列は0.125秒間隔の16のインパルスからなるものとし，このインパルスの高さをインバージョンの未知数とした．破壊フロントは気象庁の震源から同心円状に速度2.9km/sで広がるものとした．基盤のS波速度は中央防災会議の解析を参考に3.82km/sとした．インバージョンには非負の最小自乗解を求めるためのサブルーチン（Lawson and Hanson，1974）を用いた．また，すべりの時空間分布を滑らかにするための拘束条件を設けた．観測波と合成波を比較する際には記録のヘッダに記載された絶対時刻の情報を用いている．

図-3にインバージョンの結果として得られた最終すべり量の分布を示す（S波速度3.82km/s，密度2.8ton/m3ですべり量に換算）． 同図に示すように，破壊開始点付近（a）およびその北西側（b）にアスペリティを有する震源モデルが得られた．余震のモーメントとして表-1の値を用いると，図-3に示す本震の最終すべり量の分布はMW=6.4に相当する．最大すべり量は0.92mである．

インバージョンに用いた観測点における観測波と合成波の比較を図-4に示す．これらの図において，横棒の部分がインバージョンに用いた部分である．この比較において，特に，震源から見て北東の方角にあたるSZOH43（清水南），SZO012（蒲原），SZO008（沼津）において揺れ始めの位相がうまく再現されていないことに気付く．それ以外の点では，波形の再現結果は概ね良好である．
上記の点に関して，防災科学技術研究所ではDD法による詳細な余震分布を求めており，本解析で対象とした北東傾斜の断層面の南側に，もう一枚，南傾斜の断層面が存在していたことを報告している（http://www.kyoshin.bosai.go.jp/kyoshin/topics/surugawan_20090811/inversion/）．本解析でうまく再現されていない北東側の観測点の初期位相は，南傾斜の断層面の破壊によるとも考えられる．その場合，本解析は，F-netのCMT解に表れた主破壊のみを対象としていることになる．
上記のような理由で，ここに得られた震源モデルは現時点で完全なものであるとは思えないが，特に被害の大きかった牧ノ原市を含む震源の南西側の観測点における0.2-2Hzの帯域の揺れはこの震源モデルでほぼ再現できる．図-3に示す震源モデルの小断層毎，時間ウインドウ毎のモーメント解放量（その要素に割り当てられた余震のモーメントに対する比）をテキストファイルに示す．

謝辞：本研究では（独）防災科学技術研究所のK-NET，KiK-netの強震記録，F-netのCMT解，気象庁の震源データを使用しています．ここに記して謝意を表します．

（グラフをクリックすると拡大します）

SZO004	SZO006
SZO008	SZO012
SZO014	SZO017
SZO018	SZO027
SZOH33	SZOH34
SZOH36	SZOH42
SZOH43

図-4　観測波と合成波の比較（黒が観測波）