气溶胶光学厚度

1、第2章气溶胶光学厚度反演的原理和方法气溶胶光学厚度(Aerosol Optical Depth)简称AOD ,定义为介质的消光系 数在垂直方向上的积分，描述的是气溶胶对光的消减作用7。它是气溶胶最重要的参数之一，表征大气浑浊程度的关键物理量，也是确定气溶胶气候效应的重要 因素。通常高的AOD值预示着气溶胶纵向积累的增长，因此导致了大气能见 度的降低。现阶段对丁 AOD的监测主要有地基遥感和卫星遥感两种方法。其中 地基遥感乂有多种形式：多波段光度计遥感、全波段太阳直接辐射遥感、激光$ 达遥感等。其中多波段光度计遥感是目前地基遥感研究中采用的最广泛的方法。美国NASA和法国LOA-PHOTONS联

2、合建立的全球地基气溶胶遥感观测网 AERONET所使用的就是多波段太阳光度计(Sun/Sky PhotomererS ,在全球共 布设1217个站点长期观测全球气溶胶的光学特性，积累了大量的AOD数据，并用作检测气溶胶光学厚度反演精度的标准。而近年来卫星遥感技术的快速发展， 多种传感器被用来研究气溶胶特性，加上经济发展带来的大气污染问题使得利用 卫星遥感资料反演AOD成为热门课题。2.1气溶胶光学厚度反演的基本原理大气光学厚度是指沿辐射传输路径单位截面上气体吸收和粒子散射产生的 总消弱，是无纲量值。在可见光和近红外波段，它可以由下列公式计算得出：()=m( ).1( ).2( ) Q) a(

3、)(2-1)其中(九)表示大气总的光学厚度，Tm(九)表示整层大气的分子散射光学厚度，崩A)表示氧气的吸收光学厚度，旧表示臭氧的吸收光学厚度，平(幻表示水汽的吸收光学厚度， 讯九)表示气溶胶光学厚度21; 22。卫星遥感反演大气气溶胶是利用卫星传感器探测到的大气顶部的反射率，也称为表观反射率，可以表示为23：*=L / 氏 Fs(2-2)其中，L表示卫星传感器探测到的辐射值，Fs表示大气上界太阳辐射通量， . 一 改表小太阳天顶角隹的余弦值。P与地表二项反射率之间的关系可以表达为：.、：( X. us , ) Fd(志)T ( ».)i，( r.，rs , ，)= :，a( r .，

4、rs ,)1 f s ( 0 , a , Pa)(2-3)1 - S其中，e*表示传感器天顶角，缶表示太阳天顶角，巾表示太阳方位角和卫 星方位角确定的相对方位角；pa v, es,。)表示由大气分子和气溶胶散射造成的 路径辐射，它与地表状况无关；Fd(岳)表示地表反射率归一化为零时总的向下辐 射通量，也可以称为总的向下透过率，由丁气溶胶粒子对太阳光的吸收和散射作 用，它的值小丁 1.0； T(e°是向上进入卫星传感器视场方向的总透过率，S是大气后向散射比。在单次散射近似中，路径辐射Pa(0v，氏，巾)与气溶胶光学厚度p和单次散射 反射率妙之间的关系如下啊：0 .a Fa( *. ,

5、K ,)，a(乂，吐，)=，m(矿，Us, )： (2-4)4 S其中电(隹，e)是分子散射造成的路径辐射，它取决与大气模式，表示传感器天顶角的余弦值，卜 表示太阳天顶角的余弦值。在式(2-3)中，Fd(&)、 ， St3Q和S取决丁灿,鬲和Pa3v，es,。)。假设地表是均匀朗伯表面，大气垂直均匀变化，将式(2-4)代入(2-3)得：-*0 .a Pa( U、.，Us ,，)"昊，.)=m( 3，H , .)4 土.*(»., Us,) Fd( US,，0 , a ,Pa) T(US,，0 ,.a, Pa)(2-5)上式中p（%, 6s, e）为假设的朗伯体特性的

6、地表反射率，卫星传感器接收到* 一»、 , 、 、的表观发射率p （8*, a ,。）既是地表反射率p （e*，%, e）的函数，乂是气浴胶光 学厚度强的函数25。表观反射率1（匕,%, e）以及太阳和传感器的几何参数 一 一 ，,一 *（土，© ）可以从卫星遥感资料中获取，假如可以得知地表反射率P （日¥, % , 0 ）, 并用气溶胶类型和大气模式来确定 必和Pa（ 6*, % ,巾）的相关参数，理论上就可 以计算得出地面上空的气溶胶光学厚度帝。反之，若已知地面上空气溶胶光学厚度伞、气溶胶类型以及大气模式，也可以反演出地表反射率pv, 0s, e ）26。由式

7、（2-3）可以得出反演气溶胶的最优条件是地表反射率低且光谱波段波 波长较短。在地表反射率角度的情况下，气溶胶散射引起的路径辐射» 、 * 一 一 、 . 一 、pa（ e,氏，e）（与地表状况无关）对表观反射率p （6*,氏，e）起王要作用,此时 反演气溶胶光学厚度误差较小；而在地表反射率较大的情况下，地表贡献项P（膈B , " Fd（心T ® 对表观反射率p*3v,%,。）影响较大，此时反演精度1 - S较低27。2.2气溶胶光学厚度反演的主要方法自20世纪70年代中期开始，利用卫星数据反演气溶胶光学厚度的研究已经 有40年的历史，反演的方法有单通道算法、多通道

8、算法、暗像元法、结构函数 法、深蓝算法、多星协同反演法、海陆对比法、多角度偏振法、热辐射对比等28。 目前有代表性的常用气溶胶光学厚度反演算法有两种：一种是通过路径辐射项求取气溶胶光学厚度的暗像元法，另一种是通过透过率求取气溶胶光学厚度的对比 法。2.2.1暗像元法暗像元算法英文全称为 Dense Dark Vegetation,简称DDV ,它是通过路径辐 射项来计算气溶胶光学厚度。由丁地表物体的复杂多样性造成反射率变化范围很 大,彳艮难从辐射值中分理处辐射项，如果想通过辐射项来获取气溶胶信息，就必须使地表辐射值较小且能确定其精确值，这样就能够最大限度的消除地表反射率 的不确定性带来的影响。

9、在卫星影像中，大量浓密植被区由丁在可见光波段反射 率极低（约为0.010.02）,它们被称作暗像元。1988年Kaufman等利用大多数地 物在红（0.600.68两蓝（0.400.48»W反射率低的特性，根据归一化植被指 数（NDVI ）或近红外通道（2.1 H小的表观反射率进行暗像元的识别，并假定已 知这些暗像元在红蓝通道的地表反射率，依据一定的关系反演气溶胶光学厚度 29。通过大量的卫星影像资料，考虑到多种地表覆盖物，拟合得到红（0.66 V小 蓝（0.47两 和中红外通道（2.1两地表反射率的关系29:匕（2-6）2. 1:blue =4对丁 2.1 H加道卫星观测表观反射率

10、几乎不受气溶胶影响， 其值接近地表反 射率，因此可以用2.1 pni通道的表观反射率代替地表反射率，根据式（2-6）计 算出红蓝通道的地表反射率。再假定合适的气溶胶模型，就可以计算出气溶胶光 学厚度。然而对丁干旱、半干旱以及冬季城市等高发射率地区，用暗像元法来反演气 溶胶光学厚度还存在不少困难，在以上地区由丁地表的非均一性使确定地表反射 率的精确值十分不易30。当地表反射率升高时，气溶胶的指示作用降低31。一 般情况下，当地表反射率较低时，传感器接收到的辐射值随着气溶胶的增多而迅 速变大，暗像元法就是利用浓密植被在红、蓝波段的地表反射了和气溶胶的这种关系来反演光学厚度；随着地表反射率的不断增大

11、，接收到的辐射值随气溶胶的 增多而增大速度减缓，当地表反射率增大到某一程度时，辐射值将不随气溶胶的 增多而增大，甚至会出现降低的趋势。为了使暗像元方法的应用更为广泛，2002 年Kaufman等通过大量的数据验证，对以前的方法进行了扩展：对丁星下点， 暗像元法适用丁中红外通道表观反射率小丁0.4的区域；对丁非星下点，需要考虑太阳和卫星的几何参数，适用范围可以扩展到中红外通道表观反射率小丁 0.125（仰+ 1人稣）的地区（其中卜为卫星天顶角的余弦值，匕为太阳天顶角的 余弦值），红蓝通道的地表反射率依然遵循式（2-6）的函数关系。2.2.2扩展的暗像元法（V5.2算法）扩展后的暗像元法是由Lev

12、y和Remer等人提出的称为V5.2的气溶胶反演算 法32; 33。与原来的算法相比，它主要有两方面的改进：首先提出VIS/SWIR （可见光和中红外通道的反射率比值） 不再是一个常量，而是一个函数关系；其次考 虑了 NDVI （植被指数）对丁地表反射率的影响，红蓝通道和中红外通道的地表 反射率关系如式（2-7）所示34:(2-7)0.66 = f( ; 2. 1 )0.47= g(0.66)该方法中的VIS/SWIR的变化跟不但跟地表植被指数有关，根据Remer和Gatebe等的研究表明它还与太阳天顶角、传感器方位角、散射角有关，尤以散射角的影响最大35; 36。它可以表示为：-cos 1(

13、 cos u°cos sin ±sin cos )(2-8)其中8表示观测天顶角,e°表示太阳天顶角，巾表示太阳方位角和卫星方位角的相对方位角根据式(2-8)可将式(2-7)扩展为:s:0. 66s:0. 47=f( ?2. 1)=a0.66)-sP2.12 slope 0.66.2.1 + y int 0.662.1、s.P0.66 slope 0.47,0.66 十 y int 0.47266(2-9)(2-10)(2-11)其中slope 0.66.2.1 = slope 岩6WR1 + 0. 00为-0.27yint o&.w = 0. 0002务

14、 + 0.003slope 0. 47 0.66 = 0. 49y i n(t47 0.66 = 0. 0 0 5NDVSwir可表示为:m 瑚、/ .m瑚、NDVSwir = (/. 2 - ? 2.1) ( :' 1.2"2.)其中，P：2和2°1分别表示MODIS第5波段和第7波段接收到的辐射值，当NDVSwir大丁 0.6时，表示该地区的植被状况较好；当NDVIswir小丁 0.2时，表示该地区植被稀疏37。当 NDVIswir < 0- 25 时，slope ；.籍雷=0. 48 ；当 NDVSwir0-75 时，slope 携愣=0. 58 ；当

15、0. 25 M NDVIswir M 0.75 时,slope N.譬於=0. 48 + 0. 2 * ( nDVSwir 0. 25)其中slope和yint为统计方程的经验系数，slope为偏移量，yint为截距。将(2-10)代入(2-9)可得出：a. 植被指数小丁 0.25的情况：ss(2-12)p0.66 = p2.1 x(0. 48 + 0. 002g 0. 27) + 0. 0002钓 + 0. 033 ssP0.47 = P0.66 X °. 49 + 0. 005b. 植被指数大丁 0.75的情况：ss(2-13)P0.66 = P2.1 x(0. 58 + 0.

16、002。- 0. 27) + 0. 0002的 + 0. 033 ss_-_ _ _：0.47 = ：".660.490. °。5c. 植被指数介丁 0.25和0.75之间的情况：ss(2-14)P0.66 = P2.1 X 0. 48 + 0. z NDVIswir - 0. 25) + 0. 002。- 0. 270. 00025臼0. 033ss_-_ _"47 = ；0.660. 490. 0052.2.3结构函数法此时应用暗对丁旱季和中高纬度地区的冬季，植被覆盖少，地物反射率高， 目标法会造成一定的误差，因此，针对陆地亮地表的情况，发展了结构函数法来 反

17、演气溶胶光学厚度。结构函数法也称为对比法，是早期用来研究陆地污染气溶胶采用的卫星遥感 算法。它跟暗像元法一样采用的是红、 蓝通道数据。对比法需要假定在同一地区 同一时间段内地表反射率不变38。然后找出污染极小或无污染的一天作为活洁日”，以它的气溶胶数据作为参考反演 污染日”的大气气溶胶光学厚度。Tanr £和 Holben在采用对比法进行气溶胶反演的时候，假定在同一卫星观测时刻，对丁 两个相邻的像元大气是均一的，并在反演的过程中引入了结构函数的概念，使得地表反射率的确定变得容易39; 40。该算法主要利用的是表观反射率的地表贡献 项来反演气溶胶光学厚度。在暗像元方法不适用的地区，结构

18、函数法为气溶胶光 学厚度的反演提供了一条新的途径。在反演过程中，需要通过一段时间内卫星观测数据的分析， 选定其中无污染 （气溶胶光学厚度极小）的一天作为活洁日”，然后通过地面观测或者其他途径 确定这一天的气溶胶光学厚度值，作为背景的气溶胶信息。假定地表目标无变化， 通过透射函数的变化就能获取其它 污染日”的气溶胶光学厚度。由丁总透射函数 对丁相函数细节不敏感，结构函数法对散射相函数有较大的独立性，获取气溶胶光学厚度主要根据单次散射反照率和不对称因子。由式（2-3）得出相邻两点的表观反射率差值为：(2-15)* R, j) = Ri, j)听 S)( 0 V)地表结构函数为40:2 /M(d)=

19、n m_d2n(md)H(n，j E(2-16)改进后的地表结构函数为41:21n mH222M (Cl) 一 3(n - d)(m _ d) i J( i，j " ： i，j d) ( %j " d，j)(;?u _;?i，j d)(2-17)在卫星高度上的结构函数为:(2-18)M2(d) = M2(d)T2(us)T2(血)假定 消洁日”仕和 行染日” 2地表特征保持不变,则 M（d,t） = M（d, t2）,观测结构函数关系根据式(2-12)得：*(2-19)M (d, t 1)mt)，(t1)， a(t1)*M (d, t 2)EusO, v(t2), a(t2

20、)由上式可得在已知 消洁日”的气溶胶光学厚度的情况下， 污染日”的气溶胶光学 厚度就可以由卫星观测值获取。2.2.4双星协同反演算法近年来随着卫星遥感技术的不断发展，对地观测卫星携带的传感器种类不断 增多，这对丁不同平台之间的结合研究起到了推动作用。 为了解决气溶胶光学厚 度反演中未知参数过多导致气溶胶光学厚度信息难以提取的问题,唐家奎等提出了基丁Terra和Aqua双星的MODIS遥感数据协同反演算法17。该算法具有以下几种优点：（1）不需要事先假设气溶胶的类型等参数。（2）不需要假设或利用其它统计资 料预先估计地表反射率，双星协同算法是将真实的表观反射率作为反演气溶胶光 学厚度的变量之一与

21、其同时反演。（3）不收地表反射类型的影响，可应用高反 射率的陆地气溶胶反演。（4）可以同时获取两次过境时间的气溶胶信息，有利 丁分析气溶胶的动态变化规律17。基本原理如下：在整个大气层范围内对微元体积大气的辐射性质进行积分称为地球一大气 系统的辐射传输模型。根据以前学者的研究分析42,地球一大气系统之间的辐 射传输方程可以表达为：cos I 二:dz 4 二I (z, ) (z,r ,r)d , -(k；)I (z, r)(2-20)其中H表小传感器天顶角，P表小空气密度，°表小散射系数，z表小周度,?（z, r；r）表示散射相函数，G表示固体角,r表示方位（天顶角，方位角）,I z

22、,r）表示在z高度r方向的辐射强度。方程（2-20）中不能够得到 伙乙r）的 解析解，只能够获得方程的近似解，Xue等通过用多个不同变量将方程简化后得 到地表反射率L和表观反射率L之间的关系如下43：(a -b)八sec(2-21)： (Lb - a) a(1 - L)e 0L =,(a-b)【 vsec -i(Lb - a)妇- L)e 0其中，表示后向散射系数（通常取值为 0.1） , a = sec 6b=2,击表示表示气溶胶光学厚度。在双星协同反演模型中，仅考虑了空气分子和气溶胶粒子散 射两种影响因素，将大气气溶胶光学厚度分为大气分子的瑞利散射和大气气溶胶 粒子的散射两部分：瘟=盘刘+

23、 /3）（2-22）根据Linke等的研究大气分子的瑞丽散射近似解为44:项（°°） = 0. 0087-'09（2-23）利用？ngstr?m浊度公式可以求出气溶胶粒子散射光学厚度为：服）=芝（2-24）将式（2-24）和（2-23）以及（2-22）代入（2-21）中，可以得出地表反射率L以及?ngstr? m浊度系数E和波长指数a之间的关系方程。首先假设在两次卫星过境时间间隔内，地表反射率不变，且气溶胶类型保持 不变，这样波长指数"不变，只有气溶胶粒子浓度变化，即浊度系数E是变化的。然后利用0.47,0.55以及0.66Pm三个通道实现双星协同反演算法。但该算法在实际应用中存在一些制约因素：（1）反演过程中很难获得方程 的解析解，只能通过迭代法获取其近似解，计算量大，彳艮难将其业务化。（2）有些像元出现不收敛的情况误差较大.（2）不同卫星传感器之间会存在图像的配 准误差，反演过程中影响精度，导致部分反演结果误差较大。（3）