景 悅,孫艷玲*,徐 昊,陳 莉,張 輝,高 爽,付宏臣,毛 健

?

基于混合效應(yīng)模型的京津冀地區(qū)PM2.5日濃度估算

景 悅1,孫艷玲1*,徐 昊2,陳 莉1,張 輝1,高 爽1,付宏臣1,毛 健1

(1.天津師范大學(xué)地理與環(huán)境科學(xué)學(xué)院,天津 300387；2.寧夏大學(xué)經(jīng)濟(jì)管理學(xué)院,寧夏 銀川 750021)

利用2016年182d的MODIS 3km AOD數(shù)據(jù)與地面監(jiān)測數(shù)據(jù),評估了混合效應(yīng)模型不同參數(shù)組合的模擬性能,得出模型在解釋AOD-PM2.5關(guān)系時,對時間序列變異的解釋能力要比空間差異更佳.在此基礎(chǔ)上,利用混合效應(yīng)模型建立京津冀地區(qū)每日的AOD-PM2.5關(guān)系,模型擬合2為0.92,交叉驗證調(diào)整2為0.85,均方根誤差(RMSE)為12.30μg/m3,平均絕對誤差(MAE)為9.73μg/m3,說明模型擬合精度較高.基于此模型估算的2016年京津冀地區(qū)年均PM2.5濃度為42.98μg/m3,暖季(4月1日~10月31日)為43.35μg/m3, 冷季(11月1日~3月31日)為38.52μg/m3,與同時期的地面監(jiān)測數(shù)據(jù)差值分別為0.59,0.7,5.29μg/m3.空間上,京津冀地區(qū)的PM2.5濃度呈現(xiàn)南高北低的特征,有一條明顯的西南-東北走向的高值區(qū).研究結(jié)果表明,基于每日混合效應(yīng)模型可以準(zhǔn)確評估京津冀地區(qū)的地面PM2.5濃度,且模型估算的PM2.5濃度分布狀況為區(qū)域大氣污染防治提供了基礎(chǔ)的數(shù)據(jù)支撐.

MODIS AOD；PM2.5；混合效應(yīng)模型；京津冀地區(qū)

自2013年國家建立環(huán)境監(jiān)測站點以來,分析區(qū)域尺度的PM2.5污染狀況主要是利用站點獲取的數(shù)據(jù)進(jìn)行插值得到.由于環(huán)境監(jiān)測站點離散分布,且站點之間空間距離過大,會導(dǎo)致PM2.5濃度的空間分布存在較大的偏差,嚴(yán)重影響PM2.5數(shù)據(jù)的準(zhǔn)確性[1].而衛(wèi)星遙感覆蓋范圍廣,空間分辨率高,通過反演可以準(zhǔn)確獲取無站點分布區(qū)的PM2.5濃度,彌補(bǔ)了地面監(jiān)測空間上的不足[2].衛(wèi)星遙感反演主要使用的是氣溶膠光學(xué)厚度 (AOD),其定義為介質(zhì)的消光系數(shù)在垂直方向上的積分,是衡量氣溶膠的重要因子,可以用來表征大氣渾濁度或氣溶膠的總含量,有效的反映區(qū)域空氣污染程度[3].

目前,國內(nèi)外已有不少學(xué)者開展了利用AOD估算地面PM2.5濃度的研究.研究表明AOD-PM2.5之間存在明顯的正相關(guān)關(guān)系[4],在此基礎(chǔ)上,有學(xué)者通過對AOD進(jìn)行垂直訂正和濕度訂正,進(jìn)一步提高了兩者的相關(guān)性[5-8].近幾年,不同的學(xué)者分別采用多元線性回歸模型[9]、人工神經(jīng)網(wǎng)絡(luò)模型[10-11]、地理加權(quán)模型[12-14]、土地利用回歸模型[15-17]、貝葉斯統(tǒng)計模型[18]等方法利用AOD估算地面PM2.5濃度,均得到了較好的估算結(jié)果,證明了AOD估算PM2.5的可行性.2011年Lee等[19]首次使用混合效應(yīng)模型對美國東北部地區(qū)地面PM2.5濃度進(jìn)行了估算,模型擬合2高達(dá)0.97,充分證明了混合效應(yīng)模型有較好的模擬能力.后來Lee和其他學(xué)者又將此模型運用到其他地區(qū),均得到了很好驗證[20-21].以往研究中,AOD數(shù)據(jù)主要是采用分辨率10km的MODIS產(chǎn)品,很難滿足用戶高分辨率的需求. 2014年末,美國國家航空航天局 (NASA)發(fā)布了MODIS 3km AOD產(chǎn)品,為反演更高精度的地面PM2.5濃度提供了可能. Xie等[22]、Ma等[23]分別利用此數(shù)據(jù)產(chǎn)品對北京和長江三角洲地區(qū)進(jìn)行了PM2.5濃度估算,都得到了高精度的模擬結(jié)果.

自2013年國家大氣污染防治計劃實施以來,全國空氣質(zhì)量總體好轉(zhuǎn).從區(qū)域分布來看,京津冀地區(qū)大氣環(huán)境質(zhì)量雖同比有所改善,但仍是我國大氣污染最重的區(qū)域.其中,2016年統(tǒng)計的空氣質(zhì)量較差的10大城市中,6個就位于該區(qū)域,且重污染天氣主要發(fā)生在冷季.因此,本文以京津冀地區(qū)為研究區(qū),選用MODIS 3km AOD產(chǎn)品,采用混合效應(yīng)模型建立了每日變化的AOD-PM2.5關(guān)系,并使用交叉驗證檢驗?zāi)Ｐ偷臏?zhǔn)確性,在此基礎(chǔ)上估算了2016年京津冀地區(qū)年均、暖季及冷季的PM2.5濃度,為研究京津冀地區(qū)大氣污染治理提供數(shù)據(jù)支撐.

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

1.1 MODIS AOD數(shù)據(jù)

MODIS AOD數(shù)據(jù)來自于Terra衛(wèi)星,該衛(wèi)星是美國宇航局地球觀測系統(tǒng)計劃(EOS)中的第一顆衛(wèi)星,于1999年12月18日發(fā)射升空,2000年2月24日開始提供全球遙感監(jiān)測資料[6],其上搭載的中等分辨率光譜儀(MODIS),每天覆蓋全球一次,過境時間為地方時10:30am,掃描寬度2330m,星下點分辨率250m,有36個光譜通道,可反演云屬性、氣溶膠、溫度、臭氧、陸地覆被等一系列產(chǎn)品,該數(shù)據(jù)產(chǎn)品可供全球用戶免費使用.NASA提供的標(biāo)準(zhǔn)level 2級AOD產(chǎn)品空間分辨率為10km,對全球空氣污染監(jiān)測有很大應(yīng)用價值[24], 2014年發(fā)布的3km AOD產(chǎn)品已被證明其2/3的數(shù)據(jù)都在允許誤差范圍內(nèi),可以被用于全球范圍的氣溶膠變化研究[25].張西雅等[26]利用AERONET數(shù)據(jù)對Terra 3km衛(wèi)星遙感氣溶膠產(chǎn)品在京津冀地區(qū)的適用性進(jìn)行了驗證,表明該產(chǎn)品滿足區(qū)域環(huán)境空氣質(zhì)量評估的需要.本文使用的Terra MOD04_3k數(shù)據(jù)產(chǎn)品是從美國國家航空航天局的戈達(dá)德宇宙飛行中心(LAADS) (https://ladsweb. nascom.nasa.gov/)網(wǎng)站上下載得到,除去2016年2月19~27日的缺失數(shù)據(jù),共得到357d AOD數(shù)據(jù).

1.2 地面監(jiān)測PM2.5數(shù)據(jù)

PM2.5監(jiān)測數(shù)據(jù)由中國環(huán)境監(jiān)測總站(http:// 113.108.142.147:20035/emcpublish)提供.地面PM2.5數(shù)據(jù)采用的是微量震蕩天平法(TEOM)或者β射線法監(jiān)測得到,自動監(jiān)測設(shè)備實施全面質(zhì)量控制,保證了數(shù)據(jù)的準(zhǔn)確性.本文利用京津冀地區(qū)80個環(huán)境監(jiān)測站點進(jìn)行研究,各站點分布情況如圖1所示.研究收集了2016年1月1日~12月31日每天10:00和11:00的PM2.5監(jiān)測數(shù)據(jù),為匹配Terra衛(wèi)星過境的時間,將每個站點2個時段的數(shù)據(jù)平均,最終得到366d PM2.5_10:30數(shù)據(jù).

圖1 京津冀地區(qū)環(huán)境監(jiān)測站點分布

1.3 模型建立與檢驗

首先建立線性回歸模型來模擬AOD-PM2.5的關(guān)系.公式如下:

PM2.5=+′AOD (1)

式中:、為固定的截距與斜率,線性回歸模型是利用最小二乘法對一個自變量與因變量建立簡單的線性擬合關(guān)系.而實際上,由于每天的氣象條件、衛(wèi)星反演環(huán)境、局部污染情況的不同,AOD-PM2.5關(guān)系可能存在時間與空間上的變化.基于此,本文建立了2個獨立混合效應(yīng)模型來更清楚的反映AOD-PM2.5關(guān)系在時間序列及空間位置的變化.方程如下:

PM2.5ij=(+U)+(+V)′AOD+(2)

PM2.5ij=(+U)+(+V)′AOD+ (3)

(U,V)~[(00),S1](U,V)~[(00),S2]

式中:公式(2)是AOD隨每日變化的混合效應(yīng)模型,公式(3)是站點影響下的混合效應(yīng)模型.、為固定的截距與斜率,隨機(jī)參數(shù)U、V是每天變化的隨機(jī)截距與隨機(jī)斜率,U、V則是由站點之間差異導(dǎo)致的隨機(jī)截距與斜率,ε為每日的隨機(jī)殘差,ε是每個站點的隨機(jī)殘差,∑1、∑2分別為每日及站點變化下的隨機(jī)效應(yīng)對應(yīng)的協(xié)方差矩陣.通過公式(2)可以得出混合效應(yīng)模型對時間序列的解釋能力,也進(jìn)一步得出AOD-PM2.5關(guān)系有多少是由時間變化(如每天氣象條件差異等)所引起的.通過公式(3)則可進(jìn)一步解釋站點差異(地形、土地利用狀況、局部污染物排放差異)對AOD-PM2.5關(guān)系的貢獻(xiàn)度.模型擬合的效果通過計算模擬值與地面監(jiān)測值的2、平均絕對誤差(MAE)、均方根誤差(RMSE)來進(jìn)行比較,并與線性回歸模型進(jìn)行對比驗證.

交叉驗證(CV)可以用來檢驗?zāi)Ｐ偷倪^度擬合,其原理是將樣本數(shù)據(jù)隨機(jī)均勻的分成組,其中1組留作測試集,1組作為訓(xùn)練集,交叉驗證重復(fù)次,保證每一組樣本都被驗證一次,將次交叉驗證的結(jié)果取平均得到最后的結(jié)果.其中,10折交叉驗證最為常用.本文也采用10折交叉驗證的方法檢驗線性回歸模型與混合效應(yīng)模型的性能.

2 結(jié)果與分析

2.1 描述性統(tǒng)計

將研究區(qū)分為3km′3km的網(wǎng)格,用每一個網(wǎng)格去匹配其上的PM2.5站點數(shù)據(jù),如果網(wǎng)格內(nèi)對應(yīng)多個站點,則取站點間PM2.5濃度的平均值,以保證每一個AOD網(wǎng)格對應(yīng)唯一的PM2.5值.剔除每日匹配少于2對的數(shù)據(jù),最終得到182d 2841對數(shù)據(jù).AOD- PM2.5匹配結(jié)果的描述性統(tǒng)計見表1,監(jiān)測PM2.5日均濃度與AOD日均值的相關(guān)性為0.53(圖2).另外,本文將2016年的4~10月劃分為暖季, 11~3月劃分為冷季,其中暖季有2505對數(shù)據(jù),冷季有336對數(shù)據(jù).

2016年京津冀地區(qū)地面監(jiān)測的年均PM2.5濃度為43.57μg/m3(表1),暖季的PM2.5濃度范圍為1~ 201.5μg/m3,平均值為44.05μg/m3,冷季PM2.5濃度范圍為3~151.5μg/m3,平均值為33.23μg/m3. 京津冀地區(qū)2016年AOD總體平均值為0.74(表1),暖季的AOD范圍為0~4.74,平均值為0.76,冷季的AOD范圍為0~1.31,平均值為0.28.地面監(jiān)測站點插值的PM2.5空間分布特征(圖3)與2016年同時期AOD的空間分布特征類似,整體上呈現(xiàn)東南高西北低的特點.

圖2 2016年京津冀地區(qū)每日平均AOD-PM2.5匹配

表1 AOD-PM2.5匹配的描述性統(tǒng)計

2.2 模型擬合及交叉驗證

利用不同參數(shù)組合評估了混合效應(yīng)模型的擬合效果(表2).首先基于站點(每日)混合效應(yīng)模型內(nèi)部對比可知,混合效應(yīng)模型中僅考慮單因子時,隨機(jī)斜率的擬合效果比僅考慮隨機(jī)截距要好,兩者都考慮可以更進(jìn)一步提高精度,其次,對比公式(2)、(3)建立的每日混合效應(yīng)模型與站點混合效應(yīng)模型可知,混合效應(yīng)模型在解釋AOD-PM2.5關(guān)系時,對時間序列變異解釋能力要比空間差異更佳.進(jìn)一步得出AOD-PM2.5的關(guān)系其實更多的是隨時間變化,而不是空間上變化,佐證了本文進(jìn)行日估算的合理性.對比線性回歸模型與混合效應(yīng)模型,可以看出任一種混合效應(yīng)模型都比線性回歸模型擬合效果好,尤其考慮每日的混合效應(yīng)模型最佳.

圖4中,線性回歸模型與站點混合效應(yīng)模型的2分布區(qū)間為0.2~1,均出現(xiàn)了2個離群值,中位數(shù)2約為0.53左右,說明模型擬合效果一般,且模擬效果極不穩(wěn)定.相反,基于每日混合效應(yīng)模型擬合2分布區(qū)間為0.73~1,出現(xiàn)了3個小于第一個四分位數(shù)的離散值,中位數(shù)2高達(dá)0.92,尾長較短,證明模型擬合較好.同樣在MAE、RMSE分布中,線性回歸模型與站點混合效應(yīng)模型不僅中位值高,而且分布區(qū)間跨度較大,每日混合效應(yīng)模型的MAE與RMSE值均較低且分布集中,表明其誤差較小,證明基于每日混合效應(yīng)模型擬合效果最佳.

表2 模型擬合結(jié)果對比

對于不同模型交叉驗證的結(jié)果可見 (圖5),線性回歸模型與基于站點混合效應(yīng)模型的CV2為0.30(±0.03),相比驗證前2降低了0.23,說明兩模型存在嚴(yán)重的過度擬合.基于每日混合效應(yīng)模型CV2為0.85,相比驗證之前的2降低了0.07,雖然也存在部分過度擬合現(xiàn)象,但交叉驗證調(diào)整后的相關(guān)性依然很高,說明使用3km AOD產(chǎn)品結(jié)合每日混合效應(yīng)模型算法在京津冀地區(qū)可以很好的估算地面PM2.5濃度.

在模型驗證的基礎(chǔ)上,最終利用公式(2)建立一個既考慮AOD對隨機(jī)斜率也考慮其對隨機(jī)截距影響的每日混合效應(yīng)模型,得到2016年京津冀地區(qū)2841個模擬值及每天對應(yīng)的182個AOD-PM2.5關(guān)系方程.固定斜率與截距分別為37.71 (<0.0001)和19.17(<0.0001),每日變化的隨機(jī)斜率與隨機(jī)截距標(biāo)準(zhǔn)差分別為29.83,17.87,隨機(jī)殘差的標(biāo)準(zhǔn)差為12.94.整個研究區(qū)模型擬合2為0.92, RMSE為12.30μg/m3,MAE為9.73μg/m3,交叉驗證調(diào)整R為0.85.說明模擬的PM2.5濃度可以解釋平均85%的地面PM2.5濃度.模型擬合效果以不同城市進(jìn)行統(tǒng)計(表3),不同城市的2范圍為0.86~0.96,RMSE為11.20~13.31μg/m3,MAE為6.05~11.49μg/m3,交叉驗證調(diào)整后的2集中在0.70~0.91,可以看出各城市的模型擬合都有一個高的相關(guān)性和較低的誤差,其中保定市模型擬合效果最好,張家口市與承德市效果較差,但CV2依然高達(dá)0.70以上.

圖4 基于站點統(tǒng)計的不同模型擬合的R2、MAE、RMSE

圖5 不同模型交叉驗證結(jié)果

表3 基于城市統(tǒng)計的每日混合效應(yīng)模型擬合效果

2.3 京津冀地區(qū)PM2.5估算及大氣污染分布特征

利用每日混合效應(yīng)模型建立的AOD-PM2.5關(guān)系估算3km分辨率下的地面PM2.5濃度,如圖6a所示, 2016年京津冀地區(qū)年均PM2.5濃度為42.98μg/m3,值域范圍在0~528μg/m3之間,對比同時段的PM2.5監(jiān)測站點的插值圖(圖3a),可以看出,兩者的平均濃度近似均為(43±0.57)μg/m3,但估算的PM2.5濃度值域范圍更大,尤其對無站點地區(qū)得出了更為詳細(xì)的PM2.5空間信息.同時,本文還利用AOD產(chǎn)品對研究區(qū)暖季和冷季的地面PM2.5濃度進(jìn)行估算(圖6),進(jìn)一步探討模型模擬的季節(jié)適應(yīng)性問題.京津冀地區(qū)2016年暖季PM2.5濃度范圍為0~528μg/m3,平均濃度為43.35μg/m3.冷季PM2.5濃度范圍在0~267μg/m3之間,平均濃度為38.52μg/m3.與暖季、冷季站點平均PM2.5濃度(表1)對比,暖季相差+0.7μg/m3,冷季相差+5.29μg/m3,結(jié)果偏差較小,說明對研究區(qū)整體上進(jìn)行暖季、冷季PM2.5濃度估算精度較高.從空間分布上看,研究區(qū)2016年P(guān)M2.5濃度呈現(xiàn)南高北低的特征,其中暖季呈現(xiàn)高度的空間聚集性,西南-東北走向有大片PM2.5濃度高值區(qū),西北部為低值區(qū),冷季紅色高值區(qū)向北擴(kuò)散,說明污染范圍進(jìn)一步擴(kuò)大.另外,從暖季、冷季PM2.5濃度的空間分布 (圖6)可知,暖季的模擬結(jié)果空間覆蓋效果較好,冷季由于北京市城區(qū),天津市城區(qū)及張家口山區(qū)等地區(qū)AOD數(shù)據(jù)缺失,導(dǎo)致模型估算的地面PM2.5濃度出現(xiàn)部分缺失.

雖然2016年京津冀地區(qū)整體上空氣質(zhì)量整體超出了國家空氣質(zhì)量二級標(biāo)準(zhǔn),但不同城市的空氣污染狀況還是有所差異.基于城市統(tǒng)計的模型估算PM2.5濃度狀況如表4所示,其中空氣質(zhì)量較好的兩大城市為張家口市、承德市,其PM2.5濃度分別為28.01,27.04μg/m3.其余10個城市的年均PM2.5濃度都高于35μg/m3.暖季PM2.5污染狀況與全年類似.冷季北京市及天津市郊區(qū)的PM2.5濃度分別為31.41,32.01μg/m3,張家口市大部分區(qū)域冷季PM2.5濃度為34.21μg/m3.邯鄲市、邢臺市、衡水市、滄州市、石家莊市、保定市、廊坊市的空氣污染較嚴(yán)重,其PM2.5濃度都大于35μg/m3,承德市、唐山市、秦皇島市冷季的空氣質(zhì)量較好.不同城市空氣污染差異很大,表明京津冀地區(qū)在大氣污染防治中,不同城市應(yīng)因地制宜,并著重治理空氣污染嚴(yán)重地區(qū).PM2.5的年均及暖季、冷季的污染濃度大不相同,說明PM2.5濃度的季節(jié)差異很明顯,進(jìn)一步表明了空氣污染分季節(jié)探討及治理的必要性與合理性.

表4 基于城市統(tǒng)計的模型估算PM2.5濃度(μg/m3)

3 討論

3.1 模型構(gòu)建

本文采用的混合效應(yīng)模型是在固定效應(yīng)模擬AOD-PM2.5關(guān)系的基礎(chǔ)上,進(jìn)一步用每日的隨機(jī)斜率與隨機(jī)截距來代表影響每天AOD-PM2.5關(guān)系變化的眾多不確定因素(如每日氣象條件、衛(wèi)星反演條件等),并使用每日隨機(jī)殘差項訂正系統(tǒng)模擬誤差,因此擬合精度較高.本文建立的混合效應(yīng)模型未考慮站點的隨機(jī)截距項,主要由于站點數(shù)量少,并不能分配到每一個網(wǎng)格,大量網(wǎng)格并沒有站點覆蓋,而且站點的效應(yīng)僅可以糾正站點這一點的偏差,并不能糾正3km范圍的偏差,對于無站點的地區(qū)而言,站點效應(yīng)并沒有起到訂正的效果,反而引起了一些不確定性誤差.其次,混合效應(yīng)模型考慮站點與不考慮站點擬合效果近似,這與Xie等[22]的結(jié)論一致.本文進(jìn)一步通過對比不同混合效應(yīng)模型效果,證明了混合效應(yīng)模型對時間序列的模擬最好,但對空間差異的解釋力較差,因此采用日估計并未考慮站點影響來盡可能減少不必要的誤差.

3.2 模型表現(xiàn)對比

本文模型的交叉驗證調(diào)整2為0.85, 比郝靜等[27]利用10km產(chǎn)品對京津冀內(nèi)陸平原的模擬結(jié)果(CV2=0.7)好,說明3km數(shù)據(jù)產(chǎn)品要比10km數(shù)據(jù)更加精確.與Xie等[22]利用同樣數(shù)據(jù)對2013年北京地區(qū)的反演結(jié)果(CV2=0.75)相比,本文結(jié)果較好.但是與使用混合效應(yīng)模型對美國地區(qū)的擬合結(jié)果[19-21](CV2=0.88,0.92)相比,京津冀地區(qū)的模擬結(jié)果要差一些,這可能與我國地面監(jiān)測數(shù)據(jù)較少有關(guān).與其他學(xué)者采用不同的方法[12,14,18]對京津冀地區(qū)地面PM2.5濃度的反演結(jié)果對比,每日混合效應(yīng)模型同樣得出較好的模擬結(jié)果.

本文通過使用最新的MODIS 3km AOD產(chǎn)品估算了高精度的地面PM2.5濃度,為京津冀地區(qū)大氣污染防治提供了更詳細(xì)空間信息.盡管之前也有學(xué)者使用了1km分辨率的AOD產(chǎn)品數(shù)據(jù)估計地面PM2.5濃度,但其對MAIAS 1km AOD的反演結(jié)果與MODIS AOD結(jié)果作對比顯示MODIS數(shù)據(jù)反演結(jié)果與地面監(jiān)測的PM2.5濃度有更高的相關(guān)性[28].另一方面,MAIAS 1km AOD產(chǎn)品并未開放,基本很難獲取,相比之下,高質(zhì)量的MODIS 3km AOD產(chǎn)品對全球公開且可在NASA官網(wǎng)免費下載,這為我國學(xué)者研究AOD-PM2.5關(guān)系提供了巨大便利.

3.3 模型的不確定性

每日混合效應(yīng)模型不需要輔助多種要素,而是僅用AOD來建立模型,將氣象、區(qū)域差異都考慮在隨機(jī)效應(yīng)中,用更少的數(shù)據(jù)也可以得到同樣高的模擬效果,這是其他模型不可比擬的優(yōu)勢.但也有一定的局限性,隨機(jī)效應(yīng)直接模擬不確定因素的結(jié)果與在模型中加入輔助因子模擬的結(jié)果是否一致,未來需要進(jìn)一步的驗證.

另外,在相關(guān)研究中有些學(xué)者并沒有使用MODIS 3km AOD產(chǎn)品,而依舊使用10km AOD產(chǎn)品[28].這主要原因在于3km AOD目前只發(fā)布了暗像元算法(DT)產(chǎn)品,并沒有深藍(lán)算法(DB)產(chǎn)品,DT算法對植被覆蓋地區(qū)反演效果最好,對冬季城市、雪等亮地表及重污染天氣反演較差[29],使得冬季AOD數(shù)據(jù)缺失嚴(yán)重,導(dǎo)致本文冷季的AOD-PM2.5匹配數(shù)量少且部分區(qū)域估算結(jié)果缺失.這也是目前使用3km產(chǎn)品研究的局限[23],但吳健生等[14]發(fā)現(xiàn),對中國地區(qū)而言,使用不同模型的反演或者在模型中加入氣象等輔助因子都沒有改善AOD數(shù)據(jù)分辨率對反演結(jié)果影響大,這也是本文選擇高分辨率3km AOD產(chǎn)品的原因.未來需要考慮使用MODIS的雙星協(xié)同數(shù)據(jù),并對3km的暗像元算法與10km的深藍(lán)算法進(jìn)一步做數(shù)據(jù)融合,得到更多的匹配數(shù)據(jù),也會使得模型得到進(jìn)一步的提升.

4 結(jié)論

4.1 對比混合效應(yīng)模型不同參數(shù)組合的模擬結(jié)果發(fā)現(xiàn),混合效應(yīng)模型對時間變化解釋更好,對站點變化解釋力較差,其中考慮AOD對隨機(jī)截距與隨機(jī)斜率共同影響的每日混合效應(yīng)模型最佳.本文發(fā)展的2016年京津冀地區(qū)的每日混合效應(yīng)模型擬合2為0.92,交叉驗證調(diào)整2為0.85,RMSE為12.30μg/m3, MAE為9.73μg/m3,均證明模型有很好的模擬能力且估算結(jié)果與地面監(jiān)測PM2.5有很高的相關(guān)性.

4.2 利用每日混合效應(yīng)模型建立的AOD-PM2.5關(guān)系估算的2016年京津冀地區(qū)年均PM2.5濃度為42.98μg/m3,超出國家空氣質(zhì)量二級標(biāo)準(zhǔn)(35μg/m3)的22%.從空間分布來看,存在明顯的空間聚集性,有一條西南-東北走向的PM2.5高值區(qū),西北區(qū)為低值區(qū).對比同時段的監(jiān)測站點PM2.5濃度空間分布狀況,可知利用AOD估算的PM2.5空間分布狀況可以反映更詳細(xì)的空間信息,對地面PM2.5濃度的空間異質(zhì)性解釋更佳.暖季,冷季PM2.5平均濃度分別為43.35, 38.52μg/m3.季節(jié)PM2.5濃度的空間分布特征與年均基本一致.

4.3 同一個時期不同城市之間的PM2.5濃度差異明顯,同一個城市不同季節(jié)的PM2.5濃度變化較大,因此,針對京津冀地區(qū)PM2.5污染治理,應(yīng)該因地制宜,因時施宜.

[1] Hu X, Waller L A, Lyapustin A, et al. Estimating ground-level PM2.5, concentrations in the Southeastern United States using MAIAC AOD retrievals and a two-stage model [J]. Remote Sensing of Environment, 2014,140(1):220-232.

[2] 王銀牌,喻 鑫,謝廣奇.中國近15年氣溶膠光學(xué)厚度時空分布特征[J]. 中國環(huán)境科學(xué), 2018,38(2):426-434.

[3] 劉 浩,高小明,謝志英,等.京津冀晉魯區(qū)域氣溶膠光學(xué)厚度的時空特征[J]. 環(huán)境科學(xué)學(xué)報, 2015,35(5):1506-1511.

[4] Gupta P, Christopher S A, Wang J, et al. Satellite remote sensing of particulate matter and air quality assessment over global cities [J]. Atmospheric Environment, 2006,40(30):5880-5892.

[5] 王家成,朱成杰,朱 勇,等.北京地區(qū)多氣溶膠遙感參量與PM2.5相關(guān)性研究[J]. 中國環(huán)境科學(xué), 2015,35(7):1947-1956.

[6] 李成才,毛節(jié)泰,劉啟漢.利用MODIS遙感大氣氣溶膠及氣溶膠產(chǎn)品的應(yīng)用[J]. 北京大學(xué)學(xué)報(自然科學(xué)版), 2003,39(s1):108-117.

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

[8] 陳 輝,厲 青,王中挺,等.利用MODIS資料監(jiān)測京津冀地區(qū)近地面PM2.5方法研究[J]. 氣象與環(huán)境學(xué)報, 2014,30(5):27-37.

[9] 李嘯天,吳紹華,徐于月,等.江蘇省PM2.5質(zhì)量濃度的時空變化格局模擬[J]. 環(huán)境監(jiān)測管理與技術(shù), 2017,29(2):16-20.

[10] 王 敏,鄒 濱,郭 宇,等.基于BP人工神經(jīng)網(wǎng)絡(luò)的城市PM2.5濃度空間預(yù)測[J]. 環(huán)境污染與防治, 2013,35(9):63-66+70.

[11] 郭建平,吳業(yè)榮,張小曳,等. BP網(wǎng)絡(luò)框架下MODIS氣溶膠光學(xué)厚度產(chǎn)品估算中國東部PM2.5[J]. 環(huán)境科學(xué), 2013,34(3):817-825.

[12] Karimian H, Li Q, Li C C, et al. Daily estimation of fine particulate matter mass concentration through satellite based Aerosol Optical Depth [J]. 2017,IV-4/W2:175-181.

[13] 陳 輝,厲 青,張玉環(huán),等.基于地理加權(quán)模型的我國冬季PM2.5遙感估算方法研究[J]. 環(huán)境科學(xué)學(xué)報, 2016,36(6):2142-2151.

[14] 吳健生,王 茜,李嘉誠,等. PM2.5濃度空間分異模擬模型對比:以京津冀地區(qū)為例[J]. 環(huán)境科學(xué), 2017,38(6):2191-2201.

[15] 許 剛,焦利民,肖豐濤,等.土地利用回歸模型模擬京津冀PM2.5濃度空間分布[J]. 干旱區(qū)資源與環(huán)境, 2016,30(10):116-120.

[16] 吳健生,廖 星,彭 建,等.重慶市PM2.5濃度空間分異模擬及影響因子[J]. 環(huán)境科學(xué), 2015,36(3):759-767.

[17] 吳健生,謝舞丹,李嘉誠.土地利用回歸模型在大氣污染時空分異研究中的應(yīng)用[J]. 環(huán)境科學(xué), 2016,37(2):413-419.

[18] Lv B, Hu Y, Chang H H, et al. Daily estimation of ground-level PM2.5concentrations at 4km resolution over Beijing-Tianjin-Hebei by fusing MODIS AOD and ground observations [J]. Science of the Total Environment, 2017,580:235.

[19] Lee H J, Liu Y, Coull B A, et al. A novel calibration approach of MODIS AOD data to predict PM2.5concentrations [J]. Atmospheric Chemistry & Physics, 2011,11(11):9769-9795.

[20] Lee H J, Chatfield R B, Strawa A W. Enhancing the Applicability of Satellite Remote Sensing for PM2.5Estimation Using MODIS Deep Blue AOD and Land Use Regression in California, United States [J]. Environmental Science & Technology, 2016,50(12):6546.

[21] Kloog I, Chudnovsky A A, Just A C, et al. A New Hybrid Spatio-Temporal Model For Estimating Daily Multi-Year PM2.5Concentrations Across Northeastern USA Using High Resolution Aerosol Optical Depth Data [J]. Atmospheric Environment, 2014,95(1):581-590.

[22] Xie Y, Wang Y, Zhang K, et al. Daily estimation of ground-level PM2.5concentrations over Beijing using 3km resolution MODIS AOD [J]. Environmental Science & Technology, 2015,49(20):12280.

[23] Ma Z, Liu Y, Zhao Q, et al. Satellite-derived high resolution PM2.5, concentrations in Yangtze River Delta Region of China using improved linear mixed effects model [J]. Atmospheric Environment, 2016,133:156-164.

[24] Chu D A, Kaufman Y J, Zibordi G, et al. Global monitoring of air pollution over land from the Earth Observing System-Terra Moderate Resolution Imaging Spectroradiometer (MODIS) [J]. Journal of Geophysical Research Atmospheres, 2003,108(D21):4661.

[25] Munchak L A, Levy R C, Mattoo S, et al. MODIS 3km aerosol product: applications over land in an urban/suburban region [J]. Atmospheric Measurement Techniques, 2013,6(7):1747-1759.

[26] 張西雅,扈海波.京津冀地區(qū)氣溶膠時空分布及與城市化關(guān)系的研究 [J]. 大氣科學(xué), 2017,41(4):797-810.

[27] 郝 靜,孫 成,郭興宇,等.京津冀內(nèi)陸平原區(qū)PM2.5濃度時空變化定量模擬[J]. 環(huán)境科學(xué), 2018,39(4):1-13.

[28] Zongwei Ma, Xuefei Hu, Andrew M. Sayer, et al. Satellite-Based Spatiotemporal Trends in PM2.5Concentrations: China, 2004~2013 [J]. Environmental Health Perspectives, 2016,124(2):184.

[29] Ma X, Wang J, Yu F, et al. Can MODIS AOD be employed to derive PM2.5in Beijing-Tianjin-Hebei over China? [J]. Atmospheric Research, 2016,181:250-256.

Daily estimation of PM2.5concentrations based on mixed effects model in Beijing-Tianjin-Heibei region.

JING Yue1, SUN Yan-ling1*, XU Hao2, CHEN Li1, ZHANG Hui1, GAO Shuang1, FU Hong-chen1, MAO Jian1

(1.College of Geography and Environment Science, Tianjin Normal University, Tianjin 300387, China；2.College of Economics and Management, Ningxia University, Yinchuan 750021, China)., 2018,38(8)：2890~2897

The predictive performance of mixed effects model with different combinations of parameters was evaluated using the data of 182-day MODIS 3km AOD and ground monitoring concentration of PM2.5in 2016 year. The explanation capacity was better for temporal variations when explaining the relationship between AOD and PM2.5than for the spatial variations. Daily AOD-PM2.5relationship in Beijing-Tianjin-Hebei region was established based on the mixed effects model. The model predictions2, cross-validations2, RMSE and MAE were 0.92,0.85,12.30μg/m3, and 9.73μg/m3, respectively, indicated that the model showed good performance. The annual average PM2.5concentration in Beijing-Tianjin-Hebei region was 42.98μg/m3in 2016 based on the proposed model. The values for April-October and November-March were 43.35μg/m3and 38.52μg/m3. The differences were 0.59,0.7,5.29μg/m3, respectively, comparing with the ground monitoring PM2.5data at the corresponding period. PM2.5concentrations were higher in the south area and lower in the north area in Beijing-Tianjin-Hebei region with higher concentrations from southwestern to northeastern direction. PM2.5concentrations in the Beijing-Tianjin-Hebei region could be accurately evaluated based on the daily mixed effects model. The distribution of PM2.5concentrations estimated by the model could provide basic data support for the prevention and control of regional air pollution.

MODIS AOD；PM2.5；mixed effects model；Beijing-Tianjin-Hebei region

X513

A

1000-6923(2018)08-2890-08

景 悅(1994-),女,山西運城人,天津師范大學(xué)碩士研究生,主要從事大氣環(huán)境遙感研究.

2018-01-03

國家重點研發(fā)計劃青年項目(2016YFC0201700);天津市科技計劃項目(16YFXTSF00330);天津市應(yīng)用基礎(chǔ)與前沿技術(shù)研究計劃青年基金項目(16JCQNJC08600)

* 責(zé)任作者, 副教授,flyling99@163.com