夜雨聆风
你有没有想过一件事。
卫星在天上往下照，它看到的真的是海面吗。
当然不是。它看到的是一层又一层被扭曲过的影像，就像你透过好几块脏玻璃在看东西。大气在折射信号，电离层在拖慢脚步，海浪在骗你的眼睛，潮汐在把海面一会儿推高一会儿压低，连气压都在偷偷改变海面的形状。
我们之前聊过那场持续了四十五年的宇宙接力赛，从 Seasat 到 SWOT，一代代卫星上天，精度从十厘米磨到亚厘米级。但有一个关键问题我一直没展开讲，那就是，这十厘米到亚厘米之间，到底发生了什么。
今天这篇文章，就是来拆这件事的。
一张最简单的账单
测高的基本逻辑其实就一个公式。
是我们最终想要的真实海面高度，
是卫星的轨道高度，
是雷达测到的距离，
是所有需要扣除的校正项之和。
看起来简单对吧。但问题是这个 
里面塞了多少东西。
我给你报个数你就知道了。干对流层延迟大约 2.3 米，湿对流层零到五十厘米，电离层最高能到 50 厘米，海洋潮汐能达到数米量级，逆气压效应正负 15 厘米，海况偏差零到六厘米高海况下甚至能到 10 厘米。
也就是说，卫星原始测量出来的那个数字，跟真实海面高度之间的差距，最大可能有好几米。而我们要从中提取出厘米级甚至亚厘米级的海面变化信号。
这就好比你要在一堆噪音里听清一个人轻轻说了一句"海平面今年涨了三毫米"。你得先把所有噪音一个个地抠掉，而且每个都扣得足够干净，剩下的那一点点才是真货。
所以这一整套校正流程，本质上就是在做一件事，把传感器「看到的那个扭曲的世界」，一层一层还原成「真实的世界」。
先把地基打牢
在聊任何一项具体校正之前，我得先说一个容易被忽略但极其致命的前提。
公式里的 
，也就是卫星轨道高度，是所有计算的基准。如果 
本身就不准，后面所有校正做得再精致都是白搭。你想啊，
的误差会一比一地传递到最终的 
里面去。
现代测高任务是怎么定轨的呢。不是靠一颗 GPS 就完事了，而是用了多重定位约束的组合拳。
DORIS 多普勒测速系统、GPS/GNSS 导航星座、SLR 激光测距，三套系统同时工作，互相交叉验证。再加上动力学模型，要把地球重力场、大气阻力、太阳光压这些扰动因素全部算进去。
以 Jason 系列为例，它的轨道精度（
）已经逼近 2 到 3 厘米。这个数字意味着什么，意味着光是把「卫星在哪」这个问题回答清楚，人类就已经花了四十年的时间在迭代了。
地基打稳了，接下来才开始拆那些滤镜。
第一道滤镜，天空的减速带和迷雾
电磁波从卫星飞向海面的时候，它不是在真空里跑的。它得穿过大气层，而大气层会让它变慢。
这道滤镜分两层。一层叫干对流层（Dry Troposphere），一层叫湿对流层（Wet Troposphere）。
干对流层是由干燥空气引起的折射效应。它是所有校正项里量级最大的一个，大概 2.3 米。但好消息是它特别稳定，因为干空气在全球范围内的变化不大。科学家用 ECMWF 这种全球数值天气预报模型提供的气压场数据就能算得很准，校正完之后残余误差可以控制在 1 厘米以内。
湿对流层就不一样了。它是由大气中的水汽引起的，量级相对小一些，0 到 50 厘米，但它的空间和时间变化极其剧烈。热带地区水汽多，延迟就大；干旱地区水汽少，延迟就小；同一个地方今天下雨明天晴天，数值能差很多。
这个项怎么校呢。靠卫星自己带的微波辐射计（MWR）。MWR 实时测量大气的亮温，反演出水汽含量，然后做校正。但即便如此，湿对流层依然是整个校正链条里残余误差最大的重灾区之一，校正后还有 1 到 2 厘米的不确定性。
我有时候会觉得这两层挺像人的性格的。干对流层量大但老实，你摸透了它的规律它就很听话。湿对流层量小但脾气古怪，你永远不知道它下一刻要搞什么鬼。
驱散电离层的幽灵
穿过对流层之后，信号还会遇到第二道关卡，电离层。
电离层是地球大气层上部的一个区域，那里太阳辐射把气体分子电离成了自由电子和离子。当微波信号穿过这片充满了自由电子的区域时，电子会跟电磁波发生相互作用，导致信号速度变慢、路径弯曲。
在太阳活动高峰年，这项延迟误差可以达到 50 厘米。五十厘米，对于追求亚厘米精度的测高来说，这是不能容忍的。
但这里有个很妙的物理特性。电子对电磁波的延迟作用跟频率有关，频率越高影响越小。而且这个关系是可以精确描述的。
于是 TOPEX/Poseidon 做了一个天才般的设计决策，双频高度计。它同时用 
波段和 
波段发射两个不同频率的雷达波。这两个频率受到的电离层延迟不一样，但你只要测出两者的距离差 
，就可以直接在物理层面把电子的干扰给算出来并扣除掉。
经过双频校正之后，电离层的残余误差被无情压缩到了大约 0.5 厘米的极值水平。
这就是物理学的魅力。你不是在跟误差硬碰硬，你是利用误差本身的规律来消灭误差。
海面的哈哈镜
好了，信号终于穿透了所有大气层，到达了海面。你以为这下能看到真实的海面了吧。
并没有。因为海面本身就在骗你。
雷达照射到海面上的时候，海面不是一面平静的镜子。它有波浪，有波峰有波谷。而波峰和波谷对雷达波的反射特性是完全不同的。波峰像一个凸面镜，能把雷达波强烈地反射回卫星。波谷则倾向于把信号散射掉。
结果就是，雷达接收到的回波重心，会被波峰「拉低」。卫星感知到的海面高度，比真实的平均海面要偏低。
这种现象叫做海况偏差（Sea State Bias，SSB）。偏差量在 0 到 6 厘米之间，在高海况下甚至能达到 10 厘米。而且它跟有效波高（SWH）和风速高度相关，浪越大风越猛，这个偏差就越严重。
怎么修呢。科学家用复杂的二维查找表或者波谱模型，根据同步测量的波高和风速数据，把这个偏差从原始数据中剥离掉。校正后的残余误差大约 1 到 3 厘米。
说实话我觉得这个现象起名起得特别好，「哈哈镜」。你照哈哈镜的时候觉得自己变形了对吧，海面照雷达也是一样的道理。只是海面不会觉得好笑，它只是在忠实地按照物理定律反射信号而已。
地球在呼吸
你以为海面是一个静止的容器吗。不是的，海面每时每刻都在被各种力量揉捏。
最显眼的当然是海洋潮汐。月球和太阳的引力驱动海水周期性地涨落，幅度从零到数米不等。这个大家都有生活经验，海边的人都知道每天有涨潮落潮。但在测高数据里，潮汐不是一个可以忽略的背景噪声，它是一个必须用高精度模型（比如 FES 或者 GOT）精确剔除的大项。剔除之后残余 1 到 3 厘米。
但很多人不知道的是，除了海洋潮汐之外，还有三种更隐蔽的「形变」也在改变海面的高度。
固体潮（Solid Earth Tide）。不只是海水会被引力拉扯，整个固体地球也会在引力的作用下发生弹性拉伸。幅度大概 ±30 厘米。你没看错，脚下的陆地本身就在呼吸，只是你感觉不到而已。
负荷潮（Load Tide）。巨大的潮汐水体压在海床上，会导致地壳产生局部的形变。幅度 ±10 厘米左右。
极潮（Pole Tide）。地球自转轴在发生微小的漂移（极移），这种漂移会引起地球形状的细微变化，进而影响海面高度。幅度 ±5 厘米。
四类潮汐加在一起，从海洋到固体地球再到自转轴，整个地球系统都在参与这场周期性的形变舞蹈。每一项都必须用对应的模型精确计算并剔除。
有时候想到这些我会觉得有点震撼。我们以为我们在测量一个静止的海面，但实际上我们是在测量一个活的、呼吸着的星球。每一次测量都是在跟地球的脉搏打交道。
大气的隐形重量
还有一个很容易被忽视的效应，逆气压效应（Inverse Barometer）。
大气不仅有质量，而且这个质量是不均匀分布的。高压系统的空气更密更重，它会像一只隐形的手一样把海面向下压。低压系统则相反，海面会被抬升。
定量来说，海面气压每增加 1 hPa，海面就会被压低约 1 厘米。全球范围内这项效应的总波动范围大约 ±15 厘米。
修正方式是用 ECMWF 的全球气压场数据来做。不过有意思的是，在某些特定的海洋学研究中，科学家会有选择地保留或不剔除逆气压信号。为什么，因为他们想观测的是真实的海面对大气压力的动力响应，也就是风海流或者海平面上升的本来面目。这时候逆气压效应就不是噪声而是信号了。
这让我想到一个事，所谓的误差和信号，有时候只取决于你在研究什么。同一个物理过程，在这个问题里是需要被剔除的噪声，在那个问题里恰恰就是你想要测量的对象。科学里这种视角的转换特别常见。
边缘地带的手术刀
前面说的所有校正，都有一个默认前提，那就是卫星飞在开阔大洋上空，下方是一片均匀的海面。
但现实没那么乖。卫星经常会飞越南极海域的内陆湖泊、靠近海岸线的浅水区、或者覆盖着海冰的极地洋面。在这些区域，陆地的强反射会严重污染雷达回波，标准的波形模型直接失效。
怎么办呢。用**波形重定（Retracking）**技术。
标准算法用的是 Brown 模型来拟合回波波形，但在近岸和冰区，波形早就不是那个漂亮的钟形曲线了，它变得杂乱无章、充满毛刺。科学家开发了专门的算法，比如 ALES 和 SAMOSA+，它们不依赖完整的波形模型，而是直接提取波形的前沿部分，重新拟合出一个可靠的距离值。
通过这种方式，原本要被当作坏数据丢弃掉的近岸和极地观测点，被抢救回来了。这对于研究海岸带变化、极地冰盖融化这些前沿课题来说，简直是救命的技术。
我把波形重定称为边缘地带的「手术刀」。它不够优雅，但它能在常规方法宣告死亡的地方把数据救活。
一本总账
说了这么多零零碎碎的校正项，我来给你拉一张总账本。
你看这张表。初始破坏力从 2.3 米的干对流层到 ±15 厘米的逆气压，每一项单拎出来都比我们要测量的目标信号大几个数量级。但经过各自对应的「克制法宝」之后，残余误差全部被压缩到了 1 到 3 厘米的水平。
这就是人类遥感史上最严谨的除错工程之一。从米的偏差到厘米级的较量，每一个校正项背后都是几十年的理论积累、模型迭代和交叉验证。
从毛坯到精品
一张原始的雷达回波，到你电脑屏幕上那张漂亮的全球海面异常图，中间经历了什么。
L1A 级，原始回波。包含最原始的雷达反射数据和仪器底层标定信息。这是刚出炉的毛坯。
L1B 级，波形校准。物理学家的专属领地，完成仪器校准和波形重定。把原始回波变成可用的距离测量值。
L2 级，沿轨数据。这是核心产品，包含了真实海面高度以及所有校正项，供海洋学家直接调用。
L3/L4 级，全球网格产品。最终交付品，经过融合和交叉校准的无缝网格产品，比如著名的 DUACS 数据集。你在论文里用的那种彩色等值线图，大部分都是从这里来的。
整个过程就像一条精密的流水线，每一级都在为下一级提供更干净的输入。
还原真实的地球
我想用最后一点篇幅回到最初的问题。为什么要费尽心思去计算每一厘米的误差。
因为只有剥离了大气的折射、压平海浪的虚像、抚平潮汐的隆起，深藏在几千千米水面下的大洋环流，才会真正浮出水面。
未校正的仪器读数是一团混沌。剥掉大气和海况之后，你开始看到轮廓。再继续剥，绝对动力学的海面才显露出来。那才是真正有价值的东西，那是携带了热量、盐度、动量信息的海洋的脉搏。
校正的本质，是将人类的科学认知，转化为丈量宇宙的精度。
以上，既然看到这里了，如果觉得不错，随手点个赞、在看、转发三连吧，如果想第一时间收到推送，也可以给我个星标⭐～
谢谢你看我的文章，我们，下次再见。