观潮 · 模式预报与气象 AI 观察
声重力波反演绕开不确定性原理：三个半波定出局地波参数，DE2 第 8132 圈极区剖面验证
> 低轨卫星沿轨道测到的温度、密度、垂直速度扰动，是声重力波诊断的关键数据——但传统参数反演要么用平面单色波近似，要么用小波变换做局地谱分析。小波谱图擅长给出"哪段空间里有波"，却受不确定性原理的硬约束：谱分辨率提上去，空间定位必然变差；而单色波近似又只适用于理想情形，对波包这种紧凑非单色结构，反问题几乎没人做过。这篇换了个思路：不再做积分变换，直接把连续波剖面切成一段段"波段"，每个波段只看中心半波加两侧反相半波，用这三个极值点的位置和振幅解析解出波包的全部局地参数。用于 DE2 卫星三圈极区轨道的实测剖面，重构波形与观测吻合；相位分析还发现，分量之间的真实物理相位差是"经典单色相位差"与"波形几何相位差"的加权和——后者最大能贡献 10—15 度的修正。三圈轨道的相位点全部落在相位诊断图的重力波区，对应以小角度近水平传播的中尺度重力波。
缺的是哪一块
低轨卫星在 250—500 km 高度上的轨道速度远超典型声重力波的相速，测到的时间序列因此可以可靠地当成波场的一维瞬时空间快照——这是卫星资料做 AGW 诊断的物理基础。但要把它用足，现有的两套工具各有短板：
方法 · 能给什么 · 卡在哪
平面单色波近似 · 反问题已成熟：由两个相位差即可完整重建谱参数 · 只对自由传播的单色波成立，处理不了波包
小波 / 加窗 Fourier · 局地谱成分的高效估计 · 不确定性原理：谱分辨率与空间定位此消彼长；交叉小波给出的是宽而模糊的相干区，相位差在有限空间域上被平均
交叉小波 / 相干性 · 尺度依赖的相位差与相干估计 · 同上，是平均量而非局地量
还有一个结构性问题：完整刻画波包必须同时分析多个波场分量（温度、密度、速度分量），传统框架下要先比对各自的小波谱图确认谱成分一致，再用交叉小波算相位——每一步都在做空间平均。
反问题的现状是：从空间剖面直接重建波包的完整谱成分、传播方向和相位结构，此前只在平面单色波近似下被求解过；对更真实的紧凑非单色波包，文献几乎空白。
三个半波怎么换出五个参数
方法的第一步是分段。对每个剖面，先找全局极值作为第一个局部波包的主半波，再检查两侧是否各有一个反号的等距伴半波——等距性按后文给出的定量对称判据核验，不靠目测。三个关键半波的空间坐标和振幅记录下来后，对剖面剩余未分段区间循环重复。在 DE2 第 8132 圈的实例里，这样切出了 4 个物理波包。
第二步是把波包写成小波形表示：把空间坐标原点移到波包中心。这个平移会改变波相位，得到一个直接反映小波内部对称性的"几何相位"。第三步用高斯函数近似缓变包络——这是描述真实大气波包包络的标准选择——由此建立到修正 Morlet 小波的解析桥。
于是整个反演只需要三个主导半波的位置和振幅，就能解出波包的五个参数：
参数 · 性质 · 是否随分量变化
载波波数 · 谱参数 · 同一波事件的所有分量必须相同
谱带宽参数 · 谱参数 · 同上
包络空间中心 · 几何参数 · 逐分量确定
包络峰值振幅 · 几何参数 · 逐分量确定
几何相位 · 几何参数 · 逐分量确定
两个谱参数刻画所在波事件的波谱，必须对所有分量一致——这本身就是"这些分量来自同一个波事件"的判据。三个几何参数刻画局地形状，逐分量单独确定。
值得强调的是方法不锁死在高斯包络上：任何有明确物理动机的钟形包络都可以替换进去，五个参数的集合不变，只是反演公式不同。换句话说，"从三个主导半波重建波包参数"这套几何做法是普适的，差异只在解析表达式。
相位偏移要修一项几何修正
到此为止各分量是独立分析的。要拿相位差做物理诊断，必须从几何相位回到初始参考相位。结论是：分量之间的真实物理相位偏移是两项的加权和——
- 经典单色相位差：两个分量中心极值之间的空间距离除以波长再乘系数，单色极限下权重趋于 1；
- 几何相位差：各小波对称性之间的相对偏移，波形完全相同时这一项为零。
经典项占主导，几何项给出一个中等但可观的修正，最大到 10—15 度。换句话说，直接从观测波形上"看出来的"相位差，与真实的物理相位差略有不同——只要分析是建立在分量间相位差上的，这项修正就不能省。
修正后的相位差才有资格当诊断输入。对理想分层大气（等温、无风、无耗散、单原子气体近似——用于卫星高度上原子氧占主导的中间层顶以上）中的自由传播单色 AGW，作者此前的工作已证明：只用垂直速度与温度（或密度）相对涨落之间的两个相位差，就能完整重建水平波长、垂直波长和内频；且这两个相位差不依赖坐标系选取。据此画出的相位诊断图，把声波区与重力波区、向上与向下传播分在平面上的不同区域——算出两个相位角，在图上落一个点，波型和垂直传播方向立即读出，其余谱参数再解析重建。
结果：三条极区轨道都指向近水平传播的中尺度重力波
验证用 DE2 卫星三圈过极轨道：8132、4820、8129。8132 圈直接把算出的相位差代进解析公式；另两圈利用"两个剖面近乎同步、另一个相位差接近某定值"的观测特征简化公式。
轨道 · 相位点位置 · 结论
8132 · 四个相位点全部落入重力波区下半部 · 中尺度重力波，与水平面成中等角度传播
4820 · 落在同一重力波区几何中心附近 · 中尺度重力波，近水平传播
8129 · 同上 · 中尺度重力波，近水平传播
重构精度用"逐段重构再求和"的方式核验：对每个孤立波段取回波包参数、代入理论公式重建局地波形，所有波段的贡献求和后与原始观测剖面比对。第 8132 圈的三个分量（垂直速度 Vz、密度相对涨落、温度相对涨落）的重构曲线与观测吻合，全部振荡都被包络在理论导出的包络之内。
极区 AGW 的这类谱特征此前已有报道——Johnson 等、Innis 与 Conde、Fedorenko 等都用 DE2 极轨资料分析过，Fedorenko 等把极区观测特征和大振幅归因于非均匀风场对波谱的过滤。本文的结论与之一致，但方法路径不同：那些工作全部在单色平面波近似下进行，本文不依赖该近似。
原文图
第 8132 圈实测剖面与分段结果：上图为 900T′/T 与垂直速度 Vz，下图为 300ρ′/ρ 与 Vz；圆点标出各波段的中心极值与两侧反相伴极值，方框编号 ①—④ 为识别出的四个物理波包，横轴为世界时 UT。
波传播几何示意：波矢 k 与卫星轨迹夹角 θ 决定了测得波长被拉伸的程度，λ 为真实波长、Δλ 为几何拉伸量——这就是"沿轨迹一维剖面"必须做几何修正的原因。
修正 Morlet 小波与其三个主导半波：中心极值 0 与两个反号伴极值 −1、+1 的位置和振幅，唯一确定小波的五个参数；虚线为高斯包络。
第 8132 圈三个分量（Vz、300ρ′/ρ、900T′/T）的观测（实线）与理论重构（虚线）比对，粉色为理论包络；重构曲线贴合观测，全部振荡落在包络内。
第 4820 圈（左）与 8129 圈（右）的原始剖面：400ρ′/ρ、1000ρ′/ρ 与 Vz、1000T′/T 的叠加显示，两圈上密度与温度分量近乎同步，Vz 与它们之间存在明显的相位差。
第 4820 圈（左，500ρ′/ρ 与 500T′/T）与 8129 圈（右，1000ρ′/ρ 与 1000T′/T）的重构（虚线）与观测（实线）比对，横轴为磁地方时；两圈的重构同样贴合观测。
作者承认的限制
- 复杂失真段必须剔除。空间上复杂或高度畸变的波段被排除在分析之外；在实例里它们位于整个区间的外边界，但作者明确指出：实际观测资料里这类失真段可能出现在剖面的任何位置。
- 失真区里方法用不了。这些区间内波包可能重叠，信号也可能被局地噪声或其他物理性质的波过程严重污染——这种物理上的模糊性限制了方法的直接应用，需要预处理滤波或替代处理技术。
- 时间畸变没有消除，只有界定。低轨卫星轨道速度远超 AGW 相速这一前提，把时间序列当作瞬时空间快照是近似；波沿卫星运动方向还是逆着传播，会造成不同的时间畸变，作者只说"必须计入"，没有给出修正量。
- 理论模型是理想化大气的。等温、常化学成分、无风、无耗散、单原子气体近似——所有谱参数重建都建立在这套 Hines 模型框架上，实测大气不满足这些条件。
- 只有一维切面。方法给出的是沿卫星轨道方向的谱成分；垂直结构要靠相位诊断图反演，而不是直接测得。
- 后续应用尚未发表。把这套方法用到其他波过程的分析，作者写明将在后续论文中详述——本文只处理了理想分层大气这一种理论框架。
一句话判断
把"积分变换找谱"换成"形态学读波形"，代价是必须先信任分段判据，收益是绕开了不确定性原理的硬约束、拿到真正的局地参数——对卫星这种天生只有一维切面的观测，这笔账是划算的；但复杂失真段和重叠波包仍然要用别的方法先处理掉。