企业官网建设流程全解析

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通量则像是考虑了声源在风中运动，或者声音在温度梯度（导致声速变化）的大气中传播的通用公式，适用于更复杂的非均匀环境。

在实际科研和业务中：

1. 对于大尺度、准地转的波动诊断（如平流层动力学研究），E-P通量及其散度仍是黄金标准。
2. 对于天气尺度扰动的传播路径的快速、直观诊断（如看一波槽脊系统如何向下游传播），Plumb通量因其计算简便而广受欢迎。
3. 当研究焦点位于西风急流核心区、风暴轴区域，或基本流有强切变时，必须使用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。

1. 计算扰动场：ψ‘ = (Φ - [Φ]) / f，其中[Φ]是纬向平均位势高度。
2. 计算空间导数：需要∂ψ‘/∂x（经向风扰动相关）和∂ψ‘/∂y（负的纬向风扰动相关）。使用中心差分法，注意处理东西边界（通常采用循环边界）和南北边界（可能需要用单边差分或设为NaN）。
3. 组装通量分量： [ 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上的值）进行归一化，使箭头大小易于可视化。
4. 可视化：使用箭矢图quiver绘制(F_x, F_z)。通常还会叠加扰动位势高度场Φ‘的等值线，以观察波动槽脊位置与通量箭头的关系。

### 3.3 Takaya-Nakamura通量（二维）计算步骤

T-N通量计算更复杂，因为它需要基本气流的二阶导数。我们考虑在球坐标下的压力坐标系（常用）中的形式。

1. 准备基本流场：这里的基本流[u],[v]是纬向平均流。但注意，T-N通量允许[v] ≠ 0，即平均经圈环流。
2. 计算绝对涡度：η = f + ζ，其中相对涡度ζ = ∂[v]/∂x - ∂[u]/∂y（这里用的是基本流）。需要计算f和β。
3. 计算扰动流函数ψ‘：同上，ψ‘ = Φ‘ / f。
4. 计算关键中间量：这是最容易出错的一步。需要计算：
   • 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。
5. 组装通量向量：根据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等压面进行分析。

1. 数据准备：提取位势高度Z、纬向风U、经向风V、温度T。进行10天高通时间滤波（保留10天以下波动），得到天气尺度扰动场Z‘,U‘,V‘。计算纬向平均得到基本流[U],[V]。
2. 分别计算三种通量的水平分量：
   • E-P通量：使用滤波后的扰动场U‘,V‘,T‘，按照公式计算F_x和F_y。这里F_y主要体现动量通量U‘V‘。
   • Plumb通量：由扰动位势高度Φ‘ = g * Z‘计算扰动流函数ψ‘，进而计算其空间导数，得到通量。
   • T-N通量：使用纬向平均基本流[U],[V]和扰动流函数ψ‘，严格按照公式计算包含基本流切变项的通量。
3. 可视化与解读：
   • 第一张图：叠加300hPa扰动位势高度Z‘（填色）和Plumb通量矢量（箭头）。可以清晰看到，从北大西洋阻塞高压下游（东侧），有一支强烈的波作用通量箭头呈大圆路径向东南方向传播，经地中海、西亚，最终指向东亚沿岸。这与观测到的长波槽东移路径高度吻合。通量箭头在波槽区（Z‘负值区）呈辐合状，在波脊区（Z‘正值区）呈辐散状，符合理论预期。
   • 第二张图：在同一张图上对比Plumb通量（红色箭头）和T-N通量（蓝色箭头）。在中纬度西风急流轴（由[U]等值线标出）附近，两者出现了显著差异。Plumb通量箭头在急流轴上相对平直东传；而T-N通量箭头在急流入口区（北大西洋）显示出更强的向极地分量，在急流出口区（东亚沿岸）则显示出更强的向赤道分量。这表明，在急流的强切变区域，波动传播路径发生了折射，T-N通量更好地捕捉了这一物理事实。
   • 第三张图：绘制E-P通量的水平分量及其散度∇·F。可以看到，在东亚大槽加深的区域（寒潮关键区），∇·F呈现显著的负值（蓝色区域），意味着波活动在此辐散，对西风气流施加了一个向西的负强迫。这正好对应了槽后西北气流加强、引导冷空气南下的动力过程。而在大西洋阻塞高压区上游，∇·F为强正值（红色区域），表示波活动汇聚，对西风施加正强迫，有利于高压脊的维持。

### 4.3 诊断结论与预报启示

通过这个案例，我们验证了：

1. Plumb通量能清晰、直观地展示天气尺度波动的能量传播大圆路径，是追踪波列传播的优良工具，计算也相对简单。
2. Takaya-Nakamura通量在急流等强切变区域提供了更精确的波传播路径诊断，揭示了波动受基本流折射的细节，这对于理解波动如何与急流相互作用、影响下游天气至关重要。
3. 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通量进行复核和深入分析，两者结合能给出更可靠的物理图像。