[Boundary-Layer Meteorology] 海浪如何“厚此薄彼”：动量与感热交换系数脱节的物理起源

TL;DR

在海气交互研究中，一个长期存在的现象是：随着风速增加，海洋变得越来越“粗糙”（阻力系数 $C_D$ 飙升），但其交换热量的能力（Stanton Number $C_H$）却表现得异常迟钝。本文通过构建一个考虑波浪能量输入的大气表面层 e-ε 湍流模型，揭示了这一现象背后的物理本质：动量可以直接通过海浪的“形阻”传递，而热量交换则必须卑微地依赖于被波浪变相改造后的湍流扩散。

背景定位：由于物理机制不同导致的“脱节”

在经典的气象学 bulk parameterization 中，动量应力 $au$ 和感热通量 $H$ 通常被写成：

• $au={\rho }_a{C}_D(\Delta u{)}^2$
• $H={\rho }_a{c}_p{C}_H\Delta u\Delta T$

实验数据（如 HEXMAX, Anderson 1993）清晰显示，$C_D$ 随风速增加 100%-150% 是常态，而 $C_H$ 却往往被视为常数。这种脱节在数值预报中导致了巨大的不确定性，尤其是在高风速（> 15 m/s）环境下。

核心动机：挑战传统的“固体边界”假设

作者指出，现有的研究大多将海面看作一个具有特定粗糙度高度（$z_0$ 和 $z_{0t}$）的固体表面。然而，海浪不仅仅是粗糙元，它们还是动态的能量汇。

1. 动量分配：总应力 $au$ 是由湍流应力 $a{u}^t$ 和波导应力 $a{u}^w$ 共同承担的。在强风下，$a{u}^w$（Form Drag）可以占据总应力的 95%。
2. 热量单一性：海浪由于其运动的准周期性，本身并不直接贡献垂直热通量，热量完全由湍流扩散（Turbulence Diffusivity）承担。

方法论： e-ε 闭合模型夺回“失地”

为了更准确地描述海浪上方的湍流结构，作者摒弃了简单的混合长度模型（Mixing-length Scheme），转而采用两方程模型。

1. 动量平衡方程

$\frac{\partial }{\partial z}(a{u}^t+a{u}^w)=0$ 其中，$a{u}^w$ 根据波谱 $S(k,heta)$ 和增长率参数 $\beta$ 计算。这意味着随着高度增加，波导应力迅速衰减，而湍流应力 $a{u}^t$ 必须迅速增加以保持总应力恒定。

2. TKE 的产生项重构

在波浪边界层（WBL）内，湍流动能 $e$ 的产生项 $P$ 被重新定义为：

• $P^t$：平均风速切变产生的湍流。
• $P^w$：核心创新点，即由波导运动耗散转化而来的湍流动能。

图 1: 在风速 15 m/s 时，波导应力 $a{u}^w$ 在表面占据主导，导致湍流应力 $a{u}^t$ 随高度剧增。

实验与结果：Stanton Number 的真相

模型计算显示，尽管海浪大幅改变了阻力系数，但由于热交换不具备“形阻”机制，它只能通过 $P^w$ 项对湍流扩散率 $K$ 的微弱调制来响应风速变化。

• 阻力系数对比：模型结果与 Anderson (1993) 的回归曲线高度吻合，体现了海浪作为动量汇的显著作用。
• Stanton Number 稳定性：在 3-10 m/s 风速下，$C_H$ 紧贴 $1.0imes10^{-3}$。即使在 20 m/s 强风下，其增幅也仅为 10%，而非传统粗糙度关系预测的 30% 以上。

图 2: 散点为实验观测，实线为本模型。可见模型在巨大的测量噪声中精准捕捉到了 $C_H$ 的“冷淡”趋势。

深度洞察：为什么混合长度理论会失效？

作者在消融分析中发现，如果使用传统的混合长度模型，会低估近海面的湍流扩散。这是因为该模型忽略了波浪能量向湍流的转化（即 $P^w$ 项）。采用 e-ε 模型后，由于海浪诱导的 TKE 在表面层附近增加了 40%，有效增强了热量的混合效率，从而使模型更贴近真实观测。

总结与局限 (Conclusion & Limitations)

本文通过精细的物理建模证明了：海洋动量交换是力学驱动的（含形阻），而感热交换是扩散驱动的（不含形阻）。

局限性：

• 模型尚未考虑 海碎屑 (Sea Spray) 的蒸发效应。在风速 > 15-20 m/s 时，喷溅的海滴可能通过蒸发大幅增加感热和潜热，这是该领域未来的关键课题。
• 假设为稳态 conditions，未考虑波浪破碎（Wave Breaking）的强非线性过程。

启示：这一结论提醒气候学研究者，在开发海气耦合模型时，绝不能简单地通过调整 $z_0$ 来同步缩放 $C_D$ 和 $C_H$。