３．１．５

ひずみ集中帯発生にかかわる地殻構造の研究

目

次

(1)

業務の内容・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・132(a) 業務題目(b) 担当者(c) 業務の目的(d) ６か年の年次実施計画（過去年度は、実施業務の要約）
1) 平成１９年度
2) 平成２０年度
3) 平成２１年度
4) 平成２２年度
5) 平成２３年度
6) 平成２４年度(e) 平成２３年度業務目的(2)

平成２３年度の成果・・・・・・・・・・・・・・・・・・・・・・・・・・・133(a) 業務の要約(b) 業務の実施方法及び成果
1) 石狩平野東縁断層周辺での自然地震観測a) 自然地震観測の継続
b) 自然地震観測データ処理c) 予察的地震波速度構造解析
2) 石狩低地東縁断層帯周辺での比抵抗構造探査
3) 樽前山での地震観測a) 自然地震観測の継続
b) 自然地震観測データ処理・解析
4) 樽前山での地殻変動観測a) 繰り返し GPS 観測（第 3 回）b) 干渉 SAR 解析
5) 十勝岳地殻変動
a) 繰り返し GPS 観測（第 2 回）(c) 結論ならびに今後の課題(d) 引用文献(e) 成果の論文発表・口頭発表等(f) 特許出願，ソフトウエア開発，仕様・標準等の策定(3)

平成２４年度業務計画案・・・・・・・・・・・・・・・・・・・・・・・・・147

131

(1)

業務の内容

(a) 業務題目

ひずみ集中帯発生にかかわる地殻構造の研究(b) 担当者
所属機関

役職

国立大学法人北海道大学 教授

氏名

メールアドレス茂木 透

mogitisv@mail.sci.hokudai.ac.jp
mccopy_mm@mail.sci.hokudai.ac.jp

大学院理学研究院
同

教授

村上 亮

同

教授

谷岡 勇市郎 tanioka@mail.sci.hokaudai.ac.jp

同

准教授

高橋 浩晃

hiroaki@mail.sci.hokudai.ac.jp

同

准教授

橋本 武志

hashimoto@mail.sci.hokudai.ac.jp

同

准教授

大島 弘光

oshima@mail.sci.hokudai.ac.jp

同

准教授

勝俣 啓

kat@mail.sci.hokudai.ac.jp

同

助教

森 済

hmori@mail.sci.hokudai.ac.jp

同

助教

青山 裕

aoyama@mail.sci.hokudai.ac.jp

同

技術専門職員 高田 真秀

takada@sci.hokudai.ac.jp

同

技術専門職員 一柳 昌義

ichimasa@mail.sci.hokudai.ac.jp

同

技術専門職員 鈴木 敦生

suzuki@sci.hokudai.ac.jp

同

技術職員

山口 照寛

teruhiro@sci.hokudai.ac.jp

同

嘱託職員

前川 徳光

maekawa@sci.hokudai.ac.jp

同

博士研究員

奥山 哲

okuoku@mail.sci.hokudai.ac.jp

同

補佐員

黒井 和典

kuroi@mail.sci.hokudai.ac.jp(c) 業務の目的
北海道中軸帯では、島弧衝突の影響による活断層が分布しており、そこでのひずみ集中が測地学的に観測されている。その発生機構を解明するため、稠密な自然地震観測や地磁気地電流法(MT)観測を実施し、詳細な震源分布や地殻比抵抗構造等を把握する。それらの関連性を調査することで、地震活動と地殻不均質構造との関係や、地震発生に大きくかかわっていると考えられている地殻流体分布等を調査し、ひずみ集中発生機構に関する新たな知見を得る。
(d) ６か年の年次実施計画（過去年度は、実施業務の要約）
1) 平成１９年度：ひずみ集中帯にかかわる地殻構造の研究に用いる自然地震観測システムを整備するとともに、ターゲットである石狩東縁断層帯においてそれを効果的に運用できるような地震観測点配置等の検討を行った。
2) 平成２０年度：ひずみの集中が観測されている石狩低地帯周辺での地震活動を把握する
ため自然地震観測を開始した。地殻構造や流体の分布を調査するために、同地域で比抵抗構造探査を実施した。
3) 平成２１年度：自然地震観測を継続し、地震観測データの蓄積を図った。火山での地震、
132

地殻変動観測を継続した。比抵抗構造の空間密度を上げるための電磁気観測を実施した。
4) 平成２２年度：自然地震観測を継続し、地震観測データの蓄積を図った。火山での地震、地殻変動観測を継続した。比抵抗構造への海の影響について調べた。
5) 平成２３年度：自然地震観測を継続し、地震観測データの蓄積を図るとともに、予備解析を実施し、最適な震源決定手法について検討する。火山での地震、地殻変動観測を継続する。比抵抗構造について三次元性も考慮した予備的な解析を実施し、初期構造モデルを得る。
6) 平成２４年度：自然地震観測を継続し、かつ、比抵抗構造の空間密度を上げる電磁気観測を実施するとともに、地震データ・電磁気データの解析を実施し、詳細な震源分布や比抵抗構造等を得る。火山での地震、地殻変動観測データを総合的に解析する。これらのデータを総合的に検討し、ひずみ集中機構に関する知見を得る。
(e) 平成２３年度業務目的
ひずみの集中が観測されている石狩低地帯周辺での地震活動を把握するため、
自 然 地 震 観 測 を 継 続 す る と と も に 、 予備解析を実施し、最適な震源決定手法について検討する。ま た 、ひ ず み 集 中 と 火 山 活 動 と の 関 連 性 を 明 ら か に す る こ と を 目 指 し 、
樽 前 山 と そ の 周 辺 地 域 に お い て 、地 殻 変 動 源 の 深 さ 分 解 能 向 上 に 影 響 を 与 え る 観
測 網 の 空 間 配 置 を 最 適 化 す る た め 、 GPS、 合 成 開 口 レ ー ダ ー 、 そ の 他 の 地 殻 変 動
観 測 等 を 総 合 的 に 組 み 合 わ せ た 地 殻 変 動 観 測 手 法 を 開 発 し 、火 山 を 含 む 他 の ひ ず
み 集 中 帯 へ の 適 用 を 検 討 す る 。地 殻 構 造 や 流 体 の 分 布 を 調 査 す る た め に 、石 狩 地
域 、 比抵抗構造の空間密度を上げるための電 磁 気 観 測 （ マ グ ネ ト テ ル リ ク 法 探 査 等 ）を 実 施 す る と と も に 、 予備的な解析を実施し、初期構造モデルを得る。
(2)

平成２３年度の成果(a) 業務の要約
自 然 地 震 観 測 を 継 続 し 、ひ ず み の 集 中 が 観 測 さ れ て い る 石 狩 低 地 帯 周 辺 で の 地
震 活 動 を 把 握 す る た め 、 20 0 8 年 1 2 月 か ら 2 0 1 2 年 1 月 ま で の 間 に 1 2 5 2 個 の 震 源を 決 め と こ ろ 、石 狩 低 地 帯 で は 、北 部 地 域 及 び 中 南 部 で や や 地 震 活 動 が 活 発 で ある こ と が 分 か っ て き た 。ま た 、速 度 構 造 で は 石狩低地帯は大局的には低速度で特徴付
けられるとともに、その西部では高速度、東部ではより低速度の層が確認された。これらの結果を用いて最適な箴言決定手法についても研究を継続している。石 狩 地 域 の M T 法 に
よ る 比 抵 抗 構 造 探 査 で は 、周 囲 の 海 の 影 響 を 評 価 し 、そ の 影 響 を 受 け た デ ー タ を
除 い て 二 次 元 比 抵 抗 構 造 を 再 解 析 し 、三 次 元 モ デ ル の 初 期 構 造 も 作 成 し た 。石 狩
平野東縁断層帯の北部地域では断層帯の下部に周囲より比抵抗低い場所があることが推定された。
ひ ず み 集 中 と 火 山 活 動 と の 関 連 性 に つ い て 明 ら か に す る こ と を 目 指 し 、樽 前 山
及 び 十 勝 岳 地 域 に お い て 、地 殻 変 動 を 調 べ る た め 繰 り 返 し G P S 繰 り 返 し 観 測 及 び
干 渉 SAR 解 析 を 行 っ た 。GPS の 繰 り 返 し 観 測 は 順 調 に 進 め ら れ た が 、干渉 SAR 解析
は昨年まで解析に利用していた ALOS/PALSAR が 2011 年 4 月に観測を終了したため、計画を133

変更し、PS-InSAR による時系列解析を試みた。ノイズを除去しよいデータを得るためにはパラメータの調整が必要であることが分かってきた。これらのデータを解析し、総合的解析手法の開発を継続している。
(b) 業務の実施方法及び成果
1) 石狩平野東縁断層帯周辺での自然地震観測a) 自然地震観測の継続
本年度も引き続き平成 20 年度に石狩低地帯南部に新たに設置及び機器の高度化を行った 10 か所の自然地震観測点において地震観測を継続した。観測点には、高感度短周期微小地震計及び広帯域強震計が併設されており、これらのデータは 100 Hz の 24 bit でデジタ
ルサンプリングされた後、NTT 東日本フレッツ網を用いた IP 回線網を経由して北海道大学
へリアルタイムで伝送されハードディスクへ格納されている。今年度も特に大きなトラブ
ルもなく順調にデータが取得された。
b) 自然地震観測データ処理
本研究で設置された地震観測点から伝送されてくる地震波形データに加え、北海道大学の運用する既設の地震観測点、札幌市の運用する地震観測点、全国地震波形流通システム
を介してリアルタイムデータ交換を実施している気象庁、防災科学技術研究所の運用する地震観測点のデータを一元化したうえで、振幅の長時間平均と短時間平均の比からトリガーをかけてイベントの検出を行っている。トリガーされたイベント波形に対して、ノイズによるものを目視で除去し、地震によるデータについて P 相及び S 相の会話検測作業を随時実施している。
データの取得が始まった 2008 年 12 月から 2012 年 1 月までに読み取りが行われ震源が求まった地震の数は 1252 個であった。本研究では最終的にこの地域の地震波速度構造の推定を予定している。このため、上部マントルや下部地殻の地震波速度構造での波線も必要となるため、当該地域の沈み込んだスラブ内で発生するやや深発地震についても走時データ取得のため読み取り作業を行っている。震源決定の初期 P 波速度構造には北海道大学がルーチンの震源決定作業で用いている構造(笠原・他，2005 5) )を利用しているが、当該地域
では厚い第四紀堆積層が発達しているために、その影響を考慮するため本研究では表層部には 3.8 km/s の層を設定している。この構造を用いて予備震源を決定し、その結果得られる各観測点での走時残差を観測点補正値として与えて震源再計算を行って、これを初期震源とした。表層の堆積層に相当する遅い速度を持つ層を与えることで、すべての地震において地中に収束した震源を得ることができた。
図 1 には上記の方法を用いて震源決定されたすべての深さの震源を示している。やや深発地震は昨年までの結果と同様に対象領域の南半分にほぼ均等に発生している。本研究では特に石狩低地東縁断層帯南部をターゲットとしており、やや深発地震から得られる波線は当該地域の上部マントル・下部地殻域での走時データの蓄積に貢献するものと思われる。
図 2 には、石狩平野の浅発地震の活動を見るために、深さ 40 km 以浅の地震のみを示した。震源決定方法は図 1 と同じである。地震活動は、従前と同じく石狩低地東縁断層南部の断層ラインよりも東側で活発である。深さ分布を図 3 に示したが、この場所での震源の134

深さは 10 km よりも深く 30 km を超える場所でも発生が見られる。この活動は、石狩平野の東側で見られる東北日本弧の上部地殻が東傾斜に傾き下がる構造（Iwasaki et al., 2004
3)）と関連している可能性がある。

図1

自然地震観測網により得られた震央分布

図2

深さ 40 km 以浅の震源分布
135

図3

図 2 を東西方向に切った場合の地震の深さ分布c) 予察的地震波速度構造解析
昨年度に引き続き予察的な地震波速度構造の解析を実施した。用いたデータセットは昨年同様に、図 1 で示した初期震源のうち、深さ 40 km 以浅で北緯 42 度から 43 度 38 分、東経 141 度から 143 度 38 分の範囲に 2010 年 7 月までに発生した地震 166 個、P 波走時 2484個である。これに対し、Thurber(1983) 11) の手法を用いて P 波地震波速度構造と震源を同時に推定した。昨年度は 10 km の深さについて速度構造の検討を行ったが、今年度は深さ0 km 及び 20 km についても構造を推定した。また、チェッカーボードレゾリューションテ
ストを行った(図 4)。その結果、深さ 10 km において最も解の再現性が高く、深さ 0 km や
20 km では解の再現性が悪い結果となった。これらの結果は、今後蓄積されるデータや、
やや深発地震を利用することにより改善されると考えられる。
以上のような前提条件のもとに、図 5-1 にはそれぞれの深さにおける P 波速度構造のマップビューを、図 5-2 にはその東西断面を示した。これによると、石狩低地帯は大局的には低速度で特徴付けられるとともに、その西部では高速度、東部ではより低速度の層が確認される。これらの特徴は、屈折及び反射法地震探査より明らかにされた速度構造（Iwasaki
et al., 2004 3) ） と 調 和 的 で あ る 。 石 狩 低 地 帯 断 層 南 部 の 東 側 で 発 生 し て い る 地 震 は 、Iwasaki et al. (2004) 3) との比較では千島弧に押し込まれた東北日本弧上部地殻の上部付近で発生しているようにも見えるものの、今後震源の深さを含めた精査が必要である。
また、苫小牧市付近には周辺よりも高速度の異常が見られる。この付近での深層ボーリン
グでは、深さ 4500 m 付近に基盤をなす花崗岩類が認められ（Hanagata, 2002 1) ）、それと関連して当該海域は先第三系の変成岩・基盤岩からなる苫小牧リッジ構造となっていることから（国安，2010 6) ）、これらの基盤構造の高まりを反映していることも考えられる。
今後は、やや深発地震のデータを用いてより深い部分の速度構造の推定を実施し、これらの構造の深部での様態を明らかにする予定である。図の作成の一部に SEIS-PC（石川・中村，2005 8) ）を利用した。
136

図4

チェッカーボードリゾリューションテストの結果

137

図 5-1

深さ 0、10、20 km における P 波

図 5-2

地震波速度構造。

図 4-1 の東西断面における P 波速度構造。

2) 石狩低地東縁断層帯周辺での比抵抗構造探査
20 年度から 22 年度に行った広帯域 MT 法探査で、断層帯周辺の 40 の地点（図 6）でデータが得られた。これまでに、これらと既存の観測点のうち 4 本の測線上に属するデータを使用して、二次元インバージョンによる比抵抗断面の推定を行った。しかしながら、周辺の海と石狩低地帯に厚く分布する堆積層の低比抵抗によって生ずる異常な電流系による
MT データへの影響が懸念された。この場合、二次元解析では偽像を引き起こす可能性がある。そこで、本年度はこの影響を評価するために、海の分布と堆積層の低比抵抗を仮定した三次元モデル計算を行った。図 7（左）に示すように海水の比抵抗として 0.3 Ωm、堆積層には先の二次元構造を参考にして 3 及び 30 Ωm を与え、Fomenko and Mogi（2002）の三138

次元比抵抗構造に対する MT レスポンスを計算するコードにより、見かけ比抵抗、位相及びインダクションベクトルを計算した。その結果、インダクションベクトルについては 0.01
Hz 以下で観測値と同様の傾向を示し、構造の仮定が十分現実的であることが確認された（図 7 中・右）。このとき、見かけ比抵抗及び位相については、二次元構造の TE モードに相当する成分において、仮定した海や低比抵抗構造の影響が大きく現れた。海が存在する
場合と無い場合で比較すると、0.03 Hz 以下の周波数帯で、見かけ比抵抗は対数の 30%、位
相は 20º の最大変化が見られた（図 8）。したがって、堆積層の低比抵抗と周辺に三次元的に分布する海の低比抵抗が MT データに影響していることが明らかとなった。

図6

MT 観測点の分布。ISK が H20～22 年度に行った観測点を示す。茂木・他（2002） 7)による観測点（HDK）及び Yamaya (2008) 12) による観測点（TRM）も含む。Line A～D は二次元比抵抗断面の位置を示す。活断層線は中田・今泉編（2002）による。

図 7 （左）三次元モデル計算で仮定した海及び堆積層の深さ、比抵抗分布。
（中・右）0.01
Hz 及び 0.001 Hz のインダクションベクトル（実部）。赤が観測値、黒が計算値を示す。
139

図8

海の存在を仮定した場合と仮定しない場合のレスポンスの比較。上は見かけ比抵抗の対数の変化率（％）。下はインピーダンス位相の差（度）。

上記の結果は、TE モードの 0.03Hz 以下のデータを使用した二次元解析では、正しい構造が得られないことを示唆している。したがって、これに該当するデータを除去して、再び二次元インバージョンを行った。TM モードは 320-0.00034 Hz、TE モードは 320-0.03 Hzの見かけ比抵抗及び位相を使用し、Ogawa and Uchida (1996) 9) の二次元インバージョン
コードにより 4 つの比抵抗断面を得た（図 9）。約 7 km 以浅の構造は、平野部の厚い低比抵抗層（100 Ωm 以下）が特徴である。この層は断層帯の西側でより比抵抗が低く（10Ωm以下）顕著であるが、断層帯の東側では比較的高比抵抗である（10-100Ωm）。この比抵抗境界は断層帯表層から東に低角で傾斜する構造境界を成し、約 7 km の高比抵抗基盤まで延長する。この構造の特徴は、これまでに推定されている地震波速度構造（産業技術総合研究所，2007 10) ）とも矛盾せず、比抵抗構造の境界は観測地域東方の日高衝突帯より西進するスラスト構造（例えば，伊藤，2000 4) ）のデタッチメントと解釈できる。それより深部の構造は、数 100 Ωm 以上の高比抵抗で構成されているが、Line A の西端では低抵抗体が140

推定された。これは位置的に支笏カルデラの火山活動に関連した構造である可能性がある。
一方で、その他の場所は全体的になめらかで顕著な構造は推定されておらず、地殻内の流体を示すような低比抵抗体は見つからない。ただし、解析においてデータを除去したこと
により、深部構造の情報の一部も除去されていることになり、本来あるべき構造（深部の低比抵抗体など）が推定されていない可能性がある。この問題を解決するためには、海の構造を仮定した上で三次元インバージョンによる構造の推定を行う必要がある。

図9

二次元インバージョンにより得られた比抵抗断面。ITFZ：石狩低地東縁断層帯主部、ITFZ-S：石狩低地東縁断層帯南部、NPFZ：野幌丘陵断層帯。
1413) 樽前山での地震観測a) 自然地震観測の継続
平成 20 年度に樽前山南東域で機器の高度化を行った 2 か所の自然地震観測点（NSK 及びKRS）において観測を継続した。観測は広帯域地震計で行っており、順調にデータが所得されている。データは 400 MHz 帯の無線通信及び NTT フレッツ網を用いた IP 回線網を経由して北海道大学へリアルタイムで伝送されハードディスクに格納されている。
b) 自然地震観測データ処理・解析伝送された地震波形データを、北海道大学の運用する既存地震観測点のデータとあわせて検測作業を行った。2009 年 1 月から 2011 年 12 月までに樽前山周辺で発生した地震の震源分布を図 10 に示す。樽前山の山頂溶岩ドーム付近の地下 1 km 深を中心として、火山性の微小地震が発生している。樽前山北方の風不死岳直下（図 10 の KTR 北西側）でも時折群発活動が見られるほか、樽前山の南西域で深さ 6～8 km 程度の地震活動がほぼ定常的に発生していることが分かる。なお 2011 年 7 月に山頂北東のモーラップ付近（図 10MRP の南西側）を震源とする地震がまとまって発生した。

図 102009 年 1～12 月のデータを用いて決定された樽前山周辺の地震の震源分布。

4) 樽前山での地殻変動観測a) 繰り返し GPS 観測（第 3 回）
樽前山の地殻変動を調査するため、2011 年 10 月 12 日から 11 月 5 日にかけて第 3 回の繰り返し GPS 観測を行った。観測点は昨年と同じく北海道大学の管理する 14 点、国土地理142

院設置の三角点・水準点に各 1 点である。図 11 に 2 回の観測から得られた水平変動図を示す。昨年同様観測前に周囲の樹木の枝払いを実施したため、2009～2010 年のような異常は
TRE 以外では認められない。点 TRE は枝払いで視界を確保できる限度を超えて周囲の樹木が成長しており、移設を考慮すべき時期に来ていると考えられる。

図 11

繰り返し GPS 観測から得られた 2010 年 9 月～2011 年 10 月の樽前山周辺の地殻変動。

変位量は NSK（錦岡）を基準とした基線解析によるもの。背景の描画には SRTM3 DEM を用いた。
b) 干渉 SAR 解析
昨年まで解析に利用していた ALOS/PALSAR が 2011 年 4 月に観測を終了したため、計画を変 更 し 、 PS-InSAR に よ る 時 系 列 解 析 を 試 み る こ と に し た 。 PS （ Permanent/Persistent
Scatterer）-InSAR とは、植生等の観測時期によって変化するターゲットからの後方散乱を排除し、常に後方散乱パターンが一定であるターゲットのみを用いて干渉処理を行う技術である。解析には StaMPS（Hooper et al., 2007 2) ）を利用した。
2007 年 9 月 15 日の画像をリファレンスとした時の解析領域の干渉画像を図 12 に、基線情報の誤差，大気遅延等による誤差成分を分離した PS の位相を図 13 に示す。図 13 では
2008 年以降の画像にほぼ同様の位相変化が確認でき、明らかにリファレンス画像に含まれる誤差成分が分離できていない。これは PS の密度が全ピクセル数の 0.2 ％と非常に低く、PS 間の距離が大きすぎることが原因だと考えられる。樽前山の山頂部には殆ど植生は存在しないため、この PS 密度の低さの原因としては積雪が最も疑わしい。パラメータをこのような条件に最適化することにより結果を改善できるかもしれない。
143

図 12

図 13

2007 年 9 月 15 日を基準とした干渉画像。

StaMPS による解析の結果。画像の並びは図 2 に同じ。
144

5）十勝岳地殻変動
a) 繰り返し GPS 観測（第 2 回）
十勝岳の地殻変動を調査するため、2011 年 9 月 20 日から 10 月 7 日にかけて第 2 回の繰り返し GPS 観測を行った。観測点は昨年と同じく北海道大学の管理する 10 点、国土地理院
設置の三角点 2 点の計 12 点である。図 14 に 2 回の観測から得られた水平変動図を示す。2cm 程度の変位を示す点が多いが、山頂直下にあると考えられる圧力源の活動によるものとは考えにくい。

図 14

繰り返し GPS 観測から得られた 2010 年 7 月～2011 年 9 月の十勝岳周辺の地殻変動。

変位量は TRK（大麓山）を基準とした基線解析によるもの。背景の描画には SRTM3 DEM を用いた。

(c) 結論並びに今後の課題
ひずみ集中が見られる石狩東縁断層帯で継続されている稠密な地震観測を継続し、断層帯中南部と北部とでは地震活動が活発であることが分かってきた。更に観測を継続することによりこの傾向を確かめていく必要があろう。また、観測された地震を基にした P 波速度構造も明らかになってきたので、地震の観測を継続することにより更に精度の高い構造が得られるであろう。石狩低地東縁断層帯付近の比抵抗構造も明らかになってきており、周囲の海の影響を除いた比抵抗構造により、中部と南部での断層の活動性と構造との関連が見られる可能性もある。
樽前山の地殻変動観測は昨年度より観測が開始されたが、顕著な変動は観測されておらず、GPS の繰り返し観測及び干渉 SAR 観測による総合的な解析により、地殻変動源推定のための最適な観測点配置の検討を進めていくことになろう。しかし、ALOS/PALSAR が 2011年 4 月に観測を終了したため、計画を変更し、PS-InSAR による時系列解析を試みることになった。十勝岳 62－Ⅱ火口付近で顕著な地殻変動が観測され、火口付近が膨張していることが認められたが、その傾向は最近鈍ってきていることも分かってきた。引き続き地殻変動観測を継続する必要がある。
145

(d) 引用文献
1) Hanagata, S., Eocene shallow marine foraminifera from subsurface sections in he
Yufutsu-Umaoi district, Hokkaido, Japan, Paleontol. Res., 6, 147-178, 2002.
2) Hooper, A, P. Segall, and H. Zebker, Persisitent scatterer InSAR for crustal
deformation analysis, with application to Volcan Alcedo, Galapagos, J. Geophys. Res.,
112(B07407), doi:10.1029/2006JB004763, 2007.
3) Iwasaki, T, K. Adachia, T. Moriyab, H. Miyamachi, T. Matsushima, K. Miyashita,
T. Takeda, T. Taira, T. Yamadaa, K. Ohtake, Upper and middle crustal deformation ofanarc–arc

collisionacross,Hokkaido,

Japan,

inferredfromseismic

refraction/wide-angle reflection experiments, Tectonophys., 388 , 59– 73, 2004.
4) 伊藤谷生，日高衝突帯−前縁褶曲・衝上断層帯の地殻構造，石油技術協会誌，65，103-109，
2000．
5) 笠原

稔・小平秀一・本谷義信：1993年釧路沖地震の余震活動とそれ以前の北海道周辺の地震活動, 1993年釧路沖地震による被害の調査研究, 平成5年度文部省科学研究費突発災害調査研究成果報告, 27-34, 1993.
6) 国安稔，日高勇払地域，日本地方地質誌，北海道地方，631pp，2010．
7) 茂木透・日高 2000MT 探査グループ，日高地域での広帯域 MT 探査，月刊地球，24，485-487,
2002.
8) 中村浩二，石川有三：卓上の地震活動解析 SEIS-PC for Windows のこの8年,日本地震学会 2005 年秋季大会講演予稿集, 153, 2005.
9) Ogawa, Y. and T. Uchida, A two-dimensional magnetotelluric inversion assuming
Gaussian static shift, Geophys. J. Int ., 126, 69-76, 1996.
10) 産業技術総合研究所，石狩低地東縁断層帯の活動性及び活動履歴調査，「基盤的調査観
測対象断層帯の追加・補完調査」成果報告書 No.H18-8, 35pp., 2007.
11) Thurber, C. H., Earthquake locations and three-dimensional crustal structure in
the Coyote Lake area, central CaliforniaDetail Only Available, Jour. Geophys. Res.,88, 8226-36, 1983.
12) Yamaya, Y., Three dimensional resistivity structure of Tarumai Volcano by themagnetotelluric method including the effect of regional structure, Ph. D. thesis,
Hokkaido University, 110pp., 2008.
(e) 成果の論文発表・口頭発表等
著者

題名

発表先

発表年月日

山谷祐介，茂木透，本多

石狩低地東縁断層帯周辺の千葉市（日本地

平成 23 年 5亮，長谷英彰，鈴木敦生， 比抵抗構造(2)（ポスター）

球惑星科学連合

月 27 日橋本武志

2011 年大会）

Yamaya, Y., T. Mogi, H.

Resistivity imaging by

IUGG General

平成 23 年 7Hase, R. Honda, T.

Magnetotelluric method in

Assembly 2011月2日

Hashimoto and A. Suzuki

Ishikari- teichi-toen

146

fault zone (Poster)
山谷祐介，茂木透，長谷

石狩低地東縁断層帯周辺の英彰，本多亮，鈴木敦生， 比抵抗構造（ポスター）

2011 年 CA 研究

平 成 23 年

会

12 月 13 日

橋本武志(f) 特許出願，ソフトウエア開発，仕様・標準等の策定
1) 特許出願
なし
2) ソフトウエア開発
なし
3) 仕様・標準等の策定なし(3)

平成２４年度業務計画案

ひずみの集中が観測されている石狩低地帯周辺での地震活動を把握するため、
自 然 地 震 観 測 を 継 続 す る と と も に 、得 ら れ た 読 み 取 り 値 を 用 い て 地 震 波 速 度 構 造
解析を実施しそれを用いた高精度震源決定を行うことでひずみ集中帯での地震
活 動 特 性 を 明 ら か に す る 。ま た 、ひ ず み 集 中 と 火 山 活 動 と の 関 連 性 を 明 ら か に す
る こ と を 目 指 し 、ひ ず み 集 中 帯 の 西 縁 部 に 位 置 す る 樽 前 火 山 に お い て 、活 動 を 把
握 す る た め 地 震 観 測 を 継 続 す る と と も に 、 GPS の 繰 り 返 し 測 量 を 実 施 し 、 こ れ ま
で の 地 殻 変 動 観 測 デ ー タ を 総 合 的 に 解 析 す る 。ひ ず み 集 中 帯 の 東 端 に 位 置 す る 十
勝 岳 に お い て 、地 殻 変 動 源 の 深 さ 分 解 能 向 上 に 影 響 を 与 え る 観 測 網 の 空 間 配 置 を
最 適 化 す る た め 、 G PS、 合 成 開 口 レ ー ダ ー 、 そ の 他 の 地 殻 変 動 観 測 等 を 総 合 的 に組 み 合 わ せ た 地 殻 変 動 観 測 手 法 を 開 発 し 、火 山 を 含 む 他 の ひ ず み 集 中 帯 へ の 適 用
法 を 示 す 。地 殻 構 造 や 流 体 の 分 布 を 調 査 す る た め に 石 狩 地 域 で 行 っ た 電 磁 気 観 測
デ ー タ を 解 析 し 、三 次 元 比 抵 抗 構 造 を 求 め 上 で 、ひ ず み 集 中 の 地 下 構 造 の 特 徴 を明ら かに する 。国 内の 学会 等に おい て、今ま でに 得た 研究 成果 を発 表し 、新 た な専門的知見を得、成果の取りまとめに反映させる。

147