Abstract
The free-air anomalies are computed using a data set from various types of gravity measurements in the Korean Peninsula area. The gravity values extracted from the Earth Gravitational Model 2008 are used in the surrounding region. The upward continuation technique suggested by Dragomir is used in the computation of the external free-air anomalies at various altitudes. The integration radius 10 times the altitude is used in order to keep the accuracy of results and computational resources. The direct geodesic formula developed by Bowring is employed in integration. At the 1-km altitude, the free-air anomalies vary from -41.315 to 189.327 mgal with the standard deviation of 22.612 mgal. At the 3-km altitude, they vary from -36.478 to 156.209 mgal with the standard deviation of 20.641 mgal. At the 1,000-km altitude, they vary from 3.170 to 5.864 mgal with the standard deviation of 0.670 mgal. The predicted free-air anomalies at 3-km altitude are compared to the published free-air anomalies reduced from the airborne gravity measurements at the same altitude. The rms difference is 3.88 mgal. Considering the reported 2.21-mgal airborne gravity cross-over accuracy, this rms difference is not serious. Possible causes in the difference appear to be external free-air anomaly simulation errors in this work and/or the gravity reduction errors of the other. The external gravity field is predicted by adding the external free-air anomaly to the normal gravity computed using the closed form formula for the gravity above and below the surface of the ellipsoid. The predicted external gravity field in this work is expected to reasonably present the real external gravity field. This work seems to be the first structured research on the external free-air anomaly in the Korean Peninsula area, and the external gravity field can be used to improve the accuracy of the inertial navigation system.
주변 해역을 포함한 한반도 일원에서 측정된 중력자료로부터 상층중력의 고도이상(free-air anomaly)을 계산하였다. 주변 영역에서는 인접국가가 발표한 중력자료가 있는 경우 발표된 자료를 이용하였으며, 없는 경우 EGM2008(Earth Gravitational Model 2008)로부터 계산한 고도이상을 이용하였다. 중력의 상향연속은 Dragomir가 제안한 방법으로 계산하였다. 상층중력 고도이상 계산의 정확성과 계산 속도를 고려하여 적분반경은 계산 고도의 10배로 하였다. 적분에 필요한 측지선의 거리는 Bowring이 개발한 공식을 사용하였다. 위도 $33^{\circ}N{\sim}43^{\circ}N$, 경도 $124^{\circ}E{\sim}131^{\circ}E$에서 계산된 고도이상은 고도 1 km에서 -41.315에서 189.327 mgal까지 변화하고 표준 편차는 22.612 mgal이다. 고도 3 km에서는 -36.478에서 156.209 mgal까지 변화하고 표준 편차는 20.641 mgal이다. 고도 1,000 km에서는 3.170에서 5.864 mgal까지 변화하고 표준 편차는 0.670 mgal이다. 3 km 고도에서 계산된 고도이상을 같은 높이에서 측정한 항공 중력 고도이상과 비교하였다. 이들의 rms 오차는 3.88 mgal로 나타났다. 항공 중력 측정 교차점오차가 2.2 mgal 임을 고려하면 이들 오차에 의미를 부여할 수 없으며, 원인으로는 이번 연구에서 발생한 계산상 오차와 함께/또는 발표된 항공중력의 보정오차에 기인하는 것으로 사료된다. 상층중력 고도이상에 완전식으로 계산한 지구타원체 외부의 정규중력을 더하여 상층중력을 예측하였다. 이번 연구에서 국내 최초로 계산한 고도에 따른 상층중력 고도이상은 한반도 일원의 상층중력장을 잘 표현하고 있는 것으로 보이며, 상층중력장은 관성항법장치의 정확도 향상 등에 이용될 수 있을 것이다.