[论文解读] Neutral stability height correction for ocean winds

Ross N. Hoffman|arXiv (Cornell University)|Jul 7, 2011

Ocean Waves and Remote Sensing参考文献 2被引用 3

一句话总结

本文提出了一种简单而精确的解析近似方法，用于在中性稳定度条件下估算海洋表面粗糙度长度（$z_0$），从而实现对海洋风速观测的精确中性稳定度高度校正。该方法通过将 $y = \log(z/z_0)$ 与 $\gamma = (ak^2V^2)/(gz)$ 的关系进行线性拟合，得到 $z_0$ 的闭式表达式，该表达式在风速 1–30 m/s 范围内具有高精度，并可作为牛顿迭代法的优良初始猜测，将高度校正误差降低至百分之几。

ABSTRACT

Adjusting ocean wind observations to a standard height, usually 10 m, requires the use of a boundary layer model, and knowledge of the thermodynamical variables. Height adjustment is complicated by the fact that a necessary parameter, the roughness height, cannot be given in a closed form solution. If only the wind and reporting height are known, the best that can be done is to assume neutral stability. The determination of roughness height is analyzed and a simple approximation (used by Atlas et al. 2011) is derived in detail. This approximation is accurate for winds in the range of 1 - 30 m/s for neutral stratification and would be an excellent initial estimate for a Newton iteration to determine the roughness height precisely, whether or not neutral stability is assumed.

研究动机与目标

• 开发一种在仅知风速和测量高度时，可靠且闭式的海洋表面粗糙度长度（$z_0$）近似方法。
• 通过用与风速相关的 $z_0$ 估计替代恒定 $z_0$ 假设，提高海洋风速观测中性稳定度高度校正的精度。
• 提供一个 $z_0$ 的初始估计，使牛顿-拉夫森迭代法能快速收敛，实现对 $z_0$ 的精确求解，即使在非中性稳定度条件下亦可。
• 通过考虑 $z_0$ 随风速变化的特性，将风速高度校正中的系统性误差从先前的 5–15% 降低至更小。

提出的方法

• 利用查诺克公式和中性阻力系数 $C_{dn} = \left[ \frac{k}{\log(z/z_0)} \right]^2$，推导出关系式 $z_0 = \frac{a}{g} C_{dn} |V|^2$，将 $z_0$ 与风速 $|V|$ 和高度 $z$ 关联起来。
• 引入辅助变量 $\gamma = \frac{ak^2}{gz} |V|^2$，通过 $\gamma = y^2 e^{-y}$ 将 $z_0$ 与可测量量关联，其中 $y = \log(z/z_0)$。
• 采用线性拟合 $y = c_0 + c_1 \log(\gamma)$，其中 $c_0 = 3.7$，$c_1 = -1.165$，适用于 $|V| \leq 40$ m/s 的 $\gamma$ 值，该拟合基于 $z_0$ 值在对数空间 $j = 6$ 到 $30$ 的数值数据推导得出。
• 构建最终近似式 $\hat{z}_0 = z \exp[-(c_0 + c_1 \log(\gamma))]$，实现从 $|V|$ 和 $z$ 直接、非迭代地估算 $z_0$。
• 通过将 $\hat{z}_0$ 与精确值对比验证该方法，结果显示 $z_0$ 的误差极小，对风速比值（如表中所示 $|V_{10}/V_{19.5}|$）的影响可忽略不计（仅相差 0.001）。

实验结果

研究问题

• RQ1 能否推导出一种简单且闭式的海洋表面粗糙度长度（$z_0$）近似方法，使其依赖于中性稳定度下的风速和高度？
• RQ2 与精确的迭代解相比，该近似方法在典型海洋风速（1–30 m/s）范围内的准确性如何？
• RQ3 该近似方法能否作为牛顿-拉夫森迭代法求解 $z_0$ 的有效初始猜测，即使在非中性稳定度条件下也能实现精确求解？
• RQ4 使用与风速相关的 $z_0$ 对 10 米中性风速校正精度的影响如何？与假设 $z_0$ 恒定相比，其改进效果如何？

主要发现

• 所提出的近似式 $\hat{z}_0 = z \exp[-(3.7 - 1.165 \log(\gamma))]$，其中 $\gamma = \frac{ak^2}{gz} |V|^2$，在中性稳定度下对风速 1 至 30 m/s 范围内的 $z_0$ 提供了高度精确的估计。
• 通过考虑 $z_0$ 随风速的变化，该方法将高度校正误差从使用恒定 $z_0$ 时的最高 15% 降低至仅百分之几。
• 该近似方法得到的风速比值（如 $|V_{10}/V_{19.5}|$）与精确值的差异小于 0.001，证实了其在存在微小 $z_0$ 误差时仍具有实际精度。
• 该近似方法足够稳健，可作为牛顿-拉夫森迭代法的优良初始估计，该方法可在 2–4 步内收敛至机器精度。
• 推导出的系数（$c_0 = 3.7$，$c_1 = -1.165$）与 $z$、$a$ 或 $g$ 无关，使得该方法在已知 $|V|$ 时可适用于任意高度，具有普适性。
• 在缺乏热力学稳定度数据的操作场景中，该方法尤为宝贵，因为它可在极少假设下实现精确的中性稳定度校正。