山丘地形的陆面过程及边界层特征的模拟-CAPS.PDFVIP

第 16 卷 1 期 　 应 用 气 象 学 报 Vol. 16 , No. 1 　　 　 　　　 2005 年 2 月 　 JOURNAL OF APPL IED METEOROLO GICAL SCIENCE 　February 2005 山丘地形的陆面过程及边界层特征的模拟 胡小明　刘树华 (北京大学物理学院大气科学系 ,北京 100871) 摘 　　要 将模式 NP89 的陆面过程参数化方法应用到北京大学的三维复杂地形中尺度数值模 式中 ,得到了一个较理想的三维陆面过程及边界层模式 ,利用这个改进的三维模式对 20 km ×20 km 范围的山丘地形的陆面过程及边界层特征进行了数值模拟。模拟结果表明 , 由于 地形阻挡所造成山后的湍流较山前强 ,进而造成近地面温度梯度和感热支出小 ,最终造成山 后的温度比山前的温度明显偏高 ;而且随着山高的增加 ,这种现象更加明显 , 即该模式对山 丘地形条件下的陆面过程和大气边界层特征具有较强的模拟能力 ;模拟结果合理 ,对研究过 山气流形成机制、起伏地形大气边界层物理特征和污染物的扩散具有理论和应用价值。 关键词 : 陆面过程参数化方法 　三维大气边界层模式 　地气相互作用 　过山气流 引 　言 研究复杂地形特别是山区丘陵地形下的陆面过程和大气边界层特征 ,对大气边界层 物理、区域气候和污染物的扩散输送研究具有重要的理论和应用意义。长期以来许多学 者通过不同的方法对复杂地形的陆面过程和大气边界层的特征进行了大量的观测和理论 研究 , 比如 1971 年在美国 Colorado 进行的背风波观测实验[1 ] 。Durran [2 ]建立了两层大气 模式 ,研究了地形起伏对大气波的影响 ,模拟结果表明 ,地形波与整层大气情况下得到的 重力内波形态上完全不同。刘红年等[3 ]采用二阶矩湍流闭合方案 ,分别就非静力和准静 力两种条件建立 PBL 数值模式 , 计算了二维山体对气流平均场和湍流场的影响 ,并分析 [4 ] 了两种模式的计算结果的差别。齐瑛 建立了一个自由大气与对流边界层相耦合的解 析模式 ,利用该模式系统地讨论了对流边界层对过山气流的影响。陈明等[5 ] 根据太行山 区多年的气候统计资料 ,分析了太行山区焚风的气候特征 ,并利用数值模式模拟了地形的 热力强迫效应和动力影响机制对焚风形成和发展的影响。刘辉志等[6 ] 采用 ARPS 4. 0 中 尺度气象模式模拟了对流边界层中气流过山引起的地形波 ,模拟结果表明 ,当大气边界层 是对流边界层时 ,气流过山引起的地形强迫 ,仍能在上部稳定层结中造成足够的垂直扰 动 ,产生向上传播的重力内波。 在对复杂地形的陆面过程和大气边界层特征的数值模拟研究中 , 由于地形起伏,模式 坐标的设计不同于平坦地区的模式 ,在模式中需要做地形坐标变换 ,因此给模式的设计带 ( ) 国家自然科学基金项目 资助。 收到 收到再改稿。 © 1995-2005 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 　　　　　　　　　　　　　　　　　应 　用 　气 　象 　学 　报 　 　 　　　　　　　　　　16 卷 　 14 来一定的困难。本文模式边界层方程组的处理参考了北京大学的三维复杂地形中尺度数 值模式[7 ] ,陆面过程模式参考了 Noilhan 和 Planton 于 1989 年提出的陆面过程的参数化 方法 ,即NP89 模式[8 ] ,将 20 km ×20 km 范围的