第 卷第 期 年 月 ' 杨学胜 沈元芳 徐国强 辐射方案对 (!)*& 全球模式的影响 大气科学 +,-,& +, -., / #$%" 0" % # " (!)*&, 1# 0 辐射方案对 全球模式的影响 杨学胜 沈元芳 徐国强 中国气象科学研究院灾害天气国家重点实验室 数值模式创新基地 北京 摘 要 辐射过程是中期数值天气预报模式和气候模式中一个重要的物理过程 本文针对 (!)*& 1(,! #" 0) *0%" & " # 全球模式非静力半隐式半拉格朗日的特点 进行了不同辐射方案的敏感性比较试验 以期找到一个适合于 (!)*& 全球模式的辐射方案 用到的方案包括 4(5 模式中的 ((/ 长波辐射方案!&! 的 00 0 短波辐射方案以及 * 45 的长短波辐射方案 试验结果表明 * 45 的辐射方案更适合于 (!)*& 模式 天平均的北半球 ) 高度场距平相关系数较 4(5 的辐射方案提高了近 6 另外 针对实际运行中存在的辐射计算不稳定问题 对辐射方案进行了气压订正和温度订正 并对 年 月影响我国的海棠台风进行了模拟试验 气压订正后的 小时模式预报在对副热带高压和台风的描述上更接近于分析 可用预报时效也有了明显提高 关键词 辐射方案 非静力 全球模式 距平相关系数 气压订正温度订正 文章编号 中图分类号 )7 文献标识码!!" "# # $% $&!$ +! -,&8* +,0-9.,!"#$ % &'"&"" ' % ' (0" $ % " # "#$ ""$ %$ # " :" % # " # 0 #, 0 %#" # 0 %," % %" "% 0 ""%0 # ;,,%% " " "3 2$ # " 1 % 0 %" 0 0 " 0# "1 0" % # " (!)*&, 1# 0 / 0 0" % # % 0 ((/, 3 0 00 0 " 3 % # #" 4(5 4 " ( %05 % ",# 0 * 45 0" % # / " ""* 45 0" % # # "1 " " (!)*&, 1# 0 0 0 3, 0 # % " % % " ), $ " ","" " 8 # $ % 1 "6% #$ 0" " 4(5 % #! $ 0" #$ " % %" % # 3 1 % $ " 0" " 0" $ % 0 " % " % #$ "" "1" / % "" " " $ 8",0%" """ 1" $%,"0" $ 0 % 1 01 " % %" % # % " " 0 0 3, 0 # % " % % " ), $ " ","" " 8 # $ #$ 3 0 13 '! 0" $ % 0 ""%, 1# 0 # % " % % "$ 0" #$ " % %" 收稿日期 收修定稿 资助项目 国家科技支撑计划项目!! 计划!!! 作者简介 杨学胜 男 年出生 研究员 主要从事全球数值预报模式研究 * #,2 %#, 3 %
7 卷 ' 引言 辐射过程是中期天气预报模式和气候模式中最重要的物理过程参数化方案之一 最近二十多年来 随着大气遥感探测技术的发展和天气预报业务及气候变化模拟的需要 辐射过程的研究取得了很大的成效 相继开发了多个辐射参数化方案并在预报模式中得到广泛应用 比较著名的包括欧洲中期天气预报中心 * 45* $ " 0 #, 4 " 5 % " 的长短波辐射 % " 5 ; $ %5 0 #% ;1 " 的辐射方案 ;% " 7 &% ": $ " 美国大气研究中心!( " "!"#$ $ %( % 气候模式 ## " #" 0 中的大气模式! ## "!"# $ % 0 的辐射方案 英国气象局 80 气候预报与研究中心的方案 *0 0 " 等 目前 国内数值预报模式和气候模式的物理过程参数化方案 4 " 4, " 7 宇如聪等 7 大多引自国外的模式 如 4(5 4 " ( %05 % ", 模式 % 0 & ## " #" & " # 0 和!()&!03% 0(, ) 0%" & " # 等 而国家气象中心的全球中期预报业务模式 / ; 采用的则是 * 45 的物理过程 陈起英等 7 (!)*& 1(,! #" 0) *0%" & " # 全球模式是近几年发展的一个全球非静力半隐式半拉格朗日中期数值预报模式 庄世宇等 陈德辉等 +, " 采用的是 4(5 的物理过程参数化方案 但 4(5 的物理过程参数化方案主要是针对中尺度数值预报模式 对全球模式是否适合尚不得而知 因此 面对如此众多的参数化方案 我们需要找出一个更适合于自己模式的物理过程参数化方案 同时还需要针对模式运行过程中出现的问题进行优化和改进处理 本文主要对全球模式预报影响较大的辐射过程进行研究 关于辐射过程对模式预报的影响 我国气象学者在这方面也做了大量的研究工作 张华等 沈元芳等 7 王曼等 张凤 张凤等 刘玉芝等 王宏等 石广玉等 其中沈元芳等 年 研究了初始场和长波辐射对气候模拟的影响 其试验表明 不同长波辐射算法模式运行之间的变率比没有长波辐射算法变化的初始场作用大得多 并且 还发现 ## " #" 0 的算法和 * 45 的算法比较接近 而与 *) 的差别较大 但所有结果在全球许多地方和不同尺度上的差别都超过 6 远远大于初始场对气候模拟的作用 国内外对各种辐射方案都做过详细的介绍 这里我们就不再重复 本文的目的 一是找出适合于 (!)*& 全球模式的最佳辐射方案 二是针对 (!)*& 全球模式框架与物理过程的特点 以及在实际运行中出现的问题 对辐射方案中的具体计算细节 如接口 温度的计算等进行优化改进 以提高模式的预报精度 全球模式简介 (!)*& 非静力全球模式动力框架 +, " 使用一套含完全动力学过程的全可压方程组 水平方向上采用经纬度球面坐标 垂直方向采用高度地形追随坐标 积分方案使用两时间层的半隐式半拉格朗日方案 放在极地的水平! 跳点格式 垂直方向为 )$ 变量跳层分布 使用广义共轭余差方法求解关于气压扰动的三维 8 # ": 方程等 其中参考大气廓线考虑了地形的影响 模式预报变量包括扰动气压 扰动位温 水平风 和 垂直速度 和水汽 不同辐射方案的比较试验 (!)*& 全球模式的物理过程参数化方案主要来自于 4(54(5 提供了如下几种辐射方案可供选择 短波方案包括!&! "! " 0&$%!0# " " 的 00 0 " 0 0 等 长波辐射方案包括 ((/ ($0(0" / 0 " 5 ; $ %5 0 # % ;1 " 的长波辐射方案等 关于 4(5 不同辐射方案的比较试验 在 (!)*& 中尺度模式中已做过比较工作 发现
期 杨学胜等 辐射方案对 (!)*& 全球模式的影响 +! -, " / =#$%" (0" &% # " (!)*& 1 0 后的距平相关系数演变情况 试验中 (!)*& 全球模式的分辨率为 垂直方向取 层 时间步长为 可以看出 * 45 的辐射方案的可用预报时效明显优于 4(5 距平相关系数提高了 6 以上 这说明 * 45 的辐射方案更适合于 (!)*& 全球模式 ( 辐射方案的气压与温度订正 图 采用 * 45 辐射方案和 4(5 辐射方案得到的 年 月份 ) 高度场 天平均的距平相关系数演变 5, *3 " 0 3, 0 # % " % % "! )," "" * 45 0" % # * 45 (! 0 " 4(5 0" % # 4(5 (! 4(5 中的几个辐射方案效果基本相当 全球模式中的试验结果也与中尺度模式基本一致 这里就不再赘述 考虑到 4(5 的物理过程主要是针对中尺度模式 而辐射作为全球模式中一个最重要的参数化方案之一 对全球模式是否合适 为此 我们将 * 45 的辐射方案 % " 引入到 (!)*& 全球模式中 并进行了预报试验 为了检验两个辐射方案的预报效果 模式选用 *) 层的分析资料作为初始场进行预报 积分时间为 天 试验时段为 年 月 日至 月 日共计 天 从 月 日起对模式预报与 *) 的分析资料进行逐日的评分检验 比如说 对 月 日的分析进行检验时 对应的预报需要与这天的分析来比较 因此需要用到的预报包括 月 日的 天预报 月 日的 天预报 月 日的 7 天预报 月 日的 天预报 月 日的 天预报 月 日的 天预报 依此类推 最后得到对应于 月 日 日分析的检验 然后对这 天的 7 77 小时的检验评分进行平均生成 天平均后的距平相关系数 试验中 4(5 辐射方案选取 ((/ 长波辐射方案!&! 的 00 0 短波辐射方案 * 45 的辐射方案为 % " 方案 图 给出了 4(5 与 * 45 上述两种不同辐射方案在 年 月份 ) 高度场 天平均 (!)*& 模式在实际运算中 有时在积分至 7 天时会出现辐射过程计算不稳定的现象 这种情况主要发生在 天预报的全球中期预报模式中 而在 天预报以内的有限区模式中则不会出现 这与大气中当气团局地增暖时 空气柱会向上膨胀 单位气压高度差增大 使得该地区上层气压高于四周地区的气压 如下层的气压 而在 * 45 和 4(5 的 ((/ 辐射方案中 正好需要用到相邻两层的气压差的对数 当出现上层气压高于下层气压的情况时 就会出现负的对数 导致计算出现问题 因此 需要对方案进行订正 在实际大气中 除了有强烈上升运动的区域外 空气的垂直运动是较小的 因此在辐射过程中可近似的将大气当作处于静力平衡状态来处理 基于此 模式中对辐射方案采用了如下气压修正方案 假定环境大气满足静力平衡 同时保证气压是随高度单调变化的 ( ) *+ + ' 则各模式层上的气压为 (, *( +',+ ), 其中 ( 为地面气压 (, 为模式整层气压 模式垂直方向为 - 层 由于辐射是在整层上计算 还需要对 式所求的气压进行如下的订正处理 以保证整层的厚度 * ( +(-. ( +(- 其中 (,< 为模式半层气压 订正后的气压为 (,. *(,+' + )/ 7 另外 (!)*& 全球模式采用 )$ 垂直跳层分布 辐射方案中计算通量时需要用到各 陈德辉 薛纪善 杨学胜 等 新一代静力 非静力多尺度通用数值预报动力模式 半隐式半拉格朗日 理论方案设计及试验 技术报告文集 中国气象科学研究院 数值预报中心内部交流材料
卷 ' 图 *) 分析和 1% (!)*& 全球模式 小时预报的 年 月 日 9/) 高度场 1 气压订正方案 % 无订正 5, ), $ " ",",$# # *) 01% % " (!)*&, 1# 0 1 "0% " "" $ % %" % #
期 杨学胜等 辐射方案对 (!)*& 全球模式的影响 +! -, " / =#$%" (0" &% # " (!)*& 1 0 整层边界上的温度 为此 我们采用了一种气压权重方案将整层温度插值到半层上 0,. *0, (, (,. +(,. (,. (,. +(,. 0,. (,. (,. +(, (,. (,. +(, / 为了检验气压和温度修订正后的预报效果 我们进行了预报试验 其中辐射方案选的是 * 45 的方案 图 给出了原方案 即无订正方案 订正方案对 年影响我国的海棠台风 ) 高度场 小时预报及对应的分析场 模式起报时间为 年 月 日 9/ 模式初值为 *) 的 层的分析资料 从图 中可以看出 原方案对副热带高压 线的预报过强 中心强度也过强 台风位于副热带高压内 这对台风路径的预报会产生一定影响 而修正方案在对 线的位置 强度的描述上 尤其是在我国境内与分析很接近 且明显优于原方案 台风的位置也比原方案更接近于分析 图 为原方案和气压订正方案 年 月 天平均的 ) 高度场距平相关系数演变情况 其中模式初始场 *) 的 层分析资料 积分时间为 天 模式试验时段为 年 月 日至 月 日共计 7 天 从 月 日起对模式预报与 *) 的分析资料进行逐日的评分检验 然后再对 天的检验评分进行平均 可以看出 订正方案在预报的中期比原方案有明显的提高 可用预报时效达到 天 需要特别指出的是 虽然在短时间积分的有限 区模式中 上述不稳定情况较少出现 但如果模式采用 4(5 的 ((/ 辐射方案 也应该对气压进行修正 这也是需要使用该方案的研究人员引起特别注意的 ) 总结和讨论 (!)*& 作为新开发的全球非静力半隐式半拉格朗日模式 目前已在中国气象局的 = 巨型机上每天实时运行 结果提供给国家气象中心试用 该模式主要采用 4(5 物理过程参数化方案 但 4(5 的物理过程参数化方案主要是针对中尺度模式 对全球模式是否适合需要通过试验来验证 本文针对 (!)*& 全球模式的特点 对影响中期预报模式和气候模式中的一个重要物理参数化方案 辐射过程进行试验 以期找出一个适合于 (!)*& 全球模式的辐射方案 为此 我们将 * 45 的辐射方案引入到 (!)*& 全球模式中 并进行了 4(5 辐射方案与 * 45 辐射方案的比较试验 结果表明 * 45 的辐射方案更适合于 (!)*& 模式 在北半球 天平均的 ) 高度场距平相关系数提高了近 6 另外 针对 4(5 和 * 45 辐射方案中长时间积分时存在的计算不稳定问题 提出了气压订正和温度订正方案 并对影响我国的海棠台风进行了预报试验 订正后的方案在对副热带高压的描述上和台风中心的位置预报上更接近于分析 北半球 天平均的 ) 高度场距平相关系数在预报的中后期也有一定程度的提高 模式计算稳定性也有了明显提高 从辐射方案的影响试验来看 由于 * 45 的物理过程参数化方案主要针对全球中期预报模式 所以更适合于 (!)*& 全球模式 因此 有必要在今后几年将 * 45 的其他物理方案引入到 (!)*& 全球模式中来 以进一步提高模式的预报性能 参考文献 '# 图 年 月 ) 高度场 天平均的距平相关系数演变 虚线 原方案 实线 气压修正方案 5, *3 " 0 3, 0 # % " % % " )," 陈德辉 沈学顺 新一代数值预报系统 (!)*& 研究进展 应用气象学报 & - ( % "$, (!)*& %0$$%"
卷 '!$$ 0 ",%&% % 陈起英 姚明明 王雨 7 国家气象中心新一代业务中期预报模式 / ; 的主要特点 气象 >,+,#,4,+ 7!, " $ " # 0 #, " % "# 0 / ; " ",% " ",% & :! 0" $ # " :" ;= (! &4 "# $ % " 0 (!&! / % # 7 7 4 / % ## " %#" " # # 0 3 & #" 7 0 # % " 0 % 3 %" 1 3 00," " # 2$ # ", # % " 0# # 0!"# &% 7 *0 0 &,! &" 0 " 21 0" % 0 =,%, ", % # 0 > " ( " & % ;%!!8 * 7!$ # " :" " 1 $ " 0" " " "# $!"# &% 刘玉芝 石广玉 赵剑琦 一维辐射对流模式对云辐射强迫的数值模拟研究 大气科学 7 77 ; + : &,?.! " 0 " 0"3 %, % 0 1, 0# 0"3 % 3 %"3 # 0 7 77 * / 1#& ) " (0"3 " #, "# $ ((/ 30" 0% " 0,# 0 ", 3 $ ( 7 % " =#$%" %, " " 0" " $ # " :" $ % 0 $"%$ $ " " * 45# 0 4 ( 3 7 % " (0" 0% 0 0"3 $ $ " " * 45 % ", " # $ ( &% : $ 5 & / #$ 0 2%, # " 0 3 " 0!%% " $0# " 0 % #$ "" 0 %, " 0 2 $ ( 沈元芳 5 王超 初始场和长波辐射对气候模拟的影响 应用气象学报 7 & +, 54, =#$%" " 00, 3 0 " %#" # 0,!$$ 0 ",% &% % 7 沈元芳 黄丽萍 许国强 等 7 长波辐射对大气变化的敏感性和在 4(5 模式中的模拟检验 气象学报 & +,8,;$,-., " 7 / "3", 3 0" " "# $ %%, 0" # " 4 " ( %05 % "4(5# 0!%" ",%&% 沈元芳 胡江林 (!)*& 模式中的坡地辐射方案及其对短期天气过程模拟的影响 大气科学 & +,8, & $ 0% % # (!)*&0" %" # " ", " $ % 石广玉 王标 张华 等 大气气溶胶的辐射与气候效应 大气科学 7 7 &, 4,?,8 " / 0"3 0%#"% %" " # $ % 7 7 4, 48?? " 7, 0 #% %, 1"# $ %# 0 1 0 # % " # %# " 0 &% % &!7 & $$ 7 王宏 石广玉 王标 等 中国沙漠沙尘气溶胶对沙漠源区及北太平洋地区大气辐射加热的影响 大气科学 4,8,&, 4, " / #$%" 0 " #0 " " 0"3 ", " 3 0 " % 0" ")%%, 王曼 沈元芳 韩慎友 等 两组辐射方案对中国短期天气过程影响的比较和探讨 大气科学 7 7 4, & +,8& " =#$%" ", $ 0" $ # " :" % # # " ", " $ %!"# $ %&% % 7 7 4 -; 8? + "! "# $ %, % % " # 0 0" $ #%!03%!"# $ %&% % +,-, 18, " ) 0 % ":" (!)*&, 1 0 ""% # #$%" # ;,,# 0 &% % & 宇如聪 徐幼平 7!(* 及其对 年汛期降水的模拟 气象学报 7 + ( %,- + $, 7!(* 0" # " " 0 ##!%" " %&% 7 张凤 =!)! 中短波辐射方案的改进研究 = 引入 5 ; 短波辐射方案 气候与环境研究 7?,5, / " 0 " #$ 3 # " " 0" " % # =!)! ) "==% $ ",5 ; % 0 " =!)! #"%0*3 # "( % 7 张凤 曾庆存 =!)! 中短波辐射方案的改进研究 = 短波辐射方案的改进 气候与环境研究?,5,?,>,% / " 0 " #$ 3 # "
期 杨学胜等 辐射方案对 (!)*& 全球模式的影响 +! -, " / =#$%" (0" &% # " (!)*& 1 0 " 0" " % # =!)! ) "=/ #$ 3 # " 0" % # #"%0*3 # "( % 张华 石广玉 一种快速高效的逐线积分大气吸收计算方法 大气科学 7?,8 &,! "0 % " 1 %% " # " 0 " # $ %1 $"!"# $ %&% % 7 庄世宇 薛纪善 朱国富 等 (!)*& 全球三维变分同化系统 基本设计方案与理想试验 大气科学 7?, & -? " (!)*&, 1 ' " # % % # 0,0, 1 3" " " 7 大气和海洋科学快报 征稿启事 中国科学院大气物理研究所主办的英文双月刊 大气和海洋科学快报 英文刊名为!"# $ %0 @% %&% % ; "!@&; 已于 年 月创刊 快报 主要刊登大气科学和海洋物理学方面的创新研究成果 本刊将充分发挥 快报 ; " 类期刊 短平快 的特点 努力把稿件的发表周期限制在半年以内 以保证专家学者的最新研究成果得到及时发表!@&; 处于草创阶段 恳切希望得到广大科学家的支持 欢迎大大家踊跃投稿 由于是 快报 类期刊!@&; 对稿件的长度将有严格的限制 要求稿件不超过 英文单词 包括插图 一幅插图按 7 单词计算 投稿请登录 "$ 0. 2. % % %0 2 $2 注册并提交稿件 联系地址 北京 7 信箱邮政编码 电话 A 传真 A * # # $ % % 4 1 " 0. 2. % % $B % % 大气和海洋科学快报 编辑部