许冬梅 3,，沈菲菲 ，李泓，刘瑞霞，王易，束艾青

(1.南京信息工程大学 气象灾害教育部重点实验室/气候与环境变化国际合作联合实验室/气象灾害预报预警与评估协同创新中心，江苏 南京 210044；2.中国气象局上海台风研究所，上海 200030；3.中国气象局沈阳大气环境研究所，辽宁 沈阳 110166；4.高原与盆地暴雨旱涝灾害四川省重点实验室，四川 成都 610225；5.国家卫星气象中心 中国气象局中国遥感卫星辐射测量和定标重点开放实验室，北京 100081；6.江苏省气象台，江苏 南京 210008)

1 引言

西太平洋台风对我国的影响十分严重，平均每年有7～8 个台风登陆[1]，并带来暴雨、风暴潮灾害以及洪涝灾害等，造成重大的人员伤亡和财产损失。在过去的10 年，有多个超强台风对我国造成了诸如特大风暴潮、风灾、特大海潮、暴雨等多种自然灾害[2]。因此，如何有效地对台风进行预报对于防灾减灾至关重要。近年来，我国对如何提高台风路径的预报准确度做了大量的研究，统计了多种定位方法、路径预报方法、强度预报方法等，并发现数值模式对于改进台风的预报起了重要作用[3]。李秋阳等[4]研究了全球预报系统（GFS）再分析资料和日本气象厅（JMA）区域客观分析资料初始场对台风“桑美”的内部结构和移动路径的影响，发现二者均能较好地把握台风整体的流场形式，且JMA 资料在700 hPa 风场、近地面风速、降水等多个方面有更好的表现。由于台风的生命史中绝大多数时间都在宽阔的大洋面上，且其结构变化受到能量输送、中尺度强对流活动及其分布以及中纬度环流系统变形等的巨大影响。为提高数值模式的准确度，有必要获得较为精确的模式初始场，而精确的模式初始场的获得又依赖于资料同化。然而常规大洋观测系统，如浮标和船舶等，十分稀少，难以满足提升初始场的数据需求。同时，岸基雷达的有效观测半径有限，无法观测到远在大洋深处的台风。一系列的研究通过使用各种各样先进的卫星遥感仪器，已经改进了预报的精度。例如，Miyoshi 和Kunii[5]运用自适应膨胀方法对大气红外探测器（AIRS）数据进行增强，改善了大气温度和湿度数据在同化中的应用。Zou 等[6]利用GSI 分析系统，将GOES-13 和GOES-15成像仪的红外波段辐射强度观测数据进行直接同化，并在飓风天气研究预报模型（HWRF）对称涡旋初始化的基础上，增加了一个非对称分量，改善了墨西哥湾的热带气旋“德比”的轨迹和强度预报能力。俞兆文等[7]基于WRF-ARW 和WRFDA3.8 同化了AMSR2卫星微波成像资料，通过质量控制和偏差订正实验，改进了对2014 年9 号台风“威马逊”的预报。

卫星观测资料已成为数值天气预报中常规观测资料的补充，对改善模式初值条件和提高预报精度有很大的作用。地球同步卫星由于相对于地球表面静止且距离地球表面距离较远，可以对地球表面进行大面积连续监测。国内外科学家已经对静止卫星的应用开展了一系列的研究。耿晓雯等[8]在WRFDAv3.91中搭建先进的静止卫星扫描辐射成像仪AGRI（Advanced Geosynchronous Radiation Imager）同化接口，分析FY-4A 卫星天顶角对AGRI 资料偏差订正的影响，发现卫星天顶角与偏差关系不明显，不必在预报因子中进行考虑，提升了AGRI 红外通道辐射率资料在中尺度模式中的应用效率。Wang 等[9]使用日本葵花8 号（Himawari-8）卫星的辐射率资料，对2016 年7 月18−21 日的华北特大暴雨个例做了关于AHI（Advanced Himawari Imager）水汽通道辐射率资料同化特点的研究，发现引入新一代气象观测遥感器AHI 资料能够有效改善模拟的湿度场、风场和预报场，并对同化前6 h 和强度超过100 mm 的降水预报改进明显。Honda 等[10]以2015 年超强台风“苏迪罗”为研究对象，同化了Himawari-8 卫星辐射率资料，发现10 min一次的高分辨率的Himawari-8 卫星辐射率资料同化可以很好地改进核心区和外部雨带地区的热带气旋结构分析，并由此提高了对台风“苏迪罗”强度和路径预报的精确度。事实上，地球同步卫星可以清晰地观测到和热带气旋有关的旋转的对流性云并在分析热带气旋的最佳观测路径中起到重要作用。然而，云的信号影响了遥感信息的质量，降低了数据的利用率[11]。近年来，国内外专家学者在云检测方面做了卓有成效的工作。Goldberg 等[12]提出了适用于AIRS 的NESDIS-Goldberg 云检测方案，这是一种将AIRS 通道和微波通道进行经验组合，改善云检测能力的方式。官莉[13]研究了用空间匹配的高空间分辨率的成像仪的云检测信息以及AIRS 探测器序列定位误差对其观测值、晴空订正辐射强度和云特性反演的影响。官莉和王振会[14]通过运用与AIRS 空间匹配的MODIS（Moderate-Resolution Imaging Spectro-radiometer）产品确定受云污染的数据。朱文刚[11]结合GPAPES-3DVAR系统和AIRS，建立了适合于GRAPES-3DVAR 模式的晴空通道云检测方案，大大增加了资料，尤其是中高层资料的使用率，并改善了中高层的分析场。然而，云解析模型的计算需要大量资源，包括高度非线性的过程[15]，同时可能会受到参数化方案选择的影响[16]。因而，剔除云信息，仅保留晴空资料进行同化是目前解决云污染的有效方法。

静止气象卫星位于较高的地球同步轨道，因此视野广阔，可观测范围约1.7 亿km2，占地球表面的1/3。同时，它可以连续观测同一地区，一般可提供30 min一次的全盘扫描资料，或是3～5 min 提供一次特定地区的观测资料，并连续监视天气系统的演变。因此，它在对台风预报方面具有明显优势。2015 年7 月，日本气象厅开始正式运行葵花8 号地球同步卫星[16]，这是一系列第三代地球同步气象卫星中的第一颗。它成功发射于2014 年10 月7 日，并在试运行9 个月后开始业务观测。2015 年8 月，葵花8 号卫星成功地获得了台风“苏迪罗”快速增强的过程[10]。然而值得指出的是，Honda 等[10]是针对超强台风“苏迪罗”在海上的快速加强阶段进行同化研究，同化窗区间的下垫面为均一的开阔洋面。而适用于我国近海快速增强台风的且针对葵花8 号辐射率资料同化研究工作还开展较少。另外Honda 等[10]采用的是计算代价成本相对较高的集合同化方法，对我国近海快速增强的台风而言，为了赢得宝贵的预警时效，有没有可能采用一种基于快速高效的GPAPES-3DVAR 同化算法来分析葵花8 号辐射率资料并考察其对近海台风分析和预报的影响？这也恰恰是气象业务部门亟需解决的难题之一。鉴于现阶段葵花8 号资料在数值预报中的应用研究还不多，本文以2017 年13 号超强台风“天鸽”为例，同化葵花8 号地球静止卫星的红外辐射率资料，通过改进初始场，讨论其对台风路径、强度的模拟和预报能力的提升作用。

2 观测资料及同化系统

2.1 Himawari-8 晴空红外辐射率资料介绍

葵花8 号卫星搭载了先进的静止轨道成像仪AHI，它具有16 个观测通道（表1），其8、9、10（波长为6.2 μm、6.9 μm 和7.3 μm）3 个通道能够对水汽进行探测，因此有较强的对云探测的能力。表1 给出了各个通道透过率随着高度的变化，即各个通道的权重函数（dTau/dlnP）。AHI 能提供10 min 一次的16 通道的全圆盘高分辨率的观测资料，其高时间分辨率的资料能够有效改进资料同化的准确性[17]。此外，其空间分辨率也大幅提高，达到了1 km左右。已有研究揭示了葵花8 号卫星观测数据的特点以及其对降水预测的作用[18–19]。

表1 AHI 16 个通道的中心波长Table 1 The central wavelength of 16 AHI bands

本实验采用通道8、9、10 这3 个水汽通道，它们是红外波段，中心波长如表1 所示，分别为6.2 μm、6.9 μm和7.3 μm。由于本文同化的观测资料主要是针对水汽通道的AHI 辐射率资料，首先采用晴空条件下的模式温度廓线和湿度廓线作为大气辐射传输模式（Community Radiative Transfer Model，CRTM）来输入计算权重函数（图1）。由图1 可以发现葵花8 号3 个水汽通道的权重函数峰值位于对流层的中高层位置，从而为对流层的AHI 辐射率资料同化提供丰富的大气水汽信息。水汽通道8、9、10 的权重函数峰值位于300～400 hPa 附近的对流层中高层，其中通道10 较通道8 和通道9 要低，通道8 峰值最高，因此减小了地面带来的影响。在全空同化红外辐射率资料过程中，由于辐射率资料和模式变量间复杂的非线性以及观测误差和系统偏差不确定性较大等问题，本实验只同化水汽通道的晴空辐射率资料。

图1 通道8、9、10 的权重函数Fig.1 The weight function of bands 8,9 and 10

2.2 WRFDA 系统及AHI 同化模块构建

2.2.1 WRFDA 同化系统

WRF（Weather Research and Forecasting Model）模式是一种可压缩、非静力的中尺度天气预报模式。而WRFDA 是WRF 模式的同化系统，由美国国家大气研究中心开发和维护。本研究基于其中的3DVAR 同化方法，通过分析多源观测资料来改善初始场[20–22]。

2.2.2 AHI 同化模块

WRFDA 系统目前可以同化多种常规观测和非常规观测，本研究采用系统自带的CRTM（Community Radiative Transfer Model）作为观测算子来对AHI 辐射率资料进行模拟计算。为了去除受云和降水污染的格点，引入了非常严格的云检测方案[22]。具体步骤包括：通过通道14 亮温的空间和时间变化、通道14 和通道15 的相关差异、通道14 和水汽通道10 之间的空间相关以及通道7 的辐射率值来检测云。

2.2.3 偏差订正模块

卫星资料同化的过程中包含了来自于观测本身和模式以及同化系统的误差。本文利用WRFDA 系统中的回归偏差订正法来修正卫星观测算子，公式为

被修正的观测算子H˜ 在以往的观测算子H(x)基础上额外加入了总偏差常数项β0和各预报因子pi及其系数项 βi。在本研究中，首先通过离线的方式统计得到第一个同化时次所需的初始偏差订正系数，并在同化的过程中得以更新。后续的同化时次用到的偏差订正系数为前一个同化时刻更新过的系数。偏差订正使用到的预报因子包括：1 000～300 hPa 和200～50 hPa 层结厚度、地表温度以及水汽总含量。

2.3 台风“天鸽”个例

台风“天鸽”是2017 年13 号台风，于2017 年8 月20 日14 时在20.4°N，118°E 处生成，中心附近最大风力为8 级，中心最低气压为1 000 hPa。之后“天鸽”的风力逐渐加强并向西移动，其观测路径如图2a 所示。值得注意的是，“天鸽”在8 月22 日08 时和15 时风速连升两级，迅速升级为台风，并向西偏北方向移动。此后，“天鸽”的风速继续增强，移速不断增大，到8 月23 日05 时，“天鸽”移动至21.3°N，115.2°E，中心附近最大风力加强到13 级，中心气压为955 hPa；并在2 h 后进一步加强到14 级，中心气压为950 hPa；到8 月23 日10 时，其中心附近最大风速达到了15 级，中心气压为945 hPa。“天鸽”能够在5 h 内迅速增强3 级，可能是近海浅水区域（深度<30 m）反常温暖的结果。8 月23 日13 时，台风“天鸽”登陆，其登陆地点为22°N，113.2°E，登陆时等级为强台风级，并在24 日20 时被停止编号。台风“天鸽”具有以下4 个特点：第一，它的移动速度快，其预报移速一度达到30 km/h，导致应急准备严重不足；第二，它的强度变化大，在8 月23 日05 时到10 时的5 h 内中心附近最大风力由13 级迅速增加到15 级；第三，“天鸽”登陆时期恰逢天文大潮，使得风暴潮增水严重；第四，它的影响面广、影响大，按照美国联合台风预警中心（JTWC）的统计，台风“天鸽”严重摧毁了城市，造成了10 人死亡，仅澳门的经济损失估计超过15 亿美元。事实上，台风“天鸽”总共造成24 人死亡[23]，并对港珠澳一带造成严重的影响。因此，台风委员会宣布将“天鸽”除名。图2b 是8 月22 日18 时葵花8 号卫星13 通道（10.4 μm）的亮度温度图，从图中可以看到在“天鸽”中心附近有强烈的冷中心，其最低亮温可达190 K，图2b 也较好地反映出了台风螺旋雨带的结构。

图2 台风“天鸽”的观测路径（a）以及2017 年8 月22 日18 时Himawari-8 AHI 通道13 的亮度温度（b）Fig.2 The observation track of Typhoon Hato (a) and the brightness temperature of band 13 of the Himawari-8 AHI at 18:00 UTC on August 22,2017 (b)

2.4 实验设计

本实验采用WRF 模式对台风“天鸽”进行模拟。如图3 所示，模拟的中心区域为20°N，115°E。并采用两层网格嵌套方案，内外层网格格点数分别为401×329、481×430，模式顶层为10 hPa，垂直层数共61层。本实验采用了WRF 单参6 类微物理方案（WRF Single-Moment 6-class scheme）、RRTM（Rapid Radiative Transfer Model）长波辐射方案、Dudhia 短波辐射方案、MM5 similarity 近地面层方案、Noah Land Surface Model 陆面方案、Yonsei University 行星边界层方案以及Kain-Fritsch 积云参数化方案。

图3 WRF 模拟区域Fig.3 The simulation area of WRF

本次实验流程如图4，为了说明AHI 辐射率资料对同化的改进作用，设置了控制实验（CTNL，下同）和AHI 资料同化实验（AHI_DA，下同）。2017 年8 月22 日00 时（世界时）的初始场和侧边界条件由美国国家环境预测中心（National Centers for Environmental Prediction,NCEP）的0.25°×0.25°的GFS（Global Forecasting System)再分析资料提供。经过12 h 的Spinup 后得到第一次同化所需的背景场。并以此进行了连续12 h（每小时一次）的循环同化及预报实验。其中，CTNL实验采用常规观测资料进行同化实验，而AHI_DA 不仅采用常规观测资料，同时还引入AHI 观测资料进行同化。实验结果（表2）将在第3 部分进行分析。

图4 实验流程设置Fig.4 The experimental flowchart

表2 模拟实验方案设置Table 2 Simulation experiment scheme setting

3 实验结果分析

3.1 分析效果检验

3.1.1 预报因子分析

图5 给出了晴空条件下观测亮温减去模拟亮温（OMB）和各个预报因子分布的散点图，图中的阴影的颜色对应着分布密度。同时图5 列出了相应的偏差bias 和预报因子P的相关函数。通过图5 中的统计结果可以得到4 个预报因子偏差订正系数。该系数给出了不同预报因子相关的偏差系数，利用这些系数以及每次同化时背景场中的预报因子可以得到初始的偏差大小。基于初始的偏差经过同化过程进一步进行偏差订正。

图5 OMB 和预报因子50～200 hPa 层结厚度（a）、300～1 000 hPa 层结厚度（b）、总水汽含量（c）和地表温度（d）的分布散点图Fig.5 Scatter plots of OMB and predictors of 50～200 hPa layer thickness (a),300～1 000 hPa layer thickness (b),total moisture content (c),and surface temperature (d)

图6 为2017 年8 月22 日12 时AHI资料通道8、9、10 亮温观测残差频数分布直方图。图6a 黑色直方图为偏差订正前的OMB（观测亮温减去模拟亮温）频数分布。白色直方图为偏差订正处理后的结果，而灰色直方图为对分析场进行偏差订正后的OMA（观测亮温减去分析场亮温）频数分布。以通道8（图6a）为例，相较于图6a 中OMB 频数分布，图6a 中的OMA频数分布更加集中在0 K 附近，峰值位于0 K，频数达到3 500，而偏差订正前后的OMB 频数峰值分别为2 600 和2 500，且相较于OMB 更加接近于正态分布，说明分析场更加接近观测，同化效果较好。而通道9、10（图6b、图6c）也有类似的效果，且峰值频数更大，说明通道9、10 的分析场更加接近于观测场。

图6 2017 年8 月22 日12 时AHI 资料通道8（a）、通道9（b）、通道10（c）观测与偏差订正前背景场模拟的辐射率值差异分布直方图（黑色）、观测与偏差订正后背景场模拟的辐射率值差异分布直方图（白色）以及观测与分析场模拟的辐射率值差异分布直方图（灰色）Fig.6 Histogram of AHI observation minus simulated radiance from the background without bias correction (black),observation minus simulated radiance from the background with bias correction (white),and observation minus simulated radiance from the analysis (gray) for band 8 (a),band 9 (b) and band 10 (c) at 12:00 UTC on August 22th,2017

图7 为2017 年8 月22 日12 时到8 月23 日00 时每小时一次（共13 次）的亮温偏差时间序列，图7a 至图7c 依次为通道8、9、10，从上到下分别为观测数目、亮温偏差平均值和亮温偏差标准差。以通道8（图7a）为例，图中第一行为观测数目随时间的变化，整体呈减少趋势，由8 月22 日12 时的10 500 最终减少到8 月23 日00 时的7 500，这是进行了严格的质量控制的结果，云和降水点被严格剔除掉。图中第2 行为亮温偏差平均值随时间的变化，可见在偏差订正前，13 次同化的OMB 平均值大致在0.8 K 附近浮动，其中最大值为1.1 K，最小值为0.35 K，而偏差订正后，13 次同化的OMB 平均值位于0 K 附近，最大值0.35 K，最小值–0.4 K，说明偏差订正使得背景场模拟亮温更加接近于观测亮温，证明了偏差订正的有效性。而相对于OMB，OMA 的平均偏差基本同x轴重合，其平均偏差最大值为0.1 K，说明分析场基本同观测场相同，AHI 辐射率资料同化能够显著改善初始场分析。图中第3 行为亮温偏差标准差随时间的变化。偏差订正后和偏差订正前的OMB 标准差曲线重合。其中OMB 的亮温偏差标准差远大于OMA的亮温偏差标准差，除第一次同化外，其余12 次同化结果均位于1.2 K 附近，而OMA 的亮温偏差标准差在13 次同化实验中均位于0.4 K 附近。通道9（图7b）和通道10（图7c）均有类似的效果，说明AHI 辐射率资料可以显著改善亮温分析。值得注意的是，3 个通道的OMB 标准差均在8 月22 日12 时达到最大（图7a为2.0 K，图7b 为1.75 K，图7c 为1.35 K），而之后迅速减小，而OMA 标准差则相对较小，说明其改进模式背景场效果最好。

图7 2017 年8 月22 日12 时到8 月23 日00 时通道8（a）、通道9（b）、通道10（c）逐小时（共13 次）的观测数目（上），偏差平均（中）和偏差的标准偏差（下）Fig.7 Time series of hourly data counts (top),the averaged bias (middle) and the stdv of the bias (bottom) for band 8 (a),band 9 (b),band 10 (c) from 12:00 UTC August 22,2017 to 00:00 UTC August 23,2017 for total 13 cycles

3.1.2 风场和高度场分析

图8 是2017 年8 月22 日12 时00 分CTNL 实 验和AHI_DA 实验分析的500 hPa 高度场。图8a 为CTNL实验，图8b 为AHI_DA 实验。在图8a 中，闭合等压线的数目为4，最低等高线为5 780 位势米，位于20°～21°N，116°～117°E，与位于19.5°N，117.5°E 的台风中心偏离，未能模拟出台风“天鸽”的涡旋中心。而图8b中，AHI_DA 实验模拟出的闭合等压线数目增加到7，最低等高线为5 720 位势米，且其位置和台风中心重合。由图可见，AHI_DA 实验不仅较好地模拟出了台风“天鸽”的涡旋环流形势，而且其模拟出的台风涡旋强度也比CTNL 实验要强很多。

图8 2017 年8 月22 日12 时00 分500 hPa 位势高度场对比Fig.8 500 hPa geopotential height field at 12:00 UTC on August 22,2017

图9 是2017 年8 月22 日12 时00 分CTNL 实 验和AHI_DA 实验分析的海平面气压场和近地面风场的合成图。图9a 为CTNL 实验，图9b 为AHI_DA 实验。此时，台风的中心位于19.5°N，117.5°E 的海面。由图9a 可见，CTNL 实验的闭合等压线数目为3，且彼此间距离较大，海平面气压的最低值为990 hPa，说明CTNL 实验模拟的台风比较弱。在中心北偏西方向有风速的最大值，约为32 m/s，台风的风场较气压场向东北方向偏移。而AHI_DA 实验的闭合等压线数目增加到5，且等压线密度更大，反映出同化了AHI 辐射率资料后台风中心气压快速下降的特点，海平面气压最低值为978 hPa，气旋强度得到了显著增强。而且台风的大风带也由图9a 中位于台风中心北偏西方向变为图9b 中基本环绕台风中心，风速最大值位于中心偏南，并提高到了36 m/s，且台风风场与气压场基本匹配。由图可见，同化过程使得气旋强度和风速都得到了有效的增强，且使得台风的风场和气压场更加匹配。

图9 2017 年8 月22 日12 时00 分海平面气压场风场对比Fig.9 Sea level pressure and wind at 12:00 UTC on August 22,2017

3.1.3 均方根误差（RMSE）分析

在长达12 h 的同化区间内，分析场与观测场之间的误差一般可以用来评定同化系统的优劣。本文将同化区间的分析场与常规观测资料作检验。通过均方根误差平均值的廓线分布可以看出分析场和观测的接近程度。图10a、图10b 分别为同化区间的均方根误差平均值的垂直廓线分布图，由图可以发现两组分析的速度分量u和v在100 hPa 高度均方根误差最大，但是整体看来，AHI_DA 实验分析的风速的RMSE要显著小于CTNL 实验分析的RMSE。这也进一步表明通过同化AHI 辐射率资料可以通过多变量相关显著改进风场预报。图10c 表明对于温度场而言，AHI_DA 实验分析的温度场在1 000～400 hPa 之间相比CTNL 实验有显著的改进作用。通过图10d 可以进一步发现AHI_DA 实验在925～300 hPa 之间分析的RMSE 要显著低于CTNL 实验分析的湿度场RMSE，并且随着高度的降低，改进效果更加明显。表明AHI 辐射率资料中的水汽信息可以有效延伸到对流层较低层，并对背景场起到一定的修正作用。值得指出的是，湿度场在925～1 000 hPa 表现不好可能是由于地表发射率存在较大的误差以及本文所同化的通道的权重函数位于对流层中上层等原因造成的。

图10 同化分析场相对于常规观测资料的RMSE 均值的垂直廓线Fig.10 The averaged RMSE of the analyses against the conventional data

3.2 预报效果检验

3.2.1 降水效果

图11 为8 月22 日18 时 到8 月23 日00 时 的6 h累积降水，此时，台风尚未登陆。由图11a 中可以看出，在这6h内，主要降水区域集中在22°～26°N，114°～116°E 范围的沿海区域，还有一降水最大值区域位于111°E 以西的内陆区域。最大降水强度中心位于24°N，115°E 附近，6 h 累积降水量接近50 mm。然而，由图11b 可见，CTNL 实验预报得到的主要降水区域为112°～118°E 的沿海及内陆区域，且最大降水中心降水强度超过120 mm，而位于111°E 以西的降水区域则被漏报。图11c 为AHI_DA 实验预报得到的6 h 累积降水，主要降水区域为22°～31°N，111°～119°E 的沿海区域，降水区域的面积和降水强度均有所减小。可见AHI_DA 实验得到的6 h 累积降水图更加接近于实际的观测情况。

图11 2017 年8 月22 日18 时到23 日00 时6 h 累积降水Fig.11 Accumulated 6 hour precipitation initialized from 18:00 UTC on August 22,2017

本文进一步对各组实验6 h 累计降水量计算FSS（Fraction Skill Score）[24]评分。定义FSS 评分为

式中，Pfi和Poi分别是第i个网格点处超过给定阈值的预报和观测降雨区域的概率；N是验证区域中的网格点总数。FSS 的值在0 到1 之间，0 表示预报和观测的降雨区域没有重叠，1 表示预报和观测的降雨区域完全重叠。由图12 可以发现，资料同化实验明显优于控制实验，得到的FSS 显著提高。对于小于25 mm阈值的评分更加优于25 mm 阈值的评分。

图12 各量级6 h 累计降水量的FSS 评分Fig.12 FSS scores of accumulated 6 hour precipitation for different thresholds

3.2.2 路径和强度预报

为了评估AHI 辐射率资料同化对台风“天鸽”预报的影响，本文在同化后作了2017 年8 月23 日00 时至24 日00 时的24 h 确定性预报。图13a是台风“天鸽”的路径预报效果图。在预报起始的6 h 内两组实验预报的路径误差均相对较小，但是在起始预报的6 h之后，CTNL 实验预报的台风路径逐渐向偏南方向转折，并逐渐产生了较大的误差。而AHI_DA 实验虽然在起始预报的12～18 h 内也有部分向南转折的趋势，但在起始预报的18～24 h 内，误差又逐渐减小。总体而言，在24 h 的确定性预报中，相比于CTNL 实验，AHI_DA 实验预报的路径更加接近于实际观测路径。图13b 是两组实验预报的台风路径误差图。从图13b 中可以发现，初始时刻，同化实验和控制实验的预报误差分别为33 km 和42 km，而经过模式6 h 预报之后分别为33 km 和37 km，此时，两组预报实验的路径误差均相对较小。之后，控制实验的误差迅速增加，到起始预报时刻24 h 后，控制实验的路径误差达到最大，为123 km。而同化实验在起始预报的12～18 h 内也有迅速增加的趋势，在起始预报18 h 后，误差也达到了70 km；但在起始预报18 h 后，又逐渐下降，到起始预报24 h 后，路径误差下降到56 km。综合图13a 和图13b 来看，在24 h 的预报过程内，同化实验预报的路径误差均明显小于控制实验预报的路径误差。两组实验预报的台风路径均相对于最佳观测路径偏南，且随时间推移，预报的路径会更加偏南，路径误差总体也在逐步增加，24 h 内同化实验和控制实验预报的平均路径误差分别达到了约45 km 和67 km。总体而言，AHI 同化实验预报的台风路径相对于控制实验改进显著。从前文可以看出，同化AHI 辐射率资料能一定程度上对台风的动力场、热力场、质量场有所改善，因而图13b 和图13d 给出了两组实验在台风中心最低气压和近中心最大风速方面的预报效果。从图13c 可以看出，AHI 同化实验的中心最低气压始终和最佳观测更为接近，控制实验预报的台风强度偏弱。图13d 中两组实验对台风近中心最大风速的预报结果和气压场结论较为一致，控制实验预报的台风近中心最大风速偏弱，而AHI 同化实验和最佳观测的台风近中心最大风速误差始终小于2 m/s。

图13 2017 年8 月23 日00 时后24 h 确定性预报的实际路径和预报路径（a），路径误差随预报时间变化分布（b），最低海平面气压随预报时间变化分布（c）和最大风速随预报时间变化分布（d）Fig.13 24-hour deterministic forecast at 00:00 UTC on August 23,2017,the best track and forecast track (a),the track error with the forecast leading time (b),minimum sea level pressure with the forecast leading time (c) and maximum wind speed with the forecast leading time (d)

4 结论

本文基于WRFDA 同化系统，引入Himawari-8 的AHI 红外辐射率资料，在同化过程中采用严格的云检测方式剔除云和降水的影响，通过设置AHI_DA 和CTNL 两个实验，来检验AHI 晴空红外辐射率资料同化是否对台风“天鸽”的分析预报有改进作用，并可以得到以下结论：

（1）引入AHI 红外辐射率资料能使模拟亮温的偏差减小，经过偏差订正后，OMB 分布基本接近均值为0 的正态分布，但相对OMB 来说，OMA 的偏差标准差和平均值均大幅下降，分布更加收敛，说明同化情况较好。

（2）AHI 晴空红外辐射率资料的引入能够有效地改进台风的海平面气压场，近地面风场和500 hPa高度场的预报，经过同化之后，台风的中心气压和最大风速都得到了改进，且台风风场和流场也更加匹配。

（3）相对于CTNL 实验，AHI_DA 实验得到的纬向风速、经向风、温度、湿度等要素的均方根误差均有所减小，更加符合实际情况。

（4）AHI_DA 实验能够改进对6 h 累计降水的预报，相对于CTNL 实验有效减少了预报的降水区域和降水强度，使之更加接近于观测值。

（5）AHI_DA 实验和CTNL 实验在前6 h 内对台风路径的预报效果均较好，但6 h 后AHI_DA 实验的预报效果明显优于CTNL 实验。

本文实现了Himawari-8 卫星红外辐射率资料同化实验，但是由于本文检验的台风个例只有1 个，且对台风累计降水和路径的预报仍有偏差。未来将采用更先进的同化方式，如四维变分同化，集合–变分同化等进一步针对更多的台风个例进行研究；同时针对集合同化方法中水平局地化距离的调整等方案展开细致深入的研究，最后，本文去除了云和降水区域的观测资料，只进行了AHI 资料的晴空同化，因此实现AHI 资料的全天空同化也是本研究参与者的下一个主要研究方向，同时也包括联合其他卫星遥感资料进行多源同化尝试[25]。鉴于本文尝试利用日本葵花8 号高时空分辨率卫星水汽通道数据开展同化实验，得到了较好的同化和预报效果，对未来在业务化预报进一步使用新一代静止卫星葵花8 号辐射率资料提供了重要参考，同时也为我国自主静止卫星资料的同化应用奠定了一定基础。由于静止卫星的时空连续性，下一步我们将开展长时间的同化实验，进一步探讨静止卫星在改善台风的路径和强度预报中的作用。

致谢：本论文的数值计算得到了南京信息工程大学高性能计算中心的计算支持与帮助，在此表示感谢！