基于克里格插值法的MODIS全球氣溶膠光學厚度研究

魏云霞

(西安科技大學 陜西 西安 710000)

?

基于克里格插值法的MODIS全球氣溶膠光學厚度研究

魏云霞

(西安科技大學 陜西 西安 710000)

基于MODIS全球氣溶膠光學厚度AOD大氣產品數據，利用IDL8.5編程語言基于變程差值的思想對全球氣溶膠光學厚度數據進行克里格插值處理，然后利用地基太陽輻射計觀測網站點實測數據(SONET/AERONET)對插值后的結果中的中國區域進行驗證分析，得出了基于克里格插值方法的全球氣溶膠光學厚度分析方法，該方法對研究大尺度觀測范圍的氣溶膠光學厚度具有一定的參考意義，并可以通過此方法來推測出未知區域的氣溶膠光學厚度(AOD)值，從而彌補站點數據缺無法實現大范圍、宏觀的檢測不足。

氣溶膠光學厚度；克里格插值；站點實測數據

一、引言

大氣氣溶膠是指以大氣與懸浮在其中的固體和液體共同組成的多相體系，粒子的直徑一般在0.001μm～10μm之間，主要來源于自然產生和人類活動，如海鹽微粒、陸地塵土、煙塵等。氣溶膠對氣候、人類健康、遙感定量研究等都有不同程度的影響。2002年我國有統計數據的343個城市中，由于空氣污染導致死亡的人數達到5萬多人，有近40萬人因為空氣污染導致呼吸系統疾病。因此對氣溶膠各種性質的研究具有重要的意義。

氣溶膠光學厚度(Aerosol Optical Depth,AOD)是大氣氣溶膠的重要光學特征之一，它是指沿輻射傳輸路徑，單位截面上所有吸收和散射氣溶膠產生的總削弱，是無云大氣條件下，單位鉛直氣柱中氣溶膠散射造成的消光系數在垂直方向上的積分(無量綱量)，反映氣溶膠對光的衰減作用，是推算氣溶膠含量、評估大氣污染程度、研究氣溶膠氣候效應的關鍵因子，因而是研究大氣氣溶膠含量的最常用指標。

國內外學者目前對于氣溶膠的研究主要是地基遙感和衛星遙感兩種方式，本文是通過全球MODIS氣溶膠光學厚度數據，采用地統計學中的克里格插值方法反演全球范圍內的氣溶膠光學厚度，然后就中國區域與真實的站點數據的AOD進行比較，探究出基于克里格插值法的氣溶膠光學厚度分析方法對于研究氣溶膠光學厚度具有一定的意義，因此在以后研究氣溶膠光學厚度時，可以通過若干個已知點的氣溶膠光學厚度值，利用克里格插值法可以求出未知區域的氣溶膠光學厚度。

二、數據處理

數據的獲得到使用分別經歷了MODIS數據下載，數據拼接，數據重采樣，本文所采用的MODIS數據是通過NASA的官網Http://ladsweb.nascom.nasa.gov/，下載的2015年全年的全球MODIS氣溶膠光學厚度數據，下載的數據為HDF格式。然后對下載的數據進行拼接和重采樣，最后數據的像元分辨率重采樣為0.05°，3600列*7200行。

三、克里格插值

由于云覆蓋或者茂密植被像元缺失等原因，MODIS氣溶膠產品含有大量的無效值像元，本文基于地統計學中克里格插值方法，對全球大氣氣溶膠光學厚度產品進行插值處理，采用變程插值思想，利用ENVI/IDL編程語言來實現。具體思路如圖1

圖1 克里格插值流程圖

四、氣溶膠光學厚度驗證

(一)定義投影信息

由于經過克里格插值后的氣溶膠光學厚度時沒有地理坐標信息，需要為其添加投影信息，最后2015年所有的氣溶膠光學厚度產品統一添加WGS84經緯度坐標。

(二)地基太陽輻射計觀測網站點實測數據

太陽/天空輻射計(Sun/sky-radiometer)是測量太陽直射輻射信號和天空多角度散射輻射信號的遙感儀器，根據地面觀測的太陽直射及散射信號能推斷出氣溶膠的光學特性、成分含量等信息。目前，通過地基太陽輻射計觀測網站點觀測獲取的大氣參數精度最高，參數信息最全面，是驗證其他大氣產品，特別是大氣遙感衛星獲取的氣溶膠、水汽等的基準。

AERONET(Aerosol Robotic Network)是由美國宇航局(NASA)聯合全球多個區域太陽輻射計觀測網建立的地基氣溶膠監測網，在全球范圍內擁有500多個站點，是目前全球分布最廣泛的觀測網。旨在研究全球范圍內的氣溶膠傳輸、輻射效應，監測氣溶膠特性以及驗證衛星遙感氣溶膠產品。AERONET觀測網統一采用法國CIMEL公司生產的全自動CE318太陽光度計為觀測儀器，利用不同的站點的實測數據可反演氣溶膠光學厚度、柱水汽含量等大氣參數產品。AERONET數據主要分為三個處理等級：1.0級數據為未經處理的原始數據；1.5級數據為經過云掩膜處理的數據；2.0級數據是做過云掩膜處理且經過人工審查的高質量數據。本文采用AERONET 2.0級氣溶膠光學厚度產品對MODIS氣溶膠產品進行驗證。

(三)精度驗證

由于MODIS氣溶膠產品提供470、550、660 nm 3個波段處的氣溶膠光學厚度，而地基太陽輻射計觀測網AERONET/SONET中使用的CE318太陽光度計觀測的太陽輻射可以反演340、380、440、500、670、870、1020 nm 7個波段的氣溶膠光學厚度，兩者沒有相對應的波段；且地基觀測網的氣溶膠產品是空間上某些點固定間隔時間的連續觀測，而MODIS大氣產品是空間上以10 km×10 km為觀測單位的面上瞬時觀測，兩者具有不同的時空尺度；因此為了利用地基觀測網的實測數據對MODIS大氣產品進行比較，需要對地基觀測網實測氣溶膠產品進行波長插值。

1.波長插值

1929年，?ngstr?m提出了?ngstr?m指數，建立了氣溶膠光學厚度與波長之間的?ngstr?m關系式，即在沒有水汽影響的波段上,氣溶膠粒子的譜分布滿足Junge分布：

τa(λ)=βλ-α

3-11

其中，τa(λ)表示波長為λ的氣溶膠光學厚度，β表示?ngstr?m渾濁度系數，與氣溶膠粒子總數和譜分布等有關；α表示?ngstr?m波長指數，其數值與氣溶膠平均半徑有關，氣溶膠粒子越大，α越小，參數α和β可以通過2個不同波長處的氣溶膠光學厚度計算得出。本文基于500和670nm處氣溶膠光學厚度利用式3-12、3-13插值計算550nm處的氣溶膠光學厚度：

3-12

3-13

2.時空匹配

如前所述，AERONET地基觀測的氣溶膠光學厚度是以觀測站點為觀測單位的點上連續觀測，而MODIS氣溶膠光學厚度是以10km×10km為觀測單元的面上瞬時觀測。為了使兩者具有可比性，本文采用NASA檢驗全球陸地上空MODIS氣溶膠產品的辦法，即選取以站點為中心50km×50km范圍內的MODIS氣溶膠產品的空間平均與衛星過境前后30min的氣象站點實測數據時間平均進行匹配。

按照上述方法選取MODIS大氣產品和地基實測數據進行相關性分析，驗證MODIS氣溶膠產品在中國區域的適用性。收集全球范圍內的AERONET/SONET氣象站點數據，經過波長插值、時空匹配以及MOIDS氣溶膠產品無效值去除后，獲取時間和空間相匹配的氣溶膠樣然后對樣本進行線性擬合分析來驗證全球氣溶膠產品

τMODIS=aτAERONET/SONET+b

其中τMODIS為MODIS大氣產品，τAERONET/SONET為AERONET實測的大氣產品

圖中實線表示兩者擬合線，虛線表示1:1線，從散點圖中可看出，MODIS大氣產品與氣象站點實測數據擬合線接近于1:1線，兩者具有良好的統計相關關系，氣溶膠產品的相關系數R2達到0.77，RMSE為0.125；水汽產品相關系數R2達到0.98，RMSE為0.26；。除個別站點(如民勤)外，MODIS大氣產品在中國區域與氣象站點實測數據一致性較好。

四、結語

本文研究利用了MODIS的氣溶膠光學厚度AOD值數據，基于地統計學中的克里格插值的方法對全球氣溶膠光學厚度AOD的值進行了反演，對反演的結果的中國區域與真實的站點數據進行了擬合分析，結果表明，插值后MODIS氣溶膠光學厚度在中國區域與氣象站點實測數據一致性較好.，從而可以證實基于克里格插值法的氣溶膠光學厚度分析法對于研究氣溶膠具有一定的意義，該方法可用于彌補MODIS數據由于云覆蓋等原因引起的AOD值的缺失的區域，反演出的的AOD值可以代替空缺地區的AOD值進行研究分析和大范圍大氣監測。

魏云霞(1993-)，女，漢族，山西大同人，理學碩士，西安科技大學測繪科學與技術學院，研究方向：地圖學與地理信息系統。