基于WRF台风过程模拟研究.docVIP

基于WRF台风过程模拟研究 　　【摘 要】WRF （Weather Research and Forecast） 模式作为新一代的中尺度气象模式，在预测和分析中尺度气象方面具有良好的适用性。台风过程作为极端恶劣的中尺度气象现象，往往会造成巨大的人身和财产的损失。本文利用WRF模式对2011年的某次台风过程进行了模拟，分析模拟区域的选择，Nudging方法对台风路径模拟的影响，并和实际结果进行对比。 　　【关键词】WRF；台风路径；Nudging方法 　　0 引言 　　在大气动力学分析中，通常会将气象现象按照不同的时间和空间尺度进行尺度分析，从而在动力学方程中，采用不同的形式表征这些运动要素的特征和影响[1]。在中尺度气象模式中，大气运动具有较强的竖向非静力平衡性，尤其对于强烈的对流天气结构，浮力可以使气块产生较强的垂直加速度，导致猛烈天气的生成。 　　中尺度大气数值模式在 20 世纪九十年代已有相当发展，如美国国家环境预报中心用于业务预报的ETA模式，美国宾夕法尼亚大学和国家大气研究中心合作研制的MM5模式，科罗拉多州立大学研发的区域大气模拟系统――RAMS等中尺度模式，已经发展得非常成熟。WRF-ARW模式拟重点解决分辨率为1～10Km、时效为60h以内的有限区域天气预报和模拟问题。 　　对于台风路径的研究有着不同研究方法，采用数值模式的方式可以对台风进行诊断分析以及预测分析。对于预报模式可以通过确定的边界条件和数值及物理条件的设置模拟一个台风过程发展和移动过程。虽然方法为预测方法，如果在全球模拟中模拟了整个时间段的气象场，中尺度可以利用全球模拟的结果作为边界条件和初始条件以模拟更小尺度下的大气运动过程，如台风过程。利用中尺度模式更加精细的分辨率，对研究台风对地面结构的风荷载有重大意义。 　　1 WRF模式 　　WRF 模式系统[2]是 1997 年由美国国家大气研究中心（NCAR） 中小尺度气象处（MMM）、国家环境预报中心（NCEP）的环境模拟中心（EMC）、预报系统试验室的预报研究处（FRD）和俄克拉荷马大学的风暴分析预报中心（CAPS）四部门联合开发的新一代中尺度预报模式和同化系统。在其发展过程中，分别形成了一个是在NCAR的MM5模式基础上发展的ARW（Advanced Research WRF），另一个是在NCEP的Eta模式上发展而来的NMM（Nonhydrostatics Mesoscale Model），其动力框架如表1所示。本文中采用WRF-ARW模式。 　　表1 WRF动力框架 　　WRF的控制方程的通量形式如下所示： 　　WRF模式的初始条件和边界条件由大气气象资料提供，通常可以通过下载NCEP或者ECWMF的气象再分析资料获得，日本气象厅也提供了10km分辨率的分析资料GPV。 　　在WRF模式，有多种相关的物理参数化过程方案，如微物理过程方案，卷积云方案，表面层方案，大气边界层方案，大气辐射方案等等。这些方案的选择对中尺度气象现象的模拟有重要的影响。 　　FDDA（Four-dimensional data assimilation），也代指Nudging方法，可以在各时间步的模拟中使模拟结果向分析或者观测资料靠近。格点Nudging方法指的是在各个格点上使模拟结果靠近分析资料，观测Nudging方法指的是使模拟结果靠近观测数据。本文采用的是格点Nudging方法，其作用方程如下所示， 　　其中，θ代指的是Nudging方法作用的变量，F（θ）指的是不采用Nudging方法得到的变量的计算结果，Gθ指时间尺度对应的Nudging强度，Wθ指附加的调整的Nudging强度，0指Nudging目标文件在时间和空间上的插值结果。 　　2 模式设置 　　2.1 基础算例设置 　　本文在利用WRF模式模拟的台风过程是2011第12号台风Roke，主要影响的区域为日本本岛，该台风于日本静冈县登陆，而后穿越日本本州岛由宫城县再度进入太平洋。 　　模式的气象初始场采用的是日本气象局再分析资料GPV，其分辨率为5km，其气象要素按压力层分布包括高度，水平风速，温度，相对湿度，和海面更正气压值。边界条件也由该数据提供，每3小时更新一次。地表温度、湿度及海面温度由欧洲中尺度气象中心再分析资料提供（ECWMF-interm），其分辨率为0.7°，每6小时更新一次。 　　模式的地表静态资料采用的美国国家地质勘探局USGS资料，主要包括地面高程和地表覆盖信息，其最高分辨率为30s。 　　在本文的台风过程模拟中，笔者将着重研究不同模式区域的选取对台风过程模拟结果的影响，然后考虑Nudging方法作用的影响。为了排除多变量的复合影响，笔者先设置一个基础算例。算例的相关方案如表2所示。