李正強(qiáng)，許 華，張 瑩，李 莉，李東輝，侯偉真，呂 陽(yáng)，顧行發(fā)，陳興峰，陳 澄，張鳳霞，馬

中國(guó)科學(xué)院遙感與數(shù)字地球研究所國(guó)家環(huán)境保護(hù)衛(wèi)星遙感重點(diǎn)實(shí)驗(yàn)室，北京100101

灰霾污染已經(jīng)成為中國(guó)社會(huì)發(fā)展過(guò)程中的重大環(huán)境問(wèn)題。2013年1月發(fā)生的嚴(yán)重灰霾污染，波及華北及華東大部分地區(qū)，覆蓋區(qū)域超過(guò)140萬(wàn)平方公里，直接影響人口數(shù)量超過(guò)8億，其中受到嚴(yán)重污染影響的約2.5億，其覆蓋范圍、持續(xù)時(shí)間、污染程度等方面均屬全球罕見(jiàn)，引起了國(guó)內(nèi)外的廣泛關(guān)注。2013年9月，國(guó)務(wù)院頒布《大氣污染防治行動(dòng)計(jì)劃》，提出了10條共35項(xiàng)措施推動(dòng)大氣污染治理和空氣質(zhì)量改善，并由此對(duì)大氣污染、尤其是灰霾污染監(jiān)測(cè)提出了迫切需求。

霾是大氣污染的一種現(xiàn)象，主要由大氣中的懸浮顆粒物質(zhì)造成，通常顆粒會(huì)通過(guò)吸水而膨脹，導(dǎo)致體積增大，對(duì)太陽(yáng)光的衰減(消光)增強(qiáng)，使能見(jiàn)度顯著降低。在城市區(qū)域，由于人為顆粒物排放較多，為霧轉(zhuǎn)化為霾提供了有力條件，霾常常與霧同時(shí)或相繼發(fā)生，因此籠統(tǒng)上稱(chēng)為霧霾［1］?；姻渤Ｌ刂篙^高濃度的大氣細(xì)顆粒物污染現(xiàn)象，由于細(xì)顆粒物散射較短波長(zhǎng)太陽(yáng)輻射的能力較強(qiáng)，往往導(dǎo)致天空呈現(xiàn)橙灰等顏色，故此被稱(chēng)為灰霾。學(xué)術(shù)上，將大氣中懸浮顆粒統(tǒng)一稱(chēng)為氣溶膠，它特指1納米～100微米之間的大氣顆粒物質(zhì)，1納米以下常被認(rèn)為是氣－粒轉(zhuǎn)換的邊界尺寸［2］，而100微米以上的大氣懸浮顆粒物質(zhì)(雨、雪、雹等)則會(huì)快速降落至地面。粒徑在0.1～10微米的典型氣溶膠粒子在大氣光學(xué)、大氣輻射、大氣化學(xué)、大氣污染和云物理學(xué)等方面具有重要影響，是大氣中最重要的成分之一，其分布可隨時(shí)空有較大變化，且具有明顯的空間垂直分布特性，大部分顆粒物都聚集在大氣邊界層之內(nèi)［3］，濃度高時(shí)可造成顯著的灰霾污染。遙感手段通過(guò)觀測(cè)大氣顆粒物與輻射(例如太陽(yáng)光)的相互作用，可以了解氣溶膠在大氣層中的含量等參數(shù)，是獲得灰霾時(shí)空分布信息的有力工具。

1 基礎(chǔ)數(shù)據(jù)

1.1 衛(wèi)星數(shù)據(jù)

灰霾遙感的基礎(chǔ)數(shù)據(jù)目前主要是衛(wèi)星觀測(cè)到的氣溶膠光學(xué)厚度(AOD)，AOD指垂直大氣柱中氣溶膠消光系數(shù)隨高度(從地面到大氣頂)的積分。這里，我們以較易獲取的中分辨率成像光譜儀(MODIS)的AOD數(shù)據(jù)進(jìn)行說(shuō)明。MODIS傳感器搭載在美國(guó)宇航局于1999年發(fā)射的TERRA(過(guò)境時(shí)間為地方時(shí)10：30 am)和2002年發(fā)射的AQUA衛(wèi)星(過(guò)境時(shí)間為地方時(shí)13：30 pm)上，并通過(guò)星上轉(zhuǎn)發(fā)器實(shí)時(shí)向地面直接廣播，全球可以免費(fèi)接收數(shù)據(jù)，因此在許多國(guó)家和地區(qū)得到了廣泛的使用。

MODIS傳感器軌道高度距離地球705 km，最大空間分辨率為250 m，最大掃描寬度2 330 km。MODIS有36個(gè)波段，波長(zhǎng)范圍為0.4～14.4 μm，它采用暗目標(biāo)法［4］等原理在陸地上空反演氣溶膠光學(xué)厚度［5］。該方法通過(guò)2.12 μm波段來(lái)識(shí)別植被、土壤等暗目標(biāo)，并認(rèn)為暗目標(biāo)在紅光(0.66 μm)和藍(lán)光(0.47 μm)波段地表反射率較低，衛(wèi)星入瞳處的信號(hào)在這種情況下與氣溶膠光學(xué)厚度線性相關(guān)，并據(jù)此反演0.55 μm波長(zhǎng)的AOD。由于需要通過(guò)多像素平均來(lái)提高信噪比，因此目前的AOD產(chǎn)品空間分辨率為10 km。在進(jìn)行氣溶膠光學(xué)厚度反演前，需要進(jìn)行云檢測(cè)，因此在判定有云的地表上空將無(wú)AOD數(shù)據(jù)。此外，由于暗目標(biāo)法原理的限制，在冰、雪等亮地表上空，將無(wú)法進(jìn)行反演。在這種情況下，MODIS提供了基于深藍(lán)算法［6］的氣溶膠光學(xué)厚度數(shù)據(jù)。MODIS的氣溶膠光學(xué)厚度產(chǎn)品經(jīng)過(guò)了幾次重大版本修改后，已較為穩(wěn)定，并在全球范圍內(nèi)開(kāi)展了較為系統(tǒng)的檢驗(yàn)。驗(yàn)證結(jié)果顯示，與地基氣溶膠自動(dòng)觀測(cè)網(wǎng)(AERONET)［7］比較，陸地上空氣溶膠光學(xué)厚度的反演誤差約為 ±(0.05 ±0.15%)［8］。

1.2 氣象資料

由于目前衛(wèi)星傳感器資料的限制，研究灰霾污染遙感監(jiān)測(cè)仍需要用到一些輔助的同步氣象資料。通過(guò)模式精細(xì)化模擬后，氣象資料具有較好的時(shí)空分辨率可以與衛(wèi)星數(shù)據(jù)配合使用。為了獲得大氣邊界層高度(PBLH)和相對(duì)空氣濕度(RH)，利用全球氣象預(yù)報(bào)系統(tǒng)同化國(guó)際氣象組織基準(zhǔn)站的觀測(cè)數(shù)據(jù)，獲得再分析氣象場(chǎng)。隨后，利用中尺度預(yù)報(bào)模式可將再分析資料的時(shí)空分辨率精細(xì)化，獲得配合衛(wèi)星遙感使用的氣象場(chǎng)資料。

2 遙感監(jiān)測(cè)方法及系統(tǒng)設(shè)計(jì)

衛(wèi)星遙感獲得的基礎(chǔ)數(shù)據(jù)是氣溶膠光學(xué)厚度(AOD)，是整層大氣顆粒物含量的光學(xué)參數(shù)。通過(guò)AOD可進(jìn)一步獲得灰霾污染狀況的定量數(shù)值，并估計(jì)出空氣質(zhì)量指數(shù)。首先介紹遙感光學(xué)參數(shù)與近地面環(huán)境監(jiān)測(cè)參數(shù)(灰霾指數(shù)和空氣質(zhì)量指數(shù))之間的轉(zhuǎn)化關(guān)系，并建立相關(guān)估算方法;隨后，介紹灰霾污染遙感監(jiān)測(cè)系統(tǒng)的搭建，以及遙感影像的數(shù)據(jù)處理與監(jiān)測(cè)應(yīng)用。

2.1 衛(wèi)星遙感灰霾指數(shù)

衛(wèi)星遙感灰霾指數(shù)定義為0～4的整數(shù)，分別代表無(wú)灰霾、輕微灰霾、輕度灰霾、中度灰霾、重度灰霾。依據(jù)氣象學(xué)標(biāo)準(zhǔn)，灰霾指數(shù)與近地面層的大氣能見(jiàn)度相關(guān)：能見(jiàn)度在5～10 km為輕微霾，3～5 km為輕度霾，2～3 km為中度霾，2 km以下為重度霾(《霾觀測(cè)和預(yù)報(bào)等級(jí)》(QX/T 113—2010))。通常，除了云和霧等情況外，大氣氣溶膠的消光(吸收和散射)對(duì)能見(jiàn)度有決定性影響，因此理論上，在給定的大氣模式條件下，大氣消光與氣溶膠光學(xué)厚度有較好相關(guān)性。Baumer等人［9］根據(jù)氣溶膠垂直消光特性與近地面能見(jiàn)度的觀測(cè)給出了氣溶膠光學(xué)厚度與能見(jiàn)度和邊界層高度的經(jīng)驗(yàn)公式。本文基于該經(jīng)驗(yàn)公式獲得AOD至能見(jiàn)度的轉(zhuǎn)化關(guān)系，并考慮大氣邊界層的變化情況，模擬了氣溶膠光學(xué)厚度與灰霾指數(shù)的對(duì)應(yīng)關(guān)系，如圖1所示。在不同邊界層高度(PBLH)情況下，灰霾等級(jí)有著不同的變化，但灰霾指數(shù)總隨光學(xué)厚度增長(zhǎng)而增加。

圖1 灰霾指數(shù)與氣溶膠光學(xué)厚度的關(guān)系

2.2 衛(wèi)星遙感估算空氣質(zhì)量指數(shù)

空氣質(zhì)量指數(shù)(AQI)可定義為0～500之間的整數(shù)，分別代表優(yōu)(0～50)、良(51～100)、輕度污染(101～150)、中度污染(151～200)、重度污染(201～300)和嚴(yán)重污染(301～500)。2013年12月，環(huán)保部通報(bào)中國(guó)中東部大面積灰霾污染其首要污染物均為 PM2.5(http：//news．hangzhou．com．cn/gnxw/content/2013 －12/04/content_4990452．htm)，而PM2.5包含于PM10中，因此嚴(yán)重灰霾污染時(shí)空氣質(zhì)量指數(shù)主要由顆粒物濃度決定。氣溶膠光學(xué)厚度與顆粒物濃度PM10(空氣動(dòng)力學(xué)直徑小于10 μm的顆粒物)有較好的相關(guān)性?；贙oelemeijer等人［10］的研究結(jié)果，AOD與PM10線性相關(guān)的典型斜率為214 μg/m3。此外，不同天氣情況下，大氣層結(jié)構(gòu)情況不同，該關(guān)系還受氣溶膠垂直分布和吸濕特性的影響［11-12］。研究假設(shè)氣溶膠垂直分布為e指數(shù)形式，并滿(mǎn)足Kotchenmther等人［13］提出的氣溶膠消光系數(shù)與環(huán)境相對(duì)濕度的關(guān)系，來(lái)校正上述影響。已有的研究表明，以上2個(gè)假設(shè)基本適用于中國(guó)地區(qū)［14］。最后，依據(jù)《環(huán)境空氣質(zhì)量評(píng)價(jià)技術(shù)規(guī)范》(試行)(HJ 633—2013)，嚴(yán)重灰霾污染時(shí)可用PM10換算AQI。因此與圖1類(lèi)似，模擬了氣溶膠光學(xué)厚度與嚴(yán)重污染時(shí)的空氣質(zhì)量指數(shù)等級(jí)之間的關(guān)系。

圖2給出了不同空氣濕度條件下的空氣質(zhì)量指數(shù)等級(jí)。模擬結(jié)果表明，對(duì)于給定的邊界層高度，在相同空氣質(zhì)量指數(shù)的情況下，RH越高，AOD越大。這可由著名的吸濕增長(zhǎng)特性解釋?zhuān)磭?yán)重灰霾污染時(shí)空氣質(zhì)量指數(shù)由顆粒物干成分質(zhì)量決定，隨著空氣濕度增加，干成分吸水導(dǎo)致整個(gè)顆粒物體積增大，消光增強(qiáng)，表現(xiàn)為AOD增大。

圖2 空氣質(zhì)量指數(shù)等級(jí)與氣溶膠光學(xué)厚度的關(guān)系

2.3 灰霾污染遙感監(jiān)測(cè)系統(tǒng)設(shè)計(jì)

灰霾污染遙感監(jiān)測(cè)系統(tǒng)需要具有自動(dòng)化的數(shù)據(jù)處理能力和可視化的監(jiān)視、分析環(huán)境，因此系統(tǒng)結(jié)構(gòu)設(shè)計(jì)和系統(tǒng)配置需要考慮多方面的需求。基于MODIS數(shù)據(jù)實(shí)例，簡(jiǎn)要介紹灰霾遙感監(jiān)測(cè)系統(tǒng)的設(shè)計(jì)方案。

如圖3所示，系統(tǒng)的處理流程可分為數(shù)據(jù)源、預(yù)處理、產(chǎn)品生產(chǎn)和共享發(fā)布：1)MODIS上午星(Terra)和下午星(Aqua)可實(shí)現(xiàn)中國(guó)大部分地區(qū)每日兩次的遙感監(jiān)測(cè);使用MOD04大氣數(shù)據(jù)產(chǎn)品，不僅包含覆蓋陸地和海洋的氣溶膠光學(xué)厚度數(shù)據(jù)集，還提供了云掩碼數(shù)據(jù)。2)產(chǎn)品生產(chǎn)之前需要經(jīng)過(guò)預(yù)處理步驟，每景MOD04數(shù)據(jù)的觀測(cè)時(shí)間大約為5 min，而感興趣區(qū)可能出現(xiàn)在好幾景影像上，因此，利用拼接、裁剪、幾何投影和格式轉(zhuǎn)化等自動(dòng)化工具實(shí)現(xiàn)影像的流水線處理。由于灰霾污染監(jiān)測(cè)中云的影響較大，預(yù)處理中也包含云掩膜的功能。3)系統(tǒng)針對(duì)氣溶膠光學(xué)厚度信息和氣象分析資料，基于上述灰霾污染遙感模型，計(jì)算感興趣區(qū)的大氣污染狀況(灰霾指數(shù)和空氣質(zhì)量指數(shù)等級(jí))，同時(shí)生產(chǎn)產(chǎn)品的元數(shù)據(jù)。4)產(chǎn)品成果的共享發(fā)布以圖件和簡(jiǎn)報(bào)的形式體現(xiàn)。系統(tǒng)具有自動(dòng)運(yùn)行能力，并可利用GIS信息完成初步制圖和簡(jiǎn)報(bào)生成功能，最后，經(jīng)過(guò)專(zhuān)業(yè)人員修改和審核的監(jiān)測(cè)結(jié)果可由Web形式發(fā)布共享。

系統(tǒng)實(shí)現(xiàn)上，以一種低成本、高可靠、按需擴(kuò)展的分布式存儲(chǔ)為核心，構(gòu)建灰霾污染遙感監(jiān)測(cè)流水線系統(tǒng)。如圖4所示，該系統(tǒng)包括：1)分布式存儲(chǔ)集群。每個(gè)Linux存儲(chǔ)節(jié)點(diǎn)以高速局域網(wǎng)互聯(lián)，節(jié)點(diǎn)可動(dòng)態(tài)加入存儲(chǔ)集群，實(shí)現(xiàn)可持續(xù)擴(kuò)展，性能可隨著系統(tǒng)擴(kuò)展而增強(qiáng)。2)Linux計(jì)算節(jié)點(diǎn)。將高性能的服務(wù)器部署在高速局域網(wǎng)上，實(shí)現(xiàn)數(shù)據(jù)下載、影像預(yù)處理。3)產(chǎn)品生產(chǎn)。任務(wù)調(diào)度使用系統(tǒng)的守護(hù)進(jìn)程進(jìn)行數(shù)據(jù)監(jiān)控或定時(shí)器觸發(fā)。4)圖形工作站。將高性能的圖形顯示工作站部署在局域網(wǎng)(互聯(lián)網(wǎng))上，進(jìn)行自動(dòng)化與人工相結(jié)合的產(chǎn)品制圖和簡(jiǎn)報(bào)生產(chǎn)。5)產(chǎn)品提交給分布式存儲(chǔ)集群后，可通過(guò)Web服務(wù)共享。流水線采用了一種基于異構(gòu)分布式網(wǎng)絡(luò)環(huán)境的系統(tǒng)架構(gòu)，集群存儲(chǔ)空間和計(jì)算資源被虛擬成為一系列的服務(wù)，具有經(jīng)濟(jì)、靈活、穩(wěn)定、擴(kuò)展性強(qiáng)等技術(shù)特點(diǎn)，可以滿(mǎn)足衛(wèi)星遙感監(jiān)測(cè)灰霾的業(yè)務(wù)化需求。

圖3 灰霾污染衛(wèi)星遙感監(jiān)測(cè)系統(tǒng)數(shù)據(jù)流程

圖4 灰霾污染衛(wèi)星遙感監(jiān)測(cè)系統(tǒng)架構(gòu)

3 結(jié)果及驗(yàn)證

3.1 研究區(qū)域及遙感結(jié)果

基于上述方法，選擇2013年10月5—7日的MODIS數(shù)據(jù)開(kāi)展監(jiān)測(cè)實(shí)驗(yàn)，獲得了華北地區(qū)陸地上空灰霾指數(shù)以及重度來(lái)霾時(shí)的空氣質(zhì)量指數(shù)等級(jí)(圖5)。需要指出的是，在解譯遙感分布圖時(shí)，應(yīng)首先注意圖上空白區(qū)域是有云－霧覆蓋的區(qū)域，屬于衛(wèi)星數(shù)據(jù)目前無(wú)法反演AOD的范圍。通常，可將云－霧分為受污染型和未受污染型。由于灰霾、云－霧的空間連續(xù)分布特性，云－霧的受污染類(lèi)型在遙感尺度下不太可能出現(xiàn)突變，因此對(duì)于空白區(qū)域，可以依據(jù)周?chē)姻驳燃?jí)和梯度變化情況作出估計(jì)。例如，10月5日，河北、北京、河南北部等空白區(qū)域應(yīng)是受到污染的云－霧(從同期地面觀測(cè)可證實(shí)，AQI較高)。10月6日，山西區(qū)域的空白應(yīng)是受污染的云－霧。10月7日河北南部及其與河南、山東交界處的空白區(qū)域應(yīng)是受污染的云－霧，從空白區(qū)域中心處殘留(河北的東南邊界)的衛(wèi)星遙感結(jié)果看，該區(qū)域污染程度較高，延續(xù)了6日的重度污染。

圖5 灰霾指數(shù)(左)和空氣質(zhì)量指數(shù)(右)衛(wèi)星遙感圖(上：2013-10-05;中：2013-10-06;下：2013-10-07)

從圖5上衛(wèi)星遙感有值的區(qū)域可以看到：1)衛(wèi)星遙感可以很好地反映灰霾污染程度的變化。例如，10月6日，山東西北部衛(wèi)星遙感顯示出污染等級(jí)的連續(xù)梯度變化，在自西北到東南方向約100 km的距離內(nèi)，從嚴(yán)重灰霾污染逐漸變化為無(wú)灰霾，同時(shí)，從重度或嚴(yán)重污染變?yōu)閮?yōu)。2)衛(wèi)星遙感有助于發(fā)現(xiàn)可能的傳輸通道。在10月5、7日?qǐng)D上，從山西省的灰霾污染狀況可明顯看出，其主要污染帶是中部(由縱貫?zāi)媳钡拇笸?、忻州、太原、臨汾、運(yùn)城5大盆地組成的狹長(zhǎng)地帶)，當(dāng)該縱貫線上為輕或中度灰霾時(shí)(10月5日)，與其相鄰的山西東部和西部無(wú)灰霾，當(dāng)其為重度灰霾時(shí)(10月7日)，與其相鄰的區(qū)域?yàn)檩p度灰霾。若有合適的風(fēng)向、風(fēng)力等條件，將有可能形成灰霾向西北或東南傳輸?shù)耐ǖ馈?)衛(wèi)星遙感可獲得宏觀的灰霾空間分布概況。從該研究區(qū)域在這個(gè)時(shí)段的監(jiān)測(cè)看，河北(尤其南部)、山西(尤其南北縱貫盆地)、河南(北部)、山東(西北部)及其交界區(qū)域是灰霾的重污染區(qū)域，而內(nèi)蒙、遼寧、山東東部等區(qū)域的灰霾污染則較少。

3.2 衛(wèi)星面覆蓋觀測(cè)與地面點(diǎn)監(jiān)測(cè)結(jié)果對(duì)比

為驗(yàn)證衛(wèi)星遙感結(jié)果，選擇該區(qū)域內(nèi)12個(gè)環(huán)保監(jiān)測(cè)站同期的AQI與衛(wèi)星估算的AQI對(duì)比(表1)，選擇城市區(qū)域的多個(gè)監(jiān)測(cè)站點(diǎn)計(jì)算該城市平均空氣質(zhì)量指數(shù)。很多情況下，這些城市監(jiān)測(cè)站分布在城市中心區(qū)域約20～30 km的范圍內(nèi)，由于衛(wèi)星遙感產(chǎn)品的空間分布率約為10 km，因此衛(wèi)星估計(jì)的空氣質(zhì)量指數(shù)選取3×3像元(約30×30 km)，并取均值進(jìn)行對(duì)比。表1中列出了衛(wèi)星遙感監(jiān)測(cè)與環(huán)保站點(diǎn)監(jiān)測(cè)的AQI等級(jí)對(duì)比結(jié)果。對(duì)灰霾污染期間空氣質(zhì)量指數(shù)等級(jí)的相對(duì)誤差進(jìn)行統(tǒng)計(jì)，有效樣本數(shù)19個(gè)，結(jié)果顯示衛(wèi)星遙感與地面監(jiān)測(cè)站的AQI等級(jí)監(jiān)測(cè)結(jié)果之間平均相對(duì)誤差為16%。在88%的樣本中，衛(wèi)星遙感能夠較好地反映空氣質(zhì)量指數(shù)等級(jí)，誤差不大于1級(jí)，這說(shuō)明衛(wèi)星遙感能夠較真實(shí)地獲得灰霾污染狀況。此外，基于2013年10月整個(gè)華北地區(qū)的衛(wèi)星遙感結(jié)果與12個(gè)地面監(jiān)測(cè)站進(jìn)行了系統(tǒng)對(duì)比驗(yàn)證(樣本數(shù)160個(gè))，結(jié)果表明，衛(wèi)星遙感獲得灰霾污染情況下的AQI等級(jí)與站點(diǎn)觀測(cè)值的平均偏差為0.7級(jí)。

表1 衛(wèi)星遙感AQI與地面站AQI監(jiān)測(cè)結(jié)果對(duì)比

從表1還可以看出，衛(wèi)星遙感監(jiān)測(cè)灰霾污染情況下的AQI等級(jí)與地面監(jiān)測(cè)站結(jié)果的差異普遍在±1級(jí)之間，但通常低于站點(diǎn)監(jiān)測(cè)結(jié)果，其原因主要是衛(wèi)星遙感與地面監(jiān)測(cè)之間的時(shí)空代表性有所不同。首先，采用的衛(wèi)星數(shù)據(jù)大概在每天10：30 am過(guò)境，衛(wèi)星數(shù)據(jù)僅代表過(guò)境時(shí)刻的瞬時(shí)值，而地面監(jiān)測(cè)站公布的是空氣質(zhì)量指數(shù)的日均值，考慮到灰霾污染程度的日變化(一般夜間嚴(yán)重)，雖然10：30 am有較好的代表性，但不可避免仍可能與日均值存在差異。其次，衛(wèi)星監(jiān)測(cè)的面平均數(shù)據(jù)與環(huán)境監(jiān)測(cè)站的點(diǎn)監(jiān)測(cè)數(shù)據(jù)之間存在區(qū)域代表性的問(wèn)題，尤其是目前的城市監(jiān)測(cè)站多布置在污染程度相對(duì)較高的城市中心，對(duì)整個(gè)區(qū)域的代表性較弱。例如，從空間分布上來(lái)看(如圖6)，濟(jì)南市南面臨山，環(huán)保監(jiān)測(cè)主要在市區(qū)進(jìn)行，而衛(wèi)星監(jiān)測(cè)明顯體現(xiàn)出其南部山區(qū)空氣質(zhì)量較好，在約30 km的范圍內(nèi)，AQI等級(jí)從4級(jí)下降至1級(jí)，因此空間平均后衛(wèi)星監(jiān)測(cè)AQI低于站點(diǎn)監(jiān)測(cè)值。

衛(wèi)星監(jiān)測(cè)唐山地區(qū)10月6日結(jié)果較差的主要原因是，當(dāng)天遙感手段僅在唐山地區(qū)東南區(qū)域獲得一個(gè)像素的有效值，導(dǎo)致唐山市污染觀測(cè)不完全，造成監(jiān)測(cè)值偏低。有云區(qū)域無(wú)法獲得衛(wèi)星監(jiān)測(cè)結(jié)果，是當(dāng)前灰霾衛(wèi)星遙感的不足之處。

圖6 濟(jì)南市衛(wèi)星遙感空氣質(zhì)量指數(shù)等級(jí)分布圖(2013-10-06)

4結(jié)論

衛(wèi)星遙感可以獲得大范圍內(nèi)灰霾污染的定量監(jiān)測(cè)結(jié)果，具有客觀、統(tǒng)一、快速區(qū)域覆蓋等優(yōu)勢(shì)。介紹了基于氣溶膠光學(xué)厚度結(jié)合氣象數(shù)據(jù)的灰霾污染遙感監(jiān)測(cè)方法，包括衛(wèi)星遙感灰霾指數(shù)的算法，以及衛(wèi)星遙感估算灰霾污染情況下空氣質(zhì)量指數(shù)等級(jí)的算法。隨后，以中分辨率成像光譜儀MODIS為例，提出了灰霾污染遙感監(jiān)測(cè)系統(tǒng)的設(shè)計(jì)方案，介紹了數(shù)據(jù)流程和系統(tǒng)架構(gòu)。基于所發(fā)展的方法和系統(tǒng)，選擇華北地區(qū)，在2013年10月開(kāi)展了監(jiān)測(cè)驗(yàn)證，通過(guò)對(duì)比地面監(jiān)測(cè)站的空氣質(zhì)量指數(shù)，顯示灰霾污染時(shí)88%的衛(wèi)星遙感監(jiān)測(cè)結(jié)果誤差不大于1個(gè)AQI等級(jí)，平均的AQI偏差為0.7級(jí)，能夠比較真實(shí)地反映灰霾污染狀況。在此基礎(chǔ)上分析認(rèn)為，衛(wèi)星遙感可以很好地反映灰霾污染程度的變化，有助于發(fā)現(xiàn)可能的傳輸通道，并獲得灰霾的空間分布狀況。此外，站點(diǎn)觀測(cè)數(shù)據(jù)僅代表觀測(cè)點(diǎn)處局部的情況，而遙感方法可基于統(tǒng)一標(biāo)準(zhǔn)提供面覆蓋的灰霾污染狀況，從這一點(diǎn)上來(lái)說(shuō)，遙感方法比現(xiàn)有地面監(jiān)測(cè)更適合于大范圍內(nèi)灰霾污染的對(duì)比評(píng)價(jià)。這些宏觀信息的獲取能夠?yàn)榛姻脖O(jiān)測(cè)、預(yù)警、防控等，尤其是全國(guó)范圍內(nèi)的灰霾污染治理效果評(píng)價(jià)提供重要支撐。

［1］吳兌．近十年中國(guó)灰霾天氣研究綜述［J］．環(huán)境科學(xué)學(xué)報(bào)，2012，2：257-269.

［2］Kulmala M，Kontkanen J，Junninen H，et al．Direct observations of atmospheric aerosol nucleation［J］．Science，2013，339(6 122)：943-946.

［3］張婉春，張瑩，呂陽(yáng)，等．利用激光雷達(dá)探測(cè)灰霾天氣大氣邊界層高度［J］．遙感學(xué)報(bào)，2013，17(4)：981-992.

［4］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)：1 357-1 381.

［5］Levy R C，Remer L A，Mattoo S，et al．Secondgeneration operational algorithm：Retrieval of aerosol properties overland from inversion ofModerate Resolution Imaging Spectroradiometer spectral reflectance［J］．Journal of Geophysical Research：Atmospheres，2007，112(D13).

［6］Hsu N C，Tsay S C，King M D，et al．Aerosol properties over bright-reflecting source regions［J］．Geoscience and Remote Sensing，2004，42(3)：557-569.

［7］Holben B N，Eck T F，Slutsker I，et al．AERONET．A federated instrument network and data archive for aerosolcharacterization［J］．Remote Sensing of Environment，1998，66(1)：1-16.

［8］Levy R C，Remer L A，Dubovik O．Global aerosol optical properties and application to Moderate Resolution Imaging Spectroradiometer aerosol retrieval over land［J］．Journal of Geophysical Research：Atmospheres，2007，112(D13).

［9］Baumer D，Vogel B，Versick S，et al．Relationship of visibility，aerosol optical thichness and aerosol size distribution in an ageing air mass over South-West Germany［J］．Atmos Environ，2008，42：989-998.

［10］Koelemeijer R B A，Homan C D，Matthijsen J．Comparison of spatial and temporal variations of aerosol optical thickness and particulate matter in Europe［J］．Atmos Environ，2006，40：5 304-5 315.

［11］鄭卓云，陳良富，鄭君瑜，等．高分辨率氣溶膠光學(xué)厚度在珠三角及香港地區(qū)區(qū)域顆粒物監(jiān)測(cè)中的應(yīng)用研究［J］．環(huán)境科學(xué)學(xué)報(bào)，2011，6：1 154-1 161.

［12］何秀，鄧兆澤，李成才，等．MODIS氣溶膠光學(xué)厚度產(chǎn)品在地面PM10監(jiān)測(cè)方面的應(yīng)用研究［J］．北京大學(xué)學(xué)報(bào)：自然科學(xué)版，2010，2：178-184.

［13］KotchenrutherR A，HobbsP V，Hegg D A．Humidification factors for atmospheric aerosols off the mid-Atlantic coast of the United States［J］．Journal of Geophysical Research-Atmospheres，1999(D2)：2 239-2 251.

［14］李成才，毛節(jié)泰，劉啟漢，等．MODIS衛(wèi)星遙感氣溶膠產(chǎn)品在北京市大氣污染研究中的應(yīng)用［J］．中國(guó)科學(xué)(D 輯：地球科學(xué))，2005，S1：177-186.