基于ysu的高效空气质量模式在高效空气质量检测中应用.docxVIP

基于ysu的高效空气质量模式在高效空气质量检测中应用 摘要：基于WRF-Chem空气质量模式采用YSU、MYJ、MYNN2、ACM2,4种边界层方案和基于ACM2的敏感性实验方案ACM2R2(模式最低6层的湍流扩散系数阈值改为2m PM国内学者对不同边界层方案在不同地区的污染物浓度模拟表现开展了一系列研究.徐敬等国外,Cuchiara等综上,不同边界层方案在不同地区的表现不一,针对特定地区深入研究不同边界层方案的模拟特性十分必要1 资料与方法选取2016年12月28日～2017年1月5日之间南京冬季一次较为典型污染过程进行研究.期间,南京地区受一弱高压中心影响,风速较小,盛行西南风,出现了多次中度污染天气,且于2016年12月31日中午出现了短暂的重度污染天气,至2017年1月5日空气质量转为良好.PM1.1 观测资料本文使用的观测资料包括地面站点资料、探空廓线资料和ERA5再分析资料.地面站点资料包括时间分辨率为1h的气象要素及PM1.2 WRF-Chem参数设置WRF-Chem模式在美国国家大气研究中心(NCAR)研发的天气预报模式WRF基础上加入了化学模块,并在时间和空间上将化学传输模块与气象模块完全耦合.其不仅可用于天气预报或气候预测,还可模拟颗粒物和臭氧等污染物的浓度、扩散和相互作用等模式主要设置如下:模式模拟的区域如图1所示,空间分辨率为15km,具有119×112个单元网格.其气象边界由1°×1°的NCEP全球再分析资料FNL提供;针对化学边界,本文依据Du等模式中微物理过程参数化方案采用Morrison two-moment本文采用4种边界层方案,即YSU、MYJ、MYNN2和ACM2方案,并增加2个敏感性实验方案ACM2R2和ACM2R5.其中ACM2R2在ACM2方案的基础上,将模式底部6层(约500m)的湍流扩散系数最小阈值设定为2mYSU方案是非局地闭合K理论方案.在湍流局地扩散的基础上增加了对流性大尺度湍涡导致的非局地混合作用,同时考虑边界层顶的夹卷过程.MYJ方案是高阶局地闭合方案,通过湍流动能来闭合湍流方程,求解边界层内的湍流通量.MYJ考虑的物理过程复杂,不能满足假设条件时模拟结果有一定偏差,如在某些区域容易低估大气边界层高度MYNN2方案具有预报次网格尺度湍流动能的优势.MYNN2方案将气压相关项进行了参数化处理,同时考虑浮力的影响,使得其边界层高度模拟具有显著改善ACM2方案特点是通过调节湍流扩散项和非局地项之间的比例系数,实现由稳定条件下的局地扩散算法到不稳定条件下的局地扩散项加非局地扩散项算法的平缓过渡.1.3 模式气象初始场资料及排放源清单模式气象初始场资料采用分辨率为1°×1°的NCEP全球再分析资料FNL,其时间分辨率为6h.模拟时间段为北京时间2016年12月28日00:00～2017年1月5日00:00.模式外部采用准全球模拟的人为排放是从第2版半球污染传输(HTAPv2)1.4 误差分析为验证模式模拟结果的准确性,本文采用平均偏差(MB),均方根误差(RMSE)和相关系数(R)评估模式模拟结果.2 结果与讨论2.1 PM由图2可见,6种方案均较好地模拟出了PM在整个时段中,YSU、MYJ、MYNN2和ACM2,4种方案对于南京地区PM对整个时段而言,YSU、MYJ、MYNN2和ACM2,4种方案中,ACM2的平均偏差MB和均方根误差RMSE最小,而R最大,说明其在此次污染事件的模拟中对于3个统计指标均表现最优.而另外3种方案中,YSU的均方根误差最大,MYJ的平均偏差MB最大,MYNN2的R最小.此外,经过进一步对湍流扩散系数的修订,ACM2R2和ACM2R5在ACM2的基础上进一步提升了这3个统计指标与观测值的接近程度,其中平均误差由24.7μg/m2.2 2m温度与10m风速近地层温度与风速是地-气相互作用中重要的大气下边界信息,其体现了地面过程与大气过程之间相互作用的结果,其中2m温度和10m风速通常作为代表性指标量由表2可知,各方案的平均偏差均为负数,其中MYJ偏差最小为-0.3℃,而经过修订的ACM2R5偏差最大达-1.65℃,这是增加的湍流扩散更利用夜间近地层大气与地表进行热交换,从而使得降温更为剧烈.均方根误差中仍是MYJ最小为1.48℃,ACM2R5最大为1.98℃.因为温度变化趋势的模拟与观测十分吻合,6种方案的R都接近1,MYJ为0.97,其余皆为0.96.由图3(b)可以看出,本次污染过程10m风速较小且呈现明显的日变化特征,白天风速较夜间大,中午风速最大一般超过1m/s但最大不超过2m/s,夜间风速极小,至午夜时分基本接近或低于