Page 1 Development of Precise Geoid Model for the Establishment

，

， ， ，

Ｃｈｉｎａ Ｏｃｅａｎ Ｅｎｇｉｎｅｅｒｉｎｇ Ｖｏｌ． ２３ Ｎｏ． ４ ｐｐ． ６７９ － ６９４

?

２００９ Ｃｈｉｎｅｓｅ Ｏｃｅａｎ Ｅｎｇｉｎｅｅｒｉｎｇ Ｓｏｃｉｅｔｙ ＩＳＳＮ ０８９０５４８７

Ｄｅｖｅｌｏｐｍｅｎｔ ｏｆ Ｐｒｅｃｉｓｅ Ｇｅｏｉｄ Ｍｏｄｅｌ ｆｏｒ ｔｈｅ Ｅｓｔａｂｌｉｓｈｍｅｎｔ ｏｆ Ｃｏｎｓｉｓｔｅｎｔ Ｈｅｉｇｈｔ Ｓｙｓｔｅｍ ｉｎ Ｇｅｏｇａ Ｇｒａｎｄ Ｂｒｉｄｇｅ Ｃｏｎｓｔｒｕｃｔｉｏｎ Ａｒｅａ

，

ＤｏｎｇＨａ ＬＥＥ ＨｏｎｇＳｉｃ ＹＵＮ１ ａｎｄ Ｈｅ ＨＵＡＮＧ

，Ａｒｃｈｉｔｅｃｔｕｒａｌ ａｎｄ Ｅｎｖｉｒｏｎｍｅｎｔａｌ Ｓｙｓｔｅｍ Ｅｎｇｉｎｅｅｒｉｎｇ，Ｓｕｎｇｋｙｕｎｋｗａｎ Ｕｎｉｖｅｒｓｉｔｙ， Ｓｕｗｏｎ ４４０ ７４６，Ｋｏｒｅａ （Ｒｅｃｅｉｖｅｄ １ Ｓｅｐｔｅｍｂｅｒ ２００８；ｒｅｃｅｉｖｅｄ ｒｅｖｉｓｅｄ ｆｏｒｍ ２４ Ａｐｒｉｌ ２００９；ａｃｃｅｐｔｅｄ １３ Ｊｕｎｅ ２００９）

Ｄｅｐａｒｔｍｅｎｔ ｏｆ Ｃｉｖｉｌ

ＡＢＳＴＲＡＣＴ Ｔｈｅ ｎａｔｉｏｎａｌ ｂｅｎｃｈｍａｒｋｓ ｏｎ ｉｓｌａｎｄｓ ｗｅｒｅ ｍｏｓｔｌｙ ｅｓｔａｂｌｉｓｈｅｄ ｂｙ ｔｒｉｇｏｎｏｍｅｔｒｉｃ ｌｅｖｅｌｉｎｇ ｉｎ Ｋｏｒｅａ． Ｔｈｉｓ ｍｅｔｈｏｄ ｒｅｓｕｌｔｓ ｉｎ ｉｎａｃｃｕｒａｃｙ ｗｈｉｃｈ ｉｓ ａ ｓｅｒｉｏｕｓ ｐｒｏｂｌｅｍ ｉｎ Ｇｅｏｇａ Ｇｒａｎｄ Ｂｒｉｄｇｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ｗｏｒｋ ｔｈａｔ ｔｒｉｅｄ ｔｏ ｌｉｎｋ ｔｈｅ ｍａｉｎｌａｎｄ ａｎｄ ｔｈｅ ｉｓｌａｎｄｓ． Ｔｈｅ Ｇｅｏｇａ Ｇｒａｎｄ Ｂｒｉｄｇｅ ＰｕｓａｎＧｅｏｊｅ ｆｉｘｅｄ ｌｉｎｋ ｐｒｏｊｅｃｔ ｗａｓ ｓｅｌｅｃｔｅｄ ａｓ ｔｈｅ ｓｔｕｄｙ ａｒｅａ ａ ｈｕｇｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ｗｏｒｋ ｉｎ Ｋｏｒｅａ ｔｈａｔ ｗｉｌｌ ｃｏｎｎｅｃｔ ｔｈｅ ｍａｉｎｌａｎｄ Ｐｕｓａｎ ａｎｄ ａｎ ｉｓｌａｎｄ Ｇｅｏｊｅ ｉｓｌａｎｄ ． Ｈｏｗｅｖｅｒ ｔｈｅ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔｓ ｉｓｓｕｅｄ ａｔ ｂｅｎｃｈｍａｒｋｓ ＪＩＮＨ ａｎｄ ＧＯＥＪ ｗｅｒｅ ｎｏｔ ｃｏｎｓｉｓｔｅｎｔ ｂｅｃａｕｓｅ ｔｈｅｙ ｄｉｄ ｎｏｔ ｒｅｆｅｒ ｔｏ ｔｈｅ ｓａｍｅ ｚｅｒｏ ｐｏｉｎｔ ｗｈｉｃｈ ｗｏｕｌｄ ｍａｋｅ ｔｈｅ ｌｉｎｋｉｎｇ ｏｆ ｔｈｅ ｓｅｃｔｉｏｎｓ ｐｒｏｂｌｅｍａｔｉｃ． Ｔｈｉｓ ｐａｐｅｒ ｉｎｔｒｏｄｕｃｅｓ ｔｈｅ ｐｒｅｃｉｓｅ ｌｏｃａｌ ｇｅｏｉｄ ａｓ ａ ｖｅｒｔｉｃａｌ ｄａｔｕｍ ｆｏｒ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ｉｎ ｏｒｄｅｒ ｔｏ ｅｓｔａｂｌｉｓｈ ａ ｃｏｎｓｉｓｔｅｎｔ ｈｅｉｇｈｔ ｓｙｓｔｅｍ． Ｔｏ ｄｅｔｅｒｍｉｎｅ ｔｈｅ ｐｒｅｃｉｓｅ ｌｏｃａｌ ｇｅｏｉｄ ｆｏｒ ｔｈｅ ｃｏｎ ｓｔｒｕｃｔｉｏｎ ａｒｅａ ｗｅ ｆｉｒｓｔｌｙ ｄｅｖｅｌｏｐｅｄ ａ ｐｒｅｃｉｓｅ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ｆｏｒ Ｋｏｒｅａ ａｎｄ ｉｔｓ ａｄｊｏｉｎｉｎｇ ｓｅａｓ ａｓ ａ ｗｈｏｌｅ． Ｔｈｉｓ ｇｒａｖｉｍｅｔ ｒｉｃ ｇｅｏｉｄ ｗａｓ ｄｅｖｅｌｏｐｅｄ ｂｙ ｕｓｅ ｏｆ ａｌｌ ａｖａｉｌａｂｌｅ ｇｒａｖｉｔｙ ｄａｔａ ｉｎｃｌｕｄｉｎｇ ｓｕｒｆａｃｅ ａｎｄ ｓａｔｅｌｌｉｔｅ ｄａｔａ ｏｎ ｌａｎｄ ａｎｄ ｏｎ ｔｈｅ ｏｃｅａｎ． Ｔｈｅ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ｗａｓ ｃｏｍｐｕｔｅｄ ｂｙ ｓｐｈｅｒｉｃａｌ ｆａｓｔ ｆｏｕｒｉｅｒ ｔｒａｎｓｆｏｒｍ ｗｉｔｈ ｍｏｄｉｆｉｅｄ Ｓｔｏｋｅｓ′ ｋｅｒｎｅｌｓ． Ｔｈｅ ｒｅｍｏｖｅｒｅｓｔｏｒｅ ｔｅｃｈｎｉｑｕｅ ｗａｓ ｕｓｅｄ ｔｏ ｅｌｉｍｉｎａｔｅ ｔｈｅ ｔｅｒｒａｉｎ ｅｆｆｅｃｔｓ ｂｙ ｕｓｅ ｏｆ ｔｈｅ ＲＴＭ ｒｅｄｕｃｔｉｏｎ ａｎｄ ｔｏ ｄｅｔｅｒｍｉｎｅ ｔｈｅ ｒｅｓｉｄｕａｌ ｇｅｏｉｄ ｂｙ ｃｏｍｂｉｎｉｎｇ ｔｈｅ ＧＧＭ０２Ｓ ／ ＥＧＭ９６ ｇｅｏｐｏｔｅｎｔｉａｌ ｍｏｄｅｌ ｆｒｅｅａｉｒ ｇｒａｖｉｔｙ ａｎｏｍａｌｉｅｓ ａｎｄ ｈｉｇｈｒｅｓｏｌｕｔｉｏｎ ＤＥＭ ｄａｔａ． Ｆｉｎａｌｌｙ ｔｈｅ ｇｒａｖｉｍｅｔｒｉｃ ｍｏｄｅｌ ｗａｓ ｆｉｔｔｅｄ ｔｏ ｔｈｅ ｇｅｏｉｄ ｈｅｉｇｈｔｓ ｏｂｔａｉｎｅｄ ｆｒｏｍ ＧＰＳ ａｎｄ ｔｉｄｅ ｏｂｓｅｒｖａｔｉｏｎｓ Ｎ ＧＰＳ ／ Ｔｉｄｅ ｂｙ ｌｅａｓｔ ｓｑｕａｒｅ ｃｏｌｌｏｃａｔｉｏｎ ｔｏ ｐｒｏｖｉｄｅ ｔｈｅ ｆｉｎａｌ ＧＰＳｃｏｎｓｉｓｔｅｎｔ ｌｏｃａｌ ｐｒｅｃｉｓｅ ｇｅｏｉｄ ｍｏｄｅｌ． Ｔｈｅ ｐｏｓｔｆｉｔ ｅｒｒｏｒ ｓｔｄ． ｄｅｖ． ｏｆ ｔｈｅ ｆｉｎａｌ ｇｅｏｉｄ ｔｏ ｔｈｅ Ｎ ＧＰＳ ／ Ｔｉｄｅ ｄｅｒｉｖｅｄ ｆｒｏｍ ＧＰＳ ａｎｄ ｔｉｄｅ ｏｂｓｅｒｖａｔｉｏｎｓ ｗａｓ ± ２ ． ２ ｃｍ ｆｏｒ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ． Ｗｅ ｓｏｌｖｅｄ ｔｈｅ ｈｅｉｇｈｔ ｉｎｃｏｎｓｉｓｔｅｎｃｙ ｐｒｏｂｌｅｍ ｂｙ ｃａｌｃｕｌａｔｉｎｇ ｔｈｅ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔ ｏｆ ｔｈｅ ｂｅｎｃｈｍａｒｋｓ ａｎｄ ｔｈｅ ｃｏｎｔｒｏｌ ｐｏｉｎｔｓ ｕｓｉｎｇ ｔｈｅ ｆｉｎａｌ ｇｅｏｉｄ ｍｏｄｅｌ． Ａｌｓｏ ｔｈｅ ｈｉｇｈｌｙ ａｃｃｕｒａｔｅ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔ ｗａｓ ｅｓｔｉｍａｔｅｄ ｔｈｒｏｕｇｈ ｔｈｅ ＧＰＳ ／ ｌｅｖｅｌｉｎｇ ｔｅｃｈｎｉｑｕｅ ｂｙ ａｐｐｌｙｉｎｇ ｔｈｅ ｄｅｖｅｌｏｐｅｄ ｌｏｃａｌ ｐｒｅｃｉｓｅ ｇｅｏｉｄ． Ｔｈｅｒｅｆｏｒｅ ｔｈｅ ｐｒｅｃｉｓｅ ｌｏｃａｌ ｇｅｏｉｄ ｉｓ ｅｘｐｅｃｔｅｄ ｔｏ ｉｍｐｒｏｖｅ ｔｈｅ ｑｕａｌｉｔｙ ｏｆ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ｐｒｏｃｅｄｕｒｅ ｏｆ ｔｈｅ Ｇｅｏｇａ Ｇｒａｎｄ Ｂｒｉｄｇｅ．

，

（

（

）

）

（ ）

，

，

（

）

，

，

，

，

，

（ （

，

，

：

）

，

）

，

；

；

；

；

Ｋｅｙ ｗｏｒｄｓ ｐｒｅｃｉｓｅ ｌｏｃａｌ ｇｅｏｉｄ ｃｏｎｓｉｓｔｅｎｔ ｈｅｉｇｈｔ ｓｙｓｔｅｍ ｇｒａｖｉｔｙ ｏｂｓｅｒｖａｔｉｏｎｓ ＧＰＳ ／ Ｔｉｄｅ ｏｂｓｅｒｖａｔｉｏｎｓ ｌｅａｓｔ ｓｑｕａｒｅｓ ｆｉｔ ｔｉｎｇ

１ ． Ｉｎｔｒｏｄｕｃｔｉｏｎ

，

Ｉｎ Ｋｏｒｅａ ｎａｔｉｏｎａｌ ｂｅｎｃｈｍａｒｋｓ ｗｅｒｅ ｅｓｔａｂｌｉｓｈｅｄ ａｎｄ ａｒｅ ｍａｉｎｔａｉｎｅｄ ｂｙ ｔｈｅ Ｎａｔｉｏｎａｌ Ｇｅｏｇｒａｐｈｉｃ

（ ） ｐｏｓｅｓ：ｓｕｒｖｅｙｉｎｇ，ｍａｐｐｉｎｇ，ｒｅｓｅａｒｃｈ ａｎｄ ｃｏｎｓｔｒｕｃｔｉｏｎ ．

Ｉｎｆｏｒｍａｔｉｏｎ Ｉｎｓｔｉｔｕｔｅ ＮＧＩＩ ． ＮＧＩＩ ｉｓｓｕｅｓ ｔｈｅ ｅｘａｃｔ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔｓ ｏｆ ｂｅｎｃｈｍａｒｋｓ ｆｏｒ ｖａｒｉｏｕｓ ｐｕｒ

，ｔｈｅ ｂｅｎｃｈｍａｒｋｓ ａｒｅ ｉｍｐｏｒｔａｎｔ ｐｏｉｎｔｓ ｓｉｎｃｅ ｔｈｅｙ ａｒｅ ｕｓｅｄ ａｓ ａ ｒｅｆｅｒｅｎｃｅ ｏｆ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔ ｉｎ ｃｏｎｓｔｒｕｃｔｉｏｎ ｗｏｒｋ ． Ｔｈｅｒｅｆｏｒｅ，ｔｈｅ ａｌ ｔｉｔｕｄｅ ｏｆ ｔｈｅｓｅ ｂｅｎｃｈｍａｒｋｓ ｓｈｏｕｌｄ ｂｅ ｍｅａｓｕｒｅｄ ｅｘａｃｔｌｙ ｂｙ ｓｐｉｒｉｔ ｌｅｖｅｌｉｎｇ ． Ｈｏｗｅｖｅｒ，ｔｈｅ ｂｅｎｃｈｍａｒｋｓ ｌｏ Ｅｓｐｅｃｉａｌｌｙ

ｃａｔｅｄ ｏｎ ｉｓｌａｎｄｓ ｗｅｒｅ ｍｏｓｔｌｙ ｅｓｔａｂｌｉｓｈｅｄ ｂｙ ｔｒｉｇｏｎｏｍｅｔｒｉｃ ｌｅｖｅｌｉｎｇ ｍｅｔｈｏｄ ． Ｔｈｉｓ ｍｅｔｈｏｄ ｙｉｅｌｄｓ ｉｎａｃｃｕ

：

１ Ｃｏｒｒｅｓｐｏｎｄｉｎｇ ａｕｔｈｏｒ． Ｅｍａｉｌ ｙｈｓ＠ ｇｅｏ． ｓｋｋｕ． ａｃ． ｋｒ

６８０

， （），２００９，６７９ － ６９４

ＤｏｎｇＨａ ＬＥＥ ｅｔ ａｌ ． ／ Ｃｈｉｎａ Ｏｃｅａｎ Ｅｎｇｉｎｅｅｒｉｎｇ ２３ ４

， ｃａｌ ａｎｇｌｅ （ｔｈｅ ｉｎｃｉｄｅｎｔ ａｎｇｌｅ ｏｆ ｐｒｉｓｍ）ｏｆ １０° （Ｒｈｏ ｅｔ ａｌ ． ，２００１）． Ｔｈｉｓ ｑｕａｎｔｉｔｙ ｃａｎ ｂｅ ａ ｓｅｒｉｏｕｓ ｒａｔｅ ｌｅｖｅｌｉｎｇ ｒｅｓｕｌｔｓ ｕｐ ｔｏ ５ ． ６ ｍｍ ／ ｋｍ ｗｈｅｎ ｕｓｉｎｇ ＥＤＭ ａｎｄ ｏｂｓｅｒｖｉｎｇ ｓｉｍｕｌｔａｎｅｏｕｓｌｙ ｗｉｔｈ ｔｈｅ ｖｅｒｔｉ ｐｒｏｂｌｅｍ ｉｎ ｔｈｅ ｂｒｉｄｇｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ｗｏｒｋ ．

（ ）

（ ）

Ｆｉｇ． １ ． Ｒｅｌａｔｉｏｎｓｈｉｐ ｂｅｔｗｅｅｎ ｔｈｅ ｈｅｉｇｈｔ ｓｙｓｔｅｍｓ ｏｎ ｔｈｅ ｍａｉｎｌａｎｄ ＪＩＮＨ ａｎｄ ｏｎ ｔｈｅ ｉｓｌａｎｄ ＧＥＯＪ ．

，

（

）

Ｔｈｅ ｏｎｇｏｉｎｇ ｃｏｎｓｔｒｕｃｔｉｏｎ ｐｒｏｊｅｃｔ Ｇｅｏｇａ Ｇｒａｎｄ Ｂｒｉｄｇｅ ＰｕｓａｎＧｅｏｊｅ ｆｉｘｅｄ ｌｉｎｋ ｐｒｏｊｅｃｔ ｔｏ ｃｏｎ

（ ） （ ） ， ｔｉｏｎｓ：４ ． ５ ｋｍ ｃａｂｌｅｓｔａｙｅｄ ｇｉｒｄｅｒ ｂｒｉｄｇｅ，３ ． ７ ｋｍ ｉｍｍｅｒｓｅｄ ｔｕｎｎｅｌ，ａｎｄ ２５ ． ７ ｋｍ ｃｏｎｎｅｃｔｉｏｎ ｒｏａｄ ． Ｆｏｒ ｌｉｎｋｉｎｇ ｏｆ ｔｈｅ ｓｅｃｔｉｏｎｓ，ｃｏｎｔｒｏｌ ｐｏｉｎｔｓ ｗｅｒｅ ｉｎｓｔａｌｌｅｄ ａｔ ｅａｃｈ ｓｅｃｔｉｏｎ ａｎｄ ｔｈｅｉｒ ｈｅｉｇｈｔｓ ｗｅｒｅ ｃｏｍｐｕｔ ｅｄ ｂｙ ａｐｐｌｙｉｎｇ ｔｈｅ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔ ｏｆ ｔｈｅ ｎａｔｉｏｎａｌ ｂｅｎｃｈｍａｒｋｓ，ＪＩＮＨ ａｎｄ ＧＥＯＪ ． Ｈｏｗｅｖｅｒ，ｓｏｍｅ ｄｉｆｆｅｒｅｎｃｅ ｏｖｅｒ ３０ ｃｍ ｉｎ ｔｈｅ ｈｅｉｇｈｔ ｓｙｓｔｅｍ ａｓ ｓｈｏｗｎ ｉｎ Ｆｉｇ ． １ ｂｅｔｗｅｅｎ ｔｈｅ ｍａｉｎｌａｎｄ （ＪＩＮＨ）ａｎｄ ｔｈｅ ｉｓｌａｎｄ （ ＧＥＯＪ）ｗａｓ ｒｅｖｅａｌｅｄ ｉｎ ｔｈｅ ｐｒｏｃｅｓｓ ｏｆ ｔｈｅ ｒｏｕｔｅ ｓｕｒｖｅｙｉｎｇ ｆｏｒ ｒｏａｄ ｄｅｓｉｇｎ ｂｙ Ｄａｅｗｏｏ Ｅ＆Ｃ Ｃｏ ． Ｌｔｄ ． （ Ｄａｅｗｏｏ Ｅ＆Ｃ Ｃｏ ． Ｌｔｄ ． ，２００３）． Ｔｈｉｓ ｄｉｆｆｅｒｅｎｃｅ ｒｅｇａｒｄｅｄ ａｓ ｔｈｅ ｈｅｉｇｈｔ ｅｒｒｏｒ ｂｅｔｗｅｅｎ ｔｈｅ ｂｅｎｃｈ ｍａｒｋｓ （ＪＩＮＨ ａｎｄ ＧＥＯＪ）ｗｅｒｅ ｃａｕｓｅｄ ｂｙ ｉｎａｃｃｕｒａｔｅ ｌｅｖｅｌｉｎｇ ｒｅｓｕｌｔ ｏｆ ｂｅｎｃｈｍａｒｋ （ＧＥＯＪ）ｏｎ ｔｈｅ ｉｓ ｌａｎｄ，ｗｈｉｃｈ ｍｅａｎｓ ｔｈａｔ ｔｈｅ ｈｅｉｇｈｔ ｓｙｓｔｅｍ ｂｅｔｗｅｅｎ ｔｈｅ ｍａｉｎｌａｎｄ ａｎｄ ｔｈｅ ｉｓｌａｎｄ ｗｅｒｅ ｎｏｔ ｃｏｎｓｉｓｔｅｎｔ ． Ｉｎ ｏｔｈｅｒ ｗｏｒｄｓ，ｔｈｅ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔｓ ｉｓｓｕｅｄ ａｔ ｂｅｎｃｈｍａｒｋｓ （ＪＩＮＨ ａｎｄ ＧＥＯＪ）ｗｅｒｅ ｎｏｔ ｒｅｆｅｒｒｅｄ ｔｏ ｔｈｅ ｓａｍｅ ｚｅｒｏ ｐｏｉｎｔ，ｗｈｉｃｈ ｗｏｕｌｄ ｍａｋｅ ｔｈｅ ｌｉｎｋｉｎｇ ｏｆ ｔｈｅ ｓｅｃｔｉｏｎｓ ｐｒｏｂｌｅｍａｔｉｃ ． Ｄｉｆｆｅｒｅｎｔ ａｐｐｒｏａｃｈｅｓ ｈａｖｅ ｂｅｅｎ ｓｔｕｄｉｅｄ ｆｏｒ ｆｉｘｉｎｇ ｔｈｅ ｈｅｉｇｈｔ ｉｎｃｏｎｓｉｓｔｅｎｃｙ ｐｒｏｂｌｅｍ ａｍｏｎｇ ｄｉｆｆｅｒｅｎｔ ｈｅｉｇｈｔ ｓｙｓｔｅｍｓ ｉｎ ｔｈｅ ｃｏｎｔｅｘｔ ｏｆ ｈｅｉｇｈｔ ｓｙｓｔｅｍ ｕｎｉｆｉｃａｔｉｏｎ （Ｚｈａｎｇ ｅｔ ａｌ ． ，２００９）． Ｔｈｅ ｍａｉｎ ｉｓｓｕｅ ｏｆ ｈｅｉｇｈｔ ｓｙｓｔｅｍ ｕｎｉｆｉｃａｔｉｏｎ ｉｓ ｔｏ ｄｅｔｅｒｍｉｎｅ ｔｈｅ ｐｏｔｅｎｔｉａｌ ｄｉｆｆｅｒｅｎｃｅ ａｍｏｎｇ ｄｉｆｆｅｒｅｎｔ ｈｅｉｇｈｔ ｓｙｓｔｅｍｓ （ Ｒｕｍ ｍｅｌ，２０００）． Ｇｅｏｄｅｔｉｃ ｌｅｖｅｌｉｎｇ ａｎｄ ｏｃｅａｎ ｌｅｖｅｌｉｎｇ ｍｅｔｈｏｄ ａｒｅ ｗｅｌｌｋｎｏｗｎ ｆｏｒ ｔｈｅ ｈｅｉｇｈｔ ｓｙｓｔｅｍ ｕｎｉｆｉｃａ ｔｉｏｎ ． Ｔｈｅ ｇｅｏｄｅｔｉｃ ｌｅｖｅｌｉｎｇ ｉｓ ａ ｖｅｒｙ ｈｉｇｈｌｙ ａｃｃｕｒａｔｅ ｍｅｔｈｏｄ ｔｏ ｄｅｔｅｒｍｉｎｅ ｔｈｅ ｐｏｔｅｎｔｉａｌ ｄｉｆｆｅｒｅｎｃｅ，ｂｕｔ ｉｔ ｃａｎｎｏｔ ｂｅ ｕｓｅｄ ｔｏ ｕｎｉｆｙ ｔｈｅ ｄｉｆｆｅｒｅｎｔ ｈｅｉｇｈｔ ｓｙｓｔｅｍｓ ｓｅｐａｒａｔｅｄ ｂｙ ｏｃｅａｎｓ ． Ａｎｏｔｈｅｒ ｍｅｔｈｏｄ ｉｓ ｏｃｅａｎ ｌｅｖｅｌｉｎｇ （Ｅｋｍａｎ，１９９９）ｔｈａｔ ｃａｎ ｂｅ ｕｓｅｄ ｔｏ ｄｅｔｅｒｍｉｎｅ ｔｈｅ ｐｏｔｅｎｔｉａｌ ｄｉｆｆｅｒｅｎｃｅ ｂｅｔｗｅｅｎ ｄｉｆｆｅｒｅｎｔ ｈｅｉｇｈｔ ｓｙｓｔｅｍｓ ａｃｒｏｓｓ ｏｃｅａｎ，ｂｕｔ ｉｔ ｈａｓ ａ ｌｏｗ ａｃｃｕｒａｃｙ ｗｉｔｈ ｒｅｓｐｅｃｔ ｔｏ ｔｈｅ ｇｅｏｄｅｔｉｃ ｌｅｖｅｌｉｎｇ ｍｅｔｈｏｄ （Ｚｈａｎｇ ｅｔ ａｌ ． ，２００９）． Ａ ｎｅｗ ｍｅｔｈｏｄ ｆｏｒ ｕｎｉｆｙｉｎｇ ｄｉｆｆｅｒｅｎｔ ｈｅｉｇｈｔ ｓｙｓｔｅｍｓ ａｓ ｐｒｏｐｏｓｅｄ ｂｙ Ｂｕｒａ ｅｔ ａｌ ． （ １９９９） ｉｓ ｔｏ ｕｓｅ ｔｈｅ ＧＰＳ ／ ｌｅｖｅｌｉｎｇ ｄａｔａ ｆｏｒ ａ ｌｏｃａｌ ａｒｅａ ｉｎ ｏｒｄｅｒ ｔｏ ｄｅｔｅｒｍｉｎｅ ｔｈｅ ｐｏｔｅｎｔｉａｌ ｄｉｆｆｅｒｅｎｃｅ ｂｅｔｗｅｅｎ ｌｏｃａｌ ｐｒｅｃｉｓｅ ｇｅｏｉｄ ａｎｄ ｍｅａｎ ｓｅａ ｓｕｒｆａｃｅ，ｗｈｉｃｈ ｉｓ ｄｅｆｉｎｅｄ ｂｙ ａ ｔｉｄａｌ ｇａｕｇｅ ｏｂｓｅｒｖａｔｉｏｎｓ （Ｂｕｒａ ｅｔ ｎｅｃｔ ｔｈｅ ｍａｉｎｌａｎｄ Ｐｕｓａｎ ａｎｄ ｔｈｅ ｉｓｌａｎｄｓ Ｇｅｏｊｅ ｉｓｌａｎｄ ｉｎ Ｋｏｒｅａ ｃｏｎｓｉｓｔｓ ｏｆ ｔｈｒｅｅ ｄｉｆｆｅｒｅｎｔ ｓｅｃ

， （），２００９，６７９ － ６９４

ＤｏｎｇＨａ ＬＥＥ ｅｔ ａｌ ． ／ Ｃｈｉｎａ Ｏｃｅａｎ Ｅｎｇｉｎｅｅｒｉｎｇ ２３ ４

， ；

， ）

６８１

，

，

ａｌ ． ２００４ Ｚｈａｎｇ ｅｔ ａｌ ． ２００９ ． Ｉｎ ｐｒｉｎｃｉｐｌｅ ｉｆ ｗｅ ｃａｎ ｆｉｎｄ ｔｈｅ ｌｏｃａｌ ｍｅａｎ ｓｅａ ｓｕｒｆａｃｅ ｗｈｉｃｈ ｉｓ

， ｍｅａｎ ｓｅａ ｓｕｒｆａｃｅ ｆｏｒ ｄｅｔｅｒｍｉｎｉｎｇ ｔｈｅ ｐｏｔｅｎｔｉａｌ ｄｉｆｆｅｒｅｎｃｅｓ ａｍｏｎｇ ｄｉｆｆｅｒｅｎｔ ｒｅｇｉｏｎｓ （Ｆｅｉ ａｎｄ Ｓｉｄｅｒｉｓ， ２００１；Ｌｉ ｅｔ ａｌ ． ，２００３；Ｚｈａｎｇ ｅｔ ａｌ ． ，２００９；Ａｍｏｓ ａｎｄ Ｆｅａｔｈｅｒｓｔｏｎｅ，２００９）． Ｔｈｅｒｅｆｏｒｅ，Ｄａｅｗｏｏ Ｅ＆Ｃ Ｃｏ ． Ｌｔｄ ． ｒｅｃｏｇｎｉｚｅｄ ｔｈａｔ ｔｈｅ ｈｅｉｇｈｔ ｄｉｆｆｅｒｅｎｃｅｓ ｏｆ ｂｏｔｈ ｂｅｎｃｈｍａｒｋｓ （ＪＩＮＨ ａｎｄ ＧＥＯＪ）ｗｅｒｅ ｃａｕｓｅｄ ｂｙ ｔｈｅ ｉｎｃｏｎｓｉｓｔｅｎｔ ｈｅｉｇｈｔ ｓｙｓｔｅｍ ｐｒｏｂｌｅｍ ａｎｄ ｉｔ ｉｓ ｎｅｃｅｓｓａｒｙ ｔｏ ｄｅｔｅｒ ｉｒｒｅｌｅｖａｎｔ ｔｏ ｔｈｅ ｒｅｆｅｒｅｎｃｅ ｈｅｉｇｈｔ ｓｙｓｔｅｍ ｗｅ ｃａｎ ｍａｋｅ ｔｈｅ ｃｏｍｐａｒｉｓｏｎ ｏｆ ｌｏｃａｌ ｐｒｅｃｉｓｅ ｇｅｏｉｄ ｗｉｔｈ ｔｈｅ

ｍｉｎｅ ｔｈｅ ｌｏｃａｌ ｐｒｅｃｉｓｅ ｇｅｏｉｄ ｍｏｄｅｌ ｆｏｒ ｔｈｅ ｃｏｎｓｉｓｔｅｎｔ ｈｅｉｇｈｔ ｓｙｓｔｅｍ ｏｆ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ｏｆ ｔｈｅ Ｇｅｏｇａ Ｇｒａｎｄ Ｂｒｉｄｇｅ ． Ｆｉｇ ． ２ ｓｈｏｗｓ ｔｈｅ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ ｂｅｎｃｈｍａｒｋｓ ａｎｄ ｃｏｎｔｒｏｌ ｐｏｉｎｔｓ ｉｎ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ｏｆ ｔｈｅ Ｇｅｏｇａ Ｇｒａｎｄ Ｂｒｉｄｇｅ ．

（）

Ｆｉｇ． ２ ． Ｌｏｃａｔｉｏｎ ｏｆ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ｌｅｆｔ ａｎｄ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ ｂｅｎｃｈｍａｒｋｓ

（ ）

ａｎｄ ｃｏｎｔｒｏｌ ｐｏｉｎｔｓ ｉｎ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ｒｉｇｈｔ ．

（

）

Ｔｈｉｓ ｐａｐｅｒ ｄｅｓｃｒｉｂｅｓ ｔｈａｔ ａ ｎｅｗ ｖｅｒｔｉｃａｌ ｄａｔｕｍ ｔｈｅ ｃｏｎｓｉｓｔｅｎｔ ｈｅｉｇｈｔ ｓｙｓｔｅｍ ｆｏｒ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ｗａｓ ｅｓｔａｂｌｉｓｈｅｄ ｂｙ ｕｓｅ ｏｆ ｐｒｅｃｉｓｅ ｌｏｃａｌ ｇｅｏｉｄ ｍｏｄｅｌ ｔｏ ｍｉｎｉｍｉｚｅ ｔｈｅ ｃｏｓｔ ａｎｄ ｔｉｍｅ ｃｏｎｓｕｍｐｔｉｏｎ ｐｒｏｂｌｅｍ ｄｕｅ ｔｏ ｔｈｅ ｈｅｉｇｈｔ ｅｒｒｏｒｓ ｆｒｏｍ ｔｈｅ ｔｒｉｇｏｎｏｍｅｔｒｉｃ ｌｅｖｅｌｉｎｇ ｍｅｔｈｏｄ ． Ｔｏ ｄｅｔｅｒｍｉｎｅ ａ ｐｒｅｃｉｓｅ ｌｏｃａｌ

，

ｇｅｏｉｄ ｆｏｒ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ｗｅ ｆｉｒｓｔｌｙ ｄｅｖｅｌｏｐｅｄ ａ ｐｒｅｃｉｓｅ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ｆｏｒ Ｋｏｒｅａ ａｎｄ ｉｔｓ ａｄ

，

ｊｏｉｎｉｎｇ ｓｅａｓ ａｓ ａ ｗｈｏｌｅ ． Ｔｈｉｓ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ｗａｓ ｄｅｖｅｌｏｐｅｄ ｂｙ ｕｓｅ ｏｆ ａｌｌ ｔｈｅ ａｖａｉｌａｂｌｅ ｇｒａｖｉｔｙ ｄａｔａ

ｉｎｃｌｕｄｉｎｇ ｓｕｒｆａｃｅ ａｎｄ ｓａｔｅｌｌｉｔｅ ｄａｔａ ｏｎ ｌａｎｄ ａｎｄ ｏｃｅａｎ ａｒｅａ ． Ｔｈｅ ｅｘｉｓｔｉｎｇ ｓｕｒｆａｃｅ ｇｒａｖｉｍｅｔｒｉｃ ｄａｔａ ａｒｅ

， ，

ｎｏｔ ｓｕｆｆｉｃｉｅｎｔｌｙ ｄｅｎｓｅ ｆｏｒ ｔｈｅ ｄｅｔｅｒｍｉｎａｔｉｏｎ ｏｆ ａ ｐｒｅｃｉｓｅ ｌｏｃａｌ ｇｅｏｉｄ ｔｈｕｓ ｗｅ ｐｅｒｆｏｒｍ ａ ｒｅｌａｔｉｖｅ ｇｒａｖｉ ｍｅｔｒｉｃ ｓｕｒｖｅｙ ａｔ ｓｉｘｔｙｏｎｅ ｐｏｉｎｔｓ ｗｉｔｈ ２ ｋｍ ｓｐａｃｅｓ ａｐｐｒｏｘｉｍａｔｅｌｙ ｉｎ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ． Ｔｈｅ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ｗａｓ ｃｏｍｐｕｔｅｄ ｂｙ ｓｐｈｅｒｉｃａｌ ｆａｓｔ ｆｏｕｒｉｅｒ ｔｒａｎｓｆｏｒｍ

（

， ）

（ＦＦＴ）ｗｉｔｈ ｍｏｄｉｆｉｅｄ

Ｓｔｏｋｅｓ′ ｋｅｒｎｅｌｓ Ｆｏｒｓｂｅｒｇ ｅｔ ａｌ ． １９９３ ． Ｔｈｅ ｒｅｍｏｖｅｒｅｓｔｏｒｅ ｔｅｃｈｎｉｑｕｅ ｗａｓ ｕｓｅｄ ｔｏ ｅｌｉｍｉｎａｔｅ ｔｈｅ ｔｅｒ

（ ） （ ， ） ｔｈｅ ｒｅｓｉｄｕａｌ ｇｅｏｉｄ ｂｙ ｃｏｍｂｉｎｉｎｇ ａ ｇｅｏｐｏｔｅｎｔｉａｌ ｍｏｄｅｌ，ｆｒｅｅａｉｒ ｇｒａｖｉｔｙ ａｎｏｍａｌｉｅｓ，ａｎｄ ｈｉｇｈｑｕａｌｉｔｙ ｒｅｓ

ｒａｉｎ ｅｆｆｅｃｔｓ ｏｎ ｌａｎｄ ｂｙ ｔｈｅ ｒｅｓｉｄｕａｌ ｔｅｒｒａｉｎ ｍｏｄｅｌ ＲＴＭ ｒｅｄｕｃｔｉｏｎ Ｆｏｒｓｂｅｒｇ １９８５ ａｎｄ ｔｏ ｄｅｔｅｒｍｉｎｅ ｏｌｕｔｉｏｎ ＤＥＭ ｄａｔａ ．

，

Ｆｏｒ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ｔｈｅ ａｃｃｕｒａｃｙ ｏｆ ｔｈｅ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ｃｏｍｐｕｔｅｄ ｆｒｏｍ ｔｈｅ ｄｉｆｆｅｒｅｎｃｅ ｂｅ

， （），２００９，６７９ － ６９４

ＤｏｎｇＨａ ＬＥＥ ｅｔ ａｌ ． ／ Ｃｈｉｎａ Ｏｃｅａｎ Ｅｎｇｉｎｅｅｒｉｎｇ ２３ ４

６８２

（ ） （ ） ｗａｓ ｄｅｔｅｒｍｉｎｅｄ ． Ｔｈｅ ｇｅｏｉｄ ｄｉｆｆｅｒｅｎｃｅ （Ｎ － Ｎ ）ｈａｄ ａ ｍｅａｎ ｖａｌｕｅ ｏｆ ０ ． ５８６ ｍ ａｎｄ ａ ｓｔａｎｄａｒｄ ｄｅｖｉａｔｉｏｎ ｏｆ ４ ． ８ ｃｍ；ｔｈｅｓｅ ｒｅｓｕｌｔｓ ｗｅｒｅ ｑｕｉｔｅ ｓａｔｉｓｆａｃｔｏｒｙ ｔｏ ｇｉｖｅ ｔｈｅ ｌａｒｇｅ ｅｘｐｅｃｔｅｄ ｃｈａｎｇｅｓ ｉｎ ｔｈｅ ｇｒａｖｉｔｙ ｆｉｅｌｄ，ｔｈｅ ｌａｃｋ ｏｆ ａｃｃｕｒａｔｅ ｇｒａｖｉｔｙ ｄａｔａ ｏｆｆｓｈｏｒｅ，ａｎｄ ｔｈｅ ｎｏｉｓｅ ｉｎ ｔｈｅ ｔｉｄａｌ ｈｅｉｇｈｔ ｔｒａｎｓｆｅｒ ． Ａｓ ｓｕｃｈ，ｔｈｉｓ ｇｒａｖｉｍｅｔｒｉｃ ｍｏｄｅｌ ｗａｓ ｓｕｂｓｅｑｕｅｎｔｌｙ ｆｉｔｔｅｄ ｔｏ ｇｅｏｉｄａｌ ｈｅｉｇｈｔｓ ｆｒｏｍ ＧＰＳ ／ ｌｅｖｅｌｉｎｇ ａｎｄ ｔｉｄｅ ｏｂｓｅｒｖａｔｉｏｎｓ ｂｙ ｕｓｅ ｏｆ ｌｅａｓｔ ｓｑｕａｒｅ ｃｏｌｌｏｃａｔｉｏｎ，ｔｏ ｐｒｏｖｉｄｅ ｔｈｅ ｆｉｎａｌ ＧＰＳｃｏｎｓｉｓｔｅｎｔ ｌｏｃａｌ ｐｒｅｃｉｓｅ Ｋｕｒｏｉｓｈｉ ｅｔ ａｌ ． ，２００２）． Ｔｈｅ ｐｏｓｔｆｉｔ ｅｒｒｏｒ （ｓｔｄ ． ｄｅｖ ． ）ｏｆ ｔｈｅ ｆｉｎａｌ ｇｅｏｉｄ ｔｏ ｔｈｅ Ｎ ｇｅｏｉｄ ｍｏｄｅｌ （ ｄｅｒｉｖｅｄ ｆｒｏｍ ＧＰＳ ａｎｄ ｔｉｄｅ ｏｂｓｅｒｖａｔｉｏｎｓ ｗａｓ ± ２ ． ２ ｃｍ，ａｎｄ ｔｈｅ ｆｉｎａｌ ｇｅｏｉｄ ｃｏｕｌｄ ｂｅ ｕｓｅｄ ｔｏ ｏｂｔａｉｎ ｃｏｎｓｉｓｔｅｎｔ ｈｅｉｇｈｔｓ ｆｏｒ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ． Ｔｈｅｒｅｆｏｒｅ，ｔｈｅ ｉｎｃｏｎｓｉｓｔｅｎｔ ｈｅｉｇｈｔ ｓｙｓｔｅｍ ｐｒｏｂｌｅｍ （ａｓ ａ ｌｉｎｋｉｎｇ ｐｒｏｂｌｅｍ）ｗａｓ ｓｏｌｖｅｄ ｂｙ ｃａｌｃｕｌａｔｉｎｇ ｔｈｅ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔ ｏｆ ｔｈｅ ｂｅｎｃｈｍａｒｋｓ ａｎｄ ｔｈｅ ｃｏｎｔｒｏｌ

ｔｗｅｅｎ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄａｌ ｈｅｉｇｈｔ Ｎ ｇｒａｖ ａｎｄ ｇｅｏｉｄａｌ ｈｅｉｇｈｔ ｆｒｏｍ ＧＰＳ ａｎｄ ｔｉｄｅ ｏｂｓｅｒｖａｔｉｏｎｓ Ｎ ＧＰＳ ／ Ｔｉｄｅ ＧＰＳ ／ Ｔｉｄｅ

ｇｒａｖ

ＧＰＳ ／ Ｔｉｄｅ

ｐｏｉｎｔｓ ｕｓｉｎｇ ｔｈｅ ｆｉｎａｌ ｇｅｏｉｄ ｍｏｄｅｌ ａｎｄ ｅｌｌｉｐｓｏｉｄａｌ ｈｅｉｇｈｔ ｆｒｏｍ ＧＰＳ ｏｂｓｅｒｖａｔｉｏｎｓ ．

２ ． Ｇｒａｖｉｍｅｔｒｉｃ Ｇｅｏｉｄ ｏｆ Ｋｏｒｅａ ２．１

Ｃｏｍｐｕｔａｔｉｏｎａｌ Ｆｏｒｍｕｌａｅ ｆｏｒ Ｇｒａｖｉｍｅｔｒｉｃ Ｇｅｏｉｄ Ｄｅｔｅｒｍｉｎａｔｉｏｎ

，

Ｉｎ ｔｈｅ ｒｅｍｏｖｅｒｅｓｔｏｒｅ ｔｅｃｈｎｉｑｕｅ ｔｈｅ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄａｌ Ｎ ｈｅｉｇｈｔ ｉｓ ｓｐｌｉｔ ｉｎｔｏ ｔｈｒｅｅ ｐａｒｔｓ Ｎ ＝ Ｎ ＧＭ ＋ Ｎ ＲＴＭ ＋ Ｎ ｒｅｓ

，

（１）

ｗｈｅｒｅ Ｎ ＧＭ ｉｓ ｔｈｅ ｌｏｗｆｒｅｑｕｅｎｃｙ ｐａｒｔ ｏｆ ｔｈｅ ｇｅｏｉｄ ｏｂｔａｉｎｅｄ ｆｒｏｍ ａ ｇｅｏｐｏｔｅｎｔｉａｌ ｍｏｄｅｌ ． Ｆｏｒ ｃｏｍｐｕｔａｔｉｏｎ

， （ ，２００５） ｍｏｄｅｌ ａｎｄ ｔｈｅ ＥＧＭ９６ （Ｌｅｍｏｉｎｅ ｅｔ ａｌ ． ，１９９６）ｍｏｄｅｌ ｗａｓ ｕｓｅｄ ． ＥＧＭ９６ ｗａｓ ｄｅｔｅｒｍｉｎｅｄ ｆｒｏｍ ａ ｃｏｍ ｂｉｎａｔｉｏｎ ｏｆ ｓａｔｅｌｌｉｔｅ ｔｒａｃｋｉｎｇ ｄａｔａ，ｓａｔｅｌｌｉｔｅ ａｌｔｉｍｅｔｒｙ ｉｎ ｔｈｅ ｏｃｅａｎｓ ａｎｄ ｍｅａｎ ｇｒａｖｉｔｙ ａｎｏｍａｌｉｅｓ ｏｎ ｌａｎｄ ． Ｔｈｅ ｎｅｗ ｓａｔｅｌｌｉｔｅ ＧＲＡＣＥ ＧＧＭ０２Ｓ ｆｉｅｌｄ ｗａｓ ｕｓｅｄ ｕｐ ｔｏ ９０ ｄｅｇｒｅｅｓ ｔｏ ｉｍｐｒｏｖｅ ｔｈｅ ＥＧＭ９６ ｍｏｄｅｌ，ｓｏ ｔｈａｔ ＧＧＭ０２Ｓ ｃｏｕｌｄ ｂｅ ｕｓｅｄ ｅｘｃｌｕｓｉｖｅｌｙ ｔｏ ９０ ｄｅｇｒｅｅｓ，ａｎｄ ＥＧＭ９６ ａｂｏｖｅ ９０ ｄｅｇｒｅｅｓ，ｗｉｔｈ ｔｈｅ ｌｉｎｅａｒ ｍｅｒｇｉｎｇ ｏｆ ｓｐｈｅｒｉｃａｌ ｈａｒｍｏｎｉｃ ｃｏｅｆｆｉｃｉｅｎｔｓ ． Ｇｒｉｄｓ ｏｆ ｇｒａｖｉｔｙ （ ｇ ）ａｎｄ ｇｅｏｉｄ （Ｎ ）ｆｒｏｍ ｔｈｅ Δ ＧＧＭ０２Ｓ ／ ＥＧＭ９６ ｍｏｄｅｌ ｗｅｒｅ ｃｏｍｐｕｔｅｄ ｆｏｒ ｔｈｅ ａｒｅａ ３２° ～ ４０°Ｎ，１２３° ～ １３２°Ｅ． Ｎ ，ｔｈｅ ｈｉｇｈｆｒｅｑｕｅｎｃｙ ｐａｒｔ ｏｆ ｔｈｅ ｇｅｏｉｄ，ｉｓ ｔｈｅ ｔｅｒｒａｉｎ ｅｆｆｅｃｔ ｏｎ ｔｈｅ ｇｅｏｉｄ ｇｅｎｅｒａｔｅｄ ｂｙ ｔｈｅ ｒｅｓｉｄｕａｌ ｔｅｒｒａｉｎ ｍｏｄｅｌ （ ＲＴＭ）ｍｅｔｈｏｄ （Ｆｏｒｓｂｅｒｇ ｅｔ ａｌ ． ，１９８１）． Ｔｈｅ ＲＴＭ ｔｅｒｒａｉｎ ｒｅｄｕｃｔｉｏｎ ｅｖａｌｕａｔｅｓ ｔｈｅ ｇｒａｖｉｔａｔｉｏｎａｌ ｅｆｆｅｃｔｓ ｏｆ ｔｈｅ ｍａｓｓ ａｎｏｍａｌｉｅｓ ｒｅｌａｔｉｖｅ ｔｏ ａ ｍｅａｎ ｅｌｅｖａｔｉｏｎ ｓｕｒｆａｃｅ （Ｆｏｒｓｂｅｒｇ，１９８５）． Ｔｈｅ ｍｅａｎ ｅｌｅｖａｔｉｏｎ ｓｕｒｆａｃｅ ｉｓ ｄｅｔｅｒｍｉｎｅｄ ｂｙ ｍｏｖｉｎｇａｖｅｒａｇｉｎｇ ｆｉｌｔｅｒｉｎｇ ｂｙ ｌｏｃａｌ ＤＥＭ ｄａｔａ （Ｔｃｈｅｒｎｉｎｇ ｅｔ ａｌ ． ，１９９４）． Ｔｈｅ ｍｅａｎ ｅｌｅｖａｔｉｏｎ ｓｕｒｆａｃｅ ｃｏｒｒｅｓｐｏｎｄｓ ｉｎ ｐｒｉｎｃｉｐｌｅ ｔｏ ｔｈｅ ｔｏｐｏｇｒａｐｈｉｃ ｓｉｇｎａｌ ａｌｒｅａｄｙ ｐｒｅｓｅｎｔ ｉｎ ｔｈｅ ｓｐｈｅｒｉｃａｌ ｈａｒｍｏｎｉｃ ｒｅｆｅｒｅｎｃｅ ｍｏｄｅｌ ． Ｉｎ ｔｈｅ ｃａｓｅ ｏｆ Ｋｏｒｅａ，ｔｈｅ ｍｅａｎ ＤＥＭ ｓｕｒｆａｃｅ ｏｆ ａｐ ｐｒｏｘ ． ７５ ｋｍ （ ４２′ × ５８′）ｒｅｓｏｌｕｔｉｏｎ ｗａｓ ｕｓｅｄ ． Ｆｏｒ ｔｈｅ ｃａｌｃｕｌａｔｉｏｎ ｏｆ ｔｈｅ ｔｅｒｒａｉｎ ｅｆｆｅｃｔ ｂｙ ＲＴＭ，ｗｅ ｕｓｅｄ ｌｏｃａｌ ＤＥＭ ｏｎ ａ １００ ｍ ｒｅｓｏｌｕｔｉｏｎ ｇｒｉｄ， ｄｅｒｉｖｅｄ ｆｒｏｍ ｏｒｉｇｉｎａｌ ＤＥＭ ｏｎ ａ ３０ ｍ ｇｒｉｄ （ｓｅｅ Ｆｉｇ ． ３）． Ｔｈｅ ｏｒｉｇｉｎａｌ ＤＥＭ ｗａｓ ｐｒｏｄｕｃｅｄ ｆｒｏｍ ｄｉｇｉｔａｌ ｔｏｐｏｇｒａｐｈｉｃ ｍａｐｓ ｄｒａｗｎ ｏｎ ｓｃａｌｅ １： ２５０００ ａｎｄ ｐｒｏｖｉｄｅｄ ｄｉｒｅｃｔｌｙ ｂｙ Ｅｎｖｉｒｏｎｍｅｎｔａｌ Ｇｅｏｇｒａｐｈｉｃ Ｉｎｆｏｒｍａ ，Ｋｏｒｅａ （ｓｅｅ，ｈｔｔｐ：／ ／ ｅｇｉｓ ． ｍｅ ． ｇｏ ． ｋｒ ／ ）． Ｔｈｅ ｈｅｉｇｈｔ ｖａｌｕｅ ｉｎ ｌｏｃａｌ ＤＥＭ ｄｉｓｔｒｉｂｕｔ ｔｉｏｎ Ｓｙｓｔｅｍ （ＥＧＩＳ） ｏｆ ｔｈｅ ｒｅｆｅｒｅｎｃｅ ｇｅｏｉｄ ｕｎｄｕｌａｔｉｏｎ Ｎ ＧＭ ｔｈｅ ｃｏｍｂｉｎｅｄ ｍｏｄｅｌ ｏｆ ｔｈｅ ＧＧＭ０２Ｓ Ｔａｐｌｅｙ ｅｔ ａｌ ．

ＧＭ

ＧＭ

ＲＴＭ

ｅｄ ｆｒｏｍ ０ ｔｏ １９５０ ｍ ａｎｄ ｔｈｅ ｓｔａｎｄａｒｄ ｄｅｖｉａｔｉｏｎ ｗａｓ ａｂｏｕｔ ± １２ ． ０ ｍ ｂｙ ｃｏｍｐａｒｉｎｇ ｗｉｔｈ ｔｈｅ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔｓ ｏｆ ＧＰＳ ／ ｌｅｖｅｌｉｎｇ ｄａｔａ ．

， （），２００９，６７９ － ６９４

ＤｏｎｇＨａ ＬＥＥ ｅｔ ａｌ ． ／ Ｃｈｉｎａ Ｏｃｅａｎ Ｅｎｇｉｎｅｅｒｉｎｇ ２３ ４

６８３

（

， ）：

Ｔｈｅ ｅｘｐｒｅｓｓｉｏｎ ｆｏｒ ｔｈｅ ｇｅｏｉｄ ｅｆｆｅｃｔ ｏｆ ｔｈｅ ｔｏｐｏｇｒａｐｈｙ ｉｓ ｇｉｖｅｎ ａｓ Ｏｍａｎｇ ｅｔ ａｌ ． ２０００

（ ，） ｚ ＝ ｈ （ｘ ，ｙ ） （ 槡 ｘＱ

ρ ∫∫∫ γ

） （２） ｗｈｅｒｅ Ｇ ｉｓ ｔｈｅ ｇｒａｖｉｔａｔｉｏｎａｌ ｃｏｎｓｔａｎｔ ａｎｄ ρ ｉｓ ｔｈｅ ｄｅｎｓｉｔｙ ｃｏｎｓｔａｎｔ ｏｆ ｔｏｐｏｇｒａｐｈｉｃ ｍａｓｓ ． （ｘ ，ｙ ，ｚ ） ｉｓ ｔｈｅ ｃｏｏｒｄｉｎａｔｅ ｏｆ ＤＥＭ ｇｒｉｄ ｐｏｉｎｔ （Ｑ ）ａｎｄ （ｘ ，ｘ ，ｘ ）ｉｓ ｔｈｅ ｃｏｏｒｄｉｎａｔｅ ｏｆ ｇｒａｖｉｔｙ ｓｔａｔｉｏｎ （Ｐ ）． Ｉｎ ｐｒａｃｔｉｃｅ，Ｎ ｃａｎ ｂｅ ｃｏｍｐｕｔｅｄ ｉｎ ｔｈｅ ｆｒｅｑｕｅｎｃｙ ｄｏｍａｉｎ ｂｙ Ｆｏｕｒｉｅｒ ｍｅｔｈｏｄｓ，ｕｓｉｎｇ ａ ｌｉｎｅａｒ ｔｅｒｍ （ｃｏｒｒｅｓｐｏｎｄｉｎｇ ｔｏ ｔｈｅ ｒｅｓｉｄｕａｌ ｔｏｐｏｇｒａｐｈｙ ｂｅｉｎｇ ｃｏｎｄｅｎｓｅｄ ｏｎ ｔｈｅ ｇｅｏｉｄ）（Ｆｏｒｓｂｅｒｇ，１９８５） Ｇｐ １ （３） Ｎ ＝ γ（ｈ － ｈ ） ｓ ， ｗｈｅｒｅ ｓ ｉｓ ｔｈｅ ｈｏｒｉｚｏｎｔａｌ ｄｉｓｔａｎｃｅ ｂｅｔｗｅｅｎ ｔｈｅ ｇｒａｖｉｔｙ ｓｔａｔｉｏｎ ａｎｄ ｔｈｅ ＤＥＭ ｇｒｉｄ ｐｏｉｎｔ，ａｎｄ  ｄｅｎｏｔｅｓ ｔｈｅ ２Ｄ ｃｏｎｖｏｌｕｔｉｏｎ，ｆｏｒ ｄｅｔａｉｌｓ ｉｎ Ｆｏｒｓｂｅｒｇ （ １９８５）． Ｔｈｅ ｃｏｒｒｅｓｐｏｎｄｉｎｇ ＲＴＭ ｅｆｆｅｃｔｓ ｏｎ ｇｒａｖｉｔｙ ｄａｔａ （Δｇ ）ｗｅｒｅ ｃｏｍｐｕｔｅｄ ｄｉｒｅｃｔｌｙ ｂｙ ｓｐａｃｅ ｄｏｍａｉｎ ｐｒｉｓｍ ｉｎｔｅｇｒａｔｉｏｎ ｕｓｉｎｇ ｔｈｅ ｄｅｎｓｅ ｈｅｉｇｈｔ ｄａｔａ ． Ｎ ＲＴＭ ＝

Ｇ

∞

∞

－∞

－∞

ｄｘＱ ｄｙＱ ｄ ｚＱ

ｚ＝ｈ ｘ ｙ ｒｅｆ

－ ｘＰ

） ＋ （ｙ ２

Ｑ

－ ｙＰ

） ＋ （ｚ ２

Ｑ

－ ｚＰ

２

Ｑ

Ｐ

Ｐ

Ｑ

Ｑ

Ｐ

ＲＴＭ

ＲＴＭ

ｒｅｆ

ＲＴＭ

Ｆｉｇ． ３ ． Ｌｏｃａｌ ＤＥＭ ｄａｔａ ｗｉｔｈ １００ ｍ ｒｅｓｏｌｕ

（ ：）

ｔｉｏｎ ｕｎｉｔ ｍ ．

，

，

Ｎ ｒｅｓ ｔｈｅ ｍｅｄｉｕｍｆｒｅｑｕｅｎｃｙ ｐａｒｔ ｏｆ ｔｈｅ ｇｅｏｉｄ ｉｓ ｔｈｅ ｒｅｓｉｄｕａｌ ｇｅｏｉｄ ｃｏｍｐｕｔｅｄ ｆｒｏｍ

，

， Ｒ ， ｇ Ｓ（ ＝ Δ Ψ）ｄσ  ４π γ σ

ｇ Δ

ｒｅｓ

ｂｙ ｕｓｅ

ｏｆ Ｓｔｏｋｅｓ ｉｎｔｅｇｒａｔｉｏｎ ｅｘｔｅｎｄｉｎｇ ｉｎ ｐｒｉｎｃｉｐｌｅ ｔｈｅ ｉｎｔｅｇｒａｔｉｏｎ ｔｏ ｃｏｖｅｒ ｔｈｅ Ｅａｒｔｈ ． Ｎ ｒｅｓ ｗｈｅｒｅ

ｇ Δ

ｒｅｓ

（４）

ｒｅｓ

ｉｓ ｔｈｅ ｒｅｓｉｄｕａｌ ｆｒｅｅａｉｒ ａｎｏｍａｌｙ ｔｈａｔ ｒｅｍａｉｎｓ ｉｎ ｔｈｅ ｇｒａｖｉｔｙ ｄａｔａ ａｆｔｅｒ ｔｈｅ ｃｏｎｔｒｉｂｕｔｉｏｎｓ ｏｆ

ｔｈｅ ｒｅｓｉｄｕａｌ ｔｅｒｒａｉｎ ｅｆｆｅｃｔ

ｇ Δ

ＲＴＭ

ａｎｄ ｔｈｅ ｇｌｏｂａｌ ｆｉｅｌｄ

ｇ Δ

ＧＭ

ａｒｅ ｓｕｂｔｒａｃｔｅｄ ． Ｔｈｅ ｉｎｄｉｒｅｃｔ ｅｆｆｅｃｔ ｏｎ ｇｒａｖｉ

（ ） ｉｔｓ ｅｆｆｅｃｔ ｉｓ ｏｎｌｙ ａ ｆｅｗ ｃｅｎｔｉｍｅｔｒｅｓ ｏｎ ｔｈｅ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ｃｏｍｐｕｔａｔｉｏｎ （Ｃｏｒｃｈｅｔｅ ｅｔ ａｌ ． ，２００５）． Ｔｈｕｓ，ｉｎ ｔｈｉｓ ｓｔｕｄｙ，ｔｈｅ ｓｅｃｏｎｄａｒｙ ｉｎｄｉｒｅｃｔ ｅｆｆｅｃｔ ｈａｓ ｂｅｅｎ ｎｅｇｌｅｃｔｅｄ ． Ｔｈｅ ｆｕｎｃｔｉｏｎ Ｓ （ Ψ）ｉｓ Ｓｔｏｋｅｓ ｆｕｎｃｔｉｏｎ，Ψ ｂｅｉｎｇ ｔｈｅ ｓｐｈｅｒｉｃａｌ ｄｉｓｔａｎｃｅ ． ）． （５） Ｓ（ ｌｏｇ（ｓｉｎ Ψ＋ ｓｉｎ Ψ ＋ １ － ５ｃｏｓΨ － ３ｃｏｓΨ Ψ）＝ ｓｉｎ（１Ψ）－ ６ｓｉｎ Ψ ２ ２ ２ ２ Ｔｏ ｐｒｅｖｅｎｔ ｔｈｅ ｉｎｆｌｕｅｎｃｅ ｏｆ ｌｏｃａｌ ｇｒａｖｉｔｙ ｏｎ ｔｈｅ ｌｏｎｇｅｓｔ ｗａｖｅｌｅｎｇｔｈｓ （Ｖａｎíˇｃｅｋ ａｎｄ Ｆｅａｔｈｅｒｓｔｏｎｅ， ） １９９８ ，ｔｈｅ ｍｏｄｉｆｉｅｄ Ｓｔｏｋｅｓ ｋｅｒｎｅｌ Ｓ （ Ψ）ｉｓ ｕｓｅｄ ． Ｔｈｅ ｔｅｃｈｎｉｑｕｅ ｉｓ ｔｈｅ ｍｏｄｉｆｉｅｄ ＷｏｎｇＧｏｒｅ ｔｙ ｏｒ ｓｅｃｏｎｄａｒｙ ｉｎｄｉｒｅｃｔ ｅｆｆｅｃｔ ｉｓ ｎｏｔ ｃｏｎｓｉｄｅｒｅｄ ｉｎ ｍｏｓｔ ｏｆ ｔｈｅ ｃａｓｅｓ ｂｅｃａｕｓｅ ｉｔ ｈａｓ ｓｍａｌｌ ｖａｌｕｅｓ ａｎｄ

２

ｍｏｄ

， （），２００９，６７９ － ６９４

ＤｏｎｇＨａ ＬＥＥ ｅｔ ａｌ ． ／ Ｃｈｉｎａ Ｏｃｅａｎ Ｅｎｇｉｎｅｅｒｉｎｇ ２３ ４

６８４

（

， ），ｉｎ ｗｈｉｃｈ ｔｈｅ ｍｏｄｉｆｉｅｄ ｋｅｒｎｅｌ Ｓ （Ψ）ｈａｓ ｔｈｅ ｆｏｌｌｏｗｉｎｇ ｅｘｐｒｅｓ

ｍｅｔｈｏｄ Ｗｏｎｇ ａｎｄ Ｇｏｒｅ １９６９

ｍｏｄ

ｓｉｏｎ ． Ｎ２

（Ψ）＝ Ｓ（Ψ）－ ∑α （ｎ）２ｎｎ －＋ １１ Ｐ ｃｏｓ（Ψ）， （ｎ）ｉｎｃｒｅａｓｅ ｌｉｎｅａｒｌｙ ｆｒｏｍ ０ ｔｏ １ ｂｅｔｗｅｅｎ ｄｅｇｒｅｅｓ Ｎ ａｎｄ Ｎ ｗｈｅｒｅ ｔｈｅ ｃｏｅｆｆｉｃｉｅｎｔｓ α Ｓ ｍｏｄ

（６）

ｎ

ｎ＝２

１

１

２

ｆｏｒ ２  ｎ  Ｎ １

  Ｎ ２ － ｎ ＝  Ｎ２ － Ｎ１   ０

（ｎ） α

，… ，Ｎ

ｆｏｒ

Ｎ１  ｎ  Ｎ２ ． ｎ ＝ ２

ｆｏｒ

Ｎ２  ｎ

（７） ，

Ｔｈｅ ｃｏｅｆｆｉｃｉｅｎｔｓ Ｎ １ ａｎｄ Ｎ ２ ａｒｅ ｄｅｔｅｒｍｉｎｅｄ ｂｙ ｔｒｉａｌ ａｎｄ ｅｒｒｏｒ ｏｒ ｂｙ ｅｘｐｅｒｉｅｎｃｅ ａｎｄ ｒｅｐｒｅｓｅｎｔ ａ ｔｒａｄｅｏｆｆ ｂｅｔｗｅｅｎ ｆｕｌｌ ｕｓｅ ｏｆ ｔｈｅ ｓａｔｅｌｌｉｔｅ ｆｉｅｌｄｓ ｆｏｒ ｎ ＞ Ｎ １ ａｎｄ ｆｕｌｌ ｕｓｅ ｏｆ ｔｈｅ ｌｏｃａｌ ｇｒａｖｉｔｙ ｄａｔａ ｎ ＞

（

）

）

（

Ｎ ２ ． Ｔｈｅ ｔａｓｋ ｏｆ ｏｐｔｉｍｉｚａｔｉｏｎ ｉｓ ｔｙｐｉｃａｌｌｙ ｔｏ ｆｉｎｄ ａ ｓｅｔ ｏｆ Ｎ １ ａｎｄ Ｎ ２ ｆｏｒ ｗｈｉｃｈ ｔｈｅ ｆｉｔ ｂｅｔｗｅｅｎ ｔｈｅ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ａｎｄ ｔｈｅ ｇｅｏｉｄａｌ ｈｅｉｇｈｔｓ ｄｅｒｉｖｅｄ ｆｒｏｍ ＧＰＳ ／ ｌｅｖｅｌｉｎｇ ｄａｔａ ｉｓ ｔｈｅ ｂｅｓｔ ｉｎ ｔｈｅ ｌｅａｓｔ ｓｑｕａｒｅｓ

（

， ）

，

ｓｅｎｓｅ Ｉｌｉｆｆｅ ｅｔ ａｌ ． ２００３ ． Ｂｙ ｔｈｅｓｅ ｏｐｔｉｍｉｚａｔｉｏｎ ｔｈｅ ｖａｌｕｅｓ ｏｆ Ｎ １ ＝ ８０ ａｎｄ Ｎ １ ＝ ９０ ａｒｅ ｄｅｔｅｒｍｉｎｅｄ ｆｏｒ ｔｈｅ ｆｉｎａｌ ｃｏｍｐｕｔａｔｉｏｎ ｏｆ ｒｅｓｉｄｕａｌ ｇｅｏｉｄ ．

（

Ｔｈｅ ｒｅｓｉｄｕａｌ ｇｒａｖｉｔｙ ｉｓ ｔｒａｎｓｆｏｒｍｅｄ ｉｎｔｏ ｔｈｅ ｒｅｓｉｄｕａｌ ｇｅｏｉｄ ｂｙ ｍｕｌｔｉｂａｎｄ ｓｐｈｅｒｉｃａｌ ＦＦＴ Ｆｏｒｓｂｅｒｇ

， ），ｗｈｉｃｈ ｐｒｏｖｉｄｅｓ ａ ｖｉｒｔｕａｌｌｙ ｅｘａｃｔ ｉｍｐｌｅｍｅｎｔａｔｉｏｎ ｏｆ Ｓｔｏｋｅｓ ｆｏｒｍｕｌａ ｏｎ ａ ｓｐｈｅｒｅ ． Ｔｈｅ ｇｅｏｉｄ ｉｓ ｏｂｔａｉｎｅｄ ｂｙ ａ ｎｕｍｂｅｒ ｏｆ ｂａｎｄｗｉｓｅ ｏｐｅｒａｔｉｏｎｓ ｏｆ ｔｈｅ ｆｏｒｍ ，Δ ）［Δｇ （φ，λ ）ｓｉｎφ］＝ Ｆ ［Ｆ（Ｓ）Ｆ（Δｇ ｓｉｎφ）］ （８） Ｎ ＝ Ｓ（ Δ λ φ ｅｔ ａｌ ． １９９３

ｒｅｓ

－１

ｒｅｓ

ｒｅｓ

ｗｈｅｒｅ Ｆ ｉｓ ｔｈｅ ｔｗｏｄｉｍｅｎｓｉｏｎａｌ ｆｏｕｒｉｅｒ ｔｒａｎｓｆｏｒｍ ｏｐｅｒａｔｏｒ ． Ｆ

（Δｇ）＝Δｇ（ｘ ，ｙ）ｅ （

－ ｉ ｋｘｘ ＋ ｋｙｙ

）ｄｘｄｙ ，

（９）

ｗｈｅｒｅ ｋｘ ａｎｄ ｋｙ ａｒｅ ｔｈｅ ｗａｖｅ ｎｕｍｂｅｒｓ ａｎｄ １００％ ｚｅｒｏ ｐａｄｄｉｎｇ ｉｓ ｕｓｅｄ ｉｎ ｔｈｉｓ ｓｔｕｄｙ ． Ｔｈｅ ｇｅｎｅｒａｌ ｄｅ

（ ） Ｔｈｅ ｏｕｔｃｏｍｅ ｏｆ ｔｈｅ ｒｅｍｏｖｅｒｅｓｔｏｒｅ ｔｅｃｈｎｉｑｕｅ ｉｓ ａ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ，ｒｅｆｅｒｒｉｎｇ ｔｏ ａ ｇｌｏｂａｌ ｈｅｉｇｈｔ ｄａｔｕｍ ａｐｐｒｏｘｉｍａｔｅｌｙ，ｗｈｉｃｈ ｆｉｔｓ ｔｈｅ Ｋｏｒｅａｎ ｖｅｒｔｉｃａｌ ｄａｔｕｍ，ａｎｄ ｍｉｎｉｍｉｚｅｓ ｔｈｅ ｐｏｓｓｉｂｌｅ ｓｙｓｔｅｍａｔｉｃ ｅｒ ｒｏｒｓ ｉｎ ｌｏｎｇ ｗａｖｅｌｅｎｇｔｈ （ｌｏｗｆｒｅｑｕｅｎｃｙ）ｐａｒｔ ｏｆ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ｄｕｅ ｔｏ ｅｒｒｏｒｓ ｏｆ ｇｅｏｐｏｔｅｎｔｉａｌ ｍｏｄｅｌ ａｎｄ ／ ｏｒ ｔｒｕｎｃａｔｉｏｎ ｐｒｏｃｅｄｕｒｅｓ （ Ｋｕｒｏｉｓｈｉ ｅｔ ａｌ ． ，２００２）． Ｔｈｉｓ ｇｅｏｉｄ ｍｕｓｔ ｂｅ ｆｉｔｔｅｄ ｔｏ ｔｈｅ ＧＰＳ ｃｏｎｔｒｏｌ ｉｎ ，ｄｅ ｔｈｅ ｆｉｎａｌ ｇｅｏｉｄ ｄｅｔｅｒｍｉｎａｔｉｏｎ ｓｔｅｐ ． Ｔｈｅ ｓｏｆｔｗａｒｅ ｐａｃｋａｇｅ ＧＲＡＶＳＯＦＴ （Ｔｓｃｈｅｒｎｉｎｇ ｅｔ ａｌ ． ，１９９４） ｖｅｌｏｐｅｄ ａｔ ｔｈｅ Ｄａｎｉｓｈ Ｎａｔｉｏｎａｌ Ｓｐａｃｅ Ｃｅｎｔｅｒ ａｎｄ Ｕｎｉｖｅｒｓｉｔｙ ｏｆ Ｃｏｐｅｎｈａｇｅｎ，ｉｓ ｕｓｅｄ ｆｏｒ ａｌｌ ｃｏｍｐｕｔａ ｔｉｏｎｓ ｉｎ ｔｈｉｓ ｓｔｕｄｙ （Ｆｏｒｓｂｅｒｇ ｅｔ ａｌ ． ，２００３）． ｔａｉｌｓ ｏｆ ｔｈｅ ＦＦＴ ｍｅｔｈｏｄｓ ｗｅｒｅ ｇｉｖｅｎ ｂｙ Ｓｃｈｗａｒｚ ｅｔ ａｌ ． １９９０ ．

２．２

Ｇｒａｖｉｔｙ Ｄａｔａ ｉｎ Ｋｏｒｅａ

Ｓｅｖｅｒａｌ ｇｒａｖｉｔｙ ｄａｔａ ｓｏｕｒｃｅｓ ｆｏｒ ｏｎｓｈｏｒｅ Ｋｏｒｅａ ｆｒｏｍ ｖａｒｉｏｕｓ ｎａｔｉｏｎａｌ ａｎｄ ｉｎｔｅｒｎａｔｉｏｎａｌ ａｇｅｎｃｉｅｓ

，

ｗｅｒｅ ｅｖａｌｕａｔｅｄ ｉｎ ｐｒｅｐａｒａｔｉｏｎ ｆｏｒ ｔｈｅ ｇｅｏｉｄ ｄｅｔｅｒｍｉｎａｔｉｏｎ ． Ｆｒｏｍ ＮＧＩＩ ２８５７ ｐｏｉｎｔｓ ｇｒａｖｉｔｙ ｄａｔａ ｗｅｒｅ

（ ）

ｏｂｔａｉｎｅｄ ａｎｄ ２０００ ｐｏｉｎｔｓ ｇｒａｖｉｔｙ ｄａｔａ ｆｒｏｍ Ｂｕｒｅａｕ Ｇｒａｖｉｍｅｔｒｉｑｕｅ Ｉｎｔｅｒｎａｔｉｏｎａｌ ＢＧＩ ｗｅｒｅ ｕｓｅｄ ． Ｔｈｅ

，

，

ｏｃｅａｎ ｇｒａｖｉｔｙ ｄａｔａ ｄｅｒｉｖｅｄ ｆｒｏｍ ＥＲＳ ＧＥＯＳＡＴ ａｎｄ ＴＯＰＥＸ ｓａｔｅｌｌｉｔｅ ａｌｔｉｍｅｔｒｙ ｏｆ ｔｈｅ ｏｃｅａｎ ａｒｅａｓ ｆｒｏｍ

（

， ），ｗｅｒｅ ｕｓｅｄ ． Ｔｈｅｒｅ ｗｅｒｅ ｕｎｃｅｒｔａｉｎｔｉｅｓ ｉｎ ｔｈｅ ｄａｔｕｍ ａｎｄ ｓｙｓｔｅｍｓ ｏｆ ｔｈｅ ｒｅｆｅｒｅｎｃｅ ｆｏｒ ｓｏｍｅ ｏｆ ｔｈｅｓｅ ｇｒａｖｉｔｙ ｄａｔａ，ａｎｄ ｅｒｒｏｒｓ ｗｅｒｅ ｃｏｎｆｉｒｍｅｄ ｂｙ ｂｉａｓ ｃｏｍｐａｒ

ＫＭＳ２００２ ｇｌｏｂａｌ ｓｏｌｕｔｉｏｎ Ａｎｄｅｒｓｅｎ ｅｔ ａｌ ． ２００５

， （），２００９，６７９ － ６９４

ＤｏｎｇＨａ ＬＥＥ ｅｔ ａｌ ． ／ Ｃｈｉｎａ Ｏｃｅａｎ Ｅｎｇｉｎｅｅｒｉｎｇ ２３ ４

６８５

ｉｓｏｎｓ ｗｉｔｈ ｔｈｅ ＧＧＭ０２Ｓ ／ ＥＧＭ９６ ｒｅｆｅｒｅｎｃｅ ｍｏｄｅｌ ．

，

， ｗａｓ ｃａｒｒｉｅｄ ｏｕｔ ． Ｔｈｉｓ ｓｕｒｖｅｙ ｕｓｅｄ ａ Ｌａｃｏｓｔｅ ａｎｄ Ｒｏｍｂｅｒｇ Ｇ２００ ｇｒａｖｉｍｅｔｅｒ，ｔｙｉｎｇ ｔｈｅ ｒｅｌａｔｉｖｅ ｇｒａｖｉｔｙ ｍｅａｓｕｒｅｍｅｎｔｓ ｉｎｔｏ ａｎ ａｂｓｏｌｕｔｅ ｓｔａｔｉｏｎ ＤＡＥＪ ｗｉｔｈ ｔｈｅ ｒｅｆｅｒｅｎｃｅ ｇｒａｖｉｔｙ ｖａｌｕｅ，ｇ ：９７９８３２ ． ４４４ ｍｇａｌ ． Ｔｈｅ ｇｒａｖｉｍｅｔｒｉｃ ｓｕｒｖｅｙ ｔｏｏｋ ｐｌａｃｅ ｄｕｒｉｎｇ ｔｈｅ ｐｅｒｉｏｄ １２ ～ １６ Ｓｅｐ ． ２００５，ｏｂｓｅｒｖｉｎｇ ａ ｔｏｔａｌ ｏｆ ７０ ｏｂｓｅｒｖａｔｉｏｎｓ ａｔ ６１ ｐｏｉｎｔｓ ． Ａ ｌｅａｓｔｓｑｕａｒｅ ａｄｊｕｓｔｍｅｎｔ ｗａｓ ｃａｒｒｉｅｄ ｏｕｔ ｏｎ ｔｈｅ ｏｂｓｅｒｖａｔｉｏｎｓ，ｇｉｖｉｎｇ ａ ｓｔａｎ ｄａｒｄ ｅｒｒｏｒ ｏｆ ０ ． ０７３ ｍｇａｌ ｆｏｒ ａ ｓｉｎｇｌｅ ｏｂｓｅｒｖａｔｉｏｎ ． Ａ ｓｉｎｇｌｅ ｉｎｓｔｒｕｍｅｎｔ ｔａｒｅ ｏｆ ０ ． ２４ ｍｇａｌ ｗａｓ ｄｅｔｅｃｔｅｄ ｏｎ １３ Ｓｅｐ ． ｂｅｔｗｅｅｎ ｔｗｏ ｐｏｉｎｔｓ，ｓｏ ｔｈｅ ｐｏｉｎｔｓ ｉｎ ｔｈｉｓ ｓｅｑｕｅｎｃｅ ｍｉｇｈｔ ｈａｖｅ ａｎ ｅｒｒｏｒ ｏｆ ｔｈｉｓ ｍａｇｎｉｔｕｄｅ ． Ｈｏｗｅｖｅｒ，ｔｈｉｓ ｍａｇｎｉｔｕｄｅ ｏｆ ｅｒｒｏｒ ｈａｓ ｎｅｇｌｉｇｉｂｌｅ ｉｎｆｌｕｅｎｃｅ ｏｎ ｔｈｅ ｇｅｏｉｄ ｒｅｓｕｌｔｓ，ａｎｄ ａｌｌ ｏｂｓｅｒｖａｔｉｏｎｓ ａｒｅ Ｈｏｗｅｖｅｒ ｔｈｅ ｅｘｉｓｔｉｎｇ ｓｕｒｆａｃｅ ｇｒａｖｉｍｅｔｒｉｃ ｄａｔａ ｗｅｒｅ ｎｏｔ ｓｕｆｆｉｃｉｅｎｔｌｙ ｄｅｎｓｅ ｓｏ ａ ｇｒａｖｉｍｅｔｒｉｃ ｓｕｒｖｅｙ

ｕｓｅｄ ｆｏｒ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ｃｏｍｐｕｔａｔｉｏｎ ．

Ｆｉｇ． ４ ． Ｌｏｃａｔｉｏｎ ｏｆ ｎｅｗ ｇｒａｖｉｔｙ ｐｏｉｎｔｓ ａｒｏｕｎｄ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ． Ａｌｌ ｐｏｉｎｔｓ ｗｅｒｅ ｓｕｒｖｅｙｅｄ ｂｙ ｐｒｅｃｉｓｅ ＤＧＰＳ ｒｅｌａｔｉｖｅ ｔｏ ｔｈｅ ｆｕｎｄａｍｅｎｔａｌ

（ ）

ｂｅｎｃｈｍａｒｋ ａｔ Ｐｕｓａｎ ＪＩＮＨ ．

（λ ， ）ａｎｄ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔｓ （Ｈ ）ｏｆ ｇｒａｖｉｔｙ ｐｏｉｎｔｓ ａｐｐｒｏｘｉｍａｔｅｌｙ，ｗｅ ｐｅｒｆｏｒｍｅｄ ａ ｄｉｆｆｅｒｅｎｔｉａｌ ＧＰＳ （ ＤＧＰＳ）ｓｕｒｖｅｙ ｒｅｌａｔｉｖｅ ｔｏ ｔｈｅ ｆｕｎｄａｍｅｎｔａｌ ｂｅｎｃｈｍａｒｋ ａｔ Ｐｕｓａｎ （ＪＩＮＨ）ｗｉｔｈ ｇｒａｖｉｍｅｔｒｉｃ ｓｕｒｖｅｙ ｓｉｍｕｌｔａｎｅｏｕｓｌｙ ． ＧＰＳ ｏｂｓｅｒｖａｔｉｏｎ ｐｅｒｉｏｄ ｆｏｒ ｄａｔａ ｓｅｔｓ ａｒｅ ｂｅ ｔｗｅｅｎ １ ～ ２ ｈｏｕｒｓ ． Ｔｒｉｍｂｌｅ Ｇｅｏｍａｔｉｃｓ Ｏｆｆｉｃｅ （ＴＧＯ）ｓｏｆｔｗａｒｅ ｗａｓ ｕｓｅｄ ｆｏｒ ｂａｓｅｌｉｎｅ ｐｒｏｃｅｓｓｉｎｇ ｏｆ ｔｈｅ ｗｈｏｌｅ ＧＰＳ ｄａｔａ ａｎｄ ｎｅｔｗｏｒｋ ａｄｊｕｓｔｍｅｎｔ ． Ｔｈｅ ｎｅｔｗｏｒｋ ａｄｊｕｓｔｍｅｎｔ ｉｓ ｂａｓｅｄ ｏｎ ａ ｌｏｃａｌ ｄａｔｕｍ ａｓ ｔｈｅ ｆｉｘｅｄ ｈｅｉｇｈｔ ｏｆ ＪＩＮＨ ａｎｄ ｇｅｏｇｒａｐｈｉｃａｌ ｃｏｏｒｄｉｎａｔｅｓ ｏｆ ｔｗｏ ＧＰＳ ｐｅｒｍａｎｅｎｔ ｓｔａｔｉｏｎｓ （ＪＩＮＪ，ＰＵＳＮ）ｌｏｃａｔ ｅｄ ｎｅａｒ ｔｈｅ ｓｔｕｄｙ ａｒｅａ ． Ａｃｃｏｒｄｉｎｇ ｔｏ ｔｈｅ ａｄｊｕｓｔｍｅｎｔ ｒｅｓｕｌｔｓ，ｔｈｅ ｈｏｒｉｚｏｎｔａｌ （ｏｒ ｇｅｏｇｒａｐｈｉｃａｌ）ａｎｄ ｖｅｒｔｉ ｃａｌ （ｏｒ ｏｒｔｈｏｍｅｔｒｉｃ）ｃｏｏｒｄｉｎａｔｅｓ ｏｆ ｔｈｅ ｇｒａｖｉｔｙ ｐｏｉｎｔｓ ａｒｅ ｄｅｔｅｒｍｉｎｅｄ ｗｉｔｈ ｔｈｅ ｐｒｅｃｉｓｉｏｎ ｓｍａｌｌｅｒ ｔｈａｎ １ ｃｍ ａｎｄ ５ ｃｍ，ｒｅｓｐｅｃｔｉｖｅｌｙ ． Ｔｈｅ ｇｒａｖｉｔｙ ｄａｔａ ｗａｓ ｅｖｉｄｅｎｔｌｙ ａｆｆｅｃｔｅｄ ｂｙ ｓｏｍｅ ｇｒｏｓｓ ｅｒｒｏｒｓ ． Ｏｎｌｙ ａ ｆｅｗ ｌａｒｇｅ ｏｕｔｌｉｅｒ ｖａｌｕｅｓ ｗｅｒｅ ｒｅｍｏｖｅｄ ｗｉｔｈ ９５％ ｓｉｇｎｉｆｉｃａｎｃｅ ｌｅｖｅｌ，ｃｏｒｒｅｓｐｏｎｄｉｎｇ ｔｏ ｔｈｅ ａｒｅａｓ ｍａｉｎｌｙ ｃｌｏｓｅ ｔｏ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ Ｔｏ ｄｅｔｅｒｍｉｎｅ ｔｈｅ ｇｅｏｇｒａｐｈｉｃａｌ ｃｏｏｒｄｉｎａｔｅｓ

ａｎｄ ｉｎ ｔｈｅ ＢＧＩ ａｎｄ ＫＭＳ ｄａｔａ ． Ｔｈｅ ｇｒａｖｉｔｙ ｄａｔａ ｗｅｒｅ ｒｅｄｕｃｅｄ ｆｏｒ ｔｈｅ ｒｅｆｅｒｅｎｃｅ ｆｉｅｌｄ ａｎｄ ＲＴＭ ｔｅｒｒａｉｎ ｅｆｆｅｃｔｓ ． Ｔｅｒｒａｉｎ ｅｆｆｅｃｔｓ ｗｅｒｅ ｃｏｍｐｕｔｅｄ ｆｒｏｍ ｔｈｅ １００ ｍ ｒｅｓｏｌｕｔｉｏｎ ｈｅｉｇｈｔ ｄａｔａ ｕｓｉｎｇ ｐｒｉｓｍ ｉｎｔｅｇｒａｔｉｏｎ ． Ｔｈｅ ｓｔａｔｉｓｔｉｃｓ ｏｆ ｔｈｅ ｇｒａｖｉｔｙ ｒｅ ｄｕｃｔｉｏｎｓ ａｒｅ ｓｈｏｗｎ ｉｎ Ｔａｂｌｅ １ ． Ｉｔ ｉｓ ｕｎｃｅｒｔａｉｎ ｗｈｅｔｈｅｒ ｔｈｅ ｂｉａｓｅｓ ｉｎ ｔｈｅ ＮＧＩＩ ｄａｔａ ａｒｅ ｉｎｄｉｃａｔｉｖｅ ｏｆ ｒｅｓｉｄｕａｌ ｓｙｓｔｅｍａｔｉｃ ｅｒｒｏｒｓ ．

， （），２００９，６７９ － ６９４

ＤｏｎｇＨａ ＬＥＥ ｅｔ ａｌ ． ／ Ｃｈｉｎａ Ｏｃｅａｎ Ｅｎｇｉｎｅｅｒｉｎｇ ２３ ４

６８６

Ｆｉｇ． ５ ． Ｌｏｃａｔｉｏｎ ｏｆ ｇｒａｖｉｔｙ ｐｏｉｎｔｓ ｕｓｅｄ ｆｏｒ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ｄｅｔｅｒｍｉｎａｔｉｏｎ． Ｂａｃｋｇｒｏｕｎｄ ｉｓ ｔｈｅ ｒｅｄｕｃｅｄ ｇｒａｖｉｔｙ ｄａｔａ ｗｉｔｈ ＧＧＭ０２Ｓ ／ ＥＧＭ９６ ａｎｄ ＲＴＭ．

Ｔａｂｌｅ １

Ｓｔａｔｉｓｔｉｃｓ ｏｆ ｄａｔａ ｒｅｄｕｃｔｉｏｎｓ ｆｏｒ ｇｒａｖｉｔｙ ｄａｔａ ａｆｔｅｒ ｒｅｍｏｖａｌ ｏｆ ｇｒｏｓｓ ｅｒｒｏｒ Ｄａｔａ ｓｅｔ

Ｎｏ． ｏｆ ｐｔｓ．

：

ｕｎｉｔ ｍｇａｌ

Ｍｅａｎ

Ｓｔｄ． ｄｅｖ．

Ｍｉｎ．

Ｍａｘ．

２１ ． ８６

２３ ． １７

－ ２６ ． ０５

１７６ ． １０

－ ９ ． ３７

２１ ． １２

－ ６５ ． １１

１４１ ． ９３

Ｒｅｄｕｃｅｄ ｆｏｒ ＧＧＭ０２Ｓ ／ ＥＧＭ９６ ａｎｄ ＲＴＭ

－ ３ ． ５１

１１ ． ７４

－ ５２ ． ２３

９８ ． ２０

ＢＧＩ ｄａｔａ

２４ ． ９８

１６ ． ９６

－ １６ ． ６０

１４７ ． ２０

－ ９ ． ７５

１６ ． ９２

－ ４４ ． ９９

１２２ ． ９８

Ｒｅｄｕｃｅｄ ｆｏｒ ＧＧＭ０２Ｓ ／ ＥＧＭ９６ ａｎｄ ＲＴＭ

０ ． ０４

１２ ． ２１

－ ３９ ． ５５

１２３ ． ６７

Ｎｅｗ ｇｒａｖｉｔｙ ｄａｔａ

２２ ． ９５

４ ． ６７

４ ． ８３

３１ ． ９８

－ ０ ． ９４

４ ． ６４

－ １１ ． ２２

９ ． ２３

３ ． ４６

４ ． １２

－ ９ ． ０３

９ ． ７８

ＮＧＩＩ ｄａｔａ Ｒｅｄｕｃｅｄ ｆｏｒ ＧＧＭ０２Ｓ ／ ＥＧＭ９６

２７５９

Ｒｅｄｕｃｅｄ ｆｏｒ ＧＧＭ０２Ｓ ／ ＥＧＭ９６

２０００

Ｒｅｄｕｃｅｄ ｆｏｒ ＧＧＭ０２Ｓ ／ ＥＧＭ９６

６１

Ｒｅｄｕｃｅｄ ｆｏｒ ＧＧＭ０２Ｓ ／ ＥＧＭ９６ ａｎｄ ＲＴＭ

２．３

Ｐｒａｃｔｉｃａｌ Ｃｏｍｐｕｔａｔｉｏｎ ｏｆ Ｇｒａｖｉｍｅｔｒｉｃ Ｇｅｏｉｄ

Ｔｈｅ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ｗａｓ ｃｏｍｐｕｔｅｄ ｆｒｏｍ ｔｈｅ ｒｅｄｕｃｅｄ ｇｒａｖｉｔｙ ｄａｔａ ｂｙ ｇｒｉｄｄｉｎｇ ｗｉｔｈ ｌｅａｓｔ ｓｑｕａｒｅｓ

，

（ ）

， ｔｉｏｎ ｏｆ ｔｈｅ ｔｅｒｒａｉｎ ｅｆｆｅｃｔｓ ａｎｄ ＧＧＭ０２Ｓ ／ ＥＧＭ９６ ｍｏｄｅｌ ｅｆｆｅｃｔｓ ． Ｄａｔａ ｗｅｒｅ ｇｒｉｄｄｅｄ ｉｎ ｔｈｅ ａｒｅａ ３３ ～ ３９° Ｎ，１２４ ～ １３１°Ｅ，ａｔ ａ ｂａｓｉｃ ｇｒｉｄ ｓｐａｃｉｎｇ ｏｆ ０ ． ０１２５° × ０ ． ０１６６６７° （ ０ ． ７５′ × １′）ｉｎ ｌａｔｉｔｕｄｅ ａｎｄ ｌｏｎｇｉ ｔｕｄｅ ． Ｔｈｅ ｃｏｍｐａｒｉｓｏｎ ｏｆ ｏｕｒ ｒｅｓｕｌｔｓ ｗｉｔｈ ｔｈｅ ｒｅｓｕｌｔｓ ｏｆ Ｙｕｎ （ １９９５）ｉｎｄｉｃａｔｅｓ ｔｈａｔ ２Ｄ ｍｕｌｔｉｂａｎｄ ｓｐｈｅｒ ｃｏｌｌｏｃａｔｉｏｎ ｔｗｏｄｉｍｅｎｓｉｏｎａｌ ２Ｄ ｍｕｌｔｉｂａｎｄ ｓｐｈｅｒｉｃａｌ ＦＦＴ ｔｒａｎｓｆｏｒｍａｔｉｏｎ ｆｏｌｌｏｗｅｄ ｂｙ ｔｈｅ ｒｅｓｔｏｒａ

ｉｃａｌ ＦＦＴ ｗｉｔｈ ４ｂａｎｄｓ ｈａｓ ｇｉｖｅｎ ｔｈｅ ｂｅｓｔ ｓｏｌｕｔｉｏｎ ｆｏｒ Ｋｏｒｅａ ． Ｔｈｅ ＦＦＴ ｗａｓ ｐｅｒｆｏｒｍｅｄ ｂｙ ｕｓｅ ｏｆ １００％ ｚｅｒｏｐａｄｄｉｎｇ ｔｏ ｌｉｍｉｔ ｔｈｅ ｐｅｒｉｏｄｉｃｉｔｙ ｅｆｆｅｃｔｓ ． Ｔｈｅ １００％ ｚｅｒｏｐａｄｄｉｎｇ ｐｕｔｓ ｚｅｒｏｓ ａｒｏｕｎｄ ｔｈｅ ｖａｌｕｅｓ ｏｆ

（

），ｐｒａｃｔｉｃａｌｌｙ ｄｏｕｂｌｉｎｇ ｔｈｅ ｄｉｍｅｎｓｉｏｎｓ （Ｙｕｎ，１９９９）．

ｔｈｅ ｏｒｉｇｉｎａｌ ｆｉｅｌｄ ｉｎｐｕｔ ｍａｔｒｉｘ

Ｔｈｅ ｆｉｎａｌ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ｉｎ Ｆｉｇ ． ６ ｗａｓ ｃｏｍｐｕｔｅｄ ｂｙ ２Ｄ ｓｐｈｅｒｉｃａｌ ＦＦＴ ｗｉｔｈ ４ｂａｎｄｓ ａｎｄ ８０ ～

，

９０ ｍｏｄｉｆｉｅｄ ＷｏｎｇＧｏｒｅ ｋｅｒｎｅｌ ． Ｆｉｇｓ ． ７ ａｎｄ ８ ｓｈｏｗ ｔｈｅ ｇｅｏｉｄ ｅｆｆｅｃｔｓ ｏｆ ｔｈｅ ｌｏｃａｌ ｇｒａｖｉｔｙ ｄａｔａ ａｎｄ ｔｈｅ ｅｆｆｅｃｔ ｆｒｏｍ ｔｈｅ ｔｏｐｏｇｒａｐｈｙ ．

， （），２００９，６７９ － ６９４

ＤｏｎｇＨａ ＬＥＥ ｅｔ ａｌ ． ／ Ｃｈｉｎａ Ｏｃｅａｎ Ｅｎｇｉｎｅｅｒｉｎｇ ２３ ４

６８７

Ｆｉｇ． ６ ． Ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ｏｆ Ｋｏｒｅａ ｉｎ ｇｌｏｂａｌ ｖｅｒｔｉｃａｌ ｒｅｆｅｒｅｎｃｅ ｓｙｓｔｅｍ

）

（ｕｎｉｔ：

ｍ ． Ｌｏｃａｔｉｏｎ ｏｆ ＧＰＳ ｌｅｖｅｌｉｎｇ ｄａｔａ

（

ｓｈｏｗｎ ｗｉｔｈ ｃｒｏｓｓ ａｌｌ ＧＰＳ ／ Ｌｅｖｅｌｉｎｇ

）ａｎｄ ｔｒｉａｎｇｌｅ （ＧＰＳ ／ Ｌｅｖｅｌｉｎｇ ）

ｄａｔａ

ｄａｔａ ａｒｏｕｎｄ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ．

Ｆｉｇ． ７ ． Ｇｅｏｉｄ ｅｆｆｅｃｔｓ ｆｒｏｍ ｔｈｅ ｔｒａｎｓｆｏｒｍａｔｉｏｎ ｏｆ ｔｈｅ

Ｆｉｇ． ８ ． Ｒｅｓｔｏｒｅｄ ｇｅｏｉｄ ｅｆｆｅｃｔｓ ｏｆ ｔｈｅ ｒｅｓｉｄｕａｌ

ｃｏｍｂｉｎｅｄ ｒｅｓｉｄｕａｌ ｇｒａｖｉｔｙ ｄａｔａ （ｕｎｉｔ：ｍ）．

（ ：）

ｔｅｒｒａｉｎ ｍｏｄｅｌ ｕｎｉｔ ｍ ．

Ｔａｂｌｅ ２ ｓｈｏｗｓ ｔｈｅ ｓｔａｔｉｓｔｉｃｓ ｏｆ ｃｏｍｐａｒｉｓｏｎ ｏｆ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ａｎｄ ｇｅｏｉｄ ｏｂｔａｉｎｅｄ ｆｒｏｍ ｔｗｏ ＧＰＳ ／ ｌｅｖｅｌｌｉｎｇ ｄａｔａ ｓｅｔｓ ｉｎ Ｆｉｇ ． ６ ． Ｔｈｅ Ｋｏｒｅａｎ ｇｅｏｉｄ ｍｏｄｅｌ ｈａｓ ａ ｒｅａｓｏｎａｂｌｅ ｂｕｔ ｎｏｔ ｔｈｅ ｂｅｓｔ ｆｉｔ ｔｏ ＧＰＳ ． Ｔｈｉｓ ｉｓ ｐｒｏｂａｂｌｙ ｄｕｅ ｔｏ ｔｈｅ ｓｙｓｔｅｍａｔｉｃ ｅｒｒｏｒｓ ａｎｄ ｏｕｔｌｉｅｒｓ ｉｎ ｂｏｔｈ ｇｒａｖｉｔｙ ａｎｄ ＧＰＳ ／ ｌｅｖｅｌｉｎｇ ｄａｔａ ． Ａｌｓｏ

，

ｓｈｏｗｎ ｉｓ ｔｈｅ ｆｉｔ ｉｎ ａ ｍｏｒｅ ｌｏｃａｌｉｚｅｄ ａｒｅａ ａｒｏｕｎｄ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ｕｓｅｄ ｉｎ ｔｈｅ ｆｉｎａｌ ｆｉｔ ｏｆ ｔｈｅ ｇｅｏｉｄ

（

）

ｔｏ ＧＰＳ ａｎｄ ｔｉｄｅ ＧＰＳ ／ Ｔｉｄｅｇａｕｇｅ ｏｂｓｅｒｖａｔｉｏｎｓ ｉｎ ｔｈｅ ｆｏｌｌｏｗｉｎｇ ｓｅｃｔｉｏｎ ． Ｔａｂｌｅ ２

：

Ｓｔａｔｉｓｔｉｃｓ ｏｆ ｃｏｍｐａｒｉｓｏｎ ｏｆ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ａｎｄ ｇｅｏｉｄ ｆｒｏｍ ＧＰＳ ／ ｌｅｖｅｌｉｎｇ ｕｎｉｔ ｍ ＧＰＳ ｄａｔａ ｓｅｔ

ＧＰＳ ／ Ｌｅｖｅｌｉｎｇ ｄａｔａ （ａｌｌ） ＧＰＳ ／ Ｌｅｖｅｌｉｎｇ ｄａｔａ （ａｒｏｕｎｄ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ）

Ｎｏ． ｏｆ ｐｔｓ

Ｍｅａｎ

Ｓｔｄ． ｄｅｖ．

Ｍｉｎ．

Ｍａｘ．

４４７

０ ． ８５

０ ． １６

－ ０ ． ２１

２ ． ００

７４

０ ． ８７

０ ． １２

０ ． ５３

１ ． ０９

３ ． ＧＰＳ ａｎｄ Ｔｉｄｅ Ｏｂｓｅｒｖａｔｉｏｎｓ Ｔｈｅ ｔｉｄｅ ｏｂｓｅｒｖａｔｉｏｎｓ ｗｅｒｅ ｃａｒｒｉｅｄ ｏｕｔ ｐｒｉｍａｒｉｌｙ ｔｏ ａｃｈｉｅｖｅ ａｎ ｉｎｄｅｐｅｎｄｅｎｔ ｃｏｎｔｒｏｌ ｏｆ ｔｈｅ ｇｒａｖｉｍｅｔ

， （），２００９，６７９ － ６９４

ＤｏｎｇＨａ ＬＥＥ ｅｔ ａｌ ． ／ Ｃｈｉｎａ Ｏｃｅａｎ Ｅｎｇｉｎｅｅｒｉｎｇ ２３ ４

６８８

，

ｒｉｃ ｇｅｏｉｄ ｈｅｉｇｈｔｓ ａｃｒｏｓｓ ｔｈｅ ｓｔｒａｉｔ ． Ｔｈｅ ｍｅａｓｕｒｅｍｅｎｔｓ ｗｅｒｅ ｄｏｎｅ ｏｖｅｒ ａ ｆｅｗ ｄａｙｓ ｕｓｉｎｇ ｒｅｐｅａｔｅｄ ｌｏｃａｌ ｍｅａｓｕｒｅｍｅｎｔｓ ｔｏ ｔｈｅ ｉｎｓｔａｎｔａｎｅｏｕｓ ｓｅａ ｓｕｒｆａｃｅ ｉｎ ｃｏｎｎｅｃｔｉｏｎ ｗｉｔｈ ｔｈｅ ｌａｙｏｕｔ ｏｆ ｔｈｅ ｔｉｄｅ ｇａｕｇｅｓ ． Ｔｈｒｅｅ

（ ，

，

）

（ ）

ｔｅｍｐｏｒａｒｙ ｔｉｄｅ ｇａｕｇｅｓ ＪＩＮ１ ＧＥＯ１ ＧＥＯ２ ａｎｄ ｏｎｅ ｐｅｒｍａｎｅｎｔ ｔｉｄｅ ｇａｕｇｅ ＧＡＤＯ ｗｅｒｅ ｕｓｅｄ ． Ｆｉｇ ． ９ ｓｈｏｗｓ ｔｈｅ ｌｏｃａｔｉｏｎ ｏｆ ｔｈｅ ｔｉｄｅ ｏｂｓｅｒｖａｔｉｏｎ ｐｏｉｎｔｓ ．

Ｆｉｇ． ９ ． Ｌｏｃａｔｉｏｎ ｏｆ ｔｉｄｅ ｏｂｓｅｒｖａｔｉｏｎ ｐｏｉｎｔｓ．

ＧＡＤＯ ｉｓ ａ ｐｅｒｍａｎｅｎｔ ｔｉｄｅ ｇａｕｇｅ ｆｏｒ ｔｈｅ ｈｙｄｒｏｇｒａｐｈｉｃ ｓｅｒｖｉｃｅ ． Ｔｈｉｓ ｔｉｄｅ ｇａｕｇｅ ｐｒｏｖｉｄｅｓ ａｎ ｕｎｉｎ ｔｅｒｒｕｐｔｅｄ ｔｉｍｅ ｓｅｒｉｅｓ ｆｏｒ ｔｈｅ ｄｕｒａｔｉｏｎ ｏｆ ａｌｌ ｏｔｈｅｒ ｔｉｄｅ ｇａｕｇｅｓ ． Ａ ２４ｈｒ ｅｒｒｏｒ ｉｎ ｔｈｅ ｔｉｍｅ ｔａｇｇｉｎｇ ｗａｓ ｏｂ

，

ｖｉｏｕｓ ａｎｄ ｃｏｒｒｅｃｔｅｄ ｐｒｉｏｒ ｔｏ ｔｈｅ ｕｓｅ ｏｆ ｄａｔａ ． ＪＩＮ１ ＧＥＯ１ ａｎｄ ＧＥＯ２ ｗｅｒｅ ｏｂｓｅｒｖｅｄ ｂｙ ａ ｔｅｍｐｏｒａｒｙ ｔｉｄｅ

（

），ａｓ ｔｈｅ ｓａｍｅ ｐｅｒｉｏｄ ２８０ ． ０ ～ ２８６ ． ０ ｉｎ Ｊｕｌｉａｎ ｄａｙ （ＪＤ）． ＪＩＮ１ ｐｒｏｖｉｄｅｄ ａｎ ｕｎｉｎｔｅｒｒｕｐｔｅｄ ｔｉｍｅ ｓｅｒｉｅｓ ｉｎ ｔｈｅ ｐｅｒｉｏｄ，２８１ ． ０９１ ～ ２８５ ． １７４ ｉｎ ＪＤ． Ｔｈｅｓｅ ｎｏｉｓ ｅｓ ｉｎ ｔｈｅ ｔｉｍｅ ｓｅｒｉｅｓ ｏｆ ＪＩＮ１ ｗｅｒｅ ｄｕｅ ｔｏ ｗａｖｅｓ ａｎｄ ｂｏａｔｓ ． ＧＥＯ１ ａｐｐａｒｅｎｔｌｙ ｍｏｖｅｄ ａｔ ｔｗｏ ｐｅｒｉｏｄｓ，ＪＤ： ２８２ ． ２２ （＃ １）ａｎｄ ２８３ ． ６２ （＃ ２）． Ｔｈｅ ａｍｏｕｎｔ ｏｆ ｍｏｖｅｍｅｎｔ ｗａｓ ｎｏｔ ｗｅｌｌ ｄｅｆｉｎｅｄ，ｔｈｅ ｄａｔａ ｗｅｒｅ ｏｎｌｙ ｕｓｅｄ ｉｎ ｔｈｅ ｐｅｒｉｏｄ，ＪＤ：２８０ ． ２４５ ～ ２８２ ． ２００ ． ＧＥＯ２ ｗａｓ ａｌｓｏ ｍｏｖｅｄ ａｔ ｔｗｏ ｐｅｒｉｏｄｓ，ＪＤ：２８２ ． ２４ （＃ １）ａｎｄ ２８４ ． ９５ （＃ ２）． Ｔｈｅｓｅ ｍｏｖｅｍｅｎｔｓ ｗｅｒｅ ｌｉｋｅｌｙ ａｔｔｒｉｂｕｔｅｄ ｔｏ ｃｕｒｉｏｕｓ ｐｅｏｐｌｅ ｔａｋｉｎｇ ｔｈｅ ｓｅｎｓｏｒ ｏｕｔ ｏｆ ｔｈｅ ｗａｔｅｒ ａｎｄ ｔｈｒｏｗｉｎｇ ｉｔ ｂａｃｋ ａｇａｉｎ ． Ｔｈｅ ｆｉｒｓｔ ｂｒｅａｋ ｗａｓ ｃｌｅａｒｌｙ ｆｉｘａｂｌｅ （ｏｆｆｓｅｔ ｏｆ １９ ． ８ ｃｍ）． Ｔｈｅ ｕｓｅｆｕｌ ｄａｔａ ｆｏｒ ｔｈｉｓ ｔｉｄｅ ｇａｕｇｅ ｗｅｒｅ ｔｈｅｒｅｆｏｒｅ ＪＤ：２８１ ． ２３８ ～ ２８４ ． ９５２ ． ｇａｕｇｅ ｄｕｒｉｎｇ ６ ｄａｙｓ ７ Ｏｃｔ ． ２００５ ｕｎｔｉｌ １２ Ｏｃｔ ． ２００５

Ｔｈｅ ｈｅｉｇｈｔ ｏｆ ｔｈｅ ＧＰＳ ｒｅｆｅｒｅｎｃｅ ｐｏｉｎｔ ａｂｏｖｅ ｔｈｅ ｍｅａｎ ｓｅａ ｌｅｖｅｌ ｆｏｒ ａ ｔｉｍｅ ｉｎｔｅｒｖａｌ

ｃａｎ ｂｅ ｅｘｐｒｅｓｓｅｄ ａｓ：

ｔ＝ ｔ Δ

，

－ ｔ１

（１０）

ＨΔ ｔ ＝ Ｈ ｍ ＋ ｔ ｍ － ｔΔ ｔ

； ，

２

ｗｈｅｒｅ ＨΔ ｔ ｉｓ ｔｈｅ ｍｅａｎ ｓｅａ ｌｅｖｅｌ ｄｕｒｉｎｇ ｔｈｅ ｐｅｒｉｏｄ Ｈ ｍ ｔｈｅ ｍｅａｓｕｒｅｄ ｈｅｉｇｈｔ ｄｉｆｆｅｒｅｎｃｅ ｆｒｏｍ ｔｈｅ ＧＰＳ

；，

；

，

ｒｅｆｅｒｅｎｃｅ ｐｏｉｎｔ ｔ ｍ ｔｈｅ ａｖｅｒａｇｅ ｔｉｄａｌ ｒｅａｄｉｎｇｓ ｃｏｒｒｅｓｐｏｎｄｉｎｇ ｔｏ ｔｈｅ ｍｅａｓｕｒｅｍｅｎｔｓ Ｈ ｍ ａｎｄ ｔΔ ｔ ｔｈｅ ａｖｅｒａｇｅ ｔｉｄｅ ｇａｕｇｅ ｒｅａｄｉｎｇｓ ｄｕｒｉｎｇ ｔｈｅ ｌｏｎｇｅｒ ｐｅｒｉｏｄ

；

，

ｔ ． Ｔｈｅ ｈｅｉｇｈｔ ＨΔｗｉｌｌ ｂｅ ａｆｆｅｃｔｅｄ ｂｙ ｌｏｎｇｐｅｒｉ Δ ｔ

，

ｏｄ ｔｉｄｅｓ ｈｏｗｅｖｅｒ ｔｈｅ ｈｅｉｇｈｔ ｄｉｆｆｅｒｅｎｃｅｓ ｂｅｔｗｅｅｎ ｓｔａｔｉｏｎｓ ｃｏｒｒｅｓｐｏｎｄ ｔｏ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔ ｄｉｆｆｅｒｅｎｃｅｓ ｓｉｎｃｅ ｔｈｅ ｌｏｎｇｐｅｒｉｏｄ ｔｉｄｅｓ ｍａｙ ａｓｓｕｍｅｄｌｙ ｂｅ ｃｏｍｍｏｎ ｔｏ ａｌｌ ｓｔａｔｉｏｎｓ ．

，

，

Ｂａｓｅｄ ｏｎ ｔｈｅ ｕｓｅｆｕｌ ｔｉｄａｌ ｉｎｔｅｒｖａｌｓ ｗｉｔｈ ｔｈｅ ＧＡＤＯ ｔｉｄｅ ｇａｕｇｅ ｕｓｅｄ ａｓ ａ ｒｅｆｅｒｅｎｃｅ ａｎｄ ｔｈｅ ｐａｉｒｓ ｏｆ ｔｈｅ ｈｅｉｇｈｔｓ ａｂｏｖｅ ｓｅａ ｌｅｖｅｌ ｆｏｒ ｔｈｅ ｒｅｌｅｖａｎｔ ｉｎｔｅｒｖａｌｓ

ｔ ｗｅｒｅ ｄｅｔｅｒｍｉｎｅｄ，ａｓ ｓｈｏｗｎ ｉｎ Ｔａｂｌｅ ３ ． Ｔｈｅ Δ

， （），２００９，６７９ － ６９４

ＤｏｎｇＨａ ＬＥＥ ｅｔ ａｌ ． ／ Ｃｈｉｎａ Ｏｃｅａｎ Ｅｎｇｉｎｅｅｒｉｎｇ ２３ ４

６８９

Ｆｉｇ． １０ ． Ｓｅａ ｌｅｖｅｌ ｈｅｉｇｈｔ ｖａｒｉａｔｉｏｎ ｆｒｏｍ ｔｉｄａｌ ｏｂｓｅｒｖａｔｉｏｎｓ．

（）

Ｄｔｅｒｍ ｉｎ Ｔａｂｌｅ ３ ｉｎｄｉｃａｔｉｎｇ ｔｈｅ ｔｅｒｍ Ｄ ＝ Ｈ ｍ ＋ ｔ ｍ ｉｎ Ｅｑ ． １０ ｓｈｏｕｌｄ ｂｅ ｃｏｎｓｔａｎｔ ． Ｔａｂｌｅ ３

Ｈｅｉｇｈｔｓ ａｂｏｖｅ ｌｏｃａｌ ｓｅａ ｌｅｖｅｌ ｆｏｒ ｓｐｅｃｉｆｉｃ ｐｅｒｉｏｄ ｉｎｔｅｒｖａｌｓ Ｎｏ． ｏｆ

Ｍｅａｎ ｏｆ

Ｓｔｄ． ｄｅｖ． ｏｆ

Ｈ ｍ ｍｅａｓ

Ｄｔｅｒｍ

Ｄｔｅｒｍ

ＧＡＤＯ

１６

３ ． ３４

０ ． ０１４

ＧＡＤＯ





ＧＡＤＯ



ＪＩＮ１

Ｓｔａｔｉｏｎ

（）

Ｉｎｔｅｒｖａｌ ＪＤ

Ｎｏ． ｏｆ

ＨΔ ｔ

ｔｉｄａｌ ｍｅａｓ

（ｍ）

Ｃｏｍｍｅｎｔｓ

２８１ ． ０９１ ～ ２８５ ． １７４

５８６０

２ ． ３１６

Ｕｓｅ ｗｉｔｈ ＪＩＮ１



２８０ ． ２４５ ～ ２８２ ． ２００

２８０６

２ ． ２８５

Ｕｓｅ ｗｉｔｈ ＧＥＯ１





２８１ ． ２３８ ～ ２８４ ． ９５２

５３３０

２ ． ３２３

Ｕｓｅ ｗｉｔｈ ＧＥＯ２

６

３ ． １７５

０ ． ０２２

２８１ ． ０９１ ～ ２８５ ． １７４

２９４０

１ ． ９３３

ＧＥＯ１

６

４ ． ４８０

０ ． ０４９

２８０ ． ２４５ ～ ２８２ ． ２００

２８１

３ ． ３０８

ＧＥＯ２

６

４ ． ５５２

０ ． ０２３

２８１ ． ２３８ ～ ２８４ ． ９５２

５３４３

１ ． ７８１

Ｔｈｅ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔｓ ｉｎ ｔｈｅ Ｋｏｒｅａｎ ｈｅｉｇｈｔ ｓｙｓｔｅｍ ｃａｎ ｓｕｂｓｅｑｕｅｎｔｌｙ ｂｅ ｆｏｕｎｄ ｆｒｏｍ ｔｈｅ ｈｅｉｇｈｔ

，

ｄｉｆｆｅｒｅｎｃｅｓ ｏｆ Ｔａｂｌｅ ３ ． Ｆｏｒ ｔｈｉｓ ｐｕｒｐｏｓｅ ｔｈｅ ｈｅｉｇｈｔｓ ｗｅｒｅ ｔｒａｎｓｆｅｒｒｅｄ ｆｒｏｍ ｔｈｅ ｆｕｎｄａｍｅｎｔａｌ ｂｅｎｃｈｍａｒｋ ＪＩＮＨ ｔｏ ＪＩＮ１ ｂｙ ｕｓｅ ｏｆ ｔｈｅ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ． Ｔｈｉｓ ｔｒａｎｓｆｅｒ ｗａｓ ａｓｓｕｍｅｄ ｔｏ ｂｅ ａｃｃｕｒａｔｅ ｓｉｎｃｅ ｔｈｅ ｂａｓｅ

，

ｌｉｎｅ ｌｅｎｇｔｈ ｗａｓ ｖｅｒｙ ｓｍａｌｌ ． Ｔｈｅ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔ ｏｆ ＪＩＮ１ ｗａｓ ｓｕｂｓｅｑｕｅｎｔｌｙ ｔｒａｎｓｆｅｒｒｅｄ ｔｏ ＧＡＤＯ ＧＥＯ１ ａｎｄ ＧＥＯ２ ｂｙ ｕｓｅ ｏｆ ｔｈｅ ｒｅｌｅｖａｎｔ ｈｅｉｇｈｔ ｄｉｆｆｅｒｅｎｃｅｓ ｉｎ Ｔａｂｌｅ ４ ．

，

Ｔｈｅ ＧＰＳ ｃｏｏｒｄｉｎａｔｅ ｏｆ ＪＩＮＨ ｄｅｆｉｎｅｓ ｔｈｅ ｄａｔｕｍ ｏｆ ｔｈｅ ＧＰＳ ｓｙｓｔｅｍ． Ｓｉｍｉｌａｒｌｙ ｔｈｅ Ｎ ＧＰＳ ／ Ｔｉｄｅ ｄｅ

，

，

，

ｆｉｎｅｓ ｔｈｅ ｌｏｃａｌ ｇｅｏｉｄ ｄａｔｕｍ ｉ ． ｅ ． ｔｈｅ ｇｅｏｉｄ ｖａｌｕｅ ｗｈｅｎ ｓｕｂｔｒａｃｔｅｄ ｆｒｏｍ ＧＰＳ ｗｉｌｌ ｇｉｖｅ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔｓ ｉｎ ｔｈｅ Ｋｏｒｅａｎ ｎａｔｉｏｎａｌ ｈｅｉｇｈｔ ｓｙｓｔｅｍ． Ｔｈｅ ｃｏｍｐａｒｉｓｏｎ ｏｆ Ｔａｂｌｅｓ ３ ａｎｄ ４ ｓｈｏｗｓ ｔｈａｔ ｔｈｅ ｌｏｃａｌ ｓｅａ ｌｅｖｅｌ ｉｎ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ｉｓ ｃｌｏｓｅ ｔｏ ｔｈｅ ｚｅｒｏ ｈｅｉｇｈｔ ｏｆ Ｋｏｒｅａ ． ＧＡＤＯ ｈａｓ ｔｈｕｓ ａｎ ｏｒｔｈｏｍｅｔｒｉｃ

； ， ２ ． ３２３ ｍ． Ｔｈｅｒｅｆｏｒｅ，ｔｈｅ ｄｉｆｆｅｒｅｎｃｅ ｉｎ ｈｅｉｇｈｔ ｗａｓ ２５ ｃｍ，ｓｍａｌｌｅｒ ｔｈａｎ ｔｈｅ ｔｙｐｉｃａｌ ｓｅａｓｏｎａｌ ｖａｒｉａｔｉｏｎｓ ｏｆ ｈｅｉｇｈｔ ｏｆ ２ ． ５７１ ｍ ｈｏｗｅｖｅｒ ｉｔ ｈａｓ ａ ｈｅｉｇｈｔ ａｂｏｖｅ ｔｈｅ ｍｅａｎ ｓｅａ ｌｅｖｅｌ ｆｏｒ ｆｏｕｒ ｄａｙｓ ｉｎ Ｏｃｔｏｂｅｒ ２００５ ｏｆ

， （），２００９，６７９ － ６９４

ＤｏｎｇＨａ ＬＥＥ ｅｔ ａｌ ． ／ Ｃｈｉｎａ Ｏｃｅａｎ Ｅｎｇｉｎｅｅｒｉｎｇ ２３ ４

６９０

ｓｅａ ｌｅｖｅｌ ． Ｔｈｅ ｌａｓｔ ｃｏｌｕｍｎ ｉｎ Ｔａｂｌｅ ４ ｇｉｖｅｓ ａｎ ｅｓｔｉｍａｔｅ ｏｆ ｔｈｅ ｇｅｏｉｄ ｄｉｆｆｅｒｅｎｃｅ ｂｅｔｗｅｅｎ ｔｈｅ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ

（

）

Ｎ ｇｒａｖ ａｎｄ ｔｈｅ ｉｎｄｅｐｅｎｄｅｎｔ ｇｅｏｉｄ ｉｎｆｏｒｍａｔｉｏｎ Ｎ ＧＰＳ ／ Ｔｉｄｅ ｄｅｒｉｖｅｄ ｂｙ

（

Ｎ ＧＰＳ ／ Ｔｉｄｅ ＝ ｈ ＧＰＳ ／ Ｔｉｄｅ － Ｈ ｌｅｖｅｌｌｉｎｇ ．

）

（１１）

，

Ｔｈｅ ｇｅｏｉｄ ｄｉｆｆｅｒｅｎｃｅ Ｎ ＧＰＳ ／ Ｔｉｄｅ － Ｎ ｇｒａｖ ｈａｓ ａ ｍｅａｎ ｖａｌｕｅ ｏｆ ０ ． ５８６ ｍ ａｎｄ ａ ｓｔａｎｄａｒｄ ｄｅｖｉａｔｉｏｎ ｏｆ ４ ． ８ ｃｍ． Ｔｈｅ ｓｔａｎｄａｒｄ ｄｅｖｉａｔｉｏｎ ｉｓ ｑｕｉｔｅ ｓａｔｉｓｆａｃｔｏｒｙ ｔｏ ｇｉｖｅ ｔｈｅ ｌａｒｇｅ ｅｘｐｅｃｔｅｄ ｃｈａｎｇｅｓ ｉｎ ｔｈｅ ｇｒａｖｉｔｙ

，

，

ｆｉｅｌｄ ｔｈｅ ｌａｃｋ ｏｆ ａｃｃｕｒａｔｅ ｇｒａｖｉｔｙ ｄａｔａ ｏｆｆｓｈｏｒｅ ａｎｄ ｔｈｅ ｎｏｉｓｅ ｉｎ ｔｈｅ ｔｉｄａｌ ｈｅｉｇｈｔ ｔｒａｎｓｆｅｒ ． Ｔｈｅ ｍｅａｎ

（

ｖａｌｕｅ ｏｆ ｔｈｅ ｇｅｏｉｄ ｄｉｆｆｅｒｅｎｃｅ ｅｘｐｒｅｓｓｅｓ ｔｈｅ ｄｉｆｆｅｒｅｎｃｅ ｂｅｔｗｅｅｎ ｔｈｅ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ｒｅｆｅｒｒｉｎｇ ｔｏ ａ ｇｌｏｂ

）

（）

ａｌ ｄａｔｕｍ ａｎｄ ｔｈｅ ＧＰＳ ｇｅｏｉｄ ｆｒｏｍ Ｅｑ ． １１ ｗｈｉｃｈ ｂｙ ｄｅｆｉｎｉｔｉｏｎ ｒｅｆｅｒｓ ｔｏ ａ ｌｏｃａｌ ｓｅａｌｅｖｅｌ ｄａｔｕｍ． Ｔａｂｌｅ ４ ｓｈｏｗｓ ｔｈａｔ ｔｈｅ ｂｅｎｃｈｍａｒｋ ＧＥＯＪ ｏｎ ｔｈｅ Ｇｅｏｊｅ ｉｓｌａｎｄ ｄｏｅｓ ｎｏｔ ｆｉｔ ｔｈｅ ｇｅｏｉｄ ａｎｄ ｔｉｄａｌ ｈｅｉｇｈｔ

；

，

（

）

ｔｒａｎｓｆｅｒ ｔｈｅｒｅｆｏｒｅ ｔｈｅｒｅ ｉｓ ｌｉｋｅｌｙ ａｎ ｅｒｒｏｒ ｉｎ ｏｎｅ ｏｆ ｔｈｅ ｈｅｉｇｈｔｓ ｅｉｔｈｅｒ ｌｅｖｅｌｅｄ ｈｅｉｇｈｔ ｏｒ ＧＰＳ ｈｅｉｇｈｔ

，

ｆｏｒ ｔｈｉｓ ｓｔａｔｉｏｎ ． Ｔｈｅｒｅｆｏｒｅ ｔｈｅ Ｎ ＧＰＳ ／ Ｔｉｄｅ ｖａｌｕｅ ｏｆ ＧＥＯＪ ｉｓ ｎｏｔ ｕｓｅｄ ｉｎ ｔｈｅ ｆｉｔｔｉｎｇ ｆｏｒ ｆｉｎａｌ ｇｅｏｉｄ ｄｅｔｅｒ ｍｉｎａｔｉｏｎ ． Ｔａｂｌｅ ４

， ｅｌｌｉｐｓｏｉｄａｌ ｈｅｉｇｈｔｓ ｆｒｏｍ ＧＰＳ ｓｕｒｖｅｙ ａｎｄ ｄｅｒｉｖｅｄ ｇｅｏｉｄ ｈｅｉｇｈｔｓ （Ｎ

Ｈｅｉｇｈｔｓ ｉｎ ｔｈｅ Ｋｏｒｅａｎ ｈｅｉｇｈｔ ｓｙｓｔｅｍ ｏｆ ｔｈｅ ｔｉｄｅ ｇａｕｇｅ ｐｏｉｎｔｓ ｗｉｔｈ ｃｏｒｒｅｓｐｏｎｄｉｎｇ

）

ＧＰＳ ／ Ｔｉｄｅ

Ｎ ＧＰＳ ／ Ｔｉｄｅ

Ｎ ＧＰＳ ／ Ｔｉｄｅ

＝ ｈ－ Ｈ

－ Ｎ ｇｒａｖ

２８ ． ４１４

２８ ． ９５２

０ ． ５３８

２ ． １８８

２８ ． ３９４

２８ ． ９３２

０ ． ５３８

３１ ． ５６２

２ ． ５７１

２８ ． ３９２

２８ ． ９９１

０ ． ５９９

１２８ ． ６９６８０

３２ ． ５６９

３ ． ５９４

２８ ． ３２７

２８ ． ９７５

０ ． ６４８

３４ ． ８９０６１

１２８ ． ６９９１８

３０ ． ９８２

２ ． ０２９

２８ ． ３９３

２８ ． ９５３

０ ． ５６０

３４ ． ８８９４０

１２８ ． ６９０４９

４８ ． ０９０

１８ ． ７６３

２８ ． ３９９

２９ ． ３２７

０ ． ９２８

Ｅｌｌｉｐｓｏｉｄ

Ｏｒｔｈｏｍ．

ｈｅｉｇｈｔ ｈ

ｈｅｉｇｈｔ Ｈ

１２８ ． ８１７２８

３６ ． ４８９

７ ． ５３７

３５ ． ０９３１０

１２８ ． ７８９３２

３１ ． １２０

ＧＡＤＯ

３５ ． ０２４１６

１２８ ． ８１０８７

ＧＥＯ１

３４ ． ９８０６４

ＧＥＯ２

ＧＥＯＪ

Ｐｏｉｎｔ

Ｌａｔ．

Ｌｏｎ．

ＪＩＮＨ

３５ ． １０００１

ＪＩＮ１

，

，

Ｎ ｇｒａｖ

Ｃｏｍｍｅｎｔｓ

Ｄａｔｕｍ ｐｏｉｎｔ Ｏｒｔｈｏ ｈ ｆｒｏｍ ＪＩＮＨ ｂｙ ｇｅｏｉｄ Ｏｒｔｈｏ ｈ ｆｒｏｍ ＪＩＮ１ ｂｙ ｔｉｄａｌ Ｏｒｔｈｏ ｈ ｆｒｏｍ ＧＡＤＯ Ｏｒｔｈｏ ｈ ｆｒｏｍ ＧＡＤＯ Ｇｅｏｊｅ ｉｓｌａｎｄ Ｂｅｎｃｈ Ｍａｒｋ

４ ． Ｆｉｎａｌ Ｇｅｏｉｄ Ｍｏｄｅｌ ｆｏｒ ｔｈｅ Ｃｏｎｓｉｓｔｅｎｔ Ｈｅｉｇｈｔ Ｓｙｓｔｅｍ

，

Ｔｈｅ ｆｉｎａｌ ｇｅｏｉｄ ｏｆ ｔｈｅ ａｒｅａ ｗａｓ ｄｅｔｅｒｍｉｎｅｄ ｂｙ ｆｉｔｔｉｎｇ ｔｈｅ Ｎ ＧＰＳ ／ Ｔｉｄｅ ｖａｌｕｅｓ ｏｆ Ｔａｂｌｅ ４ ｔｏｇｅｔｈｅｒ ｗｉｔｈ

，

ａ ｓｕｂｓｅｔ ｏｆ ｔｈｅ ＧＰＳ ／ ｌｅｖｅｌｉｎｇ ｄａｔａ ｐｏｉｎｔｓ ａｓ ｓｈｏｗｎ ｉｎ Ｆｉｇ ． ７ ． Ｂｙ ｕｓｅ ｏｆ ｔｈｅ ｇｅｏｉｄ ｉｎｆｏｒｍａｔｉｏｎ ｆｒｏｍ

，

ＧＰＳ ／ ｌｅｖｅｌｉｎｇ ｏｒ ｔｈｅ ＧＰＳ ／ Ｔｉｄｅ ｏｂｓｅｒｖａｔｉｏｎ ｌｉｎｋｓ ｔｏｇｅｔｈｅｒ ｗｉｔｈ ｇｅｏｉｄ ｆｒｏｍ ｇｒａｖｉｍｅｔｒｙ ｔｈｅ ｌｏｎｇ ｗａｖｅ

，

ｌｅｎｇｔｈ ｇｅｏｉｄ ｅｒｒｏｒｓ ｃａｎ ｂｅ ｓｕｐｐｒｅｓｓｅｄ ａｎｄ ｔｈｅ ｉｎｈｅｒｅｎｔ ｄａｔｕｍ ｄｉｆｆｅｒｅｎｃｅｓ ｅｌｉｍｉｎａｔｅｄ ． Ｉｔ ｉｓ ｅｓｓｅｎｔｉａｌ

，

ｗｈｅｎ ｃｏｍｐｕｔｉｎｇ ＧＰＳ ｇｅｏｉｄ ｈｅｉｇｈｔｓ ｔｈａｔ ｂｏｔｈ ｌｅｖｅｌｉｎｇ ａｎｄ ＧＰＳ ｈｅｉｇｈｔｓ ａｒｅ ａｓ ｅｒｒｏｒ ｆｒｅｅ ａｓ ｐｏｓｓｉｂｌｅ

， （），２００９，６７９ － ６９４

ＤｏｎｇＨａ ＬＥＥ ｅｔ ａｌ ． ／ Ｃｈｉｎａ Ｏｃｅａｎ Ｅｎｇｉｎｅｅｒｉｎｇ ２３ ４

６９１

ｏｔｈｅｒｗｉｓｅ ｔｈｅｓｅ ｅｒｒｏｒｓ ｗｉｌｌ ｃｒｅｅｐ ｉｎｔｏ ｔｈｅ ｆｉｔｔｅｄ ｇｅｏｉｄ ． Ｔｈｅ ｆｉｔｔｉｎｇ ｏｆ ａ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ｔｏ ａ ｓｅｔ ｏｆ ＧＰＳ ｇｅｏｉｄ ｈｅｉｇｈｔｓ ｅｎｔａｉｌｓ ｍｏｄｅｌｌｉｎｇ ｔｈｅ ｄｉｆｆｅｒｅｎｃｅ ｓｉｇｎａｌ

ε＝ Ｎ

， ε

ＧＰＳ ／ Ｔｉｄｅ

（１２）

－ Ｎ ｇｒａｖ

ａｎｄ ａｄｄｉｎｇ ｔｈｅ ｍｏｄｅｌｌｅｄ ｃｏｒｒｅｃｔｉｏｎ ｓｉｇｎａｌ ｔｏ ｔｈｅ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ． Ｉｎ ｔｈｉｓ ｗａｙ ａ ｎｅｗ ｇｅｏｉｄ ｇｒｉｄ ｔｕｎｅｄ ｔｏ ｔｈｅ ｌｅｖｅｌｉｎｇ ａｎｄ ＧＰＳ ｄａｔｕｍ ｃａｎ ｂｅ ｏｂｔａｉｎｅｄ ． Ｔｈｅ ｍｅｔｈｏｄ ｏｆ ｌｅａｓｔ ｓｑｕａｒｅｓ ｃｏｌｌｏｃａｔｉｏｎ ｉｓ ｕｓｅｄ ｆｏｒ ｅｓ

，

ｔｉｍａｔｉｎｇ ｔｈｅ ｔｒｅｎｄ ａｎｄ ｍｏｄｅｌｌｉｎｇ ｔｈｅ ｒｅｓｉｄｕａｌｓ ． Ｆｏｒ ｔｒｅｎｄ ｅｓｔｉｍａｔｉｏｎ ｗｅ ｕｓｅｄ ｔｈｅ ４ｐａｒａｍｅｔｅｒ ｍｏｄｅｌ

（Ｈｅｉｓｋａｎｅｎ ａｎｄ Ｍｏｒｉｔｚ，１９６７）ｏｆ ｔｈｅ ｆｏｒｍ， （１３） ｃｏｓλ ｓｉｎλ ′， ε＝ ｃｏｓφ α＋ ｃｏｓφ α＋ ｓｉｎφ α＋ α＋ ε ｗｈｅｒｅ，ｐａｒａｍｅｔｅｒｓαｔｏαｃｏｒｒｅｓｐｏｎｄ ｔｏ ｔｈｅ ｇｅｏｉｄ ｅｆｆｅｃｔｓ ｏｆ ａ Ｈｅｌｍｅｒｔ ｔｒａｎｓｆｏｒｍａｔｉｏｎ ． ε ′ ｉｓ ｔｈｅ ｒｅｓｉｄ ｕａｌ ｇｅｏｉｄ ｅｒｒｏｒｓ ｍｏｄｅｌｅｄ ｂｙ ｌｅａｓｔ ｓｑｕａｒｅ ｃｏｌｌｏｃａｔｉｏｎ ｏｆ ｔｈｅ ｆｏｒｍ， ＾ ＝ Ｃ Ｃ ｛ （１４） ′ ｝． ε ε Ｉｎ ｔｈｅ ｃｏｌｌｏｃａｔｉｏｎ ｐｒｏｃｅｓｓ，ａ ｃｏｖａｒｉａｎｃｅ ｆｕｎｃｔｉｏｎ ｍｕｓｔ ｂｅ ａｓｓｕｍｅｄ ｆｏｒ ｔｈｅ ｒｅｓｉｄｕａｌ ｇｅｏｉｄ ｅｒｒｏｒｓε ′ ａｓ ａ ｆｕｎｃｔｉｏｎ ｏｆ ｄｉｓｔａｎｃｅ ｓ ． Ｗｅ ｕｓｅｄ ｔｈｅ ｓｅｃｏｎｄｏｒｄｅｒ Ｍａｒｋｏｖ ｃｏｖａｒｉａｎｃｅ ｆｕｎｃｔｉｏｎ ｏｆ ｔｈｅ ｆｏｒｍ （Ｉｌｉｆｆｅ ｅｔ ， ａｌ ． ，２００３） （１５） Ｃ （ｓ ）＝ Ｃ （ １ ＋α ｓ ）ｅ α． Ｓｕｃｈ ａ ｃｏｖａｒｉａｎｃｅ ｆｕｎｃｔｉｏｎ ｉｓ ｃｈａｒａｃｔｅｒｉｚｅｄ ｂｙ ｔｈｅ ｖａｒｉａｎｃｅ Ｃ ａｎｄ ｃｏｒｒｅｌａｔｉｏｎ ｌｅｎｇｔｈ ｓ ，ｗｈｉｃｈ， ｉｎ ｔｕｒｎ，ｄｅｔｅｒｍｉｎｅ ｔｈｅ ｄｅｇｒｅｅ ｏｆ ｆｉｔ ａｎｄ ｔｈｅ ｓｍｏｏｔｈｎｅｓｓ ｏｆ ｔｈｅ ｉｎｔｅｒｐｏｌａｔｅｄ ｇｅｏｉｄ ｅｒｒｏｒ ． Ｔｈｅ ｃｏｎｓｔａｎｔ α ｉｓ ｔｈｅ ｏｎｌｙ ｑｕａｎｔｉｔｙ ｔｏ ｂｅ ｓｐｅｃｉｆｉｅｄ，ｗｉｔｈ Ｃ ａｕｔｏｍａｔｉｃａｌｌｙ ａｄａｐｔｅｄ ｔｏ ｔｈｅ ｄａｔａ ． Ｉｎ ｔｈｅ ｓｅｌｅｃｔｉｏｎ ｏｆ ｔｈｅ ｃｏｒｒｅｌａｔｉｏｎ ｌｅｎｇｔｈ ａｎｄ ｎｏｉｓｅ ｏｆ ｏｂｓｅｒｖｅｄ ｅｒｒｏｒｓ，ｔｈｅ ｕｓｅｒ ｈａｓ ａ ｗｉｄｅ ｒａｎｇｅ ｏｆ ｓｅｌｅｃｔｉｏｎ ｏｐｔｉｏｎｓ ． Ｅｉｔｈｅｒ ａ ｓｔｒｏｎｇ ｆｉｔ ｔｏ ｔｈｅ ＧＰＳ ｄａｔａ，ｏｒ ａ ｍｏｒｅ ｒｅｌａｘｅｄ ｆｉｔ，ｃａｎ ｄｉｍｉｎｉｓｈ ｔｈｅ ｉｍｐａｃｔ ｏｆ ｐｏｓｓｉｂｌｅ ｅｒｒｏｒｓ ｉｎ ｔｈｅ ＧＰＳ ｄａｔａ （Ｆｏｒｓｂｅｒｇ ｅｔ ａｌ ． ，２００３）． Ｗｅ ｄｅｔｅｒｍｉｎｅｄ ａ ｃｏｒｒｅｌａｔｉｏｎ ｌｅｎｇｔｈ ｏｆ ４０ ｋｍ ｂｙ ａ ｓｔａｎｄａｒｄ ＬＳＣ ｐｒｏｃｅ ｄｕｒｅ （ Ｍｏｒｉｔｚ，１９８０）ａｎｄ １ ｃｍ ａ ｐｒｉｏｒｉ ＧＰＳ ｎｏｉｓｅ ｗａｓ ａｓｓｕｍｅｄ ． Ｔｈｅ ｇｅｏｉｄ ｃｏｒｒｅｃｔｉｏｎ ｓｉｇｎａｌ ｉｓ ｓｈｏｗｎ １

１

２

３

４

４

－１ ｘｘ ｓｘ

－ ｓ

０

０

１／２

０

ｉｎ Ｆｉｇ ． １１ ．

（

）

Ｔｈｅ ｐｏｓｔｆｉｔ ｅｒｒｏｒ ｓｔｄ ． ｄｅｖ ． ｏｆ ｔｈｅ ｆｉｎａｌ ｇｅｏｉｄ ｔｏ ｔｈｅ Ｎ ＧＰＳ ／ Ｔｉｄｅ ｄｅｒｉｖｅｄ ｆｒｏｍ ＧＰＳ ａｎｄ ｔｉｄｅ ｏｂｓｅｒ ｖａｔｉｏｎ ｉｓ ± ２ ． ２ ｃｍ ｆｏｒ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ａｎｄ ± ４ ． ３ ｃｍ ｆｏｒ ｔｈｅ ｏｔｈｅｒ ａｒｅａｓ ａｒｏｕｎｄ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ． Ｔａｂｌｅ ５ ｓｈｏｗｓ ｔｈｅ ｐｏｓｔｆｉｔ ｓｔａｔｉｓｔｉｃｓ ｆｏｒ ｔｈｅ ｌｏｃａｌ ｄａｔｕｍ． Ｔｈｅ ｆｉｎａｌ ｌｏｃａｌ ｐｒｅｃｉｓｅ ｇｅｏｉｄ ｉｓ ｓｈｏｗｎ ｉｎ Ｆｉｇ ． １２ ． Ｔａｂｌｅ ５

：

Ｐｏｓｔｆｉｔ ｓｔａｔｉｓｔｉｃｓ ｆｏｒ ｔｈｅ ｌｏｃａｌ ｄａｔｕｍ Ｄａｔａ ｓｅｔ

Ｎ Ｔｉｄｅ ｄｅｒｉｖｅｄ ｆｒｏｍ ＧＰＳ ／ ｔｉｄｅ ｇａｕｇｅ ｆｉｘ

（

）

ｐｏｉｎｔｓ ｉｎ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ

Ｎ ＧＰＳ ｄｅｒｉｖｅｄ ｆｒｏｍ ＧＰＳ ／ Ｌｅｖｅｌｉｎｇ ｄａｔａ

（ａｒｏｕｎｄ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ）

ｕｎｉｔ ｍ

Ｎｏ． ｏｆ ｐｔｓ

Ｍｅａｎ

Ｓｔｄ． ｄｅｖ．

Ｍｉｎ．

Ｍａｘ．

５

０ ． ００１

± ０ ． ０２２

－ ０ ． ０２３

０ ． ０２４

７４

－ ０ ． ００２

± ０ ． ０４３

－ ０ ． １１２

０ ． １１８

５ ． Ｒｅｓｕｌｔｓ ａｎｄ Ｃｏｎｃｌｕｓｉｏｎｓ

，

Ｉｎ ｔｈｅ ｐｒｅｓｅｎｔ ｓｔｕｄｙ ｗｅ ｄｅｖｅｌｏｐｅｄ ａ ｌｏｃａｌ ｐｒｅｃｉｓｅ ｇｅｏｉｄ ｍｏｄｅｌ ｔｏ ｓｏｌｖｅ ｔｈｅ ｐｒｏｂｌｅｍ ｏｆ ｔｈｅ ｉｎｃｏｎ

， （），２００９，６７９ － ６９４

ＤｏｎｇＨａ ＬＥＥ ｅｔ ａｌ ． ／ Ｃｈｉｎａ Ｏｃｅａｎ Ｅｎｇｉｎｅｅｒｉｎｇ ２３ ４

６９２

Ｆｉｇ． １１ ． Ｇｅｏｉｄ ｃｏｒｒｅｃｔｉｏｎ ｓｉｇｎａｌ ｆｒｏｍ ｌｏｃａｌ ＧＰＳ ／ ｌｅｖｅｌ

（ ：）

Ｆｉｇ． １２ ． Ｆｉｎａｌ ｇｅｏｉｄ ｆｏｒ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ｏｆ Ｇｅｏｇａ

ｉｎｇ ／ ｔｉｄｅ ｇａｕｇｅ ｄａｔａ ｕｎｉｔ ｍ ．

，

Ｇｒａｎｄ Ｂｒｉｄｇｅ ｆｒｏｍ ｇｒａｖｉｔｙ ＧＰＳ ａｎｄ ｔｉｄｅ ｇａｕｇｅ

，

ｍｅａｓｕｒｅｍｅｎｔｓ ｗｉｔｈ ｃｏｎｔｏｕｒ ｉｎｔｅｒｖａｌ ｏｆ ２ ｃｍ．

ｓｉｓｔｅｎｔ ｈｅｉｇｈｔ ｓｙｓｔｅｍ ｏｃｃｕｒｒｉｎｇ ｉｎ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ｏｆ ｔｈｅ Ｇｅｏｇａ Ｇｒａｎｄ Ｂｒｉｄｇｅ ｃａｕｓｅｄ ｂｙ ｔｈｅ ｌａｒｇｅ

，

ｄｉｆｆｅｒｅｎｃｅ ｏｖｅｒ ３０ ｃｍ ｉｎ ｔｈｅ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔ ｏｆ ｔｈｅ ｂｅｎｃｈｍａｒｋｓ ＪＩＮＨ ａｎｄ ＧＥＯＪ ． Ｆｏｒ ｔｈｅ ｄｅｔｅｒｍｉ

， ， ｆｏｒ Ｋｏｒｅａ ａｎｄ ｉｔｓ ａｄｊｏｉｎｉｎｇ ｓｅａｓ ａｓ ａ ｗｈｏｌｅ ． Ｔｈｉｓ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ｗａｓ ｄｅｖｅｌｏｐｅｄ ｂｙ ｕｓｉｎｇ ａｌｌ ａｖａｉｌ ａｂｌｅ ｇｒａｖｉｔｙ ｄａｔａ，ｉｎｃｌｕｄｉｎｇ ｓｕｒｆａｃｅ ａｎｄ ｓａｔｅｌｌｉｔｅ ｄａｔａ ｏｎ ｌａｎｄ ａｎｄ ｏｃｅａｎ ａｒｅａ ． Ｔｈｅ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ｗａｓ ｃｏｍｐｕｔｅｄ ｆｒｏｍ ｔｈｅ ｒｅｄｕｃｅｄ ｇｒａｖｉｔｙ ｄａｔａ ｂｙ ｇｒｉｄｄｉｎｇ ｗｉｔｈ ｌｅａｓｔ ｓｑｕａｒｅｓ ｃｏｌｌｏｃａｔｉｏｎ，２Ｄ ｍｕｌｔｉｂａｎｄ ｓｐｈｅｒｉｃａｌ ＦＦＴ ｔｒａｎｓｆｏｒｍａｔｉｏｎ，ａｎｄ ｒｅｓｔｏｒｅ ｏｆ ｔｈｅ ｔｅｒｒａｉｎ ｅｆｆｅｃｔｓ ａｎｄ ＧＧＭ０２Ｓ ／ ＥＧＭ９６ ｍｏｄｅｌ ｅｆｆｅｃｔｓ ． Ｄａｔａ ｗｅｒｅ ｇｒｉｄｄｅｄ ｉｎ ｔｈｅ ａｒｅａ ３３ ～ ３９° Ｎ，１２４ ～ １３１° Ｅ，ａｔ ａ ｂａｓｉｃ ｇｒｉｄ ｓｐａｃｉｎｇ ｏｆ ０ ． ０１２５° × ０ ． ０１６６６７° （ ０ ． ７５′ × １′）ｉｎ ｌａｔｉｔｕｄｅ ａｎｄ ｌｏｎｇｉｔｕｄｅ ． Ｔｈｅ ｆｉｎａｌ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ｗａｓ ｃｏｍｐｕｔｅｄ ｂｙ ２Ｄ ｓｐｈｅｒｉｃａｌ ＦＦＴ ｗｉｔｈ ４ｂａｎｄｓ ａｎｄ ＷｏｎｇＧｏｒｅ ｍｏｄｉｆｉｃａｔｉｏｎ ｔｏ ｄｅｇｒｅｅ ｂａｎｄ ８０ ～ ９０ ． Ｆｉｎａｌｌｙ，ｔｈｉｓ ｇｒａｖｉｍｅｔｒｉｃ ｍｏｄｅｌ ｗａｓ ｓｕｂｓｅｑｕｅｎｔｌｙ ｆｉｔｔｅｄ ｔｏ ｔｈｅ ＧＰＳ ／ ｌｅｖｅｌｉｎｇ ａｎｄ ｔｉｄｅ ｏｂｓｅｒｖａｔｉｏｎｓ ｂｙ ｌｅａｓｔ ｓｑｕａｒｅ ｃｏｌｌｏｃａｔｉｏｎ，ｔｏ ｐｒｏｖｉｄｅ ｔｈｅ ｆｉｎａｌ ＧＰＳｃｏｎｓｉｓｔｅｎｔ ｌｏｃａｌ ｐｒｅｃｉｓｅ ｇｅｏｉｄ ｍｏｄｅｌ ． Ｆｏｒ ｔｈｅ ｃｏｌ ｌｏｃａｔｉｏｎ ｐｒｏｃｅｓｓ，ｗｅ ｄｅｔｅｒｍｉｎｅｄ ａ ｃｏｒｒｅｌａｔｉｏｎ ｌｅｎｇｔｈ ｏｆ ４０ ｋｍ ｆｏｒ ａ ｓｅｃｏｎｄｏｒｄｅｒ Ｍａｒｋｏｖ ｃｏｖａｒｉａｎｃｅ ｆｕｎｃｔｉｏｎ ａｎｄ １ ｃｍ ａ ｐｒｉｏｒｉ ＧＰＳ ｎｏｉｓｅ ｗａｓ ｕｓｅｄ ． Ｔｈｅ ｐｏｓｔｆｉｔ ｅｒｒｏｒ （ｓｔｄ ． ｄｅｖ ． ）ｏｆ ｔｈｅ ｆｉｎａｌ ｇｅｏｉｄ ｔｏ ｔｈｅ ｎａｔｉｏｎ ｏｆ ａ ｐｒｅｃｉｓｅ ｌｏｃａｌ ｇｅｏｉｄ ｏｆ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ｗｅ ｆｉｒｓｔ ｄｅｖｅｌｏｐｅｄ ａ ｐｒｅｃｉｓｅ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ

Ｎ ＧＰＳ ／ Ｔｉｄｅ ｄｅｒｉｖｅｄ ｆｒｏｍ ＧＰＳ ａｎｄ ｔｉｄｅ ｏｂｓｅｒｖａｔｉｏｎｓ ｗａｓ ± ２ ． ２ ｃｍ ｉｎ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ｏｆ ｔｈｅ Ｇｅｏｇａ

，

Ｇｒａｎｄ Ｂｒｉｄｇｅ ． Ｂａｓｅｄ ｏｎ ｔｈｅ ｆｉｎａｌ ｇｅｏｉｄ ｔｈｅ ｎｅｗ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔ ｗａｓ ｃｏｍｐｕｔｅｄ ａｔ ｔｈｅ ｂｅｎｃｈｍａｒｋｓ ａｎｄ ｃｏｎｔｒｏｌ ｐｏｉｎｔｓ ｉｎ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ． Ｔａｂｌｅ ６ ｓｈｏｗｓ ｔｈｅ ｒｅｓｕｌｔｓ ｏｆ ｔｈｅ ｎｅｗ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔｓ ａｎｄ ｔｈｅ ｄｉｆｆｅｒｅｎｃｅｓ ． Ｔｈｅ ｄｉｆｆｅｒｅｎｃｅｓ ｂｅｔｗｅｅｎ ｔｈｅ ｏｌｄ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔ ｆｏｒ ｔｈｅ ｂｒｉｄｇｅ ｄｅｓｉｇｎ ａｎｄ ｔｈｅ ｎｅｗ ｏｒｔｈｏｍｅｔｒｉｃ

，

，

ｈｅｉｇｈｔ ｆｒｏｍ ｔｈｅ ｆｉｎａｌ ｇｅｏｉｄ ａｒｅ ０ ． ００７ ｍ ａｔ ＪＩＮＨ ｗｈｉｃｈ ｉｓ ｔｈｅ ｍａｉｎｌａｎｄ ｂｅｎｃｈｍａｒｋ ａｎｄ ０ ． ３４２ ｍ ａｔ

，

ＧＥＯＪ ｗｈｉｃｈ ｉｓ ｔｈｅ ｉｓｌａｎｄ ｂｅｎｃｈｍａｒｋ ． Ｔｈｅ ｍｏｄｉｆｉｃａｔｉｏｎ ｏｆ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔｓ ｆｏｒ ｂｒｉｄｇｅ ｄｅｓｉｇｎ ｉｓ ｎｅｃｅｓｓａｒｙ ｆｏｒ ａｃｃｕｒａｔｅ ａｎｄ ｅｆｆｉｃｉｅｎｔ ｃｏｎｓｔｒｕｃｔｉｏｎ ｗｏｒｋ ．

， （），２００９，６７９ － ６９４

ＤｏｎｇＨａ ＬＥＥ ｅｔ ａｌ ． ／ Ｃｈｉｎａ Ｏｃｅａｎ Ｅｎｇｉｎｅｅｒｉｎｇ ２３ ４

Ｔａｂｌｅ ６

Ｎｅｗ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔｓ ｏｆ ｂｅｎｃｈｍａｒｋｓ ａｎｄ ｃｏｎｔｒｏｌ ｐｏｉｎｔｓ ｆｒｏｍ ｔｈｅ ｆｉｎａｌ ｇｅｏｉｄ

Ｏｒｔｈｏｈｅｉｇｈｔ Ｐｏｉｎｔｓ Ｂ． Ｍ．

Ｃｏｎｔｒｏｌ Ｐｏｉｎｔ

Ｂ． Ｍ．

６９３

ＪＩＮＨ

Ｌａｔ．

Ｌｏｎ．

Ｏｒｔｈｏｈｅｉｇｈｔ ｆｏｒ

（Ｄｅｇ． ）

（Ｄｅｇ． ）

ｂｒｉｄｇｅ ｄｅｓｉｇｎ Ａ

３５ ． １０００１

１２８ ． ８１７２８

８ ． ６６７５ ｍ

Ｎｅｗ ｏｒｔｈｏｈｅｉｇｈｔ

Ｄｉｆｆｅｒｅｎｃｅ

８ ． ６７１ ｍ

０ ． ００７ ｍ

（ ） ｆｒｏｍ ｆｉｎａｌ ｇｅｏｉｄ （Ｂ） （Ａ － Ｂ）

ＣＰ１０１

３５ ． ０２８８２

１２８ ． ８１９０２

２４ ． １２１ ｍ

２４ ． ４３８ ｍ

０ ． ３１７ ｍ

ＣＰ１１６

３５ ． ０２１４６

１２８ ． ８０６１２

１５ ． ８７３ ｍ

１６ ． ２１８ ｍ

０ ． ３４５ ｍ

ＣＰ１１０

３５ ． ０１６０４

１２８ ． ７６２４７

９ ． ６２８ ｍ

９ ． ９８９ ｍ

０ ． ３６１ ｍ

ＣＰ１０４

３５ ． ０１５０６

１２８ ． ７１７５８

８０ ． ３１８ ｍ

８０ ． ６６５ ｍ

０ ． ３４７ ｍ

ＣＰ１００

３５ ． ０２１１７

１２８ ． ７１７７０

１２ ． ６４９ ｍ

１２ ． ９７５ ｍ

０ ． ３２６ ｍ

ＧＥＯＪ

３４ ． ８８９４０

１２８ ． ６９０４９

１８ ． ７６３０ ｍ

１９ ． １０５ ｍ

０ ． ３４２ ｍ

，

Ｉｎ ｔｈｅ ｐｒｅｓｅｎｔ ｓｔｕｄｙ ｗｅ ｓｏｌｖｅ ｔｈｅ ｉｎｃｏｎｓｉｓｔｅｎｔ ｈｅｉｇｈｔ ｓｙｓｔｅｍ ｐｒｏｂｌｅｍ ｂｙ ｃａｌｃｕｌａｔｉｎｇ ｔｈｅ ｅｘａｃｔ ｏｒ ｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔｓ ｏｆ ｔｈｅ ｂｅｎｃｈｍａｒｋｓ ａｎｄ ｔｈｅ ｃｏｎｔｒｏｌ ｐｏｉｎｔｓ ｂｙ ｕｓｉｎｇ ｆｉｎａｌ ｇｅｏｉｄ ｍｏｄｅｌ ． Ｔｈｅ ｆｉｎａｌ ｇｅｏｉｄ

，

ｃａｎ ｂｅ ｕｓｅｄ ｔｏ ｇｉｖｅ ｃｏｎｓｉｓｔｅｎｔ ｈｅｉｇｈｔｓ ｉｎ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ． Ａｌｓｏ ｔｈｅ ｐｒｅｃｉｓｅ ｏｒｔｈｏｍｅｔｒｉｃ ｈｅｉｇｈｔ ｉｓ

，

ｅｓｔｉｍａｔｅｄ ｔｈｒｏｕｇｈ ｔｈｅ ＧＰＳ ／ ｌｅｖｅｌｉｎｇ ｔｅｃｈｎｉｑｕｅ ｂｙ ａｐｐｌｙｉｎｇ ｔｈｅ ｆｉｎａｌ ｇｅｏｉｄ ｍｏｄｅｌ ． Ｔｈｅｒｅｆｏｒｅ ｔｈｅ ｐｒｅ

，

， ｐｒｏｃｅｄｕｒｅ ｏｆ ｔｈｅ Ｇｅｏｇａ Ｇｒａｎｄ Ｂｒｉｄｇｅ ． Ｈｏｗｅｖｅｒ，ｔｈｅ ｆｉｎａｌ ｇｅｏｉｄ ｍｏｄｅｌ ｃａｎｎｏｔ ｂｅ ｕｓｅｄ ｄｉｒｅｃｔｌｙ ｆｏｒ ａｎ ａｃｃｕｒａｔｅ ｂｒｉｄｇｅ ｃｏｎｓｔｒｕｃｔｉｏｎ，ｓｕｃｈ ａｓ ｔｈｅ ｐｒｅｃｉｓｉｏｎ ｌｏｆｔ ｏｆ ｂｒｉｄｇｅ ｐｉｅｒ ａｎｄ ｐｉｅｒ ｃａｐ ｅｔｃ ． ，ｂｅｃａｕｓｅ ｉｔ ｈａｓ ａ ｌｏｗ ａｃｃｕｒａｃｙ （± ２ ． ２ ｃｍ）ｗｉｔｈ ｒｅｓｐｅｃｔ ｔｏ ｔｈｅ ａｌｌｏｗａｂｌｅ ａｃｃｕｒａｃｙ ｏｆ ａｃｃｕｒａｔｅ ｂｒｉｄｇｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ． Ｉｎ ｏｒｄｅｒ ｔｏ ｕｓｅ ａ ｌｏｃａｌ ｇｅｏｉｄ ｍｏｄｅｌ ｄｉｒｅｃｔｌｙ ｆｏｒ ｔｈｅ ａｃｃｕｒａｔｅ ｂｒｉｄｇｅ ｃｏｎｓｔｒｕｃｔｉｏｎ，ｉｔ ｉｓ ｎｅｃｅｓｓａｒｙ ｔｏ ｅｎｈａｎｃｅ

ｃｉｓｅ ｌｏｃａｌ ｇｅｏｉｄ ｗｈｉｃｈ ｉｓ ｄｅｖｅｌｏｐｅｄ ｉｎ ｔｈｉｓ ｓｔｕｄｙ ｉｓ ｅｘｐｅｃｔｅｄ ｔｏ ｓｕｃｃｅｓｓ ｔｈｅ ｑｕａｌｉｔｙ ｏｆ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ

ｔｈｅ ａｃｃｕｒａｃｙ ｏｆ ｇｅｏｉｄ ｍｏｄｅｌ ｕｐ ｔｏ ｔｈｅ ｓｕｂｍｍ ｌｅｖｅｌ ｂｙ ｌｏｎｇｔｅｒｍ ｔｉｄａｌ ｏｂｓｅｒｖａｔｉｏｎｓ ａｎｄ ａｄｄｉｔｉｏｎａｌ ｇｒａｖ ｉｔｙ ｏｂｓｅｒｖａｔｉｏｎｓ ｉｎ ｔｈｅ ｃｏｎｓｔｒｕｃｔｉｏｎ ａｒｅａ ｏｆ ｔｈｅ ｂｒｉｄｇｅ ． Ｒｅｆｅｒｅｎｃｅｓ

，

， ， ： ， ， （）： ， ， ， ， ， （ ） ，： ， ， ， ， ， ， ，， ， ， ， ， ， ， ， ， ， ， ， ， （）： ， ， ，， ，， ， ， ， ， ， ， ， （）： ， ， ， ， ， ： ， ， （）： ， ， ， （ ） ， ， ， ， （）： ， ， ， ，： ， ， ， ， ， ，

Ａｍｏｓ Ｍ． Ｊ． ａｎｄ Ｆｅａｔｈｅｒｓｔｏｎｅ Ｗ． Ｅ． ２００９ ． Ｕｎｉｆｉｃａｔｉｏｎ ｏｆ Ｎｅｗ Ｚｅａｌａｎｄ′ｓ ｌｏｃａｌ ｖｅｒｔｉｃａｌ ｄａｔｕｍｓ ｉｔｅｒａｔｉｖｅ ｇｒａｖｉｍｅｔｒｉｃ ｑｕａｓｉｇｅｏｉｄ ｃｏｍｐｕｔａｔｉｏｎｓ Ｊｏｕｒｎａｌ ｏｆ Ｇｅｏｄｅｓｙ ８３ １ ５７ ～ ６８ ． Ａｎｄｅｒｓｏｎ Ｏ． Ｂ． Ｋｎｕｄｓｅｎ Ｐ． ａｎｄ Ｔｒｉｍｍｅｒ Ｒ． ２００５ ． Ｉｍｐｒｏｖｅｄ ｈｉｇｈ ｒｅｓｏｌｕｔｉｏｎ ａｌｔｉｍｅｔｒｉｃ ｇｒａｖｉｔｙ ｆｉｅｌｄ ｍａｐｐｉｎｇ ＫＭＳ２００２ ｇｌｏｂａｌ ｍａｒｉｎｅ ｇｒａｖｉｔｙ ｆｉｅｌｄ ． ＩＡＧ ｓｙｍｐｏｓｉａ ｉｎ Ａ Ｗｉｎｄｏｗ ｏｎ ｔｈｅ Ｆｕｔｕｒｅ ｏｆ Ｇｅｏｄｅｓｙ Ｎｅｗ Ｙｏｒｋ Ｓｐｒｉｎｇｅｒ Ｂｅｒｌｉｎ Ｈｅｉｄｅｌｂｅｒｇ １２８ ３２６ ～ ３３１ ． Ｂｕｒａ Ｍ． Ｋｏｕｂａ Ｊ． Ｋｕｍａｒ Ｍ． Ｍüｌｌｅｒ Ａ． Ｒａｄěｊ Ｋ． Ｔｒｕｅ Ｓ． Ａ． Ｖａｔｒｔ Ｖ． ａｎｄ Ｖｏｊｔíｋｏｖá Ｍ． １９９９ ． Ｇｅｏｉｄａｌ ｇｅｏｐｏｔｅｎｔｉａｌ ａｎｄ ｗｏｒｌｄ ｈｅｉｇｈｔ ｓｙｓｔｅｍ Ｓｔｕｄ ． Ｇｅｏｐｈｙｓ ． Ｇｅｏｄ ． ４３ ４ ３２７ ～ ３３７ ． Ｂｕｒａ Ｍ． Ｋｅｎｙｏｂ Ｓ． Ｋｏｕｂａ Ｊ． íｍａ Ｚ． Ｖａｔｒｔ Ｖ． ａｎｄ Ｖｏｊｔíｋｏｖá Ｍ． ２００４ ． Ａ ｇｌｏｂａｌ ｖｅｒｔｉｃａｌ ｒｅｆｅｒｅｎｃｅ ｆｒａｍｅ ｂａｓｅｄ ｏｎ ｆｏｕｒ ｒｅｇｉｏｎａｌ ｖｅｒｔｉｃａｌ ｄａｔｕｍｓ Ｓｔｕｄ ． Ｇｅｏｐｈｙｓ ． Ｇｅｏｄ ． ４８ ３ ４９３ ～ ５０２ ． Ｃｏｒｃｈｅｔｅ Ｖ． Ｃｈｏｕｒａｋ Ｍ． ａｎｄ Ｋｈａｔｔａｃｈ Ｄ． ２００５ ． Ｔｈｅ ｈｉｇｈｒｅｓｏｌｕｔｉｏｎ ｇｒａｖｉｍｅｔｒｉｃ ｇｅｏｉｄ ｏｆ Ｉｂｅｒｉａ ＩＧＧ２００５ Ｇｅｏ ｐｈｙｓ ． Ｊ ． Ｉｎｔ ． １６２ ７ ６７６ ～ ６８４ ． Ｄａｅｗｏｏ ｅｎｇｉｎｅｅｒｉｎｇ ａｎｄ ｃｏｎｓｔｒｕｃｔｉｏｎ Ｃｏ． ｌｔｄ． ２００３ ． Ｉｎｓｐｅｃｔｉｏｎ ｓｕｒｖｅｙ ｒｅｐｏｒｔ ｆｏｒ ＰｕｓａｎＧｅｏｊｅ ｆｉｘｅｄ ｌｉｎｋ ｐｒｏｊｅｃｔ Ｔｅｃｈｎｉ ｃａｌ Ｒｅｐｏｒｔｓ． ｉｎ Ｋｏｒｅａｎ Ｅｋｍａｎ Ｍ． １９９９ ． Ｕｓｉｎｇ ｍｅａｎ ｓｅａ ｓｕｒｆａｃｅ ｔｏｐｏｇｒａｐｈｙ ｆｏｒ ｔｈｅ ｄｅｔｅｒｍｉｎａｔｉｏｎ ｏｆ ｈｅｉｇｈｔ ｓｙｓｔｅｍ ｄｉｆｆｅｒｅｎｃｅｓ ａｃｒｏｓｓ ｔｈｅ Ｂａｌｔｉｃ ｓｅａ Ｍａｒｉｎｅ Ｇｅｏｄｅｓｙ ２２ １ ３１ ～ ３５ ． Ｆｅｉ Ｚ． Ｌ． ａｎｄ Ｓｉｄｅｒｉｓ Ｍ． Ｇ． ２００１ ． ＧＰＳ ｌｅｖｅｌｉｎｇ ａｎｄ ｔｈｅ ｓｅｃｏｎｄ ｇｅｏｄｅｔｉｃ ｂｏｕｎｄａｒｙ ｖａｌｕｅ ｐｒｏｂｌｅｍ ｉｎ Ｇｒａｖｉｔｙ Ｇｅｏｉｄ ａｎｄ Ｇｅｏｄｙｎａｍｉｃｓ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ａｓｓｏｃｉａｔｉｏｎ ｏｆ Ｇｅｏｄｅｓｙ Ｓｙｍｐｏｓｉａ Ｖｏｌ． １２３ Ｂｅｒｌｉｎ ＳｐｒｉｎｇｅｒＶｅｒｌａｇ ３４１ ～ ３４６ ．

６９４

， （），２００９，６７９ － ６９４

ＤｏｎｇＨａ ＬＥＥ ｅｔ ａｌ ． ／ Ｃｈｉｎａ Ｏｃｅａｎ Ｅｎｇｉｎｅｅｒｉｎｇ ２３ ４

， ， ， ， （ ）： ， ， ， ， ， ， ， （）： ， ， ， ，， ， ， ， ， ， ， ， ， ， ， ， ， ， ， （ ）： ， ， ， ， ， ， （）： ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， ， （）： （ ） ， ， ， ， ， ， ， ， ， （）： ， ， ， ， ， ， （）： ， ， ，： ， ， ， ， ， ， ， ， ， ， ， （）： ， ， ，， ，， ， ， ， ， ，， ， ， ， ， ， ， ， ， ， （）： ， ， ， ， ， ， ， ， ， ， ， （ ）： ， ， ， ， （）： ， ， ， ， ， ， ， （）： ， ，， ， ， ， ， ，

，

Ｆｏｒｓｂｅｒｇ Ｒ． ａｎｄ Ｔｓｃｈｅｒｎｉｎｇ Ｃ． Ｃ． １９８１ ． Ｔｈｅ ｕｓｅ ｏｆ ｈｅｉｇｈｔ ｄａｔａ ｉｎ ｇｒａｖｉｔｙ ｆｉｅｌｄ ａｐｐｒｏｘｉｍａｔｉｏｎ ｂｙ ｃｏｌｌｏｃａｔｉｏｎ Ｊ ． Ｇｅｏｐｈｙｓ ． Ｒｅｓ ． ８６ Ｂ９ ７８４３ ～ ７８５４ ． Ｆｏｒｓｂｅｒｇ Ｒ． １９８５ ． Ｇｒａｖｉｔｙ ｆｉｅｌｄ ｔｅｒｒａｉｎ ｅｆｆｅｃｔ ｃｏｍｐｕｔａｔｉｏｎｓ ｂｙ ＦＦＴ Ｊｏｕｒｎａｌ ｏｆ Ｇｅｏｄｅｓｙ ５９ ４ ３４２ ～ ３６０ ． Ｆｏｒｓｂｅｒｇ Ｒ． ａｎｄ Ｓｉｄｅｒｉｓ Ｍ． Ｇ． １９９３ ． Ｇｅｏｉｄ ｃｏｍｐｕｔａｔｉｏｎｓ ｂｙ ｔｈｅ ｍｕｌｔｉｂａｎｄ ｓｐｈｅｒｉｃａｌ ＦＦＴ ａｐｐｒｏａｃｈ Ｍａｎｕｓｃｒｉｐｔａ Ｇｅｏｄａｅｔｉｃａ １８ ２ ８２ ～ ９０ ． Ｆｏｒｓｂｅｒｇ Ｒ． Ｔｓｃｈｅｒｎｉｎｇ Ｃ． Ｃ． ａｎｄ Ｋｎｕｄｓｅｎ Ｐ． ２００３ ． Ａｎ ｏｖｅｒｖｉｅｗ ｍａｎｕａｌ ｆｏｒ ｔｈｅ ＧＲＡＶＳＯＦＴ Ｋｏｒｔ ＆ Ｍａｔｒｉｋｅｌ ｓｔｙｒｅｌｓｅｎ． Ｈｅｉｓｋａｎｅｎ Ｗ． Ａ． ａｎｄ Ｍｏｒｉｔｚ Ｈ． １９６７ ． Ｐｈｙｓｉｃａｌ Ｇｅｏｄｅｓｙ Ｓａｎ Ｆｒａｎｓｉｓｃｏ Ｗ． Ｈ． Ｆｒｅｅｍａｎ ａｎｄ Ｃｏ． Ｉｌｉｆｆｅ Ｊ． Ｃ． Ｚｉｅｂａｒｔ Ｍ． Ｃｒｏｓｓ Ｐ． Ａ． Ｆｏｒｓｂｅｒｇ Ｒ． Ｓｔｒｙｋｏｗｓｋｉ Ｇ． ａｎｄ Ｔｓｃｈｅｒｎｉｎｇ Ｃ． Ｃ． ２００３ ． ＯＳＧＭ０２ ａ ｎｅｗ ｍｏｄｅｌ ｆｏｒ ｃｏｎｖｅｒｔｉｎｇ ＧＰＳｄｅｒｉｖｅｄ ｈｅｉｇｈｔｓ ｔｏ ｌｏｃａｌ ｈｅｉｇｈｔ ｄａｔｕｍｓ ｉｎ Ｇｒｅａｔ Ｂｒｉｔａｉｎ ａｎｄ Ｉｒｅｌａｎｄ Ｓｕｒｖｅｙ Ｒｅｖｉｅｗ ３７ ２９０ ２７６ ～ ２９３ ． Ｋｕｒｏｉｓｈｉ Ｙ． Ａｎｄｏ Ｈ． ａｎｄ Ｆｕｋｕｄａ Ｙ． ２００２ ． Ａ ｎｅｗ ｈｙｂｒｉｄ ｇｅｏｉｄ ｍｏｄｅｌ ｆｏｒ Ｊａｐａｎ ＧＳＩＧＥＯ２０００ Ｊｏｕｒｎａｌ ｏｆ Ｇｅｏｄｅｓｙ ７６ ８ ４２８ ～ ４３６ ． Ｌｅｍｏｉｎｅ Ｆ． Ｇ． Ｓｍｉｔｈ Ｄ． Ｓ． Ｔｏｒｒｅｎｃｅ Ｍ． Ｈ． Ｗｉｌｌｉａｍｓｏｎ Ｒ． Ｇ． Ｃｏｘ Ｃ． Ｍ． Ｒａｃｈｌｉｎ Ｋ． Ｅ． Ｍｉｎｇ Ｙ． Ｍ． Ｋｅｎｙｏｎ Ｓ． Ｃ． Ｓａｌｍａｎ Ｒ． Ｔｒｉｍｍｅｒ Ｒ． Ｒａｐｐ Ｒ． Ｈ． ａｎｄ Ｎｅｒｅｍ Ｒ． Ｓ． １９９６ ． Ｔｈｅ ｄｅｖｅｌｏｐｅｍｅｎｔ ｏｆ ｔｈｅ ＮＡＳＡ ＧＳＦＣ ａｎｄ ＮＩＭＡ ｊｏｉｎｔ ｇｅｏｐｏｔｅｎｔｉａｌ ｍｏｄｅｌ Ｉｎｔ ． Ｓｙｍｐｏｓｉｕｍ ｏｎ Ｇｒａｖｉｔｙ Ｇｅｏｉｄ ａｎｄ Ｍａｒｉｎｅ Ｇｅｏｄｅｓｙ Ｔｏｋｙｏ ４６１ ～ ４６９ ． Ｌｉ Ｆ． Ｃｈｅｎ Ｗ． ａｎｄ Ｙｕｅ Ｊ． Ｌ． ２００３ ． Ｏｎ ｓｏｌｕｔｉｏｎ ａｎｄ ａｐｐｌｉｃａｔｉｏｎ ｏｆ ＧＰＳ ／ Ｇｒａｖｉｔｙ ｂｏｕｎｄａｒｙ ｖａｌｕｅ ｐｒｏｂｌｅｍ Ｃｈｉ ｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ｇｅｏｐｈｙｓｉｃｓ ４６ ５ ５９５ ～ ５９９ ． ｉｎ Ｃｈｉｎｅｓｅ Ｍｏｒｉｔｚ Ｈ． １９８０ ． Ａｄｖａｎｃｅｄ Ｐｈｙｓｉｃａｌ Ｇｅｏｄｅｓｙ Ｈｅｒｂｅｒｔ Ｗｉｃｈｍａｎ Ｖｅｒｌａｇ Ｋａｒｌｓｒｕｈｅ ＆ Ａｂａｃｕｓ Ｐｒｅｓｓ Ｔｕｎｂｒｉｄｇｅ Ｗｅｌｌｓ Ｋｅｎｔ． Ｏｍａｎｇ Ｏ． Ｃ． Ｄ． ａｎｄ Ｆｏｒｓｂｅｒｇ Ｒ． ２０００ ． Ｈｏｗ ｔｏ ｈａｎｄｌｅ ｔｏｐｏｇｒａｐｈｙ ｉｎ ｐｒａｃｔｉｃａｌ ｇｅｏｉｄ ｄｅｔｅｒｍｉｎａｔｉｏｎ ｔｈｒｅｅ ｅｘａｍ ｐｌｅｓ Ｊｏｕｒｎａｌ ｏｆ Ｇｅｏｄｅｓｙ ７４ ６ ４５８ ～ ４６６ ． Ｒｏｈ Ｔ． Ｈ． Ｓｅｏ Ｄ． Ｊ． ａｎｄ Ｌｅｅ Ｊ． Ｃ． ２００１ ． Ａｎａｌｙｓｉｓ ｏｆ ｉｎｃｉｄｅｎｃｅ ａｎｇｌｅ ｕｓｉｎｇ ｔｏｔａｌ ｓｔａｔｉｏｎ ｉｎ ｌｅｖｅｌｉｎｇ Ｋｏｒｅａｎ Ｊｏｕｒｎａｌ ｏｆ Ｇｅｏｍａｔｉｃｓ １ １ ２１ ～ ２６ ． Ｒｕｍｍｅｌ Ｒ． ２０００ ． Ｇｌｏｂａｌ ｕｎｉｆｉｃａｔｉｏｎ ｏｆ ｈｅｉｇｈｔ ｓｙｓｔｅｍ ａｎｄ ＧＯＣＥ ｉｎ Ｇｒａｖｉｔｙ Ｇｅｏｉｄ ａｎｄ Ｇｅｏｄｙｎａｍｉｃｓ Ｉｎｔｅｒｎａｔｉｏｎ ａｌ Ａｓｓｏｃｉａｔｉｏｎ ｏｆ Ｇｅｏｄｅｓｙ Ｓｙｍｐｏｓｉａ Ｖｏｌ． １２３ Ｂｅｒｌｉｎ ＳｐｒｉｎｇｅｒＶｅｒｌａｇ １５ ～ ２０ ． Ｓｃｈｗａｒｚ Ｋ． Ｐ． Ｓｉｄｅｒｉｓ Ｍ． Ｇ． ａｎｄ Ｆｏｒｓｂｅｒｇ Ｒ． １９９０ ． Ｔｈｅ ｕｓｅ ｏｆ ＦＦＴ ｔｅｃｈｎｉｑｕｅｓ ｉｎ ｐｈｙｓｉｃａｌ ｇｅｏｄｅｓｙ Ｇｅｏｐｈｙｓｉ ｃａｌ Ｊｏｕｒｎａｌ Ｉｎｔｅｒｎａｔｉｏｎａｌ １００ ３ ４８５ ～ ５１４ ． Ｔａｐｌｅｙ Ｂ． Ｒｉｅｓ Ｊ． Ｂｅｔｔａｄｐｕｒ Ｓ． Ｃｈａｍｂｅｒｓ Ｄ． Ｃｈｅｎｇ Ｍ． Ｃｏｎｄｉ Ｆ． Ｇｕｎｔｅｒ Ｂ． Ｋａｎｇ Ｚ． Ｎａｇｅｌ Ｐ． Ｐａｓｔｏｒ Ｒ． Ｐｅｋｋｅｒ Ｔ． Ｐｏｏｌｅ Ｓ． ａｎｄ Ｗａｎｇ Ｆ． ２００５ ． ＧＧＭ０２Ａｎ ｉｍｐｒｏｖｅｄ Ｅａｒｔｈ ｇｒａｖｉｔｙ ｆｉｅｌｄ ｍｏｄｅｌ ｆｒｏｍ ＧＲＡＣＥ Ｊｏｕｒｎａｌ ｏｆ Ｇｅｏｄｅｓｙ ７９ ８ ４６７ ～ ４７８ ． Ｔｓｃｈｅｒｎｉｎｇ Ｃ． Ｃ． Ｋｎｕｄｓｅｎ Ｐ． ａｎｄ Ｆｏｒｓｂｅｒｇ Ｒ． １９９４ ． Ｄｅｓｃｒｉｐｔｉｏｎ ｏｆ ｔｈｅ ＧＲＡＶＳＯＦＴ ｐａｃｋａｇｅ Ｔｅｃｈｎｉｃａｌ Ｒｅｐｏｒｔ Ｇｅｏｐｈｙｓｉｃａｌ Ｉｎｓｔｉｔｕｔｅ Ｕｎｉｖｅｒｓｉｔｙ ｏｆ Ｃｏｐｅｎｈａｇｅｎ． Ｖａｎíˇｃｅｋ Ｐ． ａｎｄ Ｆｅａｔｈｅｒｓｔｏｎｅ Ｗ． Ｅ． １９９８ ． Ｐｅｒｆｏｒｍａｎｃｅ ｏｆ ｔｈｒｅｅ ｔｙｐｅｓ ｏｆ Ｓｔｏｋｅｓ′ｓ ｋｅｒｎｅｌ ｉｎ ｔｈｅ ｃｏｍｂｉｎｅｄ ｓｏｌｕｔｉｏｎ ｆｏｒ ｔｈｅ ｇｅｏｉｄ Ｊｏｕｒｎａｌ ｏｆ Ｇｅｏｄｅｓｙ ７２ １２ ６８４ ～ ６９７ ． Ｗｏｎｇ Ｌ． ａｎｄ Ｇｏｒｅ Ｒ． １９６９ ． Ａｃｃｕｒａｃｙ ｏｆ ｇｅｏｉｄ ｈｅｉｇｈｔｓ ｆｒｏｍ ｍｏｄｉｆｉｅｄ Ｓｔｏｋｅｓ ｋｅｒｎｅｌｓ Ｇｅｏｐｈｙｓｉｃａｌ Ｊｏｕｒｎａｌ ｏｆ ｔｈｅ Ｒｏｙａｌ Ａｓｔｒｏｎｏｍｉｃａｌ Ｓｏｃｉｅｔｙ １８ １ ８１ ～ ９１ ． Ｙｕｎ Ｈ． Ｓ． １９９５ ． Ｒｅｓｕｌｔｓ ｏｆ ｔｈｅ ｇｅｏｉｄ ｃｏｍｐｕｔａｔｉｏｎ ｆｏｒ Ｋｏｒｅａｎ ｐｅｎｉｎｓｕｌａ Ｐｈ． Ｄ Ｔｈｅｓｉｓ Ｔｅｃｈｎｉｃａｌ Ｕｎｉｖｅｒｓｉｔｙ ｏｆ Ｂｕ ｄａｐｅｓｔ Ｈｕｎｇａｒｙ． Ｙｕｎ Ｈ． Ｓ． １９９９ ． Ｐｒｅｃｉｓｉｏｎ ｇｅｏｉｄ ｄｅｔｅｒｍｉｎａｔｉｏｎ ｂｙ ｓｐｈｅｒｉｃａｌ ＦＦＴ ｉｎ ａｎｄ ａｒｏｕｎｄ ｔｈｅ Ｋｏｒｅａｎ ｐｅｎｉｎｓｕｌａ Ｅａｒｔｈ Ｐｌａｎｅｔｓ Ｓｐａｃｅ ５１ １ １３ ～ １８ ． Ｚｈａｎｇ Ｌ． Ｌｉ Ｆ． Ｃｈｅｎ Ｗ． ａｎｄ Ｚｈａｎｇ Ｃ． ２００９ ． Ｈｅｉｇｈｔ ｄａｔｕｍ ｕｎｉｆｉｃａｔｉｏｎ ｂｅｔｗｅｅｎ Ｓｈｅｎｚｈｅｎ ａｎｄ Ｈｏｎｇ Ｋｏｎｇ ｕｓ ｉｎｇ ｔｈｅ ｓｏｌｕｔｉｏｎ ｏｆ ｔｈｅ ｌｉｎｅａｒｉｚｅｄ ｆｉｘｅｄｇｒａｖｉｍｅｔｒｉｃ ｂｏｕｎｄａｒｙ ｖａｌｕｅ ｐｒｏｂｌｅｍ Ｊｏｕｒｎａｌ ｏｆ Ｇｅｏｄｅｓｙ ８３ ５ ４１１ ～ ４１７ ．

， （）：

，

，

，

， ，

： ，

， ， ， ， ， ， ， ， ， ， ：

，

， ， ， ， ，， ，， ，

，

， ，

，

， （）：