国产日韩欧美一区二区三区三州_亚洲少妇熟女av_久久久久亚洲av国产精品_波多野结衣网站一区二区_亚洲欧美色片在线91_国产亚洲精品精品国产优播av_日本一区二区三区波多野结衣 _久久国产av不卡

?

MODIS C5、C6氣溶膠產(chǎn)品驗證及區(qū)域適應性評價

2016-12-22 02:55:53楊以坤田信鵬
關鍵詞:氣溶膠反射率光學

楊以坤,孫 林,韋 晶,田信鵬,賈 臣

(山東科技大學 測繪科學與工程學院,山東 青島 266590)

?

MODIS C5、C6氣溶膠產(chǎn)品驗證及區(qū)域適應性評價

楊以坤,孫 林,韋 晶,田信鵬,賈 臣

(山東科技大學 測繪科學與工程學院,山東 青島 266590)

MODIS Collection 5(C5)、Collection 6(C6)氣溶膠產(chǎn)品是當前NASA提供的兩種逐日氣溶膠光學厚度(Aerosol Optical Depth,AOD)產(chǎn)品,可以為大氣顆粒物污染評價以及氣候效應影響研究提供重要的數(shù)據(jù)支撐,但兩種產(chǎn)品分別基于不同的氣溶膠反演方法生成,明確兩種產(chǎn)品氣溶膠反演方法的精度和區(qū)域適應性對用戶合理選擇數(shù)據(jù)產(chǎn)品具有重要意義。以分布于我國的AERONET(Aerosol Robotic Network,AERONET)數(shù)據(jù)為依據(jù),分析了兩種產(chǎn)品在幾種典型區(qū)域的精度以及區(qū)域適應性。結果表明:①C6 DT算法整體反演精度略優(yōu)于C5 DT算法;②與DT算法相比DB具有較高的整體反演精度(R=0.91,RMSE~0.166,MAE~0.116),氣溶膠有效觀測天數(shù)明顯增加,能夠降低15%左右的氣溶膠高估情況;③C6 DT算法在濃密植被地區(qū)整體精度較高,但在亮地表區(qū)域適應能力較差,產(chǎn)品缺失值嚴重;而DB算法均能實現(xiàn)暗地表和亮地表區(qū)域上的氣溶膠反演,尤其在城市、稀疏植被覆蓋氣溶膠反演精度整體較高,空間連續(xù)性得到明顯改善。

MODIS;氣溶膠光學厚度;AERONET;暗目標法;深藍算法

氣溶膠光學厚度(Aerosol Optical Depth,AOD)是消光系數(shù)在垂直方向上對整層大氣的積分,也是氣溶膠光學特性的一項重要參數(shù),反映了氣溶膠對光的衰減作用,是定量描述氣溶膠氣候輻射強迫的重要參量[1]。目前大氣氣溶膠測量方法主要包括地基監(jiān)測和衛(wèi)星遙感反演技術。地基監(jiān)測能夠實時獲取空間某點在不同波長處的詳細氣溶膠光學厚度信息,但由于觀測設備、觀測條件的限制,無法獲取廣域尺度范圍內(nèi)的氣溶膠信息,難以從宏觀上反映氣溶膠的分布狀況和變化趨勢[2]。隨著衛(wèi)星遙感技術的快速發(fā)展,氣溶膠遙感反演可以彌補氣溶膠地基監(jiān)測手段在空間觀測尺度的不足,能夠獲得大范圍的氣溶膠數(shù)據(jù)。

國內(nèi)外學者針對大氣氣溶膠衛(wèi)星遙感反演做了大量研究,逐漸發(fā)展成較為成熟的氣溶膠光學厚度反演算法,主要包括暗目標算法[3-5](Dark Target,DT)、結構函數(shù)法[6]、深藍算法[7-9](Deep Blue,DB)等。這些氣溶膠反演算法已被廣泛應用到NOAA/AVHRR[10]、VIIRS[11]、MODIS[4-5、12-13]、Landsat 8[14-15]、HJ衛(wèi)星[16-17]、FY衛(wèi)星[18]等多種傳感器,并以較高精度實現(xiàn)了不同區(qū)域上空的氣溶膠光學厚度反演。其中DT和DB算法已經(jīng)作為MODIS氣溶膠產(chǎn)品的主要反演算法,并提供了全球逐日的氣溶膠觀測產(chǎn)品(MOD04)。到目前為止,MODIS氣溶膠產(chǎn)品經(jīng)歷了多次版本更新,目前,官網(wǎng)提供了MOD04 C5和C6兩種氣溶膠產(chǎn)品版本,其中DT和DB算法為兩種產(chǎn)品的核心氣溶膠反演算法。定量分析和評價不同MODIS氣溶膠在不同地區(qū)的適用性,對用戶合理選擇和使用產(chǎn)品數(shù)據(jù)進行氣溶膠研究具有重要意義。

針對MODIS氣溶膠產(chǎn)品的驗證以及區(qū)域適應性分析,已開展大量研究工作[19]。夏祥鰲[20]運用全球AERONET地基觀測數(shù)據(jù)對MODIS氣溶膠產(chǎn)品進行驗證分析,發(fā)現(xiàn)MODIS氣溶膠產(chǎn)品在陸地上整體高于實測值;王莉莉等[21]運用太陽分光觀測網(wǎng)(CSHNET)氣溶膠光學厚度地基聯(lián)網(wǎng)數(shù)據(jù),評估MODIS氣溶膠光學厚度產(chǎn)品數(shù)據(jù)在中國不同生態(tài)類型和地理區(qū)域的適用性,發(fā)現(xiàn)在地表類型均一,植被覆蓋良好的區(qū)域具有較高的反演精度。周春艷等[22]對比分析了MOD04 C4和C5兩種產(chǎn)品及其在中國北方地區(qū)的適用性,發(fā)現(xiàn)MOD04 C4和C5產(chǎn)品在中國北方的高反射地區(qū)反演精度較差,而C5產(chǎn)品在濃密植被地區(qū)具有較高的反演精度,適用性較高;Xia等[23]利用氣溶膠地基觀測數(shù)據(jù)對我國西北干旱半干旱地區(qū)的MOD04 C4、C5和DB產(chǎn)品進行驗證,結果表明C5產(chǎn)品在反演精度上較C4產(chǎn)品有較大改善。張曉雷等[24]利用AERONET地基觀測數(shù)據(jù)對MODIS C6 3 km數(shù)據(jù)進行驗證,結果表明MODIS 3 km氣溶膠產(chǎn)品質量穩(wěn)定,整體反演精度較高,但氣溶膠高估現(xiàn)象明顯。

為了定量分析和評價當前MODIS氣溶膠產(chǎn)品的精度及其在不同區(qū)域的適應性,本文選擇高、中、低地表反射率多層次分布,地表結構從簡單到復雜的多種地表類型區(qū)域作為研究區(qū),分別選擇國內(nèi)典型的植被覆蓋區(qū)、城市高亮地區(qū)和西北荒漠地區(qū),運用AERONET地基觀測數(shù)據(jù)對MOD04 C5和C6兩種氣溶膠產(chǎn)品進行驗證對比,探討和確定MOD04兩種典型產(chǎn)品在國內(nèi)不同地區(qū)的適應性。

1 MODIS氣溶膠反演算法

MODIS是搭載在Terra和Aqua衛(wèi)星上的一個重要傳感器,分別于1999年12月18日和2002年5月4日發(fā)射成功。具有36個波段,能夠實現(xiàn)0.4~14 μm全光譜范圍覆蓋,空間分辨率包括250 m、500 m和1 km,每兩天可實現(xiàn)一次全球覆蓋。由于其較高的空間分辨率、重訪周期短的優(yōu)勢,MODIS被廣泛應用到全球地表、生物圈、大氣和海洋等研究中。MODIS提供陸地、大氣和海洋等44種標準數(shù)據(jù)產(chǎn)品,其中MOD04為NASA發(fā)布的MODIS二級大氣氣溶膠產(chǎn)品,該產(chǎn)品首次給出全球尺度逐日的氣溶膠光學厚度數(shù)據(jù),在全球范圍內(nèi)得到了廣泛應用[19]。MOD04產(chǎn)品能夠提供470、550、650 nm不同波長逐日的氣溶膠觀測數(shù)據(jù),目前可用的C5版本空間分辨率為10 km,C6版本空間分辨率為3 km和10 km,反演算法包括DT算法和DB算法。

1.1 MODIS DT陸地氣溶膠反演算法

在低地表反射率情況下,衛(wèi)星測得的反射率隨大氣氣溶膠光學厚度的增加而增加,這是暗目標法氣溶膠光學厚度反演的基本原理[25]。最初MODIS地表反射率的確定是基于Kaufman等[4]提出的假設,即在暗目標區(qū)域,近紅外波段的地表反射率與紅波段和藍波段地表反射率之間存在著固定的線性關系即:0.47 μm的地表反射率為2.1 μm地表反射率的1/4,0.66 μm的地表反射率為2.1 μm地表反射率的1/2。MODIS C5與C6 DT AOD反演時地表反射率的確定采用的新思想:在暗目標地區(qū),波段與近紅外波段地表反射率的關系是散射角(Θ)和植被指數(shù)的函數(shù),紅波段與藍波段地表反射率有固定的線性關系,從而計算出紅波段和藍波段的地表反射率,詳細氣溶膠反演算法參見文獻[13]。

與MODIS C5 DT算法相比,C6 DT反演算法在重塵(AOD?1.0)時對云掩膜算法、氣溶膠類型和產(chǎn)品質量控制等方面做了較多改進。例如,在重塵時云掩膜算法利用0.47 μm(500 m分辨率)和1.38 μm(1 km分辨率)兩個波段進行。具體如下:

(1)

其中:σ*、σ分別為3×3像元的表觀反射率絕對標準偏差和平均加權標準差。

1.2 MODIS DB陸地氣溶膠反演算法

暗目標氣溶膠反演算法只能夠反演水體和濃密植被等地表反射率較低地區(qū)的氣溶膠光學厚度,對于城市、沙漠等高地表反射率地區(qū),該方法無能為力。Hsu等[7-8]發(fā)現(xiàn)即使在城市、沙漠等高地表反射率地區(qū),MODIS深藍波段的地表反射率仍然較低,如果地表反射率能夠準確確定,即可實現(xiàn)高亮地區(qū)氣溶膠光學厚度反演。基于這一發(fā)現(xiàn),Hsu等提出一種適用于高亮地表氣溶膠反演的深藍算法,運用MODIS數(shù)據(jù)成功對撒哈拉沙漠地區(qū)進行氣溶膠反演,達到較高的反演精度。

DB算法基于大部分地物的地表反射率在較短時間保持不變這一假設,利用SeaWiFS地表反射率數(shù)據(jù)構建了季節(jié)性的地表反射率數(shù)據(jù)庫,以此提供真實地表反射率數(shù)據(jù)支持氣溶膠反演。由于地表反射率數(shù)據(jù)庫僅是季節(jié)的函數(shù),這種方法在沙漠和半干旱地區(qū)等地表反射率年際變化及受雙向反射分布函數(shù)(BidirectionalReflectanceDistributionFunction,BRDF)影響較小的地區(qū)取得較好結果;但在植被和非植被混合的像元地區(qū)反演精度較低。與MOD04C5DB算法相比,C6DB算法在地表反射率數(shù)據(jù)庫構建、氣溶膠模式選擇和云/雪掩膜算法等方面進行了改進[9]。

C6DB氣溶膠反演算法在構建地表反射率數(shù)據(jù)庫時,首先利用MODIS地表類型產(chǎn)品(MCD12C1)將地表類型劃分為四類:干旱和半干旱區(qū)域、植被區(qū)域、城市及邊緣過渡區(qū)域和冰雪;除冰雪地區(qū)外在不同地表類型區(qū)域利用不同的方法計算地表反射率。

在干旱和半干旱區(qū)域,為更好地獲取觀測統(tǒng)計數(shù)據(jù),MOD04C6DB利用2002—2009年8年的MODIS表觀反射率數(shù)據(jù)輔助地表反射率數(shù)據(jù)庫的建立,利用MODIS三個波段的表觀反射率進行云、云邊緣和厚氣溶膠檢測;同時為減小地表覆蓋類型的潛在變化,地表反射率數(shù)據(jù)庫除考慮季節(jié)變化外,利用NDVI判定季節(jié)性和暫時水體。最后將像元地表反射率值與散射角進行二次多項式擬合確定地表反射率數(shù)據(jù)庫,其中散射角(Φ)定義為:

Φ=cos-1(-cosθ0cosθ+sinθ0sinθcosφ)。

(2)

其中θ0、θ、φ分別為太陽天頂角、衛(wèi)星天頂角和相對方位角。

在植被區(qū)域,根據(jù)可見光和近紅外波段地表反射率的關系說明植物物候變化對地表反射率的影響。利用AERONET站點監(jiān)測的氣溶膠光學厚度和單次散射反照率數(shù)據(jù)支持MODIS影像進行大氣校正得到參照地表反射率。根據(jù)2.1μm通道表觀反射率(R2.1)和地表覆蓋類型估計可見光通道地表反射率,公式如下:

ESR0.65=a+b×R2.1+c×(R2.1)2,

(3)

ESR0.47=d+e×ESR0.65。

(4)

其中,ESR為地表反射率的估計值,系數(shù)a、b、c、d和e是將參考地表反射率和AERONET站點觀測值進行最小二乘擬合得到。

在城市及邊緣過渡區(qū)域,可見光波段地表反射率與2.1μm的地表反射率無線性關系,并且受BRDF影響較大;因此在地表反射率庫構建時,利用混合方法確定地表反射率,即結合深藍地表底層數(shù)據(jù)庫與地表BRDF角度效應。根據(jù)AERONET站點監(jiān)測的氣溶膠厚度(AOT<0.5)對影像進行大氣校正獲得地表反射率數(shù)據(jù),以此來降低氣溶膠光學厚度對地表反射率確定的影響;同時為減小植被變化對地表反射率的影響,根據(jù)NDVI值將地表反射率劃分成3個NDVI組:NDVI≤0.19;0.190.24。然后在不同NDVI組和不同季節(jié)對散射角度和地表反射率進行二次多項式擬合,最后根據(jù)散射角度和深藍地表底層數(shù)據(jù)庫得到地表反射率數(shù)據(jù)。

C5 云/雪掩膜算法根據(jù)0.412μm波段3×3像元表觀反射率的變化,利用閾值進行云/雪掩膜。在地表反射率變化較大的地區(qū),為充分進行云像元的掩膜,C6云/雪掩膜算法除利用0.412μm波段表觀反射率外,還利用11μm波段的亮度溫度以及11μm波段和12μm波段亮度溫度之差的關系進行云檢測;綜合使用熱紅外波段和藍光波段能夠更清楚地辨別云和下墊面特別是高亮的非均一地表。C6算法利用1.38μm表觀反射率、11μm亮度溫度、11,12μm波段亮度溫度之差和大氣可降水量作為輔助數(shù)據(jù)進行薄卷云檢測。

總體來說,MODISC6DB氣溶膠光學厚度產(chǎn)品,首先利用1km分辨率影像進行氣溶膠光學厚度反演,然后再將其合成10km分辨率的氣溶膠光學厚度產(chǎn)品,并且改進了云掩模算法和地表反射率數(shù)據(jù)庫,更新了塵埃顆粒識別的氣溶膠模型方案,并修訂了產(chǎn)品質控[26]。

2 數(shù)據(jù)源及評價方法

2.1 AERONET數(shù)據(jù)及處理

2.1.1 數(shù)據(jù)介紹

AERONET(AerosolRoboticNetwork)是由NASA和LOA-PHOTONS(CNRS)聯(lián)合建立的全球地基氣溶膠遙感觀測網(wǎng)。目前,AERONET已經(jīng)覆蓋全球主要區(qū)域,該數(shù)據(jù)利用CIMEL自動太陽光度計(SPAM)觀測得到,可以提供380,440,500,675,870,1020nm波長處的氣溶膠光學厚度,每15min提供一次觀測數(shù)據(jù)[27]。數(shù)據(jù)整體精度較高,誤差約0.01~0.02。AERONET對于研究全球氣溶膠的傳輸、氣溶膠輻射效應、驗證輻射傳輸模式以及校驗衛(wèi)星遙感氣溶膠的結果具有重要作用[28]。AERONET提供了3個質量等級的氣溶膠光學厚度數(shù)據(jù):Level1.0,未經(jīng)過濾云和最后驗證的數(shù)據(jù);Level1.5,經(jīng)過嚴格濾云但沒有最后驗證的數(shù)據(jù);Level2.0,經(jīng)過嚴格濾云和質量控制的數(shù)據(jù)。本文運用的AERONET觀測數(shù)據(jù)詳細信息如表1所示。

表1AERONET站點詳細信息

Tab.1ThedetailedinformationofAERONETstation

站點名稱經(jīng)度/(°E)緯度/(°N)高程/m級別時間地表類型Beijing116.381439.976992.002.02013城市Beijing_RADI116.378640.004859.001.52013城市Beijing-CAMS116.316739.9333106.02.02013城市XiangHe116.961539.753636.002.02013植被SACOL104.137135.94601965.802.02010植被NAM_CO90.962530.77254740.002.02010—2014荒漠

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

MODIS提供了470、550、660nm波長處的氣溶膠光學厚度,與AERONET地面觀測數(shù)據(jù)沒有對應的波段。為了實現(xiàn)運用AERONET數(shù)據(jù)對MODIS氣溶膠反演產(chǎn)品進行驗證對比,首先需要對兩者進行波長匹配,本文運用的是MODIS550nm處的氣溶膠產(chǎn)品;然后利用?ngstr?m指數(shù)對AERONET數(shù)據(jù)進行插值,得到550nm處的氣溶膠實測值。在無水汽影響的波段上,氣溶膠離子的譜分布滿足Junge分布,氣溶膠光學厚度與波長的定量關系滿足下式[29]:

τa(λ)=βλ-α。

(5)

其中,τa(λ)表示波長λ的氣溶膠光學厚度;β表示?ngstr?m渾濁度系數(shù),與氣溶膠粒子總數(shù)、粒子譜分布和折射指數(shù)有關;α表示?ngstr?m波長指數(shù),與氣溶膠的平均半徑有關。利用AERONET地基監(jiān)測站點提供的0.44,0.5,0.675μm處的氣溶膠光學厚度帶入式(5)求出β和α,然后將λ=0.55μm帶入求得0.55μm處的氣溶膠光學厚度。

MOD04氣溶膠產(chǎn)品的氣溶膠值表示空間上10×10km2的瞬時氣溶膠光學厚度,而AERONET氣溶膠光學厚度是每隔15min的瞬時觀測值;MODIS衛(wèi)星在中國地區(qū)的過境時間約為當?shù)貢r間10:30,MODIS與AERONET氣溶膠監(jiān)測的空間和時間尺度均不相同。如果取衛(wèi)星過境時間的AERONET觀測值與單個像元(10×10km)衛(wèi)星監(jiān)測氣溶膠光學厚度進行對比可信度較低,并且容易出現(xiàn)奇異值現(xiàn)象。為了提高二者對比結果的可信度,需要將衛(wèi)星和地面氣溶膠監(jiān)測數(shù)據(jù)進行時空匹配。本文采用以地基觀測站點為中心,根據(jù)周圍3×3像元(10km分辨率)標準差去除異常值后取平均值,然后與衛(wèi)星過境前后30min內(nèi)的地面觀測的氣溶膠光學厚度平均值進行對比。

2.2 評價方法

利用最小二乘原理對數(shù)據(jù)進行線性擬合估計數(shù)據(jù)集擬合的斜率和截距,同時計算均方根誤差(Root-Mean-Square-Error,RMSE)、平均絕對誤差(Mean-Absolute-Error,MAE)、MODIS氣溶膠產(chǎn)品在陸地上的期望誤差(Expected-Error,EE)和相對平均偏差(Relative-Mean-Bias,RMB)等指標,定量評價和分析氣溶膠的反演精度及不確定性,其中RMB>1.0和RMB<1.0時分別表示氣溶膠光學厚度反演值整體呈現(xiàn)高估和低估現(xiàn)象。

(6)

(7)

EE=±(0.05+0.2AODAERONET);

(8)

(9)

3 結果與分析

3.1 MOD04 C5 DT與C6 DT氣溶膠產(chǎn)品對比驗證

圖1為MOD04 C5、C6 DT氣溶膠產(chǎn)品的對比驗證圖。其中,虛線表示MODIS氣溶膠產(chǎn)品在陸地上的期望誤差;黑色實線表示Y=X;長劃線為線性回歸擬合曲線。由圖1可以看出,MODIS DT算法與AERONET地基觀測數(shù)據(jù)整體具有較高的一致性(R>0.80),其中C6 DT算法的有效觀測點對整體少于C5 DT算法,原因是C6 DT算法在云掩膜、數(shù)據(jù)質量控制上進行一定改進,減少部分異常值。MOD04 C5 DT有50%的觀測點對滿足MODIS氣溶膠反演期望誤差,而C6 DT有53.1%的觀測點對滿足期望誤差??梢哉f明,與C5 DT算法相比(R=0.814,RMSE=0.205,MAE=0.141),C6 DT算法具有更高的氣溶膠整體反演精度,其反演值與地基觀測數(shù)據(jù)具有更高的一致性(R=0.881),兩者具有較小的均方根誤差(RMSE=0.173)和平均絕對誤差(MAE=0.131)。尤其在氣溶膠值較小時(AOD<0.5),C6 DT算法的氣溶膠反演值與實測值間的吻合度整體高于C5 DT算法。

3.2 MOD04 C6 DT與C6 DB氣溶膠產(chǎn)品對比驗證

圖2(a)、2(b)分別為MOD04 C6 DT和DB氣溶膠產(chǎn)品的對比驗證圖。從圖2中可以看出,與MOD04 C6 DT算法相比(R=0.881,RMSE=0.173,MAE=0.131,RMB=1.157),DB算法的氣溶膠反演值與AERONET地基觀測數(shù)據(jù)具有更高的一致性,兩者相關系數(shù)達到0.91,均方根誤差(RMSE~0.166)和相對誤差(MAE~0.116)整體較小,有61.9%的觀測點對滿足MODIS氣溶膠產(chǎn)品誤差精度要求,能夠有效減少15%左右的氣溶膠反演高估現(xiàn)象(RMB=1.001)。其次,由于DB算法能夠實現(xiàn)暗地表和亮地表區(qū)域氣溶膠反演,DB算法的有效觀測點數(shù)(N=821)明顯高于DT算法(N=371),且氣溶膠反演值域范圍更廣(0~2.5)。當氣溶膠光學厚度較小(AOD<0.5)時,DB算法的氣溶膠反演值與地基觀測數(shù)據(jù)具有更高的吻合度;同時,DB算法能夠實現(xiàn)高污染天氣(AOD>1.5)下的氣溶膠光學厚度反演,且反演精度整體較高。因此,可以說明C6 DB算法的氣溶膠反演精度整體高于C6 DT算法。

圖1 MODIS DT C5(a)與C6 (b) AOD有效性驗證(分辨率10 km)

Fig.1 Validation of MODIS DT C5(a) and C6 (b) AOD (10 km resolution)

圖2 MODIS DT C6(a)與DB C6 (b) AOD有效性驗證(分辨率10 km)

Fig.2 Validation of MODIS DT C6(a) and DB C6 (b) AOD(10 km resolution)

3.3 MOD04 C6 DT與C6 DB氣溶膠產(chǎn)品區(qū)域適應性評價

MODIS DT與DB氣溶膠反演算法在地表反射率確定方法上存在較大差異。在城市地區(qū)BRDF和混合地表類型對地表反射率的確定影響明顯;在植被茂密地區(qū),由于植被隨季節(jié)的生長規(guī)律導致地表反射率有較大變化,不同氣溶膠反演算法在不同區(qū)域的適應性有一定差異。因此,本文分別選擇濃密植被地區(qū)、氣溶膠反演困難的城市地區(qū)以及荒漠等稀疏植被地區(qū),驗證和對比MOD04 C6兩種不同氣溶膠產(chǎn)品在不同區(qū)域上的適應性。

3.3.1 城市地區(qū)

圖3是MOD04 C6 DT與DB氣溶膠產(chǎn)品在典型城市地區(qū)對比驗證圖,Beijing、Beijing_RADI和Beijing-CAMS站點位于中國北京市,地表反射率較高,屬于典型的城市高亮地表地區(qū)。從圖3可以看出,MOD04 C6 DT算法在城市地區(qū)有效反演天數(shù)整體較少,氣溶膠反演值與AERONET氣溶膠地基觀測數(shù)據(jù)的一致性整體較差(R≈0.734~0.808),具有較大均方根誤差(RMSE≈0.230~0.273)和相對誤差(MAE≈0.171~0.213),僅有41%~51%的有效觀測點數(shù)滿足MODIS氣溶膠產(chǎn)品的誤差精度要求。與DT算法相比,MOD04 C6 DB算法有效反演天數(shù)明顯增加,氣溶膠反演值與AERONET地基觀測數(shù)據(jù)整體具有更高的一致性(R≈0.890~0.924),均方根誤差(RMSE≈0.184~0.197)和平均誤差(MAE≈0.135~0.142)整體較小,有55%~61%的觀測點對滿足誤差精度要求。其次,DB算法氣溶膠反演值域范圍更廣(0~2.5),能夠實現(xiàn)高污染天氣(AOD>1.5)下城市地區(qū)的氣溶膠光學厚度反演,且反演精度整體較高??梢哉f明,城市地區(qū)DB算法比DT算法具有更高的氣溶膠整體反演精度,穩(wěn)定性較高。

圖3 城市地區(qū)MODIS DT C6(a,c,e)和DB C6 (b,d,f) AOD有效性驗證(分辨率10 km)

Fig.3 Validation of MODIS DT C6(a,c,e) and DB C6 (b,d,f) in urban area(10 km resolution)

3.3.2 植被地區(qū)

圖4是MOD04 C6 DT與DB氣溶膠產(chǎn)品在植被地區(qū)的對比驗證圖。XiangHe和SACOL站分位于中國東部地區(qū),屬于典型的植被覆蓋地區(qū),由圖4可以看出,MODIS DT和DB算法均能以較高精度實現(xiàn)氣溶膠光學厚度反演,與DB算法相比,DT算法氣溶膠反演值與實測值存在更高的一致性(R≈0.936~0.938),具有更小的均方根誤差(RMSE≈0.050~0.152)和相對誤差(MAE≈0.042~0.110),有60%~90%的觀測點對滿足誤差精度要求。當氣溶膠光學厚度較小(AOD<0.5)時,DT算法的氣溶膠反演值與地基觀測數(shù)據(jù)具有更高的吻合度,與AERONET實測值擬合曲線更接近于Y=X直線;因此,可以說明DT算法在植被覆蓋地區(qū)較DB算法具有更高的整體反演精度。

3.3.3 稀疏植被及荒漠地區(qū)

圖5為MODIS C6 DT與DB氣溶膠產(chǎn)品在稀疏植被及荒漠地區(qū)的對比驗證結果。NAM_CO站位于中國內(nèi)陸西藏高原上,海拔較高,植被稀少,地表裸露,空氣污染較輕,氣溶膠光學厚度值整體較小。由圖5可以看出,DT算法在該地區(qū)的氣溶膠有效觀測點數(shù)較少(N=80),氣溶膠反演精度整體較差(R=0.197),存在較大均方根誤差(RMSE=0.278)和相對誤差(MAE=0.188),僅有28.7%的觀測點對滿足誤差精度要求,有70%的點出現(xiàn)高估現(xiàn)象。原因是該地區(qū)屬于典型的亮地表區(qū)域,植被覆蓋較少,地表反射率的確定較為困難,導致氣溶膠反演精度整體較低,存在明顯的高估現(xiàn)象。DB算法在該地區(qū)的有效觀測天數(shù)增加(N=118),反演值與實測之間的相關性顯著提高(R=0.630),具有較小的均方根誤差和相對誤差,有超過75%的觀測點對滿足誤差精度要求,氣溶膠反演精度整體較高??梢哉f明,在稀疏植被覆蓋等亮地表區(qū)域上,MODIS DB算法能夠以較高精度確定其地表反射率,實現(xiàn)該區(qū)域的氣溶膠光學厚度反演,能夠減少或降低氣溶膠高估現(xiàn)象,彌補DT算法在該地區(qū)氣溶膠反演的不足。

圖4 濃密植被區(qū)MODIS DT C6(a,c)和DB C6 (b,d) AOD有效性驗證(分辨率10 km)

Fig.4 Validation of MODIS DT C6(a,c) and DB C6 (b,d) in vegetation area at 10 km resolution

圖5 植被稀疏地區(qū)MODIS DT C6(a)和DB C6 (b) AOD有效性驗證(分辨率10 km)

Fig.5 Validation of MODIS DT C6(a) and DB C6 (b) in sparse vegetated area(10 km resolution)

3.4 MOD04 C6 DT和C6 DB氣溶膠產(chǎn)品空間分布對比

圖6是不同地區(qū)上空MOD04 C6 DT和DB氣溶膠產(chǎn)品空間分布圖。其中圖6(a)、6(b)分別為北京及周邊地區(qū)等典型城市C6 DT和DB氣溶膠空間分布,時間為2013年10月18日;圖6(c)、6(d)分別是中國西北荒漠及西藏高原等典型荒漠、稀疏植被地區(qū)C6 DT和DB氣溶膠空間分布,時間為2013年9月12日??梢钥闯觯琈OD04 DT算法在城市(圖6(a))、荒漠(圖6(c))等高反射率地區(qū)中氣溶膠反演無能為力,缺失值較為嚴重,空間連續(xù)性差,其氣溶膠產(chǎn)品難以滿足大氣氣溶膠相關研究需要。然而MOD04 DB算法能夠克服DT算法在城市、荒漠等高亮地表區(qū)域氣溶膠反演困難問題,實現(xiàn)該區(qū)上空氣溶膠高精度反演,數(shù)據(jù)缺失程度較小,空間連續(xù)性明顯得到改善。因此,可以說明MOD04 DB算法能夠實現(xiàn)高地表反射率地區(qū)的氣溶膠反演,其氣溶膠產(chǎn)品在反演范圍和空間連續(xù)性上明顯優(yōu)于DT算法。

圖6 MODIS C6 DT AOD與C6 DB AOD空間分布比較(分辨率10 km)

Fig.6 Spatial distribution comparison MODIS C6 DT AOD and C6 DB AOD in space(10 km resolution)

4 結論

本文詳細介紹了當前常用的MODIS DT和DB兩種氣溶膠產(chǎn)品陸地算法,選擇分布在不同下墊面上的AERONET站點氣溶膠地基觀測數(shù)據(jù),對獲取的長時間序列MODIS DT和DB氣溶膠產(chǎn)品進行精度驗證對比,并分析了不同產(chǎn)品在不同區(qū)域上的適應性。結果表明:

1) MOD04 C5 DT算法與C6 DT算法具有相似氣溶膠反演精度,有超過50%的觀測點對滿足MODIS氣溶膠誤差精度要求。整體上,C6 DT略優(yōu)于C5 DT算法,原因是C6 DT算法在云掩膜、質量控制等方面進行了一定的改進。

2) 與MOD04 C6 DT算法相比,C6 DB算法有效觀測天數(shù)明顯增加,氣溶膠反演值域范圍更廣,具有更高的氣溶膠整體反演精度(R=0.91,RMSE~0.166,MAE~0.116),能夠有效降低15%左右的氣溶膠反演高估現(xiàn)象。

3) MOD04 C6 DT氣溶膠產(chǎn)品能夠以較高精度實現(xiàn)植被地區(qū)的氣溶膠反演,但在城市和稀疏植被地區(qū)氣溶膠反演精度整體較差,有效觀測天數(shù)較少;然而DB算法能夠實現(xiàn)暗地表和亮地表區(qū)域上的氣溶膠反演,其中在城市、稀疏植被地區(qū)的氣溶膠反演精度明顯優(yōu)于DT算法,而DT算法在植被地區(qū)的反演精度和穩(wěn)定性略高于DB算法。

4)空間分布上,MOD04 DT氣溶膠產(chǎn)品在城市、稀疏植被等高亮地區(qū)中缺失值較為嚴重,空間連續(xù)性差。相反,DB氣溶膠產(chǎn)品反演范圍和空間連續(xù)性上明顯優(yōu)于DT算法,能夠克服和改善DT算法在城市、荒漠等高亮地表區(qū)域氣溶膠反演困難問題。

致 謝:感謝NASA GSFC提供的MODIS氣溶膠產(chǎn)品數(shù)據(jù)和AERONET氣溶膠地基觀測數(shù)據(jù)。

[1]齊玉磊,葛覲銘,黃建平.北方地區(qū)MODIS和MISR與AERONET氣溶膠光學厚度的比較及其時空分布分析[J].科學通報,2013,17(58):1670-1679. QI Yulei,GE Jianming,HUANG Jianping.Spatial and temporal distribution of MODIS and MISR aerosol optical depth over northern China and comparison with AERONET[J].Chinese Science Bulletin,2013,17(58):1670-1679.

[2]李成才,毛節(jié)泰,劉啟漢,等.MODIS衛(wèi)星遙感氣溶膠產(chǎn)品在北京市大氣污染研究中的應用[J].中國科學:地球科學,2005,35(S1):177-186. LI Chengcai,MAO Jietai,LIU Qihan,et al.Application of MODIS satellite remote sensing aerosol products in the study of air pollution in Beijing city[J].Science China:Earth Sciences,2005,35(S1):177-186.

[3]孫林.城市地區(qū)大氣氣溶膠遙感反演研究[D].北京:中國科學院研究生院(中國科學院遙感應用研究所),2006:4 -11.

[4]KAUFMAN Y J,TANRé D,REMER L A,et al.Operational remote sensing of troposphereic aerosol over land from EOS moderate resolution imaging spectroradiometer[J].Journal of Geophysical Research Atmospheres,1997,102067(D14):17051-17067.

[5]KAUFMAN Y J,WALD A E,REMER L A,et al.The MODIS 2.1 μm channel correlation with visible reflectance for use in remote sensing of aerosol[J].IEEE Transactions on Geoscience & Remote Sensing,1997,35(5):1286-1298.

[6]TANRé D,HOLBEN B N,KAUFMAN Y J.Atmospheric correction algorithm for NOAA-AVHRR products-theory and application[J].IEEE Transactions on Geoscience & Remote Sensing,1992,30(2):231-248.

[7]HSU N C,TSAY S C,KING M D,et al.Aerosol properties over bright-reflecting source regions[J].IEEE Transactions on Geoscience & Remote Sensing,2004,42(3):557-569.

[8]HSU N C,TSAY S C,KING M,et al.Deep blue retrievals of Asian aerosol properties during ACE-Asia[J].IEEE Transactions on Geoscience & Remote Sensing,2006,11(44):3180-3195.

[9]HSU N C,JEONG M J,BETTENHAUSEN C,et al.Enhanced deep blue aerosol retrieval algorithm:The second generation[J].Journal of Geophysical Research Atmospheres,2013,118(16):9296-9315.

[10]BRAUN M.Retrieval of aerosol optical depth over land surfaces from AVHRR data[J].Atmospheric Measurement Techniques Discussions,2014,7(1):2411-2420.

[11]蘇城林,蘇林,陳良富,等.NPP VIIRS數(shù)據(jù)反演氣溶膠光學厚度[J].遙感學報,2015,19(6):978-983. SU Chenglin,SU Lin,CHEN Liangfu,et al.Retrieval of aerosol optical depth using NPP VIIRS data[J].Journal of Remote Sensing,2015,19(6):978-983.

[12]LEVY R C,REMER L A,MATTOO S,et al.A new algorithm for retrieving aerosol properties over land from MODIS spectral reflectance[J].Journal of Geophysical Research,2006,7(1):1-53.

[13]LEVY R C,REMER L A,MATTOO S,et al.Second-generation operational algorithm:Retrieval of aerosol properties over land from inversion of moderate resolution imaging spectroradiometer spectral reflectance[J].Journal of Geophysical Research Atmospheres,2007,112(D13):3710-3711.

[14]SUN L,WEI J,BILAL M,et al.Aerosol optical depth retrieval over bright areas using Landsat 8 OLI images[J].Remote Sensing,2015,8(1):1-14.

[15]孫林,田信鵬,于會泳,等.Landsat8 OLI城市地區(qū)氣溶膠光學厚度遙感反演[J].山東科技大學學報(自然科學版),2015,34(5):47-55. SUN Lin,TIAN Xinpeng,YU Huiyong,et al.Remote sensing inversion of urban aerosol optical depth from on landsat 8 OLI[J].Journal of Shandong University of Science and Technology(Nature Science),2015,34(5):47-55.

[16]SUN K,SUN C K,LIU Q H,et al.Aerosol optical depth retrieval by HJ-1/CCD supported by MODIS surface reflectance data[J].Science China:Earth Sciences,2010,53(S1):74-80.

[17]王中挺,厲青,王橋,等.利用深藍算法從HJ-1數(shù)據(jù)反演陸地氣溶膠[J].遙感學報,2012,16(3):596-610. WANG Zhongting,LI Qing,WANG Qiao,et al.HJ-1 terrestrial aerosol data retrieval using deep blue algorithm[J].Journal of Remote Sensing,2012,16(3):596-610.

[18]周永波,白潔,周著華,等.FY-3A/MERSI海上沙塵天氣氣溶膠光學厚度反演[J].遙感學報,2014,18(4):771-787. ZHOU Yongbo,BAI Jie,ZHOU Zhuhua,et al.Aerosol optical depth retrieval from FY-3A/MERSI for sand-dust weather over ocean[J].Journal of Remote Sensing,2014,18(4):771-787.

[19]ICHOKU C,KAUFMAN Y J,REMER L A,et al.Global aerosol remote sensing from MODIS[J].Advances in Space Research,2004,34(4):820-827.

[20]夏祥鰲.全球陸地上空MODIS氣溶膠光學厚度顯著偏高[J].科學通報,2006,51(19):2297-2303. XIA Xiangao.MODIS aerosol optical depth over the global land significantly higher[J].Chinese Science Bulletin,2006,51(19):2297-2303.

[21]王莉莉,辛金元,王躍思,等.CSHNET觀測網(wǎng)評估MODIS氣溶膠產(chǎn)品在中國區(qū)域的適用性[J].科學通報,2007,52(4):477-486. WANG Lili,XIN Jinyuan,WANG Yuesi,et al.CSHNET observation network to assess the applicability of MODIS aerosol products in China region[J].Chinese Science Bulletin,2007,52(4):477-486.

[22]周春艷,柳欽火,唐勇,等.MODIS氣溶膠C004、C005產(chǎn)品的對比分析及其在中國北方地區(qū)的適用性評價[J].遙感學報,2009,13(5):854-872. ZHOU Chunyan,LIU Qinhuo,TANG Yong,et al.Comparison between MODIS aerosol product C004 and C005 and evaluation of their applicability in the north of China[J].Journal of Remote Sensing,2009,13(5):854-872.

[23]LI X,XIA X A,WANG S L,et al.Validation of MODIS and deep blue aerosol optical depth retrievals in an arid/semi-arid region of northwest China[J].Particuology,2012,10(1):132-139.

[24]孫曉雷,甘偉,林燕,等.MODIS 3 km氣溶膠光學厚度產(chǎn)品檢驗及其環(huán)境空氣質量指示[J].環(huán)境科學學報,2015,35(6):1657-1666. SUN Xiaolei,GAN Wei,LIN Yan,et al.Validation of MODIS 3 km aerosol optical depth product and its air quality indication[J].Acta Scientiae Circumstantiae,2015,35(6):1657-1666.

[25]KAUFMAN Y J,SENDRA C.Algorithm for automatic atmospheric corrections to visible and near-IR satellite imagery[J].International Journal of Remote Sensing,1988,9(8):1357-1381.

[26]MI W,LI Z Q,XIA X G,et al.Evaluation of the moderate resolution imaging spectroradiometer aerosol products at two aerosol robotic network stations in China[J].Journal of Geophysical Research Atmospheres,2007,112(D22):321-341.

[27]BILAL M,NICHOL J E.Evaluation of MODIS aerosol retrieval algorithms over the Beijing-Tianjin-Hebei region during low to very high pollution events[J].Journal of Geophysical Research Atmospheres,2015,120(5):7941-7957.

[28]HOLBEN B N,ECK T F,SLUTSKER I,et al.Aeronet-A federated instrument network and data archive for aerosol characterization[J].Remote Sensing of Environment,1998,66(1):1-16.

[29]?NGSTR?M A.The parameters of atmospheric turbidity[J].Tellus,1964,16(1):64-75.

(責任編輯:高麗華)

Verification of MODIS C5 and C6 and Their Regional Adaptability Evaluation

YANG Yikun,SUN Lin,WEI Jing,TIAN Xinpeng,JIA Chen

(College of Geomatics, Shandong University of Science and Technology,Qingdao,Shandong 266590,China)

MODIS Collection 5(C5) and Collection 6(C6), the two current daily Aerosol Optical Depth (AOD) products based on the Dark Target (DT) and Deep Blue (DB) algorithms provided by NASA,provide important data support for the evaluation of atmospheric particulate matter and the effects of climate change.Because they are based on different aerosol retrieval methods,to make clear the precision and regional adaptability of the two products is of great significance for users in selecting reasonable data products.In this paper,the accuracy and regional adaptability of the two products in several typical regions were analyzed based on the data of AERONET distibuted in China.The results show that: ①The overall retrieval accuracy of DT C6 algorithm is slightly better than that of DT C5 algorithm;②Compared with DT algorithm,DB algorithm has an overall higher retrieval accuracy (R=0.91,RMSE~0.166,MAE~0.116) with significantly increased effective observation days and can reduce approximately 15% aerosol overestimation;③With higher accuracy in vegetation areas but poorer adaptability in bright areas,C6 DT algorithm has large missing values,while DB algorithm,with the capability of retrieving aerosol over both dark and bright surfaces,has higher aerosol accuracy in urban and sparse vegetated areas,thus improving the spatial continuity significantly.

MODIS;aerosol optical depth;AERONET;dark target;deep blue

2016-03-03

國家自然科學基金項目(41171270);山東省杰出青年基金項目(2012JQB01025)

楊以坤(1992—),男,山東菏澤人,碩士研究生,主要從事大氣遙感研究.E-mail:yangykun_rs@163.com 孫 林(1975—),男,山東棗莊人,教授,博士,主要從事大氣探測方面的研究.E-mail:sunlin6@126.com

P407

A

1672-3767(2016)05-0017-11

猜你喜歡
氣溶膠反射率光學
影響Mini LED板油墨層反射率的因素
近岸水體異源遙感反射率產(chǎn)品的融合方法研究
海洋通報(2022年4期)2022-10-10 07:40:26
滑輪組的裝配
具有顏色恒常性的光譜反射率重建
光學常見考題逐個擊破
氣溶膠傳播之謎
氣溶膠中210Po測定的不確定度評定
化學腐蝕硅表面結構反射率影響因素的研究*
電子器件(2017年2期)2017-04-25 08:58:37
四川盆地秋季氣溶膠與云的相關分析
光學遙感壓縮成像技術
天全县| 海淀区| 江安县| 弥勒县| 浠水县| 洪湖市| 长治市| 尼玛县| 象山县| 金塔县| 三台县| 云和县| 阳朔县| 东乡县| 喀什市| 利辛县| 合江县| 磐石市| 济阳县| 永福县| 湛江市| 电白县| 万盛区| 岑巩县| 沈丘县| 靖江市| 治多县| 新津县| 怀远县| 德清县| 昌吉市| 余庆县| 大石桥市| 怀柔区| 卢湾区| 赤水市| 清涧县| 永清县| 化德县| 佳木斯市| 漠河县|