快速检索
  气象   2016, Vol. 42 Issue (12): 1498-1505.  DOI: 10.7519/j.issn.1000-0526.2016.12.007

论文

引用本文 [复制中英文]

宫宇, 徐珺, 代刊, 等, 2016. 人工订正位涡改进模式初始场技术研究[J]. 气象, 42(12): 1498-1505. DOI: 10.7519/j.issn.1000-0526.2016.12.007.
[复制中文]
GONG Yu, XU Jun, DAI Kan, et al, 2016. Study on Improvement on the Model Initial Field by Manual Correction of Potential Vorticity[J]. Meteorological Monthly, 42(12): 1498-1505. DOI: 10.7519/j.issn.1000-0526.2016.12.007.
[复制英文]

资助项目

公益性行业(气象)科研专项(GYHY201306002) 资助

第一作者

宫宇, 主要从事天气预报工作.Email:heygy@live.cn

文章历史

2016年3月10日收稿
2016年11月07日收修定稿
人工订正位涡改进模式初始场技术研究
宫宇 , 徐珺 , 代刊 , 杨寅 , 杨舒楠     
国家气象中心,北京 100081
摘要:面对数值天气预报模式发展中的制约,如今气象界提出了预报员数据同化(data assimilation,DA)的理念,即将预报员的能力有机地融入到数值预报模式中以便发挥各自的优势。人工订正位涡改进模式初始场技术即充分体现了这一理念,其原理是利用了在中、高纬度气旋性环流场中,位涡(potential vortictiy,PV)场和卫星水汽图像(satellite water vapor imagery,WV)较好的对应关系检验数值预报,进而利用位涡的可反演性通过人工修改位涡场来订正数值预报,并已在部分国家实践中证明有着显著的改进模式预报性能的效果。本研究将通过分析对流层中高层的动力活跃区、位涡分布和卫星水汽图像暗区的关系检验模式偏差,进而通过人工订正位涡场和位涡反演,实现对模式初始场的修改,为预报员主观订正模式初始场、改进模式预报能力提供完整的技术流程。
关键词位涡    卫星水汽图像    位涡反演    模式偏差订正    
Study on Improvement on the Model Initial Field by Manual Correction of Potential Vorticity
GONG Yu, XU Jun, DAI Kan, YANG Yin, YANG Shunan    
National Meteorological Centre, Beijing 100081
Abstract: Faced with a series of limitations in developing numerical weather prediction model, a brand new conception of forecaster data assimilation (DA) has been put forward. That is an organic combination between weather forecasters' irreplaceable ability and numerical weather prediction. The improvement on the model initial field by manual correction of potential vorticity is coincident with this conception. By utilizing the principle of correspondence between potential vorticity (PV) and water vapor image (WV) observed by satellite, forecasters can exam the model and manually correct the distribution of PV, and use PV inversion to obtain a well developed model initial field. This technology has been proved efficacious in operational work in some countries. This study takes use of the principles mentioned above to achieve the whole technology and use some experimental results to show the technical processes.
Key words: potential vorticity (PV)    satellite water vapor imagery    potential vorticity inversion    model bias correction    
引言

预报员怎样与数值天气预报指南产品相互融合,并将各自价值发挥到极致,从而使两者“合力”做出的预报预警信息有“1+1>2”的效果,是学术界及各预报中心的最新关注焦点。在2015年召开的国际数值天气预报会议上,“预报员DA”的理念被低调提出(张永,2015)。“预报员DA”借助了DA理念,将预报员的预报经验,对一些天气过程中数值天气预报(numerical weather prediction, NWP)模式容易出现的“偏激”错误进行把握。此类将预报员能力与数值预报定量信息有机结合的理念或将成为预报员与预报模式长期共处的新机制。本研究的内容——人工订正位涡改进模式预报技术便是一种典型意义上的预报员DA技术。该方法的核心特点是利用等位涡面与卫星水汽图像的对应关系及位涡的可反演性将预报员的观测和长期经验快速、合理、有效地融入到模式分析场,实现预报员主观定性的预报能力同NWP的客观定量的特点进行有机的结合,从而实现提升模式预报性能。作为对NWP预报性能的有效订正和补充,该方法既充分发挥了预报员在实时预报中的不可替代的作用,又规避了NWP模式出现的“偏激”错误,十分有助于未来天气预报服务的定量化展开。

卫星水汽观测图像的观测原理决定了其可代表高层大气的运动状况,水汽图像所看到绝大部分重要结构与导致垂直运行和气流变形场的大尺度天气过程有关系。将水汽图像应用于天气尺度诊断分析的基础便是干区和湿区以及它们的边界。这些关键特征经常和重要的高空气流相关联,例如对流层顶动力异常,即位涡的入侵等过程。在热带气旋以外的气旋性环流场中,位涡场和卫星水汽图像之间有着很好的对应关系(吴君等,2007任素玲等,2014)。借助分析1.5 PVU或者2.0 PVU等位涡面的高度,可以较好地反映动力对流层顶的演变从而描绘位涡异常的主要特征。将水汽图像上的暗区位置与位涡分布进行比较,确定水汽图像中的暗区和模式输出中位涡异常区(对流层顶动力异常区域)之间是否匹配,可以检验模式对于天气尺度系统的预报强弱位置等特征;而将水汽图像上亮区位涡分布进行比较可以确定在急流的湿侧,动力对流层顶高的地方(等位涡面高)和水汽图像上灰度亮的地方是否匹配,从而检验模式对于活跃的斜压叶状云系统的预报情况。通过以上的对应关系,则开辟了一种全新的直观快速的定性模式检验途径,从而为下一步订正改进模式预报奠定了基础(许健民等,2008蒋建莹等,2014陶祖钰等,2012赵小平等,2014; 陈秀杰等,1998)。Hoskins等(1985)从理论上全面、清晰地阐述了位涡反演理论,指出在满足一定前提条件的情况下,就能从给定的位涡质量分布反演出相应的风、压和温度场,即位涡具有可反演性。因此通过模式位涡与卫星水汽图像的对比订正模式位涡,即可实现合理改进模式温压风全部物理量场的目的,开拓了主观订正模式的新思路,目前国外这项技术也已发展得较为完善并取得较好的应用效果(Hoskins et al, 1985; Georgiev et al, 2001)。

1 方法说明 1.1 人工订正位涡改进模式预报的技术流程

人工订正位涡(PV)改进模式预报的技术流程见图 1。技术流程通过客观检验方法包括要素检验,如风场、位势高度等检验或模式风场与观测风场对比, 以及PV与卫星水汽观测图像(WV)的主观对比发现模式误差后,通过PV和WV的对应关系人工修改PV,用修改后的PV场进行位涡反演并重新计算以得到新的预报,从而达到改进模式初始场的效果。

图 1 人工订正位涡改进模式预报技术路线流程 Fig. 1 The flowchart of manual correction of potential vorticity on the model initial field

为实现以上修改位涡高度并通过位涡反演得到风场和温度场,就要涉及到两大处理步骤:

① 如何使调整后的位涡较合理地融合入三维初始场。

② 对调整后的位涡进行位涡反演。

为了满足未来业务中人工修改位涡分布的合理性和交互便捷性的要求,采用最优线性无偏估计(BLUE)方法进行处理(Partick et al,2007; Arbogast et al, 2012)。其具有计算量小,运算速度快又兼顾合理性的特征。该方法结合了历史数据将人工订正的位涡合理地覆盖整个模式气柱,相较于常规插值,BLUE方法结合了模式历史物理量场统计特征,使得修改量更合理有效地融入模式场。具体方法如下:

(1) 预报员修改1.5 PVU上的高度。

该高度并不一定发生在标准模式垂直层次上,如发生在300~200 hPa,相应的200和300 hPa会产生同样的,表示根据线性关系插值求得更改高度上的PV,其中PVclosest level abovePVclosest level below分别表示距离1.5 PVU高度最近的上下两个模式面的PV

(2) 为保证PV垂直梯度不变,对改变高度的上下各4层采用同样的ΔPV

$ \begin{array}{l} \Delta PV = 1.5{\rm{PVU}} - \\ \;\;\;\;\;\;\;\;\;\;f\left({PV, P{V_{{\rm{closest level above}}}}, P{V_{{\rm{closest level below}}}}} \right) \end{array} $ (1)

(3) 通过最优线性无偏估计(BLUE)方法将PV订正信息插值到模式各层次上。其中,改变的8层作为观测矢量,原始模式场作为背景场,则订正后的PV场可写为:

$ {P_{{\rm{new}}}} = {P_{{\rm{old}}}} + \boldsymbol{B}{\boldsymbol{H}^{\rm{T}}}{(\boldsymbol{HB}{\boldsymbol{H}^{\rm{T}}} + \boldsymbol{R})^{ - 1}}\Delta P $ (2)

式中,R为观测误差矩阵,B为背景误差协方差矩阵,为8个PV订正矢量。采用1个月的数据作为样本计算BH为观测算子。假设R值小可忽略,则

$ {P_{{\rm{new}}}} = {P_{{\rm{old}}}} + \boldsymbol{B}{\boldsymbol{H}^{\rm{T}}}{(\boldsymbol{HB}{\boldsymbol{H}^{\rm{T}}})^{ - 1}}\Delta P $ (3)

由以上计算得到了预报员主观订正后的位涡场。

采用了分部式位涡反演方法对预报员主观订正后的位涡场进行反演,从而得到相应的温压风场。

1.2 位涡反演方法

目前常用的位涡反演方法主要是具有圆形轴对称结构的孤立位涡反演、准地转位涡反演和分部(片)位涡反演(piecewise PV inversion, PPVI)三种,在上述介绍的三种反演方法中,具有圆形轴对称结构的孤立位涡反演倾向于理论研究,而准地转位涡反演是线性的,对于实际的大气运动过程,分部(片)位涡反演PPVI方法(Arbogast et al, 2008; Davis, 1992a; 1992b;Davis et al, 1991; Wang et al, 2003; Charney, 1995; 寿绍文, 2010)不仅可以较合理地反演位涡,还能定量分析各种物理过程对大气运动的影响,有助于认识大气运动规律。因此研究采用的位涡反演方法是分部(片)位涡反演方法,具体方法参照了Davis等(1991)Wang等(2003)所述的反演原理。

下面以球坐标系中的大尺度运动说明分部(片)位涡反演的方程:对于大尺度运动,散度方程可简化为平衡方程,在球坐标系中,平衡方程可以写为:

$ {\nabla ^2}\mathit{\Phi } = \nabla \cdot\left({f\nabla \phi } \right) + \frac{2}{{{a^4}{\rm{co}}{{\rm{s}}^2}\varphi }}\frac{{\alpha (\partial \phi /\partial \lambda, \partial \phi /\partial \varphi)}}{{\partial \left({\lambda, \varphi } \right)}} $ (4)

式中,Φ为等压面位势,ϕ为流函数, φ为球坐标系纬度。

Ertel位涡在球坐标系中可近似写为:

$ \begin{array}{l} q = - \frac{{gk\pi }}{P}\left( {\eta \frac{{\partial \theta }}{{\partial \pi }} - \frac{1}{{a{\rm{cos}}\varphi }}\frac{{\partial v}}{{\partial \bar a}}\frac{{\partial \theta }}{{\partial \lambda }} + } \right.\\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\left. {\frac{1}{a}\frac{{\partial u}}{{\partial \pi }}\frac{{\partial \theta }}{{\partial \varphi }}} \right) \end{array} $ (5)

式中,η为绝对涡度的垂直分量,k=Rd/cp, π=cp(p/p0)k,取静力平衡近似于水平无辐散近似则式(5) 可以改写为

$ \begin{array}{l} q = \frac{{gk\pi }}{P}\left[ {\left({f + {\nabla ^2}\phi } \right)\frac{{{\partial ^2}\mathit{\Phi }}}{{\partial {\pi ^2}}} - \frac{1}{{{a^4}{{\cos }^2}\varphi }}\frac{{{\partial ^2}\phi }}{{\partial \lambda \partial \pi }}\frac{{{\partial ^2}\mathit{\Phi }}}{{\partial \lambda \partial \pi }}} \right] - \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\frac{1}{{{a^2}}}\frac{{{\partial ^2}\phi }}{{\partial \varphi \partial \pi }}\frac{{{\partial ^2}\mathit{\Phi }}}{{\partial \varphi \partial \pi }} \end{array} $ (6)

式(4) 和式(6) 构成以位势高度和流函数为未知函数的闭合方程组,其中q作为已知函数可从式(5) 计算得到。则已知位涡q,求解诊断方程式(4) 和式(6),可得到位势和流函数,从而得到满足平衡关系的温度场和风场。

分部反演中对于扰动量常用的是按时间或空间定义的平均和扰动,对任一瞬时变量A表示为

$ A\left({\lambda, \varphi, \pi, t} \right) = \bar A\left({\lambda, \varphi, \pi, t} \right) + A\prime \left({\lambda, \varphi, \pi, t} \right) $ (7)

式中,平均量定义为$\bar A = \frac{1}{\tau }\int_0^\tau {A{\rm{d}}t} $,其中τ表示扰动的生命周期,对于大尺度运动约为5 d左右。

变量的平均量不随时间变化,是与扰动无关的,将式(7) 代入式(4) 和式(6),可得到位势和流函数的平均量满足的诊断方程:

$ {\nabla ^2}\mathit{\bar \Phi } = \nabla \cdot\left({f\nabla \bar \phi } \right) + \frac{2}{{{a^4}{\rm{co}}{{\rm{s}}^2}\varphi }}\frac{{\alpha (\partial \bar \phi /\partial \lambda, \partial \bar \phi /\partial \varphi)}}{{\partial \left({\lambda, \varphi } \right)}} $ (8)
$ \begin{array}{l} \bar q = \frac{{gk\pi }}{P}\left[ {\left({f + {\nabla ^2}\bar \phi } \right)\frac{{{\partial ^2}\mathit{\bar \Phi }}}{{\partial {\pi ^2}}} - \frac{1}{{{a^4}{{\cos }^2}\varphi }}\frac{{{\partial ^2}\bar \phi }}{{\partial \lambda \partial \pi }}\frac{{{\partial ^2}\mathit{\bar \Phi }}}{{\partial \lambda \partial \pi }}} \right] - \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\frac{1}{{{a^2}}}\frac{{{\partial ^2}\bar \phi }}{{\partial \varphi \partial \pi }}\frac{{{\partial ^2}\mathit{\bar \Phi }}}{{\partial \varphi \partial \pi }} \end{array} $ (9)

将变量的平均量和扰动量的定义式代入式(4) 和式(6),并减去平均量方程式(8) 和式(9),则对应的变量扰动量的方程是:

$ \begin{array}{l} {\nabla ^2}\mathit{\Phi '} = \nabla \cdot\left({f{\nabla ^*}\phi '} \right) + \frac{2}{{{a^4}{\rm{co}}{{\rm{s}}^2}\varphi }}\left[ {\frac{{{\partial ^2}\phi '}}{{\partial {\lambda ^2}}}\frac{{{\partial ^2}\phi '}}{{\partial {\varphi ^2}}} - {{\left({\frac{{{\partial ^2}\phi '}}{{\partial \lambda \partial \varphi }}} \right)}^2}} \right.+\\ \;\;\;\;\;\;\;\;\;\;\;\left. {\frac{{{\partial ^2}\bar \phi }}{{\partial {\lambda ^2}}}\frac{{{\partial ^2}\phi '}}{{\partial {\varphi ^2}}} + \frac{{{\partial ^2}\phi '}}{{\partial \varphi }}\frac{{{\partial ^2}\bar \phi }}{{\partial \varphi }} - 2\frac{{{\partial ^2}\bar \phi }}{{\partial \lambda \partial \varphi }}\frac{{{\partial ^2}\phi '}}{{\partial \lambda \partial \varphi }}} \right] \end{array} $ (10)
$ \begin{array}{*{20}{c}} {q' = \frac{{gk\pi }}{P}\left[ {\left({f + {\nabla ^2}\bar \phi } \right)\frac{{{\partial ^2}\mathit{\bar \Phi }}}{{\partial {\pi ^2}}} + {\nabla ^2}\phi '\frac{{{\partial ^2}\mathit{\bar \Phi }}}{{\partial {\pi ^2}}}} \right. + }\\ {\;{\nabla ^2}\phi '\frac{{{\partial ^2}\mathit{\Phi '}}}{{\partial {\pi ^2}}} - \frac{1}{{{a^2}{\rm{co}}{{\rm{s}}^2}\varphi }}\left({\frac{{{\partial ^2}\phi '}}{{\partial \lambda \partial \pi }}\frac{{{\partial ^2}\mathit{\bar \Phi }}}{{\partial \lambda \partial \pi }}} \right. + }\\ {\left. {\frac{{{\partial ^2}\bar \phi }}{{\partial \lambda \partial \pi }}\frac{{{\partial ^2}\mathit{\Phi '}}}{{\partial \lambda \partial \pi }} + \frac{{{\partial ^2}\phi '}}{{\partial \lambda \partial \pi }}\frac{{{\partial ^2}\mathit{\bar \Phi }}}{{\partial \lambda \partial \pi }}} \right) - \frac{1}{{{a^2}}}\left({\frac{{{\partial ^2}\phi '}}{{\partial \varphi \partial \pi }}\frac{{{\partial ^2}\mathit{\bar \Phi }}}{{\partial \varphi \partial \pi }} + } \right.}\\ {\left. {\left. {\frac{{{\partial ^2}\bar \phi }}{{\partial \lambda \partial \pi }}\frac{{{\partial ^2}\mathit{\Phi '}}}{{\partial \lambda \partial \pi }} + \frac{{{\partial ^2}\phi '}}{{\partial \lambda \partial \pi }}\frac{{{\partial ^2}\mathit{\bar \Phi }}}{{\partial \lambda \partial \pi }}} \right)} \right]} \end{array} $ (11)

以上就是扰动位涡、扰动位势和扰动流函数所满足的诊断方程。按实际分析的需求,需进一步导出不同动力过程所满足的方程。设q′为总位涡扰动,qn为某一动力过程相对应的位涡扰动,且有$q' = \sum\limits_{n = 1}^N {{q_n}} $,该式对位势和流函数也成立。再引入式(12) 可以得到分属于不同动力过程的位涡、位势和流函数扰动所满足的诊断方程:

$ {\left[ {} \right]^*} + \left[ {^ - } \right] + \frac{1}{2}{\left[ {} \right]^\prime } $ (12)
$ \begin{array}{*{20}{c}} {{\nabla ^2}{\mathit{\Phi }_n} = \nabla \cdot\left({f\nabla {\phi _n}} \right) + \frac{2}{{{a^4}{\rm{co}}{{\rm{s}}^2}\varphi }} \times }\\ {\left[ {\frac{{{\partial ^2}{\phi ^*}}}{{\partial {\lambda ^2}}}\frac{{{\partial ^2}{\phi _n}}}{{\partial {\varphi ^2}}} + \frac{{{\partial ^2}{\phi ^*}}}{{\partial {\lambda ^2}}}\frac{{{\partial ^2}{\phi _n}}}{{\partial {\varphi ^2}}} - 2\frac{{{\partial ^2}{\phi ^*}}}{{\partial \lambda \partial \varphi }}\frac{{{\partial ^2}{\phi _n}}}{{\partial \lambda \partial \varphi }}} \right]} \end{array} $ (13)
$ \begin{array}{*{20}{c}} {{q_n} = \frac{{gk\pi }}{P}\left[ {\left({f + {\nabla ^2}{\phi ^*}} \right)\frac{{{\partial ^2}{\mathit{\Phi }_n}}}{{\partial {\pi ^2}}} + \frac{{{\partial ^2}{\mathit{\Phi }^*}}}{{\partial {\pi ^2}}}{\nabla ^2}{\phi _n} - } \right.}\\ {\frac{1}{{{a^2}{\rm{co}}{{\rm{s}}^2}\varphi }}\left({\frac{{{\partial ^2}{\phi ^*}}}{{\partial \lambda \partial \pi }}\frac{{{\partial ^2}{\mathit{\Phi }_n}}}{{\partial \lambda \partial \pi }} + \frac{{{\partial ^2}{\mathit{\Phi }^*}}}{{\partial \lambda \partial \pi }}\frac{{{\partial ^2}{\phi _n}}}{{\partial \lambda \partial \pi }}} \right)}-\\ {\left. {\frac{1}{{{a^2}}}\left({\frac{{{\partial ^2}{\phi ^*}}}{{\partial \varphi \partial \pi }}\frac{{{\partial ^2}{\mathit{\Phi }_n}}}{{\partial \varphi \partial \pi }} + \frac{{{\partial ^2}{\mathit{\Phi }^*}}}{{\partial \varphi \partial \pi }}\frac{{{\partial ^2}{\phi _n}}}{{\partial \varphi \partial \pi }}} \right)} \right]} \end{array} $ (14)

实际诊断过程中,分部反演方法的步骤如下:

(1) 根据式(4) 和式(6),由qq计算得到Φ, ϕ以及Φ*ϕ*

(2) 根据变量的平均量和扰动量定义和式(14) 计算得到Φ′, ϕ′以及Φ*, ϕ*

(3) 根据分析需要划分qn

(4) 根据式(12) 和式(13),由qn反演得到Φnϕn,从而诊断与某一动力过程qn相对应的扰动位势场和扰动风场。

2 人工订正位涡修改模式初始场个例实验 2.1 个例简介及偏差分析

分别选取了2014年9月2日华北暴雨过程和2016年5月21日江南大部的大暴雨过程进行了试验分析,从而展示人工订正位涡改进模式初始场技术的实际应用状况。两次试验分别使用的T639模式和GRAPES-GFS模式零场数据。

对于2014年9月2日华北暴雨过程,天气分析显示,500 hPa高空槽东移加深,低层则为低涡东移影响。从卫星水汽观测图像来看,随着高空槽的东移,槽后有较明显的弓状暗区逐渐增强东移,结合高层位涡扰动和与卫星水汽图像的对应关系,将9月1日20时T639模式零场数据分析的1.5 PVU等位涡面高度分布与同时刻卫星水汽观测进行叠加(图 2),可见T639模式位涡下侵位置略落后于观测的卫星水汽图像暗区位置,模式较实况位置偏西偏北,对于河套西部的位涡入侵估计过强,而对于河套地区的位涡入侵强度估计偏弱。从图 3a3b的对比来看,对于此次天气过程的关键影响系统,T639模式零场对500 hPa高空槽位置估计偏西,强度略偏强,这也映证了位涡同水汽图像暗区的对比检验结论。因此,可根据此检验结果认为T639模式9月1日20时零场天气系统位置较实况偏慢,预报员可根据此结论将位涡入侵区域向东调整,并可依据经验适当减弱河套西部地区位涡入侵的强度,增强河套地区位涡入侵强度。

图 2 T639模式2014年9月1日20时零场1.5 PVU的高度分布(等值线)和FY2G卫星同时刻的水汽观测图像 (a)人工订正位涡前,(b)人工订正位涡后 Fig. 2 Height of 1.5 PVU calculated with T639 model initial field superposed by FY2 water vapor imagery at 20:00 BT 1 September 2014 (a) before correction, (b) after correction

图 3 (a)2014年9月1日20时500 hPa天气观测分析, (b)同时刻T639模式零场的500 hPa高度场(黑色等值线)、温度场(彩色等值线)及风场(风标), (c)人工订正位涡后进行位涡反演后的500 hPa高空形势场 Fig. 3 Synoptic condition at 500 hPa at 20:00 BT 1 September 2014 (a) observation, (b) T639 model initial field, (c) after manual correction

而对于2016年5月21日江南大部的暴雨过程,天气分析显示,500 hPa高原系统东移加强,并于20日08时在四川盆地东部深厚的高空涡旋系统。从卫星水汽观测图像来看,高原系统的东移,其后部有较明显的团状暗区随之显著增强,结合高层位涡扰动和与卫星水汽图像的对应关系,将5月20日08时GRAPES-GFS模式零场数据分析的1.5 PVU等位涡面高度分布与同时刻卫星水汽观测进行叠加(图 4),可见模式位涡下侵在观测的卫星水汽图像暗区位置明显偏弱。从图 5a5b的对比来看,对于影响此次暴雨过程的高空涡旋系统,GRAPES-GFS模式零场对500 hPa涡旋强度明显偏弱,中心位置偏北,再次映证了位涡同水汽图像暗区的对比检验方法的正确性。因此,可根据此检验结果认为GRAPES-GFS模式5月20日08时零场天气系统强度明显偏弱,位置偏北,预报员可根据此结论适当增强观测卫星水汽图像暗区的位涡入侵强度,并适当向南调整。

图 4图 2,但为GRAPES-GFS模式2016年5月20日08时 Fig. 4 Same as Fig. 2, but for calculated with GRAPES-GFS model at 08:00 BT 20 May 2016

图 5图 3,但为2016年5月20日08时,图b模式为GRAPES-GFS Fig. 5 Same as Fig. 3, but for 08:00 BT 20 May 2016, GRAPES-GFS in Fig.b
2.2 人工订正及结果分析 2.2.1 2014年9月2日华北暴雨过程

图 2b,根据以上的检验结果,采用BLUE方法将河套西部位涡入侵强区域进行适当减弱并向南调整,再将河套地区上空的位涡入侵适当加强,整体上使得位涡入侵区域范围与实况卫星水汽图像观测暗区范围接近。另从一支气柱的位涡垂直廓线来看(图 6),采用BLUE方法修改后的垂直廓线基本沿袭了修改前廓线的形态特征,强度整体有所增强,大气高层的位涡强度有所减弱,体现出了BLUE订正方法的合理性。沿105°E做位涡的垂直剖面可见(图 7),采用BLUE方法修改后的位涡场分布较合理地呈现了位涡下探的形态,约38°N以北人工订正后的等位涡面整体上移,38°N以南则整体下移,同时并未出现突变的不连续面,说明人工向南订正位涡入侵区域的结果较为合理。

图 6 2014年9月1日20时T639模式某一气柱上的位涡扩线分布 (实线为人工订正前,虚线为人工订正后,灰色点线为订正前后的位涡变化) Fig. 6 Single PV profile of T639 model on 1 September 2015 (Solid line is before correction, long dashed line is after correction and dotted line is the PV change between before and after correction)

图 7 沿105°E做剖面的位涡分布 (红色线条为人工订正前,蓝色线条为人工订正后) Fig. 7 Vertical cross-section of vorticity along 105°E (Red line is before correction and blue is after correction)

在确定人工订正位涡结果合理后,则可以开始利用该位涡场进行分部式位涡反演获取相应的温、压、风场。图 3是对此次天气过程的位涡反演前后的天气形势对比。

图 3,人工订正前,我国中高纬区域呈两槽一脊形势,西北地区为一脊,5840 gpm脊点位置可达48°N以北,内蒙古中东部到河套西部地区为一槽,槽后为较强的偏北气流,中低纬度区域副热带高压较强,脊线约位于25°N左右,西脊点位于103°E左右。对比人工订正后的模式场,可以看到,在人工订正位涡位置和强度后高空形势整体出现了改变,位于河套西部地区的低槽被减弱向东调整,相应的风场有细微减弱,温度略升高。而位于我国西北地区的西风脊则有所加强,中低纬度的副热带高压则有所增强,副热带高压面积有所增大,脊线略微北抬,西脊点西进到102°E。

2.2.2 2016年3月21日江南暴雨过程

图 4b,2.1节的云图检验结果,采用BLUE方法将四川东南部的位涡入侵区域进行适当增强并向南调整,使得位涡入侵区域范围与实况卫星水汽图像观测暗区范围接近(图 4b),进而采用位涡反演获得人工订正后的的温、压、风、高度场。如图 5b,人工订正前,我国中高纬区域呈一槽一脊形势,西北地区为一脊,中纬度地区江南西部为一发展旺盛的闭合涡旋系统,低纬度地区5840 gpm等高线北界位于我国东南沿海一线,其外围偏西气流较为强盛。订正后的高空形势场可见,高空涡旋强度明显增强,区域范围向南拓展,涡旋外围气流较订正前明显增强,系统整体与实况更为接近。

从以上个例对比可见,分部式位涡反演方法较合理成功地反映了人工订正位涡所导致的大气系统的变化,有效订正了模式分析场的系统位置和强度,同时保证了大气环流整体的平衡性,证明人工订正位涡改进模式初始场的技术已基本实现。

3 讨论与结论

以上内容已基本实现了人工订正位涡改进模式初始场的工作,基于此,工作组将继续拓展该技术的应用与开发,如数值预报模式的云图模拟技术的应用。该技术利用数值预报模式输出的大气水凝物、辐射、风场等要素模拟出各类卫星的水汽和红外亮温图像。运用该技术,将业务数值预报模式的云图模拟图像与实时卫星观测进行对比检验,预报员不仅可非常直观地检验模式初始场天气系统位置较实况偏差,更可以直接检验得到强度的偏差,弥补了之前位涡入侵强度与水汽暗区强度只能依靠预报员经验确定对应关系的状况。进而预报员可通过修改位涡的位置和强度订正模式的初始场,并再次对此改进后的模式初始场进行云图模拟和检验至与实况接近,从而更大程度上实现改进模式预报性能的目的。

人工订正位涡改进模式初始场技术以类似于资料同化的思路将预报员的预报实践经验与数值模式预报定量特性有机地进行结合,更大限度地实现预报员和数值预报模式各自的特点能力,一定程度上突破了传统预报员的角色,更加突出了预报员天气学背景在业务工作中的不可替代性的作用。该技术在天气预报业务工作中有着较广泛的应用场景,如拓展模式检验岗位业务,每日例行检验观测云图和业务数值预报模式初始场,并做出相应判断和订正,从而改进业务数值预报模式的性能,同时为值班预报员提供模式订正建议;又如预报员可通过敏感性试验,将部分位涡异常系统人工剔除进行预报试验,进而探究天气过程的关键影响系统,为值班员提供重点关注天气系统建议等。

参考文献
陈秀杰, 耿勃, 叶惠明, 1998. 一次飑线天气过程的卫星水汽图像特征[J]. 气象, 24(6): 51-54. DOI:10.7519/j.issn.1000-0526.1998.06.012
蒋建莹, 汪悦国, 2014. 卫星水汽图像上两次暴雨过程的干、湿特征对比分析[J]. 气象, 40(6): 706-714. DOI:10.7519/j.issn.1000-0526.2014.06.007
任素玲, 蒋建莹, 许健民, 2014. 卫星水汽通道探测所揭示的高空流场在南亚高压东侧强降水分析中的应用[J]. 气象, 40(6): 697-705. DOI:10.7519/j.issn.1000-0526.2014.06.006
寿绍文, 2010. 位涡理论及其应用[J]. 气象, 36(3): 9-18. DOI:10.7519/j.issn.1000-0526.2010.03.002
陶祖钰, 周小刚, 郑永光, 2012. 从涡度、位涡、到平流层干侵入--位涡问题的缘起、应用及其歧途[J]. 气象, 38(1): 28-40. DOI:10.7519/j.issn.1000-0526.2012.01.003
吴君, 汤剑平, 邰庆国, 等, 2007. 切变线暴雨过程中湿位涡的中尺度时空特征[J]. 气象, 33(10): 45-51. DOI:10.7519/j.issn.1000-0526.2007.10.007
许健民, 方宗义, 2008. 《卫星水汽图像和位势涡度场在天气分析和预报中的应用》导读[J]. 气象, 34(5): 3-8. DOI:10.7519/j.issn.1000-0526.2008.05.001
张永. 2015. 世界气象预报十大趋势, 中国气象报, 2015年12月9日第四版.
赵小平, 沈新勇, 朱晶晶, 等, 2014. 2009年豫南一次强暴雨过程的位涡方程诊断分析[J]. 气象, 40(1): 38-47. DOI:10.7519/j.issn.1000-0526.2014.01.005
Arbogast P, Maynard K, Piriou C, 2012. About the reliability of manual model pv corrections to improve forecasts[J]. Wea Forecasting, 27(27): 1554-1567.
Arbogast P, Maynard K, Crepin F, 2008. Ertel potential vorticity inversion using a digital filter initialization method[J]. Quar J Roy Meteor Soc, 134(634): 1287-1296. DOI:10.1002/qj.v134:634
Charney J, 1955. The use of the primitive equations of motion in numerical prediction[J]. Tellus, 7(1): 22-26. DOI:10.3402/tellusa.v7i1.8772
Davis C A, 1992a. Piecewise potential vorticity inversion[J]. J Atmos Sci, 49(16): 1397-1411. DOI:10.1175/1520-0469(1992)049<1397:PPVI>2.0.CO;2
Davis C A, 1992b. A potential-vorticity diagnosis of the importance of initial structure and condensational heating in observed extratropical cyclogenesis[J]. Mon Wea Rev, 120(11): 2409-2428. DOI:10.1175/1520-0493(1992)120<2409:APVDOT>2.0.CO;2
Davis C A, Emanuel K A, 1991. Potential vorticity diagnostics of cyclogenesis[J]. Mon Wea Rev, 119(8): 1929-1953. DOI:10.1175/1520-0493(1991)119<1929:PVDOC>2.0.CO;2
Georgiev C G, Martín F, 2001. Use of potential vorticity fields, Meteosat water vapour imagery and pseudo water vapour images for evaluating numerical model behaviour[J]. Meteor Appl, 8(1): 57-69. DOI:10.1017/S1350482701001050
Hoskins B J, McIntyre M E, Robertson A W, 1985. On the use and significance of isentropic potential vorticity maps[J]. Quar J Roy Meteor Soc, 111(470): 877-946. DOI:10.1002/qj.49711147002
Partick Santurette, Christo G Georgiev.2007.United States: Academic Press Inc.Applying Satellite Water Vapor Imagery and Potential Vorticity Analysis, 1-200.
Wang X, Zhang D L, 2003. .Potential vorticity diagnosis of a simulated hurricane.Part Ⅰ: Formulation and quasi-balanced flow[J]. J Atmos Sci, 60(13): 1593-1607. DOI:10.1175/2999.1