中國氣溶膠光學(xué)厚度時空演變特征分析

薛文博 武衛(wèi)玲 王金南 楊金田 雷 宇 汪藝梅

(環(huán)境保護(hù)部環(huán)境規(guī)劃院，北京 100012)

1 引言

為掌握全國PM2.5污染狀況，我國于2013年開始在74個城市開展了PM2.5的監(jiān)測與數(shù)據(jù)實時發(fā)布工作，但目前我國開展PM2.5監(jiān)測的城市較少，尚難以反映全國PM2.5污染的總體狀況。由于氣溶膠分布具有較高時空變異性，盡管PM2.5地面觀測站點能夠提供時間連續(xù)的觀測數(shù)據(jù)，但不能得到連續(xù)的空間覆蓋信息，難以有效反映氣溶膠污染的時空演變趨勢。在美國、歐洲等發(fā)達(dá)國家，由于擁有較為完善的空氣質(zhì)量監(jiān)測網(wǎng)絡(luò)和大量的PM2.5濃度監(jiān)測站點，利用 AOD產(chǎn)品反演近地面 PM2.5濃度的研究已經(jīng)廣泛開展。因此用NASA中分辨率成像光譜儀(MODIS)觀測的氣溶膠光學(xué)厚度(aerosol optical depth，AOD)來表征PM2.5污染狀況是一種有效的技術(shù)手段，可以彌補(bǔ)地面觀測手段的不足，獲得高分辨率的空間觀測信息。

大氣氣溶膠是指大氣與懸浮在其中的固體和液體微粒共同組成的多相體系，粒子的粒徑多在0.001～10μm之間，作為氣溶膠最基本的光學(xué)特性之一，氣溶膠光學(xué)厚度(AOD)被定義為介質(zhì)的消光系數(shù)在垂直方向上的積分，是描述氣溶膠對光衰減作用的主要因子。毛節(jié)泰等(2002)證明了MODIS氣溶膠產(chǎn)品與地面光度計觀測結(jié)果具有很好的相關(guān)性;Chu等(2003)利用NASA發(fā)布的全球MODIS氣溶膠產(chǎn)品數(shù)據(jù)，研究全球及區(qū)域性大氣污染狀況，證實了用氣溶膠光學(xué)厚度監(jiān)測大氣污染的可行性;李成才等(2005)利用MODIS產(chǎn)品研究了北京市大氣污染狀況;Wang等(2010)考慮大氣層邊界高度、濕度等因素的影響，建立了北京市AOD和PM2.5近地面濃度之間的回歸反演模型。

氣溶膠的消光效應(yīng)主要來源于細(xì)顆粒物的貢獻(xiàn)，粗粒子消光貢獻(xiàn)相對較低，因此通過AOD的變化趨勢可以有效反映我國PM2.5的污染狀況。為了全面了解我國陸地上空氣溶膠的時空分布及變化規(guī)律，本文利用2005-2012年8年MODIS氣溶膠光學(xué)厚度產(chǎn)品，對中國氣溶膠的時空演變特征進(jìn)行了系統(tǒng)分析，著重分析了京津冀、長三角、華中、成渝四個重點區(qū)域及北京、石家莊、保定、邢臺四個典型城市的AOD演變特征，揭示了我國近八年來氣溶膠污染的時空演變過程。

2 數(shù)據(jù)與方法

2.1 數(shù)據(jù)來源

美國國家宇航局(NASA)于1999年12月發(fā)射了Terra衛(wèi)星，于2002年5月發(fā)射了Aqua衛(wèi)星，MODIS分別搭載在Terra和Aqua兩顆衛(wèi)星上。MODIS掃描寬度為2330km，提供了從0.145—14.235μm區(qū)間36個通道的信息，探測器每天覆蓋全球1次，分辨率有250m、500m和1000m三種。本文使用的數(shù)據(jù)來源于EOSAQUA/MODIS傳感器Level 2光學(xué)厚度產(chǎn)品，采用暗像元算算法，空間分辨率可達(dá)到10km。

首先對所獲得的MODIS產(chǎn)品進(jìn)行預(yù)處理，利用IDL對遙感影像進(jìn)行投影、拼接，在ARCMAP軟件中進(jìn)行柵格統(tǒng)計分析，獲取中國陸地2005-2012年550nm的年平均AOD柵格數(shù)據(jù)。

2.2 研究方法

分析我國AOD的時空演變特征，目的是研究我國PM2.5污染的時空變化趨勢，但由于AOD受到空氣中濕度、沙塵、降雨等多重因素的影響，AOD的時空變化趨勢常存在一定的隨機(jī)性，不一定能準(zhǔn)確反映PM2.5污染特征。因此選用長時間序列、多年平均的AOD數(shù)據(jù)在一定程度上可以減少氣象等隨機(jī)性因素的影響。本文統(tǒng)計分析了全國、各省、地級及以上城市2005-2007年、2008-2010年、2011-2012年三個時段的AOD均值及歷年均值，通過對比分析不同時間段、不同地區(qū)的AOD數(shù)值，揭示我國2005-2012年氣溶膠污染的總體特征。

3 AOD時空分析

3.1 全國總體趨勢

從2005-2012年全國AOD時間演變趨勢來看，與2005-2007三年均值相比，絕大部分省份2008-2010年期間AOD下降趨勢非常明顯，全國平均降幅為9.23%。而在2011-2012年期間AOD又有所回升，相比2008-2010年期間全國平均增長了6.12%，但大部分地區(qū)2011-2012兩年平均值依然低于2005-2007三年平均值，因此2005-2007年是我國氣溶膠污染最嚴(yán)重的時期;從AOD的空間分布來看，氣溶膠污染較重的地區(qū)主要集中在京津冀、長三角、華中、成渝等人口密集和工業(yè)發(fā)達(dá)地區(qū)。

京津冀地區(qū)(北京市、天津市、河北省)2005-2007年、2008-2010年、2011-2012年三個時期AOD年均值分別為0.62、0.61、0.64，依次是全國相應(yīng)時期AOD年均值 0.54、0.49、0.52 的 1.2 倍、1.3 倍、1.2 倍;2008-2010年和2011-2012期間，在全國大部分省份AOD下降的趨勢下，北京、天津地區(qū)的AOD保持增長或不變的趨勢，穩(wěn)定在0.8左右的高值區(qū)，氣溶膠污染持續(xù)嚴(yán)重。相比2005-2007年，2008-2010年北京AOD增幅為4.39%，天津增幅為10.57%，2011-2012年北京AOD幾乎無變化，天津增幅為6.53%;相比2008-2010年，北京、天津在2011-2012年間AOD值略有減小。河北省2005-2012年期間AOD總體呈持續(xù)下降趨勢，2008-2010年、2011-2012年分別相比2005-2007年期間下降10.87%、7.03%。

長三角地區(qū)(上海市、江蘇省、浙江省)2005-2007年、2008-2010年、2011-2012年三個時期AOD年平值分別為 0.77、0.72、0.76，是全國相應(yīng)時期 AOD 年均值0.54、0.49、0.52 的 1.43 倍、1.47 倍、1.45 倍，其中上海在各時間段內(nèi)AOD年均值均超過0.9，這可能與臨海濕度或海鹽粒子有很大關(guān)系。相比2005-2007年，上海、江蘇和浙江在2008-2010年期間AOD降幅分別為4.38%、4.41%、10.14%，在2011-2012年期間三個地區(qū)AOD值又有所上升。

華中地區(qū)(湖北省)2005-2007年、2008-2010年、2011-2012年三個時期AOD年均值分別為0.74、0.67、0.7，是全國相應(yīng)時期 AOD年平均值的1.38倍、1.38倍、1.34倍;與2005-2007年相比，湖北省在2008-2010年和2011-2012年總體呈持續(xù)下降趨勢，降幅分別為9.56%、5.68%。

成渝地區(qū)(四川省、重慶市)2005-2007年、2008-2010年、2011-2012年三個時期AOD年均值分別為0.73、0.62、0.65，是全國 AOD 年均值的 1.36 倍、1.27倍、1.25倍。與2005-2007年相比，四川省2008-2010年和2011-2012年期間AOD總體顯著下降，降幅分別達(dá)15.92%和11.07%;重慶市2008-2010年AOD年均值相比2005-2007年期間下降了5.36%，2011-2012年期間又回升到2005-2007年水平。

圖1 AOD年均值分布圖

3.2 典型城市分析

2013年1月份我國發(fā)生了大規(guī)模長周期的灰霾污染事件，華北地區(qū)成為灰霾污染重災(zāi)區(qū)，其中北京、石家莊、保定、邢臺等城市是華北地區(qū)污染最嚴(yán)重的城市群，因此本文重點分析上述四個城市2005-2012年間AOD的變化趨勢。分析結(jié)果表明，與全國333個地級及以上城市AOD年平均值相比，北京、石家莊、保定、邢臺地區(qū)AOD年均值一直在0.8附近波動，遠(yuǎn)高于2005-2012年全國AOD歷年均值。2005-2012期間北京、石家莊、保定、邢臺AOD的8年平均值分別是全國平均水平的1.51倍、1.53倍、1.44倍、1.52倍。這些城市近八年AOD年均值變化趨勢較為一致，具有時間同步性。

2005-2007 年全國地級及以上城市AOD年均值總體持續(xù)增長，到2008-2010年AOD年均值又呈現(xiàn)下降趨勢，在2011-2012年期間全國城市AOD年均值又明顯回升，這種變化趨勢可能與2008-2010年期間全球經(jīng)濟(jì)危機(jī)導(dǎo)致的經(jīng)濟(jì)發(fā)展放緩，能源消費(fèi)量及各種污染物排放量下降有直接關(guān)系。

表1 全國各省AOD年均值變化情況(增幅以2005-2007年為基準(zhǔn))

圖2 全國及典型城市AOD年均值變化

4 研究結(jié)論

由于大氣氣溶膠分布具有較高時空變異性，盡管PM2.5地面觀測站點能夠提供時間連續(xù)的觀測數(shù)據(jù)，但不能得到連續(xù)的空間覆蓋信息，難以有效反映氣溶膠污染的時空演變趨勢。為此，歐美廣泛開展利用AOD產(chǎn)品反演近地面PM2.5濃度的研究。用NASA中分辨率成像光譜儀(MODIS)觀測的氣溶膠光學(xué)厚度(aerosol optical depth，AOD)來表征PM2.5污染狀況是一種有效的技術(shù)手段，可以彌補(bǔ)地面觀測手段的不足，獲得高分辨率的空間觀測信息。

大氣氣溶膠是指大氣與懸浮在其中的固體和液體微粒共同組成的多相體系，粒子的粒徑多在0.001～10μm之間，作為氣溶膠最基本的光學(xué)特性之一，氣溶膠光學(xué)厚度(AOD)被定義為介質(zhì)的消光系數(shù)在垂直方向上的積分，是描述氣溶膠對光衰減作用的主要因子。氣溶膠的消光效應(yīng)主要來源于細(xì)顆粒物的貢獻(xiàn)，粗粒子消光貢獻(xiàn)相對較低，因此通過AOD的變化趨勢可以有效反映我國 PM2.5的污染狀況。

為了全面了解我國陸地上空氣溶膠的時空分布及變化規(guī)律，本文利用2005-2012年8年MODIS氣溶膠光學(xué)厚度產(chǎn)品，對中國氣溶膠的時空演變特征進(jìn)行了系統(tǒng)分析，著重分析了京津冀、長三角、華中、成渝四個重點區(qū)域及北京、石家莊、保定、邢臺四個典型城市的AOD演變特征，揭示了我國近八年來氣溶膠污染的時空演變過程。本文利用2005-2012年MODIS氣溶膠產(chǎn)品研究了我國氣溶膠光學(xué)厚度的空間分布和時間演變特征，主要結(jié)論如下:

(1)2005-2012年期間我國AOD總體呈現(xiàn)“先降低，后升高”的趨勢，2005-2007年期間是我國氣溶膠污染最嚴(yán)重的時期。全國氣溶膠污染嚴(yán)重地區(qū)主要集中在京津冀、長三角、華中、成渝等經(jīng)濟(jì)發(fā)達(dá)，人口密集的地區(qū)。

(2)京津冀、長三角、華中、成渝四個重點區(qū)域AOD年均值在不同時期內(nèi)均遠(yuǎn)高于平均全國平均水平，2005-2012年期間總體變化趨勢與全國基本一致。

(3)2005-2012年期間北京、石家莊、保定、邢臺四個城市AOD年均值持續(xù)在0.8附近的高值區(qū)波動，變化幅度較小，AOD遠(yuǎn)高于全國及重點區(qū)域平均值。四城市在2005-2012年期間AOD年均值變化趨勢存在顯著的同步性。

利用AOD表征PM2.5的污染特征有一定的科學(xué)依據(jù)，但是AOD的大小受近地面PM2.5濃度、濕度及氣溶膠垂直分布特征等多種因素的影響，用AOD直接表征近地面PM2.5濃度的大小存在一定的不確定性。因此在下一步研究中有必要對AOD進(jìn)行濕度和垂直等因素的訂正，提高AOD表征近地面PM2.5濃度的準(zhǔn)確性。

［1］高大偉，徐宏輝，郁珍艷，等.MODIS氣溶膠光學(xué)厚度在臨安大氣顆粒物監(jiān)測中的應(yīng)用［J］.環(huán)境科學(xué)研究，2012(07):739-744.

［2］傅為，鄒崢嶸.基于MODIS氣溶膠光學(xué)厚度反演的PM_(10)濃度監(jiān)測［D］.中南大學(xué)，2009.

［3］林海峰，辛金元，張文煜，等.北京市近地層顆粒物濃度與氣溶膠光學(xué)厚度相關(guān)性分析研究［J］.環(huán)境科學(xué)，2013(03):826-834.

［4］鄭卓云，陳良富，鄭君瑜，等.高分辨率氣溶膠光學(xué)厚度在珠三角及香港地區(qū)區(qū)域顆粒物監(jiān)測中的應(yīng)用研究［J］.環(huán)境科學(xué)學(xué)報，2011(06):1154-1161.

［5］Wang J，Christopher S A.Intercomparison between satellite-derived aerosol optical thickness and PM 2.5 mass:Implications for air quality studies［J］.Geophysical Research Letters，2003，30(21):2-5.

［6］Liu Y，Sarnat J A，Kilaru V，et al.Estimating ground-level PM2.5in the eastern United States using satellite remote sensing.［J］.Environmental Science ＆Technology，2005，39(9):3269-3278.

［7］King M D，Kaufman Y J，Tanre D，et al.Remote Sensing of Tropospheric Aerosols from Space:Past，Present，and Future ［J］.Bulletin of the American Meteorological Society，American Meteorological Society，1998，80(11):2229-2259.

［8］毛節(jié)泰，李成才，張軍華，等.MODIS衛(wèi)星遙感北京地區(qū)氣溶膠光學(xué)厚度及與地面光度計遙感的對比［J］.應(yīng)用氣象學(xué)報，2002(S1):127-135.

［9］Chu D A.Global monitoring of air pollution over land from the Earth Observing System-Terra Moderate Resolution Imaging Spectroradiometer(MODIS)［J］.Journal ofGeophysicalResearch， American GeophysicalUnion， 2003， 108(D21):1-18.

［10］李成才，毛節(jié)泰，劉啟漢，等.MODIS衛(wèi)星遙感氣溶膠產(chǎn)品在北京市大氣污染研究中的應(yīng)用［A］.全國優(yōu)秀青年氣象科技工作者學(xué)術(shù)研討會［C］.中國湖南長沙:2006:12.

［11］Wang Z，Chen L，Tao J，等.Satellite-based estimation of regional particulate matter(PM)in Beijing using vertical-and-RH correcting method［J］.Remote Sensing of Environment，Elsevier Inc.，2010，114(1):50-63.

［12］關(guān)佳欣，李成才.我國中東部主要地區(qū)氣溶膠光學(xué)厚度的分布和變化［J］.北京大學(xué)學(xué)報(自然科學(xué)版)，2010(02):185-191.

［13］Lin J，Nielsen C P，Zhao Y，et al.Recent changes in particulate air pollution over China observed from space and the ground:effectiveness of emission control.［J］.Environmental science technology，2010，44(20):7771-7776.

［14］王淑蘭，云雅如，胡君等.情景分析技術(shù)在制定區(qū)域大氣復(fù)合污染控制方案中的應(yīng)用研究［J］.環(huán)境與可持續(xù)發(fā)展，2012，37(4):14-20.

［15］解淑艷，王瑞斌，李建軍等.現(xiàn)代化環(huán)境空氣質(zhì)量監(jiān)測網(wǎng)絡(luò)構(gòu)想［J］.環(huán)境與可持續(xù)發(fā)展，2012，37(4):26-31.

［16］王金南，寧淼，孫亞梅.區(qū)域大氣污染聯(lián)防聯(lián)控的理論與方法分析［J］.環(huán)境與可持續(xù)發(fā)展，2012，37(5):5-10.

［17］寧淼，孫亞梅，楊金田.國內(nèi)外區(qū)域大氣污染聯(lián)防聯(lián)控管理模式分析［J］.環(huán)境與可持續(xù)發(fā)展，2012，37(5):11-18.

［18］劉炳江.“十二五”主要大氣污染物總量減排對策措施［J］.環(huán)境與可持續(xù)發(fā)展，2012，37(4):5-10.

［19］李培，王新，柴發(fā)合等.我國城市大氣污染控制綜合管理對策［J］.環(huán)境與可持續(xù)發(fā)展，2011，36(5):8-14.

［20］王帥，丁俊男，王瑞斌等.關(guān)于我國環(huán)境空氣質(zhì)量監(jiān)測點位設(shè)置的思考［J］.環(huán)境與可持續(xù)發(fā)展，2012，37(4):21-25.

［21］嚴(yán)剛，燕麗.“十二五”我國大氣顆粒物污染防治對策［J］.環(huán)境與可持續(xù)發(fā)展，2011，36(5):20-23.

［21］薛文博，王金南，楊金田，等.國內(nèi)外空氣質(zhì)量模型研究進(jìn)展［J］.環(huán)境與可持續(xù)發(fā)展，2013，38(3):14-20.