東海地方の地殻変動
6一2

国土地理院
Geospatial Information Authority of Japan [水準測量　森～掛川～御前崎]第1〜5図は東海地方(森町〜御前崎市間)の水準測量結果である.前回の観測は半年前の2016年1月である.第1図の最上段は,最新の観測結果と前回の観測結果の差による各水準点の上下変動である.森町側に対して御前崎市側が沈降する長期的な傾向に特に変化は見られない.第2図は，掛川市（140-1）から見た御前崎市（2595）の上下変動時系列である．上のプロットが生の観測値による時系列，下のプロットが年周成分を除去した後の時系列である． 2000年夏以前のSSE開始よりも前の沈降の速度と比較して， SSE進行期にある2000年秋頃から2005年夏頃までは沈降速度が速かったが，2005年夏以降は，2000年夏よりも前の沈降速度にほぼ戻ったように見える.第3図は，第2図の観測結果について，最新の変動が従来のトレンド（傾き）上にのっているかどうか等を，できるだけ定量的に評価するための資料である． 2000年秋～2005年夏のSSE進行期とその前後の期間の3つの期間に分けて，トレンドを推定した後，年周成分を推定した.上段の時系列は，第2図の年周成分を除去していない時系列のうち1995年以降のものである.破線は，3期間に分けて推定した回帰曲線である．2段目の表に回帰モデルの数値を示した．期間（2）のSSE進行期は，傾きが約-8mm/年と沈降速度が速くなったが，その後の期間（3）については約-5mm/年と期間（1）の沈降速度に近くなっている．年周変化の振幅は，小さくなっている．同時に，回帰の標準偏差も小さくなっている．一番下の段に，期間（2）から期間（3）にかけての時系列の拡大図を示した．回帰モデルからの残差による標準偏差を細い破線で示してある．長期的な傾向に特段の変化は見られない.第4図は，森町（5268）を基準とした掛川市（140-1）と御前崎市（2595）の変動時系列グラフである.森町に対する掛川市および御前崎市の長期的な沈降傾向に特段の変化は見られない.
[水準測量　御前崎　時系列]
第5図は，掛川から御前崎検潮所に至る各水準点の時系列上下変化である．御前崎検潮所付属水準点は，2009年8月駿河湾の地震時に局所的に沈下したものと考えられる．2011年4月に御前崎先端付近でわずかな降起の傾向が見られたが,その後は従来とほぼ同じトレンドで沈降している.
[GNSSと水準測量の比較　御前崎]
第6図は，電子基準点間の比高変化について，水準測量（取付観測という）の結果とGNSS連続観測結果とを比較したものである．両者はほぼ同様の傾向を示しており，最新のデータは従来の長期的な沈降傾向に沿っている.
[GNSS 上下高精度比高観測御前崎]
第7〜10図は，掛川〜御前崎間における高精度比高観測（GNSS連続観測）の結果である.第7図は，高精度比高観測点間の比高変化について，水準測量の結果とGNSS連続観測結果を比較したものである.両者はほぼ同様の傾向を示しており，最新のデータは従来の長期的な沈降傾向に沿っている.なお, (2)のＨ下板沢－Ｈ下朝比奈２の基線では最近の数ヶ月間は隆起が見られていたが，2016年9月6日にＨ下朝比奈２観測点周辺の樹木を伐採したことにより、元に戻っている.第8〜10図は、掛川〜御前崎間における高精度比高観測（GNSS連続観測）の結果である.第8図下段の比高変化グラフには約2cmのばらつきが見られるが， H下板沢（98H023）観測点に対して御前崎側の観測点が長期的な沈降の傾向にあることが見てとれる．なお，第7図と同様， (3)H下板沢－ H下朝比奈2の基線で2016年4月以降見られていた隆起は， 9月6日のＨ下朝比奈2観測点周辺の樹木伐採により，元に戻っている.第9図に，各高精度比高観測点のH下板沢観測点に対する比高変化について，1ヶ月ごと及び10日ごとの平均値を示す.各図の右に各点の上下変動速度(マイナスは沈降)が記されている. H下朝比奈2 （98H010）を除き，特段の傾向の変化は見られない．第10図左は， H細谷（98H025）観測点に対する各点の比高について， 1ヶ月平均値と3ヶ月前の1ヶ月平均値との差を，最近3年間について示したものである． H下朝比奈２を除き，特段の傾向の変化は見られない.第10図右は，同様にＨ細谷観測点に対する各点の比高の1ヶ月平均値の前月との差を，最近1年間について示したものである．特段の傾向の変化は見られない．
[水準測量　御前崎先端部]
第11〜12図は，御前崎先端部の変動を見るために小さな環で行っている水準測量の結果である．最近は概ね半年に1回の頻度で実施している．第11図の最上段は，今回2016年7月の最新の結果と前回2016年1月の結果の差による上下変動観測結果で，特段の変化は見られない.第12図は，網平均を行った結果を最近の4つの期間について示したもので，比較のため，1977年からの上下変動の累積を比較のために最下段に示す. （4）に示した最近の短期的な傾向には，従来の傾向と特段異なる変化は見られない.
[水準測量　静岡県菊川市]
第13～15図は静岡県が実施している菊川市付近の水準測量の結果である.平成26年度からは観測の頻度がこれまでの2週間に1度から1ヶ月に1度になっており，最新のデータは2016年10月期に行われた観測結果となっている．グラフの掲載順序は，第13図に東側の路線，第14図に西側の路線の結果を掲載し，各図の最上段に，一番長い路線の結果を示してある.第13図,第14図の各図の中段にはSF2129から2602-1に至る南北の短い路線（約100m）のデータが掲載されているが，これらは独立な観測値による結果である．両者とも，2602-1で2009年8月11日の駿河湾の地震時に1mmを超える沈降を示した.また，第14図の最上段にも変化が見られることから，1033も同時に沈降した可能性がある．これら，2602-1と10333を含むグラフの近似曲線は， 2009年8月11日の駿河湾の地震前までのデータを用いて計算した．東北地方太平洋沖地震による影響や顕著な傾向の変化は見られない.第15図に示した傾斜ベクトルの時間変化には，揺らぎを伴いながらも，全体としては長期的な南南東傾斜の傾向が見られる.
[GNSS　御前崎とその周辺]
第16～20図は御前崎とその周辺のGNSS連続観測結果である．三ヶ日から榛原（はいばら）に至る東西方向の基線もあわせて示している.第17図の(4)(5)において2009年の夏に見られる跳びは，2009年8月11日に発生した駿河湾の地農に伴う御前崎Ａ観測点の変動によるものである． 2011年3月11日に発生した東北地方太平洋沖地震に伴い、第17図の(5)に地震時と地震後の基線の短縮が見られる.第18図の(8)において2009年8月頃から，掛川観測点が東向きに動いたような基線長の変化が見られた後、10月に戻った．同様の変化はピラーに内蔵された傾斜計にも見られるが， GNSSの上下成分には見られない． 2010年夏にも同様の東向きの変化が見られた後， 9月28日以降戻っている. 2009年も2010年も大雨後に戻っているが,原因は不明のままである. 2011年及び2012年にはこのような変化はなかったが，2013年以降，再び同様の変化が見られるようになった.なお(6)において， 2016年4月頃に見られる飛びは，アンテナ更新によるものである.第20図の(6)において2014年6月頃から隆起する向きの変化が見られたが， 8月に観測点周辺の樹木を伐採した後に戻っており,観測点周辺の樹木の成長に伴う受信環境の悪化による誤差であった可能性がある．また，2016年6月頃からも隆起する向きの変化が見られており，同様に樹木繁茂による影響の可能性がある.一部の観測点では2010年2～3月頃にレドームの開閉を行ったことによる見かけ上の変動が含まれている場合があるので,第16図下段の観測局情報を参照する必要がある.最近のデータには，特段の傾向の変化は見られない.
[GNSS 験河湾]
第21～23図は，駿河湾とその周辺のGNSS連続観測時系列である．2014年4月頃から，焼津Ａ (990840）が隆起・沈降を繰り返しているような変化が見られるが，原因については電波干渉の影響の可能性があり，この影響を除去するために2015年6月4日にアンテナ交換とフィルター挿入を行った.その後は隆起･沈降を繰り返す変化は見られていない.その他には，傾向に特段の変化は見られない.
[長距離水管傾斜計　御前崎・切山]
第24図は，御前崎長距離水管傾斜計の月平均結果と傾斜計端点間の水準測量結果である．長距離水管傾斜計のデータは, 2012年8月14日から10月18日までの間の機器異常による欠測と2013年1月28日から2月1日までの間に行われた機器交換の前後で変化がないものと仮定してデータをつなげている．水準測量結果では，長期的な東側隆起の傾向が見られる．上側口印の水準測量のデータ，下側の○印の水管傾斜計のデータともに2009年8月11日の駿河湾の地震時の跳び等を補正して表示している.下側の○印の水管傾斜計のデータでは, 2009年6月17日に西側局舎にエアコンを設置した効果により，最近のプロットの年周成分は小さめである．なお,東側局舎へのエアコン設置は1993年4月で，2002年の冬に行われた両局舎の建て替えによって密閉性が高まったとの記録がある. 2014年の1月頃から3月頃にかけて，水管傾斜計によって測定される傾斜の値がやや小さめであったが，その後，ほぼ元のレベルに戻っている.第25図は御前崎及び切山の長距離水管傾斜計観測値の日平均値データ及び時間平均値データである．特段の傾向の変化は見られない．
[深井戸　ひずみ　御前崎]
第26図は御前崎の深さ約800mの深井戸で実施している地殻変動（ひずみ）連続観測結果である．特段の傾向の変化は見られない．なお，同地点で観測を行っていた傾斜計は，1成分が7月23日に大きな変化を示し、機器異常と見られるため、今回から掲載を取り止めている.
[東海地方の地殻変動]
第27〜31図は，白鳥(しろとり)観測点を固定局として示した，東海地方の地殻変動である.第27図上段は最近の1年間の水平変動である．比較のために，東北地方太平洋沖地震前においてスロースリップのなかった2つの時期における変動速度を中段に，スロースリップが発生していた時期の変動速度を下段に示している.最近の東海地方の地殻変動には西向きの変動が広く見られ，スロースリップの発生していなかった時期のものに近い.第28図は，上下成分について同様の比較を示すものである．水平よりもばらつきが大きい.第29～30図は，東北地方太平洋沖地震前の2008年1月～2011年1月の期間の変動を定常変動と仮定し，それからの変動の差を非定常変動として示した図である．水平成分及び上下成分のそれぞれについて，最近の約1年間の図と3ヶ月ごとの図を示す。最近の約1年間の図では，明瞭な南東向きの変動は見られない．なお，浜名湖の東側のわずかな非定常的な隆起は引き続き見られるものの、小さくなってきている.第31図は，東海地方のGNSS連続観測点の非定常地殻変動の3成分時系列である．東北地方太平洋沖地震の余効変動の影響は小さくなってきている. (7)の榛原観測点で2016年6月頃から上下成分に見られる変化は，観測点周辺の樹木の成長に伴う受信環境の悪化による誤差の可能性がある.
[東海地方　ゆっくり滑りによるプレート境界面上の滑り分布]
第32～58図は，東海地方に見られる非定常的な地殻変動について，プレート境界面上の滑りを推定した結果である．解析方法や解析期間，形状が異なる複数のモデルで推定を行っており，概要を第32図にまとめている．なお， 2016年熊本地震に伴い固定局の三隅観測点が変動したため，これらの資料では地震時の変動を補正している.
(モデル1－1, 1－2)
第33～42図は，東海地方の非定常地殻変動（2008年1月～2011年1月を基準）からまず東北地方太平洋沖地震による余効変動の影響を除去し，その結果を用いてプレート境界面上の滑りを推定した結果である．第33～37図は解析期間が2013年3月12日～2016年10月15日，第38～42図は，最近の1年間（2015年10月15日～2016年10月15日）である.第33～34図（第38～39図）は，東北地方太平洋沖地震の余効滑りの推定結果と，それから計算される東海地方での余効変動，非定常地殻変動と，余効変動の影響を除去して得られた非定常成分である．余効滑りによる変動を除去した図では，浜名湖北西で南東向きの変動，浜名湖の東側で東向きの変動と隆起が見られる.第34～35図（第39～40図）は，非定常地殻変動から一様滑りの矩形断層を推定した結果である．また，第36～37図（第41～42図）は，非定常地殻変動からプレート境界面上での滑り分布を推定した結果である．浜名湖付近に滑りの中心が推定されている．第37図（第42図）の残差を見ると，水平変動はよく説明できているが，上下変動については説明できていない変動が残っているように見える.
(モデル2-1)
第43〜50図は，東海地方の非定常地殻変動について，時間依存インバージョンにより，プレート境界での滑り分布を時間発展で推定した結果である．このモデルでは，東北地方太平洋沖地震の余効滑りとフィリピン海プレート上面での滑りを同時に推定している.なお、この資料でいう非定常地殻変動とは，時系列データから定常的な変動と周期的な変動を除去したものである．周期的な変動については，2012年1月以降の全てのデータを用いて，多項式曲線と2013年1月以前と以後でそれぞれの周期成分（年周・半年周）を同時に推定した．定常的な変動については，2008年1月～2011年1月での平均変動速度を推定した．また，平滑化のために3日間の値の平均をとっている.第43図は，東海地方の非定常地殻変動の最近の約1年間とその3ヶ月ごとの水平変動の図である．図中に示す日付を中心として，前後4日間での9日間の平均値をとり，それぞれを比較している．東北地方太平洋沖地震の余効変動の影響が見られる．3ヶ月ごとの図及び最近の約1年間の図のいずれでも，浜名湖の北西周辺でわずかに南向きの成分が見られるが，最近の3ヶ月間ではかなり小さくなってきている.第44図は，非定常地殻変動の三成分時系列をプロットしたものである.第45図は，非定常地殻変動から東北地方太平洋沖地震の余笏滑りを計算により除去した図である．2014年10月以降について，6ヶ月間ごとの水平変動を示している．全ての時期で，浜名湖の北西では南東方向の変動，浜名湖の東側では東方向への変動が見られるが， (D)の2016年4月以降では変動が小さくなっている．なお，この図を含め，第50図までのベクトル図では，2日毎にデータを抽出し，2日前及び2日後との3点での平均値を用いている．第46図は，推定された滑りの時間変化を示す資料である．上段は滑り分布，中段は誤差分布，下段はモーメントの時間変化グラフである． 2013年半ばから一貫してモーメントの解放量が増加傾向にあったが，最近はその傾向が鈍化しているように見える．なお，2016年1月頃に見られるモーメントの急激な増加は, 2016年1月上旬に伊勢湾で発生した短期的SSEによる影響である.また, 2016年夏頃に見られる増加は，年周成分が完全には除去できていない影響の可能性がある．第47〜48図は，6ヶ月間ごとの水平及び上下変動の計算値と観測値を比較したものである．水平変動は概ね説明できている．上下変動についてはばらつきが大きいため，判断できない．第49図は，滑り領域とモーメントの時間変化について，今回と前回（2001～2005年）とを比較したものである．今回の滑りは最大で7cm程度であり，滑りの中心は前回よりも南側に位置している.第50図は、東北地方太平洋沖地震の余効変動を除去して得られた非定常変動とフィリピン海プレート上面の滑り分布を対応させて示した資料である.
(モデル2ー2)
第51～58図は，第43～50図と同様に，時間依存インバージョンによりプレート境界での滑り分布を時間発展で推定した結果である．長期的SSEの進行状況をより詳細に確認するため，2016年1月上旬に伊勢湾で発生した短期的SSEによる地殻変動を除去して推定を行っている．第54図のモーメントの時間変化グラフを見ると，2015年半ば以降は，モーメント解放速度がそれ以前に比べて小さくなっていることが分かる.
 197
菊川市
1955 999須前崎市2594 
 1950-7 -牧之原市御前崎市(1) 2016/01 — 2016/06~07掛川市cm 
 25 1 0 11 –2 cm
 255 1 0 (2) 2015/07 — 2016/01 11 21 cm 0 -1 cm 1 o -1 cm 1 o –1 (5) 2014/01 － 2014/06〜07標高断面図m 150
100 50 0 ※ 本図は地理院地図を使用して作成し
     ている。なお、海域部は海上保安庁
     海洋情報部の資料を使用している。水準測量による森町～掛川市～御前崎市間における上下変動第1図Fig. 1 Vertical crustal movement from Mori town to Omaezaki city via Kakegawa city. (4) 2014/06〜07 — 2015/01 (3) 2015/01 — 2015/07為町● :水準点
O :固定鋲
水準点140-1（掛川市）を基準とした2595（御前崎市浜岡）の高さの経年変化第2図Fig. 2  Time series of height change of BM2595 (Hamaoka) as referred to BM140-1 (Kakegawa).
•各期間の回帰モデル（1次トレンド＋年周）を実線で表示している. ・回帰モデルからの残差による標準偏差を破線で示している.第３図　水準点2595（御前崎市）の経年変化　スロースリップイベント期間で分けた回帰モデルFig. 3 Regression model for the period before, during and after the slow slip event.水準点5268（森）を基準とした140-1（掛川）と2595（御前崎市浜岡）の上下変動時系列第4図Fig. 4 Time series of height change from BM 5268 (Mori) to BM140-1 (Kakegawa) and BM2595 (Hamaoka). • 1980年11月水準点2602-1を新設第５図水準点140-1（掛川市）を基準とした掛川～御前崎間の各水準点の高さの経年変化Fig. 5 Time series of height changes of benchmarks along the route between Kakegawa and Omaezaki as referred to BM140-1 (Kakegawa).
律前崎　電子基準点の上下変動

従来の傾向に変化は見られない
〈1〉静岡線(93089) -熱川(93052)辻高1925　2006　2007　2008 2003   2010   2011נו201 2014 1000　2001　2002　2003　2004 2012 1822
10/01〜10/08の平均GNSS連続報測による結果については，水準測量の会期間との差が最小となるオフセット量を求めている

水準測量とGNSS連続規測
御前崎　高精度比高観測点の上下変動
● :水準測量O : GNSS連続観測•水準測量による結果は、最初のプロット系の値を0cmとして描画している. GNSS連続観測のプロット点は、日々の産標値〈HTI :最終解）から計算した値の月平均値 最新のブロット点は10/01～10/08の平均• GNSS連続報測による結果については、水準測量の全期間との差が最小となるオフセット量を求めている
御前崎地域の高精度比高観測GNSS観測結果（基線図及び保守状況） Results of high precision vertical GNSS measurements in Omaezaki region (Baseline map and history of maintenance). -[HTI :最終解]御前崎地域の高精度比高親測GNSS親測結果第８図b Fig. 8b Time series of the height change in precision vertical GNSS measurements sites in Omaezaki region.
最新のブロット点は、月平均値は10/01〜10/08， 10日間平均値は10/08〜10/08の平均.
最新のブロット点は、 ・月平均値は、2011年3月11日発生した東北地方太平洋沖地震前後で期間を分けて回帰モデルを作成している.御前崎地域の高精度比高観測GNSS観測結果（1ヶ月間及び10日間移動平均・時系列）第9図Fig. 9 Results of high precision vertical GNSS measurements in Omaezaki region (Time series of 1 month and 10 days running mean).
98H025 98H023 ３か月98H016 98H025 98H023 1か月98H016 98H010 98H007 98H004
固定局：98H025御前崎地域の高精度比高観測GNSS観測結果（点毎の３ヶ月間及び１ヶ月間の変動量）第10図Fig. 10 Results of high precision vertical GNSS measurements in Omaezaki region (Height change every three months and 1 month at each site).
御前崎地方の上下変動（1 ）
水準測量による御前崎先端部の上下変動（1）第11図Fig. 11 Vertical crustal deformation by the precise leveling survey around Omaezaki (1/2).頃向に変化は見られない。
水準測量による御前崎先端部の上下変動（2）第12図Fig. 12 Vertical crustal deformation by the precise leveling survey around Omaezaki (2/2).
菊川市付近の水準測量結果（1）
水準点2602-1と2601の経年変化
•水準点2602-1 &含むグラフは、2009/08/11酸河海の地震でステップが生じたため、地震前までのデータについて近似由線を計所有害者になっており、そのことを目標に取り組織で入り替えた。
算し年間成分を構正している。地震後のデータは年間成分を構正していない。 •水準点S F2 ! 29は20 ! 6  4 7 月に移転した

菊川市付近の水準測量結果（2）
水準点2602-1と10333の経年変化
静岡県による短距離水準測量結果（3）第15図Results of short distance leveling (3): Vector representations of tilt derived from leveling data in (1) and (2). Original Fig. 15 data are provided by the Prefectural Government of Shizuoka.御前崎周辺ＧＮＳＳ連続観測時系列(1)
※2003年3月5日に基準局92110(つくば1)のアンテナおよびレドームの交換を実施し、解析値に補正をしています。御前崎周辺GNSS連続観測点観測結果（基線図及び保守状況）第16図Fig. 16 Results of continuous GNSS measurements in the Omaezaki region (baseline map and maintenance history ).
御前崎周辺ＧＮＳＳ連続観測時系列(2)
ふ線変化グラフ揭覧： 1996/04/01～2016/10/22 JST
基様変化グラフ800: 2014/10/01～2016/10/22 JST

御前崎周辺ＧＮＳＳ連続観測時系列(3)
あ緑変化グラフ期間: 1995/04/01～2016/10/22 JST基線変化グラフ探閲： 2014/10/01～2016/10/22 JST
御前崎周辺GNSS連続観測点観測結果（斜距離） Results of continuous GNSS measurements in the Omaezaki region (baseline length) (2/2).
御前崎周辺ＧＮＳＳ連続観測時系列(4)
比高変化グラフ揭稿：1996/04/01～2016/10/22 JST比高変化グラフ800: 2014/10/01～2016/10/22 JST
比高変化グラフ規間： 1956/04/01～2016/10/22 JST比高変化グラフ期間: 2014/10/01～2016/10/22 JST
第20図 御前崎周辺GNSS連続観測点観測結果（比高）

御前崎周辺ＧＮＳＳ連続観測時系列(5)

駿河湾周辺GNSS連続観測点観測結果（基線図及び保守状況）第21図Fig. 21 Results of continuous GNSS measurements around the Suruga Bay (Baseline map and maintenance history).験河湾周辺   ＧＮＳＳ連続観測時系列（2）あ線変化グラフ独身の産化は見られない。
第22図 駿河湾周辺GNSS連続観測点観測結果（斜距離）

駿河湾周辺   ＧＮＳＳ連続観測時系列(3)比高変化グラフ第23図 駿河湾周辺GNSS連続観測点観測結果（比高）
(relative height).
※３ 2009年6月西側局舎にエアコンを設置※ 4  2009年8月11日駿河湾の地震による跳びを補正. ※ 5  2009年8～10月センサー感度不良のため欠測.  2009年8月11日に発生した駿河湾の地震前後のデータに連続性は無い※ 6  2012年8月14～10月18日機器異常のため欠測，期間前後のデータに連続性は無い第24図 御前崎長距離水管傾斜計による傾斜観測結果Fig. 24 Results of tilt observation measured by long water tube tiltmeter at Omaezaki.第25図御前崎長距離水管傾斜計及び切山長距離水管傾斜計による傾斜観測結果Fig.25 Results of tilt observation measured by long water tube tiltmeter at Omaezaki and Kiriyama.第26図御前崎地中地殻活動監視装置による連続観測結果（ひずみ日平均値および水平ひずみ） Fig. 26 Results of continuous measurements of tilt and strain in the Omaezaki deep borehole (Daily mean value and horizontal strain ).
(観測点名:御前崎下岬)
動速度（白鳥固定） Horizontal deformation of recent 1 year in the Tokai district based on GNSS measurements and horizontal Fig. 27 deformation rates before (middle left), during (lower) and after (middle right) the Tokai slow slip (fixed Shirotori).
GNSS観測による東海地方の最近１年間の上下変動及びスロースリップ開始前・進行期・終息後の上下変第28図動速度（白島固定） Vertical deformation of recent 1 year in the Tokai district based on GNSS measurements and vertical deformation Fig. 28 rates before (middle left), during (lower) and after (middle right) the Tokai slow slip (fixed Shirotori). 2008年1月〜2011年1月のデータから平均変動速度、年周/半年周成分を推定して、元の時系列データから除去した時系列。 GNSS 観測による東海地方の最近１年間と３ヶ月ごとの非定常地殻変動（水平変動）第29図Fig. 29 Transient horizontal deformation of recent 1 year and every 3 months in the Tokai district. . GNSS 観測による東海地方の最近１年間と３ヶ月ごとの非定常地殼変動（上下変動）第30図Fig. 30 Transient vertical deformation of recent 1 year and every 3 months in the Tokai district.東海地方の非定常地殻変動時系列第31図Fig. 31 Time series of transient deformation at selected stations in the Tokai district.
※三隅には2016年4月の熊本地震に伴う地殻変動の補正を行った.第32図 東海地方ゆっくり滑りによるプレート境界面上の滑り分布 モデル一覧Fig. 32  List of models of slip distribution on the plate interface of Tokai district.
東北地方太平洋沖地震の余効すべりの推定（暫定）【モデル1－1】
余効すべりモデルから計算される東海地方の地殻変動（暫定）
第33図東北地方太平洋沖地震の余効すべりの推定（暫定）（モデル1－1） Fig. 33 Estimated postseismic slip of the 2011 off the Pacific coast of Tohoku Earthquake (preliminary result) (model 1-1). (Upper) Slip distribution (left), horizontal deformation (center), vertical deformation (right). (Lower) Comparison of calculated and observed crustal deformation of Tokai district (preliminary result). Horizontal deformation (left), vertical deformation (right).
余効変動除去前の矩形断層モデルの推定（暫定)【モデル1ー1】
余効変動除去前の矩形断層モデルの推定（暫定）（モデル1－1）第34回Fig. 34 Estimated rectangular fault model before eliminating postseismic deformation (preliminary result) (model 1-1). (1) Observed transient deformation. Horizontal deformation (left), vertical deformation (right). (Top) (Middle) (2) Estimated rectangular fault model and comparison of transient deformation and calculated deformation. (Bottom) (3) Residual between observed transient and calculated deformation form the rectangular fault model.余効変動除去後の矩形断層モデルの推定（暫定)【モデル1－1】水平変動余効変動除去後の矩形断層モデルの推定（暫定）（モデル1－1）第35区Fig. 35 Estimated rectangular fault model after eliminating postseismic deformation (preliminary result) (model 1-1). (1) Observed transient deformation. Horizontal deformation (left), vertical deformation (right). (Top) (Middle) (2) Estimated rectangular fault model and comparison of transient deformation and calculated deformation. (Bottom) (3) Residual between observed transient and calculated deformation from the rectangular fault model.余効変動除去前のすべり分布モデルの推定（暫定）（モデル1－1）第36図Fig. 36 Estimated slip distribution on the plate interface before eliminating postseismic deformation (preliminary result) (model 1-1). (Top)       (1) Slip distribution. (Middle) (2) Comparisons of transient deformation and calculated deformation. Horizontal deformation (left), vertical deformation (right). (Bottom) (3) Residual between observed transient and calculated deformation from the slip distribution model. Horizontal deformation (left), vertical deformation (right).
余効変動除去後のすべり分布モデルの推定（暫定)【モデル1－1】
第37図余効変動除去後のすべり分布モデルの推定（暫定）（モデル1－1） Fig. 37 Estimated slip distribution on the plate interface after eliminating postseismic deformation (preliminary result) (model 1-1). (Top) (1) Slip distribution. (Middle) (2) Comparisons of transient deformation and calculated deformation. Horizontal deformation (left), vertical deformation (right). (Bottom) (3) Residual between observed transient and calculated deformation from the slip distribution model. Horizontal deformation (left), vertical deformation (right).
東北地方太平洋沖地震の余効すべりの推定（暫定）【モデル1—2】
第38図東北地方太平洋沖地震の余効すべりの推定（暫定）（モデル1－2） Fig. 38 Estimated postseismic slip of the 2011 off the Pacific coast of Tohoku Earthquake (preliminary result) (model 1-2). (Upper) Slip distribution (left), horizontal deformation (center), vertical deformation (right). (Lower) Comparison of calculated and observed crustal deformation of Tokai district (preliminary result). Horizontal deformation (left), vertical deformation (right).
余効変動除去前の矩形断層モデルの推定（暫定)【モデル1ー2】
余効変動除去前の矩形断層モデルの推定（暫定）（モデル1－2）第39回Fig. 39 Estimated rectangular fault model before eliminating postseismic deformation (preliminary result) (model 1-2). (1) Observed transient deformation. Horizontal deformation (left), vertical deformation (right). (Top) (Middle)   (2) Estimated rectangular fault model and comparison of transient deformation and calculated deformation. (Bottom) (3) Residual between observed transient and calculated deformation form the rectangular fault model.
余効変動除去後の矩形断層モデルの推定（暫定)【モデル1ー2】
余効変動除去後の矩形断層モデルの推定（暫定）（モデル1－2）第40回Fig. 40 Estimated rectangular fault model after eliminating postseismic deformation (preliminary result) (model 1-2). (1) Observed transient deformation. Horizontal deformation (left), vertical deformation (right). (Top) (Middle) (2) Estimated rectangular fault model and comparison of transient deformation and calculated deformation. (Bottom) (3) Residual between observed transient and calculated deformation from the rectangular fault model.第41図余効変動除去前のすべり分布モデルの推定（暫定）（モデル1－2) Fig. 41 Estimated slip distribution on the plate interface before eliminating postseismic deformation (preliminary result) (model 1-2). (1) Slip distribution. (Top) (Middle)  (2) Comparisons of transient deformation and calculated deformation. Horizontal deformation (left), vertical deformation (right). (Bottom)  (3) Residual between observed transient and calculated deformation from the slip distribution model. Horizontal deformation (left), vertical deformation (right).余効変動除去後のすべり分布モデルの推定（暫定）（モデル1－2）第42回Fig. 42 Estimated slip distribution on the plate interface after eliminating postseismic deformation (preliminary result) (model 1-2). (1) Slip distribution. (Top) (Middle)   (2) Comparisons of transient deformation and calculated deformation. Horizontal deformation (left), vertical deformation (right). (Bottom)  (3) Residual between observed transient and calculated deformation from the slip distribution model. Horizontal deformation (left), vertical deformation (right).東海地方の非定常水平地殻変動  [固定局：三隅] 【モデル2-1】項式曲線と2013年1月以前と以後でそれぞれ周期成分（年周・半年周）を同時に推定している。一次トレンドは, 2008年1月～2011年1月のデータから推定している. •カルマンフィルターで平滑化した値を示している. •三隅には2016年4月の新本地震に伴う地殻変動の擁正を行った.
•非定常地設変動とは，時系列データから周期成分と一次トレンドを除去したもの．2012年以降のデータに対して，多項式曲線と2013年1月以前と以後でそれぞれ周期成分（年周•半年周）を同時に推定している。一次トレンドは、 2008年1月～2011年1月のデータから推定している. • 3日移動平均の慎を示している•三隅には2016年4月の熊本地震に伴う地設変動の補正を行った
東海地方の非定常地設変動時系列 [固定局：三隅]【モデル2-1】
東海地方の非定常的な地殻変動（余効変動除去後）【モデル2-1】
• GNSS連続観測の結果から非定常地設変動の時系列データを作成し、それをカルマンフィルターで平滑化したデータから.格定された東北地方太平洋沖地震後の太平洋ブレート上面における余効すべりの影響を取り除いたものを示している•非定常地設変動とは，時系列データから周期成分と一次トレンドを除去したもの．2012年以降のデータに対して、多項式曲線と2013年1月以前と以後でそれぞれ周期成分（年周・辛年周）を同時に推定している。一次トレンドは、2008年1月～ 2011年1月のデータから推定している。 •調査情円は、非定常地設定動の時系列をカルマンフィルターで平滑化してそれから推定される調査と東北地方太平洋冷地震後の太平洋ブレート上面における余効すべりのモデルから計算される地設変動の設差から計算している.
第45図　東海地方の非定常的な地殻変動（余効変動除去後）

時間依存インバージョンによるすべり分布モデルとモーメントの時間変化【モデル2-1】
• ONSS連続観測の結果から非定常地發変動の時系列データを作成し、時間依存のインバージョンでフィリピン海プレート及び太平洋プレート上面のすべりを同時に推定している. •非定常地殺変動とは、時系列データから開幕成分と一次トレンドを除去したもの．2012年以降のデータに対して、多項式曲線と2013年1月以前と以後でそれぞれ周期成分〈年間・半年間〉を同時に推定している．一次トレンドは、2008年1月～2011年1月のデータから推定している。 •推定されたフィリビン海ブレート上面のすべり量とその設差を.上段と中段にそれぞれ増間ごとに示している.

東海地方の地殻変動の観測値と計算値の比較（水平変動）【モデル2-1】
•上段は, GNSS連続観測の結果から非定常地設変動の時系列データを作成し,それをカルマンフィルターで平滑化したデータから、雑定された東北地方太平洋沖地震後の太平洋プレート上面における余効すべりの影響を取り除いたものを示している. •未定常地設変動とは、時系列データから周期成分と一次トレンドを除去したもの、2012年以降のデータに対して、多項式曲線と2013年1月以前と以後でそれぞれ周期成分（年限•半年周）を同時に推定している．一次トレンドは、2008年1月～ 2011年1月のデータから推定している。 •中段は、推定した東海地方のザベリ分布モデルから計算される水平地設定動を示している. •下段は、上段から中段を差し引いた残差を示している. •上段は、ONSS連続観測の結果から非定常地投資動の時系列データを作成し、それをカルマンフィルターで平滑化したデータから、雑定された東北地方太平洋沖地震後の太平洋ブレート上面における余効すべりの影響を取り除いたものを示している. •未定常地殻変動とは、時系列データから周期成分と一次トレンドを除去したもの。 2012年以降のデータに対して、多項式曲線と2013年1月以前と以後でそれぞれ周期成分〈年周・半年周〉を同時に推定している。一次トレンドは、2008年1月～ 2011年1月のデータから格定している. •中段は、推定した薬落地方のすべり分考モデルから計算される上下地殻変動を示している. •下段は,上段から中段を差し引いた残差を示している.
東海地方の地殻変動の観測値と計算値の比較（水平変動）
第47図

Fig. 47 Comparisons between observed and calculated horizontal displacements in the Tokai district (model 2-1). Observation (top), calculation (middle), residual (bottom).
東海地方の地殻変動の観測値と計算値の比較（上下変動）
第48図
(モデル2ー1)
Fig. 48 Comparisons between observed and calculated vertical displacements in the Tokai district (model 2-1). Observation (top), calculation (middle), residual (bottom).
東海地方の地殻変動の観測値と計算値の比較（上下変動）【モデル2-1】
非定常地設変動から推定される累積すべり分布及びモーメントの時間変化（暫定）
【モデル2-1】
等49区
非定常地殻変動から推定される界積すべり分布及びモーメントの時間変化（暫定）

Comparisons of accumulated slip distribution and time series of estimated moment between the 2014-2015 event and the last event (2001-2005) (preliminary result) (model 2-1). Slip distribution (top), time series of estimated moment (bottom).
東海地方の非定常地殻変動及びすべり分布の時間変化【モデル2-1】
・上段及び中段は、GNSS連続観測の結果から非定常地殻変動の時系列データを作成し、それをカルマンフィルターで平滑化したデータから、推定された東北地方太平洋沖地震後の太平洋プレート上面における余効すべりの影響を取り除いたものを示している。 ・非定常地投資動とは、時系列データから周期成分と一次トレンドを除去したもの. 2012年以降のデータに対して、多項式曲線と2013年1月以前と以後でそれぞれ国規成分（年国・半年国）を同時に推定している。一次トレレンドは、2008年1月～ 2011年1月のデータから推定している。 ・時間依存のインバージョンでフィリピン海プレート及び太平洋プレート上面のすべりを同時に推定した．推定されたフィリピン海プレート上面のすべり量を下段に示している.第50図
東海地方の非定常地殻変動及びすべり分布の時間変化
(モデル2ー1)
Fig. 50 Time evolution of transient displacement and slip distribution in the Tokai district (model 2-1). Horizontal deformation (top), vertical deformation (middle), slip distribution (bottom).
•非定常地發変動とは、時系列データから南期成分と一次トレンドを発去したもの. 2012年以降のデータに対して,多項式曲線と2013年1月以前と以後でそれぞれ恩精成分（年間・半年間）を同時に推定している。一次トレンドは、 2008年1月～2011年1月のデータから推定している. • 3日移動平均の植を示している. • 2016年初頭の短期的SSE前後の地設変動のオフセットを取り除いている. •三隅には2016年4月の熊本地震に伴う地設変動の雑正を行った
東海地方の非定常地殻変動時系列  [固定局：三隅] 【モデル2-2】   オフセット除去版東海地方の非定常的な地殻変動（余効変動除去後）【モデル2–2】
オフセット除去版
• GNSS連続観測の結果から非定常地設変動の時系列データを作成し、それをカルマンフィルターで平滑化したデータから、推定された東北地方太平洋沖地震後の太平洋プレート上面における余効すべりの影響を取り除いたものを示している•非定常地設変動とは,時系列データから閉期成分と一次トレンドを除去したもの. 2012年以降のデータに対して、多項式由線と2013年1月以前と以後でそれぞれ周期成分（年周・半年周）を同時に推定している。一次トレンドは、2008年1月～ 2011年1月のデータから推定している。 •設差特円は、非定常地設定動の時系列をカルマンフィルターで平滑化してそれから推定される設差と東北地方太平洋冷地震後の太平洋ブレート上面における余効すべりのモデルから計算される地紋変動の誤差から計算している. • 2016年初頭の短期的SSE前後の地設変動のオフセットを取り除いている.第53図

時間依存インバージョンによるすべり分布モデルとモーメントの時間変化【モデル2-2】
2016年初頭の短期的SSEに伴う地格変動を除去
すべり分布
平信報(TIM : 2an • GNSS連続観測の結果から非定常地設変動の時系列データを作成し、時間依存のインバージョンでフィリピン海プレート及び太平洋プレート上面のすべりを同時に推定している. •非定常地殻変動とは、時系列データから同期成分と一次トレンドを除去したもの．2012年以降のデータに対して、多項式曲線と2013年1月以前と以後でそれぞれ燭精成分（年間•半年間）を同時に推定している。一次トレンドは、2008年1月～2011年1月のデータから推定している。 •推定されたフィリピン海プレート上面のすべり量とその誤差を、上段と中段にそれぞれ期間ごとに示している. • 2016年初頭の短期的SSE前後の地設変動のオフセットを取り除いている.

東海地方の地殻変動の観測値と計算値の比較（水平変動）【モデル2-2】
•上段は, GNSS連続観測の結果から非定常地殺愛動の時系列データを作成し,それをカルマンフィルターで平滑化したデータから、幾定された東北地方太平洋沖地震後の太平洋プレート上面における余効すべりの影響を取り除いたものを示している. •赤定常地殻変動とは、時系列データから腐期成分と一次トレンドを除去したもの. 2012年以降のデータに対して、多項式曲線と2013年１月以前と以後でそれぞれ周期成分〈年間・半年間〉を同時に推定している。一次トレンドは、2008年１月～ 2011年１月のデータから推定している。 •中段は、推定した東海地方のすべり分布モデルから計算される水平地殻変動を示している. •下段は,上段から中段を差し引いた残差を示している. • 2016年初頭の短期的SSE前後の地設変動のオフセットを取り除いている
東海地方の地殻変動の観測値と計算値の比較（上下変動）【モデル2-2】
•上段は、GNSS連続観測の結果から非定常地設変動の時系列データを作成し、それをカルマンフィルターで平滑化したデータから、雑定された東北地方太平洋沖地震後の太平洋プレート上面における余効すべりの影響を取り除いたものを示している. ·赤定常地設変動とは，時系列データから周期成分と一次トレンドを除去したもの．2012年以降のデータに対して、多項式曲線と2013年1月以前と以後でそれぞれ閲覧成分〈年間・半年間〉を同時に推定している．一次トレンドは、2008年１月～ 2011年１月のデータから推定している。 •中段は、推定した東海地方のすべり分布モデルから計算される上下地紋変動を示している.下段は、上段から中段を差し引いた残差を示している。 • 2016年初頭の短期的SSE解後の地設変動のオフセットを取り除いている.
東海地方の地殻変動の観測値と計算値の比較（水平変動）
第55図
(モデル2-2)
Fig. 55 Comparisons between observed and calculated horizontal displacements in the Tokai district (model 2-2). Observation (top), calculation (middle), residual (bottom).
第56図
東海地方の地殻変動の観測値と計算値の比較（上下変動）
(モデル2-2)
Fig. 56 Comparisons between observed and calculated vertical displacements in the Tokai district (model 2-2). Observation (top), calculation (middle), residual (bottom).
非定常地殺変動から推定される果積すべり分布及びモーメントの時間変化（暫定）
【モデル2-2】 ・黒破線は、沈み込む海側プレート上面の等深線•青破線は、東海地震の想定跟墓域第57図非定常地殻変動から推定される累積すべり分布及びモーメントの時間変化（暫定）（モデル2-2） Fig. 57 Comparisons of accumulated slip distribution and time series of estimated moment between the 2014-2015 event and the last event (2001-2005) (preliminary result) (model 2-2). Slip distribution (top), time series of estimated moment (bottom). •上段及び中段は、GNSS連続戦剤の結果から非定常地設変動の時系列データを作成し、それをカルマンフィルターで平滑化したデータから、推定された東北地方太平洋沖地震後の太平洋プレート上面における余効すべりの影響を取り除いたものを示している. •非定常地設定動とは，時系列データから周期成分と一次トレンドを除去したもの．2012年以降のデータに対して、多項式曲線と2013年１月以前と以後でそれぞれ周期成分（年周・半年周）を同時に推定している。一次トレンドは、2008年1月～ 2011年1月のデータから推定している. •時間依存のインバージョンでフィリビン満プレート及び太平洋ブレート上面のすべりを同時に推定した。推定されたフィリピン海プレート上面のすべり量を下段に示している. • 2016年初頭の短期的SSE前後の地殻変動のオフセットを取り除いている.第58図
東海地方の非定常地殻変動及びすべり分布の時間変化
(モデル2-2)
Fig. 58 Time evolution of transient displacement and slip distribution in the Tokai district (model 2-2). Horizontal deformation (top), vertical deformation (middle), slip distribution (bottom).