垂直分层加密和预报区域扩大对GRAPES_TYM台风预报的影响*

麻素红 张 进 瞿安祥 王大鹏 沈学顺

MA Suhong ZHANG Jin QU Anxiang WANG Dapeng SHEN Xueshun

国家气象中心，北京，100081

National Meteorological Centre，Beijing 100081，Beijing

1 引 言

国家气象中心于2010 年基于GRAPES-MESO（陈德辉等，2006）开始研发区域模式台风数值预报系统（GRAPES_TYM），该系统于2012 年投入业务运行并得到持续改进，包括模式对流参数化方案（张进等，2017 年）、模式参考大气廓线和涡旋初始化（麻素红等，2018；麻素红，2019；Wang，et al，2020）等，GRAPES_TYM 热带气旋路径及强度预报能力不断提升。但是，该系统仍然存在较明显的系统偏差，如对超强台风（最大风速＞51.0 m/s）的强度极值预报偏低，导致强度预报存在较明显的负偏差（麻素红等，2018）。

热带气旋强度预报很具有挑战性（DeMaria，et al，2005；Rogers，et al，2006），其中模式分辨率、模式物理过程以及模式初始场对热带气旋的描述是影响热带气旋强度预报的关键因素（Cha，et al，2013；Gopalakrishnan，et al，2012；Liu， et al，2000，2012；梁旭东等，2010）。而高分辨率模式是开展热带气旋强度及结构变化研究以及业务预报的基本条件（Gopalakrishnan，et al，2012；Tallapragada，et al， 2012； Wang， 2001， 2002； Yau， et al， 2004；Zhang，et al，2011）。研 究（Zhang， et al，2003，2015；Kimball，et al，2006）表明，模式垂直分辨率的提升对热带气旋强度预报至关重要。在模式不同高度层加密，对热带气旋强度预报的作用会有明显不同，具体表现为：对模式中低层加密，模式预报的热带气旋较强；对模式高层加密，模式预报的热带气旋较弱。如2017 年NCEP （National Center for Environment Prediction）区域模式热带气旋数值预报系统HWRF（Hurricane Weather and Research Forecast）的垂直分辨率由61 层提高到75 层后，模式对热带气旋路径及强度预报误差明显减少（Zhang B L，et al，2016）。

随着热带气旋预报业务的不断拓展，GRAPES_TYM 需要为北印度洋热带气旋业务预报提供参考产品。这就需要将该系统的预报范围由覆盖西北太平洋和中国南海扩大到覆盖北印度洋。目前基于区域模式开展不同海域热带气旋数值预报的主要策略是分区域运行，如美国NCEP 的HWRF 模式针对大西洋、东北太平洋、西北太平洋3 个海域的热带气旋预报分别独立运行3 个区域模式（Tallapragada，et al，2012，2014a，2014b，2015，2016）。但Zhang X J等（2016）基于HWRF 建立了覆盖大西洋、东北太平洋的大区域热带气旋数值预报试验系统（简称HWRF-B）。试验结果显示：当洋面上2 个及以上热带气旋同时存在时，HWRF-B 的平均路径误差比HWRF 相应的平均路径误差小（Ghassan，et al，2017）。

2019 年针对上述两个问题，国家气象中心对GRAPES_TYM 进行了如下改进：（1）在模式层顶不变的条件下，模式垂直分辨率由50 层增加到68 层；（2）预报区域由覆盖西北太平洋、中国南海扩大到覆盖北印度洋。改进后的系统于2019 年8 月投入业务运行。

文中分析了模式垂直分辨率提升以及预报区域扩大对GRAPES_TYM 路径及强度预报的影响，并与NCEP-GFS 和ECMWF 进行了对比。其结果将为GRAPES_TYM 3.0 版在业务预报中更好地发挥作用提供重要参考。

2 GRAPES_TYM 概况

GRAPES_TYM 是等经纬度格点模式，水平方向采用Arakawa-C 网格，垂直方向上采用地形追随高度坐标，模式层顶在35 km 左右。试验所采用的模式物理过程包括WSM6 微物理过程、Meso-SAS 对流参数化、YSU 边界层过程、Noah 陆面过程、莫宁-奥布霍夫相似理论表面层参数化以及Goddard 短波辐射和RRTM 长波辐射参数化。

为了改善模式初始场对热带气旋的描述能力，GRAPES_TYM 采用了初始涡旋强度调整技术，使模式初始场中涡旋的初始强度接近观测（麻素红等，2018）。

3 模式垂直分层加密及对热带气旋路径和强度预报的影响

3.1 垂直分层加密方案

模式中低层加密对热带气旋强度预报至关重要。为此在水平分辨率0.09°、模式层顶不变的条件下，采用相同的模式积分区域，将GRAPES_TYM的垂直分层由50 层提升到68 层，其中900 hPa 以下垂直层次由原来的10 层加密到19 层，模式最低层由原来的52.26 m 降低至20.68 m。图1 为50 层和68 层垂直分层dz 随高度的变化。

图1 模式垂直分层（红色为68 层，黑色为50 层，dz 为两个垂直层次间的距离）Fig. 1 Model vertical levels （Red is for 68 levels and black is for 50 levels； dz indicates the depth of each layer）

3.2 统计分析

为了评估垂直分层加密对GRAPES_TYM 热带气旋路径、强度预报的影响，利用2018 年1807—1826 号共计20 个热带气旋进行了批量试验。预报试验每天2 次 （00 和12 时，世界时，下同），总样本数为98 个，预报时效120 h。以下路径及强度误差计算均采用中国气象局上海台风研究所提供的最佳路径。

图2 为批量试验平均误差分布。图中包括了最大误差（×所示）、平均误差（曲线）、误差百分比分布（10%、25%、50%、75%、90%）。结果显示：垂直分层加密对72 h 前路径预报影响不明显，72—120 h 平均路径误差分别减小5.1%、8.4%和7.1%，120 h 最大路径预报误差由1943 km 减小到1304 km（图2a）。垂直分层加密对强度预报影响更为显著（图2b）：24—120 h 最大风速平均误差分别减小14.1%（24 h）、17.1%（48 h）、29.1%（72 h）、30.9%（96 h）和9.1%（120 h）；90%、75%、50%的样本中小误差样本明显增多，72 h 和96 h 的最大误差明显减小，其中72 h 最大误差由26.4 m/s 减小到14.4 m/s，96 h 的最大误差由19.9 m/s 减小到15.5 m/s。

由于模式垂直分层加密主要集中在模式中下层，可以改善模式对强台风的预报（Zhang，et al，2015），减小强度预报的负偏差，但强度预报正偏差有所增大（图3）。当垂直分层由50 层增加到68 层时，24、48、120 h 强度预报正偏差最大值增大（图中×所示），72、96、120 h 最大负偏差减小。在24—120 h 预报区间，从90%和75%强度偏差分布来看，68 层模式的正偏差均大于50 层模式，而从25%和10%的偏差分布来看，68 层模式的负偏差均小于50 层模式。从模式积分开始，50 层模式的平均强度预报负偏差快速增大，到60 h 增加到最大（−5.23 m/s）。60—96 h 强 度 偏 差 稳 定 在−4.72—−5.23 m/s，96 h 后负偏差有所减小。

图2 平均路径误差 （a） 和平均绝对强度误差 （b） （× 代表最大误差）Fig. 2 Mean track errors （a） and Mean absolute intensity errors （b）（× indicates the maximum error）

图3 强度偏差 （单位：m/s，× 代表最大偏差）Fig. 3 Biases of intensity（unit：m/s，× indicates the maximum bias）

68 层模式预报的平均强度误差以及强度预报负偏差的减小主要源于其对强台风以及超强台风预报能力的提高。以48 h（图4a）和96 h（图4b）为例分析68 层和50 层模式强度预报误差的差值（68 层误差−50 层误差）随风速的变化，图中正值表示68 层模式的预报误差大于50 层模式，负值表示68 层模式的预报误差小于50 层模式。可以发现：无论是48 h 预报还是96 h 预报，除个别样本外，68 层模式对最大风速大于50 m/s 的超强台风强度预报误差小于50 层模式；对最大风速大于40 m/s的强台风，68 层模式的强度预报误差大部分要小于50 层模式；而对小于40 m/s 的热带气旋的强度预报，二者预报差异不明显。

由于68 层模式对强台风以及超强台风强度预报能力的提升，预报的风压关系也明显好于50 层模式（图5）。 48 h 预报中，50 层模式预报的最大风速为55 m/s 左右，中心气压915 hPa 左右（图5a），而68 层模式预报的最大风速可以达到62 m/s 左右、中心气压910 hPa 左右（图5b）；96 h 预报中，50 层模式预报的最大风速为50 m/s左右，中心气压915 hPa 左右（图5c），68 层模式预报的最大风速可以达到60 m/s、中心气压910 hPa 左右（图5d）。68 层模式预报的风压关系更接近观测。

3.3 典型个例分析−超强台风“山竹”

2018 年9 月7 日12 时，台风“山竹”在西北太平洋生成，11 日08 时加强为超强台风，并且维持超强台风级别达96 h，最强风力超过17 级（65 m/s）。15 日，台风“山竹”在菲律宾北部登陆，16 日17 时在广东台山海宴镇登陆，登陆时中心附近最大风力14 级（45 m/s），中心最低气压955 hPa。

图6 为50 层及68 层模式对“山竹”的路径和强度预报。初始时刻为2018 年9 月8 日00 时到9 月12 日00 时，预报时效为120 h，初时时间间隔为24 h。从预报路径（图6）来看，50 层和68 层模式的结果比较接近，而平均路径预报误差显示，68 层模式24 和120 h 的结果略大于50 层模式，48、72 和96 h 的结果略小于50 层模式（图略）。除了72 h 误差差值超过30 km外，其他时刻平均误差差异均小于10 km。从10 m最大风速来看，50 层模式所有时次预报的极值在50 m/s上下，相比于观测最大值65 m/s 显著偏低；68 层模式预报的最大风速可以超过60 m/s，更接近观测。68 层模式24—120 h的平均强度误差均小于50 层模式，技巧最低为9.7%（48 h）、最高可达25.7%（96 h）（图略）。由此可见，垂直分辨率提高可以显著减小GRAPES_TYM对“山竹”的强度极值的预报误差，从而减小强度预报的平均误差。

图4 48 h （a） 和96 h （b）预报68 层和50 层模式最大风速误差绝对值差（蓝色点）（红色虚线为平均值）Fig. 4 Differences （errors of L68 minus errors of L50） in absolute error of 10 m maximum wind between L68 and L50 model （blue dots） at （a） 48 h and （b） 96 h forecast （The red dashed lines show averaged values）

图5 模式预报的风压关系 （a. 50 层模式的48 h 预报， b. 68 层模式的48 h 预报， c. 50 层模式的96 h 预报， d. 68 层模式的96 h 预报；红色是观测数据，蓝色是预报数据）Fig. 5 Relationship between Vmax and central pressure （a. 48 h for L50 model，b. 48 h for L68 model， c. 96 h for L50 model， d. 96 h for L68 model； red dots and lines are observed central pressure and Vmax，blue dots and lines are forecasted central pressure and Vmax）

4 预报区域扩大对热带气旋路径及强度预报的影响

采用水平分辨率0.09°、垂直分层为68 层的GRAPES_TYM 模式，在相同物理过程条件下，对比分析模式预报区域扩大对热带气旋路径及强度预报的影响。将覆盖西北太平洋、中国南海（0°—50.25°N，90°—171°E）的预报区域标注为D1，覆盖西北太平洋、中国南海及北印度洋（15°S—60°N，40°—180°E）的预报区域标注为D2（图7）。D2 包括了整个青藏高原，地形更加复杂。

预报区域扩大到D2，除了包含更复杂的地形外，将会有更多的天气系统进入区域模式预报范围。图8 为2018 年9 月7 日12 时D2 的初始场，图中红色矩形区域为D1。可以看到D1 预报区域向北扩展到60°N、向西扩展到40°E 时，将有更大范围的西风槽系统和大陆高压进入预报范围；预报区域南扩到15°S，会将位于40°—60°E 和110°—120°E的两支越赤道气流以及赤道辐合带纳入预报范围。一方面区域模式对上述天气系统的预报能力会影响其对热带气旋路径及强度的预报，另一方面预报区域扩大会减小侧边界对位于边界附近热带气旋预报的影响，如图中所示2018 年超强台风“山竹”，9 月7 日12 时位于（12.9°N，165.3°E），而D1的东边界为170°E，两者仅相距不到5°。

利用2016—2018 年生命期超过3 d 的热带气旋对预报区域扩大进行试验，预报时效为120 h，预报间隔为12 h（00 和12 时）。

试验结果统计分析显示，区域扩大对西北太平洋及中国南海热带气旋平均路径误差及强度预报误差无明显影响，0—120 h 平均路径预报误差改变未超过2%、强度预报平均误差改变未超过5%（图略），平均路径预报误差和强度预报误差改变最大值均出现在120 h，说明预报区域扩大对预报后期影响相对明显。

图6 2018 年超强台风“山竹”路径 （a. 50 层，b. 68 层） 及10 m 最大风速 （c. 50 层，d. 68 层）（黑色线代表观测路径，其他彩色线代表不同初始时间的预报）Fig. 6 Forecast tracks of super typhoon Mangkhut （a. 50 levels，b. 68 levels） and 10 m maximum wind （c. 50 levels，d. 68 levels）（Black line：best track data， colored lines：forecasts with different initial time as indicated in （c） and （d），such as the red is the forecast initiated at 00：00 UTC 8 September 2018）

图7 积分区域D1 和D2Fig. 7 Model domains of D1 and D2

D2 和D1 区域预报路径误差和强度误差相对较小的样本数比较接近：24—120 h 路径预报误差D2 小于D1 的样本数百分比分别是47.8%、51.33%、51.3%、47%和43.7%。强度预报误差D2 小于D1的样本数百分比分别是50.2%、52.3%、48.3%、51.6%和51%。由此可见，预报区域扩大对路径预报负影响略多，对强度预报正影响略多。

从120 h 路径预报误差地理分布可以看出，D1 和D2 预报误差较大的区域均在20°N 以北，这个区域有较多的转向热带气旋（图9a、b）。而预报误差差值（D2 的路径误差-D1 的路径误差）较大的区域和D1、D2 大误差区域相近（图9c）。预报区域扩大对20°N 以北的热带气旋路径预报影响相对明显，对20°N 以南的预报路径影响不明显。20°N 以北的转向热带气旋，影响其移动路径的天气系统比较复杂，除了副热带高压的影响外，还包括西风槽等天气系统的影响。模式预报区域越大，其所覆盖的天气系统越多，由区域模式自身特征所导致的不确定性就会增大，大区域和小区域的路径预报就会存在较大的差异。

5 与NCEP-GFS 和ECMWF 台 风 路 径 及强度预报对比分析

图8 2018 年9 月7 日12 时500 hPa 位势高度（台风符号为超强台风“山竹”的位置；红色矩形所示区域为D1，外面黑色区域为D2）Fig. 8 Geopotential height at 500 hPa at 12:00 UTC 7 September 2018 （The TC symbol shows the location of the super typhoon Mongkhut； The red rectangle shows D1 and the outer domain is D2）

图9 D1 区域 （a） 和 D2 区域 （b） 120 h 预报台风路径误差地理分布以及（c） D1 和 D2 预报的路径误差差值 （D2−D1） 的地理分布 （单位：km；Min：最小值，Max：最大值，ME：平均值） Fig. 9 Distributions of track error from D1 （a），D2 （b）and track error difference （D2−D1） （unit：km；Min：Minimum，Max：Maximum， ME：Mean）

将GRAPES_TYM 3.0 版回算的2016—2018年结果与美国NCEP-GFS 和欧洲中期天气预报中心全球模式ECMWF 的预报结果进行同样本比较（图10），24—120 h 的样本数分别为576、499、393、296、215 和150。

图10 平均路径误差（a）、平均强度误差 （b）和强度预报偏差 （c）Fig. 10 Mean track errors （a），mean intensity errors （b） and biases of forecast intensity（c）

试验结果显示，ECMWF 的平均路径误差显著小于NCEP-GFS 和GRAPES_TYM，而NCEP-GFS和GRAPES_TYM 的平均路径误差比较接近（图10a），其中ECMWF 24—120 h 的平均路径误差分别为56.23、101.51、170.2、246.97 和339 km，NCEP-GFS预报的结果分别为60.97、114.13、197.94、315.58和439.58 km，GRAPES_TYM 预报的结果分别为70.13、128.65、207、322.72 和450.54 km。GRAPES_TYM 24—96 h 平均强度预报的平均误差明显小于NCEP-GFS 和ECMWF（图10b），120 h 平均强度预报误差大于ECMWF 的预报，小于NCEP-GFS 的预报，其中ECMWF 预报的结果分别为7.0、8.46、8.47、8.16 和7.15 m/s，NCEP-GFS 分别为6.13、6.95、7.75、8.25 和8.98 m/s，GRAPES_TYM 分别为5.57、6.27、6.47、6.8 和8.08 m/s。ECMWF 和NCEP-GFS 强度预报误差偏大的一个主要原因是两个全球模式均未对分析场中的涡旋强度进行调整使其接近观测强度。从强度预报偏差来看，相比于GRAPES_TYM， ECMWF 模 式24— 120 h、 NCEP-GFS 的24—72 h 强度预报存在明显负偏差即强度预报偏弱（图10c）：ECMWF 的24—120 h 强度预报负偏差分 别 为−5.47、−6.68、−6.31、−5.31 和−3.41 m/s，NCEP-GFS 的 24—72 h 强 度 预 报 负 偏 差 分 别 为−2.91、−2.83 和−1.43 m/s。

路径预报也就是热带气旋中心位置预报，即中心的经度和纬度预报。ECMWF 预报路径的经度（图11a）存在负偏差，即相对于观测路径位置偏西，但24—120 h 最大偏差均未超过0.5°；NCEP-GFS 和GRAPES_TYM 预报路径的经度存在明显的正偏差，即相对于观测路径位置偏东，尤其是在96—120 h，NCEP-GFS 模式的120 h 预报路径经度偏差可达1.65°、GRAPES_TYM 模式的经度偏差可达1.53°。ECMWF 预报的24—120 h 纬度偏差较小，且均为正偏差，即预报路径存在微小的北偏；NCEP-GFS预报路径北偏趋势比ECMWF 大，在120 h 可达0.38°；相比于上述2 个全球模式，GRAPES_TYM预报路径在纬度方向的正偏差显著偏大，24—120 h台风中心位置纬度预报偏差分别为0.28°、0.45°、0.69°、0.97°和1.31°（图11b）。

图11 台风中心位置经度偏差 （a）和纬度偏差 （b）Fig. 11 Longitude deviation of typhoon center position （a） and latitude deviation of typhoon center position （b）

6 结 论

为了改进GRAPES_TYM 对西北太平洋及中国南海热带气旋的预报能力、为北印度洋热带气旋业务预报提供产品支持，2019 年8 月GRAPES_TYM 升级为3.0 版，包括水平分辨率由0.12°提升到0.09°，垂直分辨率由50 层增加到68 层，预报区域由（0°—50.25°N，90°—171°E）扩大到（15°S—60°N，40°—180°E）。扩大后的预报区域覆盖了西北太平洋、中国南海、北印度洋以及亚洲大部分区域。本文分析了模式垂直分辨率由50 层增加到68 层以及预报区域扩大对GRAPES_TYM 热带气旋路径及强度预报的影响。

结果显示，模式垂直分辨率提高对改善GRAPES_TYM 路径及强度预报至关重要。当模式垂直分辨率由50 层增加到68 层时，96—120 h 的平均路径误差可以减小8.4%—7.1%。提高模式垂直分辨率可以显著提高GRAPES_TYM 对强台风以及超强台风强度极值的预报能力，减小强度预报负偏差、改善模式预报的风压关系。相对于垂直分辨率50 层的预报，68 层预报的24—120 h 平均强度误差可减小13.2%、16.0%、29.1%、30.9%和9.1%。

模式预报区域扩大对平均路径和强度预报误差影响不显著。预报区域扩大后，大区域的路径及强度预报误差小于小区域的次数比率在50%附近（120 h 除外，为43.7%）。20°N 以北的预报路径对预报区域扩大更为敏感，这主要是由于20°N 以北的热带气旋影响其路径预报的天气系统比较复杂，除了副热带高压，还有西风槽、大陆高压等天气系统。模式预报区域扩大会带来更多的预报不确定性。

与NCEP-GFS 和ECMWF 模式对2016—2018年热带气旋路径及强度预报相比，GRAPES_TYM 3.0 版的平均路径预报误差与NCEP-GFS 接近，同ECMWF 相比误差较大；24—96 h 强度预报误差明显小于NCEP-GFS 和ECMWF，NCEP-GFS 和ECMWF对热带气旋强度预报存在明显的负偏差。