2. 解放军理工大学指挥军官基础教育学院, 南京 211101
2. College of Basic Education for Commanding Officers, PLA University of Science and Technology, Nanjing 211101, China
台风是发生在热带海洋上的强烈气旋性漩涡,常给受影响地区带来严重灾害.台风的生成和发展主要受不稳定大气初始扰动及海面热力条件等因素的影响,其路径和强度在短时间内均可能出现很大变化,经过地区会产生破坏力很强的大风、强降水等天气现象,因此船只、浮标等常规观测手段一般无法对台风进行捕获和跟踪观测.在气象业务中,常利用可见光和红外遥感手段,对台风温度场及云系特征进行观测和分析,估计台风发展状况和中心附近最大风速.由于云层的覆盖遮挡,可见光和红外等被动遥感均难以识别云层下海面的热力学状况,对台风在不同阶段的云水含量、降水强度、云水和水汽的分布及变化等要素探测能力有限.主动式卫星遥感是目前获取台风条件下大范围区域海面风场资料的主要手段,以QSCAT散射计为代表,可提供空间分辨率25 km,精度为1 m·s-1(风速)、15°(风向)的海面风场产品,为科学研究和天气预报服务.
全极化微波辐射计是19世纪90年代中后期才发展起来的一种新的被动微波遥感技术,它的出现 打破了被动式卫星传感器只能观测海面风速、不能观测海面风向的局限,为被动式卫星遥感开启了崭新的应用领域.2003年,在NRL(Naval Research Laboratory for the U.S. Navy)、IPO(Integrated Program Office/NPOESS)和NASA(National Aeronautics and Space Administration)等机构的共同推动下,国际上第一个载有全极化微波辐射计WindSat的试验卫星Coriolis升空,用于验证全极化微波辐射传感器的可靠性,建立相应的海面风速、风向算法,为将来全极化微波辐射计的业务化运行提供依据(Gaiser et al., 2004).Coriolis卫星轨道高度为840 km,轨道倾角98.7°,重访周期为8天,升轨过境时间为地方 时18点.WindSat在10.7 GHz、18.7 GHz和37 GHz 频率上设置了全极化通道,测量对应于目标微波辐射4个Stokes参数的亮温,其首要应用是观测海面风场;在6.8 GHz和23.8 GHz频率上的水平和垂直双极化通道主要用来观测海面温度和大气水汽含 量,为利用全极化通道反演海面风矢量提供辅助资料.
微波具有能够穿透云盖、浓雾、降雨的特点,所以利用星载微波辐射计可以探测云以下大气和海表状况,是分析台风的一种有效手段.目前为止,在无降水情况下利用被动式微波传感器进行海面风速反 演,其精度可达到1 m·s-1(Bettenhausen et al., 2006); 一旦有降水出现,由于传感器定标精度的限制和风场反演算法的不成熟,原算法无法得到稳定的回归参数,导致风场反演无法进行,因此,NRL发布的WindSat风场、海面温度等产品中未包括云层较厚、降水强度较大等区域的产品数据.为弥补这一缺陷,Thomas Meissner和Frank J. Wentz提出可利用对风速敏感、受降水影响较小的通道组合方案,对降水条件下海面风场进行反演计算(Meissner and Wentz, 2009).全极化辐射计在短时间内获得全球大面积风场的同时,还可以获得大气云水含量、水汽含量和海表温度等环境场参数,是国际上星载微波遥感技术的前沿方向.开展深入的理论和数据处理研究,对于我国将来设计和制造自己的全极化辐射计具有重要的意义.
本文通过分析WindSat观测亮温在台风区域的分布特征,对利用全极化微波辐射计进行台风区域内海面风场反演的可行性进行了理论分析和验证,并在敏感性分析的基础上确定了反演使用的辐射计通道.在此基础上,设计和构建反演算法,在NRL针对WindSat发展的前向辐射传输模型FM(Forward Model)中加入了降水对大气透过率的影响,并利用台风区域实测数据对海面发射率经验校正系数进行了计算,为进行海面风场反演奠定基础.为提高海面风向反演准确度,本文设计了环境影响经验修正函数,对第三和第四Stokes通道亮温进行环境参数影响经验修正,突出海面风向的作用.在海面风场反演中,分别使用基于统计的多元线性回归算法和最大似然估计法,实现了海面风速和风向等其他物理要素的反演计算,并将海面风场反演结果与QSCAT散射计风场产品进行了对比.
2 数据来源通过WindSat数据处理技术组(Cal/Val Science Team),我们获得了JPL(Jet Propulsion Laboratory) PO.DAAC(Physical Oceanography Distributed Active Archive Center)下发的2003年11月至2004年12 月的WindSat传感器数据SDR(Sensor Data Record)、 环境数据EDR(Environmental Data Record), 数据版本TDR(Technical Data Report)为146AFBBDA.其中,SDR为定标、定位、极化旋转校正等预处理后的传感器亮温,EDR为对应于SDR观测点的产品数据,主要包含晴空条件下海面风矢量产品,还包括其他产品如海面温度、大气柱中水汽总量、大气柱中云液态水总量和降水强度等.为了验证本文海面风场反演算法的准确性,选取PO.DAAC提供的QSCAT level 2B海面风场产品进行对比分析,两颗卫星间采用60 min的时间窗口和25 km的空间窗口进行数据匹配.
依据中国台风网报告,2004年我国以东太平洋面上生成热带风暴共30次.其中,0414号台风“云娜”风力强、雨量大、影响范围广,具有典型代表性.“云娜”于2004年8月6日上午在西太平洋关岛西北的洋面发展成热带低压,移动过程中强度不断增强,10日升级为台风,12日在浙江温岭登陆时中心气压为950 hPa,中心最大风力达12级.本文研究使用了2004年8月10日WindSat在“云娜”上方过境时的观测数据(轨道号为08262)进行个例研究,分析了反演结果准确度及误差空间分布特征.在此基础上,本文针对2004年出现的台风过程逐次进行了分析计算,并将海面风场反演结果与QSCAT产品作对比,以验证反演算法的精度及可靠性.
3 台风境内风场反演可行性分析台风区的近海面风场按风速大小可分为三个区域:外围,自台风的边缘向内直到最大风速区,风速分布呈阵性,且风速向内增大;中圈,围绕台风眼的最大风速区,通常与围绕台风眼的云墙区相重合;内圈,在台风眼区,风速向中心迅速减小.海面风速的不均匀分布会影响海面粗糙度,使海面发射率发生变化,反映在全极化微波辐射计测量结果上,表现为辐射计观测亮温的不均匀分布(图 1).另外,由于风向引起的海面粗糙度的各向异性,表现为海面反射率的各向异性,因此在不同的相对风向角条件下卫 星测量的全极化辐射值是不同的.台风的切线风向沿台风一周变化约为360°,某区域和另一区域的相对风向角差别可达180°,根据电磁辐射理论,这使得目标电 磁辐射的第三和第四Stokes参数 T3和T4 出现正负符号变化,WindSat实测数据也证实了这一点(图 1).
![]() |
图 1 WindSat 10.7GHz各极化通道亮温图像 (a)垂直极化通道亮温;(b)水平极化通道亮温;(c)第三Stokes通道亮温;(d)第四Stokes通道亮温. Fig.1 Brightness temperature at 10.7 GHz from WindSat(a) V-pol brightness temperature; (b) H-pol brightness temperature; (c) The 3rd Stokes parameter; (d) The 4th Stokes parameter. |
在图 1所示的WindSat 10.7 GHz各极化通道亮温图像中去掉了陆地、海岸和近海岸数据,仅保留下垫面标识为海洋的亮温图像,这是由于近海岸水体和远离海岸的水体的运动学特征存在差异.受海岸轮廓及海底大陆架对水体运动的影响,在近海岸处无法利用全极化微波辐射计进行海面风场反演.由图 1可以看出,WindSat水平极化通道亮温 Th 的分布特征与垂直极化通道亮温 Tv 类似,数值略小.从 Tv 和 Th 的分布图像可以清晰地看到:在台风风速逐渐增大的外围区,WindSat观测亮温随着海面风速的不连续变化而变化;在围绕台风眼的中圈,由于高风速引起的海面粗糙度很高,WindSat观测亮温值也很高;在台风眼区域,WindSat观测亮温随海面风速的减小而迅速减小,台风眼的位置在亮温图像上明确地被显示出来.由10.7 GHz交叉极化通道(第三和第四Stokes极化通道)亮温图像可以看出,随相对风向角的变化, T3和T4 的符号及大小存在周期性变化.由此可见,各极化通道的亮温特征包含了海面风速和海面风向信息,可通过选择观测通道,构建反演算法,将台风区域海面风场反演出来.
4 台风境内风场反演个例 4.1 反演流程目前为止,利用全极化微波辐射计观测亮温进行海面风场反演的算法可以分为三类:统计算法、半统计算法和物理算法.统计算法是从辐射计亮温和现场测量地物参数推导出经验关系式,例如多元线性回归算法或改进D矩阵算法(Goodberlet et al., 1990),可以对除海面风向外的其他物理因素如海面风速、海面温度、大气水汽等进行反演,为海面风向的反演提供必要的环境场.半统计算法又称物理统计法,与统计算法的不同之处在于回归时使用辐射传输模型计算来仿真亮温,例如最大似然估计法MLE(Maximum Likelihood Estimate),一般用来进行海面风向反演(Meissner and Wentz, 2005),也可同时进行海面风向和风速的联合反演,但进行联合反演时风向和风速的反演精度会互相制约.物理反演算法是Wentz(1992)提出的,其本质是非线性方程组的求解问题,反演的准确及效率取决于地物模式的准确性和方程组解法的有效性.
在以上三种方法中,统计方法在大气参数和海面温度的反演方面已经获得广泛应用.本文中将利用统计方法反演大气水汽含量、云中液态水含量、海面温度和海面风速4个物理量.在海面风向的反演方面,迄今为止多采用物理统计法进行,即在计算中利用全极化微波辐射传输正演模型来模拟亮温进行迭代.由于正演模型输入需要大气和海面各部分参数,所以在使用物理统计法反演海面风向之前,需要先对模型输入量进行反演或估计.另外,为了去除大气和下垫面对海面风向测量通道 T3和T4 的影响,提高海面风向反演精度,在进行海面风向反演之前对辐射计观测亮温加入了大气和海面环境影响经验修正.反演的整体流程如图 2所示.
![]() |
图 2 反演算法流程图 Fig.2 Flow scheme of retrieval algorithm |
具体来讲分为以下几个步骤.首先,利用统计方法反演计算海面风速、海面温度、大气水汽和云中液态水含量;第二步,对第三和第四Stokes通道亮温进行大气和海面环境影响经验修正,突出海面风向的作用;第三步,将反演结果作为正演模型输入量的一部分,利用最大似然估计法反演计算海面风向;第四步,采用中值滤波的方法,在几个最大似然解中选取风向的最优估计值;最后,将风场反演结果和大气 参数反演结果输入正演模型,结合最大似然估计法进行海面风向迭代反演,若相邻两次风向迭代反演结果小于规定阈值,即作为风场最终反演结果输出.
4.2 基于敏感性分析的通道选择方案全极化微波辐射计观测亮温包含海面风场信息,可对全极化微波辐射计相对于大气和海洋等物理参数的敏感性进行分析,考察辐射计的哪些通道对哪些物理因素较为敏感,为反演海面风矢量、海面温度、水汽等重要物理因素时通道的选取提供理论支持,并为去除大气水汽和云中液态水、降水强度等在海面风矢量反演中的不利影响,提高海面风矢量反演精度奠定基础(严卫和王蕊, 2007).
同样的背景条件下,不同观测频率上环境要素对海面风向信号的干扰程度是不同的.由于高频率通道亮温对海面风速的敏感性要高于低频通道(Wentz, 1983),所以以往的研究者通常使用18.7 GHz和37 GHz进行无降水情况下风速的反演.但是,微波辐射计亮温除了受到风引起的海面粗糙度的影响外,还受到海面温度、大气状况等因素的影响,高频通道的大气透明度明显低于低频通道.在降水条件下进行海面风场反演时,采用高频通道将带来较大反演误差.
以水汽含量对微波辐射计接收亮温的影响为 例,假设背景条件为6 mm·h-1中等强度降水,海 面温度301 K,垂直气柱上云中液态水含量为0.26 mm, 海面风速13 m·s-1,观测方向与实际风向的相对风向角为45°,当大气水汽含量发生变化时,各通道将受到不同程度的影响(图 3).
![]() |
图 3 水汽含量变化对10.7 GHz、18.7 GHz和37 GHz通道的干扰情况(a)垂直极化通道;(b)水平极化通道;(c)第三Stokes通道;(d)第四Stokes通道. Fig.3 Disturbed condition at the 10.7 GHz、18.7 GHz and 37 GHz from vertical columnar water vapor(a) V-pol; (b) H-pol; (c) The 3rd Stokes parameters; (d) The 4th Stokes parameters. |
如图 3a、3b所示,观测频率较高时, Tv和Th 随水汽含量的增加将明显增大,海面风速信号在一定程度上会被掩盖.通过对台风区域WindSat不同频率实测亮温的空间分布特征进行分析(图 4),也能得出同样的结论.
![]() |
图 4 (a)WindSat 10.7 GHz垂直极化通道亮温图像;(b)WindSat 18.7 GHz垂直极化通道亮温图像;(c)WindSat 37 GHz 垂直极化通道亮温图像 Fig.4 (a) V-pol brightness temperature at 10.7 GHz from WindSat; (b) V-pol brightness temperature at 18.7 GHz from WindSat; (c) V-pol brightness temperature at 37 GHz from WindSat |
通过对比分析WindSat 10.7 GHz、18.7 GHz和37 GHz Tv 图像,可知在云覆盖或降水条件下,高频率通道卫星观测亮温受到大气环境的影响更为明显,37 GHz上 Tv 在台风眼区周围因云和降水而几乎饱和(图 4c),若利用此通道进行海面风速反演,反演误差必然较大.因此,在进行高水汽含量的台风区域风速反演时,在同水平定标精度等情况下,应尽量避免选取使用高频通道.基于上述分析结果,本文选择利用6.8 GHz和10.7 GHz通道进行海面风速反演.
对于海面风向观测而言,虽然相对风向角未变化,但随着大气水汽含量的变化,交叉极化通道的观测亮温 T3和T4 仍会产生很大变化,并且变化幅度随观测频率的增高而增大(图 4c,3d).为了避免高频通道上云和降水带来的影响,本文利用10.7 GHz单频率交叉极化通道结合MLE法进行海面风向反演,并设计了大气和海面环境影响经验修正函数,去除该通道上环境大气的亮温贡献,突出海面风向的作用,以提高反演精度.
4.3 反演算法实现 4.3.1 扩展辐射传输正演模型基于物理统计法进行海面风向反演时,需要使用全极化微波辐射传输正演模型来模拟亮温进行迭 代,因此正演模型的准确程度直接关系到最终 反演结果的准确度.以往常用的微波辐射传输模式如MODTRAN(MODerate resolution atmospheric TRANsmission)、LBLRTM(The Line-By-Line Radiative Transfer Model)、CRTM(Community Radiative Transfer Model) 等均不包含对目标微波辐射传输的全极化描述.为了适应WindSat观测海面风场的任务需求,NRL构建了晴空条件下针对WindSat的前向辐射传输模型FM(Forward Model),计算物理统计反演法中所需的通道仿真亮温.2005年11月至今, EUMETSAT NWP SAF (Satellite Application Facility on Numerical Weather Prediction)发布了快速辐射传输模式RTTOV(Radiative Transfer for TOVS) 8.7起系列版本,包含对全极化微波辐射计WindSat 的模拟计算.本文研究使用的辐射传输正演模型基于参数化模型NRL-FM(Wang et al., 2005),为使辐射传输模型能够在非晴空条件下使用,在大气部分加入了降水对大气透过率的影响,扩展后的正演模型框架如图 5所示,输入量为海面温度、海水盐度、入射角、海面风向、海面风速、观测频率及极化方式、水汽含量、云中液态水含量、降水强度等,输出量为 星载全极化微波辐射计仿真亮温.在此基础上,利用 选用的台风区域内WindSat EDR、SDR数据及LBLRTM 模式、RTTOV模式等对模型中的部分参数重新进行了拟合.
![]() |
图 5 前向辐射传输模型框架示意图 Fig.5 Diagram of Forward Model |
按照极化辐射传输理论,星载全极化微波辐射计接收到的亮温可表示为(Wang et al, 2005):
在1~37 GHz,即使有云或中等强度的降水,大气也是相当透明的,吸收成分主要为O2和H2O,大气透过率 τ 可用垂直方向上总的氧气吸收系数 AO、 水汽吸收系数 Av、 云中液态水吸收系数 AL 和降水衰减系数 AR 在观测角 θ 路径上的投影量来计算:
eφ 是风向 φ 引起的海面发射率的变化(Wang et al., 2005):
为减小模型整体误差,参考Michael等(2004)使用2003年9月至2004年2月的WindSat亮温资料研究发展的对风成海面发射率的修正方案(Wang et al., 2005),对垂直极化和水平极化通道采用各向同性的经验修正形式,对交叉极化通道采用谐波形式的经验修正(式(13)),系数 cei(i=0,1,2,3) 使用2004年8月份的WindSat亮温数据重新进行拟合得到(表1).
![]() |
表1 海面发射率经验校正系数 Table 1 Empirical correction coefficients to sea surface emissivity |
在利用星载微波辐射计垂直和水平极化通道亮温进行大气参数的反演时,多采用线性统计回归算法(Wentz, 1999):
![]() |
表2 环境参数反演系数 Table 2 Environmental parameters inversion coefficient |
![]() |
图 6 环境参数反演结果(a)海面温度;(b)大气水汽含量;(c)云中液态水含量;(d)降水强度. Fig.6 Environmental parameters retrieval results(a) Sea surface temperature; (b) Water vapor; (c) Cloud liquid water; (d) Rain rate. |
从反演结果可以看出,台风区域海面温度因降水的作用在整体上低于周围海域,但在台风中圈存在一个明显的暖涡结构(图 6a箭头所示),并且台风周围大范围水汽含量很高(图 6b),在底层风场的作用下水汽不断地被输送到台风内部,此刻台风正处于旺盛发展阶段,强度将会进一步加强.
除用于对台风状况进行定性分析外,环境场要素反演结果的另一个重要用途是对交叉极化通道观测亮温 T3和T4 进行大气和海面环境影响经验修正.无论理论分析还是实际验证均表明,WindSat各频率交叉极化通道测量的亮温除了随海面风向的不同而发生变化以外,还随着大气水汽含量 V、云中液态水含量L、海面温度Ts、降水强度R 等不同而发生变化.因此, T3和T4 不仅包含了海面风矢量的信息,还包含了大气和海洋的其他物理信息,而后者对于海面风矢量的反演来说是一种干扰因素.为了消除其他物理因素的影响,在此设计一个修正函数 f(Ts,V,L,R),将T3和T4 加上该函数,可以达到大气和海面修正的效果.
![]() |
表3 第三和第四Stokes极化通道亮温修正系数 Table 3 The 3rd and 4th Stokes temperature correction factors |
进行海面风速反演,可考虑单层大气模型,假设 Ω=1,TC=0, 将大气温度用等效温度Teff表示,认为 Ts=Teff, 则星载全极化微波辐射计垂直通道和 水平通道观测亮温可以由方程(1)和(4)推导简化为:
![]() |
图 7 (a)海面风速反演结果;(b)反演风速与QSCAT产品比较 Fig.7 (a) Sea surface wind speed retrieval results; (b) Comparison retrieved wind speed with QSCAT |
利用WindSat亮温进行海面风向的反演一般采用物理统计法进行,在各种物理统计反演方法中,MLE法具有更深刻的理论基础.以前MLE法常用于散射计海面风场的反演(解学通等, 2005),但是通过研究发现,这种方法同样可以用来进行全极化微波辐射计风场的反演.
WindSat接收到的亮温 T *B 由两部分组成:
另一方面,采用正演模型对真实值 T 0B 进行模拟时也存在误差,记为模型误差 ε M,则 T 0B 和模型预测值 T MB 之间的关系为:
理论上, JMLE 的全局最小值应对应于真实风向的最佳估计值.但实际上,由于正演模型本身的局限性及辐射计各种测量噪声的影响,真实风向的最佳估计值可能对应于 JMLE 的局部极小值.风向反演结果去模糊考虑风场的连续性,使用中值滤波的方法进行,从若干个 JMLE 的局部极小值对应的风向中选取出风向的最佳估计.使用未经环境修正的SDR进行海面风向反演时,风向反演结果与散射计QSCAT风场产品相比均方根误差为30.71°,使用 修正后的SDR进行反演,风向反演精度提高到18.56°, 风场最终反演形态见图 8.
![]() |
图 8 (a)海面风矢量反演结果;(b)QSCAT风场 Fig.8 (a) Sea surface wind vector retrieval results; (b) QSCAT wind vector at nearly same time from same area |
通过对反演误差在台风区域的空间分布情况(图 9)进行分析,可以看出海面风速和海面风向反演误差的空间分布特征是不同的.对海面风速而言,台风中圈对应的区域反演误差较大,这主要是由于台风中圈最大风速区通常与围绕台风眼的云墙区相重合,云墙和强降水表现出的亮温增强信号与高风速引起的亮温增强信号混淆在一起,导致海面风速反演结果与实际值出现较大偏差.对海面风向而言,反演误差较大的区域主要集中在相对风向角0°~ 10°、160°~200°和350°~360°区域内,这是由于交叉极化通道亮温 T3和T4 在该角度区间内测量值一般小于0.2 K,观测值被测量噪声淹没,导致无法反演出正确的海面风向.
![]() |
图 9 (a)海面风速反演误差的空间分布情况;(b)海面风向反演误差的空间分布情况(图中红色X标志台风眼位置) Fig.9 (a) The spatial distribution of sea surface wind speed retrieval error; (b) The spatial distribution of sea surface wind direction retrieval error (The red X marks the eye of the typhoon location) |
为验证本文算法的精度及可靠性,我们依据中国台风网对2004年度台风的837个时次记录,共寻 找到1013轨WindSat过境数据(包含升轨和降轨数据),其中851轨因WindSat和QSCAT在以台风眼为中心的15°×15°台风区域内匹配观测数据量小于1000未采用,实际使用量为162轨(表4),对各台风过程海面风场反演误差的统计情况见表5.
对2004年各台风过程进行分析计算,并将反演结果与QSCAT风场产品进行对比,海面风速和海面风向反演的均方根误差分别可达1.64 m·s-1和18.02°.在参与计算的数据资料中,90%左右数据点风速反演误差小于2 m·s-1、风向反演误差小于20°.通过对比分析可知,在海面风向反演计算前,采用本文设计的经验修正函数对 T3和T4 通道观测亮温进行环境影响修正,可使风向反演精度提高10°以上.
![]() |
表4 各台风过程风场计算使用的WindSat数据来源轨道编号 Table 4 WindSat data source track ID used in each typhoon process of wind vector calculation |
![]() |
表5 各台风过程海面风场反演结果与QSCAT对比误差统计 Table 5 Root mean square errors of sea surface wind speed and direction retrieval results from each typhoon process in contrast with QSCAT wind productions |
与传统的双极化微波辐射计相比,全极化微波辐射计的突出特点在于其可以提供包含海面风向在内的海面风场矢量信息,打破了传统微波辐射计只能观测海面风速的局限,是被动微波遥感领域新的发展方向.全极化微波辐射计相对于可见光、红外遥感及传统双极化微波辐射计而言,可以穿透云盖和降水,获取海面风场资料,为在台风等恶劣天气条件下进行海洋遥感提供了有效方法.
本文通过选取6.8 GHz和10.7 GHz低频通道 组合形式,基于交叉极化通道亮温 T3和T4 的风向特征,利用全极化微波辐射计WindSat实测亮温,实现了台风区域海面风场反演,90%左右数据点能达到风速反演误差小于2 m·s-1、风向反演误差小于20°的要求.反演结果可用于对台风底层海面附近最大风速及风场分布情况进行定量分析,以确定台风强度、发展状况等信息,也可作为数值预报模式的有效输入量,弥补台风区域其他观测手段的对于海面风场要素观测的局限,提高危险天气的观测和预报能力.
相对晴空条件而言,在台风等复杂天气条件下进行气象要素的微波遥感探测是困难的.随着云滴粒子和降水粒子半径的增大,大气对微波的衰减更为明显,忽略大气散射带来的误差将越来越大,最终影响产品精度,造成台风中圈的风速反演结果与实际值出现较大偏差.另外,由于云和降水等因素的影响,在复杂天气条件下准确地进行微波辐射亮温的正演仿真是很困难的,在特定的相对风向角区间内交叉极化通道的观测值将被测量噪声和模型噪声掩盖,使海面风向反演结果误差增大.本文通过设计经验修正函数,对第三和第四Stokes通道观测亮温进行环境影响修正,在一定程度上减弱了环境因素对海面风向反演精度的影响,使风向反演精度提高10°以上,为将来提供高质量的全极化微波辐射计海面风场产品奠定了基础.
致谢 感谢PO.DAAC提供的WindSat实测数据支持了本文的研究.[1] | Bettenhausen M H, Smith C K, Bevilacqua R M, et al. 2006. A nonlinear optimization algorithm for WindSat wind vector retrievals. IEEE Transactions on Geoscience and Remote Sensing, 44(3): 597-610, doi: 10.1109/TGRS.2005.862504. |
[2] | Gaiser P W, Twarog E M, Li L, et al. 2004. The WindSat space borne polarimetric microwave radiometer: sensor description and mission overview. // Proceedings of the 2004 IEEE International Geoscience and Remote Sensing Symposium (IGARSS'04, Vol.1). Anchorage, AK: IEEE, doi: 10.1109/IGARSS.2004.1369039. |
[3] | Goodberlet M A, Swift C T, Wilkerson J C. 1990. Ocean surface wind speed measurements of the Special Sensor Microwave/Imager (SSM/I). IEEE Transactions on Geoscience and Remote Sensing, 28(5): 823-828, doi: 10.1109/36.58969. |
[4] | JPL. 2006. WindSat Data Products Users' Manual. Version3. |
[5] | Liu J Y, Wang Z Z, Yin X B, et al. 2007. An ocean wind retrieval algorithm for WindSat microwave Polarimetric radiometer. Remote Sensing Technology and Application (in Chinese), 22(2): 210-215. |
[6] | Meissner T, Wentz F. 2005. Ocean retrievals for WindSat: Radiative transfer model, algorithm, validation. // Proceedings of MTS/IEEE (OCEANS 2005). Washington, DC: IEEE, 130-133, doi: 10.1109/OCEANS.2005.1639750. |
[7] | Meissner T, Wentz F J. 2009. Wind-vector retrievals under rain with passive satellite microwave radiometers. IEEE Transactions on Geoscience and Remote Sensing, 47(9): 3065-3083, doi: 10.1109/TGRS.2009.2027012. |
[8] | Steven S P, Rose L A, Gaiser P W, et al. 2004. Recent measurements of the microwave emissivity of foam: Effects on WindSat brightness temperatures. WindSat Cal/Val and Science Meeting, Solomons, MD. |
[9] | Wang N Y, Chang P S, Bettenhausen M, et al. 2005. WindSat physically based forward model: atmospheric component. // Proc. SPIE 5656, Active and Passive Remote Sensing of the Oceans, 104, doi: 10.1117/12.578811. |
[10] | Wentz F J. 1983. A model function for ocean microwave brightness temperatures. Journal of Geophysical Research: Oceans (1978—2012),88(C3): 1892-1908, doi: 10.1029/JC088iC03p01892. |
[11] | Wentz F J. 1992. Measurement of oceanic wind vector using satellite microwave radiometers. IEEE Transactions on Geoscience and Remote Sensing, 30(5): 960-972, doi: 10.1109/36.175331. |
[12] | Wentz F J. 1999. AMSR Ocean Algorithm, Remote Sensing System. RSS Tech. Report 110398. Version 2. |
[13] | Xie X T, Fang Y, Chen X X, et al. 2005. Research on numerical wind vector retrieval algorithm based on maximum likelihood estimation. Geography and Geo-Information Science (in Chinese), 21(1): 30-33, doi: 10.3969/j.issn.1672-0504.2005.01.009. |
[14] | Yan W, Wang R. 2007. Wind vector sensitivity analysis for Polarimetric microwave radiometer (in Chinese). // Jiangsu and Anhui Provinces Academic Seminar in Atmospheric Sounding, Remote Sensing and Electronic Technology. Anqing, 157-161. |
[15] | 刘璟怡, 王振占, 殷晓斌等. 2007. 一种针对Windsat极化辐射计的海面风场反演方法. 遥感技术与应用, 22(2): 210-215. |
[16] | 解学通, 方裕, 陈晓翔等. 2005. 基于最大似然估计的海面风场反演算法研究. 地理与地理信息科学, 21(1): 30-33, doi: 10.3969/j.issn.1672-0504.2005.01.009. |
[17] | 严卫, 王蕊. 2007. 全极化微波辐射计海面风矢量敏感性分析. // 第四届苏皖两省大气探测/环境遥感与电子技术学术交流研讨会. 安庆, 157-161. |