六种非感应起电参数化方案对南京一次雷暴电荷结构的数值模拟分析

利用三维非静力中尺度数值模式Meso-NH中的六种非感应起电参数化方案,采用理想场模拟南京地区雷暴云不同阶段的电荷结构特征,重点分析了雷暴云发展期雪、霰电荷密度对电荷结构的影响,并将成熟期的电荷结构与观测资料比较。结果表明,在雷暴云发展期,电荷分布在上升气流内且高度较高,TAKAH方案为反偶极结构,GARDI方案为三极结构,SAUN、BSMP为反三极结构,SAP98方案为正偶极结构,HELFA方案呈多层结构。在成熟期,电荷分布高度降低,有正负电荷共存区产生。模拟结果中,雪与霰是主要的荷电载体,荷电量及极性与方案有关。雪的电荷密度对云顶的电荷密度影响较大,霰的电荷密度对中下层的电荷密度影响较大。成熟期的电荷结构与VHF宽带干涉仪和雷声联合定位结果对比,SAP98方案模拟结果更符合该次雷暴云电荷结构。     

冰晶核化对雷暴云微物理过程和起电影响的数值模拟研究

本文在已有的二维对流云模式中采用了一种与气溶胶有关的冰晶核化方案替代原有的冰晶核化经验公式, 并选取个例, 分别就两种方案进行了模拟对比试验。模拟结果表明:(1)新方案所得冰晶比含水量主要分布在-0.1~-7.6℃温区之间, 高于原方案所得的-50.1~-24.2℃温区;在整个雷暴云的发展过程中新方案冰晶的分布高度、温度区间以及最大浓度值均大于原方案。(2)在新方案中, 温度相对较高的过冷区产生大量冰晶, 其争食云中水汽抑制了云滴、雨滴的增长。此外, 与原方案相比, 霰增长受雨滴大幅减小的影响进一步得到限制, 导致生成的霰小于原方案, 且空间分布具有较大区别。(3)两方案在雷暴云初期形成的电荷结构不同;在发展旺盛与消散阶段新方案中电荷空间分布区域和电荷量均大于原方案, 此外, 在不同时刻主正电荷区和主负电荷区的中心高度存在差异。本文对云微物理过程及起电的分析为后继探讨气溶胶与雷暴云起放电过程、电荷结构之间的相互关系提供了有利条件。     

混合误差协方差用于集合平方根滤波同化的试验

在集合卡尔曼滤波方法中, 根据预报集合统计提供的依流型而变的预报误差协方差对同化起到决定性的作用。但在集合样本容量不足及模式存在系统误差时, 由预报集合估计的预报误差协方差会出现明显偏差。既要减小这种估计偏差对同化产生的影响而又不增加计算量, 一种可供选择的方法是将定常或准定常的高斯型预报误差协方差和由预报集合估计的预报误差协方差加权平均用于集合卡尔曼滤波同化。利用浅水方程模式, 通过观测系统模拟试验检验在不同的模式误差、 集合成员数以及观测密度条件下, 将这种混合预报误差协方差矩阵用于在集合平方根滤波的效果。试验结果表明, 当预报集合成员数较多而模式又无误差时, 不必采用混合的预报误差协方差矩阵, 否则, 采用混合的预报误差协方差矩阵都有可能改进分析和预报。混合预报误差协方差的最优的权重系数与模式误差关系密切, 模式误差越大, 定常预报误差协方差的权重越大。最优的权重系数与集合成员数及观测密度也有一定关系。     

北京地区一次空气重污染过程的气象条件模拟参数化敏感性试验

气象预报是影响大气重污染预报精度的关键所在。针对2016年12月16～21日北京市一次重污染过程,开展了中尺度气象模式WRF的参数化方案配置敏感性试验。对微物理过程、长波辐射过程、短波辐射过程、陆面过程、边界层过程、近地面过程以及积云对流参数化过程进行组合优选,共设计51组参数化方案组合,分析不同模拟方案下北京市8个气象站点温度、相对湿度、10 m风速的模拟精度及其敏感性。试验结果表明:温度模拟对长波过程参数化方案最为敏感,集合离散度达2.4～7.4°C,再次是短波过程参数化方案;相对湿度模拟也对长波过程参数化方案最敏感,再次是陆面过程;风速模拟对不同过程参数化方案的敏感性程度差异不大。通过模拟结果与观测的统计对比,优选出模拟误差最小的方案组合为Lin微物理方案、RRTMG长波方案、RRTMG短波方案、Tiedtke积云对流方案、Noah陆面方案、MYNN 3rd边界层方案和MYNN近地面方案,并将其与集合平均、基准方案进行对比。对于集合平均来说,其温度模拟与观测相关系数为0.69,高于基准方案,其模拟偏差与均方根误差比基准方案低25%和11%;集合平均的相对湿度和风速模拟相比基准方案变化较小。与集合平均相比,优选方案能同时改进温度、相对湿度和风速模拟,使温度模拟偏差和均方根误差比基准方案下降35%和17%,使相对湿度模拟偏差和均方根误差下降43%和13%,使风速模拟偏差和均方根误差下降33%和24%。以上结果表明,参数化方案的敏感性试验和优选能显著减小重污染期间气象要素的模拟误差,重污染预报改进需重点关注参数化方案模拟上的不确定性。本研究也发现MYNN3rd边界层方案在这次重污染过程的气象要素模拟上具有良好性能,可为未来重污染预报改进提供参考。     

2018年2月琼州海峡持续性海雾过程的数值模拟分析

2018年2月春节期间琼州海峡发生持续性大雾天气,造成大量船舶停航。本文结合葵花8号卫星反演海雾产品、琼州海峡沿岸站点能见度观测数据及美国国家环境预报中心NCEP（National Centers for Environmental Prediction） 提供的FNL（Final Analysis）客观分析资料,对2018年2月18～20日的大雾过程进行了天气学成因分析,并进一步利用CMA-MESO（Global and Regional Assimilation and Prediction System）高分辨率数值模式从边界层方案、模式垂直分层以及海雾能见度算法三个方面进行敏感性试验,以找出模拟效果更好的模式设置方案。研究结果表明:大雾期间华南近海海温较常年平均偏低,受地面冷高压南下补充的弱冷空气影响,偏东暖湿气流流经冷海面并快速凝结。而数值模拟对比试验显示,采用YSU（Yonsei University）边界层方案、边界层垂直层次加密及美国国家海洋大气局预报系统实验室（FSL/NOAA）的海雾诊断方案（简称FSL）对改进能见度预报效果显著:YSU边界层方案比MRF（Medium Range Forecast Model）边界层方案对该次大雾过程的分布范围和最低能见度出现的时间模拟效果更优;模式低层分层加密可更好体现出低能见度的演变过程;通过能见度算法与实况对比,基于模式预报性能较好的湿度和温度预报而来的FSL算法,其能见度预报与站点实况最为接近。     

塔克拉玛干沙漠腹地大气边界层参数化方案的模拟评估

沙尘起沙、沉降、传输均受到沙漠地区大气边界层条件的制约。沙漠地区观测资料匮乏,限制大气边界层模拟效果的检验和评估。利用WRFV3.7.1中尺度数值模式中5种边界层参数化方案(ACM2、BL、MYJ、MYNN2.5、YSU),模拟2014年4月塔克拉玛干沙漠大气边界层特征,并与塔中80 m塔及风廓线雷达晴朗天气下的观测资料对比分析。结果表明:5种方案均能模拟出近地面气温及地表温度,边界层高度,感热、潜热、地表热通量的变化趋势,但未能模拟出边界层风速的日变化趋势,温风湿廓线能较好的反映晴日沙漠地区边界层结构的变化特征,但未模拟出风速随高度变化趋势。沙漠地区下垫面干燥,热容量低,晴天极易形成对流不稳定边界层,非局地湍流参数化方案,ACM2方案是沙漠地区大气边界层模拟较为合理的选择。     

青藏高原夏季土壤有机质及砾石影响水热 传输特性的数值模拟

本文针对青藏高原部分地区土壤有机质和砾石含量较高的特点,在前人工作的基础上,发展了一个新的参数化方案以描述土壤有机质和砾石对土壤导热率、导水率的影响。通过对通用陆面模式CoLM中的土壤水、热参数化方案以及地表蒸发阻抗三方面的逐步改进,对青藏高原藏东南站和纳木错站两种不同下垫面进行单点数值模拟分析。对比原方案与最终优化方案的模拟结果表明:采用新方案的CoLM模式对藏东南站土壤湿度模拟性能明显提高,平均偏差减小到0.04,而对纳木错站浅层20 cm以上土壤湿度的模拟偏差略微增大。新方案在藏东南站对土壤内部温度的模拟改善较为显著,平均偏差减小了0.2℃;而在纳木错站40 cm以上有所改进。新参数化方案较好地模拟了两个观测站表面能量通量的时间变化,纳木错站7、8月份的潜热通量改进尤为明显,比原方案减少大约20 W m-2,与观测结果较为接近。     

积云参数化对0613号台风“珊珊”数值模拟的影响研究

使用W RF的ARW动力核心模拟台风"珊珊",以研究各种积云参数化方案对于台风路径,强度和降水分布等的影响。研究表明,环境场对积云参数化的选取不敏感,积云参数化对台风的强度和路径模拟有明显的影响,KF方案模拟的台风强度较强,增强过程明显,后期移速偏快。KF和BM J对降水的模拟较好,GD方案不够理想。     

AREM模式中引入MYJ边界层方案介绍及其试验

业务实践表明,AREM模式对中国暴雨有较好的模拟能力。当前AREM模式中可供选择的边界层参数化方案比较单一,而不同边界层方案对暴雨的模拟差异较大,所以在AREM中引入一种新的边界层方案是提高其性能的一个重要方面。本文介绍了在AREM模式中引入Mellor-Yamada-Janjic′（MYJ）边界层方案的关键技术与流程,并对MYJ方案的模拟能力进行了评估。模式对2008年四川暴雨的模拟结果表明,MYJ方案有助于改善模式预报的降雨落区位置。