３－２６ 合成開口レーダ干渉法による平成 20 年（2008 年）岩手・宮城内陸地震に伴う地殻変動
Coseismic deformation of the 2008 Iwate-Miyagi inland earthquake derived
from SAR interferometry
防災科学技術研究所
National Research Institute for Earth Science and Disaster Prevention
2008年6月14 日に発生した平成20 年（2008年）岩手・宮城内陸地震に伴う地殻変動を調査するため，陸域観測技術衛星「だいち」に搭載されている合成開口レーダ（SAR（センサ名：PALSAR））による緊急観測データを用い，SAR 干渉解析を行った.その解析結果について報告する.
北行軌道に関する干渉画像（パス 402 において 2007 年 9 月 21 日と 2008 年 6 月 23 日に観測されたデータから作成）においては，震央の東側で 20cm を超えるスラントレンジ短縮，西側で 1m 程度のスラントレンジ伸張が検出された（第 1 図(a)）.南行軌道に関する干渉画像（パス 53 において 2006 年 8 月 4 日と2008 年 6 月 24 日に観測されたデータから作成）においては，震央の東側で 50cm 程度のスラントレンジ伸張，西側で 80cm 程度のスラントレンジ短縮が検出された（第 2 図(a)）.より震央に近い領域では，より大きな地殻変動が生じていたと想像されるが，両方の画像において干渉は得られなかった.これは，SAR 干渉法による検出限界を超えるような変形が生じていたためと推測される.震央の南東においては，東北大学・岩手大学合同調査グループや産業技術総合研究所の現地調査によって地震断層が発見されている.
その周辺について，それぞれの干渉画像を拡大してみると（第 1 図(b)および第 2 図(b)），地震断層が発見された地点に沿って，干渉縞の不連続が見られる.ただし，この不連続を横切る地殻変動量は，位相差における波数の不確定性を解決（アンラッピング）できないため，見積もることはできない.
SAR 干渉解析による地殻変動，GEONET（栗駒２観測点の変位は，国土地理院ホームページで公開されている図から判読），KiK-net によって得られた変位データを用いて断層すべり分布の推定を行った.
断層の走向を F-net によって求められている N209°E に設定し，断層面の浅部延長が現地調査によって発見されている地表断層と交わるように位置を固定した.その面中に設定した長さ 50km，幅 20km の面を 2km×2km に分割し，観測された地殻変動をよく説明するように，それぞれの小断層におけるすべりベクトルを推定した.このとき，面内におけるすべり分布は滑らかであると仮定し，その滑らかさの拘束の強さは， ABICを最小にするように決定した.また，断層の傾斜を20°,30°，40°に設定した場合についてそれぞれインバージョン解析を行ったところ，傾斜を 30°に設定した場合に，観測値－計算値の残差２乗和を最小とする解が得られた.以上の解析によって得られた断層すべり分布を第 3 図に示す.
震源の南にすべり量が大きい領域が推定された.この傾向は，地震観測データを用いた震源インバージョン結果とおおよそ調和的である.余震分布と比較すると，大きなすべりが推定された領域では余震が少なく，その周辺域で余震が多く発生しているように見えることが特徴である.また，すべり領域が地表近くに達している領域において，地表断層が発見さており，この結果は，現地調査の結果とも調和的であると考えられる.ただし，本結果においては，推定モーメントが地震観測結果よりも大きい，断層北部周辺において観測値と計算値の残差が大きいという問題が残されている.
謝辞：本研究で用いた PALSAR データは PIXEL において共有しているものであり,JAXA と東京大学地震研究所との共同研究契約により提供されたものである. PALSAR データの所有権は経済産業省およびJAXA にある.本解析においては GEONET のデータを使用した.
（小澤 拓）

㪄㩷㪉㪎㪇㩷㪄╙ 1 ࿑㩷 (a)ർⴕ゠㆏䈮䈍䈇䈩䋬2007 ᐕ 9 ᦬ 21 ᣣ䈫 2008 ᐕ 6 ᦬ 23 ᣣ䈮᷹ⷰ䈘䉏䈢䊂䊷䉺䈎䉌૞ᚑ䈚䈢ᐓᷤ↹௝.⿒ᤊශ䈲ᧄ㔡䈱㔡ᄩ䉕␜䈜.(b)╙ 1 ࿑(a)䈱⎕✢䈪␜䈜㗔ၞ䈱᜛ᄢ࿑.ਣශ䈲䋬᧲ർᄢቇ䊶ጤᚻᄢቇวห⺞ᩏ䉫䊦䊷䊒䈍䉋䈶↥ᬺᛛⴚ✚ว⎇ⓥᚲ䈱⃻࿾⺞ᩏ䈮䉋䈦䈩⊒⷗䈘䉏䈩䈇䉎࿾㔡ᢿጀ䈱૏⟎䉕␜䈜䋨䊖䊷䊛䊕䊷䉳䈮䈍䈇䈩౏㐿䈘䉏䈩䈇䉎࿑䈎䉌್⺒䋩.ㅘ᣿ᄥ✢䈲⋡ⷞ䈪್⺒䈚䈢ᐓᷤ❋䈱ਇㅪ⛯✢.⿒ਃⷺශ䈲ᩙ㚤ጊ䈱૏⟎䉕␜䈜.⢛᥊䈲࿖࿯࿾ℂ㒮⊒ⴕ
50m 䊜䉾䉲䊠 DEM 䈮䉋䉎࿾ᒻ㒶ᓇ࿑.
Fig. 1㩷 (a) Interferogram generated from ascending data (21 Sep. 2007 and 23 Jun. 2008). The red star
indicates the mainshock epicenter. (b) Enlarged interferomgram of Fig. 1 (a). Image area corresponds
to the broken rectangle in Fig. 1(a). Circles show locations of surface ruptures found by Tohoku
University, Iwate University, and AIST (read from figures published on homepages). The transparent
curve shows discontinuity of fringes traced by eyes. Red triangle indicates the location of Mt.
Kurikoma. Back ground image shows the topography depicted from GSI 50m-mesh DEM.

㪄㩷㪉㪎㪈㩷㪄╙ 2 ࿑㩷 (a)ධⴕ゠㆏䈮䈍䈇䈩䋬2006 ᐕ 8 ᦬ 4 ᣣ䈫 2008 ᐕ 6 ᦬ 24 ᣣ䈮᷹ⷰ䈘䉏䈢䊂䊷䉺䈎䉌૞ᚑ䈚䈢ᐓᷤ↹௝.⿒ᤊශ䈲ᧄ㔡䈱㔡ᄩ䉕␜䈜.(b)╙2 ࿑ (a)䈱⎕✢䈪␜䈜㗔ၞ䈱᜛ᄢ࿑.ਣශ䈲䋬᧲ർᄢቇ䊶ጤᚻᄢቇวห⺞ᩏ䉫䊦䊷䊒䈍䉋䈶↥ᬺᛛⴚ✚ว⎇ⓥᚲ䈱⃻࿾⺞ᩏ䈮䉋䈦䈩⊒⷗䈘䉏䈩䈇䉎࿾㔡ᢿጀ䈱૏⟎䉕␜䈜䋨䊖䊷䊛䊕䊷䉳䈮䈍䈇䈩౏㐿䈘䉏䈩䈇䉎࿑䈎䉌್⺒䋩.ㅘ᣿ᄥ✢䈲⋡ⷞ䈪್⺒䈚䈢ᐓᷤ❋䈱ਇㅪ⛯✢.⿒ਃⷺශ䈲ᩙ㚤ጊ䈱૏⟎䉕␜䈜.⢛᥊䈲࿖࿯࿾ℂ㒮⊒ⴕ
50m 䊜䉾䉲䊠 DEM 䈮䉋䉎࿾ᒻ㒶ᓇ࿑.
Fig. 2㩷 (a) Interferogram generated from descending data (4 Aug. 2006 and 24 Jun. 2008). The red starindicates the mainshock epicenter. (b) Enlarged interferomgram of Fig. 2 (a). Image area corresponds tothe broken rectangle in Fig. 2(a). Circles show locations of surface ruptures found by Tohoku
University, Iwate University, and AIST (read from figures published on homepages). The transparent
curve shows discontinuity of fringes traced by eyes. Red triangle indicates the location of Mt.
Kurikoma. Back ground image shows the topography depicted from GSI 50m-mesh DEM.

㪄㩷㪉㪎㪉㩷㪄╙ 3 ࿑㩷 (a)SAR ᐓᷤ⸃ᨆ⚿ᨐ䋬KiK-net䋬GEONET 䈮䉋䈦䈩ᓧ䉌䉏䈢࿾Ზᄌേ䈎䉌᳞䉄䈢ᢿጀ䈜䈼䉍ಽᏓ.⿒ᤊශ䈲ᧄ㔡䈱㔡Ḯ䋨᳓ᐔ૏⟎䋩䉕␜䈜.(b)࿾࿑਄䈮ᛩᓇ䈚䈢䈜䈼䉍ಽᏓ.㤥ὐ䈲᳇⽎ᐡ৻రൻ㔡Ḯ䈮䉋䉎૛㔡ಽᏓ.✛ਣශ䈲䋬᧲ർᄢቇ䊶ጤᚻᄢቇวห⺞ᩏ䉫䊦䊷䊒䈍䉋䈶↥ᬺᛛⴚ✚ว⎇ⓥᚲ䈱⃻࿾⺞ᩏ䈮䉋䈦䈩⊒⷗䈘䉏䈩䈇䉎࿾㔡ᢿጀ䈱૏⟎䉕␜䈜䋨䊖䊷䊛䊕䊷䉳䈮䈍䈇䈩౏㐿䈘䉏䈩䈇䉎࿑䈎䉌್⺒䋩.㕍⎕✢䈲ᢿጀ䊝䊂䊦ផቯ䈮䈍䈇䈩⸳ቯ䈚䈢ᢿጀ㕙䈪䈅䉎.⢛᥊䈲࿖࿯࿾ℂ㒮⊒ⴕ 50m 䊜䉾䉲䊠 DEM 䈮䉋䉎࿾ᒻ㒶ᓇ࿑.
Fig. 3㩷 (a) Fault-slip distribution estimated from crustal deformations by InSAR, KiK-net, and GEONET. Red
star shows mainshock hypocenter. (b) Fault-slip distribution projected onto the map. Dots indicate
aftershock epicenters determined by JMA. Green circles represent locations of surface ruptures foundby Tohoku University, Iwate University, and AIST (read from figures published on homepages). Blue
broken rectangle indicates fault plane used in the estimation of fault-slip distribution. Back groundimage shows the topography depicted from GSI 50m-mesh DEM.

㪄㩷㪉㪎㪊㩷㪄