边界层参数化方案对台风“烟花”北上阶段模拟影响的差异比较.pdf

第40卷 第4期2023年8月海洋预报MARINE FORECASTSVol.40，No.4Aug.2023收稿日期：2022-04-21。基金项目：天津市自然科学基金（20JCYBJC00640、20JCYBJC00780）;天津市海洋气象重点实验室开放基金项目（2022TKLOM05）;天津市气象局科研项目（202114zdxm01、202115dgxm04）。作者简介：邢蕊（1985-），女，高级工程师，博士，主要从事台风及灾害性天气高分辨率数值模拟研究。E-mail:*通信作者：杨健博（1989-），男，高级工程师，博士，主要从事大气边界层数值模拟研究。E-mail:边界层参数化方案对台风“烟花”北上阶段模拟影响的差异比较邢蕊1,2，杨健博1,3*，庄庭4，王庆元5，邱晓滨1,3，田梦1,3（1.天津市海洋气象重点实验室，天津 300074；2.天津市滨海新区气象局，天津 300457；3.天津市气象科学研究所，天津300074；4.天津市气象探测中心，天津 300061；5.天津市气象台，天津 300074）摘要：利用中尺度数值预报模式（WRF）中的 8 种边界层参数化方案（ACM2、Boulac、GBM、MYJ、MYNN、QNSE、UW、YSU），采用高分辨率数值试验的方法研究不同边界层参数化方案对台风“烟花”（2021）北上影响环渤海区域阶段的路径、强度、降水及动热力结构等方面的模拟差异。结果表明：台风北上后的路径模拟对边界层参数化方案较敏感，其中Boulac方案模拟的路径误差最小,随着积分时间延长，各试验的路径差异也越来越显著；由于台风北上阶段强度较弱，模拟的最低气压绝对误差总体偏小，基本分布在26 hPa之间；各试验对累积降水极值大小和位置的模拟存在较大差异，通过分析24 h累积降水的TS评分可知，Boulac方案在中雨以上量级的模拟中表现最优，而ACM2方案则在大雨及暴雨以上量级的模拟中表现最优；各试验对于台风动热力结构的模拟存在差异，其中动力结构的差异更加明显，并且低层大气比高层大气差异更显著。造成这种差异的原因可能是由于各方案对边界层热通量以及动热力结构的模拟存在差异，并通过边界层顶的夹卷过程将这种差异引入高层大气。关键词：边界层参数化方案；北上台风；数值模拟；环渤海区域中图分类号：P456.7 文献标识码：A文章编号：1003-0239（2023）04-0107-15DOI:10.11737/j.issn.1003-0239.2023.04.0110引言我国是世界上受台风影响最严重的国家之一1。台风登陆过程中带来的强风暴雨常常会诱发严重的自然灾害，给人民的生命和财产安全构成巨大威胁2-3。受三维观测资料的限制，对台风的深入研究主要依靠高分辨率的数值模拟手段。在数值模式中，大气中的许多重要物理过程由于受到模式分辨率等因素的影响，一般通过参数化的方法来进行刻画，如边界层湍流参数化方案、积云对流参数化方案等，不同的参数化方案对模拟结果会产生不同的影响4-7。观测分析表明，台风边界层中存在尺度为几百米几千米的小尺度系统，近地面阵风与台风边界层的精细结构有关8，台风边界层作为台风结构的一个重要组成部分，其动力和热力过程较为复杂，可直接影响到台风的结构和强度9-10。显而易见，边界层参数化方案在台风数值模拟中具有重要作用，也是造成台风强度和结构预报不确定性的来源之一。研究表明，在数值模式中是否考虑边界层参数化方案可显著影响台风的模拟结果10-11，不同的参数化方案由于对边界层中水汽等物理量以及垂直混合等物理过程描述存在差异，也会导致台风路径和强度的模拟结果出现显著差异12-13。前人已针对边界层参数化方案对台风生命史不同阶段的影响开展了研究。例如：温晓培等14研究了边界层参数化方案对台风“SANBA”（2012）初生阶段的影响，结果表明在台风形成初期，不同边海洋预报40卷界层方案对台风的强度和路径模拟具有明显影响，其中QNSE和ACM2方案对该台风发展初期边界层过程的处理较完善，能够较好地反映边界层热带扰动的形成和发展；张福青等11利用 TL6 模式研究了成熟台风边界层的作用，他指出在 24 h 内，台风边界层内各种通量的垂直输送对其路径影响很小而对成熟台风的维持非常重要，其中潜热输送是成熟台风的主要能源，影响最为重要，而感热输送对成熟台风几乎无影响；王叶红等15研究了边界层参数化方案对台风“莫兰蒂”（2016）登陆阶段的模拟影响，定量分析了中尺度数值预报（Weather Research and Forecasting，WRF）模 式WRFv3.8中7种边界层参数化方案在台风登陆减弱阶段的移动路径、强度、结构、降水量以及近地层有关物理量场等方面的差异，并指出不同的边界层参数化方案对于台风路径、强度、降水等方面的模拟具有不同的优势。此外，还有学者针对不同季节、不同影响区域以及不同初始强度的台风展开边界层参数化方案的研究。例如：丁成慧等16利用WRF模式研究了两个 局 地 闭 合（QNSE、MYJ）和 两 个 非 局 地 闭 合（YSU、ACM2）的边界层参数化方案对南海秋季台风“莎莉嘉”（2016）模拟的影响，指出不同边界层方案对于台风路径的影响较小，但对台风强度和结构有明显的影响，非局地闭合方案明显优于局地闭合方案；徐亚钦等17利用WRF模式对强影响浙江的9个西太平洋台风进行多种微物理和边界层参数化方案的模拟，定量评估不同方案下的台风路径、强度和降水的模拟性能，并指出SBU-YLin 微物理和Boulac边界层方案相对更适合浙江台风模拟；王晨稀18以 GRAPES-TCM 为试验模式，分析了两种边界层参数化方案对1109号台风“梅花”预报的影响，并按起报时刻台风强度的不同分为“强台风”和“弱台风”，对比分析了边界层参数化方案对不同强度台风的影响，结果表明YSU方案对弱台风的路径强度预报效果更优，MRF方案对强台风的路径预报效果更好，方案的强度预报效果则与对流参数化方案的选择有关。以上研究表明边界层参数化方案对生命史不同阶段、不同影响地区、不同季节以及不同初始强度的台风模拟均具有重要影响，并且对不同台风的影响存在显著差异，因此针对更多台风以及台风生命史不同阶段开展数值试验来加深边界层参数化方案对台风模拟影响的认识具有重要意义。台风登陆北上后受中纬度斜压性复杂环境及环渤海地区海陆共存下垫面的共同影响，边界层过程较为复杂，模式预报的不确定性较大，风雨精细化预报难度大，因此，准确的边界层垂直湍流输送模拟对台风减弱及陆上结构的预报有重要影响19。然而，目前针对北上影响环渤海区域的台风开展边界层参数化方案数值试验的研究相对较少，不同边界层方案针对北上台风模拟的适用性特征尚不明确。另外，近年来台风北上的频率相对较高，往往会造成暴雨以上的量级，并且有造成严重灾害的可能（如1975年3号台风、1978年5号台风北上后引发的特大暴雨分别给河南和华北东部造成严重灾害），因此针对北上影响环渤海地区的台风开展不同边界层参数化方案影响的研究具有重要意义，可为数值模式选择合适的边界层参数化方案提供参考依据，从而提高北上台风的预报准确率，有助于政府部门开展防灾减灾决策工作。本文使用的资料包括欧洲中期天气预报中心（EuropeanCentreforMedium-RangeWeatherForecasts，ECMWF）的全球气候第五代大气再分析资料 ERA5（Fifth Generation ECMWF AtmosphericReanalysis），分辨率为0.250.25；中央气象台的台风“烟花”实时路径强度资料；中国气象局CMPA实况降水资料，分辨率为0.050.05。利用WRFv4.3模式对台风“烟花”北上影响环渤海阶段进行高分辨率数值模拟，采用8种边界层参数化方案进行数值试验，研究不同方案对台风“烟花”北上阶段路径、强度、降水、动力、热力结构的影响，选取表现最优的北上台风数值模拟边界层参数化方案，为提高北上台风风雨预报准确率及业务应用提供参考。1个例简介2021年第6号台风“烟花”于7月18日生成于西北太平洋洋面，先后以西北、偏西、西北偏北路径逼近大陆沿岸，并于25日04时30分（世界时，下同）在浙江舟山普陀沿海登陆，26日01时50分在浙江嘉108邢蕊等：边界层参数化方案对台风“烟花”北上阶段模拟影响的差异比较4期兴平湖市沿海再次登陆，是有气象记录以来首个在浙江两次登陆的台风，随后台风在陆上向西北以及偏北方向移动，先后进入浙江、江苏、安徽、山东及河北境内，并最终移入渤海湾。台风“烟花”于30日变性为温带气旋，12时停止编号，在中国陆上滞留时间长达95 h，成为1949年以来滞留时间最长的台风20。7月28日00时，台风“烟花”位于500 hPa等压面上低涡的西侧、高空槽的东南方向（图略），之后台风逐渐北移并嵌入到高空槽中，导致高空槽加深，随着台风进一步东移北上，30日00时，低层台风涡旋位于高空槽前，200 hPa有西风带环境位涡侵入到台风环流中（图略），850 hPa相当位温分布图上显示有干冷空气卷入到台风环流中（图略），以上分析表明台风“烟花”在北上过程中逐渐与西风带系统发生相互作用，并最终导致变性。受台风“烟花”影响，海河流域大部分地区出现明显降雨，流域东部普降暴雨大暴雨，其中各有809站（18.1%）出现暴雨和大暴雨，1站出现特大暴雨。7 月 27 日 23 时30 日 23 时最大降雨量为250.7 mm，出现在河北省沧州市，最大小时降雨量为54.3 mm，出现在山西省大同市。另外，经统计分析知，19492019年69月登陆北上后影响环渤海地区的台风共25例，按首次登陆点的不同，可将这25例台风分为山东及山东以北登陆（第一类，11 个）和山东以南登陆（第二类，14个）两类（图略）。第一类台风北上时路径所经下垫面为先海洋后陆地，发生变性的几率较小且北上后无加强，第二类台风路径所经下垫面为先陆地后海洋，移入环渤海区域后发生变性的几率增大，部分个例存在强度加强现象21。本文中研究的台风“烟花”首次登陆点在浙江北部地区，属于第二类，移动路径所经下垫面为先陆地后海洋，在台风移动后期发生了变性，北上过程中强度无加强，但在环渤海地区造成了区域性的暴雨天气，因此对台风“烟花”进行研究可为相似路径台风的风雨预报提供参考依据，具有一定的代表性。2模式简介及数值试验设计利用中尺度数值模式WRFv4.3对台风“烟花”进行三重双向嵌套以及追随台风中心的移动嵌套方案，模拟区域中心位于37N，120E，水平分辨率为 18 km/6 km/2 km，水平方向格点数为 414372、226226、412412，垂直方向共42层，模式顶气压为50 hPa。仅在最外层使用 Kain-Fritsch积云对流参数化方案，云微物理过程选取WSM6方案，短波辐射用Dudhia方案，长波辐射用RRTM方案，陆面过程选取Noah方案。积分时段为2021年7月28日18时30日00时，共30 h，积分时间步长为30 s，积分过程中未使用动力初始化方案及 Spectral Large-scale Nudging方法。模式采用ERA5再分析资料作为初始场和侧边界条件。模 式 分 别 采 用 ACM2、Boulac、GBM、MYJ、MYNN、QNSE、UW、YSU 8个边界层参数化方案对台风“烟花”北上影响环渤海区域阶段的路径强度及风雨过程进行数值试验。由于给定的边界层参数化方案只能与有限的地面层方案配合，为排除地面层参数化方案的影响，在试验中尽量选择较为一致的地面层方案（见表1）。为便于与ERA5资料进行对比，本文使用与ERA5资料分辨率较为接近的D01网格数据进行分析。目前，常用的边界层参数化方案主要可分为3种，即局地闭合、非局地闭合以及（局地-非局地）混合型方案。局地闭合方案假设每个格点上的脉动通量完全由该格点上物理量的平均量决定；非局地闭合方案综合考虑该格点及周边格点对脉动通量的影响；而混合型方案综合考虑了以上两种闭合机制的影响。各边界层参数化方案简要概述如下：YSU方案22：YSU方案是一阶非局地闭合方案。该方案是建立在 K 理论（类比于分子扩散理论，认为湍流通量正比于网格平均量的垂直梯度）的基础上，显示引入了对逆梯度输送过程以及夹卷通量的考虑，以表征大尺度湍涡的非局地输送作用。YSU方案是目前应用最广泛的边界层方案，但一些研究指出，该方案针对夜间边界层的模拟还存在较大偏差，并会高估稳定边界层内垂直混合作用的强度23。ACM2方案24-25：ACM2方案是局地-非局地混合型方案，该方案在显式考虑非局地传输作用的基础上，引入了一个一阶涡旋扩散项，实现了对非局地向上混合以及局地向下混合两种湍流混合机109海洋预报40卷制的综合表征，改进了对近地层气象要素垂直分布的模拟效果。对于稳定或中性边界层，ACM2方案会关闭非局地传输项，仅使用局地闭合机制。ACM2方案的缺点主要在于往往明显高估了边界层高度，这可能与该方案采用了与其他边界层方案均不相同的边界层高度计算方法有关26-27。MYJ方案28-29：MYJ方案是1.5阶局地湍流动能（TKE）闭合方案，认为湍流动量交换系数与TKE的平方根成正比，且边界层高度也由TKE的垂直廓线决定。该方案对于稳定和弱不稳定条件下的边界层都有较好的模拟效果，但在对流边界层中误差较大，容易低估对流边界层内的垂直混合作用强度，导致边界层的偏冷和偏湿。因MYJ方案具有较高的计算效率，是目前应用最广泛的局地闭合型边界层方案。MYNN 方案30：MYNN 方案是在 MYJ 方案的基础上提出的湍流主长度尺度的诊断方程，其加入了水汽凝结过程，为了确保数值积分的稳定性还加入了限定条件以及重新确定了一组更加合理的湍流闭合常数。ROMN-CASCON 等31指出，MYNN方案针对雾过程的模拟效果优于包括MYJ方案在内的其他边界层方案。Boulac方案32：Boulac方案是1.5阶局地闭合方案，该方案考虑了地形对湍流动能的激发作用以及对流边界层（包括弱稳定边界层）内热量通量的逆梯度输送作用。与其他边界层方案相比，Boulac方案对于下垫面信息更加敏感，因此适用于受地形影响较大的地区和事件的模拟。QNSE方案33：QNSE方案是1.5阶局地闭合方案。该方案在不稳定条件下采用TKE闭合方案，而在稳定条件下，采用湍流谱闭合模式发展而来的K-模式。其优势在于保留了更接近真实大气的物理过程，并且考虑了湍流和波的共同作用。徐慧燕等34指出，QNSE方案针对暴雨过程的模拟效果优于其他边界层方案。GBM方案35：GBM方案是1.5阶局地闭合方案，针对边界层顶物理过程的计算采用了夹卷闭合技术。该方案适用于干对流条件下边界层垂直结构的模拟。目前针对GBM边界层方案的应用研究开展的还相对较少。UW 方案36：UW 方案是 1.5 阶局地闭合方案。该方案是在GBM方案基础上，改进了针对气候模式较长时间步长的数值积分稳定性，并提高了计算效率，同时显式地考虑了对流层的夹卷闭合过程，并且针对所有柱状大气中的湍流层采取统一处理。UW方案对于大风条件下的边界层结构的模拟效果较好，因此多应用于风能资源的预测与评估。3各试验台风路径、强度、降水的对比分析3.1台风路径各试验都模拟出了台风北上阶段的路径，模拟的台风中心整体较实况偏东（见图1a）。29日03时前，各个边界层参数化方案的模拟路径差异十分微小，总体路径误差为 3090 km，其中 29 日 0003时误差较显著（见图1a、1b，两个时次各试验的平均表1 数值试验方案设计Tab.1 Designs of numerical experiments试验名称ACM2BoulacGBMMYJMYNNQNSEUWYSU边界层参数化方案ACM2(Pleim)schemeBougeault and Lacarrere(BouLac)TKEGBM TKE-type schemeMellor-Yamada-Janjic(Eta)TKEMYNN 2.5 level TKEQNSE-EDMFBretherton-Park/UW TKE schemeYSU scheme地面层参数化方案Revised MM5 Monin-ObukhovRevised MM5 Monin-ObukhovRevised MM5 Monin-ObukhovMonin-Obukhov(Janjic Eta)Revised MM5 Monin-ObukhovQNSERevised MM5 Monin-ObukhovRevised MM5 Monin-Obukhov110邢蕊等：边界层参数化方案对台风“烟花”北上阶段模拟影响的差异比较4期误差为74.2 km）；29日0612时，模拟与实况路径的误差为3060 km（见图1b），各试验之间的路径差异略有增大；29日12时之后，各试验之间的路径差异显著增大，其中，MYNN 和 GBM 方案在 29日12时30日00时之间的平均路径误差最大（分别为96.6 km和95.5 km），Boulac方案则最小（此时段的平均路径误差为35.5 km），同样，MYNN和GBM方案在29日00时30日00时的24 h平均路径误差依然最大（分别为76.2 km和77.7 km），而Boulac方案最小（47.5 km）。由台风路径误差的标准差时间序列图（见图1c）也可清晰地看出各试验之间路径误差差异的变化趋势：随着台风向北移动，路径误差的标准差整体随时间呈增大趋势，表明随着积分时间的延长，各试验之间的路径差异增大，在29日08时之前，各试验的路径误差标准差小于10 km，之后标准差陡增，在29日17时达到最大，为51.1 km，之后标准差又呈下降趋势，但均在20 km以上,可见随着积分时间延长，边界层参数化方案对台风路径的模拟产生影响。另外，通过引入海陆分布指数（Sea-LandRatio），即台风半径300 km范围内海洋面积所占台风环流总面积的比例37分析也可知，各试验台风环流300 km半径范围内的海洋和陆地下垫面分布比例的差异随积分时间延长逐渐增大，可见随着积分时间的延长，不同方案表现出来的差异不断叠加，导致台风路径的模拟差异逐渐增大。3.2台风强度台风“烟花”在北上影响环渤海区域阶段的强度基本维持在15 m/s（热带低压级），中心最低气压则基本维持在988 hPa左右。由于10 m风场受海陆共存下垫面影响较大，因此从最低气压的角度来分析各试验台风强度的差异。从图2可以看出，8种方案模拟的最低气压所代表的台风强度均弱于实况（中央气象台发布的台风实时路径），但各参数化方案均模拟出了较为平缓的气压变化趋势，除 GBM图1各试验台风移动路径、台风路径误差、海陆分布指数和路径误差标准差分布Fig.1Typhoon track,track error,the sea-land ratio and track error standard deviation111海洋预报40卷方案外，其余7种方案最低气压分布的差异较微小，GBM方案在29日00时30日00时的平均绝对误差最大，达到5.42 hPa，误差最大值则约为6 hPa，其余方案在这一时段的平均绝对误差基本在24 hPa左右，其中ACM2方案最小，为2.42 hPa。由于台风北上阶段强度较弱，模拟的绝对强度误差总体偏小，WRF模式中的边界层参数化方案对于台风北上阶段强度的模拟误差基本为26 hPa。3.3台风降水图 3为实况与各试验模拟的 2021年 7月 29日00时30日00时24 h累积降水分布，为便于比较，此处将实况降水插值到与模拟降水一致的 D01网格上，分辨率为18 km。由图可见，实况24 h累积强降水区呈东北西南带状分布，累积降水极值为165.5 mm，出现在天津西南。就降水落区而言，各试验大于10 mm的降水区分布形态与实况较为一致，均模拟出了东北西南带状的强降水区，但各试验在强降水极值大小与分布位置等方面存在显著差异，其中UW方案的累积降水极值大于实况值，其余方案均小于实况，UW和ACM2方案与实况最接近。图3a主要降水区（黑色方框，37.240.2N，115.5 118 E）平 均 实 况 降 水 量 为 58.1 mm，ACM2、Boulac、GBM、MYJ、MYNN、QNSE、UW 和YSU方案相应的区域平均降水量分别为54.67 mm、51.05 mm、46.30 mm、52.02 mm、52.99 mm、54.88mm、54.98 mm和48.19 mm，结合图中各试验的累积降水分布情况也可见，GBM方案模拟的降水最弱，ACM2、QNSE 和 UW 方案模拟的累积降水相对较强，但各试验模拟的累积降水均弱于实况。将实况降水插值到图3黑色方框区域内D01网格点上，并计算各试验在此黑色方框区域内各量级降水的TS评分，结果见表2。通过分析可知，当累积降水10 mm 时，Boulac 方案的 TS 评分约为0.942，表现最优；当累积降水 25 mm、50 mm、100 mm时，ACM2方案的评分均为最高。因此，对于台风“烟花”的北上阶段，Boulac方案适合于中雨以上量级的模拟，而ACM2方案适合于大雨及暴雨以上量级的模拟。4各试验台风动、热力结构的对比分析为考察各试验模拟的台风动力结构的差异，将2021年7月29日06时和18时的975100 hPa的风速在实况降水极值附近区域（116.2116.6E，37.938.3N）内做平均，得到图 4。由图 4a 可见，29日06时900 hPa以下ERA5资料的区域平均风速随高度的增加逐渐增大，最大可达 19.9 m/s（925 hPa 处）；在 900650 hPa 呈减小趋势，其中825650 hPa 的风速迅速减小，减小幅度为每100 hPa约5.5 m/s；650500 hPa风速基本维持不变；500 hPa以上风速先增大，而后维持不变，最后逐渐减小。各试验模拟的区域平均风速随高度的变化趋势与ERA5资料的基本一致，但彼此之间存在显著差异。在700 hPa以下，ACM2和GBM方案的区域平均风速变化趋势和风速增减幅度与ERA5资料最接近,ACM2方案在875 hPa处风速增加到最大，为21.8 m/s，而 GBM 方案在 875 hPa 处风速增加到最大，为 20.4 m/s，其余方案模拟的风速随高度增加的幅度更大，最大增加幅度为每100 hPa约8 m/s（风速最大约28 m/s）；在700350 hPa，除GBM和MYNN方案的区域平均风速变化趋势与ERA5资料差异较大外，其余各方案之间差别不大，变化趋势与ERA5注：括号中的数字代表29日00时30日00时的平均绝对误差图2观测值与边界层参数化方案模拟的3 h间隔的最低海平面气压Fig.2Observed and simulated minimum sea level pressurewith PBL schemes with a temporal interval of 3 h112邢蕊等：边界层参数化方案对台风“烟花”北上阶段模拟影响的差异比较4期表2 8种边界层参数化方案模拟的各量级24 h累积降水的TS评分Tab.2 Threat scores of 4 magnitudes 24 h accumulated rainfall from eight PBL scheme simulations.降水等级/mmR10R25R50R100ACM20.9290.862*0.765*0.234*Boulac0.942*0.8340.7300.129GBM0.9380.7890.5480.016MYJ0.9250.8590.6920.016MYNN0.9150.8490.7070.054QNSE0.9380.8490.6860.054UW0.9240.8420.7160.083YSU0.8430.8210.7250.104注：*表示方案中最优者注：括号中数值代表24 h累积降水极值，“”代表极值降水出现的位置，黑色方框代表主要降水区图32021年7月29日00时30日00时实况及边界层参数化方案模拟的24 h累积降水量分布Fig.324 h accumulated rainfall observation and simulations with PBL schemes during 0000UTC 29 to 0000UTC 30 July，2021113海洋预报40卷资料也较为类似；300 hPa以上各试验的区域平均风速变化趋势为先增大后减小，ERA5 资料则表现为先维持不变而后减小的趋势；另外，各试验在700 hPa 以下的低层区域平均风速的差异比高层大气更显著。由各试验与ERA5资料所有层次的平均风速差异来看，ACM2方案与ERA5资料的平均风速相差约为1 m/s，两者结果最接近。与 29日 06时相比，各试验对 29日 18时（见图4b）模拟的区域平均风速垂直廓线差异增大，其中，ACM2方案的整体变化趋势与ERA5资料最接近，尤其在825700 hPa，仅有ACM2方案模拟出了风速随高度基本维持不变而后减弱的趋势，其余方案的模拟风速均为随高度的增加先增大后减小，与ERA5资料差异较大。以上分析表明，除边界层外，在数值模拟中不同边界层参数化方案对于边界层以上大气风速的模拟也存在较明显差异，并且随着积分时间的延长，各试验之间的差异越来越显著。ACM2方案总体与 ERA5 资料的风速垂直廓线分布形式较为接近。图5为各试验在29日06时沿模拟的台风中心所在经度的全风速与经向环流剖面。由图可见，各试验在 890 hPa 以下的边界层均为由南向北的流出气流，790 hPa以上则为偏北气流流向台风中心，890790 hPa 之间气流的流入、流出较弱，气流的这种分布形式可能是由于台风中心随高度向东倾斜造成的；各试验在台风中心向北约1经纬度左右的范围内为弱风速区，风速低于 4 m/s；低层大风速区基本分布在 690 hPa 以下，其中 QNSE 方案低层风速极值中心超过 26 m/s（890 hPa 处），在所有方案中最强，除 GBM 方案风速最弱外，其余方案风速极值分布在 2226 m/s；各试验在 590hPa以上存在一个随高度向北倾斜的风速大值带，GBM 方案的风速同样是最弱的。由以上分析可见，各边界层参数化方案所模拟的风场经向剖面虽总体分布形势较为一致，但高层和低层风速的强度均存在显著差异，表明边界层参数化方案不仅对低层大气的风速产生显著影响，对高层大气同样会产生影响。为考察各试验模拟的低层大气辐合、辐散的差异，图6计算了各试验区域平均850 hPa散度随时间的变化，由于台风雨带主要出现在台风中心以北400 km半径范围内（图略），因此选取这一区域范围进行研究。由图6可见，各边界层参数化方案所模拟的区域平均850 hPa散度均为负值，表明在台风中心以北400 km半径范围内存在明显的辐合，各试验的区域平均850 hPa散度存在显著差异，其中GBM方案模拟的低层大气辐合最弱，对应累积降水也最弱。29 日 0012 时（定义为白天时段），ACM2 和UW方案的850 hPa辐合相对较强，29日12时30日00时（定义为夜间时段），QNSE方案的辐合则最强，对应UW、QNSE和ACM2方案模拟的累积降水图4各试验及ERA5资料区域平均风速廓线Fig.4Regional mean profiles of wind speed from ERA5 and simulations114邢蕊等：边界层参数化方案对台风“烟花”北上阶段模拟影响的差异比较4期均较强。各参数化方案对低层大气辐合辐散模拟的差异是造成降水模拟差异的原因之一。为对比各试验对北上台风热力结构模拟的差异，将各试验与ERA5资料的温度和相对湿度在各注：图中横坐标右侧终端点为台风中心所在纬度图52021年7月29日06时各试验沿模拟的台风中心所在经度的全风速（填色，单位：m/s）和经向环流剖面（矢量，垂直速度扩大30倍，单位：m/s）Fig.5The profiles of wind speed(shaded area,unit:m/s)and longitudinal circulation(vector,The vertical speed increases by 30times,unit:m/s)along the longitude of TC center of PBL scheme simulations on 0600 UTC 29 July,2021115海洋预报40卷等压面（116.2116.6E，37.938.3N）区域内做平均，得到温度和相对湿度随高度变化的垂直廓线图。由图 7a 和 7b 可见，各试验与 ERA5 资料的温度垂直廓线总体差异不大，温度误差绝对值分布在02.5 K 之间，温度随高度增加均为递减趋势，700 hPa 高度以下的温度递减率显著小于 700 hPa以上。同样，各试验相对湿度的变化趋势与ERA5资料较为一致，500 hPa以下相对湿度随高度变化不大，基本分布在 90%100%之间（见图 7c、7d），表明北上台风环流的湿层较为深厚。由以上分析可知，各边界层参数化方案对于北上台风热力结构模拟造成的差异小于动力结构的差异。结合图8a温度误差标准差的垂直分布可见，29日06时，随着高度的增加，各试验的温度误差的离散度逐渐减小，各参数化方案对北上台风热力结构模拟的差异在大气低层更显著；29日 18时温度误差的标准差大于06时，表明随着积分时间的延长，各试验表现出来的差异也越来越明显。对于相对湿度，700hPa 以下的低层两个时次的误差标准差均较小，700 hPa以上则逐渐增大（见图8b），表明各试验的相对湿度在高层的差异大于低层，这与图7c、7d的表现较为一致。从各试验850 hPa区域平均相当位温随时间变化的情况看（见图9），除GBM方案显著偏低、随时间呈现逐渐减小的趋势以外，其余各试验 850 hPa区域平均相当位温基本分布在 342.0346.0 K，变化趋势较为类似。29日12时之前（白天时段），各试验区域平均相当位温随时间逐渐增大，差异较小；29日12时之后（夜间时段），相当位温随时间短暂增大后维持基本不变，但各试验之间的差异增大，位温最高的MYNN方案与位温最低的Boulac方案的最大差异达2.2 K。与温度的垂直分布特征类似，随着积分时间的延长，各试验850 hPa区域平均相当位温的差异也越来越大。为分析边界层参数化方案对台风结构模拟造成差异的原因，将各试验在各时次边界层顶的向上垂直质量通量38（VMT=w，VMT代表垂直质量通量，代表空气密度，w代表垂直速度）在（116.2116.6E；37.938.3N）区域内做平均，得到图10。通过分析可知，各试验在边界层顶向上的垂直质量通量具有显著差异，其中QNSE方案的垂直质量通量明显强于其他各方案，表明该方案的垂直混合较强，边界层高度也高于其他各方案（图略），较强的垂直混合能够把边界层内的能量带到高空，从而对台风结构产生影响；相比而言，GBM方案的垂直质量通量最弱，对应边界层高度最低，表明边界层顶的夹卷和垂直混合最弱。从图4中也可看出，QNSE方案在低层的风速强于其他各方案，而GBM方案的风速最弱，表明边界层顶的垂直质量输送和夹卷过程会对台风结构产生影响。另外，感热和潜热的海气交换和动量交换对台风的发生、发展十分重要，感热和潜热输送是边界层的主要物理过程15。通过分析发现，各方案潜热通量的数值显著大于感热通量，表明潜热通量对台风发展的作用更大10。潜热通量和感热通量的分布形式较为类似，海洋上热通量均明显大于陆地，但各方案在量值上具有显著差异（图略），其中QNSE方案的感热和潜热通量明显大于其他方案。有文献指出10，成熟台风的主要能量来源于潜热释放，水汽由边界层辐合上升，释放大量的凝结潜热，为台风的维持和发展提供能量，同时在洋面上通过感热的方式加热台风，为台风继续供给能量。从最低海平面气压的时间序列图中也可以看出（见图2），QNSE方案的台风中心最低气压在所有方案中相对最低，这可能与潜热释放加热大气导致气压降图68种边界层参数化方案模拟的区域平均850 hPa散度时间序列图Fig.6Time series of zone averaged 850 hPa divergence fromeight PBL scheme tests116邢蕊等：边界层参数化方案对台风“烟花”北上阶段模拟影响的差异比较4期图7区域平均温度和相对湿度的垂直廓线分布Fig.7Regional mean profiles of temperature and relative humidity图8各试验与ERA5资料的温度误差标准差和相对湿度误差标准差的垂直廓线图Fig.8Profiles of standard deviation of temperature error and relative humidity error between ERA5 and simulations of eight PBLschemes117海洋预报40卷低有关。综上可见，由于各方案模拟的边界层热通量具有较大差异，对台风发展和结构变化产生影响，通过边界层顶的夹卷和垂直质量输送过程从而对台风整体结构造成影响，各种差异不断累积，最终导致各方案的台风结构产生显著差异。5结论与讨论本 文 利 用 中 尺 度 数 值 模 式 WRFv4.3 中 的ACM2、Boulac、GBM、MYJ、MYNN、QNSE、UW、YSU 8种不同边界层参数化方案，采用高分辨率数值试验的方法研究了边界层参数化方案对台风“烟花”北上移入环渤海区域阶段的模拟影响。结果表明，不同边界层参数化方案对台风路径、降水及动热力结构的模拟差异显著。具体表现为：台风北上后的路径模拟对边界层参数化方案较敏感，其中 Boulac 方案模拟的路径误差最小（24 h平均路径误差为47.5 km），而MYNN和GBM方案模拟的路径误差最大（24 h 平均路径误差为76.2 km和77.7 km）。随着台风向北移动，各边界层参数化方案模拟的路径差异也越来越显著，表现为路径误差标准差随时间积分呈增加趋势。除GBM方案对台风最低气压的模拟显著偏大（强度偏弱）外，其余7种方案模拟的差异较小，其中 ACM2方案在所有试验中的时间平均绝对误差最小，为2.42 hPa。由于台风北上阶段强度较弱，模拟的台风最低气压绝对强度误差总体偏小，基本在26 hPa。各试验对台风降水的模拟存在显著差异，主要体现在强降水极值大小和位置方面，其中UW和ACM2方案的累积降水极值与实况最接近。通过计算24 h累积降水的TS评分可知，Boulac方案在中雨以上量级的模拟中表现最优，而ACM2方案则在大雨及暴雨以上量级的模拟中表现最优。从区域平均风速的垂直廓线及风场经向环流剖面来看，各试验对台风动力结构的模拟具有显著差异，这种差异不仅存在于边界层中，还存在于边界层以上的自由大气中，并且随着积分时间的延长，各试验之间的差异也越来越明显；另外，相比于其他方案，ACM2方案的风速廓线分布变化趋势总体与ERA5资料更接近；各方案对于850 hPa大气辐合辐散的描述也存在较为明显的差异，辐合较强的方案对应模拟的累积降水也较强。通过分析各试验温湿特征量的水平与垂直结构可知，各边界层参数化方案对北上台风热力结图98种边界层参数化方案模拟的区域平均850 hPa相当位温时间序列图Fig.9Time series of zone averaged 850 hPa equivalentpotential temperature from eight PBL scheme tests图108种边界层参数化方案模拟的区域平均边界层顶向上的垂直质量通量时间序列图Fig.10Time series of zone averaged upward vertical massflux on the top of the boundary layer from eight PBL schemetests