Jeffrey D. Kepert就职于澳大利亚气象与气候研究中心（Centre for Australian Weather and Climate Research），其研究《为热带气旋模拟选择边界层参数化方案》（Choosing a Boundary Layer Parameterization for Tropical Cyclone Modeling）发表于《每月天气评论》（Monthly Weather Review）2012年5月号。该文针对热带气旋数值模拟中边界层参数化方案的选择问题，系统评估了多种常用方案的表现，并给出了明确的推荐意见。
从学术背景来看，热带气旋的边界层与其他大气环境中的边界层存在显著差异，这主要是因为气旋自身强烈的旋转对其动力学过程产生了深刻影响。边界层过程直接调控着海气界面的热量、水汽和动量交换，从而在很大程度上决定了风暴的强度；同时，人类与环境所受的气旋影响也主要通过边界层传递。尽管已有大量研究致力于改进极端风速条件下的海表通量参数化，例如发现拖曳系数在约30米/秒以上不再随风速无限增大，但针对边界层上部湍流通量的参数化问题，研究工作相对薄弱。已有若干模拟研究表明，热带气旋的强度和结构对边界层参数化方案高度敏感，不同方案可导致强度差异高达15米/秒、中心气压差异达16百帕。然而，此前的研究虽然指出了部分方案的缺陷，却未能给出明确的方案选择建议。因此，本文的核心目标是通过在简化诊断框架下对代表性方案进行系统比较，揭示各方案差异的物理原因，并为热带气旋模拟用户提供明确的方案选择指导。
在研究方法上，作者使用了Kepert与Wang于2001年提出的诊断性热带气旋边界层模型（KW01模型）作为基础工具，并对其数值方法进行了改进。该模型在给定代表气旋主体的气压场条件下诊断边界层流动，从而排除了边界层对气旋主体的反馈作用，使得不同参数化方案之间的直接差异能够被清晰分离。模型为三维，但本文仅应用于静止风暴，因此实际为二维（半径-高度）问题；垂直方向设有20层，最低层位于10米高度，模式顶为2.25千米，水平网格距为3千米。所有实验均采用相同的海表通量参数化：拖曳系数采用Smith和Montgomery（2008）的公式，热量交换系数取0.0011，以消除海表通量差异带来的干扰。作者在KW01模型中实现了四类代表性边界层参数化方案：一是MM5模式中的Bulk方案和Hi-Res方案，二者均采用常数混合长和基于局地Richardson数的稳定度函数；二是Louis方案，采用Blackadar混合长公式，使近地面混合长趋近于kz；三是非局地K-剖面参数化（KPP）方案，其扩散率廓线为k = ku*z(1 − z/h)^p，其中h为边界层高度、p为形状参数；四是Mellor-Yamada 2.25阶闭合方案，该方案预报湍流动能q并诊断混合长。此外，作者还使用了Louis方案和Mellor-Yamada方案的中性版本，以便分离稳定度效应的影响。为了解各方案在实际应用中的使用频率，作者还统计了2006至2010年间发表于《每月天气评论》、《大气科学杂志》和《皇家气象学会季刊》上使用MM5或WRF进行热带气旋模拟的论文，发现MM5用户中约72%使用Bulk或Hi-Res方案，而WRF用户中超过四分之三使用KPP类方案，其中以YSU方案为主。
模拟结果显示，Bulk和Hi-Res方案存在严重缺陷。这两种方案给出的扩散率在最低模式层达到最大值，导致近地面风速随高度呈近似线性变化而非观测所示的对数分布。观测表明，热带气旋边界层中对数风速廓线应至少延伸至200米高度，而这两种方案在最低数百米内均无法复现这一特征。其根本原因在于采用了不随高度变化的常数混合长，违背了近地面混合长应与距地表距离成正比的理论要求。这一缺陷导致模拟的近地面入流和超梯度风显著偏强，最大入流达15.8米/秒，最大方位角风速为46.3米/秒。由于该方案在黑adar（1976）原文中仅针对夜间边界层设计，且作者明确建议在其他条件下应采用近地面切线形式的混合长，但在MM5实现中这一建议未被采纳，因此作者建议在热带气旋模拟中不要使用Bulk和Hi-Res方案。
Louis方案与Mellor-Yamada高阶闭合方案的表现则较为理想。两者均能正确产生近地面对数风速廓线，入流层深度和超梯度风强度也较为合理。Louis方案的入流和方位角风速最大值均较Bulk方案明显减弱，最大上升流也较弱且出现位置更靠外。扩散率最大值出现在入流层中部而非近地面，垂直分布更符合理论预期。Mellor-Yamada方案的结果与Louis方案在气旋内核区域非常接近，但在外围下沉区，高阶闭合方案对稳定层结导致的湍流抑制更为敏感，入流层更浅且边界层顶更清晰。中性版本的Louis方案与完整版本在上升区差异不大，但在下沉区因忽略稳定层结作用而混合更强、入流层更深。值得注意的是，混合长的渐近值l∞对模拟结果有显著影响：增大l∞会增加中上层扩散率，使近地面方位角风增强、入流减弱、入流层加深。
非局地KPP方案的表现则高度依赖于边界层高度h的诊断。当h取1.5千米时，KPP方案的扩散率远大于Louis和高阶闭合方案，导致入流和超梯度风被过度削弱。增大h会进一步加剧这一问题，这与Braun和Tao（2000）以及Smith和Thomsen（2010）在使用MRF方案时发现的过度混合现象一致。作者认为MRF方案中用于诊断h的体积Richardson数判据在强切变环境中容易给出过高的边界层高度，而YSU方案改用位温廓线判据后有所改善。即使将h设为与其他方案入流深度一致的值，KPP方案在p=2的默认形状参数下仍产生过强的混合；只有将p增大至4并采用随半径变化的h才能接近Louis方案的结果，但这一修改过于剧烈且可能破坏该方案在其他情境下的性能。因此，作者建议对KPP类方案的使用持谨慎态度，至少应检查诊断出的边界层高度是否合理。
在讨论部分，作者通过计算各方案对径向风和方位角风模拟结果在半径-高度平面内的均方根差，进一步量化了方案间的差异。结果表明，所有满足近地面k ≈ ku*z条件的方案（简单混合长方案、Louis方案、Mellor-Yamada方案）彼此差异较小，而违反该条件的Bulk-Hi-Res方案和常数扩散率方案则与其他方案差异显著。作者还发现，在四种方案中，有三种在最大上升流位置附近出现了扩散率向上延伸的瓣状结构，这与眼墙区上升运动一致，但其产生机制并非湍流动能的平流输送，而是局地强切变与上升气流中降低的静力稳定度共同作用的结果。这一发现对Nolan等人关于湍流动能平流在眼墙区作用的解释提出了补充。
文章的结论明确指出：Bulk和Hi-Res方案由于无法复现观测到的近地面对数风速廓线且违反理论预期，不应继续用于热带气旋模拟；KPP类非局地闭合方案对边界层高度诊断敏感，在热带气旋内核中该诊断本身存在困难，因此需谨慎使用；Louis方案和Mellor-Yamada高阶闭合方案在此次评估中未发现重大问题，被推荐为适合热带气旋模拟的方案。作者特别强调，近地面风速廓线的偏差不仅影响边界层本身，还会通过海表通量直接影响风暴强度和发展过程，因此这一问题绝非仅对边界层气象学家有意义。文章还指出，Louis方案与Mellor-Yamada方案在近中性条件下具有理论上的等价性——中性Louis方案可由Mellor-Yamada 2阶方案在局地生成-耗散平衡条件下导出，因此二者不能被视为完全独立的证据。
本研究的学术价值在于，它在一个可控的简化框架内系统分离了边界层参数化方案对热带气旋模拟的直接效应，明确揭示了常用方案中广泛存在但长期未被充分认识的理论缺陷，并为模式用户提供了清晰的选择指南。同时，该研究也提醒学术界，许多已发表的热带气旋动力学研究可能因使用了有缺陷的边界层方案而高估了某些机制的强度，相关结论需要谨慎对待。从应用角度看，这一工作有助于改进热带气旋强度和结构预报，进而提升防灾减灾能力。文章还首次给出了KPP方案在热带气旋条件下混合过强的定量证据，并揭示了扩散率在眼墙区的特殊分布特征，这些发现为后续边界层参数化方案的改进提供了重要参考。