基于MODIS数据的我国东南沿海某区域海面温度反演研究.pdfVIP

基于MODIS数据的我国东南沿海某区域海面温度反演研究 宋佳李良超 (西安电子科技大学理学院，西安710071) 海表面温度(SST)是海洋的主要物理参数之一，在大气与海洋间热量、动量及水汽 交换中扮演着重要的角色，它是决定海气相互作用及气候变动的重要因素。中国是海洋 大国，有着丰富的海洋资源，获取连续的海洋参数观测数据是研究海洋的必要前提，对 于我国合理有效地开发海洋资源也有着深远的影响【l】。近20年来，热红外遥感技术的飞 速发展为快速地获取区域海面温度空间差异信息提供了新的途径【2l。本文主要研究基于 中分辨率成像光谱仪MODIS数据的海面温度反演方法并反演出我国东南沿海某区域的 海面温度。 热红外遥感是获取海面热量空间分布状况的一种非常有效的手段。本文采用劈窗算 法来反演海面温度，它的主要思想是利用两个通道对水汽吸收和发射率的差异来分别建 立方程，通过解方程组进行海面温度的反演。劈窗算法是针对AVHRR的第4和5热红 因此MODIS热红外数据的3l和32波段最适合用来进行海面温度的反演瞳3。由于大气 和海面的双重影响，热红外遥感观测到的不是海面的真实温度，而是混合有大气影响和 海面热特征影响的亮度温度。对于前者，主要有大气的吸收问题，大气的吸收作用直接 决定着海面热辐射在大气中的透射程度即大气透过率。对于后者，海面发射率是一个非 常重要的因素。因此，大气透过率和海面发射率是从热红外遥感数据中反演海面温度的 基本参数。 热辐射传输方程是热红外遥感和海面温度反演的基础，考虑到大气和海面影响因素 的复杂性和不确定性，热辐射过程可表达为【2】： (1) E(Z)=s，(8)f，(9)E(C)+扛一f，(o)11一s，(p)】r，(D)E(艺J)+1If，(臼)k(￡) 式中f，(a)代表波段f的大气透过率，￡。CO)代表波段f的在观察方向0的海面发射率， 互代表海面温度，Z代表波段f的星上亮度温度，E(C)代表海面辐射强度，瓦为大气向 ． L。2 上平均温度，乃上为大气向下的平均作用温度。Planck函数日(丁)2万南描述了 黑体的光谱发射特性，是热辐射传输方程的核心组成部分，对Planck函数进行线性简化， 是从热辐射传输方程中推导海面温度反 演方法的前提。毛克彪【2】通过分析表明，针对特定的热红外波段和温度区间，可将Planck 热辐射与温度在273～322K区间内的变化关系进行计算与分析，得到热辐射强度随温度 的变化接近于如下的线性关系【2】： (2) 1849T32—26．50036 BI(丁)=O．13787Tal一31．65677；B2(丁)=O．1 scaleB、radiance 象素点的计数值，radiance offsetB的值可在相应波段数据集的属性 域中获得【3】。根据式(2)可以分别得到MODIS第31和32波段的亮度温度。 海面的辐射能经过大气传输之后才能抵达空中的遥感器。由于大气对辐射能的削弱 作用，往往只有部分辐射能最后抵达遥感器。这种削弱作用一般用大气透过率f来描述。 在海面温度反演过程中，水汽是估计大气透过率的主要考虑因素。在MODIS的36个波 段，而2波段和5波段是大气窗口波段【21。这样就可以利用比值法来求大气水汽含量， 计算公式如下【3】【4】： 中获得【31。而MODIS的31和32波段的大气透过率与大气水汽含量之间的关系可以用 LOWTRAN、MODTRAN等大气模拟软件来模拟估计，它们之间满足线性关系【2】： (3) r3l=0．97366-0．05468w；f32=0．9621-0．08991w 因此，利用式(3)可以得到31和32波段的大气透过率。 海面发射率￡定义为在相同温度和波长下观测地球表面的光谱辐射