边界层参数化方案对台风“烟花”北上阶段模拟影响的差异比较.pdf
《边界层参数化方案对台风“烟花”北上阶段模拟影响的差异比较.pdf》由会员分享,可在线阅读,更多相关《边界层参数化方案对台风“烟花”北上阶段模拟影响的差异比较.pdf(15页珍藏版)》请在咨信网上搜索。
1、第40卷 第4期2023年8月海洋预报MARINE FORECASTSVol.40,No.4Aug.2023收稿日期:2022-04-21。基金项目:天津市自然科学基金(20JCYBJC00640、20JCYBJC00780);天津市海洋气象重点实验室开放基金项目(2022TKLOM05);天津市气象局科研项目(202114zdxm01、202115dgxm04)。作者简介:邢蕊(1985-),女,高级工程师,博士,主要从事台风及灾害性天气高分辨率数值模拟研究。E-mail:*通信作者:杨健博(1989-),男,高级工程师,博士,主要从事大气边界层数值模拟研究。E-mail:边界层参数化方案对
2、台风“烟花”北上阶段模拟影响的差异比较邢蕊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)北上影响环渤海区域阶段的路径、强度、降水及动热力
3、结构等方面的模拟差异。结果表明:台风北上后的路径模拟对边界层参数化方案较敏感,其中Boulac方案模拟的路径误差最小,随着积分时间延长,各试验的路径差异也越来越显著;由于台风北上阶段强度较弱,模拟的最低气压绝对误差总体偏小,基本分布在26 hPa之间;各试验对累积降水极值大小和位置的模拟存在较大差异,通过分析24 h累积降水的TS评分可知,Boulac方案在中雨以上量级的模拟中表现最优,而ACM2方案则在大雨及暴雨以上量级的模拟中表现最优;各试验对于台风动热力结构的模拟存在差异,其中动力结构的差异更加明显,并且低层大气比高层大气差异更显著。造成这种差异的原因可能是由于各方案对边界层热通量以及动
4、热力结构的模拟存在差异,并通过边界层顶的夹卷过程将这种差异引入高层大气。关键词:边界层参数化方案;北上台风;数值模拟;环渤海区域中图分类号:P456.7 文献标识码:A文章编号:1003-0239(2023)04-0107-15DOI:10.11737/j.issn.1003-0239.2023.04.0110引言我国是世界上受台风影响最严重的国家之一1。台风登陆过程中带来的强风暴雨常常会诱发严重的自然灾害,给人民的生命和财产安全构成巨大威胁2-3。受三维观测资料的限制,对台风的深入研究主要依靠高分辨率的数值模拟手段。在数值模式中,大气中的许多重要物理过程由于受到模式分辨率等因素的影响,一般通
5、过参数化的方法来进行刻画,如边界层湍流参数化方案、积云对流参数化方案等,不同的参数化方案对模拟结果会产生不同的影响4-7。观测分析表明,台风边界层中存在尺度为几百米几千米的小尺度系统,近地面阵风与台风边界层的精细结构有关8,台风边界层作为台风结构的一个重要组成部分,其动力和热力过程较为复杂,可直接影响到台风的结构和强度9-10。显而易见,边界层参数化方案在台风数值模拟中具有重要作用,也是造成台风强度和结构预报不确定性的来源之一。研究表明,在数值模式中是否考虑边界层参数化方案可显著影响台风的模拟结果10-11,不同的参数化方案由于对边界层中水汽等物理量以及垂直混合等物理过程描述存在差异,也会导致
6、台风路径和强度的模拟结果出现显著差异12-13。前人已针对边界层参数化方案对台风生命史不同阶段的影响开展了研究。例如:温晓培等14研究了边界层参数化方案对台风“SANBA”(2012)初生阶段的影响,结果表明在台风形成初期,不同边海洋预报40卷界层方案对台风的强度和路径模拟具有明显影响,其中QNSE和ACM2方案对该台风发展初期边界层过程的处理较完善,能够较好地反映边界层热带扰动的形成和发展;张福青等11利用 TL6 模式研究了成熟台风边界层的作用,他指出在 24 h 内,台风边界层内各种通量的垂直输送对其路径影响很小而对成熟台风的维持非常重要,其中潜热输送是成熟台风的主要能源,影响最为重要,
7、而感热输送对成熟台风几乎无影响;王叶红等15研究了边界层参数化方案对台风“莫兰蒂”(2016)登陆阶段的模拟影响,定量分析了中尺度数值预报(Weather Research and Forecasting,WRF)模 式WRFv3.8中7种边界层参数化方案在台风登陆减弱阶段的移动路径、强度、结构、降水量以及近地层有关物理量场等方面的差异,并指出不同的边界层参数化方案对于台风路径、强度、降水等方面的模拟具有不同的优势。此外,还有学者针对不同季节、不同影响区域以及不同初始强度的台风展开边界层参数化方案的研究。例如:丁成慧等16利用WRF模式研究了两个 局 地 闭 合(QNSE、MYJ)和 两 个
8、非 局 地 闭 合(YSU、ACM2)的边界层参数化方案对南海秋季台风“莎莉嘉”(2016)模拟的影响,指出不同边界层方案对于台风路径的影响较小,但对台风强度和结构有明显的影响,非局地闭合方案明显优于局地闭合方案;徐亚钦等17利用WRF模式对强影响浙江的9个西太平洋台风进行多种微物理和边界层参数化方案的模拟,定量评估不同方案下的台风路径、强度和降水的模拟性能,并指出SBU-YLin 微物理和Boulac边界层方案相对更适合浙江台风模拟;王晨稀18以 GRAPES-TCM 为试验模式,分析了两种边界层参数化方案对1109号台风“梅花”预报的影响,并按起报时刻台风强度的不同分为“强台风”和“弱台风
9、”,对比分析了边界层参数化方案对不同强度台风的影响,结果表明YSU方案对弱台风的路径强度预报效果更优,MRF方案对强台风的路径预报效果更好,方案的强度预报效果则与对流参数化方案的选择有关。以上研究表明边界层参数化方案对生命史不同阶段、不同影响地区、不同季节以及不同初始强度的台风模拟均具有重要影响,并且对不同台风的影响存在显著差异,因此针对更多台风以及台风生命史不同阶段开展数值试验来加深边界层参数化方案对台风模拟影响的认识具有重要意义。台风登陆北上后受中纬度斜压性复杂环境及环渤海地区海陆共存下垫面的共同影响,边界层过程较为复杂,模式预报的不确定性较大,风雨精细化预报难度大,因此,准确的边界层垂直
10、湍流输送模拟对台风减弱及陆上结构的预报有重要影响19。然而,目前针对北上影响环渤海区域的台风开展边界层参数化方案数值试验的研究相对较少,不同边界层方案针对北上台风模拟的适用性特征尚不明确。另外,近年来台风北上的频率相对较高,往往会造成暴雨以上的量级,并且有造成严重灾害的可能(如1975年3号台风、1978年5号台风北上后引发的特大暴雨分别给河南和华北东部造成严重灾害),因此针对北上影响环渤海地区的台风开展不同边界层参数化方案影响的研究具有重要意义,可为数值模式选择合适的边界层参数化方案提供参考依据,从而提高北上台风的预报准确率,有助于政府部门开展防灾减灾决策工作。本文使用的资料包括欧洲中期天气
11、预报中心(EuropeanCentreforMedium-RangeWeatherForecasts,ECMWF)的全球气候第五代大气再分析资料 ERA5(Fifth Generation ECMWF AtmosphericReanalysis),分辨率为0.250.25;中央气象台的台风“烟花”实时路径强度资料;中国气象局CMPA实况降水资料,分辨率为0.050.05。利用WRFv4.3模式对台风“烟花”北上影响环渤海阶段进行高分辨率数值模拟,采用8种边界层参数化方案进行数值试验,研究不同方案对台风“烟花”北上阶段路径、强度、降水、动力、热力结构的影响,选取表现最优的北上台风数值模拟边界层参
12、数化方案,为提高北上台风风雨预报准确率及业务应用提供参考。1个例简介2021年第6号台风“烟花”于7月18日生成于西北太平洋洋面,先后以西北、偏西、西北偏北路径逼近大陆沿岸,并于25日04时30分(世界时,下同)在浙江舟山普陀沿海登陆,26日01时50分在浙江嘉108邢蕊等:边界层参数化方案对台风“烟花”北上阶段模拟影响的差异比较4期兴平湖市沿海再次登陆,是有气象记录以来首个在浙江两次登陆的台风,随后台风在陆上向西北以及偏北方向移动,先后进入浙江、江苏、安徽、山东及河北境内,并最终移入渤海湾。台风“烟花”于30日变性为温带气旋,12时停止编号,在中国陆上滞留时间长达95 h,成为1949年以来
13、滞留时间最长的台风20。7月28日00时,台风“烟花”位于500 hPa等压面上低涡的西侧、高空槽的东南方向(图略),之后台风逐渐北移并嵌入到高空槽中,导致高空槽加深,随着台风进一步东移北上,30日00时,低层台风涡旋位于高空槽前,200 hPa有西风带环境位涡侵入到台风环流中(图略),850 hPa相当位温分布图上显示有干冷空气卷入到台风环流中(图略),以上分析表明台风“烟花”在北上过程中逐渐与西风带系统发生相互作用,并最终导致变性。受台风“烟花”影响,海河流域大部分地区出现明显降雨,流域东部普降暴雨大暴雨,其中各有809站(18.1%)出现暴雨和大暴雨,1站出现特大暴雨。7 月 27 日
14、23 时30 日 23 时最大降雨量为250.7 mm,出现在河北省沧州市,最大小时降雨量为54.3 mm,出现在山西省大同市。另外,经统计分析知,19492019年69月登陆北上后影响环渤海地区的台风共25例,按首次登陆点的不同,可将这25例台风分为山东及山东以北登陆(第一类,11 个)和山东以南登陆(第二类,14个)两类(图略)。第一类台风北上时路径所经下垫面为先海洋后陆地,发生变性的几率较小且北上后无加强,第二类台风路径所经下垫面为先陆地后海洋,移入环渤海区域后发生变性的几率增大,部分个例存在强度加强现象21。本文中研究的台风“烟花”首次登陆点在浙江北部地区,属于第二类,移动路径所经下垫
15、面为先陆地后海洋,在台风移动后期发生了变性,北上过程中强度无加强,但在环渤海地区造成了区域性的暴雨天气,因此对台风“烟花”进行研究可为相似路径台风的风雨预报提供参考依据,具有一定的代表性。2模式简介及数值试验设计利用中尺度数值模式WRFv4.3对台风“烟花”进行三重双向嵌套以及追随台风中心的移动嵌套方案,模拟区域中心位于37N,120E,水平分辨率为 18 km/6 km/2 km,水平方向格点数为 414372、226226、412412,垂直方向共42层,模式顶气压为50 hPa。仅在最外层使用 Kain-Fritsch积云对流参数化方案,云微物理过程选取WSM6方案,短波辐射用Dudhi
16、a方案,长波辐射用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个边界层参数化方案对台风“烟花”北上影响环渤海区域阶段的路径强度及风雨过程进行数值试验。由于给定的边界层参数化方案只能与有限的地面层方案配合,为排除地面层参数化方案的影响,在试验中尽量选择较为一致的地面层方案(
17、见表1)。为便于与ERA5资料进行对比,本文使用与ERA5资料分辨率较为接近的D01网格数据进行分析。目前,常用的边界层参数化方案主要可分为3种,即局地闭合、非局地闭合以及(局地-非局地)混合型方案。局地闭合方案假设每个格点上的脉动通量完全由该格点上物理量的平均量决定;非局地闭合方案综合考虑该格点及周边格点对脉动通量的影响;而混合型方案综合考虑了以上两种闭合机制的影响。各边界层参数化方案简要概述如下:YSU方案22:YSU方案是一阶非局地闭合方案。该方案是建立在 K 理论(类比于分子扩散理论,认为湍流通量正比于网格平均量的垂直梯度)的基础上,显示引入了对逆梯度输送过程以及夹卷通量的考虑,以表征
18、大尺度湍涡的非局地输送作用。YSU方案是目前应用最广泛的边界层方案,但一些研究指出,该方案针对夜间边界层的模拟还存在较大偏差,并会高估稳定边界层内垂直混合作用的强度23。ACM2方案24-25:ACM2方案是局地-非局地混合型方案,该方案在显式考虑非局地传输作用的基础上,引入了一个一阶涡旋扩散项,实现了对非局地向上混合以及局地向下混合两种湍流混合机109海洋预报40卷制的综合表征,改进了对近地层气象要素垂直分布的模拟效果。对于稳定或中性边界层,ACM2方案会关闭非局地传输项,仅使用局地闭合机制。ACM2方案的缺点主要在于往往明显高估了边界层高度,这可能与该方案采用了与其他边界层方案均不相同的边
19、界层高度计算方法有关26-27。MYJ方案28-29:MYJ方案是1.5阶局地湍流动能(TKE)闭合方案,认为湍流动量交换系数与TKE的平方根成正比,且边界层高度也由TKE的垂直廓线决定。该方案对于稳定和弱不稳定条件下的边界层都有较好的模拟效果,但在对流边界层中误差较大,容易低估对流边界层内的垂直混合作用强度,导致边界层的偏冷和偏湿。因MYJ方案具有较高的计算效率,是目前应用最广泛的局地闭合型边界层方案。MYNN 方案30:MYNN 方案是在 MYJ 方案的基础上提出的湍流主长度尺度的诊断方程,其加入了水汽凝结过程,为了确保数值积分的稳定性还加入了限定条件以及重新确定了一组更加合理的湍流闭合常
20、数。ROMN-CASCON 等31指出,MYNN方案针对雾过程的模拟效果优于包括MYJ方案在内的其他边界层方案。Boulac方案32:Boulac方案是1.5阶局地闭合方案,该方案考虑了地形对湍流动能的激发作用以及对流边界层(包括弱稳定边界层)内热量通量的逆梯度输送作用。与其他边界层方案相比,Boulac方案对于下垫面信息更加敏感,因此适用于受地形影响较大的地区和事件的模拟。QNSE方案33:QNSE方案是1.5阶局地闭合方案。该方案在不稳定条件下采用TKE闭合方案,而在稳定条件下,采用湍流谱闭合模式发展而来的K-模式。其优势在于保留了更接近真实大气的物理过程,并且考虑了湍流和波的共同作用。徐
21、慧燕等34指出,QNSE方案针对暴雨过程的模拟效果优于其他边界层方案。GBM方案35:GBM方案是1.5阶局地闭合方案,针对边界层顶物理过程的计算采用了夹卷闭合技术。该方案适用于干对流条件下边界层垂直结构的模拟。目前针对GBM边界层方案的应用研究开展的还相对较少。UW 方案36:UW 方案是 1.5 阶局地闭合方案。该方案是在GBM方案基础上,改进了针对气候模式较长时间步长的数值积分稳定性,并提高了计算效率,同时显式地考虑了对流层的夹卷闭合过程,并且针对所有柱状大气中的湍流层采取统一处理。UW方案对于大风条件下的边界层结构的模拟效果较好,因此多应用于风能资源的预测与评估。3各试验台风路径、强度
22、、降水的对比分析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 schemeMel
23、lor-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邢蕊等:边界层参数化方案对台风“烟花”
24、北上阶段模拟影响的差异比较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)也
25、可清晰地看出各试验之间路径误差差异的变化趋势:随着台风向北移动,路径误差的标准差整体随时间呈增大趋势,表明随着积分时间的延长,各试验之间的路径差异增大,在29日08时之前,各试验的路径误差标准差小于10 km,之后标准差陡增,在29日17时达到最大,为51.1 km,之后标准差又呈下降趋势,但均在20 km以上,可见随着积分时间延长,边界层参数化方案对台风路径的模拟产生影响。另外,通过引入海陆分布指数(Sea-LandRatio),即台风半径300 km范围内海洋面积所占台风环流总面积的比例37分析也可知,各试验台风环流300 km半径范围内的海洋和陆地下垫面分布比例的差异随积分时间延长逐渐增
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 边界层 参数 方案 台风 烟花 北上 阶段 模拟 影响 差异 比较
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。