关于西北太平洋热带气旋最大风速与最低海平面气压关系的研究报告

一、 研究作者、机构及发表信息

本研究的主要作者为 Jae Won Choi（第一作者及通讯作者）、Yumi Cha、Hae-Dong Kim 和 Riyu Lu。他们的所属机构分别是：中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室（Jae Won Choi 和 Riyu Lu）、韩国国立气象科学院（Yumi Cha）以及韩国启明大学全球环境系（Hae-Dong Kim）。该项研究成果以题为《Relationship between the maximum wind speed and the minimum sea level pressure for tropical cyclones in the western north pacific》的学术论文形式，发表于《Journal of Climatology & Weather Forecasting》期刊，并于2016年12月5日在线出版。

二、 研究的学术背景与目的

本研究属于气象学，特别是热带气旋（Tropical Cyclone, TC）强度分析与预报领域。热带气旋的强度通常由两个核心参数定义：近中心最大风速（Vmax）和中心最低海平面气压（Minimum Sea Level Pressure, MSLP）。这两者之间的定量关系（风压关系，Wind-Pressure Relationship）是业务预报和构建历史气候数据集的关键工具。在数据稀疏的洋面上，卫星估算（如Dvorak技术）是评估热带气旋强度的主要手段。然而，Dvorak技术在不同业务中心的应用实践中存在差异：有些中心先用其估算最大风速，再通过风压关系反推最低气压；有些则直接估算气压，再推算风速。因此，一个准确、普适的风压关系对于客观、一致地评估热带气旋强度至关重要。

历史上，西北太平洋地区广泛使用的风压关系是基于简化梯度风平衡（Cyclostrophic Wind Balance）的经验公式，其中最著名的是Atkinson和Holliday于1977年提出的AH1977关系式（Vmax = 3.4(1010 - MSLP)^0.644）。然而，这些传统关系式通常只提供给定气压下的单一风速值，未能考虑影响梯度风平衡的多种物理因素，导致估算存在较大离散性。近期研究（如Knaff和Zehr, 2007；Holland, 2008）指出，热带气旋更符合完整的梯度风平衡（Gradient Wind Balance），并开始将纬度、环境气压、气旋尺寸、移动速度等因素纳入风压模型，但这些研究主要基于大西洋的数据。

基于此背景，本研究旨在利用西北太平洋地区特有的、质量较高的飞机探测历史资料，重新审视该区域热带气旋的风压关系。具体研究目标包括：1）系统性地评估纬度、环境气压、气旋移动速度、强度变化趋势等多个因子对西北太平洋热带气旋风压关系的影响；2）在此基础上，发展一个适用于西北太平洋的、考虑多因子的新型风压关系模型；3）通过与经典模型（如Dvorak表、AH1977、Koba等）的对比，验证新模型的优越性。

三、 详细研究流程与方法

本研究的工作流程严谨，可分为数据准备、因子影响分析、模型构建与验证三大步骤。

第一步：数据准备与样本筛选 研究使用了三类数据：1）1949年至1987年8月15日期间中国气象局（CMA）台风年鉴中的飞机探测资料；2）CMA的最佳路径数据集；3）美国国家环境预报中心（NCEP）的日再分析资料。

飞机探测资料是核心观测数据，每次探测记录包含时间、位置、估计的最低海平面气压（MSLP）、主观估算的地面最大风速以及由多普勒雷达观测的飞行高度层最大风速。为确保数据质量，研究进行了严格筛选： 1. 风速调整：由于飞行高度层风速观测误差（<1.5 m/s）小于机上人员主观估算的地面风速，研究采用Franklin等人（2003）提出的调整因子，将飞行高度层的最大风速外推至10米高度的海面最大风速（Vmax）。这一步提升了风速数据的客观性和准确性。 2. 样本筛选：仅保留调整后的Vmax与飞机主观估算的地面最大风速差异小于10 m/s，且Vmax大于17.2 m/s（热带风暴强度）的样本。经过筛选，最终获得了1766个高质量的飞机探测样本。 3. 辅助参数计算：利用CMA最佳路径数据，通过插值得到与飞机探测时间对应的气旋中心位置和强度，进而计算出气旋的移动速度（Vt）和过去12小时的强度变化趋势（∂Vmax/∂t，代表增强或减弱趋势）。利用NCEP再分析资料，以气旋中心为原点，计算800-1000公里环状区域内的平均海平面气压，作为每个样本的环境气压（Penv）。

第二步：多因子对风压关系的影响分析 此步骤旨在定量和定性地分析各候选因子如何影响MSLP与Vmax的关系。研究方法主要是分组合成分析和多元回归建模。 1. 纬度（φ）效应：将全部样本按气旋中心纬度分为三组：<15°N，15-30°N，>30°N。绘制各组内MSLP与Vmax的散点关系图（图2）。结果显示，在相同Vmax下，纬度越高的热带气旋，其MSLP倾向于越低。这是因为根据梯度风平衡公式，科里奥利力（与纬度正弦值sinφ成正比）随纬度增加而增大，为平衡更大的科里奥利力，需要更强的气压梯度力（对应更低的中心气压）。这种效应在较弱（Vmax < 35 m/s）和较强（Vmax > 50 m/s）的气旋中更为显著。 2. 环境气压（Penv）效应：按环境气压高低将样本分为两组进行合成分析（图6）。结果表明，环境气压是影响风压关系的关键因子。在相同Vmax下，环境气压高的气旋，其中心MSLP也更高。这是因为决定风速的是中心气压与环境气压的差值（气压梯度），而非中心气压的绝对值。 3. 移动速度（Vt）效应：按气旋移动速度分组分析（图4）。发现移动较快的气旋，在相同MSLP下倾向于拥有更大的Vmax。这与先前研究指出的快速移动气旋风速略大的观测事实一致。 4. 强度变化趋势（∂Vmax/∂t）效应：按气旋在过去12小时是增强还是减弱/稳定进行分组（图5）。结果显示，正在增强的气旋，其中心气压相对于其当前风速而言倾向于偏高；而正在减弱的气旋则相反。 5. 量化比较：研究构建了一系列多元线性回归模型，以MSLP（或Penv-MSLP）为因变量，依次加入Vmax、纬度、移动速度、强度趋势、环境气压等自变量，并计算各模型的解释方差（Exp-Var）和平均绝对误差（MAE）。从表1可以看出，仅使用Vmax作为预测因子时，模型解释方差为62.5%，MAE为10.51 hPa。依次加入纬度、移动速度、强度趋势因子后，模型改进有限。然而，当加入环境气压（Penv）因子后，解释方差显著提升至65.8%，MAE降至9.99 hPa，证明了环境气压是改进模型最重要的单一因子。

第三步：新风压关系模型的开发与验证 基于第二步的分析，研究选择对模型改进贡献最大的因子来构建最终的风压关系模型。强度变化趋势因子因其贡献微小（增加解释方差小于0.01%）而被排除。最终构建的多元回归方程以（Penv - MSLP）为因变量，以Vmax、Vmax²、纬度（φ）、sinφ * Vmax（体现科里奥利力效应）以及移动速度（Vt）为自变量。通过回归分析，得到了如下经验公式（方程6）： Penv - MSLP = -33.068 + 0.003Vmax² + 2.058Vmax - 1.677sinφ * Vmax + 61.677φ - 0.119Vt

该模型基于1766个样本数据拟合而成。为评估其性能，研究将其与西北太平洋业务中常用的三种经典风压关系进行了全面对比：1）基于Dvorak（1975）技术表拟合的函数；2）Atkinson和Holliday（1977）的AH1977公式；3）Koba等人（1991）的关系式。对比指标包括解释方差、偏差（Bias）和平均绝对误差（MAE），结果详见表2。

四、 主要研究结果

本研究取得了一系列清晰且具有说服力的结果：

1. 因子影响的确证：合成分析与回归统计明确证实，纬度、环境气压、移动速度和强度变化趋势均对西北太平洋热带气旋的风压关系有系统性影响。给定相同的最大风速（Vmax）时：纬度更高的气旋中心气压更低；处于更高环境气压场中的气旋中心气压更高；移动更快的气旋风速相对其气压更高；正在增强的气旋中心气压相对其当前风速偏高。其中，环境气压和纬度是影响最为显著的两个因子。

2. 新模型的优越性能：新建立的多因子风压关系模型（方程6）在所有评估指标上均显著优于传统的经典模型。

   • 解释方差：新模型解释了67.5%的方差，高于Dvorak（62.1%）、AH1977（62.1%）和Koba（61.2%）模型。
   • 误差：新模型的平均绝对误差（MAE）为9.68 hPa，远低于Dvorak（14.50 hPa）和Koba（12.23 hPa）模型，也优于AH1977模型（11.67 hPa）。
   • 偏差：新模型的偏差（Bias）近乎为零，而其他模型均存在显著的正偏差（即系统性地高估风速或低估气压），例如Dvorak模型的偏差高达12.78 hPa。
   • 物理合理性：最重要的是，新模型不再为每个中心气压提供一个唯一的风速值，而是通过引入纬度、环境气压等变量，能够产生与观测更为吻合的、合理的风速离散分布，这比传统方法提供单一对应值更具物理意义和实用性。

3. 模型对比的量化结论：表2的数据有力地支持了新模型的先进性。它不仅具有更高的统计精度（解释方差最大，误差最小），而且消除了系统偏差。这表明，忽略环境气压和纬度等关键变量的传统模型，在应用于西北太平洋时存在固有的系统性缺陷，而新模型通过纳入这些因子有效地修正了这些缺陷。

五、 研究结论与价值

本研究基于西北太平洋1949-1987年的高质量飞机探测等数据，重新审视并发展了该区域热带气旋的风压关系。主要结论是：西北太平洋热带气旋的风压关系显著受到纬度、环境气压、移动速度及强度变化趋势的影响。其中，环境气压和纬度的影响最为关键。研究成功构建了一个包含Vmax、纬度、移动速度和环境气压因子的新型多元回归风压模型。

该研究的价值体现在以下几个方面： 1. 科学价值：深化了对热带气旋风压关系物理机制的理解，明确了梯度风平衡中科里奥利力项（通过纬度体现）和环境场背景气压在西北太平洋海域的重要性，将风压关系从简单的二元统计推进到考虑多物理影响因子的更精确阶段。 2. 应用价值：所开发的新模型公式简单，易于业务化应用。它可以直接用于改进基于卫星云图（如Dvorak技术）的热带气旋强度估计，无论是通过风速推算气压，还是通过气压推算风速，都能获得比现有方法更准确、偏差更小的结果。这对于提高西北太平洋热带气旋的实时预报精度和历史资料均一性处理具有重要意义。 3. 区域针对性：研究专门针对西北太平洋海域的数据进行建模，弥补了此前类似研究（如Knaff和Zehr， Holland）主要基于大西洋数据的不足，提供了更适用于该区域的本地化风压关系。

六、 研究亮点

1. 数据独特且质量高：研究充分利用了西北太平洋在1949-1987年间宝贵的飞机探测历史资料，并采用客观方法将飞行层风速校正为海面风速，构建了一个数量可观（1766个样本）且质量可靠的数据集，为分析奠定了坚实基础。
2. 因子分析系统全面：研究没有直接套用其他海域的模型，而是系统地、逐一量化分析了多个物理因子（纬度、环境气压、移动速度、变化趋势）对风压关系的影响，并明确了各因子的相对重要性，分析过程细致、逻辑清晰。
3. 模型简洁且改进显著：最终提出的模型在引入关键物理因子的同时保持了相对简洁的线性回归形式，便于业务应用。与广泛使用的经典模型对比，其在解释方差、绝对误差和系统偏差方面均展现出实质性的改进，说服力强。
4. 明确的业务导向：研究的出发点和最终成果都紧密围绕业务预报和气候分析的需求，提出的模型旨在直接改进业务实践中强度估算的准确性，具有很高的实用价值。

七、 其他有价值内容

研究在讨论部分还简要提及，新模型相比传统方法的另一个优势在于能为给定的最低海平面气压提供“一个真实的风速离散分布”，这更符合实际大气中热带气旋结构的多样性。此外，文中指出强度变化趋势（发展率）对模型的改进贡献极小，因此在最终模型中予以剔除，这一决策也体现了研究基于数据驱动的务实态度。