「LAGEOS」SLR データ解析による下里の移動速度の再計算

海洋情報部研究報告 第 45 号 平成 21 年 3 月 28 日
REPORT OF HYDROGRAPHIC AND OCEANOGRAPHIC RESEARCHES No.45 March, 2009
「LAGEOS」SLR データ解析による下里の移動速度の再計算†
佐藤まりこ＊１，松下 優＊２，藤田雅之＊３，仙石 新＊４
Re!estimation of the velocity of Simosato from analysis of LAGEOS SLR data†
Mariko SATO＊１, Hiroshi MATSUSHITA＊２, Masayuki FUJITA＊３ and Arata SENGOKU＊４
Abstract
Simosato Hydrographic Observatory in Wakayama Prefecture has been carrying out Satellite Laser Ranging
（SLR）observation for AJISAI, LAGEOS‐１, LAGEOS‐２, and other geodetic and Earth!observation satellites,
since１９
８２. More than２
６，
００
０passes have been obtained.
In this report, we re!estimated the velocity of Simosato from LAGEOS SLR data for１５years. As a result, the
velocity of２．
９cm/year in direction to２９
４°N relative to the Eurasian plate was obtained. This velocity is close to
the subduction velocity of the Philippine Sea plate. This indicates that the interplate coupling at this region is
strong. In addition, our result is consistent with other researches using SLR data
（e.g., ITRF２０
００，ITRF２
００
５SLR
rescaled）
.
１ はじめに
れた離島や本土の海洋測地網一次基準点等１
５ヶ
海洋測地網の本土基準点である下里水路観測所
所における可搬式 SLR 装置によるキャンペーン
（和歌山県那智勝浦町）では，１９
８２年以来，国産
観測では，島嶼等の精密な位置決定のための本土
の測地衛星「AJISAI」や米国の測地衛星「LAGEOS
基準点としての役割を果たした（例えば，Sen-
‐１」「LAGEOS‐
２」等 の 人 工 衛 星 レ ー ザ ー 測 距
goku，１
９９
９）．
（SLR : Satellite Laser Ranging）観測を行ってお
長期にわたる定常的な SLR 観測は，長波長の
り，現在までの観測回数は２
６，
０
００パスを超えて
地球重力場とその時間変動の決定に貢献し，地球
いる（Fig．
１）
．
科学の基礎となる国際地球基準座標系（ITRF : In-
これまでに，２００
２年４月に経緯度の基準を世
ternational Terrestrial Reference Frame）の構築
界測地系への移行するため，日本測地系から世界
に寄与しているとともに，グローバルな視点でプ
測 地 系 へ の 変 換 パ ラ メ ー タ を 決 定（仙 石・
レート運動を観測することを可能とし，地震の発
他，２００
０）して「海洋測地成果２００
０」として取
生メカニズムの解明等の地球科学の理解にも貢献
りまとめた．また，１９８
８年から２００
１年まで行わ
している．
†Received December２
６，２
０
０
８；Accepted January２
９，２
０
０
９
＊１ 航法測地室 Geodesy and Geophysics Office
＊２ 第二管区海上保安本部２nd Regional Coast Guard Headquarters
＊３ 技術・国際課 Technology Planning & International Affairs Division
＊４ 海洋調査課 Hydrographic Surveys Division
−1−
Mariko SATO, Hiroshi MATSUSHITA, Masayuki FUJITA and Arata SENGOKU
行ったので，その結果について報告する．
２ SLR データ解析
まずはじめに，SLR データ解析について述べ
る．
解析に用いたデータは，１９９
３年から２
００
７年ま
でに世界中の SLR 観測局で取得された LAGEOS‐
１及 び LAGEOS‐２の ノ ー マ ル ポ イ ン ト（NP）
データである．年ごとの観測局数，パス数，NP
データ数を Table．
１に示す．観測局数は年平均
Fig．
１ Satellite Laser Ranging at the Simosato Hydrographic Observatory.
３７局であり，２
００５年以降若干減少しているもの
の，パス数及び NP データ数は増加傾向にあり，
現在でも世界中で精力的に観測されていることが
SLR デ ー タ か ら プ レ ー ト 運 動 を 求 め る 研 究
分かる．
は，１９
７０年代末以来，数多く行われており（例
SLR データの解析は，NASA ゴダード宇宙飛行
えば，Smith et al., １
９７
９，Tapley et al., １９
８５な
セ ン タ ー（GSFC）で 開 発 さ れ た GEODYN−II/
ど），その大半は世界中の SLR 観測局で日々多く
SOLVE（Eddy et al.,１９
９０；Ullman,１
９９
２）を用い
のデータが得られている米国の測地衛星
た多衛星結合解析により行い，各年における各観
「LAGEOS」の SLR デ ー タ を 使 用 し た も の で あ
測 局 の 位 置（１年 解）を 推 定 し た（藤 田・仙
石，１９
９７）．GEODYN−II の解析に用いたモデル
る．
Sengoku（１
９９
８）は，１９
８６年 か ら１９
９４年 ま で
は概ね次のとおりである．詳細については，藤
の８年間の「AJISAI」SLR データを用いて，下里
田・他（１９
９８），松 下・他（２
０
０５a）を 参 照 さ れ
を含む SLR 観測局のプレート運動を求め，下里
たい．
のユーラシアプレート安定域に対する移動速度
・重力モデル：GGM０１C（Tapley et al.,２
００
３）
（２
９１°，３．
２cm/year）を報告した（以下，本稿で
・海洋潮汐荷重変形：IERS Conventions（McCarthy,１９
９６）
は，方位角を北から時計回りの角度として定義す
る）．
「AJISAI」は，「LAGEOS」と 比 べ て 衛 星 の
軌道高度が低く，また形状も大きく表面力の影響
を大きく受けるため，衛星に働く力のモデル化が
Table．
１ Summary of global LAGEOS data from１
９
９
３to
２
０
０
７．
難しく，地球基準座標系の決定には不利である
が，「LAGEOS」SLR データ等から求められるプ
レート運動と整合的な結果が得られることを示し
た．この成果は，当庁が行っている GPS 地殻変
動監視観測や海底地殻変動観測においても本土基
準点の動きとして，現在でも幅広く活用されてい
る．
今般，Sengoku（１
９９８）による移動速度の推定
か ら１
０年 が 経 過 し た こ と か ら，近 年 の
「LAGEOS」SLR データを用いて下里のユーラシ
アプレート安定域に対する移動速度の再計算を
−2−
Re!estimation of the velocity of Simosato from analysis of LAGEOS SLR data
Fig．
２ Time series of estimated coordinate of Simosato station. Large earthquakes（M７．
４，M６．
９）occurred about１
０
０km southeast of Simosato on
September５th，２
０
０
４.
・地球回転パラメータ：IERS Bulletine B の最終
得されていることから，座標を固定する局として
妥当と判断している．
値
・基準座標系：ITRF２０
００（Boucher et al.,２０
０４）
解析結果の一例として，下里の経緯度及び楕円
・解析で固定した座標及び速度場：ITRF２
００
０で
体高の推移を Fig．
２に示す．参考として１ヶ月解
与 え ら れ た 米 国 メ リ ー ラ ン ド 州 の Greenbelt
（１ヶ月ごとの座標値）も併せて示した．Fig．
２よ
（７１０
５）局 の 経 緯 度 と そ の 速 度，英 国 の
り，１ヶ月解にはばらつきが見られるが，１年解
Herstmonceux
（７８４
０）局の緯度とその速度．
では非常に安定していることが分かる．２００
４年９
Greenbelt 局 と Herstmonceux 局 は，プ レ ー ト
月５日 に は，下 里 の 南 東 約１
００km の 海 域 で M
境界から十分に離れておりプレート内の安定域に
７．
４及び M６．
９の地震（２
００
４年紀伊半島南東沖
あることからほぼ一様な動きをしているとみなす
地震）が発生し，この地震に伴う地殻変動として
ことができ，また解析期間中十分なデータ量が取
南 に２cm，１∼２cm の 沈 降 が 報 告 さ れ て い る
−3−
Mariko SATO, Hiroshi MATSUSHITA, Masayuki FUJITA and Arata SENGOKU
えられるため，本稿では各１年解に対する地心補
（松下他，２００
５）．
正は行っていない．
３
移動速度の算出
具体的には，以下の手順で移動速度を算出し
３．１ 算出方法
た．
移動速度の算出は，観測局間のアーク長の変化
!各観測局の観測期間，データ量，１年解の安定
率を使用する方法で行った．アーク長とは，地球
性等から，移動速度の算出に使用する観測局を
を半径が赤道半径に等しい球と仮定したときの観
選定する．
測局間の球面上の距離であり，次の式で表わされ
"年ごとに２局間のアーク長を全ての組み合わせ
る（Fig．
３）．
について求め，直線回帰により，それぞれの
l = aeΨg
アーク長の変化率を計算する．
＝aecos （sinφ１sinφ２＋cosφ１cosφ２cos（λ１−λ２））
−１
#アーク長の変化率から，最小二乗法により，各
ここで，ae は地球の赤道半径，Ψg は観測局間
観測局の絶対速度（ここでは，no−net−rotation
の地心角度，φ は観測局の地心緯度，λ は観測局
（NNR）系に基づく速度を便宜的にこう呼ぶ）
を推定する．
の地心経度である．通常の SLR 解析で得られる
座 標 の 時 系 列 に は，座 標 系 の 原 点（地 心），方
$#で得られた絶対速度から，地質学的に求めら
向，スケールの不一致による誤差が含まれている
れたプレート運動モデル（NNR NUVEL‐
１A :
が，アーク長は座標系の方向とスケールの誤差に
DeMets et al., １９
９４）で計算される速度を差し
よる影響を受けないため，より安定した移動速度
引く．
を求めることができる．さらに，アーク長の変化
なお，#の速度の推定においては，プレート内
率は，観測局の水平速度のみに依存し，レンジバ
の安定域に位置し，NNR NUVEL‐１A で表わされ
イアス等により生じる観測局位置誤差の鉛直成分
るプレート運動とほぼ一致した動きをしている米
の影響を受けないという利点があり，プレート運
国メリーランド州の Greenbelt
（７１
０５）
局の速度，
動のような水平速度を求めるのに適している
英国の Herstmonceux
（７８４
０）局の南北方向の速度
（Sengoku，１９９
８）．なお，第２節で得られた１年
を，同モデルに整合的な ITRF２
０００（Boucher et
ごとの座標系の地心のばらつきはたかだか１cm
al., ２０
０４）で与えられている速度に固定した．こ
程度であり，アーク長に与える影響は小さいと考
れにより，ここで求めた速度は，定義上，no-netrotation
（NNR）系に基づく速度と考えることがで
きる．その量的な整合性の評価については３．
４で
議論する．
また，$で得られる速度は，言い換えると，各
観測局の位置するプレートが仮想的に剛体だとし
たときの速度と観測された速度との差であり，
ユーラシアプレート上にある下里では，この速度
が本稿の目的であるユーラシアプレート安定域に
対する移動速度となる．
３．
２ 観測局及び観測期間の選定
第２節で得られた１
９
９３年から２
００７年までの各
Fig．
３ Definition of arc between station１and station２
（Sengoku，１
９
９
８）
.
観測局の１年解の時系列から，５年以上安定した
１年解が得られていることや観測局の幾何学的配
−4−
Re!estimation of the velocity of Simosato from analysis of LAGEOS SLR data
Table．
２ List of SLR stations used for the estimation of velocities.
置等を考慮して，移動速度の算出に使用する観測
局２
４局 を 選 定 し た（Table．
２及 び Fig．
４）．ま
た，観測が行われていてもデータが少ない等の理
由により，１年解のばらつきが大きい期間は除外
しており，Table．
２には移動速度の算出に使用し
た期間も示した．
先に述べたように，下里（７８
３８）
局の近傍では，
２０
０４年９月５日に紀伊半島南東沖地震（M７．
４，
M６．
９）が発生し，２cm 程度の地殻変動が検出さ
れた（松下他，２０
０５）
．２００
４年の１年解について
は，１年分のデータに占める地震後のデータが少
なく，２
００４年の１年解への地震の影響は小さい
と考えられることから，移動速度の算出に使用し
た．２０
０５年以降の１年解には，この地震に伴う
Fig．
４ Locations of SLR stations used for velocity estimation.
地殻変動に加え，測距値の較正方法の一時的な変
更や信号の検出方法の変更等による人為的な影響
が反映されている可能性が高いため，移動速度の
算出には使用していない．また，Arequipa
（７４
０３）
局 は，２００
１年６月２
３日 に ペ ル ー 南 部 地 震（M
−5−
Mariko SATO, Hiroshi MATSUSHITA, Masayuki FUJITA and Arata SENGOKU
８．
１）が発生したため，１年解が安定している１９
９５
一例として，下里と各観測局とのアーク長変化
年から２０
００年までの解を使用した．ただし，Har-
率を Table．
３に示す．アーク長変化率の推定誤差
tebeesthoek
（７５
０１）局は２００
０年に観測を開始して
は平均０．
０
８cm/year で，全アークの平均より小
い る た め，Arequipa
（７４０
３）局 と Hartebeesthoek
さい．McDonald
（７０
８０）
局 と Arequipa
（７
４０
３）局 は
（７５０
１）局 と の ア ー ク 長 変 化 率 を 求 め る と き の
下里とのアーク長が若干伸びているが，それ以外
み，２０
０２年 か ら２
０
０７年 ま で の１年 解 を 使 用 し
の観測局は２∼７cm/year の速度で縮んでいる．
た．
下里との変化率が最も大きいのはオーストラリア
移動速度の算出に使用した観測局はヨーロッパ
の Yarragadee（７０
９０）局 で−７．３ cm/year，最 も 小
には多く集中しているが，南米，南アフリカ，
さいのは McDonald
（７０
８０）局で０．
３cm/year であ
オーストラリアの観測局を採用し，できるだけ異
る．
なるプレート上の観測局を使うよう努めた．
３．
４ 絶対速度の推定
３．３
アーク長変化率
３．
３のアーク長の変化率から最小二乗法により
次 に，３．
２で 選 定 し た２４局 の 観 測 局 に つ い
各観測局の絶対速度（no−net−rotation（NNR）系
て，２局間のアーク長を年ごとに全ての組み合わ
の速度）を求めた．速度の推定にあたっては，米
せで求め，直線回帰により，各アーク長の変化率
国メリーランド州の Greenbelt
（７１０
５）局の速度及
を計算した．アーク長の組み合わせの数は２７
６
び英国の Herstmonceux
（７
８４
０）局の南北方向の速
で，アーク長変化率の推定誤差は平均０．
１
１cm/
度を ITRF２
００
０で与えられている速度に固定し
year である．
た．絶対速度の推定結果を Table．
４及び Fig．
５に
Table．
３ Arc length rates between Simosato and other SLR stations.
−6−
* Vn, Ve of Greenbelt and Vn of Herstmonceux are fixed to the velocities of ITRF2000.
Table．
４ Estimated velocities and the differences from velocities of ITRF２
０
０
０，ITRF２
０
０
５
（SLR rescaled）
and NUVEL‐
１A.（unit : cm/year）
Re!estimation of the velocity of Simosato from analysis of LAGEOS SLR data
−7−
Mariko SATO, Hiroshi MATSUSHITA, Masayuki FUJITA and Arata SENGOKU
Fig．
５ Velocity vectors of SLR station：（a）World，
（b）Europe，
（c）Asia. Red arrows, blue arrows,
purple arrows and green arrows show the estimated velocities, velocities of ITRF２
０
０
０，ITRF
２
０
０
５
（SLR rescaled）
and NUVEL‐
１A, respectively.
示す．なお，Table．
４及び Fig．
５には，比較のた
ての観測局について０．
１cm/year 以内である．
め，ITRF２０
００，ITRF２
００
５
（SLR rescaled, Altamimi
この結果を ITRF２
００
０及び ITRF２０
０５
（SLR res-
et al., ２０
０７）
，NUVEL‐１A による速度及びそれぞ
caled）
による速度と比較すると，一部の観測局を
れの速度と今回推定された速度との差も併せて示
除き，差は０．
１∼０．
２cm/year と非常に良い一致
している．
を示していることが分かる．Riyadh
（７８
３２）局で
絶対速度の推定誤差は，南北方向，東西方向と
は，ITRF２
００
０による速度との差が２cm/year 以
も に 平 均０．
０
３cm/year で あ り，南 半 球 に あ る
上 あ る が，ITRF２０
０５
（SLR rescaled）
で速度が改
オーストラリアの Mt.Stromlo
（７
８４
９）局とペルー
善され，我々の結果との差はほとんどなくなっ
の Arequipa
（７４
０３）局ではやや大きいものの，全
た．一方，Arequipa
（７
４０
３）局では，ITRF２０
００，ITRF
−8−
Re!estimation of the velocity of Simosato from analysis of LAGEOS SLR data
Fig．
６ Correlations between estimated velocities and（a）ITRF２
０
０
０/（b）ITRF２
０
０
５
（SLR rescaled）
.
Fig．
７ Differences between estimated velocities and velocities calculated from NUVEL‐
１A.
２０
０５
（SLR rescaled）
ともに推定された速度との差
Fig．
６に示す．いずれも相関係数は０．
９
５以上と
が 大 き い．こ れ は，本 稿 で 用 い た Arequipa
非常に強い相関を示しており，相対的には ITRF
（７４０
３）局の１年解の期間が５年と短いことや他の
２０
０５
（SLR rescaled）
の方が推定結果と合ってい
観測局に比べて１年解のばらつきが大きいことが
る．
影響している可能性があるが，今後慎重に検討し
次に，推定された絶対速度を地質学的に求めら
ていきたい．また，ITRF２０
００，ITRF２
００
５
（SLR res-
れたプレート運動モデルである NUVEL‐
１A と比
caled）
の速度と今回推定された速度との相関を
較する．絶対速度と NUVEL‐１A による速度との
−9−
Mariko SATO, Hiroshi MATSUSHITA, Masayuki FUJITA and Arata SENGOKU
差 を Fig．
７に 示 す．Table．
４及 び Fig．
７よ り，プ
レート内部に位置しているヨーロッパや南アフリ
カ，オーストラリア等の観測局では，NUVEL‐
１
A との差が０．
２cm/year 程度と小さく，SLR デー
タによる絶対速度の推定結果は NUVEL‐
１A と良
い一致を示していることが分かる．プレート境界
付近に位置する観測局では，NUVEL‐
１A による
速度との差が大きいが，この差は隣接するプレー
トとの相互作用を示唆しており，次節で取り上げ
ることとする．
以 上 よ り，推 定 さ れ た 絶 対 速 度 は，ITRF
２
００
０，ITRF２
００
５
（SLR rescaled）に よ る 速 度 及 び
NUVEL‐
１A による速度（プレート境界付近を除
く）と調和的であることが示された．
３．
５ プレート境界付近に位 置 す る 観 測 局 の
NUVEL‐１A に対する移動速度
Fig．
８ Velocity vectors of Simosato station relative to
the Eurasian plate. The solid arrow shows our
result and the open arrow shows the velocity
vector estimated from Sengoku（１
９
９
８）
.
ここでは，本稿の目的である下里のユーラシア
で西北西４∼５cm/year と推定される．この速度
プレート安定域に対する移動速度について述べ
は，今回求められた速度と方向，大きさともに近
る．繰り返しになるが，ここでいう「ユーラシア
いことから，この辺りの海底下ではプレート間の
プレート安定域に対する移動速度」とは，下里が
固着が強いことが示唆される．
位置しているユーラシアプレートが仮想的に剛体
Table．
５に，今回の結果とこれまでに報告され
だとしたときの速度，すなわち NUVEL‐１A に対
ている下里の移動速度についての他の研究結果を
する速度のことであり，この速度は隣接するプ
示 す．本 稿，Torrence（１９
９８）及 び ASI０
７L０
１
レートとの相互作用を示唆していると考えられ
（Luceri and Bianco，２００
８）は「LAGEOS」SLR
る．
データ，Sengoku（１９９
８）は「AJISAI」SLR デ ー
また，下里以外でも NUVEL‐１A による速度と
タ，ITRF ２
００
０ 及 び ITRF ２００
５（SLR rescaled）
の差が大きい観測局があるため，本論の趣旨とは
（Altamimi et al.,２０
０７）は SLR/VLBI/GPS/DORIS
やや逸れるが，そのような観測局についても簡単
の観測データから求めた速度である．全ての結果
に言及することとする．
において概ね良い一致を示しており，時期の違い
と見られる有意な変化は見られない．Sengoku
３．
５．１ 下里
（１９
９８）による移動速度は，
「AJISAI」SLR データ
Table．
４及び Fig．
７より，下里では，NUVEL‐
１
を使用していることもあり，若干推定誤差が大き
A（ユーラシアプレート）に対する動きとして
いが，LAGEOS 等から求められた速度の推定誤
２
９４°，２．
９cm/year の 速 度 が 求 め ら れ た
差 は０．
１cm/year 程 度 で あ り，方 位 角 が２９
１∼
（Fig．
８）．この速度は，隣接しているフィリピン
２
９４°，大きさが２．
６∼３．
０cm/year の範囲内にあ
海プレートの沈み込みに伴って下里が西北西方向
る．比 較 的 最 近 求 め ら れ た ITRF２０
００や ITRF
に引き込まれていることを意味している．
２０
０５
（SLR rescaled）
，１９
８４年 か ら２
０
０７年 ま で の
Seno et al．
（１９９
３）によると，南海トラフに沈
LAGEOS データを使用した ASI０７L０
１と比べる
み込むフィリピン海プレートの速度は紀伊半島沖
と，方位角は若干大きい結果となった．これらの
− 10 −
Re!estimation of the velocity of Simosato from analysis of LAGEOS SLR data
Table．
５ Velocities of Simosato station by other researches.
速度の違いは，主として移動速度の推定に用いた
ラシアプレートに対する動きとして，２
９
４°，２．
９
衛星，観測局，観測データの期間によるものと考
cm/year が得られた．当庁で長年使用してきた
えられる．
Sengoku（１９９
８）の結果と比べると，速度はやや
小さく，方位角はやや大きくなったが，SLR デー
３．
５．２ その他の観測局
タから求められた他の研究結果と比べて有意な差
中 国 の Changchun
（７２３
７）局，Beijing
（７
２４
９）局，
はなく調和的だと考えられる．今後は，GPS 観
Shanghai
（７
８３
７）局 で は，NUVEL‐１A（ユ ー ラ シ
測等，他の観測手法による結果との比較も行って
アプレート安定域）に対して東から東南東の方向
いきたい．
に約１cm/year の速度が検出された．ユーラシア
移動速度を精度良く決定するためには，長期に
プレートの東方にはユーラシアプレートよりやや
わたる安定した観測が重要であることは言うまで
東向きの速度が大きいアムールプレートの存在が
もない．アジアを代表する SLR 観測局である下
示唆されており（例えば，Heki et al., １９
９９），今
里水路観測所において，今後も継続して観測が
回の結果はそれと調和的である．
行っていくことが重要である．
また，南米の Arequipa
（７４
０３）局でも NUVEL‐
１
A（南米プレート安定域）に対して南東へ１．
２cm
謝
辞
/year の速度が求められた．NUVEL‐
１A による
観測に携わった歴代の下里水路観測所の所長及
と，アレキパ沖ではナスカプレートが約８cm/
び所員の方々のご尽力に深く感謝いたします．ま
year の速度で東北東方向に沈み込んでおり，Are-
た，本稿の作成にあたり，LAGEOS‐
１，LAGEOS
quipa
（７
４０
３）局の移動速度もナスカプレートの沈
‐２のグローバルデータを使用させていただきま
み込みの影響を受けて東北東方向に移動している
した．また，SLR 解析では，NASA/GSFC の GEO-
と推定されるが，SLR データによる速度は南東方
DYN−II/SOLVE を使用させていただきました．
向を示しており方向が異なっている．今回得られ
本稿の改訂にあたっては，匿名の査読者から大変
た速度は，３．
４で述べたように ITRF２０
００や ITRF
有益なコメントをいただきました．記して感謝い
２
００５
（SLR rescaled）
とも異なっていることから，
たします．
速度の差異については今後注意深く検討していき
たい．
要
旨
下里水路観測所では，１９
８２年以来，測地衛星
４ まとめ
「AJISAI」や「LAGEOS」等の SLR 観測を実施し
１９
９３年 か ら２
００
７年 ま で の「LAGEOS」SLR
ている．今回，１９
９３年か ら２
０
０７年 ま で の１
５年
データを用いて下里のユーラシアプレートに対す
間の「LAGEOS」SLR データを用いて，下里の移
る速度の再計算を行った．その結果，下里のユー
動速度の再計算を行った．その結果，ユーラシア
− 11 −
Mariko SATO, Hiroshi MATSUSHITA, Masayuki FUJITA and Arata SENGOKU
プレート安定域に対して，２９
４°の方向に２．
９cm
/year という速度が得られた．この速度はフィリ
daf.mt.asi.it/html_old/ASImed/solution.html.
松下
タ解析におけるモデルの更新とその評価，海洋情
ピン海プレートの沈み込み速度に近く，この辺り
のプレート間カップリングが強いことを示唆して
いる．本結果は，これまでに SLR データから求
優・藤田雅之・佐藤まりこ（２
０
０
５a）
：SLR デー
報部技報，２３，７３‐７７．
松下
優・藤田雅之・河合晃司・長岡
継（２
０
０
５b）
：
SLR グローバル解析により検出された紀伊半島南
められた他の研究結果と概ね調和的である．
東沖地震による地殻変動，日本地球惑星科学連合
２００５年大会予稿集 ．
参考文献
McCarthy, D. D.
（１
９
９
６）
，IERS Conventions
（１
９
９
６）
，IERS
Altamimi, Z., X. Collilieux, J. Legrand, B. Garayt, and C.
Technical Note,21．
Boucher（２
０
０
７）
，ITRF２
０
０
５：A new release of the
Sengoku, A（１
９
９
８）
，A plate motion study using Ajisai
International Terrestrial Reference Frame based on
SLR data., Earth Planets Space,50，６
１
１
‐
６
２
７．
time series of station positions and Earth Orienta-
Sengoku, A, M. Fujita, K. Matsumoto, K. Terai,and M.
tion Parameters, J. Geophys. Res., 112，B０
９
４
０
１，
Sasaki（１
９
９
９）
，Determination of Precise Positions
doi：１
０．１
０
２
９/２
０
０
７JB０
０
４
９
４
９．
of the Mainlands and Isolated Islands in the Japanese Territory‐A Review of Marine Geodetic Con-
Boucher, C., Z. Altamimi, P. Sillard, and M Feissel−Vernie（２
０
０
４），The
ITRF２
０
０
０，IERS
Technical
trol Network Deployed by the Hydrographic Department‐，Rep. Hydrogr. Res., 35，７
９
‐
９
５．
Note,31．
DeMets, C., R. G. Gordon, D. F. Argus, and S. Stein
仙石新・加藤剛・福良博子・小山薫・藤田雅之・矢吹
（１
９
９
４）
，Effect of Recent Revisions to the Geomag-
哲一朗（２
０
０
０）
：LageosI 及び LageosII の long−arc
netic Reversal Time Scale on Estimates of Current
解析による本土海洋測地基準点座標決定と日本測
Plate Motions, Geophys. Res. Lett., 21，２
１
９
１
‐
地系と世界測地系の変換パラメータの再決定，水
２
１
９
４．
路部研究報告，３６，９７‐１０９．
Eddy, W. F., J. J. McCarthy, D. E. Pavlis, J. A. Marshall,
Seno, T., S. Stein, and A. E. Gripp（１
９
９
３）
，A Model for
S. B. Luthke, L. S. Tsaoussi, G. Leung, and D. A.
the Motion of the Philippine Sea Plate Consistent
Williams（１
９
９
０）
，GEODYN!II System Operations
with NUVEL‐
１A and Geological Data, J. Geophys.
Manual, Vol.1-5，Contractor Report, ST Syst.Corp.,
Res., 98，１
７
９
４
１
‐
１
７
９
４
８．
Lanham, Md.
Smith, D. E., R. Kolenkiewicz, P. J. Dunn, and M. H. Tor-
藤田雅之・仙石新（１
９
９７）
：NASA の衛星データ解析
rence（１
９
７
９）
，The measurement of fault motion by
ソフトウエア「GEODYN‐II」の導入，水路部技
報，１５，５‐１０．
satellite laser ranging, Tectonophysics,52，５
９
‐
６
７．
Tapley, B. D., B. E. Schutz, and R. J. Eanes（１
９
８
５）
，Sta-
藤田雅之・久保岡俊宏・池田信広（１
９
９
８）
：GEODYN−
tion coordinates, baselines, and Earth rotation from
II/SOLVE 対話型解析システム（GDIS）による SLR
LAGEOS Laser Ranging：１
９
７
６
‐
１
９
８
４，J. Geophys.
６，１０８‐１１７．
データ解析，水路部技報，１
Res., 90，９
２
３
５
‐
９
２
４
８.
Hedin, A. E．
（１
９
８
７）
，MSIS‐
８
６Thermospheric Model,
Tapley, B. D., D. P. Chambers, S. Bettadpur, and J. C.
J. Geophys. Res., 34，１
３
‐
２
４．
Ries（２
０
０
３）
，Large scale ocean circulation from the
Heki, K., S. Miyazaki, H. Takahashi, M. Kasahara, F. Ki-
GRACE
GGM０
１geoid.
Geophys.
Res.
Lett.,
mata, S. Miura, N. Vasilenco, A. Ivashchenco and K.
30，２
１
６３，doi：１
０．１
０
２
９/２
０
０
３GL０
１
８
６
２
２．
An（１
９
９
９）
，The Amurian plate motion and current
Torrence, M（１
９
９
８）
，Tectonic Plate Motion, http : //
cddisa.gsfc.nasa.gov/９
２
６/slrtecto.html.
plate kinematics in Eastern Asia, J. Geophys. Res.,
104，２
９
１
４
７
‐
２
９
１
５
５.
Ullman, R. E.（１
９
９
２）
，SOLVE Program, NASA Contract
Luceri, V. and G. Bianco（２
０
０
８）
，Earth Orientation and
site coordinates/velocities from the Matera Space
Geodesy Center solution ASI０
７L０
１,http : //geo-
− 12 −
Report NAS５
‐
２
９
３
９
３，Task５
０
３．