基于扰动模式的四维变分资料同化系统框架的设计完善和数值试验

为了建立一个应用于区域数值预报的四维变分资料同化（4DVar）系统,在近期开发的扰动预报模式GRAPES_PF基础上,开发完善增量四维变分同化系统框架。该框架中暂不包含物理过程（长短波辐射、边界层过程、对流参数化和云微物理等）。对比业务使用的GRAPES 3DVar系统,增加了温度控制变量。将无量纲Exner气压与流函数的线性风压平衡方程直接在地形追随垂直坐标面上求解,且通过广义共轭余差法（GCR）求解扰动亥姆霍兹（Helmholtz）伴随方程。利用人造“探空”资料对2015年10月台风“彩虹”进行了理想数值试验。试验结果表明,所开发的扰动四维变分同化框架得到了预期的结果,即同化更多资料并反复受到模式约束的四维变分同化系统能有效改善初值质量,进而改善区域数值预报。建立的区域四维变分同化框架合理可行,为进一步发展包含完整物理过程的区域四维变分同化系统奠定了研究基础。      

风廓线雷达资料质量控制及其同化应用

为更好地同化风廓线雷达观测资料开展了相应的质量控制与同化应用研究。针对2013年5月广东地区13部风廓线雷达的观测数据,采用经验正交函数(EOF) 分析方法对其进行质量控制。相比原始观测,经过质量控制的风场提高(降低)了来自时空大(小)尺度的贡献,较好地滤除了小尺度高频脉动,也较好地保留了大尺度平均状态与局地中小尺度系统的共同影响,并且更加接近ECMWF再分析场。此外,还对质量控制后的数据进行了垂直稀疏化。分别计算了质量控制前、后风廓线雷达观测与NCEP 6 h预报场的差值,对比差值的特征发现,经过质量控制的数据的观测增量更好地满足了高斯分布与无偏假设。针对一个实际天气个例,基于GRAPES 3D-Var同化系统,分析了质量控制后的风廓线雷达资料对模式分析与预报的影响。试验表明,在循环同化过程中加入风廓线雷达资料可以更好地描述模式初始场低层风场的特征,从而对强降水的位置与强度做出更好的预报。针对2013年5月的批量试验表明,同化风廓线雷达资料使短期降水预报有明显的改善。     

1998年中国区域降水低频变化的传播机制的初步分析

使用美国NCEP/NCAR再分析和中国实际降水的逐日资料,研究了1998年夏季中国区域低频降水传播变化的机制.结果显示,(1)中国中东部地区低频降水经向传播都与南海季风和印度季风的变化相对应,低频降水在经向传播中存在不连续现象,是由于低频下沉气流造成的;(2)中国中部区域低频降水纬向传播特征是在不同时间低频降水的纬向传播方向有明显差别,初步认为造成这种现象的可能机制是纬向风气流的变化;(3)高原地区与平原地区之间的低频雨带活动具有不连续性,这是由于在高原东坡存在低频下沉补偿气流造成的.     

GRAPES区域扰动预报模式动力框架设计及检验

设计了适用于四维变分同化系统的扰动预报模式GRAPES_PF。根据GRAPES的地形追随坐标非静力原始方程组,采用小扰动分离方法推导微分形式的线性扰动预报方程组,并利用与GRAPES非线性模式相似的数值求解方案求解线性扰动微分方程组。在设计扰动预报模式时采用了两个时间层半隐式半拉格朗日方案对动量方程、热力学方程、水汽方程和连续方程进行时间差分,空间差分方案的变量分布水平方向采用Arakawa C跳点网格,垂直方向采用Charney/Phillips跳层。利用代数消元法进一步推导得到只包含未来时刻扰动Exner气压的亥姆霍兹方程,进而通过广义共轭余差法（GCR）求解,在此基础上得到未来时刻扰动量的预报值。基于所开发的扰动模式开展了数值试验。首先在非线性模式中施加一个中尺度初始扰动高压,得到初始扰动在非线性模式中的后续演变,然后将相同的初始扰动作为扰动模式的初值进行时间积分,将扰动模式预报的结果与非线性模式的结果做了对比。结果表明,所开发的扰动模式GRAPES_PF较好地模拟了惯性重力内波的传播过程:初始高压扰动激发了一个迅速向外传播的惯性重力内波,在气压场向风场适应的过程中,水平风场、垂直运动、位温和湿度等变量均出现了扰动增量,与非线性模式得到的结果相当接近。GRAPES_PF作为GRAPES非线性模式的合理线性模式为建立基于线性扰动预报的区域四维变分同化系统奠定了科学基础。      

弱频散条件下大气运动的非线性演变过程

研究非线性过程对大气运动影响的一个常用方法是在一定条件下,将大气运动力学方程化为一些熟知的典型非线性演变方程,再利用数学中已有的结果进行讨论。Long、Redekopp等从地转涡度方程出发,证明了迭加在带状切变气流中的长的Rossby波的演变可以用著名的Kdv方程来描述。此     

华南中尺度暴雨数值预报的不确定性与集合预报试验

利用非静力MM5模式,分析了不同积云对流参数化方案对华南暖区暴雨数值预报的不确定性影响,进行了中尺度暴雨模式扰动集合预报试验。不同对流参数化方案的对流凝结加热引起不同的局地温度扰动,通过大气内部的热力动力过程,导致垂直速度的差异,进而影响网格尺度和次网格尺度降水时间、地点和强度。后续降水再通过凝结潜热释放形成新的扰动源。不同积云对流参数化方案还可引起扰动源能量传播方式不同,最终使模拟大气的动力和热力结构有差异。针对物理过程的不确定性,使用两种模式扰动方法构造集合预报扰动模式,第一种方法是随机组合不同积云对流参数化方案和边界层方案,第二种方法是扰动Grell积云对流参数化方案中主要参数振幅。集合预报结果表明,第一种方法的集合预报效果优于第二种方法,仅扰动参数振幅值似乎还不足以反映华南暴雨预报的不确定性。单一的确定性预报在暴雨落区和强度方面的可信度不稳定,集合产品能给华南暴雨过程提供更有用价值的指导预报,具有较高的应用价值。     

热带大气运动的Rossby孤波

广东热带海洋气象研究所     

区域极轨卫星ATOVS辐射偏差订正方法研究

近年来,卫星辐射资料在数值天气预报(NWP)系统中的直接同化研究取得了长足进展。为了利用TIROS业务垂直探测器(ATOVS)的辐射资料,必须对卫星观测辐射值的系统性偏差进行订正。在ECMWF原全球TOVS辐射偏差订正方案基础上,结合ATOVS资料特征和中国的实际情况,建立了适用于区域NOAA-15/16/17极轨气象卫星ATOVS辐射资料的偏差订正方案。该方案偏差订正分两步进行:首先进行扫描偏差订正,然后进行气团偏差订正。扫描偏差是临边测量相对于星下点测量的系统偏差,统计显示该种偏差具有一定的纬度依赖性,所以订正时按每10度的纬度带分别进行订正。气团偏差订正主要就是根据当时的天气条件进行订正,而天气条件一般用预报因子来定量表示。文中从中国国家气象中心T213背景场导出预报因子:(1)1000—300 hPa的厚度,(2)200—50 hPa的厚度,(3)模式地表温度,(4)总可降水量。模式预报因子的使用从观念上将对观测值的订正变为对计算前向辐射值的订正问题。试验结果表明,订正结果显著。     

THE TEST OF APPLYING RADAR TREC WIND IN THREE-DIMENSIONAL VARIATIONAL ASSIMILATION

1 INTRODUCTIONThat the radar data are reasonably ingested into theassimilation analysis system to improve the analysisquality and prediction precision is one of the problemsmeteorologists currently try to resolve. Many scholarshave already obtained some r…