Overview of Retrieval Algorithm of Cloud-top Height Based on Satellite Infrared Remote Sensing

FAN Hongjie1,2, HUANG Yipeng1, LI Wanbiao1,†

Review FAN Hongjie, HUANG Yipeng, LI Wanbiao

1. Department of Atmospheric and Oceanic Sciences, School of Physics, Peking University, Beijing 100871; 2. 63892 Troop, Luoyang 471003; † Corresponding author, E-mail: lwb@pku.edu.cn

Abstract According to recent years cloud top-height inversion algorithms using satellite visible and infrared remote sensing data, domestic and foreign algorithm research progress and problems are summarized. Through the analysis of large number of related literatures, the inversion algorithms of the infrared window methods and CO2 absorption technique are mainly introduced. On the basis of these algorithms, the current business algorithms in countries and regions are expressed. Then the advantages and disadvantages of these methods in application are reviewed, and the cloud top-height inversion solution for regions in China is present. In the future, the method of cloud top-height lookup tables based on infrared split-window is prospective. In application and research, lookup tables based on brightness temperature difference and 11 μm channel brightness temperature are constructed using split-window data of stationary satellite and cloud profile radar data of polar orbit satellite, and cloud top-height can be acquired by checking lookup tables by split-window brightness temperature difference and 11 μm channel brightness temperature in application. Key words satellite infrared remote sensing; cloud-top height; retrieval algorithm; infrared split-window

云顶高度及云的辐射特性会影响地气辐射的收支, 对大气能量收支平衡具有显著的调节作用, 同时云顶高度等参数在航空气象保障、数值天气预报等领域也有重要体现, 因此确定云顶高度等参数对于大气物理及气候研究、气象保障等方面具有重要

的现实意义。云顶参数主要包括云顶高度、云顶气压、云顶温度、有效发射率、光学厚度和云相态等, 这些参数一般难以直接测量, 必须基于相关的可测物理量进行反演后才能得到。目前基于卫星红外遥感的云顶高度反演就是一种广泛使用的手段。

卫星遥感技术可为研究人员提供不同波段的云图数据, 包括可见、短波和长波红外、微波、水汽云图等图像数据。在利用卫星红外遥感数据在云顶高度反演方面, 国内外都进行了相关的研究工作。利用卫星红外遥感数据反演云参数的研究工作开始于 20 世纪 70 年代。国外较早地开展了利用卫星遥感数据反演云顶高度的研究, 形成一些成熟的算法, 并应用于遥感卫星产品中。国内学者也针对这些方法进行了应用研究及改进工作, 大多集中在利用现有的卫星云产品进行一些统计工作方面[15]。本文主要介绍利用卫星红外遥感数据反演云顶高度的算法原理、过程以及进展, 主要包括以红外窗区和 CO2 吸收技术为主的红外波段反演算法。在这些算法的基础上, 进一步介绍目前相关国家和地区的业务算法, 并对这些算法在实际应用中的优缺点予以评述, 提出针对中国地区云顶高度反演的解决方案。

1 红外窗区为主的算法1.1 单红外窗区法

单红外窗区法[610]主要利用大气红外窗区(3~5 μm 和 8~14 μm窗区通道)中某一通道的辐射数据来确定云顶高度。该方法假定云的发射率为 1, 即将云层当做黑体, 且云的散射及云顶上方的辐射和吸收均可忽略, 将卫星测得的辐射近似为云顶的辐射,即 Iλ=bλ(tc), 其中 Iλ为辐射亮度测量值, Bλ(tc)为温度 Tc下的黑体辐射亮度。首先根据卫星测量的辐射亮度, 利用普朗克公式解算得到云顶亮温, 然后由云顶亮温与大气温度廓线(当地测量的廓线数据或通过 NCEP/NCAR 及其他数值模式预报的廓线数据)对比得到云顶高度, 因此该方法又称为亮温比对法。常用的窗区通道为 11 μm 红外窗, 包括两个通道: 长波红外通道一(10.3~11.3 μm, 简称 11 μm 通道)和长波红外通道二(11.5~12.5 μm, 简称12 μm 通道)。在 GOES, FY-2C/E 和 Himawari 等系列卫星上均有相同或类似通道的探测器。单红外窗区法是反演云顶高度最为直接和简单的方法。由于将云层当做黑体来处理, 忽略了云的透过率及云顶大气消弱的影响, 因此该方法仅适合于较高的不透明云。对于较低的云层、半透明云及亚像元云层(云层未充满传感器像元视场)会出现较大的估计偏差, 半透明云和亚像元云的云顶温度往往会被高估, 云顶高度会被低估。

对于该方法, Menzel 等[6]、Smith 等[8]、Schreiner等[9]以及 Nieman 等[10]分别做了详细的分析, 并与其他方法进行比较, 均得到单红外窗区法对云顶平均高度的低估偏差通常在 100 hpa 以上, 其中

[9] Schreiner 等 还特别针对卷云的云顶高度进行反演, 并做统计分析, 得到卷云云顶高度的低估偏差达到 350 hpa 以上。因此, 他们指出该方法通常用于较高的不透明云(顶部气压低于 600 hpa 的不透明云)。

太阳光反射率‒红外窗区法1.2 红外窗混合算法1.2.1

‒太阳光反射率 红外窗区法[1112]是利用 0.2~4 μm 波段通道和 8~12 μm红外窗口通道反演卷云的云顶高度。单红外窗区法对于半透明云和亚像元云的云顶高度估计偏差较大, 根本原因是云层的有效‒红外窗区法,发射率无法准确估计, 为了解决这一问题, Shenk等[11]提出太阳光反射率 该方法是对单红外窗区法的修正, 由云层对太阳光的反射辐射估计云的覆盖率, 由覆盖率和红外通道测量值得到云的红外辐射和大气廓线数据来估计云顶高度。卫星探测器接收到的辐射亮度包含云自身的辐射亮度、透过云层的云下背景辐射亮度、云上大气路径辐射亮度和像元视场中非云区的晴空辐射亮度(云未充满探测器视场时需要考虑)。探测器接收到的辐射亮度 I 可用下式表示:

I = Nεbcldt(pc, 0) + N(1 − ε)i′clrt(pc, 0) + Nipath(pc, 0) +(1 − N)iclr =ε′icld +(1 − ε′)iclr = Iclr + ε′(icld − Iclr) , (1)

Icld = Bcldt(pc, 0) + Ipath(pc, 0)表示卫星探测器接收到的黑体云云顶辐射亮度加云上大气辐射亮度; Iclr = I′clrt(pc, 0) + Ipath(pc, 0)表示晴空辐射亮度; Bcld, I′clr和 Ipath(pc, 0)分别表示同云温黑体辐射亮度、云下背景辐射亮度和云上大气辐射亮度; t (pc, 0)表示云顶到大气层顶的透过率; N 表示像元中云覆盖率; ε表示云的发射率, ε′表示有效发射率, ε′ = Nε。式(1)中 I和 Iclr可用下列公式表示: cld

式中, Tg和 pg分别表示地表温度和气压; Tc 和 pc 分别表示云顶温度和气压。

Icld 和 Iclr可以通过国际上流行的大气辐射传输模式来计算, 主要有 LOWTRAN, MODTRAN 和 6S等模式, 目前常用的是MODTRAN模式。MODTRAN从 LOWTRAN 发展而来, 是一种中光谱分辨率大气辐射传输模式, 是美国空军地球物理实验(AFGL)开发的计算大气透过率及辐射的软件包。MODTRAN提高了 LOWTRAN 的光谱分辨率, 改进了多次散射的算法等, 通过对几何路径、大气模式、气溶胶模式、传感器特性和地表特性进行设置, 可用于计算大气路径透过率、路径的散射辐射和路径大气热辐射等。MODTRAN 反演云顶高度的过程如下。

1) 利用 0.2~4 μm 波段卫星图像亮度、卫星观测角和太阳天顶角计算云的覆盖率 N, 计算公式[12]如下:

式中, Isun,meas表示探测器接收到的反射太阳辐射亮度; Isun 表示大气层外太阳辐射亮度(假设在大气不同层的取值相同); ρclr 表示晴空太阳辐射双向反射率, 其值由下垫面类型决定; ρcld表示云层太阳辐射双向反射率。2) 假设云(水)的红外发射率为 0.9, 晴空辐射亮度 Iclr由探测器接收的云图数据得到, 通过式(1)、Iclr 和 N 得到 Icld, 计算公式[12]如下:

3) 根据大气廓线数据, 可以利用式(3)或者MODTRAN 大气辐射模式, 计算得到一系列辐射亮度值 I′cld(tc,i)(理论计算值, 其中 Tc,i 表示大气廓线数据分层温度)。

4) 用 χ(tc,i) = |I′cld(tc,i) − Icld|表示理论计算值与测量计算值之间的误差, 由大气廓线数据计算得到一系列 χ(tc,i)值。当 χ(tc,i)取最小值时, 对应的大气廓线分层温度 TC,I 为云顶温度, 对应的高度为云顶‒高度。太阳光反射率红外窗区法依赖于云在可见光波段的反射率和红外波段的发射率, 因此最大的缺点是只适用于白天, 不适用于夜晚, 但反演精度比

单红外窗区法高。Shenk 等[11]利用该方法得到卷云的云顶高度反演偏差在±50 hpa 范围内; 云层发射率越高, 精度越高。Reynolds 等[12]利用该方法反演云顶高度得出非卷云云层的云顶高度反演偏差为H2O‒红外窗区法0.5 km, 卷云的云顶高度反演偏差为 1.1 km。H2O‒ 1.2.2 H2O‒ 红外窗区法[6,810,13]又称线性外推法或红外窗截距法, 指利用 5.7~7.1 μm H2O 吸收通道和 10.5~12.5 μm 红外窗区通道辐射计算卷云云顶高度。将式(1)分别应用于 H2O 吸收通道和红外窗区通道:

I6.7(T6.7) = I6.7(TF) + ε′6.7[i6.7(tcld) − I6.7(TF)] , (7) I11(T11) = I11(TG) + ε′11[i11(tcld) – I11(TG)] , (8)

式中, 下角标 6.7 和 11分别表示水汽通道和红外窗区通道, T6.7 和 T11分别表示两个通道的亮温, TF 和TG 分别表示两个通道的晴空亮温, Tcld 表示云顶温度, I6.7(T6.7)和 I11(T11)分别表示两个通道的辐射亮度测量值, I6.7(TF)和 I11(TG)分别表示两个通道的晴空辐射亮度, I6.7(tcld)和 I11(tcld)分别表示两个通道的黑体云辐射亮度, ε′6.7 和 ε′11分别表示两个通道的有效发射率。两个通道的有效发射率近似相等, 即   11 , 6.7消去式(7)和(8)中 和 , 得到6.7 11

I (T6.7) = AI11(T11) + b, (9) 6.7 a = [I6.7(tcld) − I6.7(TF)] / [I11(tcld) – I11(TG)] , (10) b = [I6.7(TF)I11(TCLD) − I6.7(TCLD)I11(TG)] / [I11(tcld) – I11(TG)] 。 (11)

用 L1 表示式(9)对应的曲线, 通过测量得到同一块卷云多个区域的两通道辐射亮度值, 则可确定a,b 值, 得到 L1的表达式。另外, 当云层为黑体云(即云层的两通道有效发射率为 1)时, 从式(7)和(8)得到 T6.7 = T11 = Tcld,相应地, 两通道的辐射亮度可表示为 I6.7(tcld) 和I11(tcld), 两者之间的关系可用下式表示:

I6.7 (Tcld) = f I11(tcld) , (12)

式中, f表示两通道辐射亮度比。用 L2表示式(12)对应的曲线, 通过计算一系列 Tcld 值下的 I6.7(tcld)和I11(tcld), 即可得到 L2 曲线(Tcld 取值范围可为 200~ 300 K, 或按照大气廓线温度数据取值)。

L1 和 L2 相交处的亮温值即为实际的云顶温度,由大气廓线数据即可得到云顶高度。实际应用时,可按以下步骤计算云顶高度。

1) 根据温度廓线, 用普朗克公式计算出一系列两个通道不同温度下的辐射亮度值, 得到曲线 L2。

2) 根据两个通道实测图像数据中对应两个点或多个点的辐射亮度值得到曲线 L1。

3) 计算 L1 和 L2 的交点对应的亮温, 即为云顶温度。4) 比较大气温度廓线数据, 得到云顶高度。1981 年, Szejwach[13]首次利用 H2O 吸收通道和红外窗区通道对卷云云顶亮温进行反演, 然后与水

H2O‒红外窗区法反演结果更接近云汽通道和红外通道两个单独通道获得的亮温做比较, 结果显示顶高度处的温度, 并指出该方法在反演云顶高度的H2O‒红外窗区法做了详细的分析,应用前景。其后, Menzel 等[6]和 Nieman 等[10]分别

H2O‒对 并与其他方法比较, 均得到 红外窗区法在云顶平均高

H2O‒红外窗区法得到云顶高度度估计中的偏差通常在 60~110 hpa 之间的结论。

[9] Schreiner 等 利用

H2O‒红外窗区法优于单红外窗区法,的估计偏差在 80 hp a 以内。Smith 等[8]同样做了分析比较, 得出 ‒稍次于 CO2吸收技术。

因此可知, H2O红外窗区法反演精度在 60~ ‒红外窗区法和110 hpa 之间, 优于单红外窗区法, 稍次于太阳光反‒射率 CO2吸收技术, 但相对于太阳

H2O‒光反射率 红外窗区法, 具有不受时间限制的优点。由于 红外窗区法利用了水汽吸收通道,因此在反演较低云顶高度时会受到限制, 并且从反演步骤 2 可以看出, 反演时要求云层能够充满多个像元视场, 因此该方法不适用于尺度过小的云层。

1.3 红外分裂窗查算表算法

[1419]红外分裂窗查算表算法 是利用长波红外的分裂窗通道(11 μm 通道和 12 μm 通道)数据和主动式高精度设备测量得到的云参数数据, 联合建立云顶高度查算表。其中分裂窗通道提供云层的 11 μm 通道亮温及两通道的亮温差值(用 BTD 表示,通常用 11 μm 通道亮温减去 12 μm 通道亮温), 主动式高精度测量设备提供云层的垂直结构信息, 以获取云顶高度参考值。

利用分裂窗所基于的物理基础是: 1) 在两个通道, 低层云或厚云(如积云、积雨云和层积云)的发射率是一致的, 可以看做发射率为 1 的黑云; 2) 在 两个通道, 卷云的发射率比较小, 且存在差别; 3)在红外窗主要是水汽的少量吸收, 在 12 μm通道的吸收比 11 μm通道的吸收更大一些。两个通道的发射率有如下关系[15]:  1  (1  )1.08 , (13) 12 11其中, ε11′和 ε12′分别表示两个通道的有效发射率。利用式(1), 用 11 μm通道值减 12 μm 通道值得到:

其中, Iclr,11 和 Iclr,12, Icld,11 和 Icld,12分别表示两个通道的晴空辐射亮度和黑体云辐射亮度; ∆Iclr=iclr,11 – Iclr,12,为两通道晴空辐射亮度差; ∆Icld=icld,11 – Icld,12, 为两通道黑体云辐射亮度差;  –  , 11 12为两通道有效发射率差。可做如下分析。1) 当云层较厚(如积云、层积云和积雨云)时,可以看做黑云, 发射率相同(近似为 1), 即   11  1,  0, I  I , 在红外窗区两个通12 11-12 cld道辐射亮度近似相等, 则有 ∆I11-12 = 0。2) 对于晴空来说, ε11′ = ε12′ = 0, ∆I11-12 = ∆Iclr,在红外窗区主要是水汽的少量吸收, 在 12 μm 通道的吸收比 11 μm 通道的吸收更大一些, 则在晴空的两个通道 ∆Iclr与水汽的含量(与路径长度正相关)有关, 路径越长, ∆Iclr 越大, 因此晴空 BTD 值一般为正而不为 0。可以看出, 如积状云的厚度较大, 含水量高,透过率很小, 发射率大, 则 BTD 值小; 当为卷状云时, 透过率较大, 有效发射率较小, 则 BTD 值较大; 对于晴空来说, 云有效发射率为 0, ∆Iclr 不为 0,则 BTD 为正值。在陆地上时, 湿度的减少会使得∆Iclr 变小, 则 BTD 值变小, 相当干燥的地方接近0。这样, 不同的 BTD 值可以反映不同的云层类型,不同类型的云层具有不同的云顶高度。同样, 不同的亮温也反映不同的云顶高度, 那么云顶高度则随BTD 和 11 μm 通道亮温具有一定的分布规律。利用主动式高精度测量方法得到的云参数数据, 联合

两通道亮温差和 11 μm 通道亮温, 即可建立云顶高度查算表, 在获取云顶高度时就可以利用该表进行反演, 从分裂窗数据获取两通道亮温差值和 11 μm通道亮温值, 比对查算表, 即可获取云顶高度。红外分裂窗查算表算法反演云顶高度步骤如下。

1) 选取若干年遥感卫星分裂窗图像数据和载有主动遥感探测器的卫星数据(如 94 GHZ 的星载毫米波或星载激光云廓线雷达数据), 红外图像数据与主动遥感数据序列在时间和空间上一一对应。

2) 从主动遥感数据(云产品数据)得到云顶高度数据, 从分裂窗图像数据计算得到 BTD 值和 11 μm通道亮温值。

3) 将云高数据、BTD 值和 11 μm 通道亮温值一一对应, 进行填充处理, 建立查算表。

4) 由预处理的分裂窗图像, 计算得到 BTD 值和 11 μm 通道亮温值, 与建立的查算表比对, 得到相应的云顶高度。

红外分裂窗最先被 Inoue[1415]用于热带海洋上空云的检测分类(利用 NOAA7-AVHRR 数据)。刘

[19]诚等 也利用红外分裂窗对中纬度陆地上空云进行分类检测。反演云顶高度的红外分裂窗查算表算

[1718]法由 Hamada 等 首先提出, 他们 2008 年利用GMS-5 红外分裂窗数据及地基雷达测量的云参数建立查算表[17], 但是仅适用于非降水云; 然后在2010 年又利用 MTSAT 卫星红外分裂窗数据和Cloudsat 卫星云参数数据建立热带海洋上空的强对流云的云顶高度反演查算表[18], 利用查算表估计云顶高度偏差小于 1 km (云顶亮温在 240 K 左右)。该方法可用于全天候不同时段云顶高度的反演, 但是其精度受到纬度、昼夜、季节及卫星天顶角等因素的影响, 要提高反演云顶高度的精度, 则需要在建立查算表时最大限度地降低这些因素的影响, 可针对不同纬度、昼夜、不同季节分别建立相应的查算表。

2 CO2吸收技术

CO2 吸收技术[6,810,2040](CO2 Absorption Technique, CO2AT)是到目前为止研究及使用较广泛的反演云顶高度的方法。利用 CO2吸收带反演云顶高度有两种方法: 辐射比值法(或称 CO2切片法)和辐射余差法(包括均方根值法、Chahine 辐射余差法和MLEV 法)。使用的波带包含 14.2, 13.9 和 13.61 μm 的 CO2吸收带及 11 μm红外窗区。

2.1 辐射比值法

从探测器测量得到的云层辐射亮度中减去晴空辐射亮度, 可以利用以下公式:

I – Iclr =ε′(icld –Iclr)。(15)如果使用同一视场中两个相邻的探测器通道,云量相同, 发射率近似相等, 可得有效发射率相等,则有下式:

其中, λ1 和 λ2分别表示两个通道的响应波长, I(λ1)和I(λ )分别表示两个通道测量的云层辐射亮度, Iclr(λ1)

2和 Iclr(λ2)分别表示两个通道测量的晴空辐射亮度, ps和 pc 分别表示地面和云顶气压, t ( p,0) 表示气压 1为 p处至大气层顶的垂直方向大气透过率; f(λ1, λ2, pc)和 f(λ1, λ2, p)分别表示测量数据和理论计算数据下的比值函数。这样, 在已知两个通道的测量值及大气温度和透过率廓线的基础上, 就可以求出云顶高度。计算步骤如下。

1) 由两个通道的图像数据, 得到两个通道的观测亮度值和晴空亮度值。

2) 由两个通道的观测亮度和晴空亮度, 计算比值函数 f(λ1, λ2, pc)。

3) 由温度廓线和透过率廓线, 计算比值函数f(λ1, λ2, p)。4) 计算| f(λ1, λ2, pc) – f(λ1, λ2, p)|最小值对应的p, 即为云顶气压 pc。5) 比对大气廓线数据, 得到云顶高度。辐射比值法由 Smith[20]于 1968 年提出, 然后将其用于云顶高度等参数的反演工作[2122], 高度反演结果与激光雷达、无线电探测仪测量结果匹配较

[2325]好。Menzel 等 利用该方法反演云顶高度的误差为 50 hpa。menzel 等[6]、Smith 等[8]、Schreiner

H2O‒11等[9]以及 Nieman 等[10]都对 11 μm 窗区法、CO2 比值法和 μm 窗区法反演云顶高度做了比较,得到辐射比值法的反演误差在 60~110 hpa 之间,均得出 CO2切分法反演精度在 3 种方法中最高的结论。Baum 等[26]发现在两层云(底层为不透明云)的情况下, 利用辐射比值法反演的结果比较差。Wylie 等[27]将 HIRS 反演的云顶高度(利用辐射比值法)与 SAGE-Ⅱ卫星测量结果进行比较, 发现二者

[28]之间存在 1.6 km 的偏差。Frey 等 将 MAS 卫星反演的云顶高度(利用辐射比值法)与机载激光雷达测试结果进行比较, 得到 1/3 情况下偏差为±500 m, 2/3 情况下偏差为 1500 m。naud 等[30]对 MODIS 反演的云顶高度(利用辐射比值法)与后向散射激光雷达测量结果进行比较, 得到高云误差范围为–1.4~ 2.7 km, 低云误差范围为–1.2~1.5 km。为了进一步提高 CO2比值法在多层云顶高度的反演精度, Holz

[31]等 在 CO2比值法基础上提出 CO2 筛选法, 并与辐射比值法和 MLEV 算法结果做比较, 结果显示CO2筛选法结果稍好一些。Chang 等[3233]在 CO2比值法基础上提出基于有效背景辐射的 MCO2AT 方法, 反演结果与 CALIOP (星载云和气溶胶正交极化激光雷达)测量结果相比, 辐射比值法的偏差为2.5 km, MCO2AT 方法的偏差为 1 km。

利用辐射比值法反演云顶高度时, 如果通道的辐射与晴空辐射之间差值很小(小于仪器噪声), 该方法则不适用。这种情况常发生在低层的碎云和非常薄的卷云情况下。另外, 该方法对多层云的反演误差也较大。

2.2 辐射余差法

这一类方法包括均方根值法、Chahine 辐射余差法和 MLEV 法。均方根值法通过计算测量值与理论计算值之间的均方根来反演云顶高度, 计算公式如下:

其中, Irms(pc, ε′)表示 N 个通道的辐射亮度均方差; I i 表示第 i 个通道的辐射测量值; Ii(pc, ε′)表示第meas i 个通道的理论计算辐射值。Irms(pc, ε′)取最小值时对应的 pc为云顶气压, 与大气廓线对比即可得到云顶高度。计算步骤如下。1) 由红外图像数据得到所有通道的Ii 。meas 2) 由大气廓线数据, 利用式(2)计算得到所有 通道的I i 。clr

3) 在其中一个通道(该通道权重函数 dt/d(lnp)的峰值位置相对其他通道最低), 根据大气廓线数据, 选择一组可能的云顶气压 Tc 和 pc, 由式(3)计算得到 Icld。

4) 由同一个通道下的Ii , I i 和 Icld, 用下式clr meas计算得到有效发射率 ε′:

ε′与通道无关。5) 在其余通道同一组 Tc 和 pc 下, 由 ε′并利用

i式(1)~(3)计算得到 I (pc, ε′)。

6) 根据Ii 和 Ii(pc, ε′), 通过式(4)计算 Tc 和 pc meas下的 Irms(pc, ε′)。

7) 重复步骤 3 至步骤 6, 计算不同 Tc 和 pc 下的 Irms(pc, ε′), 取其中最小值对应的 Tc 和 pc, 比对大气廓线数据, 得到云顶高度。

Chahine 辐射余差法与均方根值法在本质上是一致的, 计算公式与式(17)稍有区别, 如下所示:

云顶高度与 χ(pc, ε′)取最小值时的云顶气压 pc 相对应。

MLEV 法由美国威斯康星大学开发, 与均方根法原理相同, 利用的是有效发射率的方差, 计算公式如下:

其中, ε′(pc )表示第 i 个通道的有效发射率,  ( p )表i c示 N 个通道有效发射率的平均值。LEV 取最小值时对应的 pc为云顶气压, 相应的高度为云顶高度。Chahine[37]1974 年首次将辐射余差法用于云顶气压及云高的反演, 得到单层云的云顶高度反演精度为 25 hpa。但是, Wielicki 等[38]对这一方法(使用均方根值法, 等同于 Chahine 辐射余差法)进行重复试验, 结果显示, 只有在对 600 hpa 高度以上云顶高度的反演中, 精度才能达到 25 hpa。官莉等[39]利用 MLEV 算法反演了云顶高度, 结果与激光雷达观测和 MODIS 业务云产品比较一致。利用辐射余差法反演云顶高度时, 存在与辐射

比值法相同的问题。同时, 在计算过程中要用到多个相邻通道的均值化处理, 通道越多, 误差越小,因此反演的精度会受到传感器在某一波长附近通道数的影响。

3一些国家和地区的业务算法3.1美国

NASA (美国国家航空航天局)的卫星系列先后有 1960—1965 年的 TIROS (泰罗斯)系列极轨卫星, 1964—1978 年的 Nimbus“雨云”系列极轨卫星, 2002年起的 A-train 系列卫星; ESSA (美国国家环境科学管理局) 1966—1969 年的系列极轨卫星; NOAA (美国国家海洋和大气管理局)的卫星系列有 1970— 2009 年的极轨卫星系列(POES, 即 NOAA 系列), 1975 年开始的静止业务环境卫星系列(GOES 系列),以及 2011 年开始的 NPOESS “联合极轨卫星系统”系列卫星。

目前关于云参数反演的焦点在 GOES 系列卫星、NOAA 系列卫星和 A-train 系列卫星。在红外波段及混合波段主要利用 VISSR 数据(GOES 1~3星)、VAS 数据(GOES 4~7 星)、Imager & Sounder数据(GOES 8~16 星)、AVHRR & HIRS 数据(NOAA系列卫星)及 MODIS 数据(Aqua 卫星)。目前在红外波段, 美国在反演云顶高度(或云顶气压)的业务H2O‒红外窗区法和上广泛应用的算法[810]有基于单红外窗区的算法、

CO2辐射比值法。3.2 欧洲欧洲 ESA (欧洲太空局)的主要气象卫星包括:静止轨道卫星主要是 1977—2012 年发射的 Meteosat系列卫星(Meteosat 1~3, Meteosat 4~7 (即 MOP 1~4)和 Meteosat 8~10 (即 MSG 1~3)); 极轨轨道卫星主要是分别于 2006 年和 2012 年发射的 Metop-a & Metop-b (均载有 HIRS/4 和 AVHRR/3 传感器)。

作为静止卫星, Meteosat 是云参数反演的重要数据源。根据其通道特点, 红外波段反演云顶参数已经业务化的算法有单红外窗法[67]和 H2O红外窗区法[10,13]。

3.3 日本

日本的主要气象卫星包括: 静止气象卫星为Himawari 系列卫星, 分别是 1977—2003 年相继发射的“地球静止气象卫星”, 即 GMS 系列卫星(GMS 1~5, 即 Himawari 1~5)、接替 GMS 的“多用途运输卫星” MTSAT 系列星(MTSAT 1/2, 即 Himawari

6/7)和 2014 年发射的 Himawari8。应用的主要对象为 MTSAT2 和 Himawari8, 主要波段为可见、水汽和分裂窗波段。

日本 Himawari 系列卫星主要由美国承包研制并发射, 业务上从美国继承了相应算法。近年来,

[1718]日本京都大学 Hamada 等 分别利用 GMS-5 和MTSAT-1R 卫星数据发展了基于分裂窗数据反演云顶高度的查算表法。

3.4 中国

中国的静止气象卫星为 FY-2 系列星(FY-2A/ B/C/D/E/F), 极轨气象卫星为 FY-1 系列星(FY-1A/ B/C/D)和 FY-3 系列星(FY-3A/B/C)。风云卫星云产品的生成, 在业务上主要利用

H2O‒ FY-2 系列卫星的红外窗和水汽通道数据, 考虑到利用红外窗反演半透明云带来的误差, 利用红外窗区法[41]得到云顶温度。

4目前存在的问题及解决方案设想4.1红外窗区方法存在的问题

1) 单红外窗区法。该方法适用于厚而稠密且完全充满视场的云层。对于部分云覆盖(亚像元)和薄云层或卷云而言, 其有效发射率不为1, 云层下方的红外辐射能透过云层, 此时卫星测得的云层辐射混叠了云层下方背景的辐射, 云顶的亮度及亮温

太阳光反射率‒红外窗区法。该方法依赖于会被高估, 云顶高度则被低估。

2)云在可见光波段的反射率和红外波段的发射率, 精度比红外窗区法高, 但是短波波段双向反射率的取值以及不同层取相同的太阳辐射值的假设都会对计算精度产生影响, 最大的缺点是只适用于白天, 不适用于夜晚, 同时当云的红外率越小时, 估计偏差‒越大。

3) H2O红外窗区法。该方法利用 H2O 吸收带通道, 水汽云图主要反映对流层上部的辐射信息。当大气中相对湿度增大时, 可观测的高度会相应上移, 在反演云顶高度时受到大气湿度廓线的影响,因此该方法对中高云云顶高度的反演比较有效, 对低云会失效, 且不适用于尺度过小的云层。

4) 红外分裂窗查算表算法。该方法目前主要应用在热带海洋地区, 由于海洋背景随时间和空间变化不大, 比较均一, 云类型阈值划分较简单一些,但是中纬度陆地上背景辐射随时间和空间比较多变, 分布不均, 阈值需要进行调整。该方法确定的

云顶高度偏差受纬度、太阳天顶角和季节等因素影响较大, 在高纬度云顶高度估值明显偏低, 而且受季节影响较大(夏季高于冬季), 在大天顶角地方云顶高度同样估值偏低。在应用中要考虑这些因素的影响。

4.2 CO2吸收技术存在的问题

1) 辐射比值法。当通道的辐射与晴空辐射差值很小(小于仪器噪声)时, 该方法则不适用, 常发生在低层的碎云、非常薄的卷云和高的不透明云等情况下。

2) 辐射余差法。该方法在计算过程中要用到多个相邻通道的均值化处理(要求云发射率在各通道近似相等), 通道越多, 误差越小, 因此反演的精度会受到传感器在某一波长附近通道数的影响。

4.3 解决方案设想

分裂窗查算表法主要利用长波红外大气窗口的两个波段, 是静止卫星通用的两个通道, 可以观测较高时间分辨率、较大地域范围和不同高度的云层, 加上主动式高精度设备测量得到的云参数数据的支持, 在热带海洋地区云顶高度反演中已经取得较好的结果。可以设想, 将分裂窗查表法应用于广泛地域空间及多个时间段上具有一定的可操作性及发展前途, 因此可以进行尝试性研究。考虑到该方法确定的云顶高度偏差受纬度、卫星天顶角和季节等因素的影响, 可以根据这些因素的影响程度建立多个区域、多个季节的云顶高度查算表。对于如何将基于分裂窗的二维查算表法很好地应用于中国大陆, 我们做如下分析及考虑。

1) 为了解决大纬度跨度及大范围卫星观测天顶角的影响问题, 同时也要考虑到中国大陆的海拔分布特点, 根据红外窗的晴空云图分析地表亮温分布特点, 兼顾纬度、卫星天顶角及海拔三方面, 对大陆进行地域上的分块划分。

2) 在地域分块划分的基础上, 根据分块地域地表亮温随季节的变化特点, 在时间尺度上进行二次划分, 不同分块地域季节划分尺度有所不同。

3) 在地域和季节的二维尺度上, 建立多个基于分裂窗亮温及亮温差的查算表。实际应用时, 根据地域和季节使用对应的查算表。

4) 有关建立查算表的方法, 可以参考 1.3 节的实施步骤。

5) 对于中国大陆, 在地域和季节二维尺度上进行划分时, 可以采用 FY-2D/E/F 静止卫星的红外分 裂窗数据和 Cloudsat 极轨卫星的云产品数据(2BGEOPROF)。应该注意的是,静止卫星和极轨卫星数据在空间和时间上一定要有很好的匹配。

5 展望

分析云顶高度反演方法的反演步骤, 可以看出,影响云顶高度反演精度最主要的两个因素是大气廓‒红线数据的准确性和云顶亮温的准确性。大气廓线数

H2O‒红外窗区法中的数据对比带来误据不准确, 会给单红外窗区法、太阳光反射率外窗区法和差, 也会在 CO2吸收技术中进行理论值计算时产生误差。云顶亮温的准确性影响表现在卫星反演的云顶亮温不等同于云顶实际温度, 这样也使得对比结果产生误差。提高云顶亮温准确性的途径之一是通过其他渠道获取更多的云层信息, 使得反演的云顶亮温更接近云顶实际温度; 途径之二是反演的云顶亮温不与廓线数据做直接的比对(如分裂窗查算表法), 但这需要可参考的高精度云顶高度值。随着不同用途的气象卫星数量及卫星探测通道的增加,联合其他波段通道数据(如可见、近红外、短波红外、中波红外及毫米波波段等)和极轨卫星数据来提高云顶高度的反演精度已经具有可行性, 因此,采用多个卫星数据及多波段通道数据联合建立云顶高度反演算法成为未来的发展方向。