지구 형상(Figure of the Earth)을 정의하는 문제는 측지학(Geodesy)의 가장 근본적인 과제 중 하나이다. 인류는 역사적으로 지구의 모양을 평면, 구, 그리고 자전에 의한 편평도를 반영한 지구 타원체(Earth Ellipsoid)로 발전시켜 이해해 왔다. 그러나 지구 내부의 밀도 분포가 균일하지 않고 지형의 기복이 존재함에 따라, 단순한 기하학적 모델인 타원체만으로는 실제 지구의 물리적 특성을 완벽히 설명할 수 없다. 이러한 배경에서 도입된 개념이 바로 지오이드(Geoid)이다. 지오이드는 지구의 중력장(Gravity field)을 기반으로 정의되는 물리적 가상면으로, 지구 전체의 질량 분포와 자전의 영향을 반영한 고유한 형상을 갖는다.

물리적 관점에서 지오이드는 지구 중력장의 등전위면(Equipotential surface) 중 하나로 정의된다. 지구상의 한 점에 작용하는 중력은 질량 간의 만유인력과 지구 자전에 의한 원심력의 벡터 합으로 나타난다. 이때 중력의 방향을 나타내는 연직선(Plumb line)에 모든 지점에서 수직인 곡면을 등전위면이라 하며, 이 중 평균 해수면(Mean Sea Level, MSL)과 가장 잘 일치하는 특정 등전위면을 지오이드라고 부른다. 이론적으로 해수가 육지 내부까지 자유롭게 흐를 수 있는 가상의 관을 가정했을 때, 해류나 조석, 기압 변화 등 외부 동적 요인이 없는 정지 상태의 해수면이 이루는 면이 곧 지오이드가 된다¹⁾.

지오이드를 수학적으로 정의하기 위해 중력 포텐셜(Gravity potential) $ W $를 도입할 수 있다. 임의의 공간 좌표 $ (x, y, z) $에서 중력 포텐셜 $ W $는 질량에 의한 인력 포텐셜 $ V $와 자전에 의한 원심력 포텐셜 $ $의 합으로 표현된다.

$$ W(x, y, z) = V(x, y, z) + \Phi(x, y, z) $$

이때 지오이드는 $ W(x, y, z) = W_0 $라는 일정한 상수 값을 만족하는 집합으로 정의된다. 여기서 $ W_0 $는 전 지구적인 평균 해수면에서의 포텐셜 값을 의미한다. 등전위면으로서의 지오이드는 그 표면 위의 모든 지점에서 중력 포텐셜이 동일하므로, 지오이드 면을 따라 질량을 이동시킬 때 중력이 하는 일의 양은 0이 된다. 이러한 특성 때문에 지오이드는 물리적 고도(Altitude) 측정의 절대적인 기준면인 수준면(Level surface)의 역할을 수행한다.

실제 지구에서 지오이드는 지구 타원체와 일치하지 않고 복잡한 굴곡을 형성한다. 이는 지구 내부의 밀도가 불균일하여 지역마다 중력의 크기와 방향이 달라지기 때문이다. 예를 들어, 지각 아래에 고밀도의 물질이 집중된 지역에서는 중력이 강하게 작용하여 등전위면이 타원체 밖으로 부풀어 오르고, 반대로 밀도가 낮은 지역에서는 안으로 함몰되는 형태를 띤다. 지오이드와 기준 타원체 사이의 수직 거리를 지오이드고(Geoid height) 또는 지오이드 기복(Geoid undulation)이라 하며, 이는 전 지구적으로 약 -106m에서 +85m 사이의 범위를 갖는 것으로 알려져 있다.

또한, 지오이드는 실제의 평균 해수면과도 미세한 차이를 보인다. 실제 해수는 해류, 염분 농도 차이에 의한 밀도 변화, 기압 배치, 해수면 온도 등 역학적인 요인에 의해 끊임없이 움직이기 때문에, 정지 상태를 가정하는 지오이드 면과 완전히 일치할 수 없다. 이러한 지오이드와 실제 평균 해수면 사이의 고도 차이를 해면 지형(Sea Surface Topography)이라 하며, 이는 해양학에서 해류의 순환을 분석하는 중요한 지표가 된다. 결과적으로 지오이드는 지구의 물리적 형상을 가장 정밀하게 대변하는 기준면으로서, 정밀 측량, 위성 항법, 지구 내부 구조 연구 등 현대 지구과학 전반의 기초를 형성한다.

지구의 물리적 형상을 정의하는 핵심 개념인 지오이드는 중력 포텐셜이 일정한 기하학적 궤적, 즉 등전위면(Equipotential surface)의 일종이다. 지구의 중력장(Gravity field) 내에서 임의의 점에 작용하는 중력(Gravity)은 질량 간의 만유인력(Gravitation)과 지구 자전에 의한 원심력(Centrifugal force)의 벡터 합으로 결정된다. 이러한 중력장의 크기와 방향을 수학적으로 기술하기 위해 도입된 중력 포텐셜(Gravity potential) $ W $는 공간상의 좌표 $ (x, y, z) $의 함수로 표현되며, 지오이드는 이 포텐셜 값이 특정한 상수 $ W_0 $를 유지하는 면으로 정의된다. 물리적 관점에서 지오이드는 평균 해수면(Mean Sea Level, MSL)을 육지 내부까지 연장하였을 때 형성되는 가상의 폐곡면을 의미하며, 이는 지구 시스템 내에서 위치 에너지의 기준이 된다.

중력 포텐셜 $ W $는 다음과 같이 만유인력 포텐셜 $ V $와 원심력 포텐셜 $ $의 합으로 나타낼 수 있다.

$$ W(x, y, z) = V(x, y, z) + \Phi(x, y, z) = \int_M \frac{G \, dm}{r} + \frac{1}{2} \omega^2 (x^2 + y^2) $$

여기서 $ G $는 중력 상수이며, $ r $은 질량 요소 $ dm $으로부터의 거리, $ $는 지구 자전의 각속도를 나타낸다. 지오이드 면 위에서는 포텐셜의 변화량 $ dW $가 0이 되므로, 해당 면을 따라 질량을 이동시킬 때 중력이 수행하는 일(Work)은 0이다. 이러한 성질은 지오이드가 물리적인 고도 체계의 기준면이 되는 근거를 제공한다. 즉, 지오이드 상의 모든 지점은 동일한 중력 위치 에너지를 가지며, 유체가 정역학적 평형 상태에 있을 때 그 표면은 반드시 지오이드와 일치하게 된다.

중력 등전위면으로서 지오이드가 갖는 가장 중요한 기하학적 특성은 연직선(Plumb line)과의 직교성이다. 중력 가속도 벡터 $ $는 중력 포텐셜의 기울기(Gradient)로 정의되며, 수학적으로 등전위면의 법선(Normal) 방향을 향한다.

$$ \mathbf{g} = \nabla W = \left( \frac{\partial W}{\partial x}, \frac{\partial W}{\partial y}, \frac{\partial W}{\partial z} \right) $$

따라서 지오이드 면 위의 모든 점에서 추를 매단 실이 가리키는 방향인 연직선은 지오이드 면과 수직을 이룬다. 이는 측량 현장에서 사용하는 수준기(Level)가 물리적으로는 지오이드의 접평면을 결정하는 장치임을 의미한다. 반면, 기하학적인 모델인 기준 타원체(Reference ellipsoid)의 법선 방향은 실제 중력 방향인 연직선과 일치하지 않는 경우가 많으며, 이 두 방향 사이의 각도 차이를 수직선 편차(Deflection of the vertical)라고 한다.

지오이드의 형상은 지구 내부의 질량 분포와 밀도 불균형에 의해 결정된다. 지각 내부에 고밀도 물질이 매장되어 있거나 해저 산맥과 같은 거대 지형이 존재하는 경우, 해당 지역에서는 추가적인 인력이 발생하여 중력 포텐셜이 국지적으로 상승한다. 이로 인해 등전위면인 지오이드는 기준 타원체보다 위로 볼록하게 솟아오르는 지오이드고(Geoid height)의 양(+)의 값을 갖게 된다. 반대로 밀도가 낮은 지역이나 해구와 같은 지형에서는 지오이드가 타원체 아래로 함몰되는 음(-)의 값을 보인다. 이러한 지오이드의 굴곡은 지구 내부의 맨틀 대류나 지각 평형(Isostasy) 상태를 연구하는 데 필수적인 정보를 제공한다.

지오이드와 기준 타원체의 물리적 및 기하학적 차이를 정리하면 다음과 같다.

결과적으로 지오이드는 지구의 실제 중력장을 반영하는 유일한 물리적 기준면이다. 위성 항법 시스템(Global Navigation Satellite System, GNSS)을 통해 얻는 높이 정보는 기하학적인 타원체고이므로, 이를 실질적인 물의 흐름이나 높이 차이를 나타내는 물리적 고도로 변환하기 위해서는 해당 지역의 정확한 지오이드 모델이 반드시 확보되어야 한다²⁾. 지오이드의 정밀한 결정은 측지학뿐만 아니라 해양학, 지구물리학 분야에서 지구 시스템의 역학적 거동을 이해하는 기초가 된다³⁾.

지오이드는 이론적으로 정지 상태에 있는 해수면과 일치하는 중력 등전위면으로 정의된다. 만약 지구상에 해류, 조석, 기압 차이 등의 외부 교란 요인이 전혀 존재하지 않는다면, 해수는 오직 중력과 원심력의 합력인 중력 포텐셜의 영향만을 받아 평형 상태를 이룰 것이다. 이때 형성되는 가상의 해수면이 바로 지오이드이다. 그러나 실제 지구의 바다는 끊임없이 움직이며, 물리적·화학적 성질의 불균일성으로 인해 지오이드와 완벽히 일치하지 않는다. 이러한 실제 평균 해수면(Mean Sea Level, MSL)과 지오이드 사이의 수직적 차이를 해면 지형(Ocean Topography) 또는 해양 동역학적 지형이라 부른다.

평균 해수면이 지오이드로부터 이탈하게 만드는 주요 원인은 해양 동역학적 요인에 있다. 해류의 흐름, 해수의 온도 및 염분 분포에 따른 밀도 차이, 그리고 대기압의 공간적 불균형 등이 해수면의 형상을 변화시킨다. 특히 대규모 해류가 흐르는 해역에서는 코리올리 효과(Coriolis effect)에 의한 수평적 압력 구배력이 발생하며, 이는 해수면의 경사를 유도하여 지형류(Geostrophic current) 평형 상태를 유지하게 한다. 또한, 수온이 높아 열팽창이 일어난 해역은 상대적으로 해수면이 높게 형성되는 반면, 차갑고 밀도가 높은 해역은 해수면이 낮아지는 경향을 보인다. 이러한 요인들로 인해 실제 평균 해수면은 지오이드로부터 전 지구적으로 약 -2m에서 +2m 사이의 편차를 나타낸다.

수학적으로 실제 해수면의 높이인 해면 고도(Sea Surface Height, SSH)와 지오이드고(Geoid height), 그리고 평균 해면 지형(Mean Dynamic Topography, MDT) 사이의 관계는 다음과 같은 선형 결합으로 표현할 수 있다⁴⁾.

$$ \zeta = h - N $$

여기서 $ $는 평균 해면 지형, $ h $는 기준 타원체(Reference ellipsoid)로부터 측정한 평균 해수면의 높이, $ N $은 타원체로부터 지오이드까지의 거리인 지오이드고를 의미한다. 현대 측지학에서는 위성 고도계(Satellite altimetry)를 통해 $ h $를 정밀하게 측정하고, 이를 그레이스 위성(GRACE)이나 고체 위성(GOCE) 등 중력 탐사 위성으로부터 얻어진 지오이드 모델 $ N $과 비교함으로써 전 지구적인 평균 해면 지형의 분포를 도출한다⁵⁾.

이러한 지오이드와 평균 해수면의 상관관계 분석은 지구물리학과 해양학의 접점에서 매우 중요한 학술적 가치를 지닌다. 지오이드를 기준으로 한 해수면의 미세한 굴곡을 파악함으로써 해류의 순환 구조와 수송량을 정량적으로 계산할 수 있으며, 이는 지구 전체의 에너지 순환을 이해하는 핵심 지표가 된다. 또한, 장기적인 해수면 상승과 같은 기후 변화 양상을 모니터링하기 위해서는 해수면의 변화가 순수한 해수량의 증가에 의한 것인지, 아니면 해양 내부의 동역학적 변화에 의한 것인지를 구분해야 하므로 지오이드와의 정밀한 비교가 반드시 선행되어야 한다⁶⁾. 결과적으로 지오이드는 실제 해수면의 거동을 해석하기 위한 절대적인 수직 기준면으로서의 역할을 수행한다.

인간이 거주하는 터전인 지구의 형상을 파악하려는 노력은 인류 역사의 지적 발전과 궤를 같이해 왔다. 초기 고대 문명권에서는 직관적인 관찰에 의존하여 지구를 평면적인 원반이나 사각형으로 인식하는 지구 평면설이 지배적이었다. 그러나 고대 그리스 시대에 이르러 월식 때 달에 비치는 지구의 그림자가 곡선이라는 점, 북쪽이나 남쪽으로 이동함에 따라 보이는 별자리가 달라진다는 점 등을 근거로 지구가 구 형태라는 지구 구체설이 제기되었다. 특히 에라토스테네스(Eratosthenes)는 하짓날 시에네와 알렉산드리아에서의 태양 남중 고도 차이와 두 지점 사이의 거리를 이용하여 지구의 둘레를 최초로 계산함으로써, 지구 형상에 대한 논의를 신화의 영역에서 측지학(Geodesy)의 영역으로 전환하였다.

17세기 과학 혁명기에 접어들면서 아이작 뉴턴(Isaac Newton)은 만유인력의 법칙과 원심력의 원리를 결합하여 지구의 형상을 새롭게 정의하였다. 뉴턴은 지구가 자전함에 따라 발생하는 원심력으로 인해 적도 부근이 부풀어 오른 회전 타원체(Oblate Spheroid)의 형상을 띨 것이라고 이론적으로 예측하였다. 반면, 프랑스의 지오반니 도메니코 카시니(Giovanni Domenico Cassini)를 필두로 한 데카르트주의자들은 지구가 극 방향으로 길쭉한 장구형 타원체일 것이라고 주장하며 대립하였다. 이 논쟁은 프랑스 과학 아카데미가 파견한 라플란드와 페루 원정대의 위도 1도당 호의 길이 측정 결과, 고위도로 갈수록 호의 길이가 길어진다는 사실이 밝혀지면서 뉴턴의 편평 타원체 가설이 승리하는 것으로 종결되었다. 이후 지구의 기하학적 형상은 기준 타원체(Reference Ellipsoid)라는 수학적 모델로 정립되었다.

19세기에 이르러 정밀한 중력 관측과 수준 측량이 수행되면서, 지구의 실제 형상이 단순한 수학적 타원체와 일치하지 않는다는 사실이 드러나기 시작하였다. 카를 프리드리히 가우스(Carl Friedrich Gauss)는 지구의 물리적 표면을 “수학적으로 정의된 면이 아니라, 중력의 방향에 수직인 평형 상태의 해수면”으로 정의하며 이를 ’지구의 수학적 형상’이라 일컬었다. 가우스의 제자였던 요한 베네딕트 리스팅(Johann Benedict Listing)은 1873년 이러한 물리적 실체를 표현하기 위해 지오이드(Geoid)라는 용어를 처음으로 제안하였다. 리스팅은 지오이드를 평균 해수면을 육지 내부까지 연장하였을 때 상정되는 등전위면(Equipotential surface)으로 정의하였으며, 이는 지구 내부의 불균일한 질량 분포를 반영하는 고유한 물리적 형상임을 명시하였다.

현대 측지학에서 지구 형상의 표현은 기하학적인 기준 타원체와 물리적인 지오이드의 결합으로 완성된다. 20세기 중반 이후 인공위성을 이용한 중력 탐사와 위성 고도계 기술이 발전함에 따라, 전 지구적 규모에서의 지오이드 기복이 정밀하게 산출되었다. 오늘날 지오이드는 단순히 지구의 모양을 설명하는 단계를 넘어, GPS를 이용한 높이 측정의 기준이 되는 정표고 체계의 근간을 형성한다. 또한, 지오이드와 타원체 면 사이의 차이인 지오이드고와 중력의 방향 차이를 나타내는 수직선 편차에 대한 연구는 지구 내부의 밀도 구조와 맨틀 대류 등 지구물리학적 현상을 이해하는 핵심적인 도구로 활용되고 있다. 이러한 역사적 변천은 지구를 단순한 기하학적 대상에서 복잡한 물리적 역학 체계로 인식해 온 인류의 학술적 여정을 반영한다.

지오이드는 지구의 물리적 형상을 정의하는 가장 정밀한 모델로, 수학적으로는 지구 중력 포텐셜(gravity potential)이 일정한 값을 갖는 등전위면으로 정의된다. 지구 내부의 질량 분포가 불균일하고 자전으로 인한 원심력이 작용하기 때문에, 지오이드는 기하학적으로 단순한 구나 타원체와는 다른 복잡한 굴곡을 형성한다. 이를 수학적으로 기술하기 위해서는 지구 중력장의 물리적 특성을 지배하는 미분 방정식과 경계값 문제를 해결해야 한다.

지구의 총 중력 포텐셜 $ W $는 질량의 끌어당김에 의한 만유인력 포텐셜 $ V $와 지구 자전에 의한 원심력 포텐셜 $ $의 합으로 표현된다. $$ W(x, y, z) = V(x, y, z) + \Phi(x, y, z) $$ 여기서 원심력 포텐셜 $ $는 지구의 자전 각속도를 $ $, 회전축으로부터의 거리를 $ p $라 할 때 $ ^2 p^2 $으로 비교적 단순하게 정의되지만, 만유인력 포텐셜 $ V $는 지구 내부의 밀도 분포에 의존하므로 매우 복잡한 형태를 띤다. 질량이 존재하지 않는 지구 외부 공간에서 만유인력 포텐셜은 라플라스 방정식(Laplace’s equation)인 $ ^2 V = 0 $을 만족하며, 이를 구좌표계에서 풀이하면 구면 조화 함수(spherical harmonics)의 급수 전개 형태로 나타낼 수 있다⁷⁾.

지구 중력장 모델링에서 널리 사용되는 구면 조화 함수 전개식은 다음과 같다. $$ V(r, \theta, \lambda) = \frac{GM}{r} \left[ 1 + \sum_{n=2}^{\infty} \sum_{m=0}^{n} \left( \frac{a}{r} \right)^n \bar{P}_{nm}(\cos \theta) (\bar{C}_{nm} \cos m\lambda + \bar{S}_{nm} \sin m\lambda) \right] $$ 이 식에서 $ G $는 중력 상수, $ M $은 지구의 총 질량, $ a $는 지구의 적도 반지름을 의미한다. $ {P}%%//%%{nm} $은 정규화된 연관 르장드르 함수이며, $ {C}%%//%%{nm} $과 $ {S}_{nm} $은 지구 내부의 질량 분포 특성을 담고 있는 구면 조화 계수(Stokes coefficients)이다. 이 계수들은 인공위성 중력 탐사를 통해 정밀하게 산출되며, 전 지구적 지오이드 모델을 구축하는 핵심 데이터가 된다.

실제 지오이드의 형상을 결정할 때는 기하학적으로 정의된 기준 타원체(reference ellipsoid)와의 차이를 분석하는 방식을 취한다. 기준 타원체상의 정상 포텐셜(normal potential)을 $ U $라 할 때, 실제 포텐셜 $ W $와의 차이를 교란 포텐셜(disturbing potential) $ T = W - U $로 정의한다. 지오이드와 타원체 사이의 수직 거리인 지오이드고(geoid height) $ N $은 브룬스 공식(Bruns’ formula)에 의해 다음과 같이 교란 포텐셜과 관계를 맺는다. $$ N = \frac{T}{\gamma} $$ 여기서 $ $는 기준 타원체에서의 정상 중력을 의미한다. 즉, 특정 지점의 교란 포텐셜을 알 수 있다면 해당 지점의 지오이드 높이를 수학적으로 산출할 수 있다.

지표면에서 관측된 중력 데이터를 이용하여 지오이드를 결정할 때는 스토크스 공식(Stokes’ formula)이 사용된다. 이는 지표면 전체의 중력 이상(gravity anomaly) 데이터를 적분하여 특정 지점의 지오이드고를 구하는 방법이다. $$ N = \frac{R}{4\pi\gamma} \iint_{\sigma} \Delta g S(\psi) d\sigma $$ 이 식에서 $ R $은 지구의 평균 반지름, $ g $는 중력 이상, $ S() $는 관측점과 적분 요소 사이의 각거리에 따른 가중치 함수인 스토크스 함수이다. 이 공식은 지구 전체에 걸친 중력 관측값이 존재한다는 가정하에 성립하며, 현대 측지학에서는 이를 보완하기 위해 위성 데이터와 지상 관측 데이터를 결합한 하이브리드 모델을 사용하여 지오이드의 정밀도를 높이고 있다⁸⁾.

지오이드를 물리적으로 정의하고 그 형상을 수학적으로 기술하기 위해서는 지구 중력장의 에너지 상태를 나타내는 중력 포텐셜(Gravity Potential) 이론에 대한 이해가 선행되어야 한다. 지구상의 임의의 점에서의 중력 포텐셜 $W$는 지구의 질량 분포에 의한 만유인력(Gravitation) 포텐셜 $V$와 지구의 자전(Rotation)으로 인해 발생하는 원심력(Centrifugal Force) 포텐셜 $\Phi$의 산술적인 합으로 정의된다. 즉, 중력 포텐셜은 다음과 같은 기본 관계식을 만족한다.

$$W = V + \Phi$$

만유인력 포텐셜 $V$는 뉴턴의 만유인력 법칙에 따라 지구 전체의 질량 분포를 적분하여 산출된다. 지구의 총 질량을 $M$, 만유인력 상수를 $G$, 지구 내부의 미소 질량을 $dm$, 그리고 관측점으로부터 미소 질량까지의 거리를 $r$이라 할 때, 만유인력 포텐셜은 다음과 같이 표현된다.

$$V = G \int_{M} \frac{1}{r} dm$$

이 포텐셜 함수 $V$는 지구 질량 외부의 공간에서 라플라스 방정식(Laplace’s equation)을 만족하며, 지구 내부에서는 질량 밀도 분포에 따른 푸아송 방정식(Poisson’s equation)을 따른다. 한편, 원심력 포텐셜 $\Phi$는 지구가 일정한 각속도(Angular Velocity) $\omega$로 자전한다고 가정할 때, 회전축으로부터의 수직 거리 $p$에 의존하는 함수로 정의된다.

$$\Phi = \frac{1}{2} \omega^2 p^2 = \frac{1}{2} \omega^2 (x^2 + y^2)$$

여기서 $x, y$는 지구 중심을 원점으로 하고 자전축을 $z$축으로 하는 직교 좌표계에서의 수평 좌표를 의미한다. 중력 포텐셜 $W$의 공간적 변화율, 즉 구배(Gradient)는 해당 지점에서의 중력 가속도(Gravity Acceleration) 벡터 $\vec{g}$를 형성한다.

$$\vec{g} = \nabla W = \left( \frac{\partial W}{\partial x}, \frac{\partial W}{\partial y}, \frac{\partial W}{\partial z} \right)$$

중력 가속도 벡터는 항상 중력 포텐셜이 일정한 면인 등전위면(Equipotential Surface)에 수직으로 작용한다. 측지학적 관점에서 지오이드는 이러한 무수히 많은 등전위면 중 평균 해수면(Mean Sea Level)에 가장 잘 부합하는 특정한 등전위면 $W = W_0$로 정의된다. 따라서 지오이드면 위의 모든 지점은 동일한 중력 포텐셜 값을 가지며, 해당 면에서의 중력 방향은 항상 지오이드의 법선(Normal) 방향과 일치하게 된다.

지구 내부의 질량 분포가 불균일하기 때문에 만유인력 포텐셜 $V$는 복잡한 형태를 띠게 되며, 이는 결과적으로 지오이드가 기하학적 타원체(Ellipsoid)와 달리 굴곡진 형태를 갖게 만드는 주된 원인이 된다. 중력 포텐셜 이론은 이러한 물리적 불균질성을 수학적으로 정량화하여, 지구의 실제 형상을 모델링하는 데 핵심적인 이론적 토대를 제공한다⁹⁾.

지구의 질량 분포에 의해 발생하는 만유인력 포텐셜(Gravitational Potential) $V$는 질량이 존재하지 않는 지구 외부 공간에서 특별한 수학적 성질을 갖는다. 뉴턴의 만유인력 법칙에 기초한 포아송 방정식(Poisson’s equation)에 따르면, 임의의 지점에서의 포텐셜과 밀도 $\rho$ 사이에는 $\nabla^2 V = -4\pi G\rho$라는 관계가 성립한다. 여기서 $G$는 만유인력 상수이다. 그러나 지구 표면 밖의 대기 밀도를 무시할 수 있는 진공 영역에서는 밀도 $\rho$가 0이 되므로, 위 식은 다음과 같은 라플라스 방정식(Laplace’s equation)으로 귀결된다.

$$\nabla^2 V = \frac{\partial^2 V}{\partial x^2} + \frac{\partial^2 V}{\partial y^2} + \frac{\partial^2 V}{\partial z^2} = 0$$

이 방정식을 만족하는 함수 $V$를 조화 함수(Harmonic function)라고 정의한다. 조화 함수는 영역 내부에서 극대값이나 극소값을 갖지 않으며, 경계면에서의 값에 의해 내부의 분포가 유일하게 결정되는 특성을 지닌다. 이러한 수학적 엄밀함은 지구 외부 중력장의 구조를 파악하고, 지표면 관측 데이터를 바탕으로 공간상의 포텐셜을 추론하는 이론적 토대가 된다.

지오이드의 형상을 결정하는 과정은 수학적으로 측지학적 경계값 문제(Geodetic Boundary Value Problem, GBVP)로 정형화된다. 이는 지구 표면 혹은 지표 인근에서 관측된 물리적 데이터인 중력 가속도나 포텐셜 차이를 경계 조건(Boundary condition)으로 설정하고, 라플라스 방정식을 만족하는 포텐셜 함수를 전 지구 영역에서 산출하는 문제이다. 측지학에서 다루는 경계값 문제는 크게 세 가지 유형으로 분류된다. 첫째, 경계면에서 포텐셜 값 자체가 주어지는 디리클레 문제(Dirichlet problem)이다. 둘째, 경계면에서 포텐셜의 법선 도함수, 즉 중력의 수직 성분이 주어지는 노이만 문제(Neumann problem)이다. 셋째, 포텐셜과 그 도함수의 선형 결합이 경계 조건으로 주어지는 로빈 문제(Robin problem) 혹은 혼합 경계값 문제이다.

실제 지오이드 결정에서는 지표면에서 측정된 중력 이상(Gravity anomaly)을 경계 조건으로 활용하게 되는데, 이는 포텐셜의 미분량과 포텐셜 자체의 항이 복합적으로 포함된 형태이므로 주로 로빈 문제의 범주에 속한다. 특히 조지 가브리엘 스토크스가 제안한 고전적 해법은 지구를 구로 가정하고 지표면의 중력 데이터를 적분하여 지오이드고를 구하는 방식을 취하며, 이는 현대적인 물리 측지학의 시초가 되었다.

다만 실제 지구는 자전 운동을 하므로, 순수한 만유인력 포텐셜 $V$ 외에도 자전에 의한 원심력 포텐셜 $\Phi$를 합산한 전체 중력 포텐셜 $W = V + \Phi$를 고려해야 한다. 원심력 포텐셜은 라플라스 방정식을 만족하지 않고 $\nabla^2 \Phi = 2\omega^2$ (여기서 $\omega$는 지구 자전 각속도)를 따르기 때문에, 전체 중력 포텐셜 $W$는 엄밀하게는 조화 함수가 아니다. 따라서 정밀한 지오이드 모델링 과정에서는 원심력 효과를 수학적으로 분리하거나 보정하여, 라플라스 방정식의 해법인 구면 조화 함수(Spherical harmonics) 전개를 적용할 수 있는 형태로 문제를 변환하여 해결한다. 이러한 해석적 접근을 통해 도출된 해는 지구 중력 모델(Earth Gravitational Model)의 계수로 표현되어 전 지구적인 지오이드 기복을 기술하는 데 사용된다.¹⁰⁾

지구 외부 공간에서 중력 포텐셜(Gravity Potential)은 라플라스 방정식(Laplace’s equation)을 만족하며, 이를 전 지구적 규모에서 수학적으로 엄밀하게 표현하기 위해 구면 조화 함수(Spherical Harmonics)를 이용한 무한 급수 전개 방식이 사용된다. 이 방법은 지구의 불규칙한 질량 분포에 의한 중력장의 변동을 구면 좌표계상의 직교 함수들의 합으로 분해하여 기술하는 체계이다. 구면 조화 함수 전개는 지오이드의 전 지구적 형상을 결정하고, 지구 중력장 모델(Earth Gravitational Model)을 구축하는 데 있어 핵심적인 역할을 수행한다.

지구 외부의 임의의 점 $P(r, \theta, \lambda)$에서의 만유인력 포텐셜 $V$는 다음과 같은 구면 조화 급수로 표현된다.

$$V(r, \theta, \lambda) = \frac{GM}{r} \left[ 1 + \sum_{n=2}^{\infty} \sum_{m=0}^{n} \left( \frac{a}{r} \right)^n \left( C_{nm} \cos m\lambda + S_{nm} \sin m\lambda \right) P_{nm}(\cos \theta) \right]$$

이 식에서 $G$는 중력 상수, $M$은 지구의 총 질량이며, $r$은 지구 중심으로부터의 거리, $\theta$는 여위도(Colatitude), $\lambda$는 경도를 나타낸다. $a$는 지구의 장반경으로서 기준 거리가 된다. $P_{nm}(\cos \theta)$는 $n$차(degree) $m$계(order)의 르장드르 연관 함수(Associated Legendre function)이며, $C_{nm}$과 $S_{nm}$은 지구 내부의 질량 분포 특성을 담고 있는 구면 조화 계수이다. 이 계수들은 중력 관측 데이터를 통해 결정되며, 흔히 스토크스 계수(Stokes’ coefficients)라고도 불린다.

급수의 각 항은 지구 중력장의 서로 다른 물리적 특성을 반영한다. $n=0$인 항은 지구가 점질량이라고 가정했을 때의 기본 포텐셜을 의미한다. $n=1$인 항은 좌표계의 원점이 지구의 질량 중심(Center of mass)과 일치할 경우 이론적으로 0이 되어 소거된다. $n=2$인 항은 지구의 자전에 의한 편평도와 관련이 깊으며, 특히 $C_{20}$ 계수는 지구의 관성 모멘트 차이를 나타내는 핵심 인자로 동역학적 형상 계수(Dynamic form factor)라고도 한다. $n$의 값이 커질수록 더 짧은 파장의 중력 변동을 묘사하게 되며, 이는 지각 내부의 밀도 불균형이나 국지적인 지형 기복에 의한 중력장의 세밀한 구조를 반영한다.

실제 수치 계산에서는 계수의 값이 차수가 높아짐에 따라 급격히 작아지는 문제를 해결하기 위해 완전 정규화된 구면 조화 함수(Fully normalized spherical harmonics)를 사용한다. 정규화된 계수 $\bar{C}_{nm}, \bar{S}_{nm}$와 정규화된 함수 $\bar{P}_{nm}$을 사용하면 수치적 안정성을 확보할 수 있으며, 전 지구 중력장 모델인 EGM96이나 EGM2008 등은 이러한 정규화된 계수들을 수천 차수까지 제공함으로써 지오이드의 정밀한 형상을 복원한다.

구면 조화 함수 전개를 통해 얻어진 포텐셜 함수는 기준 타원체의 포텐셜과 비교되어 지오이드고(Geoid height)를 산출하는 기초가 된다. 전개의 최대 차수 $n_{max}$는 모델의 해상도를 결정하며, 대략 $180^\circ / n_{max}$의 공간 해상도를 갖는다. 따라서 전 지구적 규모의 완만한 지오이드 굴곡은 낮은 차수의 항들에 의해 결정되고, 지각의 밀도 차이나 산맥 등과 같은 국지적 요인에 의한 미세한 굴곡은 높은 차수의 항들에 의해 표현된다. 이러한 직교 함수 전개 방식은 복잡한 지구 중력장을 체계적으로 분류하고 해석할 수 있는 강력한 수학적 도구를 제공한다.

지오이드의 형상을 결정하는 핵심적인 수치적 방법론 중 하나는 중력 이상(Gravity anomaly) 데이터를 적분하여 특정 지점의 지오이드고(Geoid height)를 산출하는 것이다. 이러한 물리적 측지학의 고전적 해법을 제시한 인물은 조지 가브리엘 스토크스(George Gabriel Stokes)로, 그는 1849년에 지구 외부의 질량이 존재하지 않는다는 가정하에 라플라스 방정식의 경계값 문제를 해결하여 스토크스 공식(Stokes’ formula)을 도출하였다. 이 공식은 지표면에서 관측된 중력값과 준거 타원체상의 이론적 중력값 사이의 차이인 중력 이상을 전 지구적으로 적분함으로써, 타원체로부터 지오이드까지의 수직 거리인 지오이드고를 직접 계산할 수 있게 한다.

스토크스 공식의 이론적 토대는 브룬스 공식(Bruns’ formula)에 기반한다. 브룬스 공식은 교란 포텐셜(Disturbing potential) $ T $와 지오이드고 $ N $ 사이의 관계를 다음과 같이 정의한다. $$ N = \frac{T}{\gamma} $$ 여기서 $ $는 정상 중력(Normal gravity)을 의미한다. 스토크스는 지구 내부의 질량을 고려하지 않는 경계 조건하에서 포텐셜 이론을 적용하여, 지표면 전체에 걸친 중력 이상의 분포로부터 교란 포텐셜을 구하는 적분 방정식을 유도하였다. 구형 근사(Spherical approximation)를 적용한 스토크스 공식의 일반적인 형태는 다음과 같다. $$ N = \frac{R}{4\pi G} \iint_{\sigma} \Delta g S(\psi) d\sigma $$ 이 식에서 $ R $은 지구의 평균 반지름, $ G $는 평균 중력, $ g $는 중력 이상, $ $는 단위 구면을 의미한다. 이때 $ S() $는 스토크스 함수(Stokes’ function) 또는 스토크스 핵(Stokes’ kernel)이라 불리며, 계산 지점과 중력 이상 관측 지점 사이의 각거리(Angular distance) $ $에 따라 가중치를 결정하는 역할을 수행한다.

스토크스 함수는 두 지점 사이의 기하학적 관계에 따라 다음과 같이 복잡한 삼각함수의 조합으로 표현된다. $$ S(\psi) = \csc\left(\frac{\psi}{2}\right) - 6\sin\left(\frac{\psi}{2}\right) + 1 - 5\cos\psi - 3\cos\psi \ln\left(\sin\frac{\psi}{2} + \sin^2\frac{\psi}{2}\right) $$ 이 함수는 각거리 $ $가 0에 가까워질수록 급격히 증가하는 특성을 가진다. 이는 계산 지점 인근의 중력 데이터가 지오이드고 결정에 가장 지대한 영향을 미친다는 것을 의미하며, 원거리로 갈수록 그 영향력이 점차 감소하다가 특정 거리에서 음의 값을 갖는 등 주기적인 변동성을 보인다. 따라서 정밀한 지오이드 산출을 위해서는 계산 지점 주변의 고해상도 중력 자료가 필수적이다.

실제 적용 과정에서 스토크스 적분은 몇 가지 전제 조건을 요구한다. 우선 지구 외부, 즉 지표면과 지오이드 사이에 질량이 존재하지 않아야 한다는 가정이 필요하므로, 지표면에서 관측된 중력값에 대기 보정 및 지형 보정을 가하여 모든 질량을 지오이드 내부로 응축시키는 중력 약약 보정(Helmert’s condensation reduction) 등의 전처리 과정이 선행되어야 한다. 또한, 이론적으로는 전 지구적 규모의 중력 이상 데이터가 필요하지만, 실제로는 데이터의 공백이나 불균일한 분포로 인해 오차가 발생할 수 있다.

현대 측지학에서는 이러한 한계를 극복하기 위해 구면 조화 함수(Spherical harmonics)를 이용한 전 지구 중력장 모델(Global Gravitational Model, GGM)과 국지적인 중력 이상 데이터를 결합하는 방식을 주로 사용한다. 장파장(Long-wavelength) 성분의 지오이드 굴곡은 위성 궤도 분석 등을 통해 얻어진 중력장 모델로 계산하고, 단파장(Short-wavelength) 성분의 미세한 변화만을 스토크스 공식을 통해 보완함으로써 계산의 효율성과 정확도를 동시에 확보한다. 이러한 결합 모델은 위성 항법 시스템(Global Navigation Satellite System, GNSS)을 이용한 고도 결정에서 지오이드고를 보정하여 정표고를 산출하는 데 핵심적인 역할을 수행한다.

현대 측지학에서 지오이드의 형상을 정밀하게 결정하는 과정은 지구 내부의 질량 분포와 외부 중력장의 특성을 물리적·기하학적으로 해석하는 복합적인 작업이다. 지오이드 결정의 핵심은 지표면이나 우주 공간에서 관측된 중력 데이터를 바탕으로, 기준 타원체로부터 지오이드까지의 수직 거리인 지오이드고(Geoid height)를 산출하는 데 있다. 이를 위해 전통적인 중력 측량 방식부터 최첨단 인공위성 탐사 기법에 이르기까지 다양한 기술이 상호 보완적으로 활용된다.

물리적 방식의 핵심인 중력 지오이드 결정은 지표에서 관측된 중력 데이터를 바탕으로 스토크스 공식(Stokes’ formula)을 적용하여 계산한다. 이 방법은 지표면 전체의 중력 이상(Gravity anomaly) 데이터를 적분하여 특정 지점의 지오이드고를 구하는 원리이다. 그러나 현실적으로 전 지구의 중력 데이터를 조밀하게 획득하기 어렵기 때문에, 대규모 질량 분포를 반영하는 글로벌 중력 모델(Global Gravity Model, GGM)과 지역적인 세부 중력 데이터를 결합하는 제거-계산-복원(Remove-Compute-Restore, RCR) 기법이 주로 사용된다. 이 과정에서 지형의 영향을 보정하기 위한 지형 보정(Terrain correction)과 밀도 불균형에 따른 응축 보정(Condensation reduction) 등의 정교한 데이터 처리가 수반된다.

해양 영역에서는 위성 고도계(Satellite Altimeter) 기술이 지오이드 결정에 결정적인 역할을 한다. 인공위성에서 해수면을 향해 마이크로파를 발사하고 반사되어 돌아오는 시간을 측정함으로써 해수면의 높이를 파악할 수 있다. 해류나 기압의 영향을 제외한 정지된 상태의 평균 해수면은 이론적으로 지오이드와 일치하므로, 위성 고도계 데이터는 해양 지오이드의 형상을 직접적으로 투영한다. GEOSAT, ERS, Jason 시리즈와 같은 위성들은 수십 년간 축적된 데이터를 통해 전 지구적인 해양 지오이드 모델의 정밀도를 비약적으로 향상시켰다.

최근에는 위성 중력 탐사(Satellite Gravimetry) 전용 미션을 통해 지구 중력장의 장파장 성분을 극도로 정밀하게 관측하고 있다. CHAMP, GRACE(Gravity Recovery and Climate Experiment), GOCE(Gravity field and steady-state Ocean Circulation Explorer) 위성 등은 위성 간 거리 측정(Satellite-to-Satellite Tracking, SST)이나 위성 중력 변화계(Satellite Gravity Gradiometry, SGG)를 이용하여 지구 중력장의 미세한 변화를 감지한다. 특히 GOCE 미션은 지구 중력장의 고주파 성분까지 포착하여 수 센티미터 오차 범위 내의 정밀한 지오이드 모델을 구축하는 데 기여하였다¹¹⁾. 이러한 위성 데이터는 구면 조화 함수(Spherical harmonic function) 형태로 전개되어 전 지구적 기준면을 정의하는 기초 자료가 된다.

기하학적 지오이드 결정법은 위성 항법 시스템(Global Navigation Satellite System, GNSS)과 전통적인 수준 측량(Leveling) 결과를 비교하는 방식이다. GNSS를 통해 얻은 타원체고(Ellipsoidal height, $h$)와 수준 측량을 통해 얻은 표고(Orthometric height, $H$)의 차이를 이용하면 해당 지점의 지오이드고($N$)를 $N = h - H$의 관계식을 통해 직접 구할 수 있다. 이 방법은 국지적인 지역에서 매우 높은 정밀도를 보장하며, 물리적 중력 지오이드 모델의 왜곡을 보정하고 지역 좌표계와 결합하여 하이브리드 지오이드(Hybrid geoid) 모델을 제작하는 데 필수적이다¹²⁾.

결과적으로 현대의 지오이드 결정은 전 지구적 중력 모델, 고해상도 지형 모델, 그리고 지상의 GNSS/수준 측량 데이터를 통합하는 방향으로 발전하고 있다. 이러한 통합 관측 기법은 국가 고도 체계의 정밀도를 높일 뿐만 아니라, 해수면 상승 모니터링이나 지구 내부의 밀도 구조 연구 등 다양한 학술 및 실무 분야에 핵심적인 기초 정보를 제공한다.

지상 중력 측량은 지표면의 특정 지점에서 지구가 작용하는 실제 중력 가속도를 정밀하게 측정하는 과정으로, 지오이드 결정의 기초 자료를 제공한다. 이는 중력계(gravimeter)를 통해 수행되며, 측정 방식에 따라 절대 중력 측정과 상대 중력 측정으로 구분된다. 절대 중력 측정은 진공 상태에서 자유 낙하하는 물체의 가속도를 레이저 간섭계 등으로 직접 측정하여 해당 지점의 절대적인 중력값을 얻는 방식이다. 반면 상대 중력 측정은 스프링의 변위 등을 이용하여 이미 중력값을 알고 있는 기지점과 미지점 사이의 중력 차이를 측정한다. 지표면에서 관측된 중력값은 지구 내부의 밀도 불균형, 지형적 기복, 그리고 관측점의 고도 등 복합적인 물리적 정보를 내포하고 있다.

지오이드를 물리적으로 정의하기 위해서는 실제 관측된 중력값과 비교할 수 있는 이론적 기준 모델이 필요하다. 이를 위해 지구를 질량이 균일하고 기하학적으로 매끄러운 회전 타원체로 가정하는데, 이 가상 모델 위에서 계산된 중력값을 정상 중력(normal gravity)이라 한다. 정상 중력은 위도에 따른 원심력의 변화와 타원체 형상을 반영한 수식에 의해 결정되며, 현대 측지학에서는 GRS80(Geodetic Reference System 1980)이나 WGS84(World Geodetic System 1984)와 같은 표준 타원체를 기준으로 산출한다. 정상 중력은 고도에 따라 감소하는 특성을 가지므로, 지표면에서 측정된 값과 직접 비교하기 위해서는 두 값을 동일한 고도상에 위치시키는 과정이 선행되어야 한다.

지표면에서 측정된 중력값을 기준면인 지오이드나 타원체상에서의 값으로 변환하는 과정을 중력 보정(gravity reduction)이라 한다. 가장 대표적인 보정은 프리 에어 보정(free-air correction)으로, 관측점과 기준면 사이에 질량이 존재하지 않는다고 가정하고 고도 차이에 의한 중력 감소분만을 보정하는 방식이다. 지표면에서의 관측 중력을 $g$, 해수면(지오이드)으로 보정된 중력을 $g_{0}$, 관측점의 고도를 $h$라 할 때, 프리 에어 보정량은 대략 $0.3086 \times h \text{ mGal}$의 비율로 더해진다. 또한 관측점 하부에 존재하는 지각의 질량 효과를 제거하는 부게 보정(Bouguer correction)과 주변 지형의 기복에 의한 인력을 고려하는 지형 보정(terrain correction) 등을 통해 중력 데이터를 정밀하게 가공한다.

중력 이상(gravity anomaly)은 이와 같이 보정된 실제 중력값과 이론적인 정상 중력값의 차이로 정의된다. 지오이드 결정에서 가장 중요한 역할을 하는 것은 프리 에어 이상(free-air anomaly)이다. 이는 프리 에어 보정만을 거친 중력값에서 정상 중력값을 뺀 것으로, 지구 내부의 질량 분포와 지오이드의 기복 정보를 직접적으로 반영한다. 프리 에어 이상 $\Delta g$는 물리적으로 브룬스 정리(Bruns’ theorem)에 의해 교란 포텐셜(disturbing potential)과 연결되며, 이는 지오이드가 기준 타원체로부터 얼마나 떨어져 있는지를 나타내는 지오이드고(geoid height)를 산출하는 핵심 변수가 된다.

최종적으로 지상 중력 측량을 통해 확보된 전 지구적 혹은 국지적 중력 이상 데이터는 스토크스 공식(Stokes’ formula)에 입력되어 지오이드의 형상을 수치적으로 결정하는 데 사용된다. 스토크스 공식은 지표면 전체의 중력 이상을 적분하여 특정 지점의 지오이드고를 구하는 수학적 틀을 제공한다. 지상 중력 데이터는 위성 중력 탐사가 포착하기 어려운 미세한 지각 밀도 변화와 급격한 지형 변화를 정밀하게 반영할 수 있다는 장점이 있다. 따라서 정밀한 국가 지오이드 모델을 구축하기 위해서는 조밀한 격자 형태로 수행된 지상 중력 측량 성과와 위성 데이터를 결합하는 과정이 필수적이다. 이러한 과정을 통해 도출된 지오이드는 수준 측량의 기준면이 되어 정밀한 고도 체계를 확립하는 토대가 된다.

현대 측지학에서 지오이드의 정밀한 결정은 지상 중력 관측의 한계를 극복하기 위해 인공위성을 활용한 원격 탐사 기술에 크게 의존하고 있다. 지상 관측은 접근이 어려운 오지나 광대한 해양 영역에서 데이터의 공백이 발생하며, 전 지구적인 통일된 기준계를 설정하는 데 어려움이 따른다. 이러한 한계를 보완하기 위해 도입된 위성 고도계(Satellite Altimetry) 기술은 위성에서 발사한 마이크로파가 해수면에 반사되어 돌아오는 왕복 시간을 측정함으로써 위성과 해수면 사이의 거리를 산출한다. 위성의 궤도 고도를 $ H $, 측정된 거리를 $ $라고 할 때, 해수면의 타원체 고도 $ h_{sea} $는 다음과 같은 기본적인 관계식으로 표현된다.

$$ h_{sea} = H - \rho - \Delta \rho $$

여기서 $ $는 대기 지연, 조석 효과, 해수면 기압 변화 등에 따른 보정 항을 의미한다. 관측된 해수면은 평균적으로 지오이드와 일치하는 경향이 있으나, 해류와 해수의 밀도 차이로 인해 발생하는 해양 역학적 지형(Mean Dynamic Topography, MDT)만큼의 차이가 존재한다. 따라서 위성 고도계 데이터를 통해 지오이드고를 정확히 도출하기 위해서는 해양 순환 모델을 이용한 MDT의 정밀한 분리가 필수적이다.

2000년대 이후에는 지구 중력장의 미세한 변화를 직접 측정하기 위한 전용 중력 탐사 위성들이 발사되면서 지오이드 모델링의 패러다임이 전환되었다. 대표적인 미션인 GRACE(Gravity Recovery and Climate Experiment)는 동일한 궤도상에서 약 220km 떨어진 두 개의 위성이 서로를 추적하는 위성 간 추적(Satellite-to-Satellite Tracking, SST) 방식을 사용한다. 지구 내부의 질량 분포 차이에 의해 앞선 위성과 뒤처진 위성 사이의 거리가 미세하게 변하는 것을 마이크로파 또는 레이저 간섭계로 측정함으로써, 지구 중력장의 시간적·공간적 변화를 고해상도로 파악한다. 특히 GRACE는 빙하의 융해나 지하수량 변화와 같은 중력장의 시간적 변동성을 관측하는 데 탁월한 성능을 발휘하였다¹³⁾.

더욱 정밀한 정적 지오이드 모델을 구축하기 위해 유럽우주국(ESA)은 GOCE(Gravity field and steady-state Ocean Circulation Explorer) 위성을 운용하였다. GOCE는 중력 그라디오메트리(Satellite Gravity Gradiometry, SGG) 기술을 도입하여, 위성 내부에 탑재된 3축 가속도계 쌍을 통해 중력의 공간 변화율인 중력 경사(Gravity gradient)를 직접 측정하였다. 이는 기존의 궤도 섭동 분석 방식보다 짧은 파장의 중력장 신호를 포착하는 데 유리하며, 이를 통해 공간 해상도 약 100km 범위에서 1~2cm 오차 범위 내의 정밀도를 가진 지오이드를 결정하는 것이 가능해졌다¹⁴⁾.

이러한 위성 탐사 데이터는 구면 조화 함수(Spherical Harmonic Function) 전개 방식을 통해 전 지구 지구 중력장 모델(Earth Gravitational Model, EGM)로 통합된다. 대표적인 모델인 EGM2008은 위성 중력 데이터와 지상 중력 관측치, 위성 고도계 데이터를 결합하여 구축되었으며, 최대 2,190차까지 전개되어 전 지구적인 지오이드 기복을 매우 상세하게 표현한다¹⁵⁾. 위성 기술을 통한 중력장 관측은 단순히 지오이드의 형상을 정의하는 것에 그치지 않고, 해수면 상승, 판 구조론에 따른 질량 이동, 그리고 전 지구적 고도 체계의 통합을 위한 물리적 근거를 제공한다는 점에서 학문적·실무적 가치가 매우 높다.

기하학적 지오이드 결정법은 지구의 물리적 특성인 중력을 직접 측정하지 않고, 서로 다른 두 기준면으로부터 측정한 고도 값의 기하학적 차이를 이용하여 지오이드고(Geoid height)를 산출하는 방식이다. 이 방법은 현대 측지학에서 위성 항법 시스템(Global Navigation Satellite System, GNSS)의 보급과 함께 그 중요성이 비약적으로 증대되었다. 기하학적 결정법의 핵심은 우주 측지 기술을 통해 얻은 수리적 위치 정보와 지표면에서의 전통적인 측량 결과를 결합하여, 이론적 모델인 기준 타원체(Reference ellipsoid)와 물리적 실체인 지오이드 사이의 간격인 지오이드고를 직접 계산해내는 데 있다.

이 방법의 논리적 기초는 고도 체계 간의 상관관계에 있다. GNSS를 활용하면 지구 질량 중심을 원점으로 하는 지구 중심 지구 고정(Earth-Centered, Earth-Fixed, ECEF) 좌표계 상에서 정밀한 위치를 파악할 수 있으며, 이를 통해 특정 지점의 타원체 고도(Ellipsoidal height)를 얻게 된다. 타원체 고도는 지표면의 한 점으로부터 기준 타원체의 법선 방향으로 잰 수직 거리를 의미하며, 순수하게 기하학적으로 정의된 값이다. 반면, 우리가 일상적으로 사용하는 표고(Orthometric height)는 지오이드로부터 지표면까지의 높이를 의미하며, 이는 수준 측량(Leveling)과 중력 보정을 통해 결정되는 물리적 고도이다.

타원체 고도, 표고, 그리고 지오이드고 사이의 기하학적 관계는 다음과 같은 기본 방정식으로 표현된다.

$ h = H + N $

여기서 $ h $는 타원체 고도, $ H $는 표고(주로 정표고), $ N $은 지오이드고를 나타낸다. 따라서 특정 지점에서 GNSS 관측을 통해 $ h $를 구하고, 직접 수준 측량을 실시하여 $ H $를 확보한다면, 해당 지점에서의 지오이드고는 $ N = h - H $라는 단순한 차이 계산을 통해 도출된다. 이러한 방식을 GNSS 수준 측량이라 하며, 물리적 중력 데이터가 부족한 지역에서도 지오이드의 형상을 파악할 수 있는 유효한 수단을 제공한다.

기하학적 지오이드 결정법을 광역적인 모델링으로 확장하기 위해서는 다수의 공통점(Common points) 확보가 필수적이다. 타원체 고도와 표고를 모두 알고 있는 점들을 격자 형태로 배치하고, 각 지점에서 산출된 지오이드고 데이터를 바탕으로 보간법(Interpolation)을 적용한다. 이때 주로 사용되는 수학적 기법으로는 최소제곱법(Least squares method)에 기반한 다항식 적합이나 크리깅(Kriging) 등이 있으며, 이를 통해 미관측 지점의 지오이드고를 추정하는 국지적 지오이드 모델을 구축한다. 이렇게 구축된 모델은 GNSS 관측만으로도 즉각적인 표고 변환을 가능하게 하여 측량 작업의 효율성을 극대화한다.

그러나 기하학적 결정법은 몇 가지 명확한 한계를 지닌다. 우선, 산출된 지오이드고의 정밀도가 GNSS 관측 오차와 수준 측량 오차에 직접적으로 종속된다는 점이다. 특히 수준 측량은 장거리를 이동하며 수행해야 하므로 오차 전파(Error propagation)의 영향이 크며, 이는 지오이드 모델의 왜곡을 야기할 수 있다. 또한, 이 방법은 관측점이 존재하는 육상 지역으로 활용 범위가 제한되며, 지형의 기복이 심하거나 중력 이상이 급격히 변하는 지역에서는 단순한 기하학적 보간만으로 지오이드의 미세한 굴곡을 완벽히 재현하기 어렵다. 따라서 현대 측지학에서는 기하학적 방법으로 얻은 국지적 데이터와 중력 관측을 통한 물리적 데이터를 결합하여 모델의 정밀도를 보완하는 방식을 주로 채택한다.

그럼에도 불구하고 기하학적 지오이드 결정법은 국가 기준점 체계의 현대화와 고도화에 있어 핵심적인 역할을 수행한다. 특히 지각 변동이나 해수면 상승으로 인해 기존의 수직 기준면이 변하는 경우, GNSS를 이용한 주기적인 기하학적 관측은 지오이드의 시계열적 변화를 모니터링하고 정밀한 지도 제작 및 토목 공사의 기준을 유지하는 데 결정적인 기여를 한다. 이는 위성 기술과 전통적 측량학이 결합하여 지구의 형상을 더욱 입체적이고 정밀하게 이해하게 된 대표적인 사례라 할 수 있다.

지구의 물리적 형상을 정의하는 데 있어 수학적 편의를 위해 도입된 기준 타원체(Reference Ellipsoid)와 지구 중력장의 물리적 실체를 반영하는 지오이드(Geoid) 사이의 기하학적 관계를 이해하는 것은 측지학 및 공간 정보 공학의 기초가 된다. 기준 타원체는 지구의 전체적인 모양을 회전 타원체로 근사한 기하학적 모델이며, 지오이드는 평균 해수면을 연장하여 지구 전체를 감싸는 중력 등전위면으로 정의되는 물리적 모델이다. 이 두 면은 지구 내부의 질량 분포가 불균일하기 때문에 일치하지 않으며, 그 차이를 정량화함으로써 지표면 위 임의의 점에 대한 정확한 위치와 높이를 결정할 수 있다.

지표면의 한 점에 대한 높이를 표현하는 방식은 기준면의 설정에 따라 크게 두 가지로 나뉜다. 위성 항법 시스템(Global Navigation Satellite System, GNSS)을 통해 직접적으로 얻어지는 높이는 기준 타원체의 법선 방향으로 측정된 타원체고(Ellipsoidal Height, $ h $)이다. 반면, 실제 지형에서 물의 흐름이나 중력의 영향을 반영하는 물리적 높이는 지오이드로부터 연직선을 따라 측정된 정표고(Orthometric Height, $ H $)이다. 이 두 높이 사이의 기하학적 차이를 지오이드고(Geoid Height) 또는 지오이드 기복(Geoid Undulation, $ N $)이라 하며, 다음과 같은 기본적인 관계식이 성립한다.

$$ h = H + N $$

위 식에서 알 수 있듯이, 지오이드가 기준 타원체보다 위쪽에 위치하면 지오이드고 $ N $은 양(+)의 값을 가지며, 반대로 타원체 아래에 위치하면 음(-)의 값을 가진다. 이러한 편차는 지구 내부의 밀도 불균형과 지형의 기복에 의해 발생한다. 예를 들어, 거대한 산맥이나 밀도가 높은 지각 물질이 존재하는 지역에서는 중력이 강해져 지오이드가 타원체 밖으로 팽창하게 된다. 따라서 정밀한 지형 측량이나 지도 제작을 위해서는 GNSS로 측정된 타원체고를 실무에서 사용하는 정표고로 변환해야 하며, 이를 위해 해당 지역의 정밀한 지오이드 모델이 필수적으로 요구된다.

지오이드와 기준 타원체 사이의 관계는 단순히 수직 거리의 차이에 그치지 않고, 방향의 차이인 수직선 편차(Deflection of the Vertical)를 동반한다. 수직선 편차는 특정 지점에서 타원체에 수직인 법선(Normal)과 실제 중력 방향인 연직선(Plumb line) 사이의 각도 차이를 의미한다. 이는 지오이드의 경사도를 나타내는 지표가 되며, 지오이드고의 공간적 변화율과 밀접한 관련이 있다. 수직선 편차 $ $는 남북 방향 성분($ $)과 동서 방향 성분($ $)으로 분해되어 표현되며, 이는 지상에서 수행되는 천문 측량 결과와 지심 좌표계 기반의 측량 결과를 연결하는 핵심 고리 역할을 한다.

현대 고도 체계에서 지오이드는 국가 수직 기준면(Vertical Datum)의 근간이 된다. 과거에는 특정 해안의 장기 평균 해수면을 0으로 설정하여 수준 측량을 통해 고도를 전달하였으나, 이는 해양 동역학적 요인에 의한 오차를 포함하는 한계가 있었다. 최근에는 전 지구적 중력장 모델과 고정밀 위성 중력 탐사 데이터를 결합하여 정의된 세계 수직 참조 체계(International Vertical Reference System, IVRS)를 통해 보다 객관적이고 통일된 고도 기준을 확립하려는 노력이 지속되고 있다¹⁶⁾. 이러한 체계 하에서 지오이드는 단순한 가상의 면을 넘어, 범지구적 차원의 지형 변화 관측과 해수면 상승 모니터링을 위한 절대적인 기준틀로서 기능한다.

지오이드고(Geoid height) 또는 지오이드 파고(Geoid undulation)는 기준 타원체(Reference Ellipsoid)의 법선 방향을 따라 타원체면에서 지오이드면까지 측정한 수직 거리를 의미한다. 이는 지구의 물리적 형상과 기하학적 모델 사이의 이격 정도를 나타내는 핵심 지표로, 측지학(Geodesy)에서 고도 체계의 상호 변환과 지구 형상 해석의 근간을 이룬다. 기하학적으로 정의된 타원체면은 질량 분포의 불균일성을 반영하지 못하므로, 실제 중력장을 반영하는 지오이드와의 사이에는 필연적으로 편차가 발생하게 된다.

임의의 지점에서의 타원체 고도(Ellipsoid height) $h$와 정표고(Orthometric height) $H$, 그리고 지오이드고 $N$ 사이에는 다음과 같은 기본적인 기하학적 관계식이 성립한다.

$$h = H + N$$

여기서 $h$는 위성 항법 시스템(Global Navigation Satellite System, GNSS)을 통해 직접적으로 관측되는 기하학적 고도이며, $H$는 지오이드를 기준으로 측정된 물리적 고도이다. 따라서 지오이드고를 정확히 파악하는 것은 GNSS 관측값을 실무에서 사용하는 표고 체계로 변환하기 위한 필수적인 전제 조건이 된다. 이때 $N$의 부호는 지오이드가 타원체보다 위쪽에 위치하면 양(+), 아래쪽에 위치하면 음(-)의 값을 가진다.

지오이드고의 공간적 분포는 지구 내부의 질량 불균형(Mass anomaly) 및 밀도 분포와 밀접하게 연관되어 있다. 지각의 밀도가 평균보다 높거나 맨틀(Mantle) 내에서 상승류가 발생하는 지역에서는 중력 포텐셜이 상대적으로 높게 형성되어 지오이드가 기준 타원체 밖으로 솟아오르는 양의 지오이드고를 나타낸다. 반대로 밀도가 낮거나 하강류가 존재하는 지역에서는 지오이드가 타원체 안쪽으로 함몰되는 음의 지오이드고가 관측된다. 전 지구적으로 지오이드고는 대략 $-107\text{m}$에서 $+85\text{m}$ 사이의 범위를 보이며, 가장 낮은 지점은 인도양 남부에, 가장 높은 지점은 북대서양과 인도네시아 인근 해역에 분포한다¹⁷⁾.

지오이드의 굴곡은 단순한 수직 거리의 차이에 그치지 않고, 수직선 편차(Deflection of the vertical)라는 기하학적 편차를 수반한다. 수직선 편차는 지오이드면의 법선 방향인 실제 중력 방향과 기준 타원체의 법선 방향 사이의 각도 차이를 의미한다. 이는 지오이드고의 수평적 변화율, 즉 지오이드의 경사와 수학적으로 연결된다. 지오이드고의 변화가 급격한 지역일수록 수직선 편차가 크게 나타나며, 이는 천문 측량과 지상 측량 결과 사이의 불일치를 야기하는 원인이 된다.

현대 측지학에서는 지구 중력장 모델(Earth Gravitational Model, EGM)을 통해 전 지구적인 지오이드고 분포를 정밀하게 기술한다. 특히 인공위성 중력 탐사 데이터를 결합한 모델들은 수 센티미터 수준의 정밀도로 지오이드고를 산출하며, 이를 통해 전 지구적 고도 기준면의 통합과 정밀 지도 제작이 가능해졌다. 이러한 지오이드고의 정밀한 결정은 해수면 변동 연구나 지각 변동 모니터링 등 지구 과학 전반의 정량적 분석에 핵심적인 기초 자료를 제공한다.

수직선 편차(Deflection of the Vertical)는 지구상의 특정 지점에서 실제 중력 방향을 나타내는 연직선(Plumb line)과 수학적으로 정의된 기준 타원체(Reference Ellipsoid)의 법선 사이의 각도 차이를 의미한다. 이는 지구의 물리적 형상인 지오이드와 기하학적 모델인 타원체가 일치하지 않기 때문에 발생하며, 근본적으로는 지구 내부 및 지표면의 불균일한 질량 분포에 기인한다. 만약 지구가 밀도가 균일한 완전한 회전 타원체라면 중력의 방향은 타원체 법선과 일치하겠으나, 실제 지구는 지각의 두께 차이, 산맥과 해구 같은 지형적 기복, 그리고 맨틀 내부의 밀도 불균질성으로 인해 복잡한 중력장을 형성한다.

지표면의 거대한 지형 지물은 수직선 편차를 유발하는 가장 직접적인 원인이다. 뉴턴의 만유인력 법칙에 따라, 거대한 산맥은 주변의 연직선을 자신 쪽으로 끌어당기는 효과를 낸다. 예를 들어, 산맥 근처에서 관측한 천문 위도와 측지 위도 사이에는 유의미한 차이가 발생하는데, 이는 산맥의 거대한 질량이 추(plumb bob)를 미세하게 편향시키기 때문이다. 반대로 해구와 같이 질량이 결손된 지역에서는 연직선이 반대 방향으로 밀려나는 현상이 관찰된다. 이러한 지형적 영향은 지형 보정(Topographic reduction)을 통해 정량적으로 분석될 수 있으며, 이는 국지적인 지오이드의 굴곡을 결정하는 핵심 요소가 된다.

그러나 지형적 요인만으로는 실제 관측되는 수직선 편차를 모두 설명할 수 없다. 지구 내부의 아이소스타시(Isostasy, 지각 균형설) 원리에 따르면, 지표의 높은 지형은 지하에 저밀도의 뿌리를 가짐으로써 중력적 평형을 이루려는 경향이 있다. 그럼에도 불구하고 실제 지구는 완전한 평형 상태에 도달하지 못하며, 이러한 잔여 불균형은 수직선 편차의 형태로 나타난다. 또한 맨틀 대류에 의한 물질의 이동이나 지각판의 섭입 등 역동적인 지질학적 과정은 지구 내부의 밀도 경계면을 왜곡시키며, 이는 지표면에서의 중력 방향 변화를 야기한다. 즉, 수직선 편차는 지표 지형뿐만 아니라 심부의 질량 분포 이상을 모두 포괄하는 결과물이다.

수학적으로 수직선 편차는 남북 성분인 $\xi$(Xi)와 동서 성분인 $\eta$(Eta)로 분해하여 표현한다. 천문 좌표($\Phi, \Lambda$)와 측지 좌표($\phi, \lambda$)의 관계를 통해 다음과 같이 정의된다.

$$ \xi = \Phi - \phi $$ $$ \eta = (\Lambda - \lambda) \cos \phi $$

이러한 성분들은 지오이드고(Geoid height) $N$의 경사도와 밀접한 관련이 있다. 지오이드의 경사가 가파를수록 수직선 편차의 크기는 커지며, 이는 다음의 미분 관계식으로 나타낼 수 있다.

$$ \xi = -\frac{1}{R} \frac{\partial N}{\partial \phi} $$ $$ \eta = -\frac{1}{R \cos \phi} \frac{\partial N}{\partial \lambda} $$

여기서 $R$은 지구의 평균 반지름이다. 이는 수직선 편차가 지오이드라는 등전위면의 기하학적 기울기를 물리적으로 나타내는 지표임을 보여준다¹⁸⁾.

결론적으로 수직선 편차는 지구 중력장의 비균질성을 지시하는 중요한 물리량이다. 이는 단순한 측정 오차가 아니라, 지구 내부의 질량 구조와 역학적 상태를 반영하는 정밀한 정보이다. 현대 측지학에서는 Vening Meinesz 공식을 활용하여 전 지구적인 중력 이상 데이터로부터 수직선 편차를 산출하며, 이는 위성 항법 시스템(GNSS) 관측을 통한 타원체 고도를 물리적 고도인 표고로 변환하는 과정에서 필수적인 보정 항목으로 사용된다.

지표상의 높이를 결정하는 문제는 단순한 기하학적 위치 결정을 넘어, 중력이라는 물리적 힘의 방향과 크기를 고려해야 하는 복잡한 과제이다. 인간의 활동과 직결되는 수문학적 현상이나 토목 공학적 설계는 중력의 영향을 받는 물의 흐름에 의존하기 때문에, 순수한 기하학적 거리보다는 물리적 의미를 지닌 고도 체계가 필요하다. 지오이드(Geoid)는 이러한 물리적 고도 체계의 출발점이자 기준면인 수직 기준계(Vertical Datum)로서 핵심적인 역할을 수행한다. 중력 포텐셜이 일정한 등전위면(Equipotential surface)인 지오이드는 이론적으로 물이 정지해 있을 때의 표면과 일치하므로, 모든 고도 측정의 물리적 기점이 된다.

고도 체계에서 가장 널리 사용되는 정표고(Orthometric height)는 지표면의 한 점으로부터 지오이드까지의 거리를 수직선(Plumb line)을 따라 측정한 값으로 정의된다. 이때 수직선은 각 점을 지나는 중력 등전위면에 직교하는 곡선 궤적을 의미한다. 현대의 위성 항법 시스템(Global Navigation Satellite System, GNSS)을 통해 얻어지는 고도는 기준 타원체(Reference Ellipsoid)를 기준으로 하는 타원체고(Ellipsoidal height)이며, 이는 지구 내부의 질량 분포와 무관한 순수 기하학적 수치이다. 따라서 실제 지형에서 물이 흐르는 방향이나 위치 에너지를 파악하기 위해서는 타원체고에서 지오이드고(Geoid height)를 보정하여 정표고를 산출하는 과정이 반드시 수반되어야 한다.

타원체고 $h$, 정표고 $H$, 그리고 지오이드고 $N$ 사이의 관계는 일차적으로 다음과 같은 수식으로 표현된다. $$h \approx H + N$$ 위 식에서 지오이드고는 기준 타원체로부터 지오이드까지의 수직 거리를 의미하며, 지구 내부의 밀도 불균형에 따라 지역마다 서로 다른 값을 가진다. 정밀한 측지학(Geodesy)적 연산에서는 지표면과 지오이드 사이의 구간 평균 중력을 산출하여 정표고를 더욱 엄밀하게 계산한다. 이러한 물리적 고도는 두 지점 사이의 고도 차이가 존재하더라도 중력 포텐셜이 같다면 물이 흐르지 않는다는 물리적 실체를 정확히 반영한다. 만약 지오이드를 기준으로 하지 않고 타원체고만을 기준으로 삼는다면, 기하학적으로는 낮은 곳에서 높은 곳으로 물이 거슬러 올라가는 것처럼 보이는 모순이 발생할 수 있다.

각 국가는 자국 영토 내에서 일관된 고도 값을 유지하기 위해 특정 지점의 평균 해수면(Mean Sea Level)을 고도 0으로 설정하여 관리한다. 대한민국은 인천만의 평균 해수면을 기준으로 수준원점을 설치하여 국가 고도 체계의 기점으로 삼고 있다. 이론적으로 이러한 평균 해수면은 지오이드와 일치해야 하지만, 실제 해양에서는 해류, 기압, 염분 차이 등으로 인해 지오이드와 수 미터 이내의 편차가 발생하는 해면 위상차(Sea Surface Topography)가 존재한다. 이에 따라 최근에는 전 지구적 중력 모델을 기반으로 한 정밀 지오이드를 구축하여, 전통적인 수준 측량(Leveling)의 한계를 극복하고 GNSS를 이용한 효율적인 고도 결정 체계로 전환하려는 노력이 전 세계적으로 진행되고 있다. 이는 국경을 초월한 고도 체계의 통합과 해수면 상승 등 기후 변화 모니터링을 위한 정밀한 기준을 제공한다는 점에서 학문적·실무적 중요성이 매우 크다.

지오이드 모델은 현대 측지학뿐만 아니라 지구물리학, 해양학, 그리고 환경과학 등 다양한 학문적·실무적 영역에서 필수적인 기준체계로 활용된다. 지오이드의 정밀한 결정은 단순히 지구의 형상을 파악하는 것에 그치지 않고, 지표면과 지구 내부에서 발생하는 물리적 현상을 정량적으로 분석하는 토대가 된다.

실무적으로 가장 광범위하게 활용되는 분야는 위성 항법 시스템(Global Navigation Satellite System, GNSS)을 이용한 정밀 측량이다. 인공위성을 통해 얻어지는 고도 정보는 기준 타원체를 기준으로 한 타원체고(Ellipsoidal height, $ h $)이다. 그러나 실제 토목 공사나 지도 제작에서 요구되는 높이는 중력 방향과 일치하는 정표고(Orthometric height, $ H $)이다. 지오이드 모델은 타원체와 지오이드 사이의 거리인 지오이드고(Geoid height, $ N $)를 제공함으로써, 다음과 같은 관계식을 통해 GNSS 관측값을 실용적인 표고로 변환할 수 있게 한다.

$$H = h - N$$

이러한 변환은 국가 수직 기준계를 확립하고, 장거리 터널 굴착이나 교량 건설 시 수평을 유지하는 데 결정적인 역할을 한다. 특히 정밀 지오이드 모델이 구축된 지역에서는 고가의 수준 측량을 대체하여 작업 효율성을 획기적으로 높일 수 있다.

지구물리학적 연구에서 지오이드는 지구 내부의 질량 분포와 역학적 상태를 진단하는 지표가 된다. 지오이드의 요철(Undulation)은 지각과 맨틀 내부의 밀도 불균형을 반영하며, 이는 맨틀 대류의 흐름이나 지각 평형(Isostasy) 상태를 연구하는 데 중요한 정보를 제공한다. 예를 들어, 지오이드가 주변보다 높게 솟아오른 지역은 지하에 고밀도 물질이 존재하거나 상승류가 발생하는 지점으로 해석될 수 있다. 이를 통해 과학자들은 보이지 않는 지구 내부의 구조를 시각화하고 역학적 모델을 검증한다.

해양학 분야에서는 해면 위상차(Sea Surface Topography, SST)를 산출하는 데 지오이드가 핵심적인 기준면이 된다. 실제 해수면은 해류, 기압, 온도 변화 등의 영향으로 지오이드와 일치하지 않고 미세한 편차를 보인다. 위성 고도계(Satellite Altimetry)로 측정한 해수면 높이에서 지오이드고를 제외하면 순수한 해양 역학적 성분만을 추출할 수 있다. 이는 지균류(Geostrophic current)의 유속과 방향을 계산하는 근거가 되며, 전 지구적 해양 순환을 이해하는 데 필수적이다¹⁹⁾.

최근에는 지구 중력장 탐사 위성(Gravity Field and Steady-State Ocean Circulation Explorer, GOCE)이나 중력 회복 및 기후 실험 위성(Gravity Recovery and Climate Experiment, GRACE) 등을 통해 지오이드의 시간적 변화를 관측함으로써 기후 변화 연구에도 기여하고 있다. 빙하의 융해나 지하수 저장량의 변화는 국지적인 중력 변화를 야기하며, 이는 곧 지오이드의 미세한 변동으로 나타난다. 이러한 시계열 분석은 해수면 상승의 원인을 규명하고 지구의 질량 이동을 모니터링하는 정밀한 수단으로 자리 잡았다²⁰⁾.

현대 측지학의 실무적 측면에서 지오이드는 국가 고도 체계의 근간이 되는 수직 기준(Vertical Datum)을 확립하는 데 결정적인 역할을 수행한다. 전통적으로 지표면의 높이를 결정하는 수준 측량(Levelling)은 평균 해수면을 연장한 가상의 면을 기준으로 수행되어 왔으나, 물리적으로 이는 중력의 등전위면인 지오이드와 긴밀하게 연결되어 있다. 국가 지도 제작과 대규모 토목 공학 프로젝트에서 정밀한 수직 위치 정보를 확보하기 위해서는 지구 타원체와 실질적인 물리적 높이 사이의 관계를 규명하는 것이 필수적이다.

위성 항법 시스템(GNSS)의 보급은 측량 패러다임을 혁신하였으나, GNSS가 제공하는 높이 정보는 수학적으로 정의된 지구 타원체로부터의 거리인 타원체고(Ellipsoidal height, $h$)이다. 반면, 실제 지형의 기복이나 물의 흐름을 결정하는 높이는 지오이드로부터 측정한 표고(Orthometric height, $H$)이다. 이 두 값 사이의 차이를 지오이드고(Geoid height, $N$)라 하며, 다음과 같은 기본적인 기하학적 관계식이 성립한다.

$$ H = h - N $$

따라서 GNSS를 이용해 정밀한 표고를 산출하기 위해서는 해당 지역의 정밀한 지오이드 모델이 반드시 선행되어야 한다. 국가 단위의 지도 제작 시에는 지상 중력 데이터, 위성 중력 데이터, 그리고 지형 모델을 결합하여 구축한 물리적 지오이드에 국가기준점에서의 GNSS/수준 측량 성과를 통합한 하이브리드 지오이드 모델(Hybrid Geoid Model)을 주로 사용한다. 대한민국에서는 국토지리정보원을 중심으로 한국형 지오이드 모델(KNGeoid)을 지속적으로 갱신하고 있으며, 최근의 모델들은 전국적으로 수 센티미터(cm) 수준의 정확도를 확보하여 직접 수준 측량을 상당 부분 대체하고 있다²¹⁾.

정밀 지오이드 모델은 대규모 토목 공사에서 구조물의 수직적 안정성을 확보하는 데 필수적이다. 교량, 댐, 터널과 같은 대형 구조물은 미세한 고도 차이에 의해 설계 하중이나 배수 체계가 결정되므로, 광범위한 지역에 걸쳐 일관된 수직 기준을 제공하는 지오이드의 역할이 중요하다. 특히 국가 간 혹은 지역 간 고도 체계의 불일치를 해소하고 통합된 공간 정보를 구축하기 위해, 전 지구적 중력장 모델을 기반으로 한 국지적 지오이드의 정밀화 작업이 전 세계적으로 진행되고 있다. 이러한 정밀 측량 기술의 발전은 단순한 지도 제작을 넘어 디지털 트윈 구축이나 자율 주행을 위한 고정밀 도로 지도 제작 등 미래 산업의 핵심 인프라로 기능한다.

지오이드의 기복, 즉 지오이드고(Geoid height)는 지구 내부의 질량 분포가 불균일하다는 것을 보여주는 직접적인 증거이다. 지구물리학(Geophysics)적 관점에서 지오이드의 굴곡을 해석하는 것은 지표면에서 관측할 수 없는 지구 심부의 밀도(Density) 구조와 동역학적 과정을 이해하는 핵심적인 수단이 된다. 지오이드의 형태는 지구 내부의 정적 밀도 분포뿐만 아니라, 맨틀 대류(Mantle convection)와 같은 동적인 흐름에 의한 지표 및 내부 경계면의 변형에 의해 결정된다.

지각 평형설(Isostasy)은 지각과 상부 맨틀의 밀도 불균형이 지오이드에 미치는 영향을 설명하는 중요한 이론적 기초를 제공한다. 에어리(Airy) 모델에 따르면, 산맥과 같이 지표 질량이 과잉된 지역은 그 하부에 저밀도의 지각 뿌리가 발달하여 부력을 얻는다. 이때 지표의 질량 과잉은 지오이드를 상승시키려 하지만, 하부의 밀도 결손은 지오이드를 하강시키려는 상반된 효과를 가진다. 지오이드고 $N$은 중력 포텐셜의 섭동인 교란 포텐셜(Disturbing potential) $T$와 기준 중력 $\gamma$의 관계로 나타나는 브룬스 공식(Bruns’ formula)에 의해 다음과 같이 정의된다.

$$ N = \frac{T}{\gamma} $$

지각 평형이 완전히 이루어진 지역에서 단파장(Short-wavelength)의 지오이드 굴곡은 지각의 두께와 밀도 차이에 민감하게 반응하며, 이를 통해 모호 불연속면(Mohorovičić discontinuity)의 깊이를 추정할 수 있다.

반면, 전 지구적 규모에서 나타나는 장파장(Long-wavelength)의 지오이드 이상은 주로 맨틀 심부의 밀도 불균형과 관련이 있다. 판 구조론(Plate tectonics)에 의해 섭입된 차가운 판은 주위 맨틀보다 밀도가 높아 중력을 증가시키고 지오이드를 상승시키는 요인이 된다. 그러나 이러한 고밀도 물질의 하강은 맨틀의 점성 흐름을 유도하여 지구 표면을 아래로 끌어당기는 동적 지형(Dynamic topography)을 형성한다. 지표면의 하강은 질량 결손을 의미하므로 지오이드를 낮추는 효과를 낸다. 따라서 실제 관측되는 지오이드는 내부 밀도 이상에 의한 직접적인 중력 효과와 동적 지형에 의한 중력 효과의 합으로 결정된다.

맨틀의 점성 구조에 따라 이 두 효과의 상대적 크기가 달라지는데, 이는 지오이드 데이터를 역산하여 맨틀의 점성도(Viscosity) 층서 구조를 규명하는 근거가 된다. 예를 들어, 하부 맨틀의 점성도가 상부 맨틀보다 현저히 높을 경우, 내부 밀도 이상에 의한 효과가 지표 변형 효과를 압도하여 고밀도 지역에서 지오이드가 상승하게 된다. 이러한 연구는 지진파 토모그래피(Seismic tomography)를 통해 얻어진 지구 내부 영상과 결합되어, 지구 내부의 열적 상태와 물질 순환 모델을 검증하는 강력한 제약 조건으로 작용한다²²⁾.

또한, 지오이드의 형태는 핵-맨틀 경계(Core-Mantle Boundary, CMB)의 기복을 파악하는 데에도 기여한다. 외핵과 맨틀 사이의 밀도 차이는 지각과 대기 사이의 밀도 차이보다 훨씬 크기 때문에, 이 경계면의 미세한 변형도 장파장 지오이드에 유의미한 영향을 미친다. 연구자들은 구면 조화 함수(Spherical harmonics)를 이용하여 지오이드 데이터를 차수별로 분리하고, 이를 각 깊이별 밀도 불균형 모델과 비교함으로써 지구 내부의 층상 구조별 기여도를 정량화한다²³⁾. 결과적으로 지오이드 연구는 단순한 지형 측정을 넘어, 거대한 열기관으로서 작동하는 지구의 내부 동력학을 규명하는 필수적인 학문적 도구라 할 수 있다.

지오이드는 지구의 평균 해수면과 일치하는 중력 등전위면으로 정의되지만, 실제 해수면은 해류, 기압 변화, 온도 및 염분 분포의 불균형으로 인해 지오이드로부터 일정 부분 벗어나게 된다. 이러한 물리적 편차를 해수면 지형(Sea Surface Topography, SST) 또는 해양 동역학적 지형(Ocean Dynamic Topography)이라 한다. 해양학과 측지학의 결합을 통해 도출되는 해수면 지형은 해양의 대순환을 파악하는 핵심 지표가 된다. 특정 지점의 실제 해수면 고도(Sea Surface Height, SSH)와 지오이드고(Geoid Height, $N$)의 관계를 수식으로 표현하면 다음과 같다.

$$ \zeta = SSH - N $$

여기서 $\zeta$는 해수면 지형을 의미한다. 위성 고도계(Satellite Altimetry)를 통해 정밀하게 측정된 해수면 고도에서 중력 관측 위성인 GOCE(Gravity field and steady-state Ocean Circulation Explorer) 등으로 구축된 고정밀 지오이드 모델을 감산하면, 순수한 해양 동역학적 요인에 의한 해수면 높이 변화를 산출할 수 있다. 이는 지균류(Geostrophic current)의 속도와 방향을 계산하는 기초 자료가 되며, 전 지구적 해양 열 수송 및 에너지 분배 과정을 이해하는 데 결정적인 역할을 한다²⁴⁾.

현대 기후 변화 연구에서 주목하는 또 다른 측면은 지오이드의 시간적 가변성이다. 지오이드는 고정불변의 면이 아니라 지구 내부 및 표면의 질량 재분배에 따라 미세하게 변화한다. 특히 기후 변화로 인한 대규모 빙하 융해와 육상의 지하수 저장량 변화는 지구 중력장의 국지적·광역적 변동을 야기한다. 2002년 발사된 GRACE(Gravity Recovery and Climate Experiment) 위성 임무는 이러한 중력장의 시계열 변화를 관측함으로써 지오이드의 변동을 추적하는 데 혁신적인 기여를 하였다. GRACE 데이터는 남극과 그린란드의 빙하 질량 손실을 정량화하고, 이것이 전 지구적 해수면 상승에 미치는 영향을 분석하는 강력한 도구로 활용된다²⁵⁾.

지오이드의 변동 관측은 해수면 상승의 원인을 규명하는 데 있어서도 필수적이다. 해수면 상승은 크게 바닷물의 열팽창에 의한 스테릭 해수면 상승(Steric sea level rise)과 빙하 융해물 유입 등에 의한 바리스태틱 해수면 상승(Barystatic sea level rise)으로 구분된다. 위성 고도계가 전체적인 해수면 높이의 변화를 측정한다면, GRACE와 같은 중력 탐사 위성은 질량 변화에 의한 바리스태틱 성분만을 분리하여 측정할 수 있다²⁶⁾. 이 두 관측치를 비교 분석함으로써 기후 시스템 내에서의 물 순환 경로와 해양 저장 에너지의 변화를 보다 정밀하게 파악하는 것이 가능하다.

결론적으로 지오이드 모델의 정밀화와 시간 가변적 지오이드의 지속적인 모니터링은 지구 시스템 과학(Earth System Science)의 핵심 과제이다. 이는 단순한 고도 기준의 설정을 넘어, 해양 순환의 변동성 파악, 빙하 질량 수지 분석, 그리고 장기적인 기후 예측 모델의 신뢰도 향상을 위한 필수적인 물리적 토대를 제공한다.