大地水准面

EGM96模型给出的全球大地水准面高分布图，红、蓝两色分别表示高于和低于椭球面

大地水准面（德語：Geoid）是指地球重力场中，与处于自由静止状态的海平面相重合或最为接近的重力等位面。大地水准面的概念最早由德国大地测量学家卡爾·弗里德里希·高斯在1828年提出，当时高斯以“地球的数学表面”来指称与重力方向相垂直、且包含了静止的平均海水面的几何表面，并提出将前述的重力等位面作为大地高程系统的基准面。其后于1873年，高斯的学生利斯廷创造出了“Geoid”一词，用以描述高斯所提出的数学表面。

在大地测量学中，大地水准面被视为是地球的物理形状和数学形状。由于自然的地形表面形态过于复杂，大地测量学通常是将重力场中整体形状与自然表面最为接近的等位面作为地球的形状进行研究。处于静力平衡状态下的平均海水面被视作是符合这一标准的重力等位面，这一假想的海水面不受潮汐、风浪及大气压变化影响，仅在地球引力和因地球自转产生的离心惯性力的作用下保持平衡。将该平均海水面所处的重力等位面延伸到陆地内部，形成的闭合曲面即为大地水准面，其所包围的形体又被称为大地体。

1849年，英国物理学家斯托克斯提出了计算大地水准面的斯托克斯方法，即将重力位分为正常重力位和扰动位两部分，再将一个理想的、与大地水准面最为密合的旋转椭球面作为正常重力位的等位面，并以大地水准面高描述大地水准面相对于该参考椭球面的起伏，通过扰动位计算出大地水准面高，从而确定出大地水准面的形状。传统的测定大地水准面的方式即是通过在地面上建立天文大地网的方式进行重力测量，再以斯托克斯方法进行计算。这种方式受到海洋和境界的分隔，往往只能求得一定区域内的地球形状。随着20世纪中后期卫星重力测量等技术逐渐发展成熟，全球范围的大地水准面形状得以被测定至分米级的精度。莫洛坚斯基、布耶哈马等大地测量学家，亦在斯托克斯方法的基础上提出了确定大地水准面的新方法。

大地水准面是测量外业所依据的基准面，在测量学中具有重要地位。在各类高程系统中，日常中称为海拔高程的正高系统是基于大地水准面建立的。如何确定大地水准面的形状，是物理大地测量学所研究的关键问题之一。

定义

经典定义

1873年，利斯廷提出了大地水准面的经典定义，即大地水准面是与平均海水面相重合的重力等位面。在1975年，美国大地测量学家理查德·拉普（英語：Richard H. Rapp）进一步完善了这一定义，他将重力位 $W$ 视为除大气层外的地球质量所引起的重力位 $W_{g}$、大气质量所引起的重力位 $W_{a}$，以及受日、月的引潮力作用而引起的重力位 $W_{t}$的和：

$W=W_{g}+W_{a}+W_{t}$

不同于 $W_{g}$，后两项与地球内部的质量分布无关，因此可以被相对精确地计算。拉普认为，定义大地水准面的重力位应当仅包含 $W_{g}$和 $W_{a}$，即大地水准面应当是使 $W_{g}+W_{a}$ 等于常数的重力等位面，而因引潮力、风和洋流等其他因素影响产生的变化则不在考虑范围之内。这种定义方式认为，通过多年的验潮数据可以得到一个平均海水面，其上的任意一点的位置与海洋面的长期平均位置相重合。

然而，20世纪60年代以来对平均海水面的研究发现，其与大地水准面间存在最大可达2米的系统性的偏移，这就使大地水准面的经典定义暴露出了缺陷。因此，马瑟（英語：Mather R. S.）等人在1978年提出，大地水准面应重新定义为与平均海水面最为密合的重力等位面。

相对论定义

1985年，瑞典大地测量学家布耶哈马基于相对论重力提出了相对论大地水准面的概念，即一个使其上所有精密时钟的运行速率相同且与平均海水面最为密合的封闭曲面。在2019年，德国不来梅大学的丹尼斯·菲利普（德語：Dennis Philipp）利用新的测量结果，以及无含时红移势能（英語：time-independent redshift potential）等新的理论，对相对论大地水准面的概念给出了数学上的表述方法。由这一方法给出的相对论水准面与经典大地水准面之间存在着毫米级的差异。

形状

表现大地水准面差距的假彩色图像，差距放大了10,000倍，地形添加了阴影效果

表现大地水准面差距的假彩色图像，差距按比例表示

从数学上看，大地水准面是一个连续但不规则的闭合曲面，它与经过这一曲面的铅垂线处处正交。地球内部质量分布，特别是外层质量分布的不均匀性，使得大地水准面的形状变得特别复杂。在发生密度突变的区域，大地水准面的曲率会表现出不连续性。因此，大地水准面并非是解析曲面。

大地水准面的形状可以用椭球面近似，该椭球面也被称为参考椭球面。以WGS84坐标系中的椭球为例，其长半轴位于赤道面上，长度约6,378,137米；短半轴指向极点，长度约6,356,752米，比长半轴短约21千米。相对于WGS84椭球，大地水准面的起伏从最高点到最低点的范围仅不到200米。在这一差距被放大的情况下，大地水准面的形状则与土豆或梨相似。尽管大地水准面是一个不规则的曲面，但它仍比真实的地形表面更加平滑。相对于大地水准面，真实地形表面的高程差跨越了近20千米：其最高点出现在珠穆朗玛峰，高度为+8,848 米；其最低点出现在马里亚纳海沟，高度为−11,034米。

从赤道观察大地水准面（红线）相对于参考椭球面（黑线）的形状，两者间的差值为同一纬线上大地水准面差距的平均值

从地极处观察大地水准面（红线）相对于参考椭球面（黑线）的形状，两者间的差值为同一经线上大地水准面差距的平均值

起伏

主条目：大地水准面高

大地水准面起伏又称大地水准面高或大地水准面差距，是指大地水准面上的某一点，与其沿法线投影到参考椭球面上的相应位置之间的距离。由于参考椭球面是人为定义的规则曲面，大地水准面的形状可以通过大地水准面高进行表达。相对于WGS84椭球，大地水准面的正极值出现在新几内亚，为+85.4米；负极值出现在印度南部，为−107.0米。

大地水准面的起伏受地球内部密度分布的直接影响。通过频谱分析，可以将大地水准面的起伏分为长波（尺度为上千公里，幅度为±100米左右）、中波（尺度为上百公里，幅度为±10米）和短波（尺度为数十公里，幅度为±1米）部分。其中，长波部分反映地幔及地球内部的密度异常分布，而中短波部分则反映了岩石圈和地形的影响。在地球物理学研究中，大地水准面的起伏具体可分为地形的不规则部分、结晶基底起伏、莫霍面起伏、岩石圈底面起伏，以及因地幔和地核部分的横向密度不均匀造成的起伏等。

EGM96地球重力场模型给出的全球大地水准面高分布图

正异常

大地水准面高在北大西洋、欧洲、非洲的西部和南部、印度洋西南部、安地斯山脈、西太平洋及太平洋的赤道部分等地表现出正值，这些地区几乎包含了世界上所有的大型地盾、显生宙地台和地质熱點。其中，规模庞大而幅度较高的有两个区域，横跨了大西洋和非洲的正异常区域被称为大西洋-非洲大地水准面高值区（英語：Atlantic-African geoid high），太平洋中西部的正异常区域则被称为赤道太平洋大地水准面高值区（英語：Equatorial Pacific geoid high）。西太平洋的正异常区域在中生代早期曾是俯冲杂岩帶所处的位置。其他与俯冲带不相关的正异常区域则可能曾是大陆板块的聚集处，或是古板块的中心等会发生地幔绝热现象的地区，地幔的绝热会使得其本身的温度升高，其产生的抬升和岩浆作用能够产生大地水准面的正异常。

负异常

大地水准面高在印度及印度洋北部，以及南极洲附近表现出负值。其中又以印度洋北部的区域表现得尤其显著，该区域也被称为印度洋大地水准面低值区（英語：Indian Ocean geoid low）。对于这些负异常形成的原因，存在有众多推测，如软流层中的浅层物质流、板块俯冲在中地幔引起的地幔对流，以及地幔和地核的边界处发生的凹陷等。但这些推测都同意，印度洋大地水准面低地产生的原因有部分是特提斯洋在中生代的俯冲带来的密度异常。

测定

在人造卫星出现以前，测定大地水准面采用的首要方法是斯托克斯方法，这一方法的依据是乔治·斯托克斯在1849年提出的斯托克斯定理：在大地水准面的外部没有质量分布的情况下，其形状和其外部的重力场可以由地球的质量、自转角速度以及大地水准面上的重力观测值唯一确定，而与其内部的质量分布无关。这一定理说明，大地水准面的形状可以通过地面上的重力观测数据确定。然而，其在应用过程中存在诸多困难：公式中的积分难以收敛、全球性的观测数据难以获得等。

现代测定大地水准面的方式主要是结合卫星和地面的重力观测数据，采用球谐函数展开的方式对地球重力场进行表达。其中，球谐函数中的高阶项主要通过地面的重力观测数据获得，而低阶项主要是通过卫星重力探测的方式测得。卫星重力探测技术至今发展为了地面跟踪卫星轨道摄动、卫星测高、卫星跟踪卫星以及卫星重力梯度测量四种模式。而诸如CHAMP、GOCE和GRACE等新一代重力测量卫星的出现，不仅使对地球重力场的研究得以以空前的精度和空间分辨率进行，还使得对地球重力场的时变监测成为可能。以GRACE卫星为例，其对长波时变重力场的时间分辨率可达10天，测定的大地水准面年变化率精度为0.01毫米每年。

斯托克斯方法

利用斯托克斯方法确定大地水准面的形状，包括两部分内容：确定大地水准面高相对于参考椭球面的高程差（即大地水准面高），以及这一高程差在参考椭球面上的曲面位置。这两部分分别可以用斯托克斯公式和维宁·曼尼兹公式进行计算。然而，实际上的重力观测数据是在地形表面，而非大地水准面上得到的。因此，在得到地表的观测数据后，还需要除去大地水准面外的质量对重力观测值的影响，并将地面上重力观测值的位置改正到大地水准面上，这一部分的内容又被称为重力归算。通过斯托克斯方法，可以求得较高精度的高程异常并随之确定相应的大地水准面高。

斯托克斯公式

斯托克斯公式给出了大地水准面高与全球大地水准面上的重力异常及其分布的数学关系：

$N={R \over 4\pi G}\iint \limits _{\sigma }\Delta g\,S(\psi )\operatorname {d} \!\sigma$

上式中，$N$ 是大地水准面高，$R$ 是地球的平均半径，$G$ 是万有引力常数。而对于积分部分，${\begin{smallmatrix}\iint \limits _{\sigma }\end{smallmatrix}}$表示在整个单位球面上的积分，$\operatorname {d} \!\sigma$ 是立体角元素。 $\Delta g$ 是单位角元素上的重力异常，由重力测量数据给出的重力位减去计算得到的正常重力位后得到。$S(\psi )$ 则为斯托克斯函数，该项是单位球面上的被计算点与重力异常观测值所在的角元素之间的夹角 $\psi$ 的函数。

维宁·曼尼兹公式

在由斯托克斯公式求得大地水准面相对于参考椭球面的高程差后，还需通过垂线偏差改正，将大地水准面上的天文坐标转换成参考椭球面上的地理坐标。这项改正能够通过通过维宁·曼尼兹公式计算：

${\begin{Bmatrix}\xi \\\eta \end{Bmatrix}}={1 \over 4\pi G}\iint \limits _{\sigma }\Delta g\,{\begin{Bmatrix}\cos {\alpha }\\\sin {\alpha }\end{Bmatrix}}\,S(\psi )\operatorname {d} \!\sigma$

上式中的 $\xi$ 表示天文坐标相对于地理坐标的纬度差， $\eta$ 表示天文坐标相对于地理坐标的经度差，$\alpha$ 则是天文坐标相对于地理坐标的方位角。

拉普拉斯方法

地球的扁率及其他形状和重力的不规则性，会对人造卫星在空间中的运动造成影响，即卫星轨道的摄动。通过对卫星轨道进行动力学分析，可以反演出地球的引力场模型。同时利用引力位在质体外部满足拉普拉斯方程的特性，可以将其展开成球谐级数进行逼近。

拉普拉斯方程的球谐级数解

通过求解球面坐标中的拉普拉斯方程 $\Delta V=0$ ，单位球外部（$r>1$）的引力位可以展开为：

$V={\sum _{n=0}^{\infty }}{\sum _{m=0}^{n}}\left[A_{nm}{R_{nm}(\theta ,\lambda ) \over r^{n+1}}+B_{nm}{S_{nm}(\theta ,\lambda ) \over r^{n+1}}\right]$

上式中各个量的含义如下：

• $(r,\theta ,\lambda )$是空间中某特定点的球坐标，$r\$是点的向径， $\theta$和 $\lambda \$分别是点的极角和方位角

• $A_{nm}$ 和 $B_{nm}$ 是通过对单位球上的观测数据积分确定的待定系数

• $R_{nm}(\theta ,\lambda )=P_{nm}(\cos \theta )\cos {m\lambda }$，$S_{nm}(\theta ,\lambda )=P_{nm}(\cos \theta )\sin {m\lambda }$，是与缔合勒让德多项式有关的面谐函数

上式中的向径 $r$ 是被计算点相对于单位球的相对半径，而非其与地心的真实距离。$R_{nm}(\theta ,\lambda )$ 和 $S_{nm}(\theta ,\lambda )$ 是与勒让德多项式 $P_{nm}(\cos {\theta })$

卫星动力学中的球谐级数

在卫星动力学中，球谐级数常写成如下形式：

$V={GM \over r}\{{1-\sum _{n=1}^{\infty }}{\sum _{m=0}^{n}}{({a \over r})}^{n}\left[J_{nm}R_{nm}(\theta ,\lambda )+K_{nm}S_{nm}(\theta ,\lambda )\right]\}$

上式中各个量的含义如下：

• $(r,\theta ,\lambda )$是空间中某特定点的球坐标，$r\$是空间中点的地心距离， $\theta$和 $\lambda \$分别是点的地心纬度的余角和经度

• $GM$ 为重力场模型的地心引力常数
• $J_{nm}$和 $K_{nm}$ 是与地球动力学有关的参数，其与 $A_{nm}$ 和 $B_{nm}$ 的关系为

${\begin{cases}A_{nm}=-GMa^{n}J_{nm}\\B_{nm}=-GMa^{n}K_{nm}\end{cases}}$

在这一表达式中，参数 $A_{nm}$、$B_{nm}$、$J_{nm}$和 $K_{nm}$均由地球本身的特性决定，包含了地球的转动惯量和动力学形状等信息，是引力场模型中的待解参数。

模型

球谐函数的示意图，地球的引力场被视作是以上级数项的叠加

重力场模型

大地水准面是通过重力场模型进行表达的，对全球的重力场进行描述的模型即为全球重力场模型。重力场模型是对真实的重力场进行近似的数学工具，利用重力场模型可以计算得到重力场中的相关物理量。在建立重力场模型的过程中，重力位 $W$ 被分成引力位 $V$ 和离心力位 $\Phi$ 两部分（即 $W=V+\Phi$）， 其中离心力位 $\Phi$ 能够根据地球的自转角速度 $\omega$ 和点在空间中的坐标 $(x,y,z)$ 直接计算得到：

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

然而，引力位 $V$ 实际上是由边界内部的质量分布所决定的，这在确定大地水准面的过程中始终是未知量。因此，需要通过在大地水准面上的观测值来求得满足一定条件的引力位函数，这样的问题也即大地测量边值问题。通常，这一问题的解法是利用引力位在质体边界面的外部满足拉普拉斯方程的性质，在边界面的外部将其展开成用球谐级数进行逼近。值得注意的是，这一过程是在假设大地水准面的外部没有质量分布的情况下进行的。如何计算大地水准面外部的质量分布对地球重力位的影响并予以去除的过程，是重力归算所研究的问题。

引力场位的确定比离心力位的确定要复杂得多，引力场模型的建立是重力场模型建立过程中的主要部分。因此，在不严谨的表述中，有时也将引力场模型直接称为重力场模型。

球谐表达式

在引力场模型中，引力位 $V$ 在地球外部（$r>a$）表达为如下形式：

$V={\frac {GM}{r}}{\sum _{n=0}^{n_{\text{max}}}}\left({\frac {a}{r}}\right)^{n}{\sum _{m=0}^{n}}{\overline {P}}_{nm}(\sin \phi )\left[{\overline {C}}_{nm}\cos m\lambda +{\overline {S}}_{nm}\sin m\lambda \right]$

上式中的各个量的含义如下：

• $(r,\phi ,\lambda )$是空间中某特定点的球坐标，$r\$是点的地心距离， $\phi \$和 $\lambda \$分别是点的地心纬度和经度

• $GM$ 为重力场模型的地心引力常数
• $n_{\text{max}}$为该重力场模型的最大阶数

• $a$ 为地球（参考椭球体）的长半轴
• ${\overline {P}}_{nm}$ 是 $n\$阶 $m\$次的完全正规化缔合勒让德多项式
• ${\overline {C}}_{nm}$ 和 ${\overline {S}}_{nm}$ 是由测量数据所确定的该重力场模型的完全正规化系数

其中，当 $n=0$ 时，${\overline {P}}_{00}(\sin \phi )=1$，${\overline {C}}_{00}=1$，则求和符号右侧所有项的和为 $1$；当 $n=1$ 时，${\overline {C}}_{1,m}={\overline {S}}_{1,m}=0$，则求和符号右侧所有项的和为 $0$. 因此，上式可进一步表达为如下形式：

$V={\frac {GM}{r}}\left(1+{\sum _{n=2}^{n_{\text{max}}}}\left({\frac {a}{r}}\right)^{n}{\sum _{m=0}^{n}}{\overline {P}}_{nm}(\sin \phi )\left[{\overline {C}}_{nm}\cos m\lambda +{\overline {S}}_{nm}\sin m\lambda \right]\right)$

在得到引力位 $V$ 和椭球面坐标为 $(r,\phi ,\lambda )$ 的关系后，再将其与离心力位模型 $\Phi$ 相加，即可得到重力位模型 $W$，最后得到表达为 $W(r,\phi ,\lambda )=W_{0}$的大地水准面的形状。 此外，由于重力位的梯度即为重力矢量（即 $\nabla W=\mathbf {g}$）， 这一模型同时也定义出了重力加速度的大小。在实际的模型建立的过程中，也常取引力位中的数阶最大项作为正常重力场中的引力位，再将余下的引力位以扰动位的形式进行表达，正常重力场中的等位面即为参考椭球面。在美国国家地理空间情报局（NGA）发布的EGM96模型中，$n_{\text{max}}=360$，即该模型给出了最高完全至360阶次的重力场模型系数，其表达的大地水准面的分辨率可达55千米或110千米（取决于对分辨率定义的不同）。而在2009年，NGA又发布了最大阶数$n_{\text{max}}=2190$，最高完全阶次为2159次的EGM08模型，其大地水准面的分辨率可达12千米，精度在全球范围类则可达±15厘米。

空间分辨率

通常来说，根据奈奎斯特采样定理，球谐函数表达式中的阶次越高，其表达的范围也就越精细。对于以球谐函数来表达的重力场模型，其空间分辨率通常以最短半波长表示。在地球表面，空间分辨率 $\psi _{\text{min}}$ 与最大阶数 $n_{\text{max}}$ 的关系为：

$\psi _{\text{min}}={\pi R \over n_{\text{max}}}\approx {20\,000\,{\text{km}} \over n_{\text{max}}}$

其中的$R$ 为地球半径。通过这一公式，可以求出现有重力场模型常采用的空间分辨率：

最大阶数	最短半波长 / km
100	200.00
200	100.00
360	55.55
2190	9.13

现有的地球重力场模型

21世纪以来，利用CHAMP、GRACE等重力卫星的测量数据，以及海洋测高重力数据和地面重力数据等多种数据源，精度达厘米级且高时空分辨率的地球重力场模型得以被建立。这些重力场模型由国际地球重力场模型中心（ICGEM）负责收集和存档，并对其进行可视化、建模、计算以及使用DOI进行标记等。截至2020年，ICGEM存储的阶数超过2000阶的超高阶静态全球重力场模型有如下几个：

模型	发布时间	最高阶数	数据来源	参考文献
EGM2008	2009年	2190	A、G、S(Grace)
EIGEN-6C4	2014年	2190	A、G、S(Goce、Grace、Lageos)
GECO	2015年	2190	EGM2008、S(Goce)
SGG-UGM-1	2018年	2159	EGM2008、S(Goce)
XGM2019e_2159	2019年	760/2190/5540	A、G、S(GOCO06s)、T
注：其中A表示测高数据，S表示为卫星测量数据，G表示地表测量数据（包括陆地、船舶和航空测量），T表示测图数据

正高系统

椭球高h、大地水准面高N和正高H的示意图，图中红色实现为大地水准面（Geoid），黑色虚线为椭球面（Ellipsoid），注意其中的N为负值

主条目：正高

在大地测量学中，高程系统是定义某点沿特定的路径到一个参考面上距离的一维坐标系统。在传统的几何水准测量中，两点之间的高差是基于两点所在的水准面相互平行的假设进行的。当水準管居中时，水準儀的竖轴与铅垂线重合，其视线方向相切于其所在的水准面。若水准面相互平行，则在两点间进行对向观测，得到的两点间的高差应当相等；在环线上进行闭合观测，其闭合差应当为零。但事实上，水准面之间是并不相互平行，由此造成的在环线观测中留下的差值被称为水准环闭合差。传统的几何测量手段不能解决水准环闭合差不为零的问题，所以高程系统的建立必须基于重力水准测量的原理。

以大地水准面为参考面建立的高程系统被称为正高系统。在正高系统中，正高被定义为地面上任意一点沿其所在的铅垂线方向到大地水准面的距离。在重力水准测量中，通过观测大地水准面上一点 $A$ 与任意一点 $B$ 之间的重力位差${\begin{smallmatrix}\int \limits _{AB}g\operatorname {d} \!h\end{smallmatrix}}$，再通过点 $B$ 沿铅垂线到大地水准面的路径上的重力平均值 $g_{\text{m}}^{B}$，可以求得点 $B$ 的正高为：

$H^{B}={1 \over g_{\text{m}}^{B}}\int \limits _{AB}g\operatorname {d} \!h$

正高与大地高的关系

在大地测量学以外的领域，正高也被称为海拔高度。而在进行空间位置的计算（如计算GPS等全球卫星导航系统中的卫星位置）时，使用的高度则是相对于参考椭球面的距离，称为大地高或椭球高。两者间的转换需要通过大地水准面高进行：

$h=H+N$

其中，$h$ 表示大地高（椭球高），$H$ 表示正高，$N$ 表示大地水准面高。

正高与正常高的关系

主条目：正常高

正高的计算需要用到重力平均值 $g_{\text{m}}^{B}$，但该值无法在地球内部的质量和重力分布未知的情况下得到，所以正高无法被精确求得。1945年，苏联大地测量学家莫洛坚斯基提出了利用正常高来代替正高的计算方式。在这一高程系统中，由地球表面向下量取正常高而获得的表面被称为似大地水准面，由正常椭球面向上量取正常高而获得的表面则被称为似地形表面，而似大地水准面和正常椭球面之间的差距又被称为高程异常。似大地水准面与水准面不同，它没有具体的物理意义，更不是重力等位面，只是用于计算的辅助面。大地水准面和似大地水准面在山区的差异可达4米，但两者在平均海水面上是重合的。

大地水准面重力位

大地水准面重力位 $W_{0}$ 与高程基准相互对应，是人为规定的参考值。这一参考值应具有不随时间变化、长期固定，且与平均海水面的高度有明确关系的特征。在2015年，国际大地测量协会（IAG）公布的大地水准面重力位参考值为$W_{0}=6.263\,685\,34\times 10^{7}\,{\text{m}}^{2}\cdot {\text{s}}^{-2}$。但在不同的区域高程系统中，其定义的高程系统仍可能对应出不同的大地水准面重力位值。以中国1985国家高程基准为例，其与全球高基准的差异为21到23厘米。

相关条目

• 布隆斯公式

• 大地测量系统
• 大地测量边值问题
• 地球坐标系统

外部链接

• EGM96地球重力场模型（英語）
• EGM08地球重力场模型（英語）
• 国际地球重力场模型中心（ICGEM）（英語）
• 国际大地水准面服务（ISG）（英語）
• 美国国家海洋和大气管理局（NOAA）（英語）
• 椭球面、大地水准面、重力、大地测量学和地球物理学教程（Tutorial Ellipsoid, geoid, gravity, geodesy, and geophysics）（英語）
• GRACE网站上对确定地球重力方法的阐述（英語）