330 likes | 562 Views
卫星海洋遥感导论 An Introduction to Satellite Oceanic Remote Sensing. 第三部分 海洋参数的遥感反演 第十章 海面风场遥感. 武汉大学 遥感信息工程学院 刘良明 教授 lm_liu69@sohu.com. 第十章 海面风场遥感.
E N D
卫星海洋遥感导论 An Introduction to Satellite Oceanic Remote Sensing 第三部分 海洋参数的遥感反演 第十章 海面风场遥感 武汉大学 遥感信息工程学院 刘良明 教授 lm_liu69@sohu.com
第十章 海面风场遥感 LiuLM Satellite Oceanic Remote Sensing 1 海面实测风场资料获取方法 1.常规观测系统 :船舶、海上浮标及沿岸和岛屿气象台站; 2.测量海面风场的传感器:微波散射计、微波辐射计、高度计和合成孔径雷达。
第十章 海面风场遥感 LiuLM Satellite Oceanic Remote Sensing 2 10.1 微波散射计测量海面风场 10.1.1 微波散射计的发展历程 1.1973 年至1974年,Skylab实验室试用了第一个太空携带式的散射计S-193 - 证明了在理论上还是技术上都是可行的 。 2. 1978 年发射的Seasat-A海洋卫星载有散射计用于测量海面风场。 3. ESCAT散射计(ERS-1 和ERS-2 卫星)、Ku波段NSCAT 散射计(ADEOS卫星 )、SeaWinds散射计 (QuikSCAT卫星 )的应用。
10.1 微波散射计测量海面风场 LiuLM Satellite Oceanic Remote Sensing 3 10.1.2 微波散射计的测风原理 微波散射计通过测量风引起的粗糙海面对微波的后向散射特性来推算风场。在海面上,毛细波叠加在重力波上,风的变化引起海表面粗糙度的变化,使接收到的后向散射随之变化。根据后向散射与风矢量之间的相关模式,经过地球物理定标后就能得出海面风场。 从散射计截面反演风场参数包括正反演两个问题。正演问题是建立雷达截面与入射角、风速、风向关系的精确地球物理模型,反演问题则是从模型反演风速和风向。
10.1.2 微波散射计的测风原理 LiuLM Satellite Oceanic Remote Sensing 4 雷达后向散射系数б°估算 散射计接收到的回波脉冲功率等于: 其中,PS+N是回波信号与噪声的混合等效功率;PN是在脉冲间对噪声单独积分后,测出的接收机噪声和无线噪声平均功率。 根据Pr,由下式可估算出雷达后向散射系数б°:
10.1.2 微波散射计的测风原理 LiuLM Satellite Oceanic Remote Sensing 5 功率谱密度 散射计发射脉冲的后向散射主要来源于那些满足Bragg共振条件的表面小尺度波,后向散射截面取决于这些小尺度波的功率谱密度。例如,相对于海面以上10米处的风速有: 根据风速与小尺度波的功率谱的关系: 导出单位面积后向散射系数的Bragg表达式为:
10.1.2 微波散射计的测风原理 LiuLM Satellite Oceanic Remote Sensing 5 后向散射系数与海面风速关系 由(10-1-5) 式可知,海面后向散射б°随摩擦风速u.线性增长。而摩檫风速与海面以上高度为z的风速U(z)的关系可通过Monin-Obukhow方程式来确定。 L是Monin-Obukhow长度,表示稳定性效应十分重要的高度。
10.1.2 微波散射计的测风原理 LiuLM Satellite Oceanic Remote Sensing 7
10.1.2 微波散射计的测风原理 LiuLM Satellite Oceanic Remote Sensing 8 后向散射系数与风向的关系 对于13.9Ghz的微波,在入射角40°时海洋单位面积的散射截面为相对于平均反射率的角度的函数。并且还可以看出后向散射系数对观测方向和雷达极化方式的灵敏响应。当观测方向与风向一致时,后向散射系数为最大值;当观测方向处于横风方向时,后向散射系数为最小值。顺风观测和逆风观测相比,后向散射系数略微减少。这种各向异性正是微波散射计获得风向数据的依据。
10.1 微波散射计测量海面风场 LiuLM Satellite Oceanic Remote Sensing 9 10.1.3 微波散射计海面风场反演 微波散射计风矢量反演过程可以分为三个主要步骤: 1.计算海面探测单元(WVC)的归一化后向散射系数,并获得不同视角上天线的对同一探测单元的观测; 2.利用风矢量与海面散射系数的关系,也即Ku或C波段的地球物理模式函数,进行风速、风向估计;由于散射的特性,这个处理会得到多个可能的风矢量解,即风矢量解的模糊性; 3.用模糊性消除算法来确定反演的风矢量场,以确保反演的风场与地球物理学上的风场相一致。
10.1.3 微波散射计海面风场反演 LiuLM Satellite Oceanic Remote Sensing 10 1. 风矢量反演模式 主要是通过对海面探测单元的不同观测与地球物理模式函数进行最小均方拟合提取风速和风向。模式函数可以表示为: 在许多实验的基础上,已获得了雷达后向散射截面б、海面参考高度处的风速U和入射波平面与风向之间夹角θ三者之间的经验关系,它们用下式表示:
10.1.3 微波散射计海面风场反演 LiuLM Satellite Oceanic Remote Sensing 11
10.1.3 微波散射计海面风场反演 LiuLM Satellite Oceanic Remote Sensing 12 式(10-1-12)中参数a和b可以从下面的不同值的Φ三组测量结果中推导出来:
10.1.3 微波散射计海面风场反演 LiuLM Satellite Oceanic Remote Sensing 13 固定风速时,б随方位角变化的示意图
10.1.3 微波散射计海面风场反演 LiuLM Satellite Oceanic Remote Sensing 14 上图给出了所测后向散射截面与θ之间的函数关系,其中θ为入射波平面与风向之间的夹角。从这图中可以看出: ①бvv总是大于 бhh; ②±Φ内对称; ③以180°为周期(bcos Φ项); ④顺风和逆风测量时有一点不对称(acos Φ项)。
10.1.3 微波散射计海面风场反演 LiuLM Satellite Oceanic Remote Sensing 15 2. 实际应用的海面风矢量反演模式 Ku波段的微波散射计模式 适用于Seasat/SASS资料的SASS-Ⅰ模式函数。利用SASS-Ⅰ模式函数,可以成功的对Seasat/SASS资料进行海面风速和风向的反演,风速精度达到2m/s、风向误差在20°范围内。 SASS-1模式函数表示为: Ku波段的Moore模型,具体形式为: SASS-Ⅱ基本表示形式如下:
10.1.3 微波散射计海面风场反演 LiuLM Satellite Oceanic Remote Sensing 16 C波段的微波散射计模式 主要工作在欧空局发射的ERS系列卫星上,与Ku波段的散射计相比,C波段散射计在海面与较长的短尺度重力波产生Bragg共振。 一、1985年提出了C波段的雷达后向散射系数与风矢量关系的第一个经验模型CMOD-1: 二、ESA根据卫星发射后所获得的散射计风场与现场资料比较,对算法进行改进和完善,又相继开发了CMOD-2、CMOD-3、CMOD-4、CMOD5-12等
10.2 卫星高度计测量海面风速 LiuLM Satellite Oceanic Remote Sensing 17 10.2.1 高度计测量海面风速的原理 卫星高度计是天底视主动式传感器,海面平静时回波信号最强。海面在风的作用下能够产生厘米尺度的波浪,从而引起海面粗糙度(海面均方斜率)的变化。海面起伏随风增大时,把信号反射回传感器的镜面面积越来越少,回波也就越来越弱。雷达高度计对于大于或等于其工作波长(2cm左右)的海面粗糙度变化有敏感反应。
10.2.1 高度计测量海面风速的原理 LiuLM Satellite Oceanic Remote Sensing 18 当海洋表面粗糙幅度超过发射脉冲宽度时,散射将首先发生在表面波峰处,然后是波谷处,从而导致回波引导沿变大,高度计接收到回波功率区域的面积增加。根据Walsh公式,对于粗糙海表面脉冲宽度有限足迹直径为: 由上式公式可知均方根波高hk越大,脉冲和海面状态共同决定的有限足迹就越大,因此海洋回波引导沿上升时间就越长。
10.2 卫星高度计测量海面风速 散射理论表明,雷达后向散射截面( )与海面均方斜率( )之间存在下列关系: 对卫星高度计而言,θ=0,则上式可变为较简单的形式: 研究表明, 与海面风速(U)近似满足线性关系: LiuLM Satellite Oceanic Remote Sensing 19 10.2.2 海面风速模型 即高度计后向散射截面和海面风速之间存在一种反比关系。
10.2.2 海面风速模型 LiuLM Satellite Oceanic Remote Sensing 20 1. 低海况时海面风速模型 根据镜面散射理论,高度计接收每个海面微小面元反射的电磁波。相对于这些海面微小面元,雷达波束的入射角θi =0。运用物理光学和电磁场的方程,Barrick推导得出 б° : 对于高度计,θ≈0。Wu和Valenzuela给出雷达后向散射系数为:
10.2.2 海面风速模型 为泡沫覆盖率。关于 的经验模式有以下两种: LiuLM Satellite Oceanic Remote Sensing 21 2. 高海况时海面风速模型 当海面风速高于临界风速时,海面会出现泡沫。风速越大,海面泡沫覆盖率越大。在高海况下海面反射率系数总的表达式为: 1) Jin Wu 模型,表达式为: 2) Monahan 模型,表达式为:
10.2 卫星高度计测量海面风速 LiuLM Satellite Oceanic Remote Sensing 22 10.2.3 高度计测量海面风速的经验模型 1. BR模式函数 Brown 等人认为,在垂直入射的情况下,光滑海面的后向散射能量与海面的平均斜方谱成反比,根据这个假定,他们得出了下列模式函数,其均方差为1.74m/s,适合范围是1~18 m/s。 2. SB模式函数 为了克服分段函数不光滑的缺点,Brown等建议在0~16 m/s范围对上式获得的风速进行多项式拟合,从而给出SB模式函数为:
10.2.3 高度计测量海面风速的经验模型 LiuLM Satellite Oceanic Remote Sensing 23 3. CM模式函数 Chelton 和McCabe提出了下列形式的模式函数: 4. GD模式函数 Goldhirsh 和Dobson对BR模式函数在雷达散射截面 空间进行5次多项式拟合,得出了GD模式函数:
10.2.3 高度计测量海面风速的经验模型 LiuLM Satellite Oceanic Remote Sensing 24 5. WC模式函数 Witter 和Chelton在假定风速概率分布函数的年际变化可忽略的前提下对CW 模式函数进行了进一步改进和推广:一是扩大了高风速部分的动态范围;二是使之适用于Geosat高度计。WC模式函数绘于10.2.2,其均方差为1.9m/s 。
10.2.3 高度计测量海面风速的经验模型 LiuLM Satellite Oceanic Remote Sensing 25 6.CC模式函数 Carter等提出了一个折线的模式函数: 7. WU模式函数 WU从镜面反射原理出发,利用实测海面斜率数据,并注意到Geosat高度计的工作波长(约2.16m)与对雷达后向散射有贡献的截止波长非常接近,导出了下列模式函数:
10.2.3 高度计测量海面风速的经验模型 LiuLM Satellite Oceanic Remote Sensing 26 9. GG模式 Glazman 和Greysukh的研究表明,雷达散射截面与海面的成熟度有关。他们通过有效波高引入波龄这一参数,并根据波龄不同分两种情况对海面风速进行反演: 1)在低波龄情况下: 2)在高波龄情况下:
10.2.3 高度计测量海面风速的经验模型 LiuLM Satellite Oceanic Remote Sensing 27 10. YG模式函数 Young通过对理论模式预报的和卫星高度计观测的热带气旋风速进行对比分析,给出了一个适用于高风速(20~40m/s)条件下的模式函数: 11. LB 模式函数 lefevre等建立了一个多项式形式的模式函数,包括后向散射截面和有效波高两个参数:
10.2.3 高度计测量海面风速的经验模型 LiuLM Satellite Oceanic Remote Sensing 28 12. FC模式函数 Freilich 和Challenor提出了一种确定高度计海面风速模式函数的新方法。该方法利用 б°随风速单值单调变化的特点,只需已知U和б°的独立分布,就可通过统计理论导出高度计海面风速模式函数,而无需采用通常靠同步数据确定回归系数的做法。FC模式函数的表达式为: 13. HT模式函数 Hwang 等从散射理论出发提出了一种考虑长波引起的海面倾斜对高度计后向散射截面影响的模式函数。他们假定倾斜海面的斜率呈正态分布,给出了如下形式的模式函数:
10.2.3 高度计测量海面风速的经验模型 LiuLM Satellite Oceanic Remote Sensing 29 适合于U>20m/s的高度计海面风速模式函数
10.3 微波辐射计测量海面风速 LiuLM Satellite Oceanic Remote Sensing 30 微波辐射计测量海面风速的原理主要是基于海面微波辐射率与海面粗糙度之间的高相关特征,而海面粗糙度又与风速有关。海面粗糙度增加,海面辐射率增加,极化特征变弱。 辐射计风速反演算法主要有两种: 一种是统计回归分析算法,主要是通过微波辐射亮温与现场风速之间的统计回归的关系; 一种是基于辐射传递的物理算法。
10.4 合成孔径雷达获取海面风场信息 LiuLM Satellite Oceanic Remote Sensing 31 合成孔径雷达(SAR)作为微波雷达也能获取海面风场的信息. SAR具有高空间分辨率(数米至数十米)测量海面风场的能力,可弥补散射计测风的不足,近年来受到人们的关注。2002年发射的欧洲Envisat卫星上装载的ASAR就兼任测量海面风场的重任,这也推动了SAR海面风场测量的研究。 SAR主动式遥感仪器,通过发射微波束、接收来自海面的后向散射获取海面信息。 SAR图像海面风速反演方法: 一、根据风速与雷达后向散射截面之间的关系计算海面风速; 二、根据SAR图像的方位向模糊原理,利用风与高波数截断波长之间的关系,计算海面风速。
10.4 合成孔径雷达获取海面风场信息 风速的反演由3个步骤组成: ①对子图像进行辐射定标,得到每个像元的后向散射截面 ②对 进行极化校正,得到 ③输入后向散射截面 ,入射角和风向值,由公式(10-4-1)反演子图像的风速 LiuLM Satellite Oceanic Remote Sensing 32