气溶胶光学厚度.docx
《气溶胶光学厚度.docx》由会员分享，可在线阅读，更多相关《气溶胶光学厚度.docx（13页珍藏版）》请在冰豆网上搜索。
气溶胶光学厚度
第2章气溶胶光学厚度反演的原理和方法
气溶胶光学厚度（AerosolOpticalDepth）简称AOD，定义为介质的消光系数在垂直方向上的积分，描述的是气溶胶对光的消减作用⑺。
它是气溶胶最重要
的参数之一，表征大气浑浊程度的关键物理量，也是确定气溶胶气候效应的重要因素。
。
通常高的AOD值预示着气溶胶纵向积累的增长，因此导致了大气能见度的降低。
现阶段对于AOD的监测主要有地基遥感和卫星遥感两种方法。
其中地基遥感又有多种形式：
多波段光度计遥感、全波段太阳直接辐射遥感、激光雷达遥感等。
其中多波段光度计遥感是目前地基遥感研究中采用的最广泛的方法。
美国NASA和法国LOA-PHOTONS联合建立的全球地基气溶胶遥感观测网AERONET所使用的就是多波段太阳光度计（Sun/SkyPhotomererS，在全球共布设1217个站点长期观测全球气溶胶的光学特性，积累了大量的AOD数据，并用作检测气溶胶光学厚度反演精度的标准。
而近年来卫星遥感技术的快速发展，多种传感器被用来研究气溶胶特性，加上经济发展带来的大气污染问题使得利用卫星遥感资料反演AOD成为热门课题。
2.1气溶胶光学厚度反演的基本原理
大气光学厚度是指沿辐射传输路径单位截面上气体吸收和粒子散射产生的总消弱，是无纲量值。
在可见光和近红外波段，它可以由下列公式计算得出：
（0=,m（0v?
<1（0V<<2（，）■.，〔（，）■.a（，）（2-1）
其中（’）表示大气总的光学厚度，,m（'）表示整层大气的分子散射光学厚度，
..1（）表示氧气的吸收光学厚度，.2（'）表示臭氧的吸收光学厚度，表示
水汽的吸收光学厚度，a（'）表示气溶胶光学厚度［21;221
卫星遥感反演大气气溶胶是利用卫星传感器探测到的大气顶部的反射率，也
称为表观反射率，可以表示为［23］:
（2-2）
其中，L表示卫星传感器探测到的辐射值，Fs表示大气上界太阳辐射通量，
展表示太阳天顶角压的余弦值。
P与地表二项反射率之间的关系可以表达为：
:
（戈..，-s,）Fd（夬）T（）,）
i■（•.，rs,'）=:
匕（*•，Vs，*（2-3）
1-SP
其中，K表示传感器天顶角，,s表示太阳天顶角，表示太阳方位角和卫
星方位角确定的相对方位角；嘉（）.’茂，）表示由大气分子和气溶胶散射造成的路径辐射，它与地表状况无关；Fd（却表示地表反射率归一化为零时总的向下辐射通量，也可以称为总的向下透过率，由于气溶胶粒子对太阳光的吸收和散射作用，它的值小于1.0；T（n）是向上进入卫星传感器视场方向的总透过率，S是
大气后向散射比。
在单次散射近似中，路径辐射嘉（）’爻，）与气溶胶光学厚度a和单次散射反射率*之间的关系如下[24]:
0.apaG.,茂，）
鼻（’，，）=：
m（丁，茂，）（2-4）
4一S
S
其中讣（K，ts,）是分子散射造成的路径辐射，它取决与大气模式，=表
示传感器天顶角的余弦值，」表示太阳天顶角的余弦值。
在式（2-3）中，Fd（ds）、
S
TO,）和S取决于-，o，a和Pa（Ts,）。
-'0.aPaG..,為，）
假设地表是均匀朗伯表面，大气垂直均匀变化，将式（2-4）代入（2-3）得：
（2-5）
P*（TV,es,©）=Pm（g，a,©）
「（J,Vs,）Fd（JS,0，.a,Pa）T（JS,，0,.a,Pa）
上式中r（戈,戈，）为假设的朗伯体特性的地表反射率，卫星传感器接收到的表观发射率r,:
1s,■-）既是地表反射率：
、（）,*,'）的函数，又是气溶胶光
学厚度谚的函数［25］。
表观反射率P*（g8s,）以及太阳和传感器的几何参数（ev,03,）可以从卫星遥感资料中获取，假如可以得知地表反射率p*^v,es^），并用气溶胶类型和大气模式来确定-.0和Pa（，,■-）的相关参数，理论上就可
以计算得出地面上空的气溶胶光学厚度ao反之，若已知地面上空气溶胶光学厚
度為、气溶胶类型以及大气模式，也可以反演出地表反射率P*（S，0s，）［26］o
由式（2-3）可以得出反演气溶胶的最优条件是地表反射率低且光谱波段波
波长较短。
在地表反射率角度的情况下，气溶胶散射引起的路径辐射
Pa（日W6s，©）（与地表状况无关）对表观反射率P（日v，Os，©）起主要作用，此时反演气溶胶光学厚度误差较小；而在地表反射率较大的情况下，地表贡献项：
（^▼茂，）Fd（「s）T（戈）对表观反射率：
、*（）,寸$「）影响较大，此时反演精度
1-s"
较低［27］o
2.2气溶胶光学厚度反演的主要方法
自20世纪70年代中期开始，利用卫星数据反演气溶胶光学厚度的研究已经有40年的历史，反演的方法有单通道算法、多通道算法、暗像元法、结构函数法、深蓝算法、多星协同反演法、海陆对比法、多角度偏振法、热辐射对比等［28］o目前有代表性的常用气溶胶光学厚度反演算法有两种：
一种是通过路径辐射项求
取气溶胶光学厚度的暗像元法，另一种是通过透过率求取气溶胶光学厚度的对比法。
2.2.1暗像元法
暗像元算法英文全称为DenseDarkVegetation,简称DDV，它是通过路径辐射项来计算气溶胶光学厚度。
由于地表物体的复杂多样性造成反射率变化范围很大，很难从辐射值中分理处辐射项，如果想通过辐射项来获取气溶胶信息，就必
须使地表辐射值较小且能确定其精确值，这样就能够最大限度的消除地表反射率的不确定性带来的影响。
在卫星影像中，大量浓密植被区由于在可见光波段反射率极低（约为0.01~0.02）,它们被称作暗像元。
1988年Kaufman等利用大多数地物在红（0.60~0.68^m蓝（0.40~0.48卩波段反射率低的特性，根据归一化植被指
数（NDVI）或近红外通道（2.1^m的表观反射率进行暗像元的识别，并假定已知这些暗像元在红蓝通道的地表反射率，依据一定的关系反演气溶胶光学厚度［29］。
通过大量的卫星影像资料，考虑到多种地表覆盖物，拟合得到红（0.66蓝（0.47和中红外通道（2.1地表反射率的关系［29］:
;'2.1
（2-6）
■2.1
4
对于2.1yn通道卫星观测表观反射率几乎不受气溶胶影响，其值接近地表反射率，因此可以用2.1通道的表观反射率代替地表反射率，根据式（2-6）计算出红蓝通道的地表反射率。
再假定合适的气溶胶模型，就可以计算出气溶胶光学厚度。
然而对于干旱、半干旱以及冬季城市等高发射率地区，用暗像元法来反演气溶胶光学厚度还存在不少困难，在以上地区由于地表的非均一性使确定地表反射率的精确值十分不易［30］。
当地表反射率升高时，气溶胶的指示作用降低［31］。
一般情况下，当地表反射率较低时，传感器接收到的辐射值随着气溶胶的增多而迅速变大，暗像元法就是利用浓密植被在红、蓝波段的地表反射了和气溶胶的这种
关系来反演光学厚度；随着地表反射率的不断增大，接收到的辐射值随气溶胶的增多而增大速度减缓，当地表反射率增大到某一程度时，辐射值将不随气溶胶的增多而增大，甚至会出现降低的趋势。
为了使暗像元方法的应用更为广泛，2002年Kaufman等通过大量的数据验证，对以前的方法进行了扩展：
对于星下点，暗像元法适用于中红外通道表观反射率小于0.4的区域；对于非星下点，需要考
虑太阳和卫星的几何参数，适用范围可以扩展到中红外通道表观反射率小于0.125（1=•1'匚°）的地区（其中丿为卫星天顶角的余弦值，%为太阳天顶角的余弦值），红蓝通道的地表反射率依然遵循式（2-6）的函数关系。
2.2.2扩展的暗像元法（V5.2算法）
扩展后的暗像元法是由Levy和Remer等人提出的称为V5.2的气溶胶反演算法［32;33］。
与原来的算法相比，它主要有两方面的改进：
首先提出VIS/SWIR（可
见光和中红外通道的反射率比值）不再是一个常量，而是一个函数关系；其次考虑了NDVI（植被指数）对于地表反射率的影响，红蓝通道和中红外通道的地表反射率关系如式（2-7）所示［34］:
（2-7）
「0.66=f（『2.1）
「0.47二9（「0.66）
该方法中的VIS/SWIR的变化跟不但跟地表植被指数有关，根据Remer和
Gatebe等的研究表明它还与太阳天顶角、传感器方位角、散射角有关，尤以散射
角的影响最大［35;36］。
它可以表示为：
1
（2-8）
其中二表示观测天顶角,
状表示太阳天顶角，表示太阳方位角和卫星方位
角的相对方位角。
根据式（
2-8）可将式（2-7）扩展为:
s
:
'0.66
s
:
'0.47
二f（‘2.1）
=g（:
'0.66）-
s.
P2.12slope0.66.21+yint0.66,21
s
Po.66slope0.47.0.66+yint
0.470.66
（2-9）
C-?
-cos—（-cosncosrsin^sincos）
其中
slope0.66.21=slopeN.D；SWR1+0.00^-0.27
（2-10）
yint0.662.〔=0.0002旳+0.003
Slope0.470.66=0.49
yin（t470.66=°.005
NDVIswir可表示为:
NDVSwir=（—.2-;?
2.1）（;?
1.^;?
2.1）（2-11）
其中，汀2和分别表示MODIS第5波段和第7波段接收到的辐射值，
当NDVSwir大于0.6时，表示该地区的植被状况较好；当NDVIswir小于0.2时，表
示该地区植被稀疏［37］。
当NDVIswirV0.25时，slope黑答=0.48；
当NDVIswirA0.75时，slope：
囂笄=0.58；
当o.25兰nDVIswir兰0.75时，slope播笄=48+0.2*（nDVSwir-0.25）
其中slope和yint为统计方程的经验系数，slope为偏移量，yint为截距。
将（2-10）代入（2-9）可得出：
a.植被指数小于0.25的情况：
ss——
（2-12）
p0.66=p2j"0.48+0.002。
-0.27）+0.0002泪+0.033
ss
P0.47=P0.66汉°.49+0.005
b.植被指数大于0.75的情况：
ss
（2-13）
P0.66=P2M"0.58+0.002。
一0.27）+0.0002対+0.033
ss
P0.47=P0.66汉0.49+0.005
c.植被指数介于0.25和0.75之间的情况：
ss
（2-14）
P0.66=P2.1汉[0.48+0.2NDVIswir-0.25）+0.002。
一0.27]
0.00025。
0.033
ss
P0.47=P0.66汉0.49+0.005
2.2.3结构函数法
此时应用暗
对于旱季和中高纬度地区的冬季，植被覆盖少，地物反射率高，目标法会造成一定的误差，因此，针对陆地亮地表的情况，发展了结构函数法来反演气溶胶光学厚度。
结构函数法也称为对比法，是早期用来研究陆地污染气溶胶采用的卫星遥感算法。
它跟暗像元法一样采用的是红、蓝通道数据。
对比法需要假定在同一地区同一时间段内地表反射率不变［38］。
然后找出污染极小或无污染的一天作为清洁
日”，以它的气溶胶数据作为参考反演污染日”的大气气溶胶光学厚度。
TanreftHolben在采用对比法进行气溶胶反演的时候，假定在同一卫星观测时刻，对于两个相邻的像元大气是均一的，并在反演的过程中引入了结构函数的概念，使得
地表反射率的确定变得容易［39;40］。
该算法主要利用的是表观反射率的地表贡献项来反演气溶胶光学厚度。
在暗像元方法不适用的地区，结构函数法为气溶胶光学厚度的反演提供了一条新的途径。
在反演过程中，需要通过一段时间内卫星观测数据的分析，选定其中无污染（气溶胶光学厚度极小）的一天作为清洁日”然后通过地面观测或者其他途径确定这一天的气溶胶光学厚度值，作为背景的气溶胶信息。
假定地表目标无变化，通过透射函数的变化就能获取其它污染日”的气溶胶光学厚度。
由于总透射函数对于相函数细节不敏感，结构函数法对散射相函数有较大的独立性，获取气溶胶光学厚度主要根据单次散射反照率和不对称因子。
由式（2-3）得出相邻两点的表观反射率差值为：
:
?
（i,j）=
二'（i,j）T（0S）（0V）
（2-15）
地表结构函数为[40]:
M（d）=
nm_d2
n（m_d）H（%「i，jd）
（2-16）
改进后的地表结构函数为[41]
nm_d
2I222
M（d）一3（n-d）（m-d）ij丿°i,j一:
i,jd>（:
?
^_:
?
i-d，j）'（A，j一;?
i，jd）]
（2-17）
在卫星高度上的结构函数为:
M2（d）二M（d）T2（7is）T2Gv）
（2-18）
假定清洁日”1t和污染日”2地表特征保持不变，则M（d,t1）=M（d,t2），观测结构函数关系根据式（2-12）得：
M（d,t1）=Tl9s（t1）,6v（t1）,Ta（t1）]（219
M（d,t2）Xds（t2）,rv（t2）,.a（t2）]
由上式可得在已知清洁日”的气溶胶光学厚度的情况下，污染日”的气溶胶光学
厚度就可以由卫星观测值获取。
224双星协同反演算法
近年来随着卫星遥感技术的不断发展，对地观测卫星携带的传感器种类不断增多，这对于不同平台之间的结合研究起到了推动作用。
为了解决气溶胶光学厚度反演中未知
参数过多导致气溶胶光学厚度信息难以提取的问题，唐家奎等提出了基于Terra
和Aqua双星的MODIS遥感数据协同反演算法［17］。
该算法具有以下几种优点：
（1）不需要事先假设气溶胶的类型等参数。
（2）不需要假设或利用其它统计资料预先估计地表反射率，双星协同算法是将真实的表观反射率作为反演气溶胶光学厚度的变量之一与其同时反演。
（3）不收地表反射类型的影响，可应用高反射率的陆地气溶胶反演。
（4）可以同时获取两次过境时间的气溶胶信息，有利于分析气溶胶的动态变化规律［17］。
基本原理如下：
在整个大气层范围内对微元体积大气的辐射性质进行积分称为地球一大气系统的辐射传输模型。
根据以前学者的研究分析［42］，地球一大气系统之间的辐射传输方程可以表达为：
COSI匚
Pdz4兀
I（z,■）（z,r,r）d•—（kOi'（z,r）
（2-20）
其中表示传感器天顶角，「表示空气密度，二表示散射系数，z表示咼度,
（z,r,r）表示散射相函数，•「表示固体角，r表示方位（天顶角，方位角），
I（z,r）表示在z高度r方向的辐射强度。
方程（2-20）中不能够得到|'（z,r）的
解析解，只能够获得方程的近似解，Xue等通过用多个不同变量将方程简化后得到地表反射率L和表观反射率L之间的关系如下［43］:
（a_b）E/sec野
（2-21）
_（Lb-a）a（1-L）e0_
（a_b）EJsec&
（Lb-a）b1-L）e0
其中■表示后向散射系数（通常取值为0.1），a二sec二'，b=2，表示表示气溶胶光学厚度。
在双星协同反演模型中，仅考虑了空气分子和气溶胶粒子散射两种影响因素，将大气气溶胶光学厚度分为大气分子的瑞利散射和大气气溶胶粒子的散射两部分：
（2-22）
根据Linke等的研究大气分子的瑞丽散射近似解为［44］:
（2-23）
4.09
•mC：
）二0.00879一
利用？
ngstr?
m浊度公式可以求出气溶胶粒子散射光学厚度为:
（2-24）
将式（2-24）和（2-23）以及（2-22）代入（2-21）中，可以得出地表反射率L以及？
ngstr?
m浊度系数一：
和波长指数:
-之间的关系方程。
首先假设在两次卫星过境时间间隔内，地表反射率不变，且气溶胶类型保持不变，这样波长指数：
•不变，只有气溶胶粒子浓度变化，即浊度系数［是变化的。
然后利用0.47,0.55以及0.66^m三个通道实现双星协同反演算法
但该算法在实际应用中存在一些制约因素：
（1）反演过程中很难获得方程
的解析解，只能通过迭代法获取其近似解，计算量大，很难将其业务化。
（2）
有些像元出现不收敛的情况误差较大.
（2）不同卫星传感器之间会存在图像的配准误差，反演过程中影响精度，导致部分反演结果误差较大。
（3）不同卫星过境时间有一定间隔，气溶胶以及大气中其它成分的性质数量会发生变化，导致部
分像元反演结果无法收敛。
2.2.5高反差地表法
20世纪70年代Kondratyev等提出了利用高反差地表发反演晴空条件下陆地上空的气溶胶光学厚度［45;46］。
利用该方法反演陆陆地上空气溶胶光学厚度需要在两个空间位置相近的区域，假设大气光学特性不变，选择明暗两种对比明显的像元进行反演。
这里要求的明暗像元是相对而言，并不局限于浓密植被覆盖的特定区域，在一定程度上解决了亮地表上空反演的问题。
但实际应用中，该方法要求传感器误差极小，且地表光谱反射值的差异要足够大，实现业务化反演比较困难。
2.3本章小结
本章节主要介绍了利用遥感影像数据反演气溶胶光学厚度的基本原理，并介绍暗像元法、V5.2算法、结构函数法、双星协同反演算法以及深蓝算法和高反差地法。
详细讲述了各自算法的原理，并描述了它们的应用范围、反演效果以及优缺点。
总体来说，暗像元法应用最为广泛，它利用红、蓝以及中红外通道信息，根据它们之间地表反射率的关系进行气溶胶光学厚度的反演，但仅适用于地表反射率较低时的情况，相对来说结构函数法在针对城市及其周边地区的气溶胶污染监测有一定的优势。
双星协同反演算法虽然不受地表类型限制，适用于城市等亮地表的气溶胶光学厚度反演，但由于该算法计算量大而且容易受到外界环境（如传感器间隔时间内温度、湿度）等的影响部分像元存在较大误差。
高反差地法获取的反演结果只是相对气溶胶光学厚度，很难实现业务化。
而V5.2算法即具有暗像元算法的优点，又在一定程度上减小了“亮”地表带来的影响。
本文的研究区域—北京市位于华北大平原北部，四季分明，在春、夏、秋季节城市及周边植被覆盖良好，冬季干旱城区地表反射率高，若利用暗像元法会存在一定误差，而扩展的暗像元法在一定程度上可以解决此问题。
升级会员 