基于多变量灰色模型的台风强度模拟方法

随着台风全路径模拟技术的不断完善,全海域、大范围的台风危险性分析方法得到发展,然而强度模型的模拟精度不能满足需求。针对台风强度演变受多因素影响、动力关系复杂及强度变化随机性强的特点,基于多变量灰色模型对台风强度模拟方法进行研究,并建立起台风多变量灰色强度模型。利用该模型还原了台风“山竹”的强度演变过程,同时对西北太平洋地区整体台风强度进行模拟及检验分析。结果表明:垂直风切变及垂直速度与台风强度的灰色关联度分别为0.771 8、0.745 1,显著高于其他3个环境因素;台风“山竹”强度模拟的后验差比值及小概率误差分别为0.349 1、0.960 8,达到了最高模拟精度;模型模拟所得西北太平洋地区整体的强度变化趋势及分布特征均与地区实际情况保持一致。

位涡源汇和位涡环流及其天气气候意义

在扼要回顾地表位涡研究进展的基础上,本文介绍了复杂地形下的位涡及位涡制造的计算及近年来关于位涡源汇和位涡环流(PVC)的研究进展,侧重介绍青藏高原表层位涡的特殊性及其对天气气候的重要影响。阐明对于绝热和无摩擦大气运动,由于位涡本身的结构重组(位涡重构)可以引起垂直涡度的发展,在夏季可以激发高原涡形成,冬季使青藏高原东部成为重要的表面涡源。基于导得的包括非绝热加热作用的、与等熵面的位移相联系的垂直运动(ωID)方程,进一步阐明青藏高原制造的正位涡沿西风气流东传会引起下游地区低空气旋性涡度、偏南风、和上升运动发展,导致位涡平流随高度增加,激发极端天气气候事件发生发展。指出青藏高原地表加热和云底的潜热释放的日变化显著地影响着地表层位涡的日变化,导致青藏高原的低涡降水系统多在午后至夜间发生发展。证明与传统的青藏高原感热加热指数相比,青藏高原地表层位涡指数能够更好地刻画关于降水的季节变化,与亚洲夏季风降水相关更密切。本文还简单介绍了PVC的概念。指出由于区域边界面的PVC的辐散辐合的变化直接与区域位涡的变化相关联,为保持北半球位涡总量的相对稳定,跨赤道面上的PVC变化与地表PVC的变化必须相互补充,因此跨赤道面上的位涡环流的变化可以成为监测近地表气候变化的窗口。近赤道的海气相互作用能够直接造成沿赤道垂直面上的纬向风垂直切变的变化,激发跨赤道位涡环流异常,从而通过大气内部PVC的变化和青藏高原的调控影响北半球近地表的气候变化。结果表明位涡环流分析为建立热带和热带外大气环流变化的联系开辟了新的蹊径,有着广阔的应用前景。

2024年2月17—23日中国大范围强寒潮雨雪冰冻强对流过程涉及的若干问题

2024年2月中下旬我国出现一次多灾种高影响天气过程。这是一次几十年一遇的过程,出现了大范围强寒潮雨雪冰冻天气并伴随强对流发生,涉及到强寒潮、沙尘、降雨、降雪、冻雨,以及强对流和伴随的大冰雹和雷暴大风,多种高影响天气在一次过程中都有所呈现,其过程之复杂异常罕见。本文针对此次过程中值得深入探讨的问题、高影响天气发生发展可能机理、相应的预报挑战等进行简要梳理,为后续对此次过程的细致和深入分析研究做一个引子。

山西北部极端降水时空演变规律及与大气环流因子的响应

基于山西北部28个国家气象站点,1972—2020年逐日降水资料,选用8个极端降水指数,采用线性回归、Pearson相关分析、连续小波和交叉小波变换分析等方法,研究了山西北部极端降水的时空变化及其与大气环流因子的相关性和周期特征。结果表明:(1)在时间上,山西北部8个极端降水指数都是在20世纪70年代后期和21世纪00年代后期到10年代,年总降水量(PRCPTOT)、中雨以上日数(R10mm)、强降水量(R95P)、极强降水量(R99P)、1 d最大降水量(Rx1day)、5 d最大降水量(Rx5day)均增多,日降水强度(SDII)显著增强,持续湿期日数(CWD)也略有增多。整个20世纪80年代降水异常偏少。(2)在空间上,极端降水指数呈从东北向西南地区逐步增加的态势。从站点趋势变化来分析,大多数站点的极端降水指数呈上升趋势,其中,上升趋势最显著的站点均位于忻州市境内西南部。朔州市境内和忻州市东南部站点PRCPTOT和SDII都呈增加趋势,但CWD则呈减少趋势,由此说明朔州市境内和忻州市东南部地区发生极端降水事件的概率较大。(3)通过小波变换分析发现,1990—2020年山西北部极端降水指数表现出约4 a左右的周期特征。在选取的大气环流指数中,北大西洋涛动指数(NAO)对山西北部极端降水影响最明显,NAO越大,PRCPTOT、R10mm、R95p、R99p、Rx1day、Rx5day、SDII则越小,同时CWD也越少,山西北部产生少雨干旱的几率较大,反之,易发生多雨、洪涝现象。研究结果可为山西北部气象灾害的防御提供科学理论依据。

葵花-8卫星AOD资料在CMA-MESO/CUACE CW 3DVar同化系统中的个例应用研究

CMA-MESO/CUACE CW化学天气耦合模式是自主研发的大气化学耦合模式,目前CMA-MESO/CUACE CW3DVar同化系统实现了地面气溶胶观测可吸入颗粒物PM_(2.5)和PM10的同化,为增强耦合同化系统非常规观测的同化能力,文中在CMA-MESO大气化学天气耦合三维变分同化框架基础上,利用查表法获得气溶胶消光系数,然后建立气溶胶光学厚度(AerosolOpticalDepth,AOD)和气溶胶组分之间关系的观测算子、切线性观测算子和伴随观测算子,实现AOD观测资料的同化应用。针对2016年12月18-20日华北、黄淮地区一次污染天气过程进行同化预报试验,试验结果表明同化葵花-8卫星(Hi‐mawari-8)气溶胶光学厚度观测后,PM_(2.5)分析的重污染区范围有所扩大,山西东南部分析与实况分布更为接近,但是山东大部地区PM_(2.5)分析偏强,与观测相比PM_(2.5)质量浓度存在高估。同时同化Himawari-8 AOD观测和地面气溶胶站点观测的PM_(2.5)分析最优,分析与观测距平相关系数最高,平均偏差、均方根误差及标准差最小。重污染区的PM_(2.5)预报检验结果表明,同化Himawari-8AOD观测对大于350µg·m^(-3)量级PM_(2.5)预报正贡献可以持续到48 h,但整体来说,同时同化Himawari-8AOD观测和地面气溶胶站点观测对各个量级的PM_(2.5)质量浓度预报质量最优。

风压平衡关系对高分辨率区域模式地面资料同化的影响分析

随着模式分辨率的提高,数值模式可分辨的短波信息增加,以地转平衡为基础的风压平衡约束关系不能完全满足低纬度地区及对流尺度高分辨率模式资料同化的需求。提出了基于多重网格方法构建1 km分辨率数值模式变分同化风压平衡约束关系的方案,基于CMA-MESO (China Meteorological Administration Mesoscale Model)三维变分(Three-Dimensional Variational Data Assimilation,简称3DVAR)和1 km分辨率中尺度区域模式开展了单点同化试验和典型个例分析。单点同化试验表明,气压和风的同化结果在没有平衡约束关系条件下比线性平衡约束下更接近观测。地面自动气象站观测资料同化试验表明新风压平衡约束关系将明显增加模式低层几十千米尺度的分析增量,其次是增加100—200 km尺度的分析增量。新风压平衡约束试验的同化分析增量相对于线性平衡约束的对照试验可以给出更多的增量中心位置,更好描述低层水汽辐合,对模式初始0—3 h降水、0—6 h地面10 m风和地面气压的模拟有改进。研究表明在大尺度风压平衡约束的基础上增加局地无风压平衡约束的订正,可以提高对流尺度天气系统地面短波信息的同化能力,改进高分辨率模式短时临近预报。

不同边界层参数化方案在中国西南地区东部高分辨率气候模拟中的敏感性研究

为了进一步评估和提高区域模式对西南地区东部高分辨率气候的模拟能力,利用WRF模式,采用多种边界层参数化方案(下称“不同方案”)对西南地区东部1998—2019年夏季降水和气温进行双重嵌套模拟(外层为D01,内层为D02)。对比不同方案模拟结果表明:多年平均降水量在D01中基本为湿偏差;D02中在四川盆地和重庆低海拔地区为干偏差,湿偏差主要位于贵州和重庆的城口、石柱和武隆一带的地形复杂区;总体上D01中ACM2方案误差最小,D02中MYJ方案误差最小。对多年平均气温的模拟在D01中除了四川盆地一带为暖偏差外其余大部地区基本为冷偏差,D02中大部地区为暖偏差;总体上D01和D02中MYJ方案误差最小,YSU方案最大。对于降水量和平均气温年际变化的模拟技巧在D01和D02中相对较高的地区均集中在重庆中西部和湖北大部地区;降水量总体为YSU方案最高,MYJ方案最低;平均气温总体为MYJ方案最高,ACM2方案最低。因此,提升模式分辨率至对流尺度后对不同气象要素模拟技巧最优的方案存在差异,需根据业务情况选择适合本地的参数化方案。

重庆两场森林火灾蔓延特征的数值模拟

为了深入认识同时段内发生的两场森林火灾蔓延特征以及大气、地形等对不同林火行为的影响,进而为多起林火扑灭火工作提供科学支撑,选取2022年8月17日发生在涪陵两场森林火灾,利用高分辨率地理信息、植被数据和气象数据对两场林火行为进行分析研究,同时运用林火大气耦合WRF-Fire(Weather Research and Forecasting Model with Fire Module)模式进行数值模拟。研究发现:①受深厚高压控制,两场林火是在持续高温、干燥背景下发生的,气象干旱等级达到重旱,火险气象等级高。②模式模拟风向风速变化与实况较为一致,基本体现出火场风向和风速突变,表明WRF-Fire模式可以较准确地再现火场蔓延情况。③基于林火蔓延特征时间和空间变化可将大梁山林火发展分为6个阶段、北山坪林火分为4个阶段,造成直线距离相隔4.4 km两个火场蔓延发展程度不同的原因在于地形差异造成的局地风场差异。

融合DEM与FY-4A数据的ECMWF预报产品深度学习订正方法

精准的数值天气预报是精细化气象公共服务和商业服务的重要前提。欧洲中期天气预报中心(European Center forMedium Weather Forecasting,ECMWF)预报产品在全球被广泛采用,但始终存在系统预报误差。针对数值天气预报中的误差和多源数据融合中的非线性映射等问题,设计了一个ECMWF数值预报产品的深度学习订正模型(Numerical Forecast CorrectionNetwork,NFC-Net)。NFC-Net引入了FY-4A卫星观测数据、数字高程模型数据(Digital Elevation Model,DEM)和ERA5历史实况数据订正预报结果,利用多源数据空间分辨率对齐模块、时空特征提取模块解决多源异构数据特征的提取与融合问题,并通过UNet网络实现ECMWF预报产品的订正。为了评估所提算法的性能,利用NFC-Net对ECMWF产品中的2 m气温和10 m风速两个天气要素开展订正试验,并将试验结果与ECMWF预报结果、ANO方法订正结果、Convlstm方法订正结果、Fuse-CUnet方法订正结果和ERA5实况进行对比。结果显示,NFC-Net模型订正的2 m气温和10 m风速的均方根误差(Root Mean Squared Error,RMSE)分别较ECMWF预报产品下降49.71%和50.86%。表明NFC-Net模型可以充分利用多源数据有效改善复杂地形条件下的订正结果。NFC-Net模型可用于订正ECMWF预报结果,显著提升数值天气预报的精度。

大气热带波动研究进展

本文简要回顾了热带波动的浅水理论解及对应的波动特征。系统介绍了从观测资料中区分和分离热带波动的信号的方法,指出近年来在热带波动观测和分离研究方面取得了明显进展,许多不同于经典波动理论的新特征被揭示。在此基础上回顾了几种影响波动演变的机制(如,波流相互作用、波动—积云对流相互作用、水汽作用等)来部分解释观测热带波动与理论的差异的可能原因。这些研究表明局地大尺度环流和湿对流等在热带波动结构演变、强度和位置中有着重要作用。文中还综述了热带波动对热带气旋生成和路径的影响及可能机制,这表明热带波动的变化能够较好解释热带天气(如,热带气旋)的变化特征。最后,大气热带波动研究依然有许多科学问题亟待解决,它对于大气科学发展具有重要科学意义。

WRF模式不同微物理方案对青海高原夏季一次降水过程模拟差异的初步探讨

为了解不同微物理方案对青海高原地区夏季降水过程模拟的影响,利用WRF模式和NCEP再分析资料,选取Lin方案、WSM6方案和New Thompson方案等3种微物理过程参数化方案,模拟了青海高原地区2017年夏季一次典型降水过程,并结合探空、降水等观测资料对模拟结果进行了探讨。结果表明:WRF模式3种微物理参数化方案对温度廓线的模拟表现出较好的一致性,且均与实况接近,但对不同区域相对湿度廓线的模拟能力有明显差异;3种微物理参数化方案均能模拟出本次降水的主要走向、雨带及强降水中心,但模拟的雨带和强降水中心位置较实况偏东且局部雨量偏大;WSM6方案对降水空间分布的模拟与实况最接近,WSM6和New Thompson方案模拟的降水偏差均略小于Lin方案;3种方案降水模拟结果中,小雨TS评分最高,其次是中雨,对大雨模拟的可信度均较差;降水模拟小雨TS评分最高的是New Thompson方案,为0.60,WSM6方案次之,Lin方案较差,为0.43;对中量降水模拟,WSM6方案略优于其他两个方案;降水模拟New Thompson和WSM6方案整体优于Lin方案。

饱和湿大气中惯性重力内波不稳定判据的探讨

为分析饱和湿斜压大气中任意频率的惯性重力内波的稳定性,利用非静力平衡Boussinesq近似的方程组,在饱和湿斜压大气中考虑湿绝热过程,着重分析讨论任意频率的惯性重力内波不稳定的充分条件(即判据),主要得出:当饱和湿斜压大气的位温经向梯度与饱和相当位温经向梯度之和大于零时,波动不稳定的必要条件是波动向北倾斜;反之,波动不稳定的必要条件是波动向南倾斜。在基流、基本温度场分布不同的饱和湿斜压环境大气中,任意频率惯性重力内波的不稳定判据呈多种不同的表达式,变化复杂。一般情况下,波动不稳定的充分条件即不稳定判据是饱和湿斜压环境大气不仅满足动热力参数Δ大于零,且波动等位相面的倾斜度必须位于与Δ有关的两个参数值A_(1)、A_(2)之间。当位温经向梯度与饱和假相当位温经向梯度的方向相反,波动相对更稳定。如果惯性重力内波发生不稳定,饱和湿斜压环境大气的等饱和相当位温面的倾斜度可能超过等绝对动量面的倾斜度,也可能小于等绝对动量面的倾斜度。

2023年6月一次北部湾低压极端暴雨过程天气气候特征

利用逐日再分析资料,对2023年6月7-10日北部湾低压持续性极端暴雨过程的天气气候特征进行研究。结果表明,(1)持续性极端暴雨前期和期间,副热带高压异常偏北偏西,位置稳定。孟加拉湾-南海-菲律宾一带热带辐合带活跃,南海西部有低压系统生成,并移入北部湾。(2)季节内振荡(MJO)处在第2、3位相,受其东传影响,中高纬度环流形成的槽脊位置有利冷空气南下影响广西。中层500 hPa我国东北地区附近有东北低涡维持,低涡后部强盛偏北风引导冷空气南下,受偏北偏西的副高阻挡,有小股冷空气渗透南下影响广西。(3)低层850 hPa印度洋越赤道气流和西南季风强盛但活动偏南,印度洋-中南半岛-南海的西南风急流受台风“古超”阻挡转向成南风,使气流在北部湾形成气旋性弯曲。广西西北部有弱东北风,增强北部湾周边气流的气旋性辐合。(4)广西南部和北部湾一带位于200 hPa东亚副热带高空急流右侧与850 hPa西南风低空急流左侧,高空辐散和低层辐合,1000~300 hPa水汽通量辐合区水汽条件充沛。中、低纬度天气系统的协同作用,为这次持续性极端暴雨提供有利的稳定大尺度环流背景,以及充分的热力、动力条件。

近40年中国冬季寒潮的气候特征及大气环流异常

利用1980—2023年高分辨率中国气象站观测数据,根据寒潮标准及强度指数,对中国寒潮事件的强度及影响区域进行客观分类,探讨全国性、区域性寒潮事件的时空变化特征和环流演变。结果表明:近40年我国冬季寒潮频次呈显著减少趋势,其影响范围扩大,而强寒潮的强度呈显著增加趋势,且年际变化幅度明显增大。中国寒潮冷空气主要来源于新地岛东南地区,路径因寒潮类型而异。分析全国型、东北华北型和西北华北型寒潮前期和同期大气环流异常特征发现:格陵兰岛的异常深厚暖高压是全国型寒潮的重要前兆,欧亚大陆对流层中高层纬向波列是其爆发的显著特征;东北华北型寒潮与冷涡在中低纬度异常高压系统阻挡下的东移有关;西北华北型寒潮与东欧平原上空暖性高压脊的发展及欧亚大陆两脊一槽的形势密切相关。所有类型寒潮爆发前均有乌拉尔阻塞高压的维持和西伯利亚地区冷空气的堆积。

基于CMA-TRAMS模式地形高度偏差的地面气温误差订正方法研究

采用一元线性方法建立南海台风模式CMA-TRAMS地形高度偏差和地面气温预报误差的回归关系,分别开展不分级、高度偏差分级和地面气温误差分级的三种订正方法的研究,并进行订正效果评估。结果表明,模式地面气温预报误差与地形高度偏差总体呈负的线性相关关系,地面气温预报绝对误差随地形高度偏差绝对值增大而增大(对模式地形高度偏低站点尤为明显),但不同时刻地面气温预报误差特征表现不同,模式对地形高度偏高(即模式地形高于测站高度)和地形高度偏差小于50m的站点,06时地面气温(世界时,下同)预报总体偏低,对地形高度偏低大于50m的站点(即模式地形低于测站高度),06时地面气温预报总体偏高;而无论站点地形高度偏差如何,模式对18时地面气温预报总体偏高。三种订正方法中地面气温误差分级法能有效地减小地面气温预报误差,该方法订正后的分析场准确率可达96%~99%,12~48小时时效预报场准确率总体可提升至90%以上,该方法具有回归关系稳定、效果显著、适用性广、简单易行等特点。

青藏高原东北部地区城市PM_(2.5)和PM_(10)时空分布特征及气象因素的影响

随着工业化进程的迅速发展,大气污染问题已不容忽视,青藏高原东北部作为中国重要的生态安全战略地区,由于其常年多风的地理特点导致土地荒漠化,而土地荒漠化会带来颗粒物污染的问题。笔者以青藏高原东北部地区逐月可吸入颗粒物(PM_(10))和细颗粒物(PM_(2.5))浓度为基础,分析了大气颗粒物PM_(2.5)和PM_(10)时空分布特征,与气象因素(降水量、气温和风速)之间的相关性及受气象因素的影响程度。结果表明:①东部人口密集和经济发达的西宁市、海东市和黄南州PM_(2.5)和PM_(10)较高,以上3个市州的PM_(2.5)平均水平分别为44.2μg/m^(3),44.7μg/m^(3)和36.5μg/m^(3),PM_(10)平均水平分别为99.1μg/m^(3),99.7μg/m^(3)和72.2μg/m^(3);2015~2019年的时间分布上各地区颗粒物浓度呈现逐年下降的趋势;空间分布表明PM_(2.5)呈现从西到东逐渐升高的趋势,PM_(10)则呈东高西低分布。②各地区气温和降水量的峰值均出现在夏季,呈现出“Λ”型的分布规律;而各地区的PM_(2.5)、PM_(10)逐月浓度变化整体呈现出“V”型的分布规律,非采暖季颗粒物浓度最低、采暖季颗粒物浓度最高。③各种气象因素的影响中,PM_(2.5)和PM_(10)与降水量、气温、风速均呈负相关,并且PM_(2.5)浓度受到风速的负向影响,而PM_(10)浓度受到风速的显著正向影响,表明风起扬尘对该区域大气污染贡献突出但风速与污染物浓度的作用机制复杂。本研究可为典型地区空气质量的改善与预测提供理论基础与参考依据。

影响华北重度和轻度高温热浪的大气背景场差异分析

首先识别了1961—2020年华北地区的高温热浪事件,在此基础上筛选出重度和轻度高温热浪事件多发年份,对比分析了重度和轻度高温热浪事件的时空变化特征以及大气背景差异。结果表明,重度和轻度高温热浪发生频次的高值区分别位于华北北部和中南部。自20世纪60年代以来,华北重度和轻度高温热浪事件频次均呈现显著的增长趋势,其中重度高温热浪事件频次的增加更为明显。对应华北重度高温热浪事件多发,东亚中纬度地区为准相当正压结构的异常反气旋性环流控制,华北地区位于异常反气旋性环流的东南部。在该异常环流背景下,华北地区上空呈现异常下沉运动、云量偏少、大气整体偏干。在异常下沉绝热增温的同时,晴朗干燥的配置有利于更多的太阳辐射到达地表,相应的向上长波辐射和感热通量增加,进一步加热近地表大气。在华北轻度高温热浪事件高发年,华北位于异常反气旋(东南侧)和气旋(西北侧)性环流之间,主要受异常西南风控制。异常西南风向华北输送暖平流加热大气的同时,增加了华北上空的水汽。偏多的水汽捕获更多的长波辐射并反射至地表,利于气温升高。

中国气象局乌鲁木齐沙漠气象研究所在次中尺度运动对稳定边界层相似性关系影响的研究中取得进展

次中尺度运动是一类大气现象,空间尺度通常介于几百米到几千米之间,时间尺度介于几分钟到几十分钟之间。常见的大气密度流、孤立波、重力波、低空急流、微冷锋、水平蜿蜒流等都可归为次中尺度运动。稳定大气边界层(SBL)的湍流统计特征会受到非平稳的次中尺度运动的影响,这限制了湍流相似理论在稳定边界层中的应用,导致在强稳定条件下近地层温度、风速等气象要素场的预报存在较大的不确定性。中国气象局乌鲁木齐沙漠气象研究所王敏仲研究员和张海亮副研究员提出了一种利用离散正交小波对稳定边界层次中尺度运动进行识别和滤波的方法,利用2021年6月29日—7月31日塔克拉玛干沙漠夜间高频三维风速和空气虚温的观测数据,分析了滤除次中尺度运动前后湍流相似性关系的变化特征,旨在深入认识次中尺度运动对稳定边界层相似性关系的影响特征和作用机制。

内蒙古巴根黑格其尔铅锌矿地质特征及找矿研究

内蒙古巴根黑格其尔铅锌矿床赋存于上石炭统宝力格庙组安山质凝灰岩中,呈透镜状、脉状。F 3断层是重要的控矿构造,二长花岗斑岩床对成矿控制作用明显,高磁异常对找矿具有直接的指示作用。本文在总结矿床地质特征、矿床成因、成矿规律及控矿因素等基础上,以地层、构造、岩浆岩、成矿空间、磁异常“五位一体”找矿模式选出两处靶区,以期实现找矿突破,为矿区增储上产提供资源基础。

适用于非静力大气模式的近似黎曼求解器应用研究

基于多矩非静力大气模式,开展了3类垂向近似黎曼求解器应用研究。多矩非静力大气模式具有高精度与数值守恒特性,其垂向采用守恒的有限差分格式进行数值离散,而网格单元边界通量计算是通过求解黎曼问题来实现的,因此采用合适的近似黎曼求解器对准确模拟非静力大气垂直运动显得十分关键。LLF(Local Lax-Friedrich)、LMARS(Low Mach Approximate Riemann Solver)和HLLC(Harten-Lax-van Leer Contact)为计算流体力学(CFD)中常用的3种近似黎曼求解器,它们的计算代价和复杂程度逐渐增加。一维标准数值试验表明:LLF计算最为经济,但具有较强的耗散;LMARS具有适用于大气流动的假设,对于数值粘性的控制较好且计算量不大;HLLC建立的三波模型可以避免对中间特征场的过度数值耗散。基于LLF近似黎曼求解器计算经济的特点,通过优化LLF近似黎曼求解器各特征波动的粘性系数,能够实现与LMARS、HLLC近似黎曼求解器相同的性能,且计算代价最小。二维非静力数值试验表明,优化的LLF近似黎曼求解器能够规避常规LLF近似黎曼求解器的数值耗散过大问题,正确模拟小尺度非静力垂直运动,达到更复杂的LMARS、HLLC近似黎曼求解器模拟效果且并未增加计算量,这为非静力大气数值模式提供了良好的参考价值。