WRF介绍-PPT

9.2WRF模式相关介绍主讲人：胡亮帆PPT制作：张瑜、陈德圆、裴坤宁9.2.1WRF模式的简介WRF（WeatherResearchandForecastingModel）模式是由美国环境预报中心（NCEP）,美国国家大气研究中心（NCAR）以及多个大学、研究所和业务部门联合研发的一种统一的中尺度天气预报模式。WRF采用Fortran90语言编写，它的特点是灵活、易维护、可扩展、有效以及适用计算平台广泛。其主要特色在于先进的数据同化技术、功能强大的嵌套能力和先进的物理过程，特别是在对流和中尺度降水处理能力方面更有优势。WRF模式适用范围很广，从中小尺度到全球尺度的数值预报和模拟都有广泛的应用.既可以用于业务数值天气预报，也可以用于大气数值模拟研究领域，包括数据同化的研究、物理过程参数化的研究、区域气候模拟、空气质量模拟、海气耦合以及理想实验模拟等。一些网址WRF主页：下载：:：模式的整体框架介绍主要由四部分组成：预处理系统（用于将数据进行插值和模式标准初始化、定义模式区域、选择地图投影方式）、同化系统（包括三维变分同化）、动力内核以及后处理（图形软件包）部分。模式的动力内核（或者框架）分为WRF—ARW（用于科学研究）和WRF—NMM（用于业务预报）两种模块。ARW和NMM均包含于WRF基础软件框架中，它们之间除了动力求解方法不同之外，均共享相同的WRF模式系统框架和物理过程模块。ARW(AdvancedResearchWRF)是在NCAR的MM5模式基础上发展起来,用于研究；NMM(NonhydrostaticMesoscaleModel)是在NCEP的Eta模式基础上发展起来。这里介绍的是ARW。WRF模式系统的主要组成模块TheWRFModelingSystemProgramComponents9.2.3WRF—ARW•ARW模块的动力框架采用完全可压缩、非静力平衡（带有静力平衡选项）欧拉模型，用具有守恒性的变量的通量形式表示。通量形式欧拉方程组：在Arekawa-C型网格点上既有风速矢量又有标量，但它们在网格点上的定义位置并不相同。水平风速的U,V分量分别定义在四方形单元格点区域的正交边界上，而温度、湿度、气压等标量则定义在四方形单元格点区域的中央。如左图的C网格分布，双向嵌套且子区域上的格距与其母区域上格距呈奇数比关系，模式在处理嵌套区域计算时，子区域从其母区域中获取初始信息，经过积分计算，模式计算后子区域上相邻9个单元格中的标量（或者相邻3个单元格边框处风速矢量）的平均值通过反馈机制返回给粗网格上重合的格点。此嵌套网格与母区域比为3:1•水平方向采用ArakawaC(荒川C)网格点（重点考虑1-10km）•垂直方向则采用地形跟随质量坐标。•时间积分方案上采用三阶或者四阶的Runge-Kutta算法。RK3方案对中央差分以及上风平流方案都具有较好的稳定性。其稳定时间步长大小比二阶蛙跃式时间步长方案要大2~3倍，可以节省机时。三阶Runge-Kutta积分方案9.2.4物理过程介绍在数值模式模拟天气过程时，往由于模式分辨率不足等原因，对次网格尺度的物理过程不能很好的描述，因此需要诸如辐射、边界层微物理等物理过程参数化来完善模拟的效果。WRF模式的主要物理过程包括：1.微物理过程2.积云对流3.辐射4.行星边界层（PBL）5.陆面过程物理过程之间的相互作用气相化学机理•CBM-IV(CarbonBondMechanismIV)•CBM-Z(CarbonBondMechanismversionZ)•RADM2(RegionalAcidDepositionMechanism)•RACM(RegionalAtmosphericChemistryMechanism)•SAPRC99(StatewideAirPollutionResearchCentermechanism)WRF-CHEM主要参数化方案—化学机理气溶胶化学机理•SectionalMOSAIC•ModalMADE/SORGAM•GOCART•MADE/VBS9.2.4.1微物理过程1-1Kessler暖云方案：该方案来源于COMMAS模式，是一个简单的暖云降水方案，考虑的微物理过程包括：雨水的产生、降落与蒸发，云水的碰并增长与自动转化，以及由凝结产生云水的过程。该方案显式预报水汽、云水和雨水，无冰相过程。1-2Purdue-Lin方案：该方案包括了对水汽、云水、雨、云冰、雪和霰6种类型水成物的处理。该方案是WRF模式中相对比较复杂的微物理方案，更适合于理论研究。1-3WRF单参数—3类水成物(WSM3)方案：该方案包括冰沉降和新的冰相参数化。与其他方案最大的不同之处在于，该方案是基于冰的质量含量而非利用温度的诊断关系计算冰粒子浓度。假设高于冰点的水成物为云水和雨，冰点以下的为云冰和雪，对包含冰过程的计算效率很高。可以对三类水成物（即水汽、云水或云冰、雨或雪）进行预报，被称为简单冰方案。但要注意的是，该方案缺少过冷水和逐步融化率过程。9.2.4.1微物理过程1-4WSM5方案：该方案与WSM3的简单冰方案类似，但由于将水汽、雨、雪、云冰和云水存储在5个不同数组，因此允许有过冷水的存在，并且允许雪下降到融化层以下进行逐步融化。该方案与PurdueLin方案不同的是，该方案对冰和水的饱和调整过程是分开处理的。另外，该方案在格距介于中尺度和可分辨云尺度的格点计算效率很高。1-5WSM6方案：该方案是对WSM5方案的补充和扩展，另外包括了霰及其相关过程。对于较粗分辨率的中尺度格点而言，WSM3、WSM5、和WSM6方案差别不大，但是对云分辨尺度它们差别显著。如果考虑计算效率和理论基础，这三种方案中WSM6方案是最适合于云分辨格点的。1-6Eta格点尺度云—降水方案：预报模式平流项中的水汽和总凝结降水变化。程序中用一个局地存储数组保存初始猜测场信息，即分解得到的云水、雨、云冰和冰相降水物密度的贡献。（冰相降水物的密度是根据存有冰的总增长信息的局地数组估算的，冰的总增长信息又与气相沉降和液相碰撞有关）9.2.4.1微物理过程1-7Thompson方案：该与早期单参数总体微物理参数化（BMP）方案相比，该方案不仅改进了物理过程和程序代码，还加入了一个查算表以提供更多详细的谱分档机制。该方案假设雪粒的尺度分布不仅与冰水含量有关，而且还与温度有关系，并进一步假设雪粒为非球形，其总体密度与粒径大小有反相关，与观测一致。（与之相反，几乎其他所有BMP方案都假设雪粒为球形，并且密度为常数）1-8Goddard积云集合（GCE）模式方案：第一，可以选择霰或者冰雹作为第三种粒子。第二，加入了新的饱和调整方法，确保（过饱和）不饱和不会发生在天气晴朗（多云）的网格点上。第三、所有不涉及融化、蒸发和凝华的微物理过程都只在一种热动力状态上进行计算，保证所有这些过程都能被平等处理。第四、某一种水成物的所有沉降过程的总和不能超过它本身的质量，确保在微物理过程计算中保持液水总量的平衡。1-9Morrison双参数方案：该方案考虑了六种水成物：水汽、云滴、云冰、雨、雪和霰/冰雹，通过自定义开关来决定霰和冰雹的选择。诊断变量主要包括云冰、雨、雪和霰或冰雹的数浓度和质量混合比，云水和水汽的混合比（总共10个变量）。这种双参数的预报能更准确地处理粒子的谱分布，而粒子的谱分布又是计算各种微物理过程发生概率和估算云/降水的关键因子9.2.4.2积云对流2-1Kain-Fritsch方案：该方案是质量通量类型，在Eta模式中进行测试调整，采用一个含有水汽上升和下降过程的简单云模式，包括卷入和卷出，以及相对粗糙的微物理过程。2-2Betts-Miller-Janjic方案：该方案是对流调整方案，其最主要的改进在于引入成云效率参数，这样就增加了一个决定大气加热和水汽目标廓线的自由度。其中，浅对流调整是参数化的重要部分。2-3Grell-Devenyi集合方案：该方案是在每个网格点上运行多个积云参数化方案，并将其结果进行集合平均，然后反馈到模式格点。2-4Grell-3方案：该方案是WRF第三版本首次引入的新方案，与Grell-Devenyi方案有很多相同之处，都是基于集合平均方法。该方案更适用于格距小于10km的情况。9.2.4.3辐射3-1快速辐射传输长波模式（RRTM)：该模式来自MM5模式，采用谱段处理方案。它是利用一个预设好的查算表准确地表示由于水汽、臭氧、二氧化碳和痕量气体（大气中含量极少的气体）引起的长波辐射过程，同时也能表示云的光学厚度。3-2GFDL长波辐射方案：该方案来自美国地球流体动力实验室，计算与二氧化碳、水汽、和臭氧相关谱段上的长波辐射，采用简化的交换方法。该方案中云的重叠是随机的。3-3CAM长波辐射方案：该方案是用于NCAR大气环流模式的谱段处理方案，主要用来进行气候模拟。它能够处理几种痕量气体，并与可分辨云和云量进行相互作用。3-4GFDL短波辐射方案：该方案同样来自于美国地球流体动力实验室，考虑了二氧化碳、水汽和臭氧的效应，采用积分时间间隔内日间平均的太阳天顶角余弦值计算段波辐射。该方案中云的重叠是随机的。9.2.4.3辐射3-5MM5Dadhia短波辐射方案：该方案来自于MM5模式，对短波辐射通量向下进行简单积分。它能够反映晴空散射、水汽吸收和云的反射与吸收。另外，还考虑了地形坡度和阴影对地表短波辐射通量影响。3-6Goddard短波辐射方案：该方案共有11个谱段，采用二流近似方法计算太阳短波辐射的散射和直接辐射分量，并且方案中考虑了已有臭氧垂直廓线的气候分布。3-7CAM短波辐射方案：该方案也是用于NCAR大气环流模式的谱段处理方案，主要用来气候模拟。它能够处理几种气溶胶和痕量气体的光学特征。该方案特别适用于区域气候模拟。9.2.4.4行星边界层4-1MRF边界层方案：该方案最大的优点在于其有效弥补了K理论的不足，解决了由大涡输送所导致的逆梯度输送问题，对于不稳定或对流这样混合很好的边界层条件，该参数化能够给出比较理想的模拟结果。其不足之处在于，当风速较大时，MRF参数化方案存在过度混合的问题，会导致对流性降水减少。4-2Yonsei大学（YSU）边界层方案：该方案是MRF边界层方案的改进版本，在MRF边界层方案中增加了对边界层顶夹卷层的显式处理，从而有效解决了MRF方案中过度混合的缺点。4-3MYJ边界层方案：该方案是基于1.5阶湍流闭合的边界层参数化模式，从2阶闭合方案简化而来。相对于完整的2阶闭合方法，它计算量小，有一定的精确性。4-4ACM2边界层方案：该方案是一种非局地闭合方案，在对流的条件下，既能模拟由大尺度的湍流涡旋产生的输送过程，又能反映次网格小尺度的湍流混合过程。因此，该方案能真实模拟边界层热通量和位温廓线。9.2.4.5陆面过程5-15层热量扩散方案：该方案基于MM5的5层土壤湿度模式，分别是1、2、4、8、16cm，在第5层底的土壤温度固定为一个气候平均值。能量收支包括辐射、感热和潜热通量，同时也允许雪盖效应。5-2Noha陆面过程模式：该模式是OSU陆面过程模式的后继版，与原先的相比，可以预报土壤结冰、积雪影响，提高了处理城市地面的能力，考虑了地面发射体的性质，这些是OSU所没有的。5-3快速更新循环（RUC）陆面过程模式：该模式为多层土壤模式（6层，可增至9层以上），各层深度为0、5、20、40、160、300cm，在土壤上层垂直分辨率较高，另外还包含多层雪模式。它考虑了土壤结冰过程、不均匀雪地、雪的温度和密度差异，以及植被效应和冠层水。5-4P