038.ニユージーランドの測地系2000
2002年11月


1.はじめに
 ニユージーランドではITRF96に立脚した新しい測地系New Zealand Geodetic Datum 2000(NZGD2000)が1999年8月25日に公開されている。これは進展する社会の要望にこたえ最新の技術によって21世紀の測地系を実現すべく、古いNew Zealand Geodetic Datum 1949(NZGD49)を更新したもので、その設定の手順や結果はニユージーランド土地情報局(Land Information New Zealand,LINZ)のHome Page(http://www.linz.gov.nz)に掲載されている。ニユージーランドの新測地系は、特定の元期における座標値とともにその変化率も与えるなど、極めて斬新な試みに充ちていて今のところ世界のうちでも最も先端的な測地系であろう。
本稿では、文献

Merrin Pearse :Realisation of the New Zealand Geodetic Datum 2000,OSG Technical Report 5,1 June 2000.

によりながら、ニユージーランド測地系2000(NZGD2000)を紹介する。
2.ニユージーランド測地系2000(NZGD2000)の特性
まずニユージーランド測地系2000(NZGD2000)のもつ特徴について述べる。
地心測地系である。以前のNZGD49はその中心が地心から200mほどずれていたがこのことが修正されている。
国際測地学協会の決定した最新の楕円体であるGRS1980楕円体を採用している。
新測地系は国際地心測地システム(International Terrestrial Reference System,ITRS)の考えを元期2000.0で国際地心測地系96(International Terrestrial Reference Frame 96,ITRF96)として実現している。ITRF96はWGS84(G873)と数cmしか異ならない。
各測地点の変動速度もITRSに立脚して算出されている。これによりユーザーは任意の元期以外の時点での最新の測地点の位置を算出して利用できる。
NZGD2000の各測地点の座標は新観測によって変えられるものである。
3.GPSデーターの集約
 1992年に、ニユージーランド測地局(後に土地情報局LINZに併合)はGPS国家網の試験観測を始めた。1993年には30点観測した。これはNZGD2000Gの0および1等点である。0等点は3点(AUCK,CHAT,WGTN)ある。1等点は28点ある。図1にその分布を示す。1994〜1998年の間も観測した。これらの観測では24時間観測を3ケ日続けている。1等点は2〜4等点でうずめられる。2等点は70km間隔、3等点は20km間隔、4等点は1km間隔である。3〜4等点は全国をカバーしているわけではない。

□ NZGD2000網のO〜1等点

4.集約したGPS観測値の整約計算
 NZGD2000確立のために1996.5年元期のITRF96座標値の計算をおこなうことをキャンベラ大学に委託した。解析はGAMIT/GLOBK/GLORGによる。座標値を計算するうえでIGS点をおおいに利用した。
 毎日のGAMIT解では観測点の位置とともに、衛星軌道、地球自転パラメーター、を算出する。これをGLOBKでまとめてひとつのシステムとする。この解をITRF96解にするためにGLORGを使う。これによりスケール、回転、移行の定誤差がなくなる。この計算のさいにIGS点の位置と移動速度の数値を使った。解析の特徴は次のごとくである。
解析ソフトとしてはGAMIT/GLOBKを使った。
すべての基線の解析は一重差によっている。
対流圏モデルは2時間ごとに設定している。
電離層のモデルは2時間ごとに設定している。
サイクルスリップの検出にCVIEWというソフトを使っている。
各観測ごとに網平均し、座標値のほかに共分散行列をフアイルとして出力している。
高さは潮汐の影響のないシステムとして算出してある。

1996.5 5元期のITRF座標

Geodctic
Code
Name X (m) Y (m) Z(m)
AUCK Whangaparaoa -5105680.990 461564.048 -3782181.773
CHAT Chatham Island -4590670.893 -275483.014 -4404596.787
WGTN Wellington Airport -4777269.350 434269.976 -4189484.652
1004 Mt York -4371448.508 950019.355 -4531598.951
1017 Hyde Rock -4408673.299 841182.570 -4518904.677
1103 Mt Horrible -4509242.347 709568.151 -4440278.862
1153 Isolated Hill -4660964.413 571445.661 -4302316.292
1181 Mt Murchison -4727392.500 622413.005 -4224184.004
1215 Eringa -4794050.724 364491.783 -4177890.168
1231 Mt Stewart -4860522.601 383529.069 -4098473.496
1259 Huirangi -4929040.037 498221.579 -4004033.171
1273 Okahuatiu -4989460.451 191252.559 -3955756.527
1305 Te Pohue -5042731.036 140230.637 -3890300.265
1314 Maketu -5039298.757 311182.284 -3884430.072
1344 Ruahine -5128814.165 401970.916 -3758286.418
1361 Mahora -5138072.408 560948.499 -3724886.327
1367 Pekearenga -5167214.344 496239.215 -3693852.201
1394 Te Paki -5222589.733 662414 397 -3589435.368
1420 Mt Greenland -4616400.196 745241.589 -4324328.213
1501 Bluff Hill -4922647.771 265115.050 -4033578.818
2085 Marotiri -4977955.116 355512.020 -3959577.453
5508 Windsor Castle -4590826.049 584593 .009 -4374757.107
5509 Three Sisters -4303267.130 894811.292 -4606633.315
6731 Haast -4521641.725 878624.013 -4396964.293
A31C Milford -4442644.982 950468.863 -4461626.486
A33D Te Herunga -5041355.010 441059.373 -3869624.530
A70X Parapara -4802026.173 617521.260 -4138427.835
B03W Puysegur Point -4305507.267 1024975.858 -4577337.960
B28C Cape Cambell -4741512.950 480470.749 -4225018.452
OUSD Otago University -4387888.550 733420.870 -4555178.594
WELL Wellington -4780648.764 436507.201 -4185440.238

5.速度モデルの決定
 地質核科学研究所(Geology and Nuclear Science Institute,GNS)に、GPSデータを使った1997GNS速度モデルを作ることを委託した。解析方法は座標変化プログラムADJCOORDによる(Bibby,1982,Crook,1992)。詳しい報告はBeavan(1998)にある。約30点について1996.5のITRF96の位置と速度場の解をもとめたが、ある点では大きな誤差が出た。NZGD2000の座標値は2000.0元期で与えてある。

@ 0〜1等点の1996.5元期の座標値を速度場モデルによって2000.0に化成した。95%信頼楕円の軸長は水平0.01m,高さ0.02mである。
A 2〜5等点は0〜1等点を与点として網平均する。SNAP(Crook,1995)による。この際の重量はそれぞれの点のアプリオリな誤差により求める。この誤差を表2に示した。観測のアプリオリ誤差は定数項と距離に依存する項とからなる。


観測のアプリオリ誤差

CONSTANT Line Length Error
E mm N mm U mm E ppm N ppm U ppm
1st Order   1.2 1.2 1.5   0.004 0.004 0.15
2st Order 1.2 1.2 1.5 0.4 0.4 1.5
3st Order 4 4 5 1.2 1.2 5
4st Order 4 4 5 4 4 15
5st Order 4 4 5 12 12 50

点数は
2等点  313点
3等点 1339点
4等点 1714点

Beaven(1998)とMorgan・Pearse(1999)の速度場の比較。
95%信頼楕円はBeavenの推定値に付してある。


元期2000.0でのITRF96でのNZGD2000座標値

Geodctic Code Name Order Latitude
(°’”)
Longitude
(°’”)
Ellipsoidal
Height(m)
AUCK Whangaparaoa Zero 36 36 10.2402 S 174 50 03.7880 E 132.711
CHAT Chatham Island Zero 43 57 20.8326 S 176 33 57.0227 W 57.9811
WGTN Wellington Airport Zero 41 19 24.4454 S 174 48 21.2186 E 26.073
1004 Mt York First 45 33 43.6111 S 167 44 20.1266 E 411.196
1017 Hyde Rock First 45 23 15.5209 S 169 11 51.7276 E 1680.809
1103 Mt Horrible First 44 24 02.0503 S 171 03 26.4405 E 397.157
1153 Isolated Hill First 42 41 14.7040 S 173 00 37.0016 E 405.505
1181 Mt Murchison First 41 43 44.6973 S 172 29 58.2847 E 1486.646
1215 Eringa First 41 10 48.5107 S 175 39 07.7910 E 590.791
1231 Mt Stewart First 40 14 24.7159 S 175 29 17.9227 E 143.609
1259 Huirangi First 39 08 02.3985 S 174 13 41.5683 E 263.041
1273 Okahuatiu First 38 34 30.5480 S 177 48 17.4692 E 323.407
1305 Te Pohue First 37 49 28.3511 S 178 24 25.5711 E 360.454
1314 Maketu First 37 45 34.0792 S 176 27 59.0730 E 95.727
1344 Ruahine First 36 19 58.9983 S 175 31 06.9717 E 438.003
1361 Mahora First 35 57 43.5863 S 173 46 09.8974 E 164.970
1367 Pekearenga First 35 37 02.0783 S 174 30 51.6930 E 174.403
1394 Te Paki First 34 27 59.7095 S 172 46 17.0730 E 354.050
1420 Mt Greenland First 42 57 11.6870 S 170 49 46.7693 E 919.307
1501 Bluff Hill First 39 28 44.3468 S 176 55 02.0846 E 119.271
2085 Marotiri First 38 36 57.7701 S 175 54 54.0969 E 760.272
5508 Windsor Castle First 43 34 53.4132 S 172 44 34.9627 E 335.355
5509 Three Sisters First 46 32 12.9470 S 168 15 12.3806 E 176.343
6731 Haast First 43 51 38.9427 S 169 00 12.9269 E 14.412
A31C Milford First 44 40 24.6238 S 167 55 26.6376 E 9.546
A33D Te Herunga First 37 35 21.7825 S 175 00 00.0836 E 318.912
A70X Parapara First 40 42 46.8026 S 172 40 19.9543 E 169.539
B03W Puysegur Point First 46 09 23.0087 S 166 36 33.5751 E 44.265
B28C Cape Cambell First 41 44 56.5656 S 173 12 49.7107 E 254.535
OUSD Otago University First 45 52 10.2057 S 170 30 39.3147 E 26.197
WELL Wellington First 41.16.29.6122 S 174 46 58.6319 E 37.686

5.まとめ
 NZGD2000の0〜1等点の初期座標値を算出し、これより2〜4等点がどのように決定されたかを述べた。座標値の元期は2000.0である。新座標値はLINZのサイトにて見出すことができる。