문서의 이전 판입니다!
지구의 형상은 지표면의 복잡한 고저 차이와 내부 질량 분포의 불균일성으로 인해 기하학적으로 단순하게 정의하기 어렵다. 그러나 정밀한 위치 결정과 지도 제작을 위해서는 지구의 형상을 수학적으로 다룰 수 있는 모델이 필수적이다. 이러한 필요에 따라 도입된 지구 타원체는 지구의 물리적 실체인 지오이드(Geoid)를 기하학적으로 가장 가깝게 근사한 회전 타원체(Ellipsoid of revolution)로 정의된다. 이는 지구의 자전축을 회전축으로 하여 타원을 회전시킨 입체 형상을 의미하며, 지구의 전체적인 기하학적 특성을 대표하는 표준 모델로 기능한다.
지구가 완벽한 구형이 아닌 타원체의 형상을 띠게 된 물리적 근거는 지구의 자전에 기인한다. 아이작 뉴턴(Isaac Newton)은 지구가 자전함에 따라 발생하는 원심력(Centrifugal force)이 지구 구성 물질을 적도 방향으로 밀어낼 것이라고 예측하였다. 자전축으로부터 수직 거리 $ r $만큼 떨어진 지점에서 각속도 $ $로 회전하는 물체가 받는 원심력 가속도 $ a_c $는 다음과 같이 표현된다.
$$ a_c = \omega^2 r $$
이 수식에 따르면 원심력은 자전축으로부터의 거리 $ r $에 비례하므로, 적도 지역에서 최대가 되고 양극점에서는 0이 된다. 결과적으로 지구의 모든 지점에서 작용하는 중력(Gravity)은 지구 질량에 의한 만유인력과 자전에 의한 원심력의 벡터 합으로 나타나며, 이는 적도 부근을 팽창시키는 동역학적 요인으로 작용한다.
지구 형성 초기, 지구가 용융 상태의 유체였다고 가정하면 지구는 유체 정역학적 평형(Hydrostatic equilibrium) 상태에 도달하게 된다. 유체 정역학적 평형 상태의 회전체는 내부의 압력 구배와 원심력, 만유인력이 평형을 이루어 표면의 모든 지점에서 중력 포텐셜(Potential)이 일정해지는 형태를 갖춘다. 이러한 물리적 메커니즘에 의해 지구는 적도 반지름이 극 반지름보다 긴 편평 타원체(Oblate spheroid)의 형상을 확립하게 되었다.
지구 타원체의 기하학적 특성을 결정하는 핵심 요소는 적도 반지름인 장반경(Semimajor axis) $ a $와 극 반지름인 단반경(Semiminor axis) $ b $이다. 이 두 반지름의 차이를 장반경으로 나눈 값인 편평률(Flattening, $ f $)은 지구의 찌그러진 정도를 정량적으로 나타내는 지표가 된다.
$$ f = \frac{a - b}{a} $$
현대 측지학에서는 인공위성 관측 데이터와 중력 측정을 바탕으로 이 수치들을 정밀하게 산출한다. 지구 타원체는 단순히 기하학적 형태를 정의하는 데 그치지 않고, 지구 내부의 밀도 분포와 회전 속도 등 물리적 상수를 포함하는 지구 참조 시스템(Earth Reference System)의 기초가 된다. 따라서 지구 타원체의 정의와 그 물리적 근거에 대한 이해는 위성 항법 시스템(Global Navigation Satellite System, GNSS)이나 정밀 지도 제작과 같은 현대 과학 기술의 논리적 토대를 형성한다.
지구의 기하학적 형상은 단순히 고정된 고체의 형태가 아니라, 자전(rotation)이라는 역학적 과정에 의해 결정되는 유동적인 평형 상태의 결과물이다. 지구가 자전함에 따라 발생하는 원심력(centrifugal force)은 적도 부근에서 가장 강하게 작용하며, 이는 지구 내부의 중력(gravity)과 상호작용하여 적도 방향으로의 질량 팽창을 유도한다. 그 결과 지구는 완벽한 구(sphere)가 아닌, 적도 반지름이 극 반지름보다 더 긴 편평 타원체(oblate spheroid)의 형상을 띠게 된다. 이러한 기하학적 특성은 지구 물리학적 관측과 측지학(geodesy)적 계산의 기초가 되며, 지구 표면의 위치를 정밀하게 규정하기 위한 수학적 모델의 도입을 필연적으로 요구한다.
기하학적 관점에서 지구 타원체는 타원을 그 단축(minor axis)을 회전축으로 하여 회전시켰을 때 생성되는 회전 타원체(ellipsoid of revolution)로 정의된다. 이를 수학적으로 표현하기 위해 3차원 직교 좌표계(Cartesian coordinate system)에서의 타원체 방정식을 도입할 수 있다. 타원체의 중심을 원점으로 하고, 회전축을 $z$축, 적도면을 $xy$평면으로 설정할 때, 타원체의 표면을 이루는 점 $(x, y, z)$들의 집합은 다음과 같은 방정식을 만족한다.
$$ \frac{x^2 + y^2}{a^2} + \frac{z^2}{b^2} = 1 $$
여기서 $a$는 타원체의 장반경(semi-major axis)으로 적도 반지름에 해당하며, $b$는 단반경(semi-minor axis)으로 극 반지름에 해당한다. 지구의 경우 $a > b$인 관계가 성립하며, 이 두 매개변수의 차이는 지구가 얼마나 평평한지를 나타내는 편평률(flattening, $f$)로 정량화된다. 편평률은 다음과 같이 정의된다.
$$ f = \frac{a - b}{a} $$
이 수치는 지구의 물리적 특성을 반영하는 핵심 상수로, 측지계(geodetic system)를 설계할 때 반드시 결정되어야 하는 기하학적 요소이다. 예를 들어 현대 측지학의 기준이 되는 세계 측지 시스템 1984(World Geodetic System 1984, WGS 84)에서는 장반경 $a$를 6,378,137.0미터로, 편평률 $f$의 역수를 298.257223563으로 정의하고 있다1).
실제 지구의 표면은 산맥, 해구, 그리고 밀도 불균형에 따른 중력 변화로 인해 매우 불규칙한 형상을 지닌다. 이처럼 복잡한 물리적 표면을 직접적인 계산의 기준으로 삼기에는 수학적 취급이 극도로 난해하다. 따라서 측지학에서는 실제 지구의 형상인 지오이드(geoid)에 가장 근사하면서도 기하학적으로 명확하게 정의되는 타원체를 설정하여 이를 수평 위치 결정의 기준으로 삼는다. 이러한 수학적 추상화는 지도 투영(map projection), 항법 시스템(navigation system), 지리 정보 시스템(Geographic Information System, GIS) 등에서 위치 좌표를 일관되게 정의하고 변환하는 데 필수적인 역할을 수행한다. 특히 현대의 위성 항법 시스템(Global Navigation Satellite System, GNSS)은 지구 중심을 원점으로 하는 세계 표준 타원체 모델을 사용하여 지구상의 모든 지점을 3차원 좌표로 정밀하게 산출한다.
지구가 완전한 구형을 이루지 않고 적도 부근이 부풀어 오른 회전 타원체(oblate spheroid)의 형상을 갖게 된 근본적인 원인은 지구의 자전에 따른 역학적 상호작용에 있다. 지구상의 모든 질점에는 두 가지 주요한 힘이 작용하는데, 하나는 지구의 질량에 의해 중심 방향으로 끌어당기는 만유인력이고, 다른 하나는 자전축을 중심으로 회전함에 따라 바깥쪽으로 밀어내는 원심력이다. 이 두 힘의 벡터 합을 중력이라 정의하며, 지구의 물리적 형상은 바로 이 중력이 지표면의 모든 지점에서 수직으로 작용하는 평형 상태를 찾아가는 과정에서 결정된다.
자전에 의한 원심력 $ _c $는 다음과 같은 수식으로 표현된다.
$$ \mathbf{F}_c = m \omega^2 \mathbf{r} $$
여기서 $ m $은 물체의 질량, $ $는 지구 자전의 각속도, $ $은 자전축에서 해당 지점까지의 수직 거리 벡터이다. 원심력의 크기는 자전축으로부터의 거리 $ r $에 비례하므로, 자전축 위에 위치한 양 극점에서는 0이 되며, 자전축에서 가장 멀리 떨어진 적도 지역에서 최대가 된다. 반면 만유인력은 지구 중심을 향해 거의 일정하게 작용하므로, 적도 지역에서는 원심력이 만유인력의 일부를 상쇄하여 극 지역보다 상대적으로 약한 중력이 형성된다.
지구의 형상을 결정짓는 핵심적인 물리 원리는 정역학적 평형(hydrostatic equilibrium)이다. 태초에 고온의 유체 상태였던 지구는 중력과 원심력의 영향 아래 물질이 재배치되는 과정을 거쳤다. 유체 내부의 압력 경사력과 외부의 중력이 균형을 이룰 때, 유체의 표면은 중력 포텐셜이 일정한 등포텐셜면을 형성하게 된다. 만약 지구가 자전하지 않는다면 등포텐셜면은 완전한 구형을 이루겠지만, 자전으로 인한 원심력이 가해지면 적도 방향으로 물질이 쏠리게 되어 적도 반지름이 극 반지름보다 긴 타원체 형상에서 평형을 이루게 된다.
이러한 물리적 관계는 알렉시 클로드 클레로(Alexis Claude Clairaut)가 제시한 클레로의 정리(Clairaut’s theorem)를 통해 수학적으로 구체화되었다. 클레로의 정리는 지구의 편평도와 위도에 따른 중력의 변화 사이의 관계를 설명하며, 지구가 내부적으로 밀도 층을 이루고 있는 회전 유체라고 가정할 때 외부 형상이 어떻게 결정되는지를 보여준다. 이 이론에 따르면, 지구의 회전 속도와 내부 질량 분포가 주어졌을 때 정역학적 평형을 만족하는 표면의 기하학적 형태는 특정 이심률을 가진 타원체로 수렴하게 된다.
결과적으로 지구는 자전이라는 동역학적 요인과 만유인력이라는 정역학적 요인이 결합하여, 적도 반지름이 극 반지름보다 약 21km 더 긴 지구 타원체의 모습을 띠게 되었다. 이는 단순히 기하학적인 가정이 아니라, 행성의 회전과 질량 분포가 만들어내는 에너지 최소화 상태의 물리적 결과물이다. 이러한 형성 원리는 지구뿐만 아니라 자전하는 모든 천체에 보편적으로 적용되는 물리 법칙의 산물이라 할 수 있다.2)
인류가 거주하는 터전인 지구의 형상을 파악하려는 시도는 고대부터 지속되었다. 초기 인류는 직관적인 관찰을 바탕으로 지구가 평면이라는 인식을 가졌으나, 천문 관측과 항해 경험이 축적되면서 구형 지구(Spherical Earth) 가설이 대두되었다. 피타고라스(Pythagoras)는 기하학적 완벽함을 근거로 지구가 구형일 것이라 추론하였으며, 아리스토텔레스(Aristotle)는 월식 때 달에 비친 지구의 그림자가 곡선이라는 점과 남북으로 이동할 때 보이는 별의 고도가 달라진다는 사실을 들어 이를 과학적으로 논증하였다. 이후 에라토스테네스(Eratosthenes)는 하짓날 시에네(Syene)와 알렉산드리아(Alexandria)에서의 태양 남중 고도 차이를 이용하여 지구의 둘레를 최초로 수치화하였다. 그는 두 지점 사이의 거리 $ s $와 중심각 $ $ 사이의 비례 관계를 다음과 같이 이용하였다.
$$ \frac{\theta}{360^\circ} = \frac{s}{C} $$
여기서 $ C $는 지구의 전체 둘레를 의미한다. 비록 측정상의 오차는 존재하였으나, 이는 지구의 크기를 수학적으로 산출한 최초의 정량적 시도로 평가받는다.
17세기에 이르러 아이작 뉴턴(Isaac Newton)은 자신의 저서인 자연철학의 수학적 원리(Philosophiae Naturalis Principia Mathematica)에서 역학적 원리를 바탕으로 지구 형상에 대한 새로운 이론을 제시하였다. 뉴턴은 지구의 자전으로 인해 발생하는 원심력(Centrifugal force)이 만유인력(Universal gravitation)과 상호작용하여, 유체 상태였던 초기 지구가 적도 방향으로 부풀어 오른 편평 타원체(Oblate spheroid)가 되었을 것이라고 예측하였다. 그는 지구가 균질한 회전 유체라고 가정하고, 적도 반지름이 극 반지름보다 약 230분의 1만큼 더 길어야 한다는 결론을 도출하였다. 반면 프랑스의 조반니 도메니코 카시니(Giovanni Domenico Cassini)와 자크 카시니(Jacques Cassini) 부자는 프랑스 내 자오선 호의 길이를 측정한 데이터를 근거로, 북쪽으로 갈수록 위도 1도에 해당하는 거리가 짧아진다는 점을 들어 지구가 남북으로 긴 장구 타원체(Prolate spheroid)라고 주장하였다.
이른바 영국과 프랑스의 과학 논쟁으로 비화한 이 대립을 해결하기 위해 프랑스 과학 아카데미(French Academy of Sciences)는 대규모 측량 원정대를 파견하였다. 1735년 피에르 부게(Pierre Bouguer)와 샤를 마리 드 라 콩다민(Charles Marie de La Condamine)이 이끄는 원정대가 적도 인근인 페루(Peru, 현재의 에콰도르)로 향하였고, 1736년에는 피에르 루이 모페르튀이(Pierre Louis Maupertuis)가 이끄는 원정대가 북극권인 라플란드(Lapland)로 파견되었다. 이들의 임무는 서로 다른 위도에서 자오선(Meridian) 호 1도의 실제 지표 거리를 정밀하게 측정하는 것이었다. 만약 뉴턴의 주장이 옳다면 고위도로 갈수록 곡률이 작아져 1도의 거리가 길어질 것이고, 카시니의 주장이 옳다면 그 반대의 결과가 나타날 것이었다.
원정대의 측정 결과, 라플란드에서의 1도 호의 길이가 페루에서의 길이보다 더 길다는 사실이 확정적으로 밝혀졌다. 이는 위도가 높아질수록 지표면이 더 평평해짐을 의미하며, 지구가 자전의 영향으로 형성된 편평 타원체임을 입증하는 결정적인 증거가 되었다. 모페르튀이는 이 공로로 ’지구를 납작하게 만든 사람’이라는 별칭을 얻기도 하였다. 이 역사적 과정은 단순한 형상 확인을 넘어, 고전 역학의 보편성을 검증하고 현대 측지학(Geodesy)이 정밀 과학으로 도약하는 계기가 되었다. 이후 과학계의 관심은 지구가 완벽한 타원체에서 얼마나 벗어나 있는지를 측정하는 문제와, 실제 중력 방향에 수직인 등포텐셜 면인 지오이드(Geoid)를 규명하는 방향으로 진화하였다. 3)
인류가 거주하는 지표면의 전체적인 형상에 대한 인식은 신화적 상상력에서 출발하여 점차 기하학적 추론과 실증적 관측의 영역으로 진화하였다. 초기 고대 문명권에서는 직관적인 경험에 의존하여 대지를 평평한 원반이나 사각형의 형태로 간주하였으나, 해양 항해의 확대와 천문 관측 기술의 발달은 지구가 구형(Spherical shape)이라는 구형설을 태동시켰다. 피타고라스(Pythagoras)는 만물의 근원을 수와 기하학적 조화에서 찾으며 구를 가장 완벽한 입체로 간주하였고, 이러한 철학적 신념을 바탕으로 지구 구형설을 최초로 주장한 것으로 알려져 있다. 이후 아리스토텔레스(Aristotle)는 경험론적 근거를 제시하며 이 가설을 과학적 이론의 단계로 격상시켰다. 그는 월식 때 달에 비치는 지구의 그림자가 항상 둥글다는 점, 북쪽이나 남쪽으로 이동함에 따라 지평선 위로 보이는 별의 고도가 변하거나 새로운 별이 나타난다는 점 등을 근거로 지구가 구체임을 역설하였다.
지구의 크기를 수학적으로 산출하려는 최초의 과학적 시도는 헬레니즘 시대의 학자 에라토스테네스(Eratosthenes)에 의해 이루어졌다. 그는 하절기 하지 정오에 이집트의 시에네(Syene, 현재의 아스완)에서는 햇빛이 깊은 우물 바닥까지 수직으로 비치지만, 북쪽으로 떨어진 알렉산드리아(Alexandria)에서는 막대기 그림자가 생겨 태양의 천정각(Zenith distance)이 발생한다는 사실에 주목하였다. 에라토스테네스는 지구가 완전한 구형이고 태양광선이 평행하게 입사한다고 가정하였으며, 두 지점 사이의 거리 $s$와 중심각 $\theta$ 사이의 관계를 이용하여 지구의 전체 둘레 $C$를 계산하였다.
두 지점 사이의 원호의 길이 $s$는 지구 반지름 $R$과 중심각 $\theta$에 대하여 다음과 같은 관계를 갖는다. $$ s = R\theta $$ 알렉산드리아에서 측정한 태양의 천정각 $\theta$가 원둘레의 50분의 1인 약 $7.2^{\circ}$였고, 두 도시 사이의 거리가 약 5,000스타디아(Stadia)임을 이용하여 그는 지구의 둘레를 약 250,000스타디아로 추정하였다. 비록 당시의 거리 측정 단위인 스타디아의 현대적 환산치에 대해서는 학설이 분분하나, 그의 측정값은 현대의 정밀 관측값과 비교했을 때 오차 범위가 약 1~15% 이내에 불과할 정도로 놀라운 정확성을 보여주었다. 이는 측지학(Geodesy)의 역사에서 기하학적 원리를 실물 지구에 적용한 기념비적인 사건으로 평가받는다.
로마 제국 쇠퇴 이후 유럽의 측지학적 탐구는 정체되었으나, 고대 그리스의 지식은 이슬람 세계로 전승되어 더욱 정교해졌다. 9세기 아바스 왕조의 칼리프 알 마문(Al-Ma’mun)은 고대 문헌의 기록을 검증하기 위해 시난(Sinjar) 평원에서 대규모 자오선 측정을 명령하였다. 이슬람 천문학자들은 위도 1도에 해당하는 자오선 호의 길이를 직접 보행하며 측정하였으며, 이를 통해 지구의 크기를 산출하였다. 이 시기의 연구는 프톨레마이오스(Ptolemy)의 우주관을 계승하면서도 실측을 통해 수치를 보정하려는 시도가 돋보였으며, 이후 대항해 시대에 유럽으로 재유입되어 지구 형상에 대한 근대적 인식의 토대가 되었다. 고전적 지구 형상론은 지구가 단순한 구체라는 전제하에 전개되었으나, 이러한 초기 측정 노력은 17세기 후반 아이작 뉴턴(Isaac Newton)에 의해 지구가 자전의 영향으로 편평해진 타원체일 것이라는 역학적 가설이 등장하기 전까지 인류의 표준적인 지구관으로 자리 잡았다.
지구의 형상이 완전한 구형이 아니라는 의구심은 17세기 후반 아이작 뉴턴(Isaac Newton)의 역학적 추론과 지오반니 도메니코 카시니(Giovanni Domenico Cassini)를 필두로 한 프랑스 천문학자들의 실측 결과가 충돌하며 과학사의 중대한 논쟁으로 부상하였다. 이른바 ’뉴턴-카시니 논쟁’은 단순히 지구의 기하학적 형태를 결정하는 문제를 넘어, 뉴턴이 제시한 고전 역학의 보편적 원리가 실제 지구라는 거대 물리계에 적용될 수 있는지를 검증하는 시금석이 되었다.
아이작 뉴턴은 1687년 저서 자연철학의 수학적 원리(Philosophiae Naturalis Principia Mathematica)에서 지구가 자전함에 따라 발생하는 원심력의 영향을 이론적으로 고찰하였다. 뉴턴은 지구가 과거에 유체 상태였다고 가정할 때, 자전축에 수직인 방향으로 작용하는 원심력이 적도 지역에서 가장 강하고 극 지점에서 영(0)이 된다는 점에 주목하였다. 이러한 힘의 불균형은 적도 부근을 부풀게 하고 극 방향을 수축시켜, 지구를 편평 타원체(Oblate Spheroid)의 형상으로 만들 것이라고 예측하였다. 뉴턴은 중력과 원심력의 평형 관계를 이용하여 지구의 편평도(Flattening)를 약 1/230로 계산하였으며, 이는 위도가 높아질수록 자오선의 곡률 반지름이 커져야 함을 의미하였다. 즉, 뉴턴의 가설이 타당하다면 고위도 지역에서 위도 1도에 해당하는 자오선 호의 길이는 저위도 지역보다 길어야 한다.
반면, 프랑스 왕립 천문대의 초대 대장이었던 지오반니 도메니코 카시니와 그의 아들 자크 카시니(Jacques Cassini)는 실제 지표면 측량 데이터를 근거로 뉴턴의 주장에 반박하였다. 카시니 가문은 1683년부터 1718년까지 프랑스를 남북으로 가로지르는 삼각 측량을 수행하였다. 그들은 파리를 중심으로 북쪽의 됭케르크와 남쪽의 콜리우르를 잇는 자오선 호의 길이를 측정하였는데, 측정 결과 위도가 높아질수록 1도 사이의 거리가 오히려 짧아진다는 결론에 도달하였다. 이는 지구가 적도보다 극 방향으로 더 길쭉한 장구 타원체(Prolate Spheroid)임을 시사하는 것이었다. 카시니의 측정치는 당시 프랑스 과학계의 실증주의적 전통과 결합하여, 영국 측의 이론적 가설인 뉴턴 역학에 대한 강력한 도전으로 받아들여졌다.
두 주장 사이의 모순은 측지학적 관점에서 곡률의 정의와 밀접하게 연관된다. 지구가 편평 타원체라면 고위도로 갈수록 표면이 완만해지므로 곡률 반지름이 커지게 되고, 결과적으로 동일한 각도에 대응하는 호의 길이는 늘어난다. 반대로 장구 타원체라면 고위도에서 표면이 더 급격하게 휘어지므로 곡률 반지름이 작아지고 호의 길이도 짧아진다. 따라서 자오선 호의 길이를 정밀하게 측정하는 것은 지구의 형상을 판별하는 결정적인 실험적 증거가 되었다.
이 논쟁은 1730년대에 이르러 프랑스 과학 아카데미(Académie des Sciences)가 주도한 대규모 원정 측정을 통해 종지부를 찍게 된다. 아카데미는 위도에 따른 호 길이의 차이를 명확히 확인하기 위해 적도 인근의 페루(현재의 에콰도르)와 북극권인 라플란드로 각각 원정대를 파견하였다. 피에르 루이 모페르튀(Pierre Louis Maupertuis)가 이끄는 라플란드 원정대와 피에르 부게(Pierre Bouguer) 등이 참여한 페루 원정대의 측정 결과, 고위도인 라플란드의 자오선 1도 길이가 프랑스나 적도 지역보다 길다는 사실이 최종적으로 입증되었다. 이로써 카시니의 측정 오류가 밝혀졌고, 지구가 편평 타원체라는 뉴턴의 역학적 예측은 실증적으로 확립되었다. 이 논쟁의 종결은 만유인력의 법칙이 지구 규모의 형상을 결정하는 근본 원리임을 증명하였으며, 현대 측지학과 지구 물리 모델링의 기초를 마련하는 계기가 되었다.
프랑스 과학 아카데미의 원정 측정을 통해 지구가 편평 타원체임이 증명되는 과정을 설명한다.
지구의 형상을 기하학적으로 정의하기 위해서는 회전 타원체(ellipsoid of revolution)를 규정하는 독립적인 매개변수가 필요하다. 일반적으로 측지학에서는 타원체의 크기를 결정하는 장반경(semi-major axis) $a$와 형상을 결정하는 편평률(flattening) $f$를 기본 상수로 채택한다. 장반경은 타원체의 중심에서 적도까지의 거리를 의미하며, 단반경(semi-minor axis) $b$는 중심에서 극점까지의 거리를 나타낸다. 이들 사이의 관계로부터 도출되는 편평률은 다음과 같이 정의된다.
$$ f = \frac{a - b}{a} $$
편평률은 지구가 자전에 의한 원심력으로 인해 적도 방향으로 얼마나 팽창했는지를 보여주는 무차원 척도이다. 이와 함께 타원의 기하학적 특성을 나타내는 이심률(eccentricity) 또한 중요한 요소이다. 측지 계산에서는 주로 제1 이심률의 제곱($e^2$)과 제2 이심률의 제곱(${e'}^2$)이 빈번하게 사용되며, 장반경 및 단반경과의 관계식은 다음과 같다.
$$ e^2 = \frac{a^2 - b^2}{a^2} = 2f - f^2, \quad {e'}^2 = \frac{a^2 - b^2}{b^2} = \frac{e^2}{1 - e^2} $$
타원체 표면상의 특정 지점에서 정의되는 곡률은 해당 지점의 위도에 따라 달라진다. 이를 기술하기 위해 자오선 곡률 반경(meridian radius of curvature) $M$과 유부선 곡률 반경(prime vertical radius of curvature) $N$을 정의한다. 자오선 곡률 반경은 남북 방향의 곡률을, 유부선 곡률 반경은 동서 방향의 곡률을 의미한다. 임의의 지리학적 위도 $\phi$에서 각 곡률 반경은 다음과 같은 수식으로 산출된다.
$$ M = \frac{a(1 - e^2)}{(1 - e^2 \sin^2 \phi)^{3/2}} $$ $$ N = \frac{a}{\sqrt{1 - e^2 \sin^2 \phi}} $$
여기서 유부선 곡률 반경 $N$은 타원체 표면의 한 점에서 법선이 자전축과 만나는 지점까지의 거리를 의미하며, 이는 측지 좌표를 3차원 직교 좌표계로 변환할 때 핵심적인 기하학적 인자로 작용한다. 타원체상의 한 점에서 모든 방향에 대한 평균적인 곡률은 위 두 반경의 기하평균인 $ $으로 정의되기도 한다.
현대 측지학에서 사용되는 수치 모델은 인공위성 관측과 중력 측정을 통해 결정된 정밀한 상수값들을 제공한다. 대표적인 모델인 지구 참조 시스템 80(Geodetic Reference System 1980, GRS80)은 국제측지학 및 지구물리학 연맹(IUGG)에서 채택한 표준으로, 장반경 $a = 6,378,137.0 \, \text{m}$와 역편평률 $1/f = 298.257222101$을 기본 상수로 정의한다4). 또한, 세계 측지 시스템 84(World Geodetic System 1984, WGS84)는 지피에스(GPS)의 기준 모델로 사용되며, GRS80과 매우 유사한 수치를 가지나 물리 상수의 미세한 차이로 인해 역편평률 값이 $1/f = 298.257223563$으로 정의된다.
이러한 수치 모델들은 단순히 기하학적 형상만을 규정하는 것이 아니라, 지구 질량, 자전 각속도, 중력 포텐셜 계수 등 지구의 물리적 특성을 종합적으로 반영한다. 따라서 타원체 모델은 지도 제작뿐만 아니라 지구 물리 데이터의 보정과 인공위성 궤도 결정 등 정밀한 과학 기술 분야의 기초가 된다.
적도 반지름인 장반경과 극 반지름인 단반경의 정의 및 편평률의 계산법을 다룬다.
타원체의 기하학적 특성을 수량화하기 위해 가장 먼저 정의되어야 할 지표는 이심률(eccentricity)이다. 이심률은 타원이 원에서 얼마나 벗어나 있는지를 나타내는 척도로서, 지구 타원체의 형상을 결정하는 핵심 매개변수이다. 측지학에서는 일반적으로 제1 이심률과 제2 이심률을 구분하여 사용한다. 타원체의 장반경을 $a$, 단반경을 $b$라고 할 때, 제1 이심률 $e$는 다음과 같이 정의된다.
$$ e = \sqrt{\frac{a^2 - b^2}{a^2}} $$
이와 더불어 단반경을 분모로 취하는 제2 이심률 $e'$은 다음과 같은 관계를 갖는다.
$$ e' = \sqrt{\frac{a^2 - b^2}{b^2}} $$
이심률은 편평률(flattening) $f = (a-b)/a$와 수학적으로 밀접한 관련이 있다. 예를 들어 $e^2 = 2f - f^2$의 관계가 성립하며, 이러한 수치들은 세계 측지 시스템 84(WGS84)나 지구 참조 시스템 80(GRS80)과 같은 표준 타원체 모델의 제원을 결정하는 기초가 된다.
지구 타원체는 구(sphere)와 달리 표면의 위치, 즉 지리 위도(geodetic latitude)에 따라 곡률이 계속해서 변화한다. 따라서 타원체상의 한 점에서의 곡률을 정의하기 위해서는 방향에 따른 곡률 반경(radius of curvature)을 개별적으로 고찰해야 한다. 측지 계산에서 가장 중요하게 다뤄지는 두 가지 곡률 반경은 자오선 곡률 반경과 동서 곡률 반경이다.
자오선 곡률 반경(meridian radius of curvature) $M$은 남북 방향인 자오선 호의 곡률을 나타내며, 위도 $\phi$의 함수로 다음과 같이 표현된다.
$$ M = \frac{a(1-e^2)}{(1-e^2 \sin^2 \phi)^{3/2}} $$
수식의 구조에서 알 수 있듯이, 자오선 곡률 반경은 적도($\phi = 0^\circ$)에서 최소값을 가지며 극지방($\phi = 90^\circ$)으로 갈수록 점차 증가한다. 이는 지구가 자전에 의한 원심력으로 인해 극 방향이 상대적으로 평평한 형태를 띠고 있음을 의미한다. 기하학적으로 평평한 면일수록 이를 근사하는 원의 반지름인 곡률 반경은 더 커지기 때문이다.
동서 곡률 반경(prime vertical radius of curvature) $N$은 자오선에 수직인 방향, 즉 동서 방향의 곡률을 나타낸다. 이는 타원체 표면의 한 점에서의 법선이 자전축과 만나는 지점까지의 거리로 정의되기도 하며, 수식은 다음과 같다.
$$ N = \frac{a}{\sqrt{1-e^2 \sin^2 \phi}} $$
동서 곡률 반경 역시 위도가 높아짐에 따라 그 값이 증가하는 특성을 보이지만, 모든 위도에서 자오선 곡률 반경보다 크거나 같은 값을 유지한다($N \ge M$). 특히 적도에서는 $N$이 장반경 $a$와 일치하게 된다. 이 두 곡률 반경은 측지 좌표를 3차원 직교 좌표계로 변환하거나, 타원체면상의 거리를 계산하는 측지선(geodesic) 산출 과정에서 필수적인 요소이다.
특정 지점에서의 평균적인 곡률 특성을 파악하기 위해 가우스 곡률 반경(Gaussian radius of curvature) $R = \sqrt{MN}$을 사용하기도 한다. 이는 타원체상의 미소 면적을 구하거나 지도 투영 시 발생하는 면적 왜곡을 분석할 때 유용하게 활용된다. 결과적으로 이심률과 위도에 따른 곡률 반경의 변화를 이해하는 것은 정밀한 측량 및 지형 공간 정보 구축의 기하학적 토대를 마련하는 과정이라 할 수 있다.
지심 위도, 지리 위도, 화법 위도 등 타원체상에서 정의되는 다양한 위도의 개념을 비교한다.
특정 지역이나 전 지구를 대표하기 위해 설정된 표준 타원체 모델들을 분류하고 특징을 분석한다.
특정 국가나 지역의 지형에 최적화하여 설정된 과거의 타원체 모델들과 그 한계를 기술한다.
인공위성 관측을 바탕으로 지구 전체의 질량 중심을 원점으로 하는 현대적 타원체 모델을 다룬다.
현재 가장 널리 사용되는 세계 표준 타원체 모델의 제원과 위성 항법에서의 역할을 설명한다.
현대 측지학의 기준이 되는 물리적, 기하학적 상수를 포함한 타원체 모델을 소개한다.
지구의 물리적 형상을 이해하기 위해서는 수학적으로 정의된 지구 타원체(Earth Ellipsoid)와 물리적 실체인 지오이드(Geoid) 사이의 상관관계를 파악하는 것이 필수적이다. 지구 타원체는 지구의 기하학적 형태를 단순화하여 좌표계의 기준으로 삼기 위해 도입된 회전 타원체인 반면, 지오이드는 지구의 중력(Gravity) 분포에 의해 결정되는 등포텐셜면(Equipotential surface) 중 평균 해수면과 일치하는 면을 의미한다. 지구 내부의 밀도 불균질성(Density inhomogeneity)과 지형의 기복으로 인해 실제 중력 방향은 타원체의 법선 방향과 일치하지 않으며, 이로 인해 두 면 사이에는 필연적인 차이가 발생한다. 이러한 두 면의 수직적 거리를 지오이드고(Geoid height) 또는 지오이드 기복(Geoid undulation)이라 하며, 이는 지구 물리적 특성을 해석하고 정밀한 위치 결정을 수행하는 데 있어 핵심적인 지표가 된다.
지오이드와 타원체의 관계는 고도 체계의 변환을 통해 구체화된다. 현대 측지학에서 주로 사용되는 세 가지 고도 개념은 타원체고(Ellipsoid height), 표고(Orthometric height), 그리고 지오이드고이다. 타원체고($h$)는 위성 항법 시스템(Global Navigation Satellite System, GNSS)을 통해 직접 측정되는 값으로, 준거 타원체 면에서 지표면의 한 점까지 법선을 따라 측정한 거리이다. 반면, 우리가 흔히 해발고도라고 부르는 표고($H$)는 지오이드면으로부터 지표면까지 연직선(Plumb line)을 따라 측정한 거리를 의미한다. 이들 사이의 기하학적 관계는 다음과 같은 수식으로 표현된다.
$$h = H + N$$
여기서 $N$은 해당 지점에서의 지오이드고를 나타낸다. 지오이드가 타원체보다 위쪽에 위치하면 지오이드고는 양(+)의 값을 가지며, 반대의 경우에는 음(-)의 값을 가진다. 이러한 관계식은 GNSS를 이용하여 얻은 기하학적 위치 정보를 실제 물이 흐르는 방향이나 지형적 높이와 연관된 물리적 고도로 변환할 때 결정적인 역할을 수행한다.
지오이드와 타원체의 편차는 지구 내부의 질량 분포 상태를 반영한다. 중력 이상(Gravity anomaly)이 발생하는 지역, 즉 주변보다 밀도가 높은 물질이 매장되어 있거나 거대한 산맥이 존재하는 지역에서는 중력이 강하게 작용하여 지오이드면이 타원체 밖으로 부풀어 오르게 된다. 반대로 질량이 결손된 해구(海溝)나 저밀도 지역에서는 지오이드면이 타원체 안쪽으로 함몰되는 경향을 보인다. 이러한 공간적 변동성은 지구 내부의 맨틀 대류(Mantle convection)나 지각 평형(Isostasy) 상태를 연구하는 데 중요한 정보를 제공한다. 따라서 정밀한 지오이드 모델의 구축은 단순히 고도 변환을 위한 도구를 넘어, 지구물리학(Geophysics)적 탐사 및 지구 구조 해석의 기초 자료로 활용된다.
최근에는 인공위성을 이용한 위성 중력 탐사(Satellite gravimetry) 기술의 발달로 전 지구적 규모의 정밀한 지오이드 모델이 산출되고 있다. 지구 참조 시스템 80(Geodetic Reference System 1980, GRS80)이나 세계 측지 시스템 84(World Geodetic System 1984, WGS84)와 같은 현대적 준거 타원체는 이러한 지오이드와의 편차를 최소화하도록 설계되었으나, 국지적인 질량 불균형에 따른 수십 미터 단위의 지오이드 기복은 여전히 존재한다. 따라서 정밀한 지도 제작이나 토목 공학(Civil engineering) 설계, 수로 조사 등에서는 단순한 타원체 모델에 의존하기보다, 해당 지역의 정밀 지오이드 모델을 결합하여 수준 측량(Leveling) 결과와 부합하는 고도 체계를 확립하는 것이 필수적이다. 이는 기하학적 형상과 물리적 역학계 사이의 간극을 보정하여 실용적인 공간 정보의 신뢰성을 확보하는 과정이라 할 수 있다.
두 면 사이의 거리 차이인 지오이드고의 개념과 그 물리적 발생 원인을 설명한다.
타원체고, 해발고도, 지오이드고의 정의와 이들 사이의 수치적 변환 관계를 다룬다.
현대 측지학에서 지구 타원체는 단순히 지구의 형상을 근사하는 모델을 넘어, 정밀 위치 결정과 지구 관측 시스템의 물리적 토대를 제공한다. 특히 위성 항법 시스템(Global Navigation Satellite System, GNSS)의 보급은 지구 타원체 모델의 실용적 중요성을 극대화하였다. 미국의 GPS가 채택하고 있는 세계 측지 시스템 84(World Geodetic System 1984, WGS 84)는 전 지구적 표준 타원체로서, 위성 궤도 계산과 수신기 위치 결정의 기준 프레임을 제공한다. 위성 항법 수신기는 위성으로부터 전송된 신호를 바탕으로 지구 중심을 원점으로 하는 3차원 데카르트 좌표계(Cartesian Coordinate System) 상의 위치 $(x, y, z)$를 산출하며, 이를 사용자가 이해하기 쉬운 경위도 좌표로 변환하기 위해 수학적 타원체 모델을 필수적으로 활용한다5).
지구 타원체상의 기하학적 좌표를 3차원 직교 좌표로 변환하는 관계식은 다음과 같이 정의된다. 임의의 지점에서 지리 위도(geodetic latitude)를 $\phi$, 경도를 $\lambda$, 타원체로부터의 높이를 $h$라 할 때, 직교 좌표 $(x, y, z)$는 아래의 수식으로 표현된다.
$$x = (N + h) \cos \phi \cos \lambda$$ $$y = (N + h) \cos \phi \sin \lambda$$ $$z = \left( N(1 - e^2) + h \right) \sin \phi$$
여기서 $N$은 해당 위도에서의 곡률 반경이며, $e$는 타원체의 이심률(eccentricity)이다. 이러한 수치 모델은 자율주행 자동차, 무인 항공기, 정밀 농업 등 고정밀 위치 정보가 요구되는 현대 기술 분야에서 위치 오차를 최소화하는 핵심적인 역할을 수행한다.
지리 정보 시스템(Geographic Information System, GIS)과 지도 제작 분야에서도 지구 타원체는 표준화된 공간 데이터 구축의 근간이 된다. 과거에는 각 국가가 자국의 지형에 최적화된 국지 준거 타원체를 사용하였으나, 현대에는 전 지구적 호환성을 위해 국제 지구 참조 시스템(International Terrestrial Reference System, ITRS)에 기반한 타원체 모델로 통합되는 추세이다. 이를 통해 서로 다른 국가에서 제작된 수치 지도가 하나의 좌표계 위에서 정합될 수 있으며, 이는 구글 어스나 네이버 지도와 같은 범지구적 위치 기반 서비스의 기술적 기초가 되었다. 특히 타원체는 평면 지도를 제작하기 위한 지도 투영법(Map Projection)의 출발점으로서, 투영 과정에서 발생하는 면적, 거리, 방향의 왜곡을 수학적으로 보정하는 기준면이 된다.
지구 물리 연구와 환경 모니터링 영역에서 지구 타원체는 지구의 미세한 변화를 감지하는 척도로 기능한다. 인공위성 레이저 측거(Satellite Laser Ranging, SLR)나 심우주 망원경(Very Long Baseline Interferometry, VLBI) 기술을 통해 결정된 정밀한 타원체 모델은 지각 변동이나 대륙 이동의 속도를 밀리미터 단위로 측정하는 기준이 된다. 또한 해면 고도계(Satellite Altimetry) 위성은 타원체를 기준으로 해수면의 높이 변화를 관측함으로써 기후 변화에 따른 해수면 상승 현상을 정량적으로 분석한다. 이때 실제 중력 방향을 반영하는 지오이드(Geoid)와 수학적 타원체 사이의 편차인 지오이드고(Geoid Height)를 정밀하게 산출하는 과정은 지구 내부의 밀도 분포와 역학적 구조를 이해하는 데 필수적인 정보를 제공한다.
국제 민간 항공 및 해양 항행 분야에서도 지구 타원체는 안전 운항을 위한 국제 표준으로 작용한다. 국제 민간 항공 기구(International Civil Aviation Organization, ICAO)와 국제 수로 기구(International Hydrographic Office, IHO)는 전 세계 항공로와 해도 제작의 기준 타원체로 WGS 84를 규정하고 있다. 이는 항공기와 선박이 국경을 넘나들 때 동일한 좌표 체계를 사용함으로써 위치 착오로 인한 사고를 방지하기 위함이다. 이처럼 지구 타원체는 순수 과학적 모형을 넘어, 현대 사회의 초연결 네트워크와 정밀 산업을 지탱하는 비보이지 않는 인프라로서 기능하고 있다.
지피에스 등의 시스템이 수신기의 3차원 좌표를 타원체 기반으로 산출하는 원리를 설명한다.
타원체상의 곡면 위치를 평면 지도로 옮기는 과정에서 발생하는 왜곡과 이를 보정하는 투영법을 다룬다.
중력 측정이나 해수면 변동 분석 시 타원체 모델을 기준으로 데이터를 보정하는 과정을 기술한다.