《气溶胶光学厚度》由会员分享,可在线阅读,更多相关《气溶胶光学厚度(11页珍藏版)》请在金锄头文库上搜索。
1、第2 章 气溶胶光学厚度反演的原理和方法气溶胶光学厚度(Aerosol Optical Depth)简称AOD,定义为介质的消光系 数在垂直方向上的积分,描述的是气溶胶对光的消减作用7。它是气溶胶最重 要的参数之一,表征大气浑浊程度的关键物理量,也是确定气溶胶气候效应的 重要因素。通常高的AOD值预示着气溶胶纵向积累的增长,因此导致了大 气能见度的降低。现阶段对于AOD的监测主要有地基遥感和卫星遥感两种方 法。其中地基遥感又有多种形式:多波段光度计遥感、全波段太阳直接辐射遥 感、激光雷达遥感等。其中多波段光度计遥感是目前地基遥感研究中采用的最 广泛的方法。美国NASA和法国LOA-PHOTON
2、S联合建立的全球地基气溶胶 遥感观测网AERONET所使用的就是多波段太阳光度计(Sun/Sky Photomerers),在全球共布设1217个站点长期观测全球气溶胶的光学特性,积 累了大量的AOD数据,并用作检测气溶胶光学厚度反演精度的标准。而近年 来卫星遥感技术的快速发展,多种传感器被用来研究气溶胶特性,加上经济发 展带来的大气污染问题使得利用卫星遥感资料反演AOD成为热门课题。2.1 气溶胶光学厚度反演的基本原理大气光学厚度是指沿辐射传输路径单位截面上气体吸收和粒子散射产生的 总消弱,是无纲量值。在可见光和近红外波段,它可以由下列公式计算得出:T(九)=工(九)+ 工(九)+ x (九
3、)+ (九)+ T (九)(2-1)m12卩a其中T(九)表示大气总的光学厚度,t (九)表示整层大气的分子散射光学厚m度,t (九)表示氧气的吸收光学厚度,T (九)表示臭氧的吸收光学厚度,T (九)12卩表示水汽的吸收光学厚度,T (九)表示气溶胶光学厚度21; 22。a卫星遥感反演大气气溶胶是利用卫星传感器探测到的大气顶部的反射率, 也称为表观反射率,可以表示为23:p* = :n:L / 卩 F(2-2)ss其中,L表示卫星传感器探测到的辐射值,F表示大气上界太阳辐射通s量,卩表示太阳天顶角e的余弦值。p*与地表二项反射率之间的关系可以表达 ss为:p (e , e , e) F (e
4、 )T (e )p*(e , e , e) = p (e , e , e) +v sa v s(2-3)vs 一 d sv1 一 s p*其中, e 表示传感器天顶角, e 表示太阳天顶角, e 表示太阳方位角和卫 vs星方位角确定的相对方位角;p (e , e , e)表示由大气分子和气溶胶散射造成a v s的路径辐射,它与地表状况无关; F (e )表示地表反射率归一化为零时总的向 ds下辐射通量,也可以称为总的向下透过率,由于气溶胶粒子对太阳光的吸收和散射作用,它的值小于1.0; T (e )是向上进入卫星传感器视场方向的总透过 v率, S 是大气后向散射比。在单次散射近似中,路径辐射p
5、 (e ,e ,e)与气溶胶光学厚度t和单次散 a v sa射反射率之间的关系如下24:0w t p (e , e , e)p (e , e , e )=p (e , e , e)+v s(2-4)a v s m v s4 卩卩vS其中p (e , e , e)是分子散射造成的路径辐射,它取决与大气模式,卩表 m v sv示传感器天顶角的余弦值,卩 表示太阳天顶角的余弦值。在式(2-3)中,Sf (e)、t (e)和 s 取决于,t 和 p (e , e , e)。d sv0 a a v s假设地表是均匀朗伯表面,大气垂直均匀变化,将式( 2-4)代入( 2-3)得:w t p (0 , e
6、, e) p*(e, 0, e)= p (0,0, e)+ a a s(2-5)v sm V s4 卩 paSp (0 ,0 ,e)F (0 ,w ,t ,P)T(0 ,w ,t ,P)+ v sdS O a. aS Oa. a1 s (w , t ,P)p*Oaa上式中 p (0 ,0 ,e )为假设的朗伯体特性的地表反射率,卫星传感器接收 vs到的表观发射率 p*(0 ,0 ,e )既是地表反射率 p (0 ,0 ,e )的函数,又是气溶v sv s胶光学厚度t的函数I。表观反射率p*(0,0 ,e)以及太阳和传感器的几何参 av s数(0 ,0 ,e )可以从卫星遥感资料中获取,假如可以
7、得知地表反射率vsp* , 0, e),并用气溶胶类型和大气模式来确定和p , 0, e)的相关参v sO a v s数,理论上就可以计算得出地面上空的气溶胶光学厚度t。反之,若已知地面a上空气溶胶光学厚度t、气溶胶类型以及大气模式,也可以反演出地表反射率 ap*(0 , 0 , e )26。vs由式(2-3)可以得出反演气溶胶的最优条件是地表反射率低且光谱波段波 波长较短。在地表反射率角度的情况下,气溶胶散射引起的路径辐射 p (0,0,e)(与地表状况无关)对表观反射率p* (0,0 ,e)起主要作用,此a v sv s时反演气溶胶光学厚度误差较小;而在地表反射率较大的情况下,地表贡献项
8、p (0v,0s,e)Fd(0s)T (0v)对表观反射率p*(0,0 ,e)影响较大,此时反演精度1 - s p*v s较低27。2.2 气溶胶光学厚度反演的主要方法自 20 世纪 70 年代中期开始,利用卫星数据反演气溶胶光学厚度的研究已 经有40 年的历史,反演的方法有单通道算法、多通道算法、暗像元法、结构函 数法、深蓝算法、多星协同反演法、海陆对比法、多角度偏振法、热辐射对比 等28。目前有代表性的常用气溶胶光学厚度反演算法有两种:一种是通过路径 辐射项求取气溶胶光学厚度的暗像元法,另一种是通过透过率求取气溶胶光学 厚度的对比法。2.2.1 暗像元法2P*2.4(2-6)暗像元算法英文
9、全称为Dense Dark Vegetation,简称DDV,它是通过路径 辐射项来计算气溶胶光学厚度。由于地表物体的复杂多样性造成反射率变化范 围很大,很难从辐射值中分理处辐射项,如果想通过辐射项来获取气溶胶信 息,就必须使地表辐射值较小且能确定其精确值,这样就能够最大限度的消除 地表反射率的不确定性带来的影响。在卫星影像中,大量浓密植被区由于在可 见光波段反射率极低(约为0.010.02),它们被称作暗像元。1988年Kaufman等 利用大多数地物在红(0.600.68 pm )蓝(0.400.48 pm)波段反射率低的特性,根 据归一化植被指数(NDVI)或近红外通道(2.1pm)的表
10、观反射率进行暗像元 的识别,并假定已知这些暗像元在红蓝通道的地表反射率,依据一定的关系反 演气溶胶光学厚度29。通过大量的卫星影像资料,考虑到多种地表覆盖物,拟 合得到红(0.66pm)蓝(0.47pm)和中红外通道(2.1pm)地表反射率的关系 29:red blue对于2.1pm通道卫星观测表观反射率几乎不受气溶胶影响,其值接近地表 反射率,因此可以用2.1pm通道的表观反射率代替地表反射率,根据式(2-6) 计算出红蓝通道的地表反射率。再假定合适的气溶胶模型,就可以计算出气溶 胶光学厚度。然而对于干旱、半干旱以及冬季城市等高发射率地区,用暗像元法来反演 气溶胶光学厚度还存在不少困难,在以
11、上地区由于地表的非均一性使确定地表 反射率的精确值十分不易30。当地表反射率升高时,气溶胶的指示作用降低 31。一般情况下,当地表反射率较低时,传感器接收到的辐射值随着气溶胶的 增多而迅速变大,暗像元法就是利用浓密植被在红、蓝波段的地表反射了和气 溶胶的这种关系来反演光学厚度;随着地表反射率的不断增大,接收到的辐射 值随气溶胶的增多而增大速度减缓,当地表反射率增大到某一程度时,辐射值 将不随气溶胶的增多而增大,甚至会出现降低的趋势。为了使暗像元方法的应 用更为广泛,2002年Kaufman等通过大量的数据验证,对以前的方法进行了扩 展:对于星下点,暗像元法适用于中红外通道表观反射率小于0.4的
12、区域;对 于非星下点,需要考虑太阳和卫星的几何参数,适用范围可以扩展到中红外通道表观反射率小于0.125(1 + 1,.必)的地区(其中卩为卫星天顶角的余弦值, 卩0为太阳天顶角的余弦值),红蓝通道的地表反射率依然遵循式(2-6)的函数 关系。2.2.2扩展的暗像元法(V5.2算法)扩展后的暗像元法是由 Levy 和 Remer 等人提出的称为 V5.2 的气溶胶反演 算法32; 33。与原来的算法相比,它主要有两方面的改进:首先提出 VIS/SWIR (可见光和中红外通道的反射率比值)不再是一个常量,而是一个函数关系; 其次考虑了 NDVI (植被指数)对于地表反射率的影响,红蓝通道和中红外
13、通 道的地表反射率关系如式(2-7)所示34:ps0.66ps0.47=f(ps )2.1=g(ps )0.66(2-7)该方法中的VIS/SWIR的变化跟不但跟地表植被指数有关,根据Remer和 Gatebe等的研究表明它还与太阳天顶角、传感器方位角、散射角有关,尤以散 射角的影响最大35; 36。它可以表示为:0 = cos ( cos A cos 0 + sin 0 sin 0 cos d)00(2-8)其中0表示观测天顶角,00表示太阳天顶角,0表示太阳方位角和卫星方 位角的相对方位角。根据式(2-8)可将式(2-7)扩展为:ps = f(ps ) = psslope + y int0
14、. 66 2. 1 2. 12 0. 66 2. 1 0. 66 2. 1 ps = g(ps ) = ps slope+ y int0.470.660.660.47 0.660.47 0.66(2-9)其中slope = slopeNDVIsWiR + 0. 0020 - 0. 270. 662. 1 0. 66 2. 1yint = 0. 000250 + 0. 0030. 66 2. 1slope = 0. 490. 47 0. 66yint = 0. 0050. 47 0. 66(2-10)NDVI SWIR可表示为:(2-11)其中,pm和pm分别表示MODIS第5波段和第7波段接收
15、到的辐射值,1.2 2.1当NDVI大于0.6时,表示该地区的植被状况较好;当NDVI小于0.2时,表SWIRSWIR示该地区植被稀疏37。当 NDVI 0.75 时,slopeNDVIsWIR 二 0. 58 ;SWIR0. 66 2. 1当 0.25 NDVI 0.75 时,slopeNDVIR 二 0.48 + 0.2 * (NDVI 0. 25)SWIR0.66 2. 1SWIR其中slope和yint为统计方程的经验系数,slope为偏移量,yint为截距。将(2-10)代入(2-9)可得出:a.植被指数小于0.25的情况:b.c.005植被指数大于 0.75 的情况:=ps x (0. 58 + 0. 0020 0. 27) + 0. 000250 + 0. 033x 0.49 + 0. 005植被指数介于 0.2