３－１３ 近地強震記録による平成 20 年（2008 年）岩手・宮城内陸地震の震源過程
Rupture process of the 2008 Iwate-Miyagi Nairiku earthquake derived fromnear-source strong-motion records
防災科学技術研究所
National Research Institute for Earth Science and Disaster Prevention
2008 年 6 月 14 日 8 時 43 分に発生した平成 20 年（2008 年）岩手・宮城内陸地震（北緯 39.0度，東経 140.9 度，深さ 8km; 気象庁）について，K-NET，KiK-net の断層近傍の強震波形記録を用いて，震源過程のインバージョンを行った．
◎断層面モデルと震源過程のパラメータ化
第１図にHi-netの再検測による本震後１日間の余震分布とF-netのモーメントテンソル逆解析から推定された震源メカニズム及び Hi-net の P 波の押し引き分布による解を示した．得られている
余震分布の走向は北東-南西で，F-net のモーメントテンソル解の走向と調和的である．インバージ
ョンに使用する断層面モデルは，破壊開始点を DD 法により再決定された位置 1)，走向を F-net のモーメントテンソル解より 209 度，傾斜は地表変状の出現位置 2)を考慮して 40 度とし，大きさは余震分布の広がりを参考に長さ 40 km，幅 18 km とした（第１図に示した長方形）．
インバージョン解析において断層面上のすべり破壊過程は，時間・空間的に離散化して表現されている 3)．空間的には 2 km 四方の小断層 180 個（20×9）に分けた．時間的には各小断層において破壊開始点から一定速度で広がる同心円が到達してから時間幅 0.8 秒のスムーズドランプ関数を0.4 秒間隔で 7 つ並べることによって表現した．各小断層からの理論地震波形は，観測点ごとに異なる地下構造を考慮した１次元成層構造モデルを仮定して，離散化波数法 4)と反射透過係数法 5)により点震源の波形を計算し，これに小断層内部の破壊伝播の効果を付加した 6)．
◎解析に用いた波形データ
防災科研 K-NET および KiK-net（YMTH10 を除き地中）の 14 観測点（第１図）で得られた加速度強震波形に，0.1 から 1.0Hz のバンドパスフィルターをかけ，積分することにより得られた速度波形の S 波到達 1 秒前から観測点ごとに異なる時間長を切り出し（9 秒から 16 秒）データとした．IWTH25 が断層面直上に位置する唯一の観測点であることを考慮し，その重みを４倍とした．
◎波形インバージョン
各小断層の各タイムウィンドウのすべり量は，観測記録と理論波形の差の最小二乗法により解いた．インバージョンには，すべりの方向をモーメントテンソル解のメカニズムのすべり方向 104 度から片側 45 度の幅の中に納める拘束条件 7)と，時間的・空間的に近接したすべりを平滑化する拘束条件 8)をかけている．平滑化の強さは ABIC9)を参考とし，第一タイムウィンドウをトリガーする同心円の伝播速度は残差を最小とするものを選んだ．
◎結果
第２図に推定されたすべり分布を，第３図に観測波形と合成波形の比較を示す．断層面全体での㪄㩷㪈㪌㪏㩷㪄地震モーメントは 2.83×1019 Nm (Mw 6.9)である．破壊は，主として破壊開始点の南側の浅い方に進展し，破壊開始点の南約 8 km の領域で最大すべり量 5.7 m が推定された．第一タイムウィンド
ウをトリガーする同心円の伝播速度は 1800 m/s である．
注：なお，本解析は暫定的なものであり，今後修正される可能性がある．
（文責：鈴木亘・青井真（防災科研）
・関口春子（京大防災研・防災科研））
http://www.k-net.bosai.go.jp/k-net/topics/Iwatemiyaginairiku_080614/inversion/

引用文献
1) 関根秀太郎・汐見勝彦・武田哲也・浅野陽一・小原一成 (2008), DD 法による 2008 年岩手・宮城内陸地震震源域の地震活動 ―複雑な余震分布と直上観測点の意義―, 2008 年日本地震学会秋季大会, X1-017.
2) 産業技術総合研究所・活断層研究センター (2008), 2008 年岩手・宮城内陸地震速報,
http://unit.aist.go.jp/actfault/katsudo/jishin/iwate_miyagi/.
3) Hartzell, S. H. and T. H. Heaton (1983), Inversion of strong ground motion and teleseismic
waveform data for the fault rupture history of the 1979 Imperial Valley, California, earthquake,
Bull. Seism. Soc. Am., 73, 1553-1583.
4) Bouchon, M. (1981), A simple method to calculate Green's function for elastic layered media,
Bull. Seism. Soc. Am., 71, 959-971.
5) Kennett, B. L. and N. J. Kerry (1979), Seismic waves in a stratified half space, Geophys. J. R.
astr. Soc., 57, 557-583.
6) Sekiguchi, H., K. Irikura, and T. Iwata. (2002). Source inversion for estimating continuous
slip distribution on the fault, —Introduction of Green's functions convolved with a correction
function to give moving dislocation effects in subfaults—, Geophys. J. Int., 150, 377-391, 2002.
7) Lawson, C. L., and R. J. Hanson (1974). Solving Least Squares Problems, Prentice-Hall, Inc.,New Jersey, 340 pp.
8) Sekiguchi, H., K. Irikura, and T. Iwata (2000), Fault geometry at the rupture termination of the1995 Hyogo-ken Nanbu earthquake, Bull. Seism. Soc. Am., 90, 117–133.
9) Akaike, H. (1980), Likelihood and the Bayes procedure, in Bayesian statics, edited by J. M.
Bernardo, M. H. DeGroot, D. V. Lindlely, and A. F. M. Smith, University Press, Valencia, Spain.

㪄㩷㪈㪌㪐㩷㪄

╙㧝࿑ ࠗࡦࡃ࡯࡚ࠫࡦ⸃ᨆߦ↪޿ߚᢿጀ㕙ߣ᷹ⷰὐಽᏓ㧚ᤊශ㧘㕍ਣߪߘࠇߙࠇᧄ㔡෸߮㧝ᣣ㑆ߦ⊒↢ߒߚ૛㔡ߩ㔡ᄩࠍ␜ߔ㧚㤥޿චሼߪ↥✚⎇ߦࠃࠅႎ๔ߐࠇߡ޿ࠆ࿾⴫ᄌ⁁ߩ૏⟎ࠍ␜ߔ㧚Hi-net ߦࠃࠆ P ᵄೋേ⸃ߣ F-net ߦࠃࠆࡕ࡯ࡔࡦ࠻࠹ࡦ࠰࡞⸃ࠍฝ਄ߦ␜ߔ㧚Fig㧚1 Fault plane model and slip distribution used for the inversion analysis. Yellow star and blue circles
indicate the epicenters of the mainshock and aftershock within the day, respectively. Black crosses
indicate the location of the surface rupture reported by AFRC/AIST. Source mechanisms determinedby the P-wave polarity analysis using Hi-net data and moment tensor inversion using F-net data arealso shown.

╙㧞࿑ ᢿጀ㕙਄ߩߔߴࠅಽᏓ࿑㧚ᤊශߪ⎕უ㐿ᆎὐࠍ␜ߔ㧚ߔߴࠅߩࠦࡦ࠲࡯ߪ 1m 㑆㓒ߢឬ߆ࠇߡ޿ࠆ㧚Fig㧚2 Slip distribution on the assumed fault plane. A star indicates the rupture starting point.
Contours are drawn at 1 m interval.

㪄㩷㪈㪍㪇㩷㪄╙㧟࿑‫ⷰޓ‬᷹ᵄᒻ㧔㤥㧕ߣวᚑᵄᒻ㧔⿒㧕ߩᲧセ㧚ฝ਄ߩᢙሼߪߘࠇߙࠇߩᦨᄢ୯㧔m/s㧕ࠍ␜ߔ㧚Fig㧚3 Comparison between the observed and synthetic waveforms. The maximum velocities
for the observation and simulation are shown upper right of each trace.

㪄㩷㪈㪍㪈㩷㪄