SAR 遥感图像的海面后向散射强度的变化很大程度上依赖于海面波谱中与SAR 工作波长相近的部分.海洋和大气现象对海面波谱的影响,在一定条件下可以成像于海面SAR 遥感图像.通常观测的大气现象包括:边界层对流､对流风暴以及由海岸地形或海气温差等引起的各种大量的中尺度环流.

千米尺度的海洋大气边界层对流是一种常见的大气现象,是在海洋表面虚位温高于底层大气虚位温的位置,由扰动动能的浮力效应形成.由于空气的垂直运动,对流涡旋中高动量的下降气体向下传输到近海面,导致近海面层的垂直风切变增加,以阵风的形式叠加在平均风场上,从而使下降气流区域的海面粗糙度增大,在SAR 遥感图像中形成了高后向散射强度区.反之,对流涡旋的上升区域对应于SAR 遥感图像中的低后向散射强度区.高分辨率的星载SAR 遥感图像足可以描述这种由浮力效应导致的水平变化的扰动尺度.

大气边界层对流的结果主要有两种形式,三维对流涡旋Cell和二维水平滚轴涡旋Roll.理论和观测结果都表明,三维对流涡旋Cell通常出现在低风速条件下,二维水平滚轴涡旋Roll在强风条件下更加普遍.Deardorff^[1]针对观测结果给出了理论解释,认为风速和海气温差对对流现象出现的形式起决定作用,在给定风速条件下,海气温差越大,就越可能形成三维对流涡旋Cell.

两种边界层对流现象都是基于高分辨率SAR气象学的重要研究目标.Sikora等^[2-4]论述了三维对流涡旋Cell在SAR 遥感图像上的斑点状后向散射特征,Mourad^[5]对二维水平滚轴涡旋Roll的线性特征SAR 遥感图像给出了解释.Young等^[6-7]利用具有三维对流涡旋Cell特征的SAR 遥感图像反演得到了MABL 高度及其浮力通量等其他扰动统计参数,Alpers 等^[8-9]利用具有二维水平滚轴涡旋Roll特征的SAR 图像反演了MABL 高度.由于气溶胶粒子主要来源于地球表面,故MABL 中富集气溶胶.Lidar可以探测到MABL 气溶胶后向散射信号和梯度变化,得到高分辨率的MABL 剖面结构特征,利用二次微分法处理探测的气溶胶后向散射信号,可以获取MABL 高度.

本文利用具有三维对流涡旋Cell和处于热不稳定状态的二维水平滚轴涡旋Roll对流特征现象的星载SAR 遥感图像,分别获取了三维对流涡旋Cell的特征长度和二维水平滚轴涡旋Roll的平均水平波长,根据边界层相似理论与对流理论得到了MABL 高度,首次对二者反演的MABL 高度进行对比,同时利用探空资料数据､Lidar 实验获取的MABL 高度进行验证分析.其结果显示二者符合较好,为后继以SAR､探空资料数据与Lidar为基础的集成技术探测MABL 特征提供了参考.

利用三维对流涡旋Cell的空间尺度和密度特征得到垂直于风向方向的平均一维水平能量谱,通过谱分析法^[10]确定小尺度峰值波长,即三维对流涡旋Cell的特征长度λCell,根据边界层相似理论提供的λCell与边界层高度Zi₁.5倍定量关系得到MABL高度^[11].

二维水平滚轴涡旋Roll通常用来求解风向,也利用二维水平滚轴涡旋Roll的空间尺度特性,得到垂直于风向方向的二维水平滚轴涡旋Roll平均水平波长λRoll,在热不稳定MABL 情况下,根据线性Rayleigh-Benars对流理论,利用2.8倍定量关系^[11],得到MABL高度.图 1为SAR 遥感图像反演MABL高度的流程.

图 1

Fig. 1

图 1 星载SAR 遥感图像反演MABL 高度流程 Fig. 1 Flowchart of the MABL depth retrieval from SAR image

基于星载SAR 遥感图像反演MABL 高度流程,首先对SAR 遥感图像进行Lee滤波以消除其固有的斑点噪声.由于海面平均风场的存在,SAR 遥感图像上呈现出具有线性排列的斑点状后向散射特征(三维对流涡旋Cell)或线性的风条纹纹理特征(二维水平滚轴涡旋Roll),因此可根据SAR 遥感图像低波数能量谱获取海面风向,但此方法获取的海面风向具有180°模糊,该海面风向的不确定性,可通过SAR 遥感图像本身信息或其他风场提供的信息加以消除^[12].

利用具有三维对流涡旋Cell特征的SAR 遥感图像反演MABL 高度时,需要对SAR 遥感图像进行辐射定标,得到雷达海面后向散射强度,以消除SAR 侧视所引起的能量值分布不均匀的现象.同时将遥感图像分辨率降为50m,消除高频信号对小尺度峰值波长信息的影响.

此外,SAR 遥感图像中不仅包含了大气对流引起的大气现象,也包含了许多大尺度､中尺度等海洋现象,可采用尺度分离的方式将三维对流涡旋Cell现象之外的特征滤除.三维对流涡旋Cell使得海表面风场模式发生变化,导致海面粗糙度发生相应变化,从而SAR 遥感图像很大程度上体现海面风场模式的变化特征,因此可利用谱相似性原理分析风速谱得到图像特征尺度分离的依据^[10].在体现大气现象的伽马中尺度和小尺度之间的谱隙对应的尺度大约是5Zi^[13].设600m 为实验海域MABL 的平均高度,则600m 与3000m 为尺度分离时的特征尺度,从而确保将大尺度和中尺度对一维水平能量谱的影响去除而不影响小尺度峰值波长信息.

最后采用双线形插值的方法将尺度分离后的SAR 遥感图像旋转到沿海面风向方向的自然坐标系,并在(K×log(S(K))､log(K))双对数坐标系下,得到垂直于风向方向的平均一维水平能量谱,确定三维对流涡旋Cell的特征长度,根据边界层相似理论得到MABL 高度.

对垂直于风向方向具有二维水平滚轴涡旋Roll特征的SAR 遥感图像要进行抽样平均,降低遥感图像分辨率,消除高频信息的影响,更好的体现条纹特征.最后滑动平均滤波后的平均一维水平灰度值确定平均水平波长.在热不稳定MABL 情况下,根据线性Rayleigh-Benars对流理论得到MABL 高度.

图 2a是青岛近海分辨率为25 m 的ENVISATASAR 遥感图像,成像时间为2005-01-1421∶41(BJT).从图中截取具有三维对流涡旋Cell特征的遥感子图像作为实例研究的数据(图 2b所示),可见截取的SAR 遥感子图像具有线性排列斑点状的海面后向散射特征.

图 2

Fig. 2

图 2 (a)ENVISAT ASAR 遥感图像;(b)具有三维对流涡旋的星载SAR 遥感子图像 Fig. 2 (a) ENVISAT ASAR image and (b) sub SAR image with three

分析截取的SAR 遥感子图像低波数能量谱,并基于分析图 2a中落山风所导致的近岸海域SAR 遥感图像呈亮色的信息确定风向为西北风(图 3 所示).图 4是进行特征尺度分离后的SAR 遥感子图像,可以看出,相比图 2b更好地体现斑点状海面后向散射特征.

图 3

Fig. 3

图 3 星载SAR 遥感子图像低波数能量图与海面风向 Fig. 3 Energnergy spectrum of sub SAR rnrnge and the wind drection

图 4

Fig. 4

图 4 特征尺度分离后的星载SAR 遥感子图像 Fig. 4 Diagnostic scale-s印amted sub SAR image

将特征尺度分离后的星载SAR 遥感子图像(图 4)进行快速傅里叶变换,得到垂直于风向方向的平均一维水平能量谱,判断三维对流涡旋Cell的特征长度为1560m(图 5a).根据边界层相似理论提供的λCell与边界层高度Zi₁.5倍定量关系,得到MABL高度为1040m.

图 5

Fig. 5

图 5 (a)垂直于风向的一维水平能量谱;(b)青岛海滨Lidar获取的MABL 高度;(c)青岛气象台探空资料数据位温剖面 Fig. 5 (a) Cross wind direction one dimensional energy spectrum； (b) MABL depth of Qingdao Seashore retrieved

SAR 遥感图像计算得到的MABL 高度结果与同步Lidar获取的MABL 高度(图 5b)及探空资料数据位温剖面(图 5c)进行验证分析.表 1列出了进行同步实验的相关地理位置和MABL 探测高度.

表 1(Table 1)

表 1 SAR、Lidar和气象探空资料获取的MABL高度的比较 Table 1 The comparisonbetween SAR，Lidar and weather data

表 1 SAR、Lidar和气象探空资料获取的MABL高度的比较 Table 1 The comparisonbetween SAR，Lidar and weather data

通过表 1可以看出,星载SAR 遥感图像反演的MABL高度比Lidar获取的MABL高度偏低,而气象探空资料估计的MABL 高度较高.λCell与Lidar获取的MABL 高度比值为1.35,接近1.5,与气象探空资料相比,SAR 遥感图像反演的MABL 高度误差达到14%,而Lidar获取的结果误差为5%.

图 6a是2006年春季SOLAS黄海实验同步ENVISAT ASAR 遥感图像,图像分辨率为25 m,成像时间是2006-04-2721∶31 (BJT).从该遥感图像左下方到左上方可以观测到二维水平滚轴涡旋Roll到三维对流涡旋Cell两种对流现象的转变,单纯从遥感图像的灰度值可以初步判断遥感图像左下方海域的风速值高于左上方海域的风速值.关于两种对流方式之间转变的定量关系,将会在以后的研究中展开.

图 6

Fig. 6

图 6 (a)ENVISAT ASAR 星载SAR 遥感图像;(b)具有三维对流涡旋Cell特征的星载SAR 遥感子图像;(c)具有二维水平滚轴涡旋Roll的星载SAR 遥感子图像 Fig. 6 (a) ENVISAT ASAR SAR image, (b) Sub SAR image with three dimensional convective cells and (c) Sub SAR image with two dimensional convective rolls

通过分析该研究海域SAR 遥感子图像低波数能量图,判断风向为西南风向.从该遥感图像中分别截取具有三维对流涡旋Cell特征的遥感子图像(图 6b)和垂直于风向方向上具有二维水平滚轴涡旋Roll特征的遥感子图像(图 6c)作为实例研究数据.基于三维对流涡旋Cell特征的SAR 遥感子图像,反演得到该海域MABL 高度为722m.

二维水平滚轴涡旋Roll在SAR遥感图像上体现为黑白相间的条纹,条纹间的间距一般在1~8km之间.基于星载SAR 遥感图像反演MABL 高度流程,将具有Roll特征的SAR 子图像平滑后的平均一维灰度值,利用其空间特性得到平均水平波长为1896m (图 7所示).

图 7

Fig. 7

图 7 滤波后的平均一维灰度值 Fig. 7 Filtered mean one dimensional intensity

根据SOLAS黄海实验现场实验数据表明(表 2所示),该时刻海洋表面虚位温高于底层大气的虚位温,即海气温差Tsea-air>0,则MABL 处于热不稳定状态.根据二维水平滚轴涡旋Roll水平波长与MABL高度在热不稳定条件下的线性Rayleigh-Benars对流理论,得到该海域MABL 为677m.

表 2(Table 2)

表 2 现场实验数据给出的海洋表面温与底层大气温度 Table 2 The potential temperature from the In-situ experiment data

表 2 现场实验数据给出的海洋表面温与底层大气温度 Table 2 The potential temperature from the In-situ experiment data

由于Lidar探测位置与SAR 遥感图像中心位置距离大约270km,所以Lidar获取的MABL 高度不能作为定量验证的标准.根据香港天文台2006-04-2708∶00提供的天气状况信息,Lidar数据可以作为SAR 遥感图像反演的MABL 高度参考值^[14].表 3给出进行同步实验的相关地理位置和MABL探测高度.

表 3(Table 3)

表 3 SAR和Udar所探测的MABL高度比较 Table 3 The comparison between SAR and Lidar

表 3 SAR和Udar所探测的MABL高度比较 Table 3 The comparison between SAR and Lidar

从表 3可以看出,利用具有三维对流涡旋Cell特征的SAR 遥感子图像反演的MABL 高度与利用具有二维水平滚轴涡旋Roll特征的SAR 遥感子图像反演MABL 高度结果均高于Lidar获取高度.将Lidar获取的MABL 高度设为SAR 成像位置MABL 的高度,得到λRoll与Lidar获取MABL 高度比值为3.01,介于实验观测比值2~6 之间^[9],与Stull^[11]提到的比值3最为接近.λCell与Lidar获取的MABL高度比值为1.7,在实验经验比值范围1~5之间^[2],接近1.5.

本文通过具有三维对流涡旋Cell和二维水平滚轴涡旋Roll特征的星载SAR 遥感图像反演了MABL 高度,首次对二者反演的MABL 结果进行比较,结果吻合较好.同时首次提出利用Lidar获取的MABL 高度验证星载SAR 遥感图像反演的MABL 高度,为SAR､Lidar与探空资料数据集成技术探测MABL 特征进行了首次尝试.

通过本文表 1与表 3可以看出,SAR 遥感图像反演的MABL 高度与Lidar获取的MABL 高度结果虽然基本一致,但有一定的出入.这可能由以下因素造成:(1)λCell/Roll与MABL 高度的定量比值要求在实测的风向条件下,而本文的海面风向是通过对SAR 遥感子图像处理得到,与实际风向会有一定的出入,影响MABL 高度反演结果.(2)对SAR 遥感图像进行快速傅里叶变换得到垂直于海面风向方向的平均一维水平能量谱时,由于风向与实际风向的出入而引入沿风向方向的剪切力对最小峰值信息造成的失真,及能量谱离散化过程中由于分辨率限制导致最小峰值的偏差.(3)尺度分离过程中特征尺度的选择,及其根据先验知识选择的图像旋转到自然坐标系下存在斑点状后向散射特征的不均匀性问题,都会引入一维水平能量谱多尺度的峰值信息.而通过平均得到的一维水平能量谱和一维灰度值,也会降低低波数能量谱的可信度和空间分布信息.最后,Lidar测量气溶胶后向散射信号及其通过探空资料数据得到MABL 高度过程中,存在的系统误差与二者测量时刻的时间延迟及空间间距也会对验证分析结果产生一定的影响.

本研究通过实例研究与验证分析展示了以高分辨率､大面积观测为特点的星载SAR 遥感图像探测MABL 的广阔前景,同时表明可利用探空资料､Lidar的高精度剖面测量MABL 结果,实时､同步地校正星载SAR 遥感图像大面积观测MABL 的定量反演结果,为后继以SAR､探空资料与Lidar为基础的集成技术探测MABL 特征提供了参考.