[JGR Oceans] 突破 MOST 限制：基于 Richardson 数的海气通量参数化新进展

TL;DR

本文基于南南海（South China Sea）观测平台两年的实测数据，深入剖析了低中风速下海气交换系数的复杂行为。研究发现，传统的 Monin-Obukhov 相似理论（MOST）在处理低风速时存在系统性偏差。作者提出了一套全新的、免迭代的参数化方案，通过引入简化的 Richardson 数，不仅提高了动量、热量和水汽通量的计算精度，还揭示了离岸风与向岸风在阻力系数上的显著差异。

背景定位

海气通量（Momentum, Sensible Heat, Water Vapor）是连接海洋与大气的桥梁。长期以来，学术界主要依赖 COARE 算法（基于 MOST 理论）。但在低风速或极稳定层结下，MOST 的迭代过程往往导致计算不稳定甚至不收敛。本工作属于典型的“基于观测的理论修补与创新”，在海气相互作用领域具有重要的实用价值。

痛点与动机：为什么 MOST 不够用了？

1. 自相关问题 (Self-correlation)：MOST 需要利用摩擦速度 $u_{\ast }$ 来计算 Obukhov 长度 $L$，而 $L$ 反过来又用于修正 $u_{\ast }$，这种闭环迭代在观测数据有噪声时会产生虚假的相关性。
2. 物理直觉缺失：在风速极小时（< 5 m/s），海水表面看起来像“光滑”的流体，此时阻力系数 $C_D$ 反而随风速降低而升高，传统线性模型难以捕捉这一非线性拐点。
3. 地理差异：近岸海域受 Fetch（风区长度）影响极大。离岸风形成的波浪由于演化时间短（波龄小），表面更粗糙，阻力更大；这一效应在通用算法中常被忽视。

核心方法论：免迭代的 Richardson 数模型

作者提出将交换系数 $C_x$ 表达为风速函数与稳定性修正函数的乘积： $C_x=C_{xu}(U)\cdot G_x(R_i)$

1. 架构解析

• 风速项 $C_{xu}(U)$：不再简单套用线性公式，而是针对向岸/离岸风分别拟合了二阶多项式，精准捕捉了低风速下的“粘性系数效应”。
• 稳定性项 $G_x(R_i)$：引入了一个简化的 Richardson 数 $R_i$。

上图展示了 $G_D,G_H,G_E$ 与 $R_i$ 的线性拟合关系，证明了该简化指标能有效替代复杂的 Obukhov 长度。

实验与结果分析

通过对 12,240 组动量通量、5,813 组感热通量和 5,637 组水汽通量的分析，研究得出了以下关键结论：

$C_D$ 的“三段式”行为

• 极低风速 (< 5 m/s)：$C_{DN}$ 随风速减小而增大，体现了空气动力学光滑表面的物理特性。
• 中等风速 (5-12 m/s)：$C_{DN}$ 基本保持在 0.001 左右，表现平稳。
• 较高风速 (> 12 m/s)：由于波浪破碎和 Fetch 限制，$C_{DN}$ 快速上升。

图 5 清晰地展示了 $C_D$ 随风速变化的 U 型趋势，以及实测数据与 COARE 3.0 的对比。

SOTA 对比

作者将新模型与标杆算法 COARE 3.0 做了同台竞技。结果显示，在使用相同的 Charnock 系数下，新模型在各个指标（RMS误差、Bias）上均达到或超过了 COARE 3.0。特别是在不需要繁琐迭代的情况下，计算效率提升巨大，极其适合嵌入到需要频繁调用的海洋气象耦合模型中。

深度洞察与总结

关键启示 (Takeaways)

• 波龄的重要性：离岸风因为 Fetch 短，生成的“年轻波”比向岸风带来的“ mature swells”更粗糙，这解释了为什么近岸风速相同但阻力不同的现象。
• 热量与水汽的恒定性：与 $C_D$ 不同，$C_H$ 和 $C_E$ 在中性条件下基本是常数（分别为 1.20 和 0.96），不随风速显著波动。

局限性与未来研究

尽管新模型在低中风速下表现卓越，但在强风（> 20 m/s）或台风条件下，由于缺乏极端天气观测数据，该模型的“饱和效应”仍有待进一步验证。未来的方向应考虑将该方法推广到深海大洋以及多模态的海气界面（如考虑海草、油污等对阻力的影响）。

关键词：Air-sea fluxes, Richardson number, Monin-Obukhov similarity theory, Exchange coefficients, South China Sea.