技術報告
海洋情報部研究報告 第 52 号 平成 27 年 3 月 2 日

REPORT OF HYDROGRAPHIC AND OCEANOGRAPHIC RESEARCHES No.52 March, 2015

海底地殻変動観測における重心推定法の効果についての考察†
秋山裕平＊ 1，渡邉俊一＊ 2
Study of stability and efficiency of the transponder array constraint method for the GPS acoustic observation†
Yuhei AKIYAMA ＊ 1 and Shun-ichi WATANABE ＊ 2Abstract
For the precise GPS acoustic seafloor positioning, we evaluated the stability and the efficiency of thetransponder array constraint method. Constraining the geometry of the seafloor transponder array via themulti epoch analysis enables us to estimate more precise position of the seafloor sites. The accuracy ofthe positioning should depend on the accuracy of the geometry of the array. In this study, we found that,though the geometry of the array converges to the probable one as the number of epoch used in the multi
epoch analysis increases, the convergence rate strongly depends on the seafloor site. It suggests that weshould estimate the accuracy of the geometry for each site. Moreover, we revealed the strong relationshipbetween the horizontal projected areas of the transponder array and the ver tical deviations of theestimated position. This suggests that the sound speed should be corrected by constraining the geometryof the transponder array. The results will lead to the evaluation of the efficiency of the transponder array
constraint method.
1

はじめになど）や，2011 年東北地方太平洋沖地震（M 9.0）

海上保安庁海洋情報部では，海洋プレートの沈に伴う海底地殻変動（Sato et al., 2011; Watanabe

み込みに伴う陸側プレートの地殻変動を監視するet al., 2014）などを検出している．

ため，GPS−音響測距結合方式による海底地殻変

本観測の技術開発においては，船底へのトラン動観測の技術開発および海底観測点の展開を行っスデューサ装備といったハードウェアの更新によている．観測点は，主に日本海溝および南海トラる観測精度・効率の向上（Sato et al., 2013）だけ

フといった日本近海の海溝・トラフ沿いの陸側斜でなく，複数エポック一括解析による重心推定法面に設置されており，測量船を用いて繰り返し観（松本・他，2008）に代表されるように，解析ス測を行っている．これまでに，宮城県沖，福島県キームを改良することによる精度向上が図られて沖等における地震間の海底地殻変動（藤田，2006

きた．本稿では，重心推定法が解析結果に与える† Received September 19, 2014; Accepted November 10, 2014
＊ 1 海洋調査課 海洋防災調査室 Geodesy and Geophysics Office, Hydrographic Surveys Division＊ 2 技術・国際課 海洋研究室
Ocean Research Laboratory, Technology Planning and International Affairs Division− 88 −

Study of stability and efficiency of the transponder array constraint method for the GPS acoustic observation

み合わせて，海中音速度構造とトランスポンダの効果について考察を行った．

位置を最小二乗法によって推定している（藤田・2海底地殻変動観測の手法他，2004）．

海底地殻変動観測システムの概念図を Fig. 1

各観測点の基準位置（以後，局位置と呼ぶ）は，に 示 す． 本 観 測 で は， マ ス ト に GNSS（Global

4 台のトランスポンダで構成されるアレイの重心Navigation Satellite System）アンテナを，船底に位置として定義される．この局位置を繰り返し測音響トランスデューサを装備した測量船を用い定し，経時変化を求めることで，海底の動きを検て，海底に設置された音響トランスポンダ（以出する．

後，トランスポンダと呼ぶ）の位置をセンチメー

局位置を求める際には，原理的には，個々のトトルの精度で求めている．トランスポンダは，特ランスポンダの 3 次元座標値をそれぞれ独立に求定の信号に応答して受信した音響波形をそのままめたうえで重心位置を算出（この手法を，以後，返信する機器である．各観測点には，その海域の独立解析と呼ぶ）できる．一方で，精度を向上さ水深と同程度の直径を持つ円周上の東西南北にせるために，推定するパラメータの数を減らす解計 4 台のトランスポンダが設置されている．これ析手法も用いられている．松本・他（2008）は，らのトランスポンダに対し，測量船のトランス局所的な地殻変動が無視できる場合に，トランスデューサを用いて音響測距観測を実施している．

ポンダの相対位置関係（以後，アレイ形状と呼ぶ）

観測は，水深の 2 倍程度の直径を持つ広さの海上を既知として拘束することにより精度向上が実現に設定された測線に沿って航走しながら実施されできるという重心推定法のアイデア（藤田・他，る．測距データは，各観測点につき一回の観測2005）を解析において適用するため，全観測エポッエポック（位置測定を実施する一回の観測単位）クに亘ってアレイ形状が不変であるという条件下当たり約 5000 ショット取得される．さらに，海で各エポックの局位置の相対位置変化を求める，中の音速度構造を把握するため，数時間ごとに複数エポック一括局位置解析を導入した（以下，CTD（Conductivity Temperature Depth profiler）
・

本稿ではアレイを拘束する手法であることを明示

XCTD（eXpendable CTD）・XBT（Expendableするため，重心推定法をアレイ固定手法と呼ぶこBathythermograph）による水温・塩分濃度プロととする）．一括解析では，アレイ形状を複数のファイルの観測を行う．これらのデータと GNSSエポックのデータを用いて決定するため，観測エによって求めた測量船のグローバルな位置とを組ポック数が増加することで，より安定した解が得られると期待される．
しかし，一括解析によるアレイ固定手法は，新たなエポックが加わるたびに，これまでに蓄積された全てのエポックのデータも含めて解析する必要があり，計算コストが大きくなる．一方で，アレイ形状はデータの増加に伴ってある一定の形状に収束すると考えられるので，アレイ形状の不定性によって生じる局位置解の不定性が，局位置解の精度の範囲内で無視できるようになるエポック

Fig. 1．Schematic picture of the GPS acoustic seafloor
geodetic observation.
図 1．GPS−音響測距結合方式による海底地殻変動観測の模式図．

数が存在すると考えられる．もしそのエポック数がわかれば，それ以降のエポックに対して，それ以前のエポックの一括解析によって得られたアレイ形状を拘束条件として解析することができ，計− 89 −

Yuhei AKIYAMA and Shun-ichi WATANABE
Table 1．Data summary used in this study.
表 1．解析に使用したデータ．

算コストを抑えることができる．石川・佐藤（2012）

ポンダを設置し，更新をしている．その際，継続は，アレイ固定手法の精度評価のためアレイ形状的な観測を行うため，新・旧トランスポンダに対の決定精度に関する考察を行い，一括解析をするする同時観測を実施することで，それらの相対的エポック数の増加に伴い，アレイ形状が安定するな位置関係を結び付けたうえで，新トランスポン傾向を見出した．この考察をさらに進めることダへの観測に移行する（石川，2011）．この更新で，アレイ形状が局位置解の精度に対して十分に作業により異なるアレイが構成されるため，更新安定するために必要なエポック数について議論で前を旧局，更新後を新局として，それぞれ別のデーきる．

タセットとして扱った．

そこで本稿では，まず石川・佐藤（2012）と同これらの条件を満たし，かつエポック数が 5 回様の手法で，エポック数の増加に伴うアレイ形状以上となるデータセット「宮城沖 1（新局）」，「宮の安定化について調べた．その上で，アレイ形状

城沖 1（旧局）」，「釜石沖 1」，「釜石沖 2」および「東の変化が局位置解に与える影響について考察を海沖 1」について解析を実施した．各データセッ行った．

トの詳細については，Table 1 にまとめる．

3

3.2

アレイ固定手法（重心推定法）の評価手法

3.1

使用データ

解析手法

まず，一括解析に使用するエポック数の増加に

本稿の目的はエポック数に対するアレイ固定手伴うアレイ形状の安定化について調べるため，石法の安定性を見積もることであるため，エポック川・佐藤（2012）と同様に，使用エポック数に対数が十分に確保できるデータセットについて解析するアレイ形状の変化について調べた．
ここで，アレイ形状を定量的に評価するためのする必要がある．
ただし，アレイ固定手法が適用できるのは，ア指標を導入する．アレイの北，東，南，西のトラレイ形状が変化していないと期待される一連のエンスポンダをそれぞれ 1，2，3，4 とする．デー

ポックに対してである．そこで，東北地方太平洋タ期間の最初のエポックから第 M エポックまで沖地震に伴う地殻変動によってトランスポンダののデータを使用して解析した際の，第 i トランス相対位置が有意に変化した日本海溝沿いの観測点ポンダの位置を，Gi,j（M）とする．j ＝ 1，2，3については，本震以降のエポックのみを使用した．

は，それぞれ各海域における東西方向（東が正），また，トランスポンダは筐体に内蔵された電池に南北方向（北が正），上下方向（上が正）を示し，より動作しているため，定期的に新たなトランス座標原点はアレイ重心とする．このアレイ形状 G− 90 −Study of stability and efficiency of the transponder array constraint method for the GPS acoustic observation

を，基準形状と比較して評価する．その基準形状

1（新局）」ではσは 2 cm 程度以内であった．「宮には，各データセットにおける最大のエポック数

城沖 1（旧局）」では M ＝ 1 のときにσが大きく

Mmax でのアレイ形状，G（Mmax）を用いた．これなっているが，これは，最初のエポックで推定さは，エポック数が多くなるほどアレイ形状を決定れたアレイ形状が適切でなく，その後のエポックするために使用できるデータ量が増加し，最も確の追加によって修正されたためと考えられる．
「釜からしい解が得られるためである．

石沖 2」では M ≧ 4 でσは 3 cm 程度以内に収束

さらに，アレイ形状の変形量を成分ごとに比較した．これらに対し，「釜石沖 1」はσが 4−7 cm と，Gi,j（M）の Gi,j（Mmax）に対する偏差の，するため，他の観測点に比べて大きかった．しかし，「釜石i ＝ 1，2，3，4 についての RMS を求め，これを沖 1」も含め，すべての観測点において，エポッσ j（M）と定義する．このσ j を，M についてク数の増加に伴ってアレイ形状が安定する傾向が評価した．

見られた．また，「釜石沖 2」の第 3 エポックの次に，アレイ形状が局位置解に与える影響につように，アレイ形状の収束する傾向から大きく外いて調べた．まず，それぞれの M ＜ Mmax についれるケースもあったが，これは別のエポックがさて，G（M）の配置でアレイ形状を固定し，M ＋らに加わっていくことで解消された．

1 番目以降のエポックのデータについて局位置を次に，各 M について得られた ΔPh（M）およ求めた．これは，第 M エポックまでのデータでび ΔPv（M）を，Fig. 3 に示す．「宮城沖 1（新局）」，アレイ形状が十分収束したと仮定し，それ以降の「宮城沖 1（旧局）」および「釜石沖 1」では，Mエポックについて G（M）を固定して解析するこが増加するにつれて ΔPh の値が減少する傾向が見とに等しい．得られた第 n エポックの局位置をられた．一方，「東海沖 1」では M の増加に伴うPj（n，G（M））とする．なお，各アレイの絶対ΔPh の値の系統的な変化は見られなかったが，全位置のゆらぎによる局位置時系列のオフセットを体を通して ΔPh が 2 cm 以下であった．「釜石沖無視するため，P は，各時系列における第 1 エポッ2」でも系統的な変化は見られず，ΔPh が 4 cm 程クからの相対変位として定義する．これは，海底度であった．ΔPv については，どの観測点につい地殻変動観測の目的が，相対変位および変位速度ても，M の増加に伴う系統的な変化は見られなの抽出であるためである．

かった．また，σが特に大きい「宮城沖 1（旧局）」

これらの局位置解について，全エポック一括解析による結果 Pj（n，G（Mmax））を基準解としてM ＝ 1 については，ΔPh および ΔPv の値も大きくなっていた．

比較する．各 M ＜ Mmax について，局位置の水平成分と上下成分の基準解に対するばらつきをそれ5

考察

ぞれ定量的に評価するため，評価パラメータ ΔPh

エポック数の増加に伴うアレイ形状の変化につP（G（M））（M），ΔP（を導入する．それぞれ，v M）
j n，いて調べた結果（Fig. 2），石川・佐藤（2012）の−Pj（n，G（Mmax））の水平成分および上下成分結果と同様，一括解析を実施するエポック数の増の絶対値の，M ＋ 1 ≦ n ≦ Mmax についての平均加に伴い，アレイ形状がある一定の形状に収束す値として求めた．

ることが確認された．それに加えて，観測海域による違いも顕著であることがわかった．また，「釜4．結果石沖 2」の結果からは，あるエポックの追加によっ

Fig. 2 に，エポック数の増加に伴うアレイ形てアレイ形状が収束の傾向から外れることも生じ状 G（M）の変化の様子と，G（M）の G（Mmax）うるが，さらなるエポックの追加で解消されるこからの偏差の RMS，σを示す．「東海沖 1」ではとが確認された．このことは，エポックごとの異いずれの M に対してもσは 1 cm 程度，「宮城沖常値に対するアレイ固定手法の有効性を示してい− 91 −

Yuhei AKIYAMA and Shun-ichi WATANABE

宮城沖１(新局)
1510EWNS

98UD

756σ(cm)Northward (cm)

101 epoch
2 epochs
3 epochs
4 epochs
5 epochs
6 epochs
7 epochs
8 epochs
9 epochs
10 epochs

054-532-10

10-15-15-10-505Eastward (cm)101152345678
M : Number of epoch used in multi-epoch analysis

910 (Mmax)

宮城沖１(旧局)
1510EWNS

98UD

756σ(cm)Northward (cm)

101 epoch
2 epochs
3 epochs
4 epochs
5 epochs

054-532-10

10-15-15-10-505
Eastward (cm)101155 (Mmax)234
M : Number of epoch used in multi-epoch analysis

釜石沖１1510EW
NS

98UD

756σ(cm)Northward (cm)

101 epoch
2 epochs
3 epochs
4 epochs
5 epochs
6 epochs
7 epochs

054-532-10

1-15

0-15-10-50
5Eastward (cm)101512345
M : Number of epoch used in multi-epoch analysis

7 (Mmax)

66σ(cm)Northward (cm)

釜石沖２Fig. 2．
Estimated transponder positions relative to the full epoch position（left）. To emphasize their shapes, the1015
1 epoch plotted to form a diamond whose diagonals are 20 cm. RMS of the relativereference transponders are
2 epochs
3 epochs
EW
94 epochs
positions of the transponders（right）.
5 epochs to the full epoch positions
6 epochs
10NS7 epochs
88 epochs
図 2．エポック数の増加に伴うアレイ形状（左）．形状の変化を強調するため，各トランスポンダを20 cm の長さUD7．
5の対角線を持つ正方形からのずれとして表した．各 G（M）における G（Mmax）からの偏差の RMS（右）
054-5

32-10

1− 92 −0-15-15-10-50
5Eastward (cm)1015

123456
M : Number of epoch used in multi-epoch analysis

78 (Mmax)

-532-101

Studyof stability and efficiency of the transponder
array constraint method for the GPS acoustic observation-150-15-10-505Eastward (cm)101512345M : Number of epoch used in multi-epoch analysis

7 (Mmax)

6釜石沖２15101 epoch
2 epochs
3 epochs
4 epochs
5 epochs
6 epochs
7 epochs
8 epochs

108UD

75Northward (cm)EW
NS

9σ(cm)

60

54-532-10

10-15-15-10-50
5Eastward (cm)1011523456
M : Number of epoch used in multi-epoch analysis

8 (Mmax)

7東海沖１15105EW
NS

98UD

76σ(cm)Northward (cm)

101 epoch
2 epochs
3 epochs
4 epochs
5 epochs
6 epochs
7 epochs
8 epochs
9 epochs
10 epochs
11 epochs
12 epochs
13 epochs

054-5

32-10

1-150-15-10-50
5Eastward (cm)

1015

1

2

3

4
56
78
910
M : Number of epoch used in multi-epoch analysis

111213 (Mmax)

Fig. 2．（continued）
図 2．（続き）(a) horizontal(b) vertical

1010

9

98

8

7

7

6

6

ΔPv (cm)

ΔPh (cm)

34.7cm(M=1)

5435.1cm(M=1)

宮城沖１（新局）
宮城沖１（旧局）
釜石沖１
釜石沖２
東海沖１543322110012345678910111213

M : Number of epoch used in mul-epoch analysis

12345678Fig. 3．Values of（a）ΔPh and（b）ΔPv.
図 3．M ＝ Mmax を基準とした局位置解の偏差の（a）水平成分 ΔPh，と（b）上下成分 ΔPv．

− 93 −

910

M : Number of epoch used in mul-epoch analysis

111213

Yuhei AKIYAMA and Shun-ichi WATANABE

直上の海面を中心とする測線に沿って実施される．
次に，アレイ形状が局位置に与える影響についる．つまり，平均的にはアレイの中心から各トラて調査した結果（Fig. 3）からは，「釜石沖 2」をンスポンダに対して測距される．解析の際に同時除き，一括解析に使用するエポック数 M が増加に推定する音速度の誤差は，それぞれの測距におすることで，局位置の水平成分が 1−2 cm 程度いて視線方向の距離変化となって表れるので，例に収束する傾向が見られた．「釜石沖 2」ではア

えば，音速度の推定値が時空間的に一律に，実際

レイ形状 G（M）は M の増加に伴って収束する

より大きく見積もられた場合には，見かけの距離傾向にあるが，ΔPh の値が大きい．この原因につが遠くなる．各トランスポンダに対する測距は，いては，本稿の解析結果が M ＝ Mmax における解アレイ中心側からのデータが多いため，この場合，を基準としているため，第 Mmax エポックの局位測線の中心からより遠いところにトランスポンダ置解に含まれる誤差を捉えたものである可能性が位置が推定される．すなわち，アレイが相似形にある．また，「宮城沖 1（旧局）」では M ＝ 1 の拡大し，局位置がより深く推定されると考えられときに ΔPh および ΔPv の値が大きくなったが，こる．

れは第 1 エポックのデータから求められたアレイ

この仮説について定量的な評価を行うため，一形状が，何らかの理由のより，他のすべてのエポッ括解析で得られたアレイ形状 G（Mmax）と局位置クから得られる形状と異なっているためと考えら

上下成分 hmulti（n）＝ P3（n，G（Mmax））を基準

れる．次のエポックが加わることで ΔPh およびとして，独立解析で各トランスポンダについて位ΔPv の値が収束する傾向がみられることから，こ置を推定した際のアレイ形状 Geach（n）と，局位のデータは本議論には影響がないといえる．

置上下成分 heach（n）と比較をした．アレイ形状

その他の観測点については，収束の傾向がそれの拡大・縮小を表現するパラメータとしては，アぞれ異なっている．特に「東海沖 1」では，エポッ

レイ G を水平面上に投影した際の面積 S（G）をク数に関わらず ΔPh が 1−2 cm 以内で安定して使用した．

おり，「宮城沖 1（新局）」でも，5 エポック目以降，その結果として，全観測点のそれぞれのエポッ2 cm 以内となる．Fig. 2 の結果を見ると，アレイクにおけるアレイ水平投影面積の変化率 SR ＝ S

形状の変形量の指標であるσも，「東海沖 1」，お（Geach）/S（G（Mmax））と，局位置の上下成分のよび「宮城沖 1（新局）」の 5 エポック目以降は，2変化率の 2 乗 hR2 ＝（heach/hmulti）とを Fig. 4 に示す．

各成分 1 cm 程度以内で概ね安定している．このFig. 4 から，すべての観測点について，アレイ面ことにより，G（M）が 1 cm 程度の精度で安定積と局位置の上下成分が同一の比率で変化するこして得られれば，G（M）の不定性による局位置とが確認された．

解のばらつきは 1−2 cm 程度に抑えられると言もしアレイ面積と局位置上下成分が前述のようえる．しかしながら，得られたデータ数は十分とな音線の幾何学的配置で決定されるのであれば，は言えないので，今後もデータ蓄積による継続的測線とトランスポンダの幾何学的配置は観測点にな評価が望まれる．

よらずほぼ相似形をなしているため，観測点によ最後に，アレイ形状と局位置の上下成分との関らず SR と hR2 は同等の関係性を持つことが期待係について考察を行う．石川・佐藤（2012）も指される．Fig. 4 で得られた，両者の比率が観測点摘しているように，音速度推定の系統的誤差は，にほとんど依存しないという結果は，この仮説のアレイ形状を相似形に保ったまま拡大・縮小する帰結と調和的である．このことは，アレイ形状の変形と，局位置の上下成分とに同時に影響を与え固定が水中音速度の一律なバイアス誤差の抑制，る．これは，音響測距時の平均的な視線方向と関ひいては局位置上下成分の補正に対して重要な役連付けて説明される．音響測距は，アレイの重心割を果たしていることを示唆している．

− 94 −Study of stability and efficiency of the transponder array constraint method for the GPS acoustic observation
1.0003

えられる．加えて，トランスポンダの更新時の新宮城沖１（新局）

1.0002

宮城沖１（旧局）

1.0001

釜石沖２

旧局同時観測における適切なアレイ固定手法について，例えば既に多くのデータが得られている旧

釜石沖１

局アレイの形状のみを固定して，データの少ない

東海沖１新局アレイの形状を求めるといった手法について

SR

1も考察し，解析効率の向上を目指していきたい．

※0.9999

60.9998

本稿では，トランスポンダアレイ形状の安定性とアレイ固定手法の効果を見積もるため，一括解0.99970.99960.99975まとめ析の使用エポック数がアレイ形状に与える影響，※ 2012/04 のエポック

0.999850.99995

hR21

およびアレイ形状が局位置解に与える影響につい1.0001て調べた．その結果，石川・佐藤（2012）の結果

2Fig. 4．SR plotted against hR .
図 4．各エポックにおけるアレイ水平投影面積の変化
率 SR と局位置の上下成分の変化率の 2 乗 hR2．と同様，使用エポック数の増加に伴いアレイ形状逆に，これらの関係性から外れているエポックの安定性には海域による依存性があることが示唆については，さらに他の要素の影響を受けていされ，解析時間を短縮する際には，海域ごとに必ると考えられる．例えば，「釜石沖 1」の 2012 年要なエポック数を指定する必要があると考えられ2が安定することが確認され，一括解析による局位置推定の精度向上が確認された．また，一括解析

4 月のエポックについては，SR と hR の比が，他る．また，アレイ形状の水平投影面積と局位置ののエポックに比して大きく異なっている．一括解

上下成分には強い相関が見られ，アレイ固定手法析によって得られた局位置の時系列については海の音速度補正への有効性が確認された．

上保安庁（2014）によって報告されているが，こ今後，本稿で得られた結果を基に，数値計算なのエポックについては時系列の傾向から大きく外どを用いたアレイ固定手法の効果の評価による解れたものであった．実際に，このエポックでは析効率の向上や，局位置の水平成分に特に影響をXBT 等の海中温度観測によって音速度構造の空及ぼすアレイ形状パラメータを抽出することによ間的不均質が確認されており，それが局位置解のるエポックごとの誤差評価や誤差要因の推定を目誤差の原因である可能性がある．

指すことで，海底地殻変動観測の精度向上に寄与

今後，局位置解に特徴的な影響を与える，アレできるであろう．

イ面積 S（G）のようなアレイ形状パラメータを謝さらに抽出することで，各エポックの局位置解の辞

KGPS 解 析 に は，NASA/GSFC の Colombo 博

誤差の評価や，個々のエポックについての誤差要

士開発のソフトウェア「IT」（Colombo, 1998）を因の特定が可能になるかもしれない．
本稿では実際のデータを用いた解析結果を示し

用い，陸上基準点として，国土地理院より電子基たが，基準とした全エポック一括解析による解自準点 1 秒データを提供いただきました．また，藤体にも誤差が含まれうるという手法上の限界があ

田雅之博士，松本良浩氏，および渡邊奈保子氏にる．今後の課題として，引き続き蓄積される実デーは，本論文の改善に資する大変有益なコメントをタによる検証だけでなく，数値計算などによる，いただきました．記して感謝します．

より制御された条件下での解析などを通じて，アレイ固定の効果を見積もることが必要であると考− 95 −

Yuhei AKIYAMA and Shun-ichi WATANABE

文

献

Yokota, N. Ujihara, and A. Asada（2014）

Colombo, O. L（1998）Long Distance Kinematic

Evidence of viscoelastic deformation

GPS, GPS for Geodesy, 2nd Edition, 537 567,

following the 2011 Tohoku oki earthquake

Springer.

revealed from seafloor geodetic observation,Geophys. Res. Lett., 41, 5789 5796,

藤田雅之・佐藤まりこ・矢吹哲一朗（2004）海底
地殻変動観測における局位置解析ソフトウェ

doi:10.1002/2014GL061134.

アの開発，海洋情報部技報，22，50 56．
要

藤田雅之・石川直史・松本良浩・望月将志・佐藤旨

まりこ・矢吹哲一朗・浅田昭（2005）宮城県

海上保安庁海洋情報部で実施している GPS−

沖海底の地殻変動と重心推定法による時系列

音響測距結合方式による海底地殻変動観測の解析の検証，日本測地学会第 104 回講演会要旨，効率を向上させるため，トランスポンダアレイの47 48．

形状を拘束する重心推定法の安定性と効果につい藤田雅之（2006）GPS/ 音響測距結合方式によるて評価を行った．使用するエポック数やアレイ形

海底地殻変動観測〜海上保安庁の取り組み

状といった解析条件が局位置解に与える影響につ（レビュー）〜，海洋情報部研究報告，42，1

いて調べた結果，一括解析で使用するエポック数14．の増加に伴い，アレイ形状がある一定の形状に収石川直史（2011）海底地殻変動観測における海束し，局位置解の水平成分も一定値に収束するこ底音響基準局の更新，海洋情報部研究報告，とが確認された．これにより，アレイ固定手法の47，44 50．

有効性が示された．また，アレイ形状の水平投影

石川直史・佐藤まりこ（2012）海底地殻変動観測

面積と局位置の上下成分に強い相関があることがにおける重心推定法の評価，海洋情報部研究確認され，アレイ固定手法の音速度補正への有効報告，48，74 84．

性が示された．

海上保安庁（2014）東北地方太平洋沖地震後の海底地殻変動観測結果，地震予知連絡会会報，
92，89 93．
松本良浩・藤田雅之・石川直史（2008）海底地殻変動観測における複数エポック一括局位置解
析手法の導入，海洋情報部技報，26，16 22．
Sato, M., T. Ishikawa, N. Ujihara, S. Yoshida, M.
Fujita, M. Mochizuki, and A. Asada（2011）
Displacement above the hypocenter of the2011 Tohoku oki earthquake, Science, 332,
1395, doi:10.1126/science.1207401.
Sato, M., M. Fujita, Y. Matsumoto, H. Saito, T.
Ishikawa, and T. Asakura（2013）Improvementof GPS/acoustic seafloor positioning precision
through controlling the ship's track line, J.
Geod., 87, 825 842, doi:10.1007/s00190 013
0649 9.
Watanabe, S., M. Sato, M. Fujita, T. Ishikawa, Y.
− 96 −