楊洵，張振超，茍青松

內(nèi)蒙古農(nóng)業(yè)大學水利與土木建筑工程學院

基于MODIS數(shù)據(jù)的呼包地區(qū)氣溶膠反演

楊洵，張振超，茍青松

內(nèi)蒙古農(nóng)業(yè)大學水利與土木建筑工程學院

氣溶膠對地球環(huán)境具有重要影響，故而備受關注。為了研究呼包地區(qū)氣溶膠特性，文章利用Aqua衛(wèi)星MODIS數(shù)據(jù)作為數(shù)據(jù)源，基于6S模型，利用暗像元法反演呼包地區(qū)氣溶膠光學厚度，并結合地面觀測數(shù)據(jù)進行精度驗證。反演結果表明，呼包地區(qū)整體空氣質(zhì)量良好，利用MODIS數(shù)據(jù)反演AOT從而監(jiān)測呼包地區(qū)環(huán)境變化趨勢較為可行，能在一定程度上為環(huán)保部門提供決策依據(jù)。

MODIS；氣溶膠光學厚度；暗像元法

引言

氣溶膠是液態(tài)或固態(tài)等多相微粒在空氣中的懸浮體系，粒徑大小范圍在10-3～10μm，有自然成因和人為成因。因其對大氣能見度[1]、空氣質(zhì)量乃至整個地氣系統(tǒng)皆有重要影響，并且會間接地改變大氣化學成分從而危害人體健康，因此，對氣溶膠進行相關研究顯得尤為重要，也備受人們重視。

雖然傳統(tǒng)基站測量監(jiān)測可溶膠數(shù)據(jù)準確、可信度較高，但不能從宏觀角度大尺度地反映氣溶膠空間分布及其運行規(guī)律，且成本較高，遙感衛(wèi)星反演正好彌補了這些缺陷。目前主要通過反演氣溶膠對光的削減作用的量度—氣溶膠光學厚度（AOT）來研究其相關特性，用氣溶膠光學厚度來描述氣溶膠消光性及其結果可以作為評估大氣污染程度指標之一。2017年內(nèi)蒙古自治區(qū)在政府工作報告中關于對大氣污染防治提出了實施工業(yè)污染源全面達標排放計劃，因此，使用MODIS數(shù)據(jù)對呼包地區(qū)進行氣溶膠反演對于監(jiān)測AOT空間分布及其變化趨勢具有相當重要的現(xiàn)實意義。

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

本文選取呼和浩特和包頭地區(qū)作為研究區(qū)域，區(qū)域經(jīng)緯度為E：109°25′～112°17′，N：39°37′～41°30′。呼包地區(qū)屬大陸性干旱氣候，晝夜溫差較大，年均降水量約400毫米，冬季漫長嚴寒，夏季短暫炎熱，該區(qū)屬典型侵蝕構造地形，地貌類型較復雜。

本文反演數(shù)據(jù)來源于美國LAADS DAAC的MODIS MYD021KM數(shù)據(jù)，由于采用采取的方法為暗像元法，因其比較適合夏季氣溶膠反演，故下載的數(shù)據(jù)時間段為7—8月。

2 AOT反演理論

2.1 反演原理

陸地上氣溶膠遙感反演發(fā)展于衛(wèi)星傳感器探測到的大氣上界的表觀反射率ρ?，它的表達式如下所示：

上式中L是大氣上界輻射,F0大氣外界太陽輻射能量，μ0太陽天頂角的余弦[2]。

ρ?與地表反射率ρ(θ,θ0,?)之間具有以下關系：

θ為衛(wèi)星天頂角,?相對方位角；ρa(θ,θ0,?)為路徑輻射,θ0為太陽天頂角；Fd(θ0)為向下透過率,其值小于1.0；T(θ)是向上總透過率,S是大氣后向散射比,ρ''是在觀測角和入射角上平均的地表反射率。

為了從表觀反射率反演氣溶膠光學厚度,需要合理假定氣溶膠模型，以提供氣溶膠散射相函數(shù)ρa(θ,θ0,?)，在確定了上述參數(shù)的輻射傳輸模型后,理論上就可以反演得到其上空的氣溶膠光學厚度τλ，這就是氣溶膠光學厚度遙感反演的基礎[3]。

2.2 DDV估算反射率

考慮到地表覆蓋的多樣性，擬合得到三個不同區(qū)域的波段與地表之間的關系，這三個波段分別為：藍光波段（0.47 μm）、紅光波段（0.66μm）和中紅外波段（2.1μm）。又因氣溶膠在2.1μm波段幾乎是透明，故在2.1μm處的表觀反射率相當于地表反射率，所以可以用2.1μm波段表示其他2個波段，關系式如下：

ρR=ρ?

2.1/2，ρB=ρ?2.1/4式中ρR是紅波段地表反射率,ρB是藍波段地表反射率,ρ?2.1是波段等于2.1μm的表觀反射率。

2.3 查找表生成

已知ρa(θ,θ0,?)和表觀反射率ρ?2.1，地表反射率ρ，S，T時，利用IDL調(diào)用6S輻射模型得到氣溶膠反演的查找表，查找表建立起來了θ，θ0，?，ρa，S，T，氣溶膠光學厚度（AOT）之間的對應關系。

2.4 反演流程

利用MODIS L1B數(shù)據(jù)結合角度數(shù)據(jù)對研究區(qū)氣溶膠厚度進行反演，在對發(fā)射率、反射率等各項數(shù)據(jù)進行預處理的基礎上，再利用IDL語言調(diào)用6s模型生成查找表，對反演結果進行地面觀測數(shù)據(jù)驗證，最后對反演結果進行時空變化規(guī)律分析。

2.5 反演結果及精度分析

在經(jīng)過各項數(shù)據(jù)處理基礎上，利用遙感技術對呼包地區(qū)7月2日、7月15日、8月1日、8月15日的MODIS數(shù)據(jù)進行AOT反演，反演結果如下圖所示：

圖1 2016年7月2日反演結果

圖2 2016年8月15日反演結果

AOT遙感反演結果需結合地面觀測數(shù)據(jù)進行驗證，氣溶膠光學厚度地面觀測站在呼和浩特和包頭地區(qū)設有AOE—Baotou站。通過上述影像，根據(jù)AOE—Baotou站的經(jīng)緯度信息讀取AOT的值與AOE—Baotou站同時間段數(shù)據(jù)對比，對比結果如表1。

表1 包頭站數(shù)據(jù)對比及誤差分析

前面通過IDL調(diào)用6S輻射模型得到氣溶膠反演的查找表，查找表中所含的AOD值只有0.0001、0.25幾個值，即當真實AOD值小于0.0001時反演結果為0，當AOD值小于0.25是反演結果為0.0001。故上述數(shù)據(jù)均在實際值內(nèi)，這表明反演結果具有可靠性。

3 結語

本文采用暗像元法對呼包地區(qū)AOT進行了反演，通過反演結果結合地面觀測數(shù)據(jù)可得出以下幾點：

（1）通過地面觀測數(shù)據(jù)驗證遙感反演的AOT精度，結果表明MODIS數(shù)據(jù)反演所得AOT在呼包地區(qū)比較符合實際，并且產(chǎn)生的誤差在允許范圍內(nèi)，故此可利用遙感反演AOT來分析其在呼包地區(qū)空間分布及變化規(guī)律。

（2）由圖1—圖2可以看出同時段大部分區(qū)域AOT值小于0.5，說明呼包地區(qū)大部分地方的AOT的值較小，城區(qū)AOT值明顯高于郊區(qū)，由此說明氣溶膠與人口密集程度、人類建設活動密切相關。

（3）如圖2所示，7月15日AOT平均值明顯低于其它時間AOT值，這說明夏季雨水對氣溶膠具有一定的吸附作用。

本文驗證了利用MODIS數(shù)據(jù)反演AOT來監(jiān)測呼包地區(qū)環(huán)境變化趨勢較為可行，反演結果可信度較高，能在一定程度上為環(huán)保部門提供決策依據(jù)。

[1]Aiimoto H.Global Air Quality and Poll-ution[J].Science，2003（302）：1716-1719.

[2]王耀庭,王橋,李炬,竇軍霞.氣溶膠光學厚度在大氣污染研究中的應用[A].中國氣象學會2007年年會,2007.

[3]劉兆亮,麻金繼.基于MODIS數(shù)據(jù)反演高反射率地區(qū)氣溶膠光學厚度的方法研究[J].大氣與環(huán)境光學學報,2012.