T639全球模式的台风初始化方案升级试验

由于洋面上缺少有效的台风观测资料,台风初始化对全球数值预报业务模式而言,仍然是一个悬而未决的问题。国家气象中心自从2009年将T639全球数值预报系统投入业务运行以来,应用了一套完整可行的台风初始化方案,其技术路线是在台风刚发生时的预报时刻,在背景场嵌入人造bogus涡旋;而在后续的循环滚动预报时刻,采取涡旋重定位和涡旋调整技术方案;因此,第一个预报时刻的涡旋形成技术直接影响到后续时刻的台风涡旋质量乃至预报效果。但是,人造bogus涡旋主要依据统计的经验模型建立,其结构与全球数值模式的动力学、物理学属性并不匹配。基于国家气象中心T639全球资料分析同化-预报循环系统,一套升级版的台风初始化方案被发展了起来,与业务方案的主要区别在于升级方案利用资料同化技术来形成初始涡旋而不是直接嵌入人造bogus涡旋;这在很大程度上减少了人为主观因素对台风初始结构的影响,而更多地是依靠变分资料同化来协调约束产生三维空间的涡旋环流,这样产生的涡旋环流不但与周围环境场比较协调,而且与模式的动力学、物理学属性也比较匹配。应用新方案,本文对生成于西北太平洋2011-2012年27个台风进行了数值试验,初步分析表明,新方案在实现台风初始涡旋环流结构方面有着不错的表现,相比于业务使用的方案而言,新方案台风路径平均预报误差有了不错幅度的降低,2~5 d预报平均路径误差普遍降低了3%~15%。     

一种基于位涡反演的变性台风初始化方案

针对变性台风的结构特征和Bogus 方案的局限性,提出了一种新的基于位涡反演技术的台风初始化方案。实际应用效果表明,该变性台风初始化方案既能使初始场中的台风中心位置和强度与观测相吻合,又能维持变性台风的结构特征(斜压性和不对称性),同时能保证初始场中各要素场之间的协调一致性。数值试验结果表明,采用该方案后,在数值模拟的前24 小时对台风中心强度和移动路径的模拟效果有较明显改进。      

CMA-GFS全球预报系统中的台风初始化

在CMA-GFS(CMA Global Forecast System)全球四维变分资料同化系统(4DVar)基础上,参照BDA(Bogus Data Assimilation)方法,建立了一个全球模式台风初始化方案.该方案通过4DVar同化窗口吸收诊断处理后的1 h间隔台风中心定位及中心气压信息,利用模式动力物理约束...     

非对称台风风场的动力初始化应用研究

在对比NCAR-AFWA三维Bogus(虚拟)风场及风压场方案的基础上,虚拟台风模型中的水平风场采用考虑台风移速、移向、摩擦等客观因素的非对称方案,垂直风场采用考虑高空反气旋环流的风廓线方案,然后利用Nudging动力初始化技术形成新的模式初始场.最后对"云娜"台风进行数值模拟,结果表明:通过合理引入客观因子及观测事实对台风Bogus模型进行构建,并利用同化方法使各个物理量相互协调,可以使初始场与实况更加接近,从而在一定程度上提高台风路径及强度的预报水平,并对降水分布的预报效果也有所改善.     

几种台风初始化方法及其对比实验

利用中尺度MM5V3模式及其伴随系统，运用四维变分同化方法，将Bogus台风作为观测资料进而优化台风初始场，对0422号台风个例进行了数值模拟和试验，并与传统Bogus方法和Nudging方法进行比较。结果表明，经过BDA方案优化的台风初始场更加接近实际台风的结构，有效地提高了台风路径和强度预报的准确率。     

台风涡旋循环初始化方法及其在GRAPES-TCM中的应用

由于观测资料的缺乏,全球数值模式中台风涡旋一般较弱.但在区域数值模式中由于能使用更高的分辨率,采用适合台风的物理过程计算方案,因此区域模式中的台风涡旋能得到更好的发展.通常情况下,由于区域模式每次预报都使用全球模式的预报结果作为背景场,因此在区域模式中得到发展的台风涡旋并不能出现在下一次预报的初始场中.本研究结合BDA和涡旋重定位方法,采用MC-3DVar技术将台风涡旋初始化与资料同化技术结合起来,提出了循环初始化技术.在该技术中,提取区域模式6 h预报的台风涡旋周围的廓线资料,并根据台风观测位置将这些廓线资料平移到正确位置进行同化,从而使得台风涡旋在区域模式中得到循环发展.个例试验的结果表明循环初始化方案能更准确地刻画初始台风的强度和结构,相比控制试验,其初始台风强度更接近实况,结构更为合理,而且与观测更为接近,相应的,台风路径和强度预报也有很大的改善.对超强台风"蔷薇"的22次预报的总体结果表明,相对控制试验,循环初始化方案使得试验的总体路径预报距离误差从24 h的113 km降低到89 km,从48 h的188 km降低到184 km.另一方面,循环初始化方案对强度预报也有很大改善,24 h中心气压均方根误差从控制试验的32 hPa降到22 hPa.2009年汛期开始,以循环初始化方案为基础的GRAPES-TCM投入准业务试验,在对2009年2、3、4和6号热带气旋的准业务运行的预报结果的分析表明,新的方案较之采用涡旋重定位的GRAPES-TCM业务系统在路径和强度预报能力方面均有较大改善.     

四维变分同化方法在台风初始化及其数值模拟中的应用--对9608号(Herb)台风的研究结果

The MM5 and its four dimensional variational data assimilation (4D-Var) system are used in this paper. Based on the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis data, the authors generate an optimal initial condition for a typhoon by using the bogus data assimilation (BDA) scheme. BDA is able to recover many of the structural features of typhoons including a warm-core vertex, the correct center position, and the strong circulation.As a result of BDA using a bogus surface low, dramatic improvement is achieved in the 72 h prediction of typhoon Herb. Through several cases, the initialization by BDA effectively generates the harmonious inner structure of the typhoon, but which is lacking in the original analysis field. Therefore the intensity forecast is improved greatly. Some improvements are made in the track forecast, but more work still needs to be done.     

反演台风海平面气压场与NCEP及报文资料的对比

利用QuikSCAT散射计获取的海面风场资料,借助UWPBL模式,对2006年3个台风个例的海平面气压场进行了反演。通过将反演结果与NCEP海平面气压场及台风报文资料的对比分析,表明反演的台风海平面气压场准确性较NCEP海平面气压场高,但与报文资料相比,中心气压明显偏高,因此这种利用散射计风场反演海平面气压的方法虽有一定的准确性,但精确度还有待提高。     

一种基于分析增量更新技术的台风初始化方案

台风初始化方案最大的困难在于其只能根据有限的观测资料（台风中心位置、最低气压、最大风速和大风半径等）来构造一个与模式动力-物理过程协调的涡旋模型。首先根据实际观测资料对背景场中的涡旋扰动进行重定位和风速调整,然后将分析增量更新（Incremental Analysis Update,IAU）技术应用于台风初始化方案中,将调整后的涡旋当作一个强迫项逐渐加入到模式预报过程中,通过模式自身调整来得到一个协调性较好的台风初始结构,从而改善模式对台风的预报性能。对“山竹”台风的多次预报结果表明:（1）台风初始化对于台风路径误差影响较小,对于强度预报改进则比较明显。特别是在台风生成初期,台风初始化技术能够有效地增强全球模式分析场中的涡旋强度,并解决预报过程中强度较弱的问题。（2）根据预先给定的三维风场,IAU技术能通过模式自身预报过程对其他变量进行调整,从而得到一个热力和动力协调的初始涡旋结构。相对于仅对初始风场的调整,它对24 h之后的路径和强度预报误差会有更进一步的改善。（3）对IAU中的松弛时间进行参数敏感性试验,发现该变量取3—6 h效果较好。      

采用预报涡旋的初始化方案对2015年台风“莲花”、“灿鸿”的试验研究

提出了一种采用预报涡旋的初始化方案,用预报涡旋代替bogus模型参与构建模式初始场,采用权重形式合成预报涡旋和分析涡旋获取台风初始涡旋。针对2015年“莲花”和“灿鸿”台风,基于该初始化方案设计了一系列对比试验进行数值模拟,并对结果进行分析。结果表明:(1)该方案得到的台风初始涡旋结构比bogus模型合理;(2)预报涡旋权重不宜取太大;(3)从长时效预报效果看,采用24 h内预报涡旋比采用长时效预报涡旋台风的路径和强度误差减小;(4)采用同一权重方案对“莲花”、“灿鸿”预报的改进效果不同,其原因与预报涡旋和分析涡旋的协调程度有关。多台风情形下可在初步评估的基础上采用不同时效的预报涡旋和不同权重方案。      

云雨条件下AMSR2微波成像资料同化试验及其在台风预报中的应用

采用WRF-ARW和WRFDA的3.8版本进行云雨条件下AMSR2(Advanced Microwave Scanning Radiometer 2)卫星微波成像资料的同化试验,并选择2014年7月第9号台风"威马逊"为研究对象,检验同化云雨条件下的AMSR2资料对台风预报结果的影响。结果表明:(1)云雨条件下的质量控制方案所保留的像元明显多于晴空条件下质量控制方案,且亮温模拟结果更加接近观测;(2)云雨影响的加入使得风场、高度场和海平面气压场的模拟结果好于晴空条件下模拟,且对台风中心有增温、增湿的作用。(3)云雨条件下模拟对台风的路径和强度预报均有改善。     

The effect of three-dimensional variational data assimilation of QuikSCAT data on the numerical simulation of typhoon track and intensity

In this paper, the three-dimensional variational data assimilation scheme (3DVAR) in the mesoscalemod el version 5 (MM5) of the US Pennsylvania State University/National Center for Atmospheric Research is used to study the effect of assimilating the sea-wind data from QuikSCAT on the prediction of typhoon track and intensity. The case of Typhoon Dujuan (2003) is first tested and the results show appreciable improvements. Twelve other cases in 2003 are then evaluated. The assimilation of the QuikSCAT data produces significant impacts on the structure of Dujuan in terms of the horizontal and vertical winds, sealevel pressure and temperature at the initial time. With the assimilation, the 24-h (48-h) track prediction of 11 (10) out of the 12 typhoons is improved. The 24-h (48-h) prediction of typhoon intensity is also improved in 10 (9) of the 12 cases. These experiments therefore demonstrate that assimilation of the QuikSCAT sea-wind data can increase the accuracy of typhoon track and intensity predictions through modification of the initial fields associated with the typhoon.     

基于KPCA的台风强度神经网络集合预报方法研究

针对影响台风强度前期预报因子较多以及因子的非线性变化特点,首先采用逐步回归方法筛选出部分预报因子,再利用核主成分分析方法在剩余的预报因子中提取包含了原数据较多信息的核主成分与前期选入的预报因子共同作为模型输入.进一步考虑到神经网络集合预报中个体的准确性和差异性的权衡问题,在不同的初始条件下生成若干组神经网络,分别选择每组中性能最优的个体,建立了一种新的非线性神经网络集合预报模型.最后以西北太平洋海域2001-2010年5-10月的台风强度为研究对象进行了预报试验.结果表明,这种神经网络集合预报模型的预报结果符合实际应用的要求,其预报平均绝对误差明显小于同等条件下的神经网络方法和逐步回归预报方法.     

台风“碧利斯”的雷达反演及改进

利用两步变分反演方法,研究2006年第4号台风“碧利斯”,对照单部雷达资料和两部雷达资料的反演结果。结果表明：利用单部雷达资料不能很好的反演出台风的风场结构特征,而两部雷达资料在第一步反演背景场时,可以看出台风的涡旋特征。通过检验发现,此次台风过程的水平切向风在最大风速半径内部与Rankine模型相近,外部与Chen 3模型接近。为了更有效地利用单部雷达的观测信息准确定位台风的移动路径,故在单部雷达资料第一步反演背景场时加入模型弱约束。反演结果表明,加入模型弱约束后,反演的背景风场就能看出台风的涡旋特征,通过第二步反演后与两部雷达资料结果相近;加模型弱约束的单部雷达反演结果与两部雷达反演结果间u、v分量的相关系数均在0.8或以上,2-7 km高度层u、v分量的均方根误差明显减小。     

台风“碧利斯”的雷达反演与改进

利用两步变分反演方法研究2006年第4号台风＂碧利斯＂,比较单部雷达资料和双部雷达资料的反演结果。结果表明：利用单部雷达资料不能很好地反演出台风的风场结构特征,而双部雷达资料在第一步反演背景场时,可以显示出台风的涡旋特征。通过检验发现,此次台风过程的水平切向风在最大风速半径内部与Rankine模型相近,外部与Chen 3模型接近。为了更有效地利用单部雷达的观测信息准确定位台风的移动路径,故在单部雷达资料第一步反演背景场时加入模型弱约束。反演结果表明,加入模型弱约束后,反演的背景风场可以看出台风的涡旋特征,通过第二步反演后与双部雷达资料结果相近;加模型弱约束的单部雷达反演结果与双部雷达反演结果间u、v分量的相关系数均为0.8或以上,2—7 km高度层u、v分量的均方根误差明显减小。     

基于ECMWF模式预报的台风降水地形订正方法

为了解决复杂地形条件下ECMWF模式预报的台风降水较实况显著小的问题,对Smith 1979年提出的地形降水方程进行改进,提出以饱和湿层高度作为方程积分上限,针对不同高度地形设定不同的降水效率;以无量纲湿弗劳德数大于1作为有、无地形降水的判据;利用ECWMF细网格预报场,通过迎风坡地形降水估算方程来订正模式预报的台风降水。用该地形降水订正方法对1617号台风“鲇鱼”的降水进行了订正预报。结果表明,虽然在一些小尺度地形区域会产生明显的空报,但是对于大尺度地形区域的强降水有显著的订正效果。对1513、1521和1614台风的订正结果进一步表明,该地形降水订正方法对改进台风极端降水预报效果显著。需要指出的是,采用的地形降水订正方法仅考虑了稳定条件下的地形降水,对于其他情形下的地形降水订正方法尚待进一步的研究。      

GRAPES区域集合预报条件性台风涡旋重定位方法研究

为了在集合预报中更合理描述台风涡旋中心定位的不确定性,采用2009—2018年中国气象局和日本气象厅台风最佳路径数据,分析台风最佳路径涡旋中心定位的不确定性特征,在此基础上设计条件性台风涡旋重定位方法(Conditional Typhoon Vortex Relocation,CTVR),构建集合成员台风涡旋中心重定位阈值条件、台风涡旋分离数学处理及涡旋重定位等数学处理过程,利用中国气象局数值预报中心区域集合预报系统(Global/Regional Assimilation and Prediciton System-Regional Ensemble System,GRAPES-REPS)对2018年西北太平洋上的3个台风(1808号"玛莉亚"、1824号"谭美"和1825号"康妮")进行轴对称结构和轴对称+非对称结构条件性台风涡旋重定位两种方案的集合预报试验和检验评估。结果表明:(1)中国气象局和日本气象厅台风最佳路径误差平均值为13.72 km,可视为台风涡旋中心定位不确定性的合理估计值;(2)统计检验结果和典型个例分析表明,采用轴对称结构和轴对称+非对称结构条件性台风涡旋重定位...     

基于降水极端预报指数的福建台风极端降水预报研究

应用1961—2017年中国气象局热带气旋最佳路径数据集、国家地面气象观测站日降水观测资料和2015年8月—2017年12月欧洲中期天气预报中心(ECMWF)集合预报系统降水极端预报指数(EFI)数据,根据百分位法定义台风影响期间福建省各站点的台风极端降水阈值,采用最小阈值法剔除台风极端降水时EFI箱线图中的异常值,保...