鄒健++崔婉哲

摘要： 氣溶膠是氣候變化和大氣污染的重要影響因素。為了研究渭南地區(qū)氣溶膠光學(xué)厚度等性質(zhì)，使用2010年7月16日MODIS-L1B數(shù)據(jù)，利用黑暗像元法，在ENVI5.0系列軟件中進(jìn)行反演操作，反演出渭南地區(qū)當(dāng)日的氣溶膠光學(xué)厚度空間分布。通過對反演結(jié)果做進(jìn)一步的空間分析。研究結(jié)果表明：渭南的平均AOD為0.45；渭南中部臨渭區(qū)的AOD高于四周；以108.7°為界限，東邊的AOD值明顯小于西邊；渭南夏季AOD偏低。

Abstract： The aerosol is the important influence factor of climatic variation and atmospheric pollution. In order to study the optical thickness and other properties of aerosol in Weinan area， the MODIS-L1B data in July 16， 2010 and the Dark pixel method are used to carry out the inversion operation in ENVI5.0 series software. The spatial distribution of aerosol optical thickness in Weinan area was obtained. The inversion results are further analyzed. Research results show that the average AOD of WeiNan is 0.45， the AOD in Linwei District of central Weinan is higher than that of the surrounding area， and based on 108.7° line， the AOD value of east is evidently less than that in west and AOD value is low in summer.

關(guān)鍵詞： 濃密植被法；氣溶膠光學(xué)厚度；MODIS；ENVI；反演

Key words： dense vegetation method；aerosol optics thickness；MODIS；ENVI；inversion

中圖分類號：X513 文獻(xiàn)標(biāo)識碼：A 文章編號：1006-4311（2017）12-0175-03

0 引言

大氣與懸浮在其中的固體和液態(tài)且小于10的微粒共同組成的多相體系稱為氣溶膠[1]。大氣氣溶膠光學(xué)厚度（Aerosol Optical Depth，AOD）的物理意義是沿輻射傳射路徑，單位截面上氣溶膠吸收和散射對太陽輻射產(chǎn)生的總削弱，它和垂直方向上的總氣溶膠濃度有關(guān)[2]。研究氣溶膠的光學(xué)厚度可以間接反映反演地區(qū)大氣質(zhì)量，對于研究陜西關(guān)中地區(qū)人們健康、地區(qū)氣候效應(yīng)以及大氣污染都有一定參考價(jià)值。首先，氣溶膠通過對太陽短波輻射的散射，以調(diào)節(jié)地—?dú)庀到y(tǒng)的反射率，致使地—?dú)庀到y(tǒng)降溫，因而對地—?dú)庀到y(tǒng)的輻射平衡產(chǎn)生間接影響；其次，氣溶膠是?！憽?dú)庀到y(tǒng)中一種重要的組成成分，是大氣輻射傳輸過程中三大活躍要素之一。它通過三個主要機(jī)制來影響氣候：①氣溶膠通過對吸收和散射太陽輻射，來影響太陽到達(dá)地球的熱量；②氣溶膠能一定程度變云屬性，進(jìn)而影響到達(dá)地球的太陽輻射量；③氣溶膠能對大氣中各類化學(xué)組成產(chǎn)生間接影響，通過改變這些化學(xué)成分及其濃度來間接影響氣候[3-5]。

傳統(tǒng)的氣溶膠測量主要是地面基站測量，這種測量方法對氣溶膠的監(jiān)測具有準(zhǔn)確、數(shù)據(jù)可信度高等特點(diǎn)。缺點(diǎn)是不具有連續(xù)性、難以從大尺度的范圍來對氣溶膠空間分布、運(yùn)動軌跡進(jìn)行研究；其次就是成本高。目前也有大量針對氣溶膠光學(xué)性質(zhì)的研究，MODIS為全球提供了針對分辨率為10km的氣溶膠產(chǎn)品，但對于研究城市范圍來說，10km分辨率的氣溶膠產(chǎn)品不能滿足需要，所以1km分辨率的衛(wèi)星遙感反演，可以更好地反映城市和小區(qū)域的氣溶膠狀況[2]。本研究使用MODIS數(shù)據(jù)對渭南地區(qū)的氣溶膠光學(xué)厚度反演，可以客觀反應(yīng)渭南地區(qū)AOD空間變化，對改善生態(tài)、促進(jìn)社會經(jīng)濟(jì)可持續(xù)發(fā)展具有重要意義[6]。

1 研究區(qū)概況及數(shù)據(jù)

渭南市，隸屬于陜西省，位于東經(jīng)108°50′-110°38′和北緯34°13′～35°52′之間，地處陜西關(guān)中渭河平原東部、關(guān)中地區(qū)，是陜西省的“東大門”?？偯娣e約13134平方千米，人口約560萬。渭南屬暖溫帶半濕潤半干旱季風(fēng)氣候，四季分明，光照充足，雨量適宜。渭南地勢以渭河為軸線，形成南北兩山、兩塬和中部平川五大地貌類型區(qū)型。

本文選取MODIS TERRA衛(wèi)星數(shù)據(jù)產(chǎn)品 MODIS L1B數(shù)據(jù)MOD021KM.A2010228.0310.005.hdf。利用經(jīng)典的黑暗像元法，對渭南市進(jìn)行1km空間分辨率氣溶膠光學(xué)厚度反演，并分析渭南市的氣溶膠空間變化特征。

2 氣溶膠光學(xué)厚度反演算法

2.1 反演原理

衛(wèi)星遙感反演氣溶膠光學(xué)厚度基本原理見式（1）：

Rref（？滋，？茲，？滋0，？茲0）=？籽a(bǔ)（？滋，？茲，？滋0，？茲0）+T（？滋0）T（？滋）？籽b/（1-？籽bs）

（1）

其中Rref是衛(wèi)星觀測到的整層反射率，是關(guān)于氣溶膠厚度的函數(shù)，？籽a(bǔ)為整層大氣層反射率，？滋為觀測天頂角的余弦，？茲為觀測天頂角的方位角，？滋0為太陽天頂角的余弦，？茲0為太陽天頂角的方位角，T（？滋0）為太陽到地面的大氣層總透射率，T（？滋）為從地面到衛(wèi)星大氣層總的透過率（直射+漫射），s為大氣的下墊面反射率，系數(shù)1/（1-？籽bs）代表地面和大氣層多次散射作用。在假定一定的大氣氣溶膠模型前提下，如果知道了下墊面的反射率？籽b，根據(jù)公式（1）即可反演得到氣溶膠的光學(xué)厚度[7]。

2.2 DDV（暗像元法）地表反射率估算

Kaufman等最先提出DDV（暗像元法），Kaufman等通過大量觀測實(shí)驗(yàn)研究指出，在清潔的大氣條件下，由于綠色植被對太陽光的強(qiáng)烈吸收，植被的反射率在紅波段通道、藍(lán)波段通道和2.1～2.6？滋m通道呈現(xiàn)出線性關(guān)系。由于氣溶膠主要影響可見光通道的圖像，而2.1～2.6？滋m通道的圖像幾乎不受氣溶膠的影響，反映的是地面的特征。根據(jù)地表植被覆蓋區(qū)的紅（0.66？滋m）、藍(lán)（0.47？滋m）通道和2.1～2.6？滋m通道的線性關(guān)系，紅、藍(lán)通道的地表反射率可以從2.1？滋m通道的地表反射率估算出來[7]，即：

p0.47=p■■/4，p0.66=p■■/2（2）

2.3 LUT查找表

它是基于6S大氣輻射傳輸模型計(jì)算的查找表（LUT），并利用地面的暗目標(biāo)自動反演陸地大氣氣溶膠光學(xué)厚度。LUT查找表是利用IDL調(diào)用6S輻射模型得到的，采用的是通用的參數(shù)，3-9月期間都可以用這一個查找表進(jìn)行氣溶膠反演[8]。

3 數(shù)據(jù)處理

利用MODIS的L1B數(shù)據(jù)反演氣溶膠的流程如圖1所示。

首先是數(shù)據(jù)預(yù)處理：幾何糾正，對MODIS數(shù)據(jù)進(jìn)行幾何糾正，具體包括發(fā)射率糾正、反射率糾正，并對發(fā)射率、反射率合成剪裁生成，以及對角度數(shù)據(jù)的的幾何糾正；云監(jiān)測，運(yùn)用IDL語言編寫的拓展工具，進(jìn)行云監(jiān)測，其基本原理是利用MODIS某對云檢測有效果的波段設(shè)置閾值把云的分布給提取出來。

其次就是數(shù)據(jù)可視化：調(diào)用通過6S模型建立的通用LUT查找表，得到包括渭南市的氣溶膠光學(xué)厚度圖，利用矢量 “渭南市行政邊界”裁出渭南地區(qū)，生成圖2。

4 反演結(jié)果與分析

4.1 渭南市光學(xué)厚度空間分析

在渭南市利用MODIS L1B數(shù)據(jù)進(jìn)行了1km空間分辨率氣溶膠光學(xué)厚度的反演，提取渭南地區(qū)的平均光學(xué)厚度得出渭南2010年7月16日11時(shí)的氣溶膠平均光學(xué)厚度為0.45。從圖2提取北緯34.5度附近提取橫向氣溶膠的光學(xué)厚度，見圖3。

從圖2和圖3我們得到：

①以109.1°為界限，東邊的AOD值明顯小于西邊。其原因和關(guān)中地區(qū)特有的地形的位置因素有關(guān)，渭南處于關(guān)中渭河平原東部，關(guān)中地區(qū)的城市呈現(xiàn)一種橫向帶狀分布狀態(tài)，渭南距離寶雞、西安等關(guān)中地區(qū)的工業(yè)發(fā)達(dá)城市有一定距離，本身化工、金屬冶煉等相比之下較不發(fā)達(dá)，所以以108.7°邊界為界限，東西兩邊AOD值具有明顯的差異。

②渭南地區(qū)7月16日AOD平均值本身較低。這是因?yàn)榇藭r(shí)渭南處于夏季本身降水較多，空氣濕潤，水汽對氣溶膠粒子具有吸附作用；太陽直射能力強(qiáng)，植被茂盛，植被也對氣溶膠粒子具有一定的吸附作用；而且周邊區(qū)域環(huán)境都較好。綜合以上幾種原因，渭南地區(qū)夏季AOD值偏低。

③渭南市內(nèi)AOD值分布不均衡，渭南中部地區(qū)的AOD值明顯高于周邊，其原因和渭南市的工業(yè)和城市分布有著極大的關(guān)系，渭南中部為臨渭區(qū)，是渭南市人口工業(yè)集中區(qū)，所產(chǎn)生的污染物及熱量比周邊多，進(jìn)而增大了渭南中部的AOD值。

4.2 反演優(yōu)缺點(diǎn)

①本文介紹了利用衛(wèi)星遙感圖像在ENVI5.0下的反演一個區(qū)域氣溶膠光學(xué)厚度的的方法，利用該方法，能較全面的反映渭南整個地區(qū)氣溶膠光學(xué)厚度的分布，相較于地面實(shí)測數(shù)據(jù)來說，更能從全局以及跨區(qū)域分析渭南氣溶膠光學(xué)厚度的分布、趨勢。

②基于Kaufman等提出的濃密植被法，具有反演原理簡單，流程較易等特點(diǎn)，但是由于該算法是基于一些列植被指數(shù)，而北方地區(qū)普遍植被情況總體較差，僅能在3-9月進(jìn)行反演操作。

③反演結(jié)果與實(shí)際的AOD值存在一定的偏差。如，利用6S輻射模型建立的LUT查找表，由于氣溶膠類型的不確定可能會帶來很大的誤差；在做反射率糾正時(shí)，氣溶膠下方的植被覆蓋區(qū)的地表的紅（0.66？滋m）、藍(lán)（0.47？滋m）通道的反射率是從2.1？滋m通道的反射率估算出來，由于是線性相關(guān)公式，本身與實(shí)際情況具有一定誤差。

5 展望與建議

為了對小區(qū)域氣溶膠時(shí)空分布、運(yùn)動軌跡等問題進(jìn)行研究，建議采用衛(wèi)星遙感反演城市和小區(qū)域氣溶膠狀況，目前這種方法已經(jīng)在我國廣泛應(yīng)用。另外，建議在全國氣象、環(huán)保等部門加深這種方法的推廣力度，也鼓勵相關(guān)專業(yè)學(xué)習(xí)這種方法。

參考文獻(xiàn)：

[1] 李成才，毛杰泰，劉啟漢.用MODIS遙感資料分析四川盆地氣溶膠光學(xué)厚度時(shí)空分布[J].應(yīng)用氣象學(xué)報(bào)，2003，14（1）：1-7.

[2]劉昆，倪長健.成都地區(qū)MODIS衛(wèi)星AOD產(chǎn)品的適用性驗(yàn)證及應(yīng)用[J].成都信息工程學(xué)院學(xué)報(bào)，2015，30（3）：298-302.

[3]王耀庭，王橋，王艷姣，等.大氣氣溶膠性質(zhì)及其衛(wèi)星遙感反演[J].環(huán)境科學(xué)研究，2005，18（6）：27-33.

[4]莫招育，陳志明，黃囧麗，等.南寧地區(qū)秋季一次連續(xù)灰霾過程的污染特征及成因分析[J].桂林理工大學(xué)學(xué)報(bào)，2015，35（3）：540-548.

[5]李成才，毛節(jié)泰，劉啟漢，等.利用MODIS衛(wèi)星和激光雷達(dá)遙感資料研究香港地區(qū)的一次大氣氣溶膠污染[J].應(yīng)用氣象學(xué)報(bào)，2004，15（6）：641-650.

[6]韓晨.基于SPOT-VGT的陜西黃河濕地保護(hù)區(qū)植被覆蓋研究[J].湖南農(nóng)業(yè)科學(xué)，2016（55）：859-862.

[7]仁佳，王振會，孫林，等.基于MODIS數(shù)據(jù)反演江浙皖地區(qū)氣溶膠光學(xué)厚度[J].環(huán)境科學(xué)與技術(shù)，2010，33（8）：167-171.

[8]王玲，田慶久，李珊珊.利用MODIS資料反演杭州市500米分辨率氣溶膠光學(xué)厚度[J].遙感應(yīng)用，2010（3）：50-54.

[9]趙英時(shí).遙感應(yīng)用分析原理與方法[M].北京：科學(xué)出版社，2012：19-28.