1. 从天气预报的“误差”说起为什么需要理解Rossby波作用通量如果你关注过中长期的天气预报可能会发现一个有趣的现象有时候数值预报模型对一周后的天气形势预测得相当准确但有时候哪怕只是三四天后的预报也会出现巨大的偏差。这种不确定性很大程度上源于大气中一种大尺度的、蜿蜒曲折的波动——Rossby波。Rossby波是大气和海洋中能量和动量传输的关键载体它决定了阻塞高压、切断低压等持久性天气系统的生消直接影响着我们的寒潮、暴雨和持续高温天气。然而仅仅知道“有Rossby波”是远远不够的就像我们知道汽车有发动机但如果不清楚发动机的功率和扭矩如何传递到车轮就无法理解汽车为何能加速或爬坡。在动力气象学的研究和业务预报中我们迫切需要一种工具来定量诊断Rossby波的能量是如何传播、在哪里汇聚或发散的。这就是“Rossby波作用通量”概念的价值所在。它不是一个单一的公式而是一组强有力的诊断工具能够将复杂的波动三维结构浓缩为直观的矢量箭头清晰地告诉我们波能量从何处来、到何处去。理解并正确计算这三种通量是深入解读数值预报产品、改进模式物理过程参数化、乃至理解气候变化中大气环流异常的核心技能。本文将从一个预报员或研究者的实操视角深入拆解大气动力学中最常用的三种Rossby波作用通量Eliassen-Palm通量、Plumb通量和Takaya-Nakamura波作用通量。我们将避开繁琐的纯理论推导聚焦于它们的物理图像、计算公式的“零件拆解”、编程实现中的关键细节以及在实际天气个例分析中如何解读这些图形让你不仅能看懂文献中的图更能自己动手算、自己分析用。2. 物理内核三种通量究竟在“通”什么在深入公式之前我们必须建立清晰的物理图像。Rossby波作用通量本质是描述“波作用量”的输送。波作用量是一个守恒量在无耗散、无强迫的基本气流中可以粗略理解为波动的“能量护照”。它的通量即代表了波能量传播的方向和强度。三种通量公式形式不同源于它们基于不同的理论框架和近似适用于不同的场景。### 2.1 Eliassen-Palm通量奠基与准地转框架下的经典Eliassen-Palm通量是历史上最早提出且应用最广泛的波作用通量。它建立在准地转理论框架下。这个框架做了关键简化过滤掉了声波和重力波专注于描述大尺度、慢变的Rossby波过程同时认为大气是层结稳定、近似处于地转平衡和静力平衡状态的。这非常符合我们对中高纬度大尺度运动的第一直觉。E-P通量是一个三维矢量F。它的三个分量有明确的物理意义F_x (纬向分量)主要与经向的热量输送相关。它可以写成 -ρ₀ *v‘θ‘的形式其中v‘是经向风扰动θ‘是位温扰动。当波动使得暖空气向北v‘0, θ‘0或冷空气向南输送时这个分量不为零。它反映了波动通过经向热输送对基本西风动量的影响。F_y (经向分量)主要与动量的经向输送相关但形式稍复杂。它正比于u‘v‘纬向风和经向风扰动的协方差。这个分量直接体现了波动动量通量的辐合辐散是导致西风急流加速或减速的直接动力因素。F_z (垂直分量)在准地转框架下它与经向风扰动和垂直速度扰动或位温扰动的协方差有关描述了波动能量在垂直方向的传播对于理解波动是否能够向上传播到平流层至关重要。注意在多数教科书的推导中E-P通量的散度 ∇·F有一个极其重要的性质它正比于波动对基本西风气流的净强迫。∇·F 0 表示波活动在此地汇聚对西风气流施加一个向东的加速作用∇·F 0 表示波活动在此地辐散对西风气流施加一个向西的减速作用。这是E-P通量图解的核心——我们不仅看箭头方向更要看箭头是“发散”还是“汇聚”。### 2.2 Plumb通量从小振幅到有限振幅的推广Plumb通量可以看作是E-P通量在小振幅波动假设下的一个特例但它更重要的价值在于其清晰的物理表述和易于计算的优点。Plumb在1985年的经典论文中从波作用量守恒方程出发在“缓变基本流”和“局地平面波”的假设下推导出了一套更直观的通量公式。Plumb通量的核心思想是对于一个单色波具有局地波数k, l, m其波作用通量方向与局地群速度方向一致并且大小正比于波作用密度。它的计算公式直接使用了扰动场形式相对简洁。以二维经度-高度剖面为例其分量常写为[ F^{(p)}_x \propto \frac{1}{2} \frac{\partial \psi‘}{\partial x} \frac{\partial \psi‘}{\partial y}, \quad F^{(p)}_z \propto \frac{1}{2} \frac{f^2}{N^2} \frac{\partial \psi‘}{\partial x} \frac{\partial \psi‘}{\partial z} ]其中ψ‘是扰动流函数。可以看到它完全由扰动量的空间导数构成计算时无需显式涉及基本气流基本流的信息隐含在扰动结构中了。这使得Plumb通量在诊断观测资料或模式输出时非常方便特别是当基本气流定义模糊比如存在强烈的水平切变时。### 2.3 Takaya-Nakamura波作用通量适用于基本气流非均匀的利器前两种通量尤其是E-P通量在基本气流如西风均匀或缓变的假设下效果最好。然而现实大气中的西风急流具有很强的水平和垂直切变是非均匀的。在急流入口区、出口区或急流轴附近波动特性会发生剧烈变化传统的E-P或Plumb通量可能无法准确描述波动的传播。Takaya和Nakamura在1997年提出的波作用通量正是为了解决这一问题。他们的推导基于线性化的原始方程并且没有假设基本气流是均匀的。T-N通量是一个三维矢量其表达式包含了基本气流的二阶导数项能够更好地捕捉在非均匀基本流中波动的折射和传播效应。T-N通量的公式较为复杂但其物理意义非常深刻它定义了在非均匀流中一个与局地群速度平行的广义波作用通量。当基本气流均匀时它会退化为类似Plumb通量的形式。因此T-N通量被认为是更通用、更稳健的诊断工具特别适用于分析中纬度天气尺度波动与急流之间的相互作用例如罗斯贝波列在穿越急流时的行为。### 2.4 三者的关系与选用指南我们可以用一个简单的类比来理解三者的关系研究声音在空气中的传播。E-P通量像是用流体力学方程严格推导出的声能通量考虑了空气的密度和状态方程非常严谨但计算复杂。Plumb通量像是在均匀、静止空气中测量声强的简化公式计算简单直观在均匀环境下很准确。T-N通量则像是考虑了声源在风中运动或者声音在温度梯度导致声速变化的大气中传播的通用公式适用于更复杂的非均匀环境。在实际科研和业务中对于大尺度、准地转的波动诊断如平流层动力学研究E-P通量及其散度仍是黄金标准。对于天气尺度扰动的传播路径的快速、直观诊断如看一波槽脊系统如何向下游传播Plumb通量因其计算简便而广受欢迎。当研究焦点位于西风急流核心区、风暴轴区域或基本流有强切变时必须使用T-N通量才能得到可靠的波传播路径诊断。否则可能会误判波能量的来源和汇。3. 从公式到代码关键参数与计算陷阱理论清晰后下一步就是动手计算。这里以最常用的Plumb通量二维和T-N通量二维为例详解编程实现中的“魔鬼细节”。我们假设你已有处理好的网格化数据如ERA5再分析资料变量包括纬向风u、经向风v、位温theta或温度T、位势高度z等。### 3.1 数据预处理滤波与基本流分解计算波作用通量的第一步永远是从总场中分离出“基本流”和“扰动”。这不是简单的气候态扣除对于天气尺度过程通常采用时空滤波。时间滤波要诊断特定时间尺度如天气尺度3-10天的波动需要对时间序列进行带通滤波。常用的是Lanczos滤波器或Butterworth滤波器。例如使用scipy.signal中的函数。关键点滤波会产生边界效应通常需要舍弃前后一定时长的数据如滤波器阶数的一半。空间滤波有时需要分离出纬向平均流和扰动。基本流[X]通常取为沿纬圈的平均扰动X‘ X - [X]。对于T-N通量这里的基本流[X]就是公式中需要的非均匀基本流。变量计算需要准备好位温θ、浮力频率N²、科里奥利参数f及其β效应β。N² (g/θ) * dθ/dz这里dθ/dz需要中心差分计算注意垂直坐标是几何高度还是气压坐标需要进行转换。### 3.2 Plumb通量二维计算步骤假设我们在一个经度-压力或高度剖面上计算变量是扰动流函数ψ‘。ψ‘可以通过扰动地转风关系从扰动位势高度Φ‘得到u‘_g - (1/f) ∂Φ‘/∂y,v‘_g (1/f) ∂Φ‘/∂x。在准地转框架下流函数ψ‘ Φ‘ / f。计算扰动场ψ‘ (Φ - [Φ]) / f其中[Φ]是纬向平均位势高度。计算空间导数需要∂ψ‘/∂x经向风扰动相关和∂ψ‘/∂y负的纬向风扰动相关。使用中心差分法注意处理东西边界通常采用循环边界和南北边界可能需要用单边差分或设为NaN。组装通量分量 [ F_x \frac{1}{2|\vec{U}_s|} \left[ \left( \frac{\partial \psi‘}{\partial x} \right)^2 - \psi‘ \frac{\partial^2 \psi‘}{\partial x^2} \right] \quad \text{(简化形式之一具体形式需查原文)} ] 更常见的操作化公式是直接使用导数乘积形式并乘以一个参考尺度的基本流速度U_s通常取剖面纬向平均风的某个特征值如300hPa上的值进行归一化使箭头大小易于可视化。可视化使用箭矢图quiver绘制(F_x, F_z)。通常还会叠加扰动位势高度场Φ‘的等值线以观察波动槽脊位置与通量箭头的关系。### 3.3 Takaya-Nakamura通量二维计算步骤T-N通量计算更复杂因为它需要基本气流的二阶导数。我们考虑在球坐标下的压力坐标系常用中的形式。准备基本流场这里的基本流[u],[v]是纬向平均流。但注意T-N通量允许[v] ≠ 0即平均经圈环流。计算绝对涡度η f ζ其中相对涡度ζ ∂[v]/∂x - ∂[u]/∂y这里用的是基本流。需要计算f和β。计算扰动流函数ψ‘同上ψ‘ Φ‘ / f。计算关键中间量这是最容易出错的一步。需要计算U [u] / cosφ(φ是纬度)V [v] / cosφ基本流的变形张量项如∂U/∂x,∂U/∂y,∂V/∂x,∂V/∂y。绝对涡度的梯度∂η/∂x,∂η/∂y。静力稳定度参数S (R/H) * (∂[T]/∂p - κ[T]/p)其中H是尺度高度κR/c_p。组装通量向量根据Takaya Nakamura (2001, JAS) 的公式(3.5)通量分量(F_x, F_y)由以下项组成 [ F_x \propto \frac{1}{2|\vec{U}_s|} \left[ U (\psi‘x^2 - \psi‘ \psi‘{xx}) V (\psi‘_x \psi‘y - \psi‘ \psi‘{xy}) \text{包含} \eta_y, S_y \text{的项} \right] ]F_y有对称的形式。编程时必须严格按照论文中的公式逐项实现并确保微分计算准确。强烈建议将计算过程封装成函数并对一个已知的简单波例如平面波进行测试验证通量方向是否与群速度方向一致。### 3.4 编程实战中的避坑指南差分方案的选择对于二阶导数中心差分是必须的。对于一阶导数在边界处可使用向前/向后差分但内部点务必用中心差分以保证精度。使用numpy.gradient函数可以方便地计算中心差分。单位制统一确保所有物理量使用国际单位制米、秒、开尔文、帕斯卡。经纬度需要转换为弧度距离需要转换为米使用地球半径。垂直坐标插值再分析资料通常是等压面坐标但公式中可能涉及对高度的导数。有时需要将数据从气压坐标插值到对数压力坐标或几何高度坐标这需要谨慎处理避免引入虚假梯度。可视化缩放通量箭头通常数量级很大直接绘制会重叠混乱。标准做法是进行归一化处理例如除以某个特征速度或通量的最大值或者使用quiver的scale和width参数精细调整。同时绘制通量矢量的散度场着色图叠加在箭头图上是解读动力强迫的关键。结果的物理合理性检查计算完成后问自己几个问题通量箭头是否大致垂直于波动等值线在波源地区如地形强迫上游是否呈现辐散状在西风急流核心区T-N通量是否显示出与Plumb通量不同的路径这些检查能帮你快速发现计算错误。4. 案例诊断一次东亚寒潮过程中的波流相互作用让我们用一个具体的例子将上述知识串联起来。分析2023年1月一次影响中国的强寒潮过程。我们使用ERA5逐日再分析资料水平分辨率0.25°等压面资料。### 4.1 天气过程简述与问题提出在寒潮爆发前约一周北大西洋上空出现了一个强烈的阻塞高压导致上游环流经向度加大。随后一个长波槽在东移过程中不断加深最终引导极地冷空气大举南下。我们的科学问题是这次寒潮相关的Rossby波能量是如何从北大西洋传播到东亚地区的三种通量给出的传播路径诊断有何异同### 4.2 计算与对比分析我们选取寒潮爆发前5天波动发展期的300hPa等压面进行分析。数据准备提取位势高度Z、纬向风U、经向风V、温度T。进行10天高通时间滤波保留10天以下波动得到天气尺度扰动场Z‘,U‘,V‘。计算纬向平均得到基本流[U],[V]。分别计算三种通量的水平分量E-P通量使用滤波后的扰动场U‘,V‘,T‘按照公式计算F_x和F_y。这里F_y主要体现动量通量U‘V‘。Plumb通量由扰动位势高度Φ‘ g * Z‘计算扰动流函数ψ‘进而计算其空间导数得到通量。T-N通量使用纬向平均基本流[U],[V]和扰动流函数ψ‘严格按照公式计算包含基本流切变项的通量。可视化与解读第一张图叠加300hPa扰动位势高度Z‘填色和Plumb通量矢量箭头。可以清晰看到从北大西洋阻塞高压下游东侧有一支强烈的波作用通量箭头呈大圆路径向东南方向传播经地中海、西亚最终指向东亚沿岸。这与观测到的长波槽东移路径高度吻合。通量箭头在波槽区Z‘负值区呈辐合状在波脊区Z‘正值区呈辐散状符合理论预期。第二张图在同一张图上对比Plumb通量红色箭头和T-N通量蓝色箭头。在中纬度西风急流轴由[U]等值线标出附近两者出现了显著差异。Plumb通量箭头在急流轴上相对平直东传而T-N通量箭头在急流入口区北大西洋显示出更强的向极地分量在急流出口区东亚沿岸则显示出更强的向赤道分量。这表明在急流的强切变区域波动传播路径发生了折射T-N通量更好地捕捉了这一物理事实。第三张图绘制E-P通量的水平分量及其散度∇·F。可以看到在东亚大槽加深的区域寒潮关键区∇·F呈现显著的负值蓝色区域意味着波活动在此辐散对西风气流施加了一个向西的负强迫。这正好对应了槽后西北气流加强、引导冷空气南下的动力过程。而在大西洋阻塞高压区上游∇·F为强正值红色区域表示波活动汇聚对西风施加正强迫有利于高压脊的维持。### 4.3 诊断结论与预报启示通过这个案例我们验证了Plumb通量能清晰、直观地展示天气尺度波动的能量传播大圆路径是追踪波列传播的优良工具计算也相对简单。Takaya-Nakamura通量在急流等强切变区域提供了更精确的波传播路径诊断揭示了波动受基本流折射的细节这对于理解波动如何与急流相互作用、影响下游天气至关重要。E-P通量及其散度则直接定量揭示了波动对平均西风气流的动力强迫位置和强度将环流变化槽脊发展与动力过程波流相互作用直接联系起来。对于预报员而言在分析中期预报时可以计算Plumb或T-N通量来预判上游波动的能量是否会向下游如我国传播。如果发现有一支强而持续的通量箭头从关键区指向预报区域就需要高度关注未来该区域环流发生剧烈调整如槽加深、脊发展的可能性这往往是寒潮、暴雨等极端天气的前兆信号。5. 进阶讨论边界、瞬变与气候尺度应用掌握了基础计算和天气尺度诊断后我们可以将视野拓展到更复杂和更宏观的场景。### 5.1 球面几何与边界条件处理全球大气是球面在极区和高纬度使用笛卡尔坐标近似会引入较大误差。因此在实际的全球诊断中所有导数都应在球坐标下进行。这意味着经向导数∂/∂y应转换为(1/R) ∂/∂φ其中R是地球半径φ是纬度弧度。纬向导数∂/∂x应转换为(1/(R cosφ)) ∂/∂λ其中λ是经度弧度。面积元在计算通量散度进行区域积分时必须考虑球面面积元cosφ dφ dλ。在南北边界通常取为85°N和85°S避免极地奇点通量的计算需要特别处理。一种常见做法是在边界附近使用单边差分公式或者直接将这些格点的通量值设为缺失值在可视化时进行掩膜。### 5.2 瞬变波与定常波分解我们之前的分析针对的是瞬变波天气尺度扰动。但在气候学中定常波由地形和海陆热力差异强迫产生的、相位固定的波动同样重要。诊断定常波的传播需要将时间平均场如冬季平均视为“基本流”将气候态场相对于纬向平均的偏差视为“定常扰动”然后代入通量公式计算。这时计算出的通量描述的是气候平均意义上波能量的输送。例如用它来诊断青藏高原对大气环流的定常强迫效应是如何向上游和下游传播的。### 5.3 在气候变化研究中的应用波作用通量是连接热带外大气环流异常与遥相关型的桥梁。在气候变化研究中常被用于诊断ENSO的影响厄尔尼诺事件导致的赤道太平洋热力异常如何通过激发罗斯贝波列如太平洋-北美型PNA影响中高纬度天气。北极放大效应北极快速增暖如何通过改变经向温度梯度进而影响西风急流和罗斯贝波的稳定性导致中纬度极端天气事件增多。外强迫响应分析气候模式模拟的温室气体增加情景下定常波通量的变化从而理解未来区域气候如北美西岸降水、东亚季风可能如何改变。在这些应用中通常计算的是气候平均场如30年冬季平均的定常波通量或者计算异常场如厄尔尼诺年减去拉尼娜年的通量以揭示外强迫导致的波传播路径的变化。### 5.4 一个常被忽视的要点通量的“非唯一性”最后必须指出一个重要的理论问题波作用通量的定义不是唯一的。从不同的波动方程准地转、原始方程、浅水方程出发或者对波作用量密度进行不同的变换都可能导出形式上不同的通量表达式。E-P、Plumb、T-N通量只是众多可能形式中经过实践检验、物理意义最清晰的几种。它们在小振幅、缓变基本流的假设下其散度项∇·F所表征的波流相互作用物理是一致的。但在有限振幅或强非线性的情况下不同通量形式给出的具体矢量方向可能会有差异。因此在报告结果时明确说明你所使用的是哪一种通量及其依据是严谨的科学态度。在实际工作中我通常会先用Plumb通量做快速诊断如果发现关键区域如急流附近的传播特征复杂则会用T-N通量进行复核和深入分析两者结合能给出更可靠的物理图像。