2002年北京地区一次强沙尘天气过程的数值模拟研究.pdfVIP

002年北京地区一次强沙尘天气过程的数值模拟研究 程丛兰张小玲王迎春刘伟东徐晓峰 (中国气象局北京城市气象研究所，北京100089) 沙尘天气是影响我国北方地区的一种非常严重的灾害性天气，扬沙和沙尘暴天气 发生时，地表大量的沙粒和尘土被大风吹起进入大气，并随着高空气流向下游作远距 输送，给途径地区的环境和生态造成很大危害。沙尘天气对太阳辐射、能见度及云物 理过程等均有重要影响，注入高层大气的沙尘还进入大气环流过程，所产生的大量悬 浮颗粒在下游地区和整个大气造成了严重的空气污染，近年来，沙尘天气引起了政府 和科学家的广泛关注。由于受观测资料限制，对沙尘暴起动、传输的微过程和分布特 征了解得很少，稀疏的观测数据，甚至卫星得到的整体图象都难以了解其内在过程， 因而用模式模拟的方法来研究和预报沙尘暴发生发展机制和潜在影响在沙尘天气监 测、预警、预报中尤为重要。有关沙尘天气和沙尘暴过程的数值模拟，国内外早已经 开始研究。研究和事实证明沙尘暴天气过程存在大、中、小尺度相互作用的复杂过程， 突发性强，存在很强的斜压不稳定，静力平衡模式已经不太适应。我们在利用高分辨 率的非静力平衡中尺度模式MM5动力框架和详细物理过程，加入了沙尘起沙、传输模 式并与之相耦合，发展了沙尘的起沙、扬沙、输送、沉降模式系统。利用该模式对 2002年3月20日影响北京的强沙尘天气过程进行数值研究。 2002年3月18一--22日全国大部分地区出现了较为严重的沙尘天气过程。从18 日开始，沙尘天气出现在西北地区的甘肃、内蒙古西部、蒙古国西部等地方，随着冷 空气向偏南方向移动，大风经过河套、华北、东北一直向南推进，于3月20日上午 影响到北京。此次强沙尘过程期间北京地区各气象站均观测到了扬沙、浮尘天气，个 别观测站达到了沙尘暴的强度，水平能见度最差时只有200m左右。我们选用19日 12时(世界时)作为初始资料，对此次影响北京沙尘天气过程进行数值研究。 1模式介绍 1．1中尺度动力模式MM5 在进行沙尘天气数值预报和沙尘预报耦合模式的研究时，我们利用PSu／NCAR开 发的最新版本的中尺度模式姗15V3版。动力模式MM5是高分辨率中尺度模拟系统，它 是一个P—o坐标系中的三维流体非静力模式，模式对于中尺度大气系统的发生和发 展至关重要的有关水汽相变，长波短波辐射、行星边界层和陆面过程等物理过程的考 虑，趋于复杂逼真，可供选择的参数化方案较多，是一个比较理想的中尺度动力模式， 352 可以用来进行理论研究和实时预报，对于不同尺度天气现象间的相互作用有较好的分 辨模拟能力，适合用于中尺度天气过程的试验。 1．2沙尘模式 沙尘预报和输送基本方程为： 华p*K—m2af a石 av a仃‘。” 一—∑掣KvK+P+PQ菩+JD吉 0仃 gw一竽娑“一掣婪v 其中：。7‘-一一P P P d％Poy p+=p。一P。，P。和P。分别代表地面气压和模式项层气压，u、v、W分别代表水平风 速u分量、v分量和垂直风速，x、Y、o代表三维坐标(水平坐标和垂直坐标)，g代 表重力加速度，p代表空气密度，Q：代表沙尘分档浓度，瞄代表沙尘沉降速度，D 代表沙尘，K代表分档，m代表地图放大因子，t代表时间，DIV由MM5输出，参看 删5。 本模式选用在天气预报和研究中广泛应用的美国PSU／NCAR发展的非静力中尺度 模式MM5作为动力模式，并发展了沙尘的起沙、扬沙、输送、干沉降、湿沉降模式， 对沙尘粒子谱进行了较细致的分档。模式在处理起沙通量时，引进了风阵性系数，比 较真实地反映了摩擦速度的阵性和地表不均一性特征，对起沙量和起沙临界速度起到 订正作用；模式中湿沉降的考虑，特别按显式和隐式分别考虑，理论上比较客观地模 拟了沙尘的湿沉降过程；沙尘的抬升分为网格尺度和次网格考虑，比较充分地考虑了 沙尘的可能抬升高度；利用遥感资料，对地表进行分类，土壤类型分类，植被覆盖情 况等，分析处理，总结出起沙系数，对起沙区域和起沙量进行订正，比较真实地模拟 出沙尘的分布特征。最后，实现了沙尘起沙传输沉降模式与非静力中尺度模式的完全 耦合，模拟研究影响北京地区沙尘的起沙、输送、沉降的物理过程和物理特征。 2结果分析 2002年3月19～22日北京地区2