大地测量与大地水准面
Geodesy and the Geoid
你站在海边说“这里海拔 0 米”,这句话听起来简单,实际上背后是一整套地球科学工程。海面会被潮汐、风暴、洋流、气压和温度扰动;陆地会因板块运动、地下水抽取、冰盖消融后的地壳回弹而升降;GPS 给出的高度也不是“离海面多高”,而是相对于一个数学椭球面的高度。
大地测量(geodesy) 正是研究地球形状、重力场和地球自转的学科。它回答三个基础问题:我在哪里?我有多高?地球本身正在怎样变形?把地球先当成球、再用相似三角形估周长,那一步见 埃拉托色尼测地球周长;本篇从椭球、大地水准面与卫星高程接着往下写。
破除误解:GPS 高度不是海拔
手机或测量接收机能用卫星定位告诉你一个高度,但这个高度通常是椭球高(ellipsoidal height):点到参考椭球面的距离。传统地图、工程和洪水风险图更关心的是正高(orthometric height):点到大地水准面的距离,也就是接近“平均海平面”的高度。两者之间差一个量:大地水准面起伏(geoid undulation)。可以用一个近似关系理解:
正高 H = 椭球高 h - 大地水准面起伏 N如果没有大地水准面模型,GPS 的“高度”可能和工程海拔差几十米。NOAA 的 GEOID18 模型就是为了把北美地区 GPS 得到的椭球高转换成 NAVD 88 正高,服务地图、测绘与工程应用。
手机屏幕上的“海拔”几乎总是把椭球高直接显示出来。接收机解出的是点到参考椭球的法向距离;WGS84 把长半径规定为 6378137 米,扁率按约定取值,好让卫星轨道和地面坐标落在同一套几何里。这块椭球光滑、对称,并不迁就真实海面,因此也不是水位零点的别名。
大地水准面起伏把两套高度拉开。全球重力模型 EGM2008 给出的 N,大约从印度洋中部低于椭球一百零七米,到冰岛与新几内亚附近高出约八十六米。差几十米是常态,不是测量噪声。要把 GNSS 高度接到洪水图、潮位零点或旧水准点,必须先减去当地的 N;少了这一步,数字再漂亮也只是几何量。
GNSS 给出的椭球高可以在同一全球框架里互比;正高却必须声明用的是哪一张大地水准面。同一套椭球高,换成另一张 N,正高就会整体平移。海拔数字能不能并排对比,取决于 N 是不是同一张图。
大地水准面是什么
大地水准面不是一个真实可见的海面,而是一个等重力位面:如果海水只受地球重力和自转影响、没有风浪洋流,它会贴合这个面。它不是光滑椭球。地球内部密度不均,山脉、海沟、地幔异常、冰盖质量都会让重力场起伏。质量多的地方,大地水准面会鼓起;质量少的地方,它会下凹。
这就是大地水准面最重要的洞察:高度不是纯几何问题,也是重力问题。测量“高”其实是在问:水会往哪里流?同一水平面上的两点,重力势能是否相同?
等势面的意思是:把一个单位质量沿这个面挪动,重力做功为零。海水若只受重力和地球自转离心力、不受风和洋流推动,就会贴着这个面停住。
真实海面还会被风应力、密度差和西边界流抬高或压低大约一两米,这叫动力海面地形。验潮站看到的“平均海平面”因此是一个地方性的、含动力扰动的水位统计,不能直接当成全球统一的重力等位面。
潮汐再叠一层周期性起伏,瞬时海面更不能当零点。万有引力给出点质量之间的力;大地水准面是把整颗非均匀地球的质量与自转积分成一张等势面。光滑椭球故意丢掉这些起伏,所以它能算轨道,却不能告诉你水往哪边流。
三套表面:地形、椭球与大地水准面
大地测量常同时处理三套表面:
- 真实地形面:山、海、城市、冰川所在的粗糙表面。
- 参考椭球面:用数学椭球近似地球整体形状,便于坐标计算。
- 大地水准面:由地球重力场定义,近似平均海平面,是高程系统的物理基准。
地图投影 解决的是球面/椭球面到平面地图的变换;大地测量则进一步追问,这个椭球到底怎样放在地球上,高度又以哪个物理面为零点。
如果说 GIS 让不同图层叠在一起,大地测量就是保证这些图层真的“对齐”的底座。
同一地点可以同时拥有三个“高度”:离地形面的相对高差、到椭球的 h、到大地水准面的 H。工程图上的等高线通常跟 H 走,卫星轨道和 GNSS 跟 h 走,登山者看到的山脊是地形面。把三者当成同一个数,坐标会在对接处突然跳开。
参考椭球还有另一层选择:它怎么嵌入地球。原点是否与地心重合、轴向是否跟着自转轴,决定了同一组经纬度在不同框架里会不会差出好几米。WGS84 与国际地球参考框架在毫米到厘米量级上接近,旧的地方独立坐标却可能差几十到上百米。先问你站在哪一套表面上,再谈投影和叠图。
为什么它重要
大地测量不是测绘行业的内部细节,它决定许多公共问题能否算准。
海平面上升:沿海城市关心的不是全球平均海平面,而是相对海平面:海面上升多少,陆地同时是否下沉。三角洲地区的地面沉降可能让洪水风险远高于全球平均数。
桥梁、隧道与大坝:长距离工程必须统一坐标与高程基准。不同施工段如果使用不同 datum,最后可能不是“误差几厘米”,而是结构无法闭合。
地震与火山监测:GNSS 连续站可以测到板块一年几毫米到几厘米的运动。地震前后的同震位移、火山膨胀、冰川卸载后的地壳回弹,都需要大地测量。
地下水与冰盖质量变化:GRACE 与 GRACE-FO 重力卫星通过测量地球重力场变化,反推陆地水储量、冰盖质量和海洋质量重新分布。
相对海平面把这套框架接到气候风险上。全球平均海面可以在涨,某段海岸却因为陆地抬升而相对下降;抽取地下水造成的沉降,则能让相对海平面涨得比全球平均数快一个数量级。海平面变化真正决定一座城市会不会被淹的,是水位相对当地陆地,而不是全球一个平均数。
从三角测量到卫星大地测量
传统大地测量依赖三角网、水准测量和天文观测。测量员在山顶架设经纬仪,一段一段把控制网铺开。二十世纪后半叶,卫星改变了尺度。VLBI、SLR、DORIS、GNSS、卫星重力任务共同构成现代大地测量系统。今天,一个全球参考框架不再是画在纸上的网,而是由遍布全球的观测站、卫星轨道、地球自转参数和重力模型共同维持。
这也是为什么 USGS 要求出版物明确记录水平与垂直 datum。没有 datum 的坐标像没有单位的数字:看似精确,实际不可复现。
今天的全球框架是运动学的。国际地球参考框架由 GNSS、VLBI、SLR、DORIS 共同维持,每个站都带速度,坐标必须写成某框架、某历元。少了框架名和日期的数字无法复现,因为板块以厘米每年挪动,任何静止坐标都只是某一天的快照。
一个工程现场会怎样用它
想象一座沿海城市要重建排水系统。工程师从 GNSS 接收机得到井盖、泵站和河堤的椭球高;水务部门手里有旧水准网的正高;洪水模型又需要统一的地形高程。如果直接把这些数字混在一起,模型会把一些本来低洼的街区看成安全区,或把安全区误判成淹没区。正确做法是先确认坐标参考框架和高程基准,再用地区大地水准面模型把椭球高转换为正高,并检查本地是否存在地面沉降。三角洲、填海区、矿区和地下水超采区尤其敏感,因为地表本身可能每年下降数毫米到数厘米。
这说明大地测量不是“地图上差一点没关系”的技术。它直接影响洪水边界、保险费率、道路坡度、污水自流方向和防潮堤高度。
混用这些数字之所以危险,不是因为少做了某道放样工序。椭球高、正高、旧潮位零点看起来都叫“米”,零点却不是同一个物理面,相减没有意义。先问高度相对于谁,比先问仪器显示了几位小数更重要。
时间维度:坐标也有日期
现代大地测量越来越强调坐标的历元(epoch)。同一 GNSS 站在 2010 年、2020 年、2030 年的位置并不完全一样。板块运动是长期背景:北美、西欧、澳大利亚、太平洋板块都在以每年厘米级速度相对移动。地震则会造成突然跳变,强震后的余滑还会持续数月到数年。冰后回弹也是典型例子。末次冰期消退后,斯堪的纳维亚和加拿大部分地区地壳仍在缓慢抬升;同时,冰盖消融让重力场和海平面相对关系发生改变。沿海海平面记录若不校正陆地升降,就会把地面运动误读成海水运动。
因此,精密坐标不只是“经纬度 + 高程”,还应包含参考框架、观测日期和速度模型。地球不是静止的测量对象,而是测量系统的一部分。
连大地水准面也有日期。冰盖消融和质量再分布会改重力场,N 会以毫米每年的量级缓慢漂移。精密高程若要跨十年对比,需要说明用的是哪一代重力模型和哪一个历元的框架。高度不是一次性量完就永远有效的标签。
板块构造在地质图上画的箭头,和 GNSS 连续站给出的速度场是同一件事,只是时间单位从百万年换成了年。框架带着速度,正是为了把这块缓慢变形的底板写进坐标。
重力异常:看不见的地下结构
重力测量还能揭示地下结构。密度较高的岩体、矿体或地幔上涌会造成正重力异常;低密度沉积盆地、盐丘或部分岩浆房会造成负异常。这不是“隔空看见地下”,而是用牛顿引力反推质量分布。单靠重力资料通常存在多解性,必须与地震波、钻孔、磁法和地质图共同约束。但在海洋、冰盖和难以进入地区,重力异常常是理解地壳结构的少数窗口。
大地水准面模型本身也依赖全球重力观测。也就是说,我们脚下“高度零点”的定义,最终连接到整颗地球的质量分布。
EGM2008 一类全球模型把卫星轨道摄动、卫星跟踪和地面重力异常拼成球谐展开。大地水准面是这张场里一张选定的等位面,长波由卫星约束,短波仍要靠航空或地面测量补。没有这张场,GPS 高度就变不成正高;没有地面短波,区域转换会在山脉和海沟处留下米级偏差。
争议与边界
大地测量的难点不在于“没有数据”,而在于所有东西都在动。板块每年移动,地震会瞬间改变地表位置,冰盖消融会让地壳回弹,地下水抽取会造成沉降。一个坐标参考框架必须不断更新,否则“同一个点”的坐标会随时间漂移。高程基准也并非全球统一。美国的 NAVD 88、欧洲各国高程系统、地方潮位站基准之间不能随意混用。全球工程、卫星遥感与海平面研究越来越需要可互操作的垂直基准。
更深的哲学问题是:我们习惯把地球看作一个固定舞台,但大地测量告诉我们,舞台本身在缓慢弯曲、漂移、升降。所谓“位置”,其实是某个时间、某个参考框架、某套物理模型下的约定。
把各国旧水准网拉到同一张全球等势面上,是正在进行的垂直基准现代化,而不是已经完成的事实。GNSS 能给出全球一致的椭球高,正高却仍取决于当地大地水准面和验潮原点。在国际高程参考框架被广泛采用之前,跨洋对比“海拔”仍要先做基准转换,不能假设所有地图上的零点已经重合。
跨域连接
- 数值方法:用重力反推地下质量分布是病态反演——同一份重力异常可由多种密度结构产生,单靠重力资料必然多解。所以勘探必须把重力与地震波、钻孔、磁法共同约束;推论:只凭重力异常就宣布"发现矿体"的报告,可信度天然有限。
- 地图投影:投影解决"球面怎么摊平",大地测量解决"这个数学椭球怎么放上地球、高度以谁为零"。三套表面——真实地形、参考椭球、大地水准面——各司其职;把它们混为一谈,是一切坐标混乱的总根源。
- 欧几里得几何:参考椭球是把地球压进规则可算几何的近似。近似就有代价:卫星给出的椭球高与工程要的"海拔"(正高)之间,差着可达几十米的大地水准面起伏。推论:不做水准面模型转换就把卫星高度当海拔,洪水风险图会系统性地画错。
- 引力物理学:高度不是纯几何问题,而是重力问题——大地水准面是等重力位面,"水往哪里流"由它决定。因此重力卫星测量重力场的微小变化,就能反推冰盖消融与地下水储量;地球的质量分布一变,"海拔零点"本身就跟着动。
- 保险与工程经济:高程基准直接连着钱:洪水边界决定保险费率,防潮堤高度决定工程投资,长距离工程依赖统一基准。坐标没有基准就像数字没有单位——推论:基准混用的工程,误差不会在图纸上显形,而会以结构无法闭合或错保的形式兑现。
参考文献
- National Geodetic Survey. GEOID18 Interactive Computation. NOAA, 2022.(GPS 椭球高与 NAVD 88 正高转换)
- USGS. Guidance on Use and Documentation of Horizontal and Vertical Datums in USGS Publications. 2019.
- Hofmann-Wellenhof, B. & Moritz, H. Physical Geodesy. 2nd ed., Springer, 2005.
- Seeber, G. Satellite Geodesy. 2nd ed., De Gruyter, 2003.
- Tapley, B. D. et al. "Contributions of GRACE to understanding climate change." Nature Climate Change, 9, 2019. doi:10.1038/s41558-019-0456-2
- Pavlis, N. K., Holmes, S. A., Kenyon, S. C., & Factor, J. K. "The development and evaluation of the Earth Gravitational Model 2008 (EGM2008)." Journal of Geophysical Research: Solid Earth 117, 2012. doi:10.1029/2011JB008916
- Altamimi, Z. et al. "ITRF2020: an augmented reference frame refining the modeling of nonlinear station motions." Journal of Geodesy 97, 2023. doi:10.1007/s00190-023-01738-w
延伸阅读
- Torge, W. & Müller, J. Geodesy. 4th ed., De Gruyter, 2012.
- UNAVCO / GAGE Facility. Geoid Height Calculator.(大地水准面起伏计算工具)