３－１７

平成 20 年（2008 年）岩手・宮城内陸地震について
On the Iwate-Miyagi Nairiku earthquake in 2008
2008 年岩手・宮城内陸地震緊急観測グループ／東北大学大学院理学研究科

Group for the Emergent Observation of the Iwate-Miyagi Nairiku Earthquake in 2008,
and Graduate School of Science, Tohoku University
2008 年 6 月 14 日 8 時 43 分に M7.2 の岩手・宮城内陸地震が発生した．この地震の直後から各研究機関による余震・GPS・電磁気に関する合同観測が実施された．
余震観測では北海道大学・弘前大学・東北大学・東京大学地震研究所・名古屋大学・京都大学防災研究所・高知大学・九州大学・鹿児島大学・防災科学技術研究所の 10 機関が合同で臨時余震観測を行った．その観測点配置を第 1 図に示す．この合同余震観測によって得られたデータを元にDD トモグラフィが実施され，高精度の地震波速度構造と震源分布が得られた．求められた震源分
布を第 2 図に，また得られた S 速度構造偏差を第 3 図に示す．余震は高速度域に集中して発生している．火山の下に高温域と考えられる低速度域があり，このためにこの領域では余震が発生していないと考えられる．余震は基本的に西に傾斜した面状に分布していると考えられるが，火山の下で余震が発生していないために，その傾向はあまり明瞭ではない．
深さ 24km における S 波速度偏差と再決定された過去の地震活動 (岡田・他，2008)1) とを比較したのが第 4 図である．この図を見ると，この領域には火山フロントに沿う南北方向の地震活動帯と，そのすぐ東の北北西-南南東方向に分布した活動帯があり，今回の地震は，その二つの活動帯があたかも収束するように見える領域付近で発生した．２つの地震帯の深部には地震波速度低速度域が見られ，それは地下深部から供給される流体の上昇経路であると考えられる．1962 年・2003 年宮城県北部地震や，1970 年秋田県南東部地震，1996 年の鬼首付近の地震の震源域の深部延長にも低速度域が分布しており，これらの地震の発生も深部から供給された流体と関係していると考えられる．
今回の地震の前駆的活動を第 5 図に示す．本震発生直前の 08:01 と 08:11 に発生した前震と 08:43に発生した本震の震源位置はほとんど同じである．また，1999 年と 2000 年にそれぞれ M4.3 とM4.9 を最大とする地震活動が今回の本震の震源の東側浅部に発生しており，その活動域では今回の地震の余震活動は低調である．つまり，今回の地震断層の浅部延長上で 1999 年と 2000 年に活動があって，それが今回の地震の発生に繋がったと考えられる．
震源域直上の一関市厳美町字祭畤（まつるべ）に東北大学が設置した GPS 連続観測点 ICNS（位置は震源の南西 2.5km 程度）における地震発生前後の変動および GPS 観測点配置図を第 6 図に示す．地震波によって得られたメカニズム解から今回の地震は東西圧縮の逆断層型であったことが明らかになっており，ICNS 観測点が北東に隆起したことは断層面が西に傾斜していることを示している．GPS 観測データを用いて 2 枚の矩形断層を仮定した場合の地震断層モデルを第 7 図に示す．
推定されたモーメントマグニチュードは２枚の矩形断層をあわせて 6.9 となった．この地震のあと
GPS 大学連合によって設置された臨時点も加えて Yagi and Kikuchi (2003)2) の時間依存インバージョン法によって推定された余効すべりの時空間分布を第 8 図に示す．余効すべりは基本的に断層の上端付近の浅部と栗駒山直下付近に推定されている．7 月 12 日までの解析で，余効すべりによるモーメント解放量はモーメントマグニチュード換算で北側断層が 5.9，南側が 6.3 と推定された．
この震源域周辺では広帯域 MT 法探査によって比抵抗構造が推定されていた．３本の測線にそっ㪄㩷㪈㪎㪉㩷㪄た２次元比抵抗断面に余震の震源分布を重ねた図を第 9 図に示す．基本的に高比抵抗域に余震が分布していることがわかる．この地震のあと，地震前後の比抵抗変化を調べるために，東北大学・秋田大学・北海道大学・東京大学地震研究所による合同 MT 観測が実施されている．
謝辞
余震の震源決定と地震波速度構造推定には，防災科学技術研究所 Hi-net・気象庁・国立天文台水沢VERA 観測所のデータ，東京工業大学が実施した臨時余震観測データ，および原子力安全基盤機構（JNES）が平成 19・20 年度に実施した内陸の活断層調査に基づく震源断層評価手法の検討事業で取得されたデータも使用させていただいた．GPS データ解析においては国土地理院 GEONET・国立天文台水沢 VERA 観測所のデータ，および原子力安全基盤機構（JNES）の上記事業で取得された GPS データも使用させていただいた．GPS 解析ソフトとしては，GpsTools ver.0.6.3 (高須・笠井，2005)3) を利用させていただいた．さらに観測にあたっては，一関市，一関市教育委員会，栗原市はじめ地元の関係諸機関・自治体に多大なるご協力をいただいた．記して感謝申し上げる．
参考文献
1) 岡田・他，東北地方中央部およびその周辺の地殻内地震波速度不均質構造と地震活動，日本地球惑星科学連合 2008 年大会予稿集，S147-008，2008．
2) Yagi and Kikuchi, Partitioning between seismogenic and aseismic slip as highlighted fromslow slip events in Hyuga-nada, Japan, GRL, 30, doi:10.1029/2002GL015664, 2003.
3) 高須・笠井, 測位衛星用高精度軌道・時計決定ソフトウェアの開発, 第 49 回宇宙科学技術連合
講演会, 広島, 2005．

㪄㩷㪈㪎㪊㩷㪄

㻌

╙1࿑วห૛㔡᷹ⷰࠣ࡞࡯ࡊߦࠃࠆ⥃ᤨ૛㔡᷹ⷰὐಽᏓ㧔8 ᦬ 12 ᣣᤨὐ㧕㧚✛ਣߪࠝࡈ࡜ࠗࡦ㧘✛྾ⷺߪࠝࡦ࡜ࠗࡦ㧘⿒ਣߪ 7 ᦬ਛᣨ߹ߢࠝࡈ࡜ࠗࡦ᷹ⷰࠍⴕߥߞߚ⥃ᤨ᷹ⷰὐࠍ␜ߔ㧚⿒྾ⷺߪ᳇⽎ᐡ࡮㒐ἴ⑼ቇᛛⴚ⎇ⓥᚲ࡮᧲ർᄢቇߩࠝࡦ࡜ࠗࡦ᷹ⷰὐࠍ⴫ߔ㧚ዊߐߥਣߣ
ᤊශߪߘࠇߙࠇ᧲ർᄢቇߩ⥄േಣℂߦࠃࠆ૛㔡ߣᧄ㔡ߩ㔡ᄩࠍ⴫ߔ㧚⿒✢ߪᵴᢿጀ㧘⿒ਃⷺߪ╙྾♿ߩἫጊࠍ␜ߔ㧚Fig.1 Distribution of temporary observation stations deployed by the Group for the aftershockobservations of the Iwate-Miyagi Nairiku Earthquake in 2008 (GIMNE2008). Greencircles and squares indicate offline and online temporary stations respectively. Redcircles denote the stations operated till the middle of July 2008. Red squares show theonline stations installed by JMA, NIED, and Tohoku University (TU). Star and small
circles indicate the epicenters of the main shock and aftershocks determined by theautomatic hypocenter determination system of TU.

㪄㩷㪈㪎㪋㩷㪄╙ 2 ࿑ วห૛㔡᷹ⷰࠣ࡞࡯ࡊߦࠃࠆ૛㔡ಽᏓ㧔6/23㨪30㧘7/2㨪7㧕㧚DD ࠻ࡕࠣ࡜ࡈࠖᴺߢ᳿ቯߐࠇߚ⚿ᨐࠍ␜ߔ㧚਄࿑㧦⷏ർ⷏-᧲ධ᧲ᣇะߩ㋦⋥ᢿ㕙࿑㧚ฝਅ࿑㧦ർർ᧲-ධධ⷏ᣇะߩ㋦⋥ᢿ㕙࿑㧚Ꮐਅ࿑㧦㔡ᄩಽᏓ㧚ᢿ㕙࿑ਛߩ⿒⦡྾ⷺߪ಴ᐫᢿጀߩ૏⟎㧘⿒✢྾ⷺߪ࿾⴫ᄌ⁁ߩ૏⟎ࠍ␜ߔ㧚08:01 ߣ 08:11 ߦ⊒↢ߒߚ೨㔡ࠍ㤥ᤊߢ㧘08:43 ߦ⊒↢ߒߚᧄ㔡ࠍ⊕ᤊߢ␜ߔ㧚ਣශ߇૛㔡ߢᷓߐࠍ⦡ߢ⴫ߔ㧚⿒ਃⷺߪ╙྾♿Ἣጊࠍ␜ߔ㧚Fig.2 Aftershock distribution determined using data obtained by GIMNE2008 and DD
tomography method. Top: WNW-ESE vertical cross sections. Bottom right: NNE-SSW
vertical cross section. Bottom left: Epicenter distribution. Color shows the depth.

㪄㩷㪈㪎㪌㩷㪄

╙ 3 ࿑ DD ࠻ࡕࠣ࡜ࡈࠖᴺߦࠃߞߡផቯߐࠇߚ㔡Ḯၞߦ߅ߌࠆ
S ᵄㅦᐲ஍ᏅಽᏓ㧚਄࿑㧦⷏ർ⷏㻌᧲ධ᧲ᣇะߩᢿ㕙㧚⊕ᤊߣ⊕ਣߪౣ᳿ቯߐࠇߚᧄ㔡ߣ૛㔡ߩ㔡Ḯࠍ␜ߔ㧚ਅ࿑㧦ᷓߐ 3kmߦ߅ߌࠆಽᏓ㧚⿒྾ⷺߪ࿾㔡ᢿጀߩ᭎⇛ࠍ⴫ߔ㧚
Fig.3 S-wave velocity perturbation distribution estimated from the temporary observation
data using DD tomography method. Top: WNW-ESE vertical cross-sections. Bottom:
Map view at the depth of 3 km.

㪄㩷㪈㪎㪍㩷㪄

㻌(a)(b)╙ 4 ࿑ ᄢ࿾㔡ߣᓸዊ࿾㔡ߣᷓㇱ᭴ㅧߩ㑐ଥ㧚ጟ↰࡮ઁ (2008)1) ߦࠃࠆ 1997 ᐕ߆ࠄ 2008 ᐕ߹ߢߩᓸዊ࿾㔡㔡ᄩಽᏓࠍ␜ߔ㧚Ꮐ࿑㧦ᓸዊ࿾㔡ߩᵴേߣㆊ෰ߩ M5.5 એ਄ߩ࿾㔡ߩᢿጀ㕙ߩ᭎⇛㧔⿒⍱ᒻ㧕㧚ฝ࿑㧦ᷓߐ 24km ߦ߅ߌࠆ S ᵄㅦᐲ஍Ꮕ㧚⊕චሼߪᷓㇱૐ๟ᵄᓸዊ࿾㔡ࠍ⴫ߔ㧚
Fig.4 Relation between the large earthquakes, microearthquakes and deep structure.
Hypocenters relocated by Okada et al. (2008)1) for the period from 1997 to 2008 areshown.
Left: Relation between the microearthquakes (black crosses) and large
earthquakes of M5.5 or larger (red rectangles).
Right: Relation between themicorearthquake and S-wave velocity perturbation distributions at the depth of 24 km.
White crosses indicate the deep crustal low-frequency earthquakes.

㪄㩷㪈㪎㪎㩷㪄

╙ 5 ࿑ ᧄ㔡ߩ೨ߩ࿾㔡ᵴേ㧚ጟ↰࡮ઁ (2008)1) ߦࠃࠆ 1997 ᐕ߆ࠄ 2008 ᐕ߹ߢߩ㔡ḮಽᏓࠍ␜ߔ㧚⊕ᤊ㧘㤥ᤊ㧘Ἧ⦡ਣߪ੹࿁ߩ࿾㔡ߩᧄ㔡㧘೨㔡㧔08:01㧘08:11㧕ߣ૛㔡㧔6/14-17㧘20㧕ࠍ␜ߔ㧚⿒⦡ߩࠪࡦࡏ࡞ߪ 1999 ᐕ㧘㕍⦡ߩࠪࡦࡏ࡞ߪ 2000 ᐕߩᵴേࠍ⴫ߔ㧚⿒⪉ဳ߇ 1999
ᐕ 4 ᦬ 19 ᣣߩ M4.3 ߩ࿾㔡㧘㕍⪉ᒻ߇ 2000 ᐕ 2 ᦬ 11 ᣣߩ M4.9 ߩ࿾㔡ࠍ␜ߔ㧚਄࿑㧦㔡ᄩಽᏓߣࡔࠞ࠾࠭ࡓ⸃㧚ਅ࿑㧦਄࿑ߩ A-A' ᢿ㕙㧚
Fig.5 Seismic activity before the main shock. Hypocenter distribution determined by Okada
et al. (2008)1) for the period from 1997 to 2008 is shown. White star, black stars andgrey circles respectively denote the main shock, foreshocks and aftershocks of the 2008
activity. Red and blue symbols show the 1999 and 2000 activities, respectively. Red andblue diamonds respectively denote the M4.3 event on April 19, 1999, and the M4.9 event
on February 11, 2000. Top: Epicenter distribution and focal mechanisms. Bottom:
Vertical cross section for the A-A' region in the top figure.

㪄㩷㪈㪎㪏㩷㪄╙ 6 ࿑ GPS ㅪ⛯᷹ⷰὐ ICNS ߇ᝒ߃ߚ࿾㔡ᤨ࿾Ზᄌേ㧚Ꮐ࿑㧦᷹ⷰὐ㈩⟎㧚ฝ࿑㧦ᄌ૏ߩᤨ㑆ᄌൻ㧚
Fig.6 Coseismic step detected by the GPS station ICNS. Left: GPS station distribution.
Right: Time histories of the displacement at ICNS.

╙ 7 ࿑ GPS ᷹ⷰ࠺࡯࠲߆ࠄផቯߐࠇߚ⍱ᒻᢿጀࡕ࠺࡞㧚(a) ࿾㔡ᤨ᳓ᐔᄌേߩ᷹ⷰ୯㧔㤥⍫ශ㧕ߣ⸘▚୯㧔⊕⍫ශ㧕ߩᲧセ㧚(b) ࿾㔡ᤨ਄ਅᄌേߩ᷹ⷰ୯㧔㤥᫔㧕ߣ⸘▚୯㧔⊕᫔㧕ߩᲧセ㧚
Fig.7 Rectangular fault model estimated from GPS data. (a) Comparison between theobservation (black arrow) and synthetic (white arrow) of the horizontal displacement.
(b) Comparison between the observation (black bar) and synthetic (white bar) of thevertical displacement.

㪄㩷㪈㪎㪐㩷㪄

(a) 140° 30’ 140° 45’ 141° 00’
06/14-06/2106/21-06/28

06/28-07/05

07/05-07/12

39° 15’39° 00’38° 45’0052008/06/14 - 2008/07/12
10(b) 39° 15’

200 400 600 800 1000
cm/yr

20HMYO0193
150OHYU39° 00’

10HNDR

30
25

520
15
100Cal.
Obs.
15cm/yr

38° 45’140° 30’

50Slip(cm )140° 45’141° 00’

141° 15’╙ 8 ࿑ Yagi and Kikuchi (2003)2) ߩᤨ㑆ଐሽࠗࡦࡃ࡯࡚ࠫࡦߦࠃࠅផቯߐࠇߚ૛ലߔߴࠅߩᤨⓨ㑆ಽᏓ㧚(a) 7 ᣣᲤߩࠬ࠽࠶ࡊ࡚ࠪ࠶࠻㧚(b) ផቯߐࠇߚ᳓ᐔᄌ૏ㅦᐲ㧔⊕⍫ශ㧕ߣ᷹ⷰ୯㧔㤥⍫ශ㧕ߩᲧセ㧚
Fig.8 Afterslip distribution using the time-dependent inversion method developed by Yagi andKikuchi (2003)2). (a) Snap shots every 7 days. (b) Comparison between the calculated
displacement velocity (white arrow) and observation (black arrow).

㪄㩷㪈㪏㪇㩷㪄

(a)

(a)

(b)

(b)

(c)

(c)

╙ 9 ࿑ 㔡ḮၞߩᲧᛶ᛫᭴ㅧ㧚ᝌ౉࿑ߩ㧟ᧄߩ᷹✢ߦߟ޿ߡ㧞ᰴరᢿ㕙ࠍ␜ߔ㧚ਣශߪวห᷹ⷰߦࠃ
ߞߡᓧࠄࠇߚ૛㔡ಽᏓࠍ␜ߔ㧚Fig.9 Resistivity distributions around the source area. Circles denote the hypocenters of theaftershocks relocated from the data obtained by GIMNE2008.

㪄㩷㪈㪏㪈㩷㪄