３ .３

ＧＰＳ観測による詳細なひずみ分布の解明

目次

(1)

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

平成２２年度の成果・・・・・・・・・・・・・・・・・・・・・・・・・・・257(a) 業務の要約(b) 業務の実施方法及び成果
1) 観測の実施
2) データ解析
3) 変形モデルの検討(c) 結論ならびに今後の課題(d) 引用文献(e) 成果の論文発表・口頭発表等(f) 特許出願，ソフトウエア開発，仕様・標準等の策定(3)

平成２３年度業務計画案・・・・・・・・・・・・・・・・・・・・・・・・・268

254

(1)

業務の内容

(a) 業 務 題 目

ＧＰＳ観測による詳細なひずみ分布の解明(b) 担当者
所属機関

氏名

メールアドレス教授

鷺谷 威

sagiya@seis.nagoya-u.ac.jp

同教授

木股 文昭

kimata@seis.nagoya-u.ac.jp

同

特任研究員

松多 信尚

matta@seis.nagoya-u.ac.jp

同

大学院生

小澤 和浩

同

大学院生

山本 淳平

同

大学院生

岡崎 仁一

同

大学院生

佐々木 朋樹

同大学院生

白坂 祥之介

同

大学院生

竹村 裕樹

同

大学院生

若杉 貴浩

国立大学法人名古屋大

技術職員

奥田 隆

okuda@seis.nagoya-u.ac.jp

同

技術職員

堀川 信一郎horikawa@seis.nagoya-u.ac.jp

国立大学法人北海道大

准教授

高橋 浩晃

hiroaki@mail.sci.hokudai.ac.jp

同

技術職員

山口 照寬

teruhiro@sci.hokudai.ac.jp

同

技術職員

一柳 昌義

ichimasa@mail.sci.hokudai.ac.jp

教授

三浦 哲

miura@aob.geophys.tohoku.ac.jp

同助教

太田 雄策

ohta@aob.geophys.tohoku.ac.jp

同

技術職員

出町 知嗣

demachi@aob.gp.tohoku.ac.jp

教授

加藤 照之

teru@eri.u-tokyo.ac.jp

同

助教

福田 淳一

fukuda@eri.u-tokyo.ac.jp

国立大学法人富山大学

教授

竹内 章

takeuchi@sci.u-toyama.ac.jp

同准教授

楠本 成寿

kusumoto@sci.u-toyama.ac.jp

国立大学法人京都大学

准教授

宮崎 真一

shinichi.miyazaki@kugi.kyoto-u.ac.jp

教授

田部井 隆雄tabei@kochi-u.ac.jp

准教授

松島 健

mat@sevo.kyushu-u.ac.jp

国立大学法人名古屋大

役職

学大学院環境学研究科

学全学技術センター

学大学院理学研究院

国立大学法人東北大学
大学院理学研究科

国立大学法人東京大学
地震研究所

大学院理工学研究部

大学院理学研究科
国立大学法人高知大学
理学部
国立大学法人九州大学
大学院理学研究院
255

准教授

中尾 茂

nakao@sci.kagoshima-u.ac.jp

室長

棚田 俊収

tanada@bosai.go.jp

同

研究員

河野 裕希

yuhki-k@bosai.go.jp

神奈川県温泉地学研究研究員

原田 昌武

harada@onken.odawara.kanagawa.jp

国立大学法人鹿児島大学理学部
独立行政法人防災科学技術研究所所(c) 業務の目的
ひ ず み 集 中 帯 に お い て 精 密 な GPSキ ャ ン ペ ー ン 観 測 を 実 施 し 、 ひ ず み 集 中 帯 を横 切 る 方 向 で の 地 殻 変 動 分 布 の 現 状 を 明 ら か に す る と と も に 、得 ら れ た 地 殻 変 動分布を説明できる物理モデルを構築する。
(d) ５ヵ年の年次実施計画（過去年度は、実施業務の要約）
1) 平 成 ２ ０ 年 度 ： ひ ず み 集 中 帯 に お け る 詳 細 な ひ ず み 分 布 を 調 査 す る た め の 準
備 と し て 、必 要 な 観 測 機 器 を 整 備 し 、国 土 地 理 院 G P S観 測 網 G E ON E T の デ ー タ を
有 効 活 用 で き る よ う に 、 ひ ず み 集 中 帯 を 横 切 る 2箇 所 の 測 線 （ 新 潟 県 上 越 市 か ら南 魚 沼 市 、 柏 崎 市 か ら 南 魚 沼 市 を 想 定 ） に 沿 っ て 、 そ れ ぞ れ 、 2 - 3km 程 度 の 間 隔で GP S の キ ャ ン ペ ー ン 観 測 用 の 観 測 点 を 設 置 し 、 年 1 回 の キ ャ ン ペ ー ン 観 測 を 実施 し た 。観 測 は 、各 観 測 点 に 設 置 さ れ た ボ ル ト に ア ン テ ナ を 固 定 し て 1 〜 2 週 間 の連続観測を実施した。年周変化の影響を避けるため、観測は秋に実施する。得られ た 観 測 デ ー タ を 周 囲 の G E ON E T 観 測 点 と 一 緒 に 解 析 し 、毎 日 の 精 密 な 座 標 値 を得た。
2) 平成２１年度：GPS観測データ解析のための体制を整備した。また，前年度設置された観測点において年1回、2ヶ月程度のキャンペーン観測を実施した。得られたデータを周囲のGEONET観測点のデータと合わせた解析を行って精密な日座標値を算出し、前年度の座標値との比較から1年分の変位ベクトルを得た。
3) 平成２２年度：平成20年度に設置された観測点において年1回のキャンペーン観測を実施し、データ解析を行なうと共に、これまでのデータから初期的なモデルを構築する。
4) 平成２３年度：平成20年度に設置された設置された観測点において年1回のキャンペー
ン観測を実施し、データ解析を行なうと共に、自然地震観測や地殻構造探査などから得られる地下構造を参考にモデルの改良を行なう。
5) 平成２４年度：平成20年度に設置された設置された観測点において年1回のキャンペー
ン観測を実施し、データ解析を行なうと共に、各サブプロジェクトから得られる成果を参考に、ひずみ集中を引き起こしている地殻内の短縮変形の分布や断層幾何形状およびすべり 分布などを推 定し、ひずみ集中帯を横切る詳細な地殻変動分布とそれを説明する物理モデルを構築する。
(e) 平成２２年度業務目的
平成 20 年度に設置された観測点 49 ヶ所において、3〜4 週間のキャンペーン観測を 1回実施する。さらに、GPS 観測データ解析のための体制を整備してデータ解析を行い、前256

年度の解析結果との比較から地殻変動分布に関する最初の結果を得る。
(2)

平成２２年度の成果(a) 業務の要約
平成 21 年 7 月から 12 月にかけて、新潟県糸魚川市・上越市・十日町市・南魚沼市・長岡市・柏崎市の合計 49 観測点において GPS 観測を実施した。このうち、9 観測点については 7 月から 11 月の 4 ヶ月間、残り 40 観測点では 9 月から 11 月の 2 ヶ月間の観測を実施し、昨年度までよりも質・量ともに上回る良好なデータを取得した。なお、1 観測点（あさひ荘）については、既存の GPS アンテナ固定ボルトが、設置した建物の改修によって使用不能となったため、新たに設置したボルトで観測を実施した。こうして取得した観測データを周囲の GEONET 観測点 28 点および IGS 観測点 8 点と合わせてデータ解析を実施し、精密な日座標値を算出した。得られた座標値のばらつきは水平成分で約 3mm、上下成分で約 13mm と、連続観測点と遜色ない結果が得られた。こうして得られた座標値と過去 2 年間の観測で得られた座標値を用いて、ひずみ集中帯の地殻変動速度分布を推定した。その結果、糸魚川市周辺、東頸城丘陵周辺に東西方向の短縮変形の集中域が見られることが明らかとなった。この地域の地殻変形モデル構築の第一段階として、2004 年新潟県中越地震、2007 年新潟県中越沖地震を含む時期の GEONET 観測網による GPS 座標値データに基づいてひずみ集中帯の変形と応力蓄積に関する考察を進めた。その結果、これら2 つの地震を境に顕著な地殻変動パターンの変化が見られること、地震の発生等に伴って、地殻深部での非地震性の断層すべりが影響を受けている可能性のあることが分かった。本
研究の成果の一部を 2010 年 8 月に台北で開催された西太平洋地球物理学会議、2010 年12 月にサンフランシスコで開催された米国地球物理学連合大会、2011 年 1 月にバンダアチェ（インドネシア）で開催された国際ワークショップで発表するとともに、これらの学会において断層変形過程のモデル化に参考となる情報を収集した。
(b) 業務の実施方法及び成果
1) 観測の実施
平成 20 年度に設置した観測点において、平成 22 年 7 月から 11 月にかけて、GPS の連続観測を実施した。観測点の座標値をより正確に得るため、昨年よりも観測期間を 2 ヶ月に延長するとともに、9 観測点においては 7 月から 11 月までの最長約 4 ヶ月の観測を実施した。うち 7 観測点においては通常 30 秒サンプリングの観測を行うところを 1Hz サンプリングとし、大地震発生時等の地震波の検出が可能となるような観測を実施した。なお、
馬場小学校観測点では、観測機器の不調によりデータが取得できなかったため、12 月 10日まで追加の観測を実施した。長岡市小国支所観測点では、従来から設置していた機器の不調のために今年度は観測データが得られなかった。また、山古志観測点では、機器の不調により観測データが 3 日分しか得られなかったが、過去 2 年分のデータを合わせて地殻変動速度を算出するには十分と考え、追加の観測は実施しなかった。
あさひ荘観測点では、昨年度に屋上工事が実施された結果、従来の観測点が使用不能となっていた。そこで、今年度は新たにアンテナ固定用のポールを設置し、新観測点におい257

てデータを取得した。
以上述べた点以外では、不順な天候のためにソーラー電源が十分な電力を供給できなくなったり、受信機のメモリが一杯になったりして、観測が途中で停止した場所はあったが、1〜4 ヶ月間の連続観測データを取得することができた。今年度の観測の様子を図 1 に、また、観測データの取得状況を表 1 にまとめる。
観測に際しては、昨年度本研究予算で購入した GPS 受信機に加え、各大学で保有する既存の機器も合わせて使用した。その際、各観測点において、昨年度使用したのと同一機種の GPS アンテナを使用するよう注意した。これは、GPS アンテナの機種により、計算される座標値にバイアスを生じることが知られているため、そうした誤差を避ける目的である。

図1

観測機器設置状況および観測作業の様子。
（左上）保倉小学校。
（右上）磯部小学校。

（下）国府小学校。
258

表1

GPS キャンペーン観測によるデータの取得状況

259

2)

データ解析
本研究で得た観測データ（49 点）に、周囲の GPS 連続観測（28 点）および IGS 観測点（7 点）を加え、合計 84 点のデータ解析を実施して、毎日の精密な座標値を推定した。
データ解析には Bernese GPS Software(version 5.0)を使用し、IGS による精密暦および地球回転パラメータを使用して解析を実施した。また、解析にあたっては、ITRF2005 座標系に準拠し、IGS 観測点（TSKB を除く 6 点）を ITRF2005 で与えられる座標値に強く拘束して座標値を得た。
その結果、過去 3 年分の座標変化を直線近似した場合の残差のばらつき（RMS）として、東西成分では 1.0-6.7mm（平均 2.3mm）、南北成分で 1.1-4.3mm（平均 1.9mm）、上下成分で 4.6-10.0mm（平均 6.7mm）という結果が得られた。積雪等の影響で観測期間中に一部系統的に座標値のずれている日があり、こうした日のデータに注意してさらに解析を行えばさらなる精度の向上が可能である。糸魚川市の連続観測点 0245 に対する座標時系列の例を図 2 に示す。固定点に近い糸魚川市の浦本小学校（URAS）や、高田平野に位置する保倉小学校（HOKS）ではあまり顕著な変動が見られないが、十日町市の旧名ヶ山小学校（MKYS）や田沢小学校（TAZS）では糸魚川に対して西向きの動き、すなわち観測網全体で東西短縮となるような変形が生じていることが分かる。観測網周辺の連続観測点については 2008 年以降の毎日のデータを使って座標値を計算しており、精密な座標値変化が得られている。上越市大潟の 0241 観測点および南魚沼市六日町の 0246 観測点の座標値変化を図 2 に示した。
これらの日座標値時系列に対して図 2 に示すように最小二乗法により直線を当てはめ、その傾きを求めることで各観測点の地殻変動速度を求めた。その結果を図 3 に示す。昨年までの結果と異なり、観測地域内でほぼ整合的な地殻変動速度分布を得ることができた。
キャンペーン観測点と連続観測点でも系統的な速度ベクトルのずれなどは見られず、精度の高い地殻変動速度分布が得られつつあることが分かる。
図 3 からも明らかなように、この地域の地殻変動は、ほぼ東西方向の短縮変形が主である。観測地域内には西から順に糸魚川—静岡構造線、高田平野の西縁および東縁、十日町盆地、六日町盆地などにほぼ南北走向の活断層の存在が知られている。これらの活断層に関連する地殻変動がどのように生じているかが本研究の最大の関心事であるが、図 3 を見る限り、最も顕著な東西短縮変形が生じているのは、活断層のトレースの周辺ではなく、上越市と十日町市の境界付近に広がる東頸城丘陵の中央部付近のように見える。こうした傾向は、地殻変動速度の東西成分の変化の様子を経度に対してプロットした図 4 においても明らかであり、東経 138.5 度付近で東西方向の速度成分が 4mm/年程度不連続に変化しているように見える。また、糸魚川市と上越市の間でも 3-4mm/年程度の短縮が生じているように見える。こちらは、東頸城丘陵のような不連続な変化ではなく、連続的な変化のようである。より詳細な議論を行うためには観測データの精度をさらに高める必要があると考えられる。

260

図2

座標時系列の例（左上）浦本小学校、（右上）保倉小学校、（左下）旧名ヶ山小学校、（右下）田沢小学校。ともに糸魚川市の 0245 観測点に対する動きを示す。赤い線は最小二乗法により求めた近似直線である。

261

図 2（続き）

座標時系列の例（左）0241（上越市大潟）、（右）0246(南魚沼市六日町)

ともに糸魚川市の 0245 観測点に対する動きを示す。赤い線は最小二乗法により求めた近似直線である。

図3

水平地殻変動速度の分布。基準点は糸魚川（0245）である。

262

図 4東西方向の地殻変動速度分布プロファイル。影をつけた部分で速度が急変しているように見える。
3) 変形モデルの検討
既に述べたように、2010 年度までの観測により、新潟県上越—中越地域のひずみ集中帯における詳細な変形の様子が明らかになりつつある。本観測地域では、上越市と十日町市の境界付近の東頸城丘陵において最も顕著な変形集中が生じている可能性が高い。
この地域は、日本海拡大に伴って生じた厚い堆積層が短縮変形を受けて褶曲していると考えられており、地表付近でも活褶曲構造が見られる。観測網に近い十日町市室野では、
図 5 に示すような泥火山が存在するが、これも東西短縮に伴う活褶曲の背斜部で生じたものと考えられ、GPS で検出された変形の集中と調和的である（新谷・田中,2009）。このこ263

とから、東頸城丘陵における集中した東西短縮変形は主として非弾性的なものであり、活断層等における地震の準備過程を直接反映したものでは無いと考えられる。

図5

十日町市室野の泥火山（写真中央）

一方、糸魚川と高田平野の間で見られる短縮変形はより広い範囲に広がっている可能性があり、高田平野西縁断層で 1751 年にマグニチュード 7 クラスの地震が発生したと考えられていることも考え合わせれば、活断層の応力蓄積を反映した短縮変形である可能性がある。同様に、十日町断層帯や六日町断層帯の南部の地殻変動の詳細は現時点において明らかとなっていないが、この地域は 2004 年新潟県中越地震で破壊しておらず、マグニチ
ュード 7 クラスの地震が起きる可能性のある場所と考えられることから、さらに地殻変動観測の精度を向上させて検討を行う必要がある。
本研究の稠密 GPS 観測により得られる詳細な地殻変動分布に基づいてひずみ集中帯の変形過程をモデル化する第一段階として、2004 年新潟県中越地震から 2007 年新潟県中越沖地震の発生へと至る過程に関して、断層間の相互作用について検討し、シナリオ構築作業を実施した。
これらの地震は、わずか 3 年間の間に 30km ほどしか離れていない 2 ヶ所で M6.8 の地震が発生した事象であり、両者の間には因果関係があると考えられる。また、GPS 観測結果によれば、2004 年以前にはこれら 2 つの地震の震源域の間で短縮変形が進行する一方で、2004 年中越地震の震源域付近では、2004 年の地震発生前に短縮変形が次第に減速していた。こうした観測結果から、2 つの震源域の間では地震発生前に逆断層が非地震性のすべりを生じていたと考えられる。このすべりに伴う短縮変形は、本研究の GPS 観測によって東頸城丘陵で見出された変形と類似のものと考えることができる。こうした非弾性変形は 2004 年の地震後に鈍化しており、地震に伴う応力変化の寄与があったと考えられる。図 6 は 2004 年中越地震によって周囲の断層運動が受ける影響をクーロン破壊関数で評価したものである。2007 年中越沖地震は主として東に傾斜する逆断層が活動したものと264考えられているが、2004 年の地震によって東に高角で傾斜する逆断層は破壊し易くなることが示されており、両者の因果関係を支持する。一方、西傾斜の逆断層は水平に近い場合を除いて破壊しにくくなるような応力変化であり、2 つの震源域の間では短縮変形が抑制された可能性がある。

図 62004 年中越地震の断層運動によるクーロン破壊関数の変化。中越沖地震の震源域（左側の楕円）では、西傾斜断層では低角の場合、東傾斜断層では高角の場合に破壊が促進されることが分かる。
(c) 結論ならびに今後の課題
新潟県上越・中越地域の 49 観測点において 1〜4 ヶ月間の GPS キャンペーン観測を実施した。このデータを周辺の GEONET 観測点のデータと合わせて解析することにより、各観測点の日座標値を精密に推定し、地殻変動の分布を推定した。その結果、東頸城丘陵中央部で集中した東西短縮変形が見出された。また、観測網西側の糸魚川市から高田平野の間でも年間 3-4mm の短縮変形がみられ、断層の応力蓄積過程と関連している可能性が考えられる。今後、観測を継続することによって精度向上を図り、その結果に基づいてひずみ集中帯の変形様式ならびに地震発生メカニズムの解明につなげていきたいと考える。
今年度のまとめを進めていた平成 23 年 3 月 11 日に平成 23 年東北地方太平洋沖地震（M9.0）が発生した。この地震は、知られている中では日本最大の地震であり、国土地理院の GPS 観測結果によれば、本研究の観測地域においても 5-10cm 程度の東向き変位が地265

震時に生じたと考えられる。また、この大地震に続いて、翌 3 月 12 日未明には、新潟・長 野 県 境 付 近 で M6.7 の 地 震 が 発 生 し た 。 こ の 地 震 で も 震 源 域 に 近 い 松 之 山 観 測 点 で39.3cm、安塚観測点で 6.8cm と顕著な地殻変動が生じた。さらに、東北地方太平洋沖地震の発生後には東日本全体で大規模な余効変動が生じており、今後数年単位で継続する可能性がある。こうした状況から、これまで本研究で調査を進めてきたひずみ集中帯の地殻変動にも顕著な変化が生じるものと予想される。こうした状況は研究開始時に想定していなかった事態であるため、平成 23 年度以降の研究の進め方については、早い段階で再検討する必要がある。まず、平成 23 年度の早い段階で一部の GPS 観測を実施し、その結果を踏まえて今後の方針について検討する必要がある。
(d) 引用文献
1) 新 谷 俊 一 ・ 田 中 和 広 , 新 潟 県 十 日 町 市 に お け る 泥 火 山 噴 出 物 の 起 源 , 地 学 雑 誌 , 118,340-349, 2009.
(e) 成果の論文発表・口頭発表等
著者

題名

発表先

発表年月日

鷺谷威，伊藤武男，高日本海東縁ひずみ集中帯におけ

千葉市（日本地

平成 22 年 5

橋浩晃，三浦哲，太田るGPS稠密観測（序報）（ポスタ球惑星科学連月雄策，加藤照之，福田ー）合 2010 年 大淳一，竹内章，楠本成会）

寿，宮崎真一，田部井隆雄，松島健，中尾茂，棚田俊收，原田昌武，奥田隆，堀川信一郎，
河野裕希
鷺谷威，大園真子，平稠 密 GPS観 測 か ら 見 る 跡 津 川 断

千葉市（日本地

平成 22 年 5

原和朗，橋本学，細善

層の応力蓄積過程（口頭）

球惑星科学連月信，和田安男，竹内章，合 2010 年 大道家涼介会）
Crustal deformation around the

台北（Western平 成 22 年 6Nishimura, and K.

northernPacific

月

Ozawa

TectonicLine

centralJapan,T.

Sagiya,T.

Itoigawa-Shizuoka
faultzone,revealedbyGeophysics
Meeting）

dense GPS observation（口頭）鷺谷威T. Sagiya

日本列島内陸部の変形に関する

広島市（日本地平 成 22年 10

諸 問 題 —測 地 観 測 の 視 点 か ら —震 学 会 2010 年

月

（口頭）

秋季大会）
inland東京（日中内陸

平 成 22年 11

earthquakes: implications from

地震ワークシ月

Crustalstrain

and266dense GPS observations（口頭） ョップ）鷺谷威

1586年 天 正 地 震 の 震 源 断 層 ： 測

名古屋市（日本平 成 22年 11

地学的な視点（口頭）

活断層学会シ月ンポジウム）
T. Sagiya, K. Ozawa,

Consistency and inconsistency

サンフランシ平 成 22年 12M.

T.

between geodetic and geologic

ス コ （ AGU

月

Y.

fault slip rates in central Japan

Fall Meeting）

Ohzono,Nishimura,andHoso（口頭）

松多信尚, 太田陽子,

台湾東海岸の津波堆積物

安藤雅孝, 西川由香,原口強,Lin

千葉市（日本地

平成 22 年 5球惑星科学連

月合 2010 年 大Cheng-Horng会）

N. Matta, Y. Ota, M.

Possibility of paleotsunami at

台北（Western

平 成 22 年 6Ando, T. Haraguti,

the Chenggong site, east coast

Pacific

月

Y. Nishikawa, W.S.

of Taiwan（口頭）Geophysics
Meeting）

Chen
N.

Matta,S.

Relationshipofthesurface

台北（Western

平 成 22 年 6月

Ishiguro, M. Murase,

deformation between long term

Pacific

W.S. Chenand short term at Yuli area

Geophysics

along the Longitudinal ValleyMeeting）

Fault in Eastern Taiwan
S.

Creeping deformation along theバンダアチェ

平 成 23 年 1

Ishiguro, M. Murase,

Longitudinal Valley fault at(Internationa月

W. S. ChenYuli

l workshop onN.

Matta,area,

eastern

Taiwan,estimated by precise leveling

Geo-hazard

survey and aerial photographic

Investigation

surveyofSumatran

fault

andSoft-launching Geophysics
department）
M.

Ohzono,T.

Strainaccumulation

process

Sagiya, K. Hirahara,around the Atotsugawa faultM.

systemHashimoto,A.

Takeuchi, Y. Hoso, Y.

intheGeophys.

J.

平成23年

Int.

Niigata-Kobe

Tectonic Zone, central Japan

Wada, K. Onoue, F.
Ohya, and R. Doke
鈴木康弘，杉戸信彦， 平 均 変 位 速 度 分 布 に 基 づ く 糸 魚
隈元崇，澤祥，渡辺満

川 ―静 岡 構 造 線 断 層 帯 北 部 の 地
267

活断層研究,

平 成 22年 1033, 1-14

月

久，松多信尚，廣内大震発生予測助，谷口薫，田力正好，石黒聡士，佐藤善輝
Y. Ota and N. Matta

The 1999 earthquake fault andTerr.

Atmos.

its repeated occurrence at theOcean

21 ，Earthquake Museum

737-741

平成22年(f) 特許出願，ソフトウエア開発，仕様・標準等の策定
1) 特許出願
なし
2) ソフトウエア開発
なし
3) 仕様・標準等の策定なし(3)

平成２３年度業務計画案

平 成 2 0 年 度 に 新 潟 県 上 越 市・南 魚 沼 市 な ど に 設 置 し た 観 測 点 4 9 ヶ 所 に お い て 、キ ャ ン ペ ー ン 観 測 を 1 回 実 施 す る 。 さ ら に 、 GPS 観 測 デ ー タ の 解 析 を 行 い 、 4 年
間 の 解 析 結 果 か ら 地 殻 変 動 分 布 の 概 要 を ま と め る と と も に 、そ の 結 果 に 基 づ い て
変形 過程 のモ デル を改 良す る。また 、研 究成 果を 論文 にま とめ ，学 術誌 に 投 稿 する と と も に 、 研 究 成 果 を 日 本 地 球 惑 星 科 学 連 合 2010 年 大 会 、 日 本 地 震 学 会 秋 季
大 会 、ア ジ ア オ セ ア ニ ア 地 球 科 学 会 、米 国 地 球 物 理 学 連 合 秋 季 大 会 等 の 学 会 で 発表する。また、これらの学会で関連する研究の情報収集を行う。
な お 、 年 度 の 早 い 時 期 に 一 部 の 観 測 点 で GP S 観 測 を 実 施 し 、 東 北 地 方 太 平 洋
沖 地 震 の 影 響 を 把 握 し 、そ の 結 果 に 基 づ い て 、必 要 で あ れ ば 研 究 方 針 の 見 直 し を行う。

268