APP下载

云偏振遥感综述

2022-06-09王宇瑶麻金继李婧晗洪津李正强

遥感学报 2022年5期
关键词:偏振气溶胶反演

王宇瑶,麻金继,李婧晗,洪津,李正强

1.安徽师范大学地理与旅游学院,芜湖241003;

2.中国科学院安徽光学精密机械研究所通用光学定标与表征技术重点实验室,合肥230031;

3.中国科学院空天信息创新研究院国家环境保护卫星遥感重点实验室,北京100101;

4.资源环境与地理信息工程安徽省工程技术研究中心,芜湖241003

1 引言

云由大气中的水蒸气经遇冷液化或凝华而形成的水滴和冰晶凝结组成,覆盖了全球67.5%的地球表面,随着时间和空间位置的变化以及下垫面的不同情况呈现出不同的分布及特征(Ilčev,2019)。IPCC 第六次评估报告指出,近年来在云的研究方面虽然取得了一些进展,但云仍然是气候变化最主要的不确定因子之一,其复杂性与多样性影响着全球气候模拟结果的准确性。(IPCC,2007,2021;Rossow 和Schiffer,1999)。同时,云作为太阳入射辐射和地球表面出射辐射的必经屏障,在大气能量交换、地球辐射收支、降水循环系统与地—气系统中也起着不可忽视的调节作用(Carney 和Vincent,1986;任建奇等,2010;Stephens 等,1990;Trenberth 等,2009);并且对地表与卫星传感器之间的信号传播具有严重的阻碍作用,会影响遥感信息的获取质量(Matsui 等,2004)。云遥感平台有地基平台、机载平台和星载平台,其中卫星观测平台因为可以进行全球尺度的云遥感研究而显得尤为重要(严卫等,2008)。由于辐射强度信号对于厚云,尤其是对流云的内部物理特征无法进行精确地探测,所以仅利用辐射强度信息无法实现高质量的云观测(杨之文等,2005);而被称为大气粒子特性的“敏感指示器”的偏振光谱信息能够探测出云粒子的尺度、形状等微物理参量,在云遥感领域展现出独特的应用优势(郝增周等,2006)。相对于传统光谱辐射云遥感而言,利用多角度偏振观测能更好地从观测数据中剥离出下垫面信息,从而获取更多更准确的大气散射信息(吴浩,2016)。自从法国发射了POLDER系列传感器以来,通过其在气溶胶和云产品方面的应用(Buriez等,1997;Goloub等,2000;Tanré 等,2011;汪杰君等,2021),国际上很多科学家意识到偏振传感器在探测气溶胶和云参量方面具有较好的优势(Nadal 和Breon,1999;Peralta 等,2007;Li 等,2018)。虽然目前国际上偏振传感器的发展很快,气溶胶偏振算法较多,但云偏振遥感的算法相对较少,本文梳理了近30年来偏振传感器的发展和云遥感算法的演变,并详细介绍了典型偏振传感器及其经典云算法的贡献;同时还对未来偏振传感器的发展方向及其对应的算法进行了展望。

2 云遥感偏振传感器

鉴于偏振传感器的蓬勃发展,基于已发射偏振仪器和待发射偏振仪器两个方面,从星载、机载和地基3 个平台进行载荷的详细介绍。表1 按首字母顺序汇总了所有载荷名称的简写、全称及中文名称,表2列举了已投入使用的偏振载荷。

表1 载荷名称汇总表Table 1 The summary table of the instrument name

表2 已投入使用的偏振载荷参数Table 2 Polarization load parameters that have been put into use

2.1 已发射偏振仪器

2.1.1 星载偏振仪器

(1)POLDER。法国空间局(CNES)为了验证偏振探测气溶胶和云的优势,发射了POLDER(POLarization and Directionality of the Earth’s Reflectances)系列传感器。POLDER-1 的发射有效验证了偏振在探测小粒子和水云方面的优势,但其空间分辨率有限;POLDER-2 更好地验证了偏振在高空间分辨率上探测气溶胶和云的优势,但由于卫星的故障而停止运行;POLDER-3 更进一步优化了POLDER-1、2 传感器的通道设置,于2004-12-18 搭载于PARASOL 卫星升空,并加入“A-Train”计划,在轨运行近10年结束观测任务,其波段信息如表3所示,结构示意图如图1所示(Andre和李玲,1997;Herman等,1997;徐飞飞,2018)。

图1 POLDER结构示意图(Andre和李玲,1997)Fig.1 The schematic diagram of POLDER(Andre and Li,1997)

表3 POLDER-3通道信息表Table 3 POLDER-3 channel information table

POLDER系列传感器通过获取对应通道顶部大气光谱辐射值从而获得辐亮度和偏振辐射值,与传统的光谱传感器相比增加了维度信息,从而更适合探测云和气溶胶的微物理特性。在气溶胶反演、地球辐射收支等研究方面作出了贡献(Bréon和Colzy,1999;Fougnie 和Bach,2009;Lacagnina 等,2017;Nadal 和Bréon,1999;Tanré 等,2011);同时可以进行云识别和更高精度的云参量反演,即云相态、云高、云滴谱、云光学厚度等(Buriez 等,1997;陈震霆等,2018;Goloub 等,1994;Herman 等,2005;李佳欣等,2020;Riedi 等,2001,2010;汪杰君等,2021;Yasumoto 等,2002)。

(2)APS。法国POLDER 系列传感器的发射及其气溶胶和云等系列产品的成功应用,让美国航天局(NASA)意识到了偏振探测在大气环境遥感领域的优势,同时也认识到低空间分辨率和低偏振精度对气溶胶反演的不确定性问题,从而NASA设计了非成像高偏振度的气溶胶偏振传感器APS(Aerosol Polarimetry Sensor)。APS 是一种摆扫式偏振仪,采用了空间、时间、光谱上同步的成像测量模式,具有很高的偏振精度(偏振度DoLP 优于0.2%),可与A-Train 系列的其他仪器协同测量(Schmunk,2018);同时APS 是沿轨道扫描的非成像仪器,这在一定程度上会影响数据使用范围,其波段设置如表4所示,APS 载荷如图2所示(Peralta等,2007)。

表4 APS通道信息表Table 4 APS channel information table

图2 荣耀卫星和ASP载荷示意图(根据(Peralta等,2007)修改)Fig.2 Illustration of the Glory satellite and ASP(Peralta et al.,2007)

APS 可以对全球云和气溶胶的分布、微物理特性、化学成分进行高精度和全范围的连续观测,以便分析气溶胶和云对全球气候变化的影响(Mishchenko 等,2007)。在云参量反演过程中,由于APS 具备高偏振精度的特点,获取的数据可用来反演云滴粒子半径、云滴数浓度、冰水路径、光学厚度等(Alexandrov 等,2015,2018;Segal-Rozenhaimer等,2018;Sinclair等,2019)。

虽然APS 于2011-03-04 发射失败,但由于它的航空样机RSP 在探测气溶胶和云参量等方面表现出良好的效果,使得APS 的设计理念具有重要的参考价值(Alexandrov 等,2015,2018;Segal-Rozenhaimer等,2018;Sinclair等,2019)。国际多个偏振计的概念设计都继承了APS的设计理念,如中国高分五号02 星搭载的POSP;乌克兰气溶胶计划即将发射的ScanPol(2.2.4 Aerosol-UA:ScanPol、MSIP)。

(3)DPC。鉴于POLDER 系列的空间分辨率较低,RSP反演气溶胶和云参量的良好性能,中国综合两个仪器各自设计理念的优点,研制了大气气溶胶多角度偏振探测仪(DPC),于2018-05-09搭载于高分五号卫星上成功发射,目前已有相关文献对高分五号载荷进行了总结(陈良富等,2021)。DPC 作为一台集多角度、多光谱和偏振探测能力于一体的星载传感器,对同一个像元,一般可获得9 个角度、3 个偏振波段(每个偏振测量波段采用3 个检偏方向获取偏振图像)、5 个非偏振波段的观测信息,共计126个有效观测值(陈斐楠等,2019;康晴等,2018;Li 等,2018;骆冬根等,2019),其波段设置见表5,DPC及GF-5示意图如图3所示(Li等,2018)。

表5 DPC通道信息表Table 5 DPC channel information table

图3 DPC和GF-5示意图(Li等,2018)Fig.3 Illustration of the DPC and GF-5 satellite(Li et al.,2018)

DPC 可用于大气气溶胶分布、地表植被分布及海水叶绿素浓度等多项研究,其在轨定标方法可为中国自主偏振仪器在轨定标提供理论基础(陈斐楠等,2019;顾行发等,2011;涂碧海等,2020;谢一凇等,2019;张苗苗等,2018);还可以通过获取云的强度辐射和偏振辐射信号,为云的识别和参量反演提供更多有效的观测信息,提高对云参量反演的精度,并利用多角度信息表征云的三维结构(常钰阳,2020;常钰阳等,2020;李超,2019;Li 等,2021;王佳佳等,2021;伟乐斯,2020;Yu等,2019)。

DPC 在大气环境探测方面具有广泛的应用前景,自2018年高分五号(GF-5)成功发射以来,目前GF-5 02星已搭载其改进版升空,同时大气环境监测卫星、高精度温室气体综合探测卫星以及陆地生态系统碳监测卫星等多个卫星均计划搭载该卫星载荷。

(4)DPC-II、POSP。自2018年高分五号成功发射以来,偏振成像仪DPC 在大气环境探测方面表现出广泛的应用前景,同时鉴于APS 高精度非成像偏振仪对提高探测精度有较大的帮助,中国高分五号的后续型号GF-5 02星将搭载两台具有偏振能力且可以进行协同观测的载荷,分别为可成像的大气气溶胶多角度偏振探测仪(DPC-II)和非成像的高精度偏振扫描仪(POSP),二者波段设置如表6所示。

表6 DPC-II与POSP波段参数设置Table 6 Parameter settings of DPC-II and POSP bands

DPC-II 在一代基础上增加了多角度观测角度数(从12 个观测角度增加至17 个观测角度)。其主要用来获取地球大气多角度多光谱偏振辐射数据,用于大气气溶胶和云的观测(包括云粒子尺寸分布、云相态、云光学厚度等),可获取全球大气气溶胶和云的时空分布信息。

高精度偏振扫描仪(POSP)采用沿轨扫描方式,9 个谱段全部为偏振通道,并将工作谱段拓展到紫外及短波红外。POSP 能通过“偏振交火”对DPC 实现精度传递(李正强等,2019)。这种设计组合为协同观测提供了非常宽的光谱和偏振范围(380—2250 nm)以及多角度能力,是现阶段气溶胶和云遥感综合观测能力最强的星载探测方案之一(杨洪春,2019)。目前两台载荷已于2021-09-07搭载升空。

(5)其他星载仪器:MAPI/MAI、CAPI、SGLI。多角度偏振成像仪MAPI 或称MAI(Multi-Angle Polarization Imager)搭载于天宫二号空间实验室(TG-2),于2016-09-15 成功发射,是国内首台应用于空间探测的多角度偏振仪器(郭俊杰等,2019a,2019b;张一鹏等,2020a,2020b)。MAI可以获取气溶胶和大气水汽含量等信息,有效探测卷云,反演云顶高度,并探究云粒子的尺寸、形状等微观特性,其外观如图4所示(景振华等,2020;殷德奎,2019;张伟,2016)。但由于受TG-2 轨道的约束,只能获取南北纬42°区间的数据(杨洪春,2019)。

图4 多角度偏振成像仪(殷德奎,2019)Fig.4 Multi-angle polarization imager(Yin,2019)

多谱段云与气溶胶探测仪CAPI(Cloud and Aerosol Polarimetric Imager)搭载于中国发射的TanSat 卫星上,具有5 个光谱通道,可以测量云、大气颗粒物等信息,从而能够为更准确地反演二氧化碳剔除干扰因素,并为研究大气污染物(如PM2.5等)的成因提供重要数据支撑(Chen 等,2017;Wang 等,2017),其可见光近红外单元及短波红外模块见图5(杨洪春,2019)。由于不具备多角度观测能力,CAPI 的探测能力受到限制(石光明等,2014)。

图5 CAPI可见近红外单元和短波红外模块(杨洪春,2019)Fig.5 CAPI’s visible and near-infrared unit and short-wave infrared module(Yang,2019)

第二代全球成像仪SGLI(Second generation Global Imager)搭载于日本全球变化观测任务的气候卫星上(GCOM-C),于2017-12-23发射(图6)。SGLI 用于长时间监测气溶胶—云相互作用的变化趋势,可利用该数据反演得到冰晶单次散射特性数据库,为全球高精度的冰云探测提供服务(Imaoka 等,2010;Letu 等,2016;李正强等,2019)。

图6 SGLI红外扫描仪和可见光和近红外辐射计Fig.6 SGLI-InfraRed Scanner and Visible and Nearinfrared Radiometer

2.1.2 机载偏振仪器

(1)RSP。为了验证高精度偏振对气溶胶和云粒子的探测优势,美国航天局设计了一种沿轨道扫描器RSP(Research Scanning Polarimeter)。RSP 是APS 的原型样机,其通道和相关参数设置和APS 完全一致(表7),相关航飞数据可为APS的在轨运行提供宝贵的参考信息(Alexandrov 等,2018)。

表7 RSP通道信息表Table 7 RSP channel information table

自1999年以来RSP 已经完成多个野外飞行实验任务,证明了高精度偏振的特殊潜能,验证了精细偏振度对云和气溶胶反演精度的提高有至关重要的作用(Alexandrov等,2020;Cairns等,2003;Chowdhary 等,2001,2005;Gao 等,2020;Stamnes 等,2018)。所获得的观测数据已广泛用于气溶胶和云参量的反演,对云微物理特性有较好的探测精度;可反演得到云光学厚度、云相态、云滴粒子有效半径、云滴粒子数浓度、冰水路径等相关云参量(Alexandrov等,2015,2018;Segal-Rozenhaimer等,2018;Sinclair等,2019)。

(2)其他机载仪器:AirHARP、SPEX airborne。AirHARP(Airborne Hyper-Angular Rainbow Polarimeter)是一种宽视场成像仪,由马里兰大学的地球与空间研究所设计,具有在670 nm通道可获取60个沿航迹视角和其他通道(440、550和870 nm)可获取20 个沿航迹视角的超角度成像能力(Fu 等,2020;Puthukkudy 等,2020)。AirHARP 可用于表征云和气溶胶的光学特性,其670 nm 通道是专门为利用偏振信号研究云虹而设计的,具有在亚千米空间分辨率下表征整视场云微物理参数的能力(McBride等,2020)。

SPEX airborne(Airborne Spectrometer for Planetary Exploration)是待发射偏振传感器SPEXone 的机载原型(2.2.4 PACE:SPEXone、HARP-2),其主要用于气溶胶探测研究(Knobelspiesse 等,2020;Van Amerongen等,2018)。

2.1.3 地基偏振仪器

太阳天空辐射计CE318 是法国CIMEL 公司生产的一种自动跟踪太阳的辐射计,可以进行云检测(李莉等,2019)。CE318 已布设在全球气溶胶自动观测网(AERONET)和中国太阳/天空辐射计观测网(SONET)等多个地基观测网络中开展长期观测(Holben等,1998)。

多角度多光谱偏振型太阳天空辐射计CE318-DP 是太阳天空辐射计中的一种类型,提供偏振测量功能。它具有的9个光谱通道和CE318仪器的中心波长相同,但与CE318 不同的是,除936 nm 通道外,其他通道均带偏振;所以与传统的强度测量相比,可以独立获得340—1640 nm 范围内8 个通道的多角度线偏振信息。

全天空成像偏振云探测器(FSIPCD)由匈牙利埃斯特拉托研发有限公司研发,该仪器可获取红光650 nm、绿光550 nm、蓝光450 nm 光谱范围内光的强度、线偏振度和可选的圆偏振度(Horváth,2014)。因为云层或多或少地改变了晴朗天空的偏振模式特征,利用全天空成像偏振云探测器探测到的偏振模式可进一步探测云层的二维分布、确定云层的相对高度(Suhai 和Horváth,2004)。这些地基偏振仪器的发展和其对应算法的开发,能为在轨偏振传感器提供产品精度验证并对气溶胶和云反演算法的改进提供服务。

2.2 计划发射偏振仪器

随着偏振数据在大气环境领域的广泛应用,科学家越来越意识到偏振传感器在探测人为气溶胶、云识别,尤其在污染水云特征识别等方面具有强大的优势。为了更好地获得多光谱高分辨率多角度偏振遥感数据,美国、中国、乌克兰等国家计划于近些年发射多台偏振载荷,从而为污染小粒子观测提供更优质的偏振数据,计划发射偏振载荷相关参量如表8所示。

表8 计划发射偏振载荷参数Table 8 Polarization load parameters to be launched

2.2.1 EPS-SG:3MI

随着POLDER 数据在大气遥感领域的广泛应用,欧盟认识到长时间偏振数据的优势,计划于EUMETSAT 极地系统—第二代(EPS-SG)任务搭载多角度—多通道—多偏振成像仪3MI(Multi-Viewing-Channel-Polarisation Imager)如图7所示。3MI 继承了POLDER(2.1.1.1 POLDER)传感器的所有优点,同时其通道分辨率和偏振通道数都有所改进,如表9(Fougnie等,2018a;Heidinger等,2016;Marbach 等,2015;Yao 等,2010)所示;此外3MI 传感器还将与EPS-SG 上搭载的其他相关仪器进行同步协同观测,可为云产品的精细反演提供更多的数据。

表9 POLDER-3与3MI波段参数对比Table 9 Comparison of band parameters between POLDER-3 and 3MI

图7 3MI整机外形图和METOP-SG Satellite A 结构示意图Fig.7 The overall appearance of 3MI and the schematic diagram of METOP-SG Satellite A

3MI能提供多光谱、多角度的大气层顶出射辐射图像,专门用于气象气候、大气成分、空气质量、气溶胶和云层特征的监测(Riedi 等,2014)。对于云,3MI 在云产品方面与POLDER 非常相似,可以获得云掩膜的产品以及云相态、云顶压强、云光学厚度和云滴谱等云参数反演产品,同时由于增强了仪器性能,从而提高了数据质量(Fougnie等,2018b)。

2.2.2 OTB-2:MAIA

由于POLDER、DPC等多角度偏振数据在细颗粒物探测方面具有广泛的应用前景,以及A-Train系列卫星的协同观测在联合反演产品方面表现出独特的优势,美国NASA提出在多角度卫星的基础上通过结合偏振测量和扩展光谱范围,将多光谱、偏振和多角度能力结合到一个单一的集成成像系统中。以Terra 卫星上的多角度成像光谱仪(MISR)为基础开发了多角度气溶胶成像仪MAIA(Multi-Angle Imager for Aerosols),提高了MISR 对大气颗粒物成分的灵敏度;MAIA 预计于2022年发射,并免费公开气溶胶和颗粒物的观测数据(Liu 和Diner,2017)。

MAIA 旨在研究可吸入颗粒物对健康的不利影响,其继承于MISR 的特点使其可以对云进行有效观测,位于强水蒸气吸收特征中的谱带可以增强卷云识别能力,氧气A 波段附近的通道用可以用于反演云顶压强,多角度观测数据可以重建云的三维形态(Diner等,1998,2005,2018a,2018b;Martonchik等,2009;Segal-Rozenhaimer等,2018)。

2.2.3 PACE:SPEXone、HARP-2

为了提高对海洋生态系统的认识,美国宇航局将于2024年发射浮游生物—气溶胶—云—海洋生态系统计划PACE(Plankton,Aerosol,Cloud and ocean Ecosystem)(Rietjens等,2019)。其上将携带两种具有偏振能力的载荷:紧凑型五角分光偏振计SPEXone 和超角度彩虹偏振仪-2 HARP-2(机载版本见2.1.2 其他机载仪器:AirHARP、SPEX airborne)。

SPEXone 是在荷兰应用科学研究组织(TNO)的支持下,由荷兰空间研究所(SRON Netherlands Institute for Space Research)和空客国防与荷兰航天公司(Airbus Defense and Space Netherlands)合作开发的。将通过5 个视角(±57°、±20°、0°)对光谱范围385—770 nm 内的辐照度和偏振度进行超光谱测量,其对线偏振度(DoLP)的测量精度为0.3%(Hasekamp等,2019)。

HARP 是一种宽视场成像仪,其一代载荷于2019-11-02搭载在CubeSat平台上发射到国际空间站,可以以低于5 km 的高度进行宽条带偏振云反演。然而由于受到航天器存储和下行链路能力的严重限制,每天只能看到一个场景。在即将发射的改进版HARP-2 中,该问题已被解决,同时HARP-2 继承了宽条带偏振云反演能力,将与SPEXone协同观测,为未来的地球科学空间任务提供高分辨率、超角度成像的偏振数据(Frouin 等,2019)。

2.2.4 Aerosol-UA:ScanPol、MSIP

虽然美国NASA 的Glory 任务发射失败了,其设计目的仍有借鉴意义,乌克兰国家天文台(MAO)基于Glory 计划启动了大气气溶胶探测计划Aerosol-UA,将多通道扫描旋光仪ScanPol 和多光谱广角成像仪MSIP 搭载于同一平台进行协同观测,目的是利用气溶胶偏振传感器监测大气中对流层与平流层中气溶胶和云的主要特征及其时空分布,该任务将于2022年发射。

ScanPol 是一种连续扫描偏振传感器,继承了APS 的设计原理,与APS 相比,ScanPol 中的光谱通道数量减少到6 个,但增加了370 nm 处的光谱通道(Milinevsky 等,2016)。ScanPol 能够在很宽的光谱范围内收集从云中散射的偏振辐射数据,从而对全球的云层分布进行遥感研究,并具有分离卷云信号的能力(Syniavskyi 等,2015)。

MSIP是一种多光谱广角偏振成像仪,可以获取测量区域内有关大气状态(云分布)和地表的信息。两个载荷的协同观测将对研究自然和人为气溶胶对气候与生态的影响作出贡献(Milinevsky等,2015)。

3 云偏振遥感经典方法

云层是太阳入射辐射和地球表面出射辐射的必经屏障,并在全球大气能量交换、水汽循环和全球辐射收支平衡中起着不可忽视的作用,进而影响着全球的气候变化。当前已有较多学者通过传统的遥感观测手段从辐射强度、空间几何特性等方面对云进行了一系列研究,然而存在着云的检测精度、物理特性与光学特性探测效果差等问题。通过对偏振数据在云检测、云物理特性、云光学特性方面展开的研究进行论述,梳理偏振仪器的云研究(表10),从而揭示云偏振遥感研究发展历程以及偏振探测的优势,探究偏振数据的应用在云相关研究中起到的重要作用。

表10 云偏振遥感的经典算法Table 10 Classical algorithms for cloud polarization remote sensing

3.1 云检测

现有的云检测算法非常多,主要基于光谱通道、多角度、几何或纹理特征进行云像元的识别(Christodoulou 等,2003;Li 等,2004,2019;Wang 等,2016),但反演结果存在问题,尤其在下垫面为亮地表时的云检测以及云边缘的检测方面存在较大问题(Mahajan 和Fataniya,2020)。由于偏振信息增加了信息的维度,所以能够提供更好的解决方案,在云检测方面更有优势(麻金继等,2021)。因为偏振载荷在波段设置上不同于传统的光学载荷,所以其云检测算法通常需要结合偏振与非偏振信息考虑(陈震霆等,2018)。POLDER 云检测算法最早由Bréon 与Buriez 提出并由大量学者持续研究发展,算法大多是通过对每一个观测方向的表观压强测试、865 nm 通道或670 nm 通道反射率测试、443 nm 通道偏振反射率测试、865 nm 偏振通道的测试来识别云像元,可充分使用云的偏振特性进行云识别(Bréon 和Colzy,1999;Buriez 等,1997)。至今,该算法相对完整,但仍然存在一些问题:首先,由于POLDER 空间分辨率较低导致云检测精度受到限制;其次,未考虑到不同的下垫面问题和时间空间变化的问题;此外,该算法中对检测结果缺失质量评价,缺乏云置信度产品。

高分五号卫星搭载的DPC 传感器具有比法国POLDER传感器更高的光谱空间分辨率和更多的观测角,从而在云遥感方面提高云检测精度并反演得到更多的云参量。国内相关学者(李超,2019;Li 等,2021)利用DPC 数据在考虑了不同的下垫面及时空变化的情况下,构建了全球的云检测算法;该算法首先使用矢量辐射传输方程模拟了不同时空、不同大气模式和不同下垫面的相关晴空和云的辐射值,正演获得云识别的动态阈值;与此同时,为了提高特殊场景及亮地表上空云检测结果的准确性,该算法还构建了全球月平均地表反照率数据库和全球冰雪覆盖数据库;此外,在上述云判识算法基础上,增加了晴空判识方法,从而增加了DPC 云判识结果的准确性评价;弥补POLDER云检测结果质量评估标准的缺失。该研究通过获取的云判识结果构建了DPC 云置信度评价体系,将获得的云识别结果分为了高置信度云、可能性云、高置信度晴空和可能性晴空4类;与相关经典传感器的产品结果对比,该方法云检测结果精度为93.06%(CALIPSO)和85%(MODIS)(李超,2019)。伟乐斯(2020)在参考POLDER云检测算法的基础上,利用表观压强、多波段反射率、NDVI、偏振反射率等信息获取DPC 的云检测结果,识别结果与MODIS 产品进行对比,云命中率平均值相较于MYD06 云掩膜产品高出了13.50%,同时错误预报率平均值相差仅为3.561%。此外伟乐斯等(2021)还使用高精度的云—气溶胶激光雷达的垂直特性掩码数据(CALIPSO-VFM)作为训练样本,针对DPC 的多偏振多角度数据,使用随机森林的方法,利用机器学习训练得到云识别结果。

此外Wang等(2017)使用CAPI数据在870 nm和670 nm 通道的反射率及其比值来进行云识别;基于MAPI 数据,Gao 等(2019)使用陆地归一化植被指数(NDVI)、670 nm 波段反射率及观测结果与模拟结果的差异作为3个判断指标来进行云识别,与POLDER 官方公布的云掩膜产品的一致性约为90%。

虽然上述卫星遥感可以获得覆盖面积广的云探测资料,但是仍受到偏振分辨率不足的限制,因此有必要利用数据质量更高、偏振精度更高的地基观测产品来补充和验证卫星遥感产品,从而验证卫星反演算法。李莉等(2019)开发了一种利用CE318-DP进行云检测的方法,使云检测的精度基本达到激光雷达进行云检测的精度,填补了主平面偏振观测无云检测方法的空白,同时也为未来空基遥感云检测算法的验证提供了一种选择。

3.2 云参量反演

3.2.1 云相态

传统非偏探测方法主要是基于仪器的红外通道利用亮温进行云相态的识别,这种方式会受到温度反演精度的影响(Baum 等,2003;Nasiri 和Kahn,2008;Platnick 等,2003);除温度外,冰云和水云的另一个显著区别是成云粒子的形状,这一差别使得用偏振特性区分云相态优势显著,即偏振观测主要依赖于粒子的形状和尺寸,表现出球形液滴的水云粒子相对于冰云在散射角140°附近的偏振观测存在明显的峰值,通常称为主虹,并在150°附近出现多个附属虹,这也是进行多角度偏振云相态判识的主要理论依据。Goloub 等(2000)应用不同相态的云在865 nm 通道偏振特征差异,基于云在散射角60°—140°范围内的斜率和散射角140°—180°范围内的标准差不同,构建了POLDER的云相态反演算法,为云相态的识别开辟了新的算法。Riedi 等(2010)综合了常用的3 种方法:多角度偏振测量、短波红外和可见光测量、热红外测量,使用POLDER 数据和MODIS 数据协同观测得到一个综合的云相态反演算法,并可以根据3种算法的识别结果获得云相态置信度因子。

郭俊杰等(2018)利用国产机载多角度偏振成像仪(MAI)的首次观测试验数据获取了水云在140°散射角附近呈现出的主虹特征,证明了国产多角度偏振成像仪具有不依赖温度等先验知识而进行云相态识别的能力。相关学者同样基于水云的虹效应,构建了国产DPC 数据的云相态判识算法模型,将识别结果与MODIS 和CALIPSO 结果进行对比,发现DPC 水云识别的准确性为93.70%,冰云识别的准确性为85%(李超,2019;Li 等,2021)。Shang 等(2020)在POLDER 云相态算法的基础上,基于不同散射角范围的偏振特征提出了一种多重测试的DPC 云相态算法,将算法运用于全球,发现该算法反演得到的水云和冰云的分布情况与MODIS 和POLDER 具有较好的一致性,但是在混合相云的分布上存在差异。

3.2.2 云顶高度

传统利用波段辐射特性反演云顶高度的方法主要有亮温法、氧气A 吸收带法、CO2薄片法。亮温反演法主要利用3—5 μm以及10—12 μm的红外窗口的辐射数据来反演云顶温度,进而根据温度廓线信息得到云顶高度,所以它受温度廓线的影响比较大(Hollars 等,2004);氧气A 吸收带法的原理是根据反演出的氧气透过率,结合大气压强廓线得到云顶压强值,该方法对薄云的反演精度比较低(Piters等,2006;张岩等,2011);CO2薄片法是利用CO2的红外吸收波段对大气高度层次变化敏感的特点得到云顶压强,对低云的反演精度较差(林琳等,2006;Platnick等,2003)。而利用POLDER 数据进行云顶气压和高度的反演时有两种方法:第1种是其具备两个氧气A 波段,可利用氧气透过率差异来反演云顶高度(Vanbauce 等,1998),该方法与传统非偏仪器的方法一致;第2种是基于偏振辐射信息的云顶瑞利压强反演云顶高度,如Goloub 等发现机载POLDER-1 数据在865 nm 观察到的偏振可以忽略不计,而450 nm 的偏振光主要是分子散射的贡献,所以可利用两通道偏振反射率的差异获得云顶气压,进而得到云顶高度(Goloub等,1994)。

乔瑞等(2021)基于DPC数据的763 nm(带宽10 nm)和765 nm(带宽40 nm)两个氧气A吸收波段,通过多项式拟合得到云顶压强和两波段表观反射率比值的经验关系进行云顶压强反演的研究,反演结果与MODIS相比平均偏差约为98.6 hPa。

3.2.3 光学厚度

传统的双光谱技术是利用410 nm 和1650 nm或2130 nm 通道同时获取云光学厚度和有效粒子半径(Nakajima和King,1990;Nakajima等,1991),尽管该方法在AVHRR和MODIS上具有悠久的应用历史(Platnick 等,2017),但在以下几个方面存在问题:首先,410 nm 通道受瑞利散射的影响显著,会导致光学厚度反演的不确定性;其次,选择的冰晶模型不是最先进的。国际上相关学者针对于传统方法出现的问题,考虑到波长对地表反照率影响较大而对云顶反照率的影响较小的特点,基于陆地和海洋分别使用POLDER 数据的670 nm和865 nm 反射率来反演云的光学厚度,从而避免410 nm 通道受瑞利散射的影响导致反演结果的不确定(Goloub 等,1999)。李树(2018)基于POLDER-3 数据利用矢量辐射传输模型RT3 分别计算在670 nm 和865 nm 大气顶的反射率、云光学厚度和观测几何一一对应构建查找表从而进行海陆上空的云光学厚度反演;并就海洋上空的反演结果,针对GHM、ASC、SC 等3 种冰云模型与POLDER 冰云产品进行了相关性分析,得出ASC模型的反演结果与POLDER-3 的结果最为相近,间接验证了ASC 模型的有效粒子半径变化对冰云光学厚度反演结果影响较小,研究表明建议采用ASC 模型作为DPC 冰云光学厚度反演的散射模型。随后,国内其他学者在针对DPC 数据的云光学厚度反演中,陆地、海洋上空分别选取了670 nm 波段、865 nm 波段,水云和冰云模型上分别选择了标准水滴模型和ASC 模型进行云光学厚度反演(李超,2019;Li等,2021)。

3.2.4 云滴谱

在卫星遥感中获取云滴尺寸分布(CDSD)参数的方法一般有两种:其中双光谱反射率法是利用测量的可见光和近红外波段的反射率与预先计算的查找表中的反射率进行比较,从而获得有效半径CDR,但无法获取有效方差EV(Nakajima 和King,1990;Nakajima 等,1991);第2 种多角度偏振反射率法是基于副虹(散射角范围:145°—165°)偏振反射率的峰值对CDR 敏感、振幅对有效方差EV 敏感的原理,来获取CDR 和EV,但是反演要求水云具有较高的均匀性分布(Bréon 和Goloub,1998;Bréon 和Doutriaux-Boucher,2005)。Shang 等(2015,2019)在前人的基础上通过加入主虹(散射角范围:135°—145°)的偏振观测,从而提升查找表检索范围,以便改进空间分辨率至40—60 km,达到降低由于水云的非均匀性导致的反演误差目的。由于利用偏振信息可以对不同尺寸半径的粒子进行很好的区分,如Xu 等(2019),Shen 等(2021)利用偏振对小粒子敏感的特性,构建了烟雾污染云和大陆性污染积云识别的方法。

对于云滴尺寸而言,机载RSP 独特的高角度分辨率允许使用在散射角范围为135°—165°的偏振反射中观察到的彩虹结构来表征液滴尺寸(Alexandrov 等,2015)。分析彩虹结构的方法有两种:第1 种是利用参数拟合法得到液滴尺寸分布(DSD)的有效半径和方差;第2 种是彩虹傅里叶变换(RFT),这是一种非参数反演方法,不需要任何关于液滴尺寸分布函数的先验知识,并且计算速度快,通常在670 nm 和865 nm 波段观察到最佳结果(Alexandrov 等,2012)。Alexandrov 等(2018)将RFT 反演结果与RSP 现场测量结果进行比较来估计反演精度,利用偏振云弓对云滴尺寸分布(droplet size distributions DSDs)的强依赖性,反演云滴大小,然后比较了由RSP 观测得到的云滴尺寸分布(DSDs)和由云滴探针(CDP)测量得到的云滴尺寸分布;这是首次通过与相关的现场数据直接比较来验证极化液滴尺寸反演的有效性;通过提出的极化云滴尺寸反演技术的验证,证明了偏振观测的价值。

除了RSP 偏振数据之外,Knobelspiesse 等(2020)利用机载AirHARP 偏振数据,使用传统的参数拟合方法反演了云滴谱分布,达到了提高相关产品精度的效果。

3.2.5 冰水路径

利用可见光—红外成像辐射仪观测的非偏振数据反演卷云的特性,需要对卷云中的冰晶颗粒物形状等特性进行相关假设,从而达到反演卷云中冰晶的数量、半径、光学厚度和冰水路径等参量;而利用RSP 传感器获得偏振数据来反演卷云相关特征时,无需对卷云中的冰晶颗粒物的形状等特性进行假设,从而减小了卷云相关参量反演结果的不确定性,以便冰水路径(IWP)反演精度更高(Cairns 等,2003),这充分说明了偏振数据在冰云参量反演方面的优势。

4 结语

通过对近30年来偏振传感器的梳理,我们发现偏振传感器从最开始的低空间分辨率、少偏振光谱、少角度和低偏振精度能力,逐渐发展到目前偏振传感器具有中高空间分辨率、多偏振光谱、多角度和高偏振能力,再到未来高空间分辨率、多光谱、多角度与高精度偏振耦合的偏振传感器。仪器的发展方向为:(1)单载荷通过提升自身能力实现高精度多角度偏振探测。例如欧盟即将发射的3MI将在POLDER的基础上进一步改进,增加短波红外信息的探测可以使得部分传统云相关研究算法得以在新数据中运用;美国将发射的MAIA是在MISR 的基础上增加了偏振观测能力,增加信息获取的维度从而提升算法精度。(2)多载荷通过协同观测实现多光谱通道、多角度、高精度偏振融合。例如中国提出的“POSP+DPC”的“偏振交火”方案、美国与欧洲提出的“APS+3MI”的联合观测方案以及乌克兰启动的大气气溶胶探测计划Aerosol-UA,都将通过偏振多角度成像仪与高精度的偏振探测仪相结合实现联合观测,该方案也是未来偏振仪器的主要发展方向之一。

在云遥感相关研究应用中,增加偏振信号能够很好地克服部分云遥感中出现的困难:在云检测中,多通道偏振信息有利于修正云边缘误分类;在云相态的识别中,利用水云865 nm 偏振通道散射角140°附近呈现出的主虹特征,可以完成更为精确地云相态区分;与此同时,使用偏振探测能更有效地分析云层内部的微物理状态,在解析云粒子形状和尺寸分布时更为有效;除此之外,进行神经网络的训练时,在输入向量中加入线性偏振度可以提高训练精度。虽然偏振观测带来的观测目标信息量增加,使得云相关研究取得了一些进步,但是还有一些问题需在未来的研究中解决,如现实中云系统组成复杂,变化较快,但利用偏振数据进行云相关研究时涉及多层云或者复杂云系统的研究较少,未来应多考虑复合云场的偏振成像;在冰云辐射模拟过程中,由于卷云中冰晶粒子的形状和构成成分的复杂性比较大,未来需要进一步针对研究区域建立冰晶粒子散射特性数据库,来提高冰云偏振辐射模拟的精确度。随着偏振仪器精度的提升及多角度偏振仪器和高精度偏振计联合观测方案的实现,希望可以在上述偏振云遥感研究领域取得一些进展。本文虽然归纳总结了偏振传感器及其偏振云遥感的相关算法的发展历程,然而由于篇幅有限,对云偏振遥感研究领域调研还不够全面;此外针对经典方法没有对其不确定性加以描述。

综上所述,随着偏振传感器的日新月异,偏振通道越来越多、偏振精度越来越高、空间时间分辨率越来越高,可获取的偏振信号越来越多越来越复杂;其设计对应的遥感算法时,假设将会越来越少,算法将会越来越严谨,其反演精度将会越来越高。考虑到随着经验的积累、偏振传感器数量和种类的增加,相信未来可以使用多种长时间序列的偏振观测数据进行对比验证和更深入的研究,期待偏振数据在优化云参量方面会有重大进展。

猜你喜欢

偏振气溶胶反演
首个使用偏振的 超快光处理器面世
美国X射线偏振测量天文台特点分析
病毒通过气溶胶传播,还能开窗通风吗
气溶胶传播病毒概率有多大
基于红外高光谱探测器的大气CO2反演通道选择
反演变换的概念及其几个性质
基于ModelVision软件的三维磁异常反演方法
计算机仿真光注入VCSEL的偏振开关及双稳特性
气溶胶科学
探讨光的偏振现象