热带气旋内核动力学 · 依据参考文献（MK97、Wang 2002、Chen & Yau 2001 等）
涡旋罗斯贝波
台风的平均位涡（PV）从中心向外减小。这个径向 PV 梯度对涡旋来说，就像地球涡度的南北梯度 β 对大气一样，是一种回复力。眼墙附近任何非对称扰动都会因此变成一种波：它比平均流转得慢，同时向外传播，被差动旋转越卷越紧，最后停在 PV 梯度消失的半径附近，并把能量和角动量交还给平均涡旋。
理论：Montgomery & Kallenbach 1997，PDF 第 1、10–11 页。这一比喻在原文中写作：行星涡度的经向导数被风暴涡度的径向导数取代。
动图 · 对照试验（线性无辐散正压模型，按 MK97 框架在本页实时计算）
同一个扰动，有回复力与没有回复力
左右两边是同一个台风、同一个初始扰动（波数 n、放在半径 r₀ 处的一圈正负相间的涡度异常）。左边保留平均 PV 梯度，扰动是涡旋罗斯贝波；右边人为删掉 PV 梯度那一项，只剩差动旋转的剪切，扰动只是被动地被卷成螺旋。白点随 r₀ 处的平均流转动，黑色三角标出 r₀ 处的波峰：左边三角逐渐落在白点后面，右边两者一直重合。
有 PV 梯度：涡旋罗斯贝波 t = 0 h
去掉 PV 梯度：只有被动剪切 t = 0 h
正涡度异常
负涡度异常
○ 白点：随 r₀ 处平均流转动
▲ 黑三角：r₀ 处的波峰
虚线圆：RMW
视野半径 100 km，北半球（逆时针旋转）
左侧试验 |ζ′| 的半径–时间图（纵轴 0–120 km，横轴 0–8 h）。亮带向上走 = 波包向外传播；走平 = 到了驻点半径。
紫：平均涡度 ζ̄(r)；灰：平均切向风 v̄(r)（都已归一化）；青：波对平均切向风的累计改变 Δv̄(r)（m/s，扰动振幅取平均涡度峰值的 10%）。
在这个平滑的单极涡旋里，波形只比平均流慢约 5–13%（波数越高越接近平均流），比观测到的 0.5–0.7 倍慢得少。按频散关系，本征传播速度正比于 ζ̄′，眼墙附近更陡的 PV 梯度会让波慢得更多；这一解释是我的推断。
这是示意性的数值试验，不是某篇论文的结果：单极涡旋取 Vmax = 45 m/s、RMW ≈ 39 km 的高斯涡度分布；环状涡旋取 32–40 km 的高涡度环，眼区涡度为环的 10%。模型只含水平平流（无对流、无摩擦、无辐散），加了很小的扩散以抑制越卷越细的结构。检验见文末。
产生原因之一 · 回复力
为什么它逆着平均流传播
把眼墙附近的一圈 PV 等值线“拉直”来看：横轴是方位角（平均流方向向右），纵轴指向中心（平均 PV 向上增大）。等值线被扰动成波状后，被推向外的地方带出了高 PV（正异常），被推向内的地方带进了低 PV（负异常）。
正异常诱生气旋式环流，负异常诱生反气旋式环流。两者在“+”的下游、“−”的上游一侧共同产生向内的流，把那里的等值线推向内，于是“−”的位置向上游移动；“+”同理。波形整体相对平均流向后传播，所以涡旋罗斯贝波转得比平均流慢。这与行星罗斯贝波西传是同一个机制，只是 β 换成了 −dζ̄/dr。示意图为概念图。
频散关系（MK97 用 WKB 方法导出）
ω = n Ω̄ + (n / R) · ζ̄′ / (k² + n²/R²)式 (12)，PDF 第 10 页。ω 为频率，n 为方位波数，k 为径向波数，R 为参考半径，Ω̄ 为平均角速度，ζ̄′ = dζ̄/dr。单极涡旋 ζ̄′ < 0，第二项为负：波相对局地角速度“倒退”。
Cpλ = R Ω̄ + ζ̄′ / (k² + n²/R²)方位相速度，式 (16)，PDF 第 11 页：等于平均流速度加上一个负的本征传播速度。
k(t) = k₀ − n t Ω̄′式 (13)，PDF 第 10 页：差动旋转（Ω̄′ < 0）让径向波数随时间增大，扰动被越卷越紧，成为拖尾螺旋。
Cgr = −2 k n ζ̄′ / [R (k² + n²/R²)²]径向群速度，式 (17)，PDF 第 11 页：拖尾螺旋 k > 0、ζ̄′ < 0 ⇒ Cgr > 0，能量向外传。k 越大 Cgr 越小（∝ k⁻³），所以波包会停下来。
rs = R + ζ̄′₀ / [R Ω̄′₀ (k₀² + n²/R²)]驻点半径（t → ∞），PDF 第 11 页。原文强调：驻点半径的存在意味着这些波被限制在涡旋附近，不能辐射到无穷远。
产生原因之二 · 谁制造了非对称扰动
扰动从哪里来
回复力只决定“扰动会怎样传播”。台风里不断有过程制造非对称 PV，它们一出现就以涡旋罗斯贝波的方式演变：
非对称对流加热
Wang 2002 的模拟中，波的上升运动在方位和径向上都超前气旋性涡度约 1/4 波长，说明与低层辐合、上升运动耦合的对流加热是驱动这些波的力量；非绝热加热对平均和扰动 PV 收支都贡献很大。Wang 2002a，PDF 第 1 页
对流注入的 PV 异常（生成期）
在准地转模式里，用参数化的对流 PV 注入代表对流，涡旋通过吸收同号 PV 异常、排出异号异常而增强；用观测估计的 PV 注入量，可在合理时间内增强到 15 m/s。Montgomery & Enagonio 1998，PDF 第 1 页
垂直风切变造成的倾斜
涡旋倾斜可以投影成涡旋罗斯贝波，这些波绕涡旋进动（Jones 1995 的结论，转引自 Jones et al. 2009）。PV 梯度和临界半径决定倾斜能否回正。Jones et al. 2009，PDF 第 1 页
环状 PV 的正压不稳定
眼墙高、眼区低的环状 PV 在 RMW 附近满足正压不稳定条件，小扰动迅速增长，涡度聚成离散的团块，形成多边形眼墙，并与眼区中尺度涡相联系；最终重新排列成近单极的涡旋。Schubert et al. 1999，PDF 第 1 页
在上面的模型里试切到“环状”涡旋：扰动不再衰减，而是指数增长；n = 4–5 增长最快（e 折时间约 20–25 分钟）。这就是正压不稳定。切回“单极”，同样的扰动只会被卷紧、向外传、衰减。
特征
怎么认出它
| 特征 | 文献中的表述或数值 | 来源 |
|---|
| 转得比平均流慢 | 波逆着平均切向流方位传播，绕眼墙一圈的周期约为随平均流运动的气块的 2 倍 | Wang 2002a，p1 |
| Elena（1985）椭圆眼墙以局地切向风的一半旋转；内雨带方位速度约 25 m/s，约为局地切向风的 68% | Corbosiero 2006，p1 |
| 向外传播 | 逆着边界层入流向外传播；Elena 内雨带平均向外 5.2 m/s | Wang 2002a p1；Corbosiero 2006 p1 |
| 观测到的内雨带径向尺度约 10 km、向外约 4 m/s（Tuttle & Gall 1995，转引） | MK97，p1 |
| 只存在于 PV 梯度足够大的地方 | 只在离中心约 60 km 内可见，更远处平均 PV 径向梯度太弱，无法支持涡旋罗斯贝波 | Wang 2002a，p1 |
| 驻点半径 | 向外传播并在 70–90 km 停滞，那里径向 PV 梯度消失或反号 | Wang 2002b，p1 |
| 空间结构 | 以波数 1、2 为主（中低层）；气旋式向内卷的螺旋，振幅在 RMW 附近最大，随高度向外倾斜 | Wang 2002a，p1 |
| 平衡运动 | 位势与风场近于平衡：辐合的气旋式流对应低位势异常；出现在平衡模式中，因此不是重力波 | Wang 2002a p1；MK97 p24 |
| 与云带耦合 | PV 带与云带强耦合，边界层内及其顶部的 PV 异常与摩擦作用产生上升运动，形成内云带，可伸到地面以上 7 km | Chen & Yau 2001，p1、p7 |
影响
它对台风做了什么
概念示意，尺度按 1 km ≈ 2.4 px 粗略画出。各要素来自 Wang 2002b（PDF 第 1 页）和 Schubert et al. 1999（PDF 第 1 页）；把它们画在同一张图里是我的综合。
眼墙形状与中尺度涡
波绕眼墙传播，使眼墙形状改变，出现气旋式旋转的多边形眼墙。Wang 2002b p1环状 PV 不稳定时，眼墙涡度聚成团块，与眼区中尺度涡相联系。Schubert 1999 p1
内、外螺旋雨带
从眼墙发出向外传播的内螺旋雨带；波在驻点半径（70–90 km）停滞，强外雨带最常出现在那里。Wang 2002b p1PV 带与云带耦合。Chen & Yau 2001 p1
强度：两种方向都有
轴对称化时能量从非对称扰动转给平均涡旋，即使强扰动也能使平均涡旋增强。Möller & Montgomery 1999 p1持续在 RMW 附近注入涡度，加上波的向外传播，构成一种增强机制。MK97 p1但在全物理模拟中，波把角动量从眼墙送进眼区，眼区加速、眼墙附近切向风减速。Wang 2002a/b p1
眼墙破碎与恢复
外雨带扰动眼墙且与眼墙中强波同相时，眼墙可能破碎，再通过轴对称化恢复；破碎（恢复）对应减弱（增强）。Wang 2002b p1高 PV 被波向内输送并混合，提高内核 PV，使台风增强。Chen & Yau 2001 p1
在模型里看 Δv̄单极涡旋、n = 2、r₀ 取 30 km 左右时，右下图的 Δv̄ 在扰动所在半径附近为正、稍外为负，加速比减速大。MK97 在 PDF 第 15 页描述了同样的结构：Δv̄ 的零点在驻点半径，最大加速大于减速。本页模型的零点位置与 WKB 驻点公式并不完全一致，因为初始扰动没有卷紧（k₀ = 0），不满足 WKB 的前提（MK97 PDF 第 14 页讨论过这一点）。
模型说明与检验
动图背后的计算
∂ζ̂/∂t = −i n Ω̄ ζ̂ + (i n / r) ζ̄′ ψ̂ + ν ∇²ζ̂， ∇²ψ̂ = ζ̂线性化的无辐散正压涡度方程，单个方位波数 n。ψ̂ 为扰动流函数，径向扰动风 u′ = −(1/r)∂ψ′/∂λ。右边第二项就是平均 PV 梯度提供的回复力；右侧对照试验把这一项设为零。ν = 20 m²/s；网格 Δr = 0.5 km，Δt = 30 s，四阶 Runge–Kutta。WKB 近似下这一方程给出 MK97 式 (12)。
∂v̄/∂t = −(1/r²) ∂(r² ⟨u′v′⟩)/∂rMK97 式 (19)，PDF 第 14 页：二阶精度的平均流变化，用于右下图的 Δv̄。
用一个有精确解的情形检验求解器：对兰金涡，边缘波（Kelvin 波）的图样角速度精确等于 Ω₀(n − 1)/n。用光滑边缘近似兰金涡，结果如下：
| 边缘宽度 | n | 本模型 | 精确值 Ω₀(n−1)/n | 相对误差 |
|---|
| 2 km | 2 | 4.96×10⁻⁴ s⁻¹ | 5.00×10⁻⁴ s⁻¹ | 0.7% |
| 2 km | 3 | 6.59×10⁻⁴ s⁻¹ | 6.67×10⁻⁴ s⁻¹ | 1.1% |
| 1 km | 3 | 6.69×10⁻⁴ s⁻¹ | 6.67×10⁻⁴ s⁻¹ | 0.4% |
剩余误差主要来自光滑边缘与理想阶跃的差别，以及初始扰动中混入的连续谱成分。
参考文献
- Montgomery & Kallenbach 1997, QJRMS — 理论、频散关系、驻点半径、平均流变化（PDF 第 1、10、11、14、15、24 页）
- Wang 2002a, JAS, Part I — 结构、相速度、60 km 范围、PV 与能量收支（PDF 第 1 页）
- Wang 2002b, JAS, Part II — 角动量输送、多边形眼墙、70–90 km 驻点、眼墙破碎与恢复（PDF 第 1 页）
- Chen & Yau 2001, JAS — PV 带与云带耦合、频散关系的三维推广（PDF 第 1、7 页）
- Möller & Montgomery 1999, JAS — 轴对称化与平均涡旋增强（PDF 第 1 页）
- Montgomery & Enagonio 1998, JAS — 对流强迫的涡旋罗斯贝波与生成（PDF 第 1 页）
- Corbosiero et al. 2006, MWR — Elena 的椭圆眼墙与内雨带观测（PDF 第 1 页）
- Schubert et al. 1999, JAS — 正压不稳定、多边形眼墙与 PV 混合（PDF 第 1 页）
- Jones et al. 2009, JAS — 倾斜投影为涡旋罗斯贝波、回正（PDF 第 1 页）