基于FY-4A大气温度廓线数据的西北太平洋低云云顶高度反演研究

李顺，衣立，徐娜，时晓曚，张苏平，刘敬武

(1.中国海洋大学深海圈层与地球系统前沿科学中心，山东 青岛 266100；2.中国海洋大学海洋与大气学院，山东 青岛 266100；3.国家卫星气象中心，北京 100081；4.青岛市气象台，山东 青岛 266003)

引言

云的特性，包括云量、云高度、垂直剖面、光学厚度、液体和冰水含量以及颗粒大小，其影响地球的辐射收支，与云-气候反馈之间存在潜在联系，并且云通过调节地球的辐射收支、水循环来影响气候变化，不同高度云的特性会产生不同的气候影响[1-3]。其中低云的特性变化对气候的影响显著，当低云云量增量达到5%时，就可以抵消CO2翻倍造成的辐射强迫[4]；但是，低云在气候模拟中的不确定性导致模拟结果的差异很大[5]。西北太平洋是低云和海雾的频发区[6-8]，也包含了强烈的海洋锋区，海洋锋区通过海表面温度、风、涡动等影响西北太平洋的海气边界层，从而影响低云[9-11]。因此，作为区分低云与中高云的宏观属性，准确的云顶高度在研究中十分重要。另外，中国重要港口和沿岸机场均在西北太平洋沿岸，当低云云底接地形成大雾时[12-14]，会产生航班延误、港口停运等不利影响。可利用云顶高度与云层厚度之差来区分低云和雾，因此云顶高度也是区分低云与雾的重要参数[15]。

目前，海上站点观测不足以支撑大范围连续监测的需要，卫星遥感是获取海上云顶高度信息的重要手段。卫星遥感反演云顶高度的技术分为被动遥感技术和主被动结合遥感技术。主被动结合遥感技术可以利用激光雷达或毫米波雷达，得到相对准确的云顶高度，然而无论是星载雷达还是地基雷达，观测范围受时空条件的限制。地基雷达可观测局地云高的连续变化；相比地基雷达，星载雷达虽然可以获取更大范围的云顶高度，但是无法对局地云高进行连续观测[16-19]。被动遥感可以利用可见光成像仪，运用云体、卫星的相对位置等空间信息和视角差异、云阴影方法等几何方法，直接估计云顶高度，也可以利用卫星传感器被动接收的各个光谱通道信息，通过不同的反演算法，获取云顶高度[20-21]。MODIS的云顶高度产品采用CO2切片法和红外窗区法计算获取云顶高度[22]，CHENG et al.[23]利用云阴影法估计北极低云的云顶高度。然而，静止卫星反演的低云云顶高度无论是在统计结果还是个例表现上均与主动遥感探测结果有显著的差距[18-19]。MARCHAND et al.[24]比较低云云顶高度的反演结果时，发现不同卫星的海上低云(尤其是层积云)云顶高度产品之间出现很大的差异，是因为反演算法中使用天气模式预报的大气温度廓线无法准确地模拟层积云顶部的盖顶逆温，从而降低低云云顶高度的反演准确性。虽然使用固定的温度递减率来代替大气温度廓线可极大地改进局地平均的层积云云顶高度，但改变观测区域后，原有固定的温度递减率不再适用[25-27]。

此前研究和现场观测表明，海上低层层状云基本位于一个盖顶逆温层之下。NICHOLLS and LEIGHTON[28]在6次层积云的飞机观测中，得到层积云平均垂直特征，发现单个云层位于一个强烈而浅薄的逆温层下，稳定的层结隔断云层与下垫面或者埃克曼(Ekman)层之间的联系，由观测结果看，逆温层底的位置与云顶高度位置基本一致。BRETHERTON et al.[29]也发现逆温层底所在高度与附近层积云的云顶高度密切相关。出现这一现象的原因是云顶受到长波辐射以及环境条件的影响。海上低层层状云的云顶会发生强烈的长波辐射冷却，它可以维持并增强云顶逆温层，使其在垂直方向几十米内温度变化10～20 K[8]；同时，在一些强烈下沉区域上，空气由于下沉至近地层形成逆温层结，逆温层结可以抑制水汽的向上输送，因此这些低层层状云常在逆温层下出现[28,30-31]。基于低层层状云的云顶经常出现强烈逆温层这一现象，YI et al.[32]利用逆温层所在高度表征黄海低云和雾的云顶高度，其方法是使用卫星反演逆温层顶强度和厚度，结合黄海沿岸几个测站的平均探空廓线，来确定低云的云顶高度。为提高低云云顶高度的反演准确性，基于层积云顶部强烈逆温层的持续存在以及其底部和层积云云顶高度的密切相关性(图1)，利用在西北太平洋的大气温度廓线，提取到温度廓线上逆温层底高度，并将其作为低云的云顶高度，这种反演低云云顶高度的方法可称为“逆温法”。

图1 逆温层底反演低云云顶高度示意图Fig.1 Schematic diagram of retrieval of low cloud CTH (cloud top height) by temperature inversion layer base

相较于气球无线电探空，高光谱红外探测仪反演的温度廓线可以得到范围大、时间间隔规律的大气温度廓线，以往这类仪器多搭载于NOAA-19、FY-3C、MetOp-A等极轨卫星上，然而对同一地区的回访周期较长[33]。FY-4A/GIIRS(Geostationary Interferometric Infrared Sounder)是第一个搭载在静止卫星的高光谱红外探测仪，利用GIIRS观测的高光谱数据，得到西北太平洋区域时空分辨率较高的大气温度廓线产品[34-35]，与欧洲极轨卫星MetOp-A/IASI(Infrared Atmospheric Sounding Interferometer)的大气温度廓线比较可见，FY-4A/GIIRS的温度反演精度明显优于IASI[36]，因此基于FY-4A/GIIRS观测资料反演的大气温度廓线具有重要的应用价值。

选取2019年1月18日—12月31日的星载激光雷达CALIPSO(Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation)低云云顶高度(cloud top height，CTH)和相应时刻FY-4A/GIIRS/AVP(Atmosphere Vertical Profiles)数据，利用“逆温法”计算得到低云的云顶高度(命名为CTHGIIRS)，并与CALIPSO观测的云顶高度、FY-4A/AGRI(Advanced Geosynchronous Radiation Imager)/CTH云顶高度产品进行比较，评估本方法在反演低云云顶高度上的效果。同时参考ISCCP(International Satellite Cloud Climatology Project)确定的低云云顶气压应高于680 hPa(约3 200 m)，并考虑到CALIPSO对多层云下的低云云高探测可能存在偏差，低云个例均选取单层低云，即在低云上无其他云或者气溶胶的遮挡[37]。

1 数据与方法

1.1 卫星数据

本文主要使用FY-4A数据和CALIPSO数据，其中FY-4A卫星搭载多种先进的遥感仪器，包括多通道扫描成像辐射计(AGRI)和干涉式大气垂直探测仪(GIIRS)[35,38-40]。目前FY-4A已有的云顶高度产品为FY-4A/AGRI/CTH，反演使用的数据是FY-4A/GIIRS/AVP，检验反演结果和现有云高产品准确性的是CALIPSO的Level 2B产品大气垂直特征分类标识(vertical feature mask，VFM)——CALIPSO/VFM。

1.1.1 反演数据——FY-4A/GIIRS/AVP

FY-4A/GIIRS/AVP是FY-4A利用干涉式大气垂直探测仪数据反演得到的大气垂直廓线产品。AVP的反演过程使用FYGAT-S(Fengyun Geostationary Algorithm Testbed-Sounder)算法，集成4种大气辐射模型，应用统计回归检索算法和非线性物理检索算法，得到大气垂直廓线AVP产品[35]。AVP产品时间分辨率为55 min，空间分辨率为16 km×16 km，同时在680 hPa 以下，AVP垂直结构一共14层，平均每层之间高度约230 m。该产品主要反映大气中温、湿度和臭氧等痕量气体的三维分布和变化，同时具有较高的时间分辨率[41]，适合对低云的大范围快速探测。选取AVP产品中几个用于“逆温法”反演低云云顶高度的数据，其中最主要的数据包括大气温度廓线(atmosphere temperature profiles，AT_Prof)以及每一层对应的气压，同时需要考虑到AT_Prof的数据质量问题，选取了AT_Prof的质量标签(quality flag of atmosphere temperature profiles，AT_Prof_Qflag)，对已选取的逆温层进行质量控制，剔除质量不合格的数据[42-43](详见http://fy4.nsmc.org.cn/portal/cn/theme/FY4A_product.html)。

1.1.2 对比数据——FY-4A/AGRI/CTH

FY-4A/AGRI/CTH是FY-4A利用多通道扫描成像辐射计数据反演得到的云顶高度产品，简称CTHAGRI。CTHAGRI的反演使用FCTHA(Fengyun Cloud Top Height Algorithm)算法反演得到，该算法是FYGAT-I(Fengyun Geostationary Algorithm Testbed-Imager)系统的一部分[41,44-45]，主要算法是红外分裂窗区法和CO2切片法，得到云顶高度CTHAGRI，其时间分辨率为15 min，空间分辨率为4 km×4 km。

1.1.3 检验数据——CALIPSO/VFM

CALIPSO/VFM用于检验反演结果和现有云高产品准确性。其数据内容是使用Feature_Classification_Flag来指定某一高度的物体分类，给出8.2 km以下垂直分辨率30 m/层，8.2～20.2 km垂直分辨率为60 m/层的云和气溶胶等目标物的分类信息以及目标物对应的海拔高度[46-47]。Feature_Classification_Flag可显示某一高度上的目标物分类为云，称之为云标签。利用云标签的垂直分布，可排除受中高云和气溶胶干扰的单层低云。确定未受中高云和气溶胶影响的单层低云后，找到单层低云最高云标签的海拔高度，即可提取到单层低云的云顶高度[37]，后文简称为CTHVFM。

1.2 匹配方法

采用时空邻近点匹配比较它们的云顶高度。为确保选取的反演数据GIIRS/AVP和比较数据CTHAGRI在时间和空间上要尽可能地接近单层低云出现的CTHVFM，计算单层低云的CTHVFM与GIIRS/AVP产品、CTHAGRI之间的空间距离与时间间隔。当满足时间间隔在1 h以内、空间距离不超过0.5°时，取与CTHVFM时间间隔最小、空间距离最近的CTHAGRI与相应的CTHVFM进行比较；取与CTHVFM时间间隔最小、空间距离最近的AVP进行反演得到CTHGIIRS，并与相应的CTHVFM进行比较。

计算CTHVFM与GIIRS/AVP产品、CTHAGRI之间的时间间隔

Δt(i)=|tVFM-tm(i)|<1 h;i=1,2,3,…。

(1)

假设当i=p,Δt(p)=min(Δt)，则取p时刻的产品计算CTHVFM与GIIRS/AVP产品、CTHAGRI之间的空间距离

<0.5°;k=1,2,3,…。

(2)

假设k=q,D(q)=min(D)，则提取q位置上的垂直廓线数据或者云顶高度数据进行后续工作，其中下标m为AVP或CTHAGRI。使用的具体时空匹配原则如公式(1)和公式(2)所示，其中，L和L′分别表示数据的经度和纬度，D是CTHVFM与其余两组数据之间的空间距离，Δt是指CTHVFM与其余两组数据之间的时间间隔。

2 低云云顶高度反演

2.1 反演方法

本文的“逆温法”参考YI et al.[32]的方法，获取大气廓线上逆温层底的位置作为低云云顶高度。尝试使用FY-4A卫星反演的大气温度廓线，直接从反演的温度廓线上确定逆温层底的位置，从而估计云顶高度。需要注意的是逆温层存在季节性变化的特点[48-50]，一般在4月逆温层强度强，7月逆温层强度弱[6]，这可能导致“逆温法”的反演效果存在季节性差异，详细内容在2.3节讨论。

利用“逆温法”反演海上单层低云云顶高度的具体方法是，当得到满足1.2节匹配方法的GIIRS大气温度廓线后，依据以下公式提取最低逆温层底对应的高度当作低层层状云的云顶高度，即CTHGIIRS：

T(k)

(3)

Q(k)≠3，4，

(4)

其中，T为温度，i为大气温度廓线的层数，即逆温层底的温度小于相邻层的温度，Q为温度廓线的质量标签。需要注意的是，FY-4A卫星反演的AVP需要进行质量控制，其中AVP质量控制标签0为“prefect”，质量控制标签1为“good”，质量控制标签2为“bad”，质量控制标签3和4为“do not use”。反演出CTHGIIRS前，需要剔除廓线质量标签为3和4的点[42-43]。

选取单层低云个例的时间范围是2019年1月18日—12月31日，在挑选单层低云完成时空临近点匹配以及对温度廓线的质量控制后，在750个GIIRS/AVP上利用“逆温法”反演单层低云的云顶高度。

2.2 统计分析

利用2.1节中所述的云高反演方法，得到CTHGIIRS。统计结果显示，在“逆温法”得出的云顶高度和VFM探测的云顶高度之间，平均高度差为-0.27 km，均方根误差为1.03 km，相关系数通过5%的显著性水平检验。如表1所示，CTHGIIRS与CTHVFM之间的云顶高度差(cloud top height difference，CTHD)(CTHDGIIRS-VFM)的平均值、均方根误差明显优于CTHAGRI与CTHVFM的云顶高度差(CTHDAGRI-VFM)，表明对低云的反演上，CTHGIIRS比CTHAGRI表现出色。

图2显示，在|VCTHD|<500 m的区间内，两组CTHD概率分布分别是PGIIRS-VFM为37.7%，PAGRI-VFM为29.1%，CTHAGRI、CTHGIIRS与CTHVFM仍有显著差距，而且在VCTHD=0线周围，CTHDGIIRS-VFM概率分布明显更多，证明在统计结果上，相比CTHAGRI，CTHGIIRS更吻合CALIPSO的观测结果。

表1 低云云顶高度差(CTHD)统计结果

图2 低云云顶高度差概率分布Fig.2 Probability density function of low cloud CTHD

两者整体概率分布(图2)和累计分布(图3)显示，CTHDAGRI-VFM和CTHDGIIRS-VFM总体分布均倾向于负值，两者的CTH反演结果均偏小，最大峰值均位于[-500，-250]m之内，其中，PAGRI-VFM分布在VCTHD<0的区域内高达71.6%。CTHDGIIRS-VFM分布为偏左单峰结构，与CTHDAGRI-VFM的分布结构较为一致，在VCTHD∈(-1 250，0]m时，两者的概率非常接近(PGIIRS-VFM=47.6%，PAGRI-VFM=47.2%)，同时两者分布函数斜率也基本一致。在最大峰值两侧的翼区，CTHDGIIRS-VFM与CTHDAGRI-VFM的差异十分显著。在VCTHD∈(0,750]m时，概率分布上，CTHDGIIRS-VFM概率24.4%远高于CTHDAGRI-VFM概率13.8%，分布函数上，CTHDGIIRS-VFM的涨幅明显大于CTHDAGRI-VFM的涨幅；而VCTHD∈(-2 250，-1 250]m，情况则恰好相反。这样的概率分布和累计分布正是CTHDGIIRS-VFM的平均值、均方根误差明显优于CTHDAGRI-VFM的原因。

图3 低云云顶高度差累积分布函数Fig.3 Cumulative distribution function of low cloud CTHD

进一步分析VFM与AGRI、GIIRS之间的对应关系。如图4所示，发现CTHAGRI的频率fAGRI主要在对角线上半区集中分布，CTHAGRI主要集中在250～1 000 m之间，与CTHVFM相比，CTHAGRI显著偏小；其频率分布显示CTHAGRI没有和CTHVFM一同变化的趋势，证明CTHAGRI与CTHVFM之间无明显的对应关系；CTHGIIRS的频率fGIIRS主要在对角线上半区有集中分布，表明整体CTHGIIRS同样偏小；同时，fGIIRS沿对角线离散分布，可见CTHGIIRS与CTHVFM之间存在一定的对应关系。频率图的分析结果与表1中统计结果一致。

2.3 结果讨论

相比FY-4A云顶高度产品CTHAGRI在低云上的效果，“逆温法”反演的低云云顶高度CTHGIIRS有一定提升，但距离CTHVFM仍存在差异，其原因可能与GIIRS/AVP时空分辨率等因素相关。

(1)GIIRS/AVP产品分辨率。AVP与CTHVFM采取最邻点匹配的方法，寻找与CTHVFM时间间隔最小、空间距离最近的AVP进行反演，但由于AVP本身的时空分辨率原因，存在一些不可避免的客观误差。此外，相较无线电探空数据10 m/层，AVP垂直分辨率相对粗糙，得到的逆温层结构可能不精准，这也对反演结果造成影响。

图4 CALIPSO/VFM的低云云顶高度(CTHVFM)与两种云顶高度产品(a. CTHGIIRS， b. CTHAGRI)之间对应关系的二维联合直方图(网格分辨率：0.25 km×0.25 km；色阶：占比，白色表示不存在)Fig.4 Two-dimensional joint histogram of the correspondence between CTHVFM and two CTH products (a. CTHGIIRS, b. CTHAGRI) (grid resolution: 0.25 km×0.25 km； color scale: percentage, white grid indicates no values)

(2)海上逆温层存在季节差异。图5显示了不同季节下两种方法反演的云顶高度与CTHVFM之间云顶高度差的平均差(mean deviation，MD)和标准差(standard deviation，STD)。春、冬季，MDGIIRS-VFM分别为-0.13 km和-0.29 km，低于其在夏、秋季节的数值，说明“逆温法”在春季和冬季的反演效果要好于夏季和秋季。同时，夏季的逆温层相对较弱，大气温度廓线难以捕捉逆温层位置信息，导致夏季低云采样比例非常小。但是，在其他季节，MDGIIRS-VFM仍明显优于MDAGRI-VFM，如图5所示。

图5 不同季节下FY-4A和VFM之间差异的平均值(实线)和标准差(虚线)(红色表示GIIRS-VFM，蓝色表示AGRI-VFM，绿色柱为采样占比)Fig.5 Mean (solid line) and standard deviation (dotted line) of differences between FY-4A and VFM in different seasons (red line indicates GIIRS-VFM, blue line indicates AGRI-VFM, and green bar represents the sampling proportion)

(3)不同云类的云顶逆温结构存在差异。如表2所示，逆温层底高度与低层层状云的云顶有十分明显的对应关系。根据ISCCP的分类，在气压大于680 hPa的低云中，不仅有层云、层积云，还有积云存在。积云的云顶与逆温层之间并没有发现明显的对应关系，而这类低云的云顶高度同样难以测准[15]。利用VFM的云类型分类(其中“low overcast，opaque”是指层云或者雾[51])，表2显示了不同云类型中两种方法反演的云顶高度与CTHVFM之间云顶高度差的平均差(MD)和标准差(STD)。正在转变的层积云和层云或者雾的CTHGIIRS平均高度差分别为-0.61 km和-0.40 km，比CTHAGRI有更出色的表现，证明该方法在低层层状云的云顶高度上有良好的效果；而因为积云不符合“云顶高度在逆温层底”这一假设，所以在积云个例上出现CTHAGRI优于CTHGIIRS的现象。

3 低云云顶高度个例分析

因为考虑有质量控制，CTHGIIRS的数目锐减，选取2个CTHGIIRS可比较数相对多的个例用于分析。

低云个例结果(图6)显示，本文方法暴露一些明显缺陷。图6a显示，CTHGIIRS零散地分布于CTHVFM和CTHAGRI之间周围，还有部分CTHGIIRS与CTHAGRI重合；除了极少数CTHGIIRS略高于CTHVFM，CTHGIIRS大多均小于CTHVFM。图6b中CTHGIIRS的表现略优于图6a，CTHGIIRS数目较多，并且在22°N以北的低云区，CTHGIIRS更接近CTHVFM，但总体来看，大部分CTHGIIRS仍低估低云的云顶高度。显而易见地，虽然相较于CTHAGRI，CTHGIIRS与CTHVFM之间的差距更小，但是CTHGIIRS的分布相对稀疏，无法完整反映低云云顶高度的连续变化。由于对温度廓线进行了质量控制，2个西北太平洋低云个例上可提供的CTHGIIRS数目也不一致。而CTHAGRI也存在明显问题，CTHAGRI数值与CTHVFM相比总体偏小，同时其刻画的云顶变化趋势，也与CTHVFM相距甚远。以上结果显示：就云顶高度的数值而言，CTHGIIRS和CTHAGRI的效果与之前的统计分析结果显示一致；就单点的云顶高度而言，两者对CTH估计均总体偏低，CTHGIIRS优于CTHAGRI。

表2 不同云类的低云云顶高度差(CTHD)统计

图6 两种云顶高度产品(CTHAGRI和CTHGIIRS)与CALIPSO/VFM的低云云顶高度(CTHVFM)对比(a. 2019年1月28日，b. 2019年2月1日)Fig.6 Comparison between two CTH products (CTHAGRI and CTHGIIRS) and CTHVFM in two low cloud cases on 28 January (a) and 1 February (b) 2019

4 结论与讨论

(1)基于逆温层与低层层状云之间的关系，尝试利用FY-4A/GIIRS/AVP，提出一种在西北太平洋区域内反演低云云顶高度的方法。对比2019年CALIPSO探测的西北太平洋单层低云云顶高度和FY-4A/AGRI/CTH、“逆温法”反演的CTH之间的差距，结果显示，CTHAGRI减去CTHVFM的平均云顶高度差CTHDAGRI-VFM为-0.49 km，CTHGIIRS减去CTHVFM的平均云顶高度差CTHDGIIRS-VFM为-0.27 km，CTHGIIRS和CTHAGRI普遍低估低云的云顶高度。同时，CTHAGRI与CTHVFM的相关性不显著，CTHGIIRS与CTHVFM存在一定相关关系，可见在统计结果上，CTHGIIRS比CTHAGRI更加接近CTHVFM。但即便“逆温法”反演的云顶高度CTHGIIRS对比已有的CTHAGRI，对低云云顶高度的反演上取得更理想的统计结果，在具体的西北太平洋低云个例上，CTHGIIRS和CTHAGRI均不能准确描述低云云顶高度的变化，因此目前卫星对低云的云顶高度反演还有很大的提升空间。

(2)对CTHGIIRS与CTHVFM之间产生误差的原因，主要从GIIRS/AVP时空分辨率、低云出现的季节和低云种类3个方面进行研究。GIIRS/AVP的时空分辨率对反演的CTHGIIRS产生一定影响。使用就近原则比较云顶高度，CTHGIIRS与CTHVFM在客观上存在一定差距；而AVP垂直分辨率为230 m/层，使用相对粗糙的垂直分辨率可能对逆温层底高度判断不准确，导致CTHGIIRS与CTHVFM的差异；同时，逆温层有季节变化的特征，在不同季节，CTHGIIRS也会受到影响；并且发现CTHGIIRS对低层层状云的反演更为准确。

(3)利用FY-4A/GIIRS/AVP和“逆温法”反演的低云云顶高度优于现有的FY-4A/AGRI/CTH，为西北太平洋低云云顶高度的反演提供了新思路。该方法反演的云顶高度在低云个例的表现上仍不理想，其准确性受到AVP产品的时空分辨率、季节和低云云种的影响。同时，云顶高度产品的时间连续性也十分重要，然而本文用于检验的CALIPSO的云顶高度数据对同一区域一天仅能扫描2次，难以满足时间连续的条件。下一步将从西北太平洋的各种气象要素等方面入手探究提高准确性的方法，因为该方法的反演结果也有可能受到海雾或者低云发生时的天气形势影响。期待在未来可以应用到时空分辨率更高、质量更优的FY-4A/GIIRS/AVP产品和FY-4A成像仪数据。

致谢：感谢国家卫星气象中心(National Satellite Meteorological Center，NSMC)提供的GIIRS AVP反演数据和AGRI CTH数据，以及美国国家航空航天局(National Aeronautics and Space Administration，NASA)兰利研究中心的大气科学数据中心(Atmospheric Science Data Center，ASDC)提供的CALIPSO数据。