WRF-Chem模式介绍完整版

1、第二章 WRF-Chem 模式介绍WRF-Chemlg式是由美国NOAA预报系统实验室（FSL开发的，气象模式（WRF）和化学模式（Chem）在线完全耦合的新一代的区域空气质量模式。图 给出了 WRF-Chem勺流程框架图。WRF-chemfe含了一种全新的大气化学模式理念。它的化学和气象过程使用相同的水平和垂直坐标系， 相同的物理参数化方案， 不存在时间上的插值， 并且能够考虑化学对气象过程的反馈作用。有别于这之前的大气化学模式，如 SAQM模式、CALGRID1式、MODEL3-CAMQ1式等，它们的气象过程和化学过程是分开的， 一般先运行中尺度气象模式， 得到一定时间间隔的气象场， 然后

2、提供给化学模式使用。这样分开处理以后，存在一些问题：首先，利用这样的气象资料驱动化学过程的时候就存在时间和空间上的插值， 而且丢失了一些小于输出间隔的气象过程，如一次短时间的降水等，而这些过程对化学过程来说可能是很重要的；其次， 气象模式和化学模式使用的物理参数化方案可能是不一样的； 再次， 不能考虑化学过程对气象过程的反馈作用。 事实上， 在实际大气中化学和气象过程是同时发生的， 并且能够互相影响， 如气溶胶能影响地气系统辐射平衡， 气溶胶作为云凝结核，能影响降水，而气温、云和降水对化学过程也有非常强烈的影响。因此，WRF-Chen#g够模拟再现一种更加真实的大气环境。最初版本的WRF-ch

3、emft2002年推出，目前的版本为（2009年4月16日），本文所采用的是WRF-chem 。图 WRF-Chem流程图(来自WRF-Chem V3用户手册)WRF ( Weather Research Forecast , Skamarock et al., 2008)式系统是美国 气象界联合开发的新一代中尺度预报模式和同化系统。WRF 奠式是一个可用来进 行1至10公里内高分辨率模拟的数值模式，同时，也是一个可以做各种不同广泛 应用的数值模式，例如：业务单位正规预报、区域气候模拟、空气质量模拟，理 想个例模拟实验等。故此模式发展的主要目的是改进现有的中尺度数值模式，例如：MM5(NCAR

4、)、ETA(NCEP/NOAA)RUC(FSL/NOAA),希望可以将学术研究以 及业务单位所使用的数值模式整合成单一系统。这个模式采用高度模块化、并行化和分层设计技术，集成了迄今为止在中尺度方面的研究成果。模拟和实时预报 试验表明，WRF真式系统在预报各种天气中都具有较好的性能，具有广阔的应用 前景。化学模式包括了污染物的传输和扩散、干湿沉降、气相化学反应、源排放、光分解、气溶胶动力学和气溶胶化学（包括无机和有机气溶胶）等，并且每一个 过程也都是高度模块化的， 有利于模式的扩展和维护， 也有利于用户选择最合适自己的方案。以下对各部分的处理作一个简单介绍。输送WRF-chem使用的是质量坐标框

5、架，平流输送的处理保持质量和标量的守恒的，空间上采用 5 阶或 6 阶差分，时间上采用 3 阶 Runge-Kutta 分裂显式差分。湍流输送使用阶的 Mellor-Yamada 闭合方案。干沉降WRF-chem中各种痕量气体和气溶胶的干沉降通量的计算使用三层阻力（空气动力学阻抗、 次表层阻抗和表面阻抗） 模型。 表面阻抗的参数化使用了 Wesely（ 1989） 提出的方案。在这种参数化方案中， 表面阻抗主要来自土壤和植被表面，植被特性由使用的下垫面类型资料和季节决定， 表面阻抗也依赖于扩散系数、 活 性气体的可溶性和化学活性。硫酸盐的干沉降使用了不同的方案。 模式中假定硫酸盐都以气溶胶态的

6、形式 存在，干沉降使用了 Erisman等（ 1994）提出的方案。k 阶距多粒径分布气溶胶的干沉降速度v） dk 可以表示为：vdk(radkrardk vGk )vGk()ra 为表面阻抗，v）Gk 为重力沉降速度，rdk 为布朗扩散率。自然源排放WRF-chem 包含了一种完全在线的自然源排放计算方案及两种半在线的方案。其中在线计算方案采用了 Guenther 等（1993, 1994） , Simpson 等（1995）和Schoenemeyer等（1997）的参数化方案，包含了异戊二烯、单菇烯、其他 VOC以及土壤氮的排放。该方案中生成的单据烯和其他VOC按RADM2的物种进行拆分。

7、森林异戊二烯的排放依赖于温度和有效光合辐射通量。 Guenther 等（ 1993）给出了异戊二烯排放的参数化方案， 其排放率正比于标准温度和标准有效光合辐射通量下的排放率。 农田和草地异戊二烯的排放仅仅考虑为温度的函数。 单萜烯、其他VOC和氮的排放也仅仅考虑为温度的函数。其他VOC的排放目前了解的比较少，因此它使用了和单萜烯相同的温度修正系数(Guenther, et al., 1994)。Grell 等 (2000)给出了标准温度和标准光合作用有效辐射通量下的排放率，其中落叶林、针叶林和混合林的排放采用 Guenther 等(1994)，农田和草地的排放采用Schoenemeyer等(1

8、997)。由于使用了 RADM2机制，所有氮的排放都当成是NO,这是一种过高估计，因为忽略了N2O的排放。光分解计算WRF-Chem目前主要包含了两种光分解计算方案：TUV和FASTJTUV (Madronich et al., 1987)光分解速率的计算方案考虑了21 个物种的光分解反应。对应于波长 、物种 i 的光分解系数Ji ，取决于反应物分子的特性，如吸收系数i ( ) 和量子产率i ( ) 。JiIA( , ) i( ) i ( )d ()式中 I A( , ) 为入射通量，它的计算使用了 Wiscombe 的辐射传输模式，考虑了02和03的吸收、雷利散射、气溶胶粒子和云的散射和吸收

9、。吸收系数i()和量子产率 M )的计算采用了 Stockwell等( 1990)给出的方案。FAST-J (Wild, et al., 2000)利用每个尺度段粒子的数浓度、折射指数以及湿半径，根据米散射理论分别计算了在300、 400、 600 和 999 nm 波长下的光学厚度、单次散射反照率0和不对称因子g (气溶胶光学特性的计算参考本文的)，然后再计算气体的光分解速率。 它充分考虑了计算的效率和精度， 可以在三维化学输送模式中使用。 Barnard 等( 2004)对在亚利桑那州、菲尼克斯、得克萨斯州和休斯敦观测的与FAST-H算的NO2光分解速率进行了比较，结果显示不论在晴空还是有

10、云大气情况下， 模拟值都在观测值的误差范围之内， 显示出了很好的一致性。气相化学机制WRF-Chem目前可选择4个光化学机制，它们分别为RADM2 (Stockwell et al.,1990), RACM, CBM4 (Gery et al., 1989和) CBMZ (Zaveri and Peters, 1999。)RADM2机制是第二版本的区域酸性沉降机制，该机制综合考虑了化学过程、预报精度和计算机资源，被广泛应用于各种大气化学模式当中。 RADM2 中无机物包括 14 个稳定物种， 4 个活泼的中间物种和 3 个足量稳定物种(氧气、 氮气和水汽) ；有机物包括26 类稳定物种和 16

11、 类过氧自由基。使用一种准稳态近似(QSSA的方法求解22个诊断物种、38个预报物种的 浓度。38个预报物种的浓度变化速率方程的计算使用后向欧拉差分方案。RACM 机制 0 0 0 0 0 0 0CBM4 机制是一种按结构集成的浓缩机制，该机制包含了 33 个反应物种以及 81 个反应。所有参与反应的物种被分为 4 大类：无机盐物种，明确对待的有机物种， 按碳碳键结构集成的有机物种和按分子量集成的有机物种。 无机盐化学是完全明确对待的， 没有物种的集成。 明确对待的有机物种包含了甲醛， 乙烯和异戊二烯。碳键集成物种描述了不同类型的化学碳键。CBMZ基于CBM4,它包含了 55个预报物种和134

12、个光化学反应。CBMZ在CBM-IV的基础上增加了过氧烷基和过氧酰基的相互作用以及它们与硝基的反应； 增加了长寿命的有机硝酸盐和过氧化氢物的反应机制； 明确对待了不活泼的烷烃如甲烷、乙烷；修订了无机化学、异戊二烯化学以及活泼的烷烃、烯烃和芳香烂的化学反应；耦合了海洋二甲基硫(DMS用排放及其化学反应。CBMZ使用一种体制依赖方法， 该方法把化学动力学机制区分为背景的、 人为源的和自然源的三种亚机制，这样可以全面缩减计算时间。Jiang等(2002)和Fast等(2002)分别从城市和区域尺度对CBMZ 的空气质量模拟效果作了评估。MADE/SORGAM 和气溶胶参数化WRF-Chem目前包含了

13、两个气溶胶计算方案，分别为MOSAIC (Zaveri, et al., 2005a,。MADE/SORGAM欧洲气溶胶动力学模式 MADE是由区域颗粒物模式PPM发展而来的，能够 提供详细的关于粒子化学组成、尺度分布以及影响粒子数浓度的动力学过程的信 息。早期版本的MADE仅限于亚微米量级的无机盐和水组成的气溶胶，进一步改 进引入了 MODELS/CMAQ奠式系统的气溶胶部分，使彳# MADE包含了粗模态粒子 及更详细的细颗粒化学组成的描述。Schell等（2001）把二次有机气溶胶（SOA） 加入到了 MADE中，即二次有机气溶胶模式（SORGAM。亚微米气溶胶的尺度分布分为两个相互重叠的

14、模态，各个模态满足对数正态分布。“、 N rn(ln dp)exp2 ln g1 (ln dpln dpg)2LAT2 ln2 g()式中N为粒子的数浓度，d p为粒子直径,dpg为中值直径，g为尺度分布的标()()二次气溶胶最重要的形成过程是在硫酸-水系统中的均相核化，计算采用准偏差。尺度分布的k阶距Mk定义为:Mkdpn(lndp)d(ln dp) 将（）式代入（）式可得卜2Mk Ndkgexpk2ln2 gKulmala等（1998）给定的方法。粒子的凝结增长分为两步，一是化学反应产生可 凝结的蒸汽，二是挥发性物种在气溶胶表面的凝结和蒸发。 在MADE中忽略了开 尔文效应，对于连续的和自

15、由分子的机制，允许 Mk随时间变化的计算。这个变率及它的标准偏差的数学表达式可参考Binkowski等（1995）。MADE中假定在粒 子碰并过程中，粒子仍然满足对数正态分布。止匕外，仅仅考虑了由于布朗运动引 起的碰并。碰并过程的数学表达式可参考 Whitby等(1997)和Binkowski等(1995)。无机盐化学基于MARSo Binkowski等(1995)对它作了改进，利用热力学平衡的方法来计算硫酸盐/硝酸盐/ 氨盐 / 水气溶胶的化学组成。根据铵离子和硫酸根离子的比值， 考虑了两种机制： 当比值小于2 时， 模式处理为氢离子摩尔浓度的三次多项式，如果有足够的水和液态氨，硝酸盐溶解计

16、算，对离子态大于 50，假定硝酸盐不存在； 当比值大于等于2 时, 所有硫酸盐都假定为硫酸铵并在有水存在的条件下计算。活度系数的计算采用 Bromley 的方法。有机化学基于SORGAM SORGAM中假定SOA之间相互作用，使用了一种准理想的处理方法。气粒转化使用了 Odum 等( 1996)的参数化方法。由于缺乏活度系数的信息，所有的活度系数都假定为 1。SORGAM分开考虑人为源和自然 源产生的前体物对SOA的贡献，它是为RACM气相化学机制设计的。如果采用 RADM2机制，自然源SOA以及对应的前体物浓度则设定为 0。MOSAICMOSAIC采用分段的方法来对待气溶胶的尺度分布，尺度段

17、的数目是灵活可 变的，WRF-Chem中目前分8个或4个谱段来表征气溶胶粒子的尺度分布，表 1 为 8 个尺度段干气溶胶的尺度范围。MOSAIC中包含了硫酸盐、硝酸盐、钱盐、钠盐、氯盐、其他无机盐、有机碳、元素碳、水以及钙盐等气溶胶，它使用气溶胶多组分平衡方法 (MESA,Zaveri, et al., 2005b)和多组分的泰勒扩展方法 (MTET, Zaveri, et al., 2005球模 拟硫酸盐、硝酸盐、铵盐、钠盐、钙盐、氯盐和水气溶胶的热力学平衡。气相到颗粒相的质量输送和凝结采用的是自适应时间分裂式欧拉方法(ASTEEM)它是一种动态的方法，不同粒子尺度段或尺度群粒子的气粒转化的

18、时间特征量相近。通常解气粒转化微分方程使用一种有效的时间分裂方法， 这种方法不需要体积平 衡假定或混合处理。MOSAIC目前没有包含粒子核化、碰并和SOA形成等过程。Zhang and Wexler (2002)和 Wexler 等 (1994) 认为在短时间尺度模拟过程中，碰并过程是不重要 的， 并且在很多其它气溶胶模式中也没有考虑碰并过程。 尽管同质核化在新气溶胶粒子形成过程中是一个很重要的过程，但是目前对同质核化的机制及其核化率 是不确定的(Kulmala, et al., 2004)表1. 8个尺度段干气溶胶的空气动力学直径范围Table 1. Particle Dry-Diamete

19、r Range for the Eight Size Bins Employed by MOSAIC 尺度段 直径下界(um)直径上界(um)12345678参考文献Schell, B., Ackermann, ., Hass, H., Binkowski, ., Ebel, A. Modeling the formation of secondary organic aerosol within a comprehensive air quality model system. Journal of Geophysical Research, 2001, 106, 2827528293.Za

20、veri, R. A., and L. K. Peters. A new lumped structure photochemical mechanism for large-scale applications. J. Geophys. Res., 1999, 104(30), 38715.Zaveri ., . Easter, and . Wexler , A new method for multi-component activity coefficients of electrolytes in aqueous atmospheric aerosols, J. Geophys. Re

21、s., 2005, 110, D02201, doi:2004JD004681.Zaveri ., . Easter, and . Peters . A computationally efficient multicomponent equilibrium solver for aerosols (MESA). J. Geophys. Res. 2005, 110, D24203,doi:2004JD005618.Barnard, J. C., E. G. Chapman, J. D. Fast, J. R. Schmelzer,J. R. Schlusser,and R. E. Shett

22、er. An evaluation of the FAST-Jphotolysis model for predicting nitrogen dioxide photolysis rates under clear and cloudy sky conditions. Atmos Environ., 2004, 38, 3393W403.Wild, O., X. Zhu, and M. J. Prather. Fast-J: Accurate simulation of in- and below cloud photolysis in tropospheric chemical models. J. Atmos. Chem., 2000, 37, 245282.Jiang, G., and J. D. Fast. Modeling the effects of VOC and NOx emission sources on ozone formation in Houston during the TexAQS2000 field campaign. Atmos. Environ., 2004, 38, 5071- 5085.Fast, J. D., R. A. Zaveri