基于MODIS數(shù)據(jù)的環(huán)渤海地區(qū)氣溶膠時空變化特征分析*

紀曉璐 廉麗姝 周甜甜 于鳳碩

(曲阜師范大學地理與旅游學院，山東 日照 276826)

基于MODIS數(shù)據(jù)的環(huán)渤海地區(qū)氣溶膠時空變化特征分析*

紀曉璐 廉麗姝#周甜甜 于鳳碩

(曲阜師范大學地理與旅游學院，山東 日照 276826)

利用2000—2015年美國國家航空航天局的Aqua衛(wèi)星 MODIS 2級氣溶膠產(chǎn)品MOD04數(shù)據(jù)，探討了環(huán)渤海地區(qū)氣溶膠光學厚度(AOD)的時空變化特征。結果表明：(1)環(huán)渤海地區(qū)的AOD年均值空間分布與海拔密切相關，高海拔地區(qū)AOD普遍較低，而低海拔地區(qū)AOD普遍較高。污染最嚴重的季節(jié)為夏季，AOD與空氣質量指數(shù)(AQI)相關性較高。(2)16年來，環(huán)渤海地區(qū)AOD年均值為0.465～0.682，2012年最高，2001年最低，平均值為0.592，遠高于全球陸上平均水平0.190。年內變化呈現(xiàn)雙峰型，最高值出現(xiàn)在7月，最低值出現(xiàn)在12月。AOD年內變化與風速呈負相關關系。

氣溶膠光學厚度 時空分布 變化特征 MODIS數(shù)據(jù) 環(huán)渤海地區(qū)

大氣氣溶膠粒子的直徑一般為0.001～100.000 μm，其組分復雜、來源多樣[1]。雖然氣溶膠粒子粒徑和質量微小，在大氣成分中所占質量分數(shù)很小，但其可以通過散射和吸收太陽輻射影響地氣系統(tǒng)的輻射平衡，通過凝結核影響云和降水過程，因此對區(qū)域乃至全球氣候都有重要影響。同時，氣溶膠還是局地光化學煙霧的重要原因，影響大氣能見度和人類健康。氣溶膠光學厚度(AOD)是描述氣溶膠對光的衰減作用的重要參數(shù)，通常AOD高預示著氣溶膠縱向積累的增長，導致大氣能見度降低。AOD在推算氣溶膠含量、評估大氣污染程度、確定氣溶膠氣候效應方面都具有重要作用[2-4]，因此其時空分布和變化特征引起了國內外學者的廣泛關注[5-8]。

環(huán)渤海地區(qū)包括遼東半島、山東半島和京津冀，同時延伸輻射到遼寧省、山東省、山西省以及內蒙古自治區(qū)中東部，是繼“珠江三角洲”、“長江三角洲”地區(qū)之后， 中國經(jīng)濟增長的“第三極”。因此，由工業(yè)生產(chǎn)、汽車尾氣排放等產(chǎn)生的空氣污染物日益增加，導致灰霾天氣頻發(fā)，對人體健康構成嚴重影響。環(huán)渤海地區(qū)關于AOD的研究主要集中在北京市及其周邊的京津冀地區(qū)[9-12]。本研究對環(huán)渤海地區(qū)的AOD研究擴大到了遼寧省、山東省和京津冀區(qū)域，利用MODIS 2級氣溶膠產(chǎn)品MOD04數(shù)據(jù)分析該地區(qū)AOD的時空變化特征，以期為環(huán)渤海地區(qū)大氣氣溶膠的區(qū)域性研究提供依據(jù)。

1 方 法

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

AOD數(shù)據(jù)來自美國國家航空航天局(NASA)網(wǎng)站(http://ladsweb.nascom.nasa.gov)發(fā)布的Aqua衛(wèi)星MODIS 2級氣溶膠產(chǎn)品MOD04數(shù)據(jù)，空間分辨率為10 km×10 km，時間范圍是從2000年3月到2016年2月(即2000—2015年)，波長取550 nm處。

數(shù)字高程模型(DEM)數(shù)據(jù)來自地理空間數(shù)據(jù)云網(wǎng)站(http://www.gscloud.cn/search)，由NASA的“奮進號”衛(wèi)星的航天飛機雷達地形測繪任務(SRTM)提供，分辨率為90 m，由此得到的環(huán)渤海地區(qū)地勢圖如圖1所示。

圖1 環(huán)渤海地區(qū)地勢圖Fig.1 Chorography of the areas around Bohai Sea region

風速數(shù)據(jù)來自中國氣象數(shù)據(jù)網(wǎng)(http://data.cma.cn/)?？諝赓|量指數(shù)(AQI)數(shù)據(jù)來自環(huán)境保護部數(shù)據(jù)中心(http://datacenter.mep.gov.cn/)。

1.2 研究方法

采用MOD04每日AOD數(shù)據(jù)，利用ENVI軟件和ArcGIS軟件對其進行預處理，作出2000—2015年環(huán)渤海地區(qū)AOD年均值空間分布圖和各季節(jié)多年平均值空間分布圖，并分析年際變化特征和年內變化特征。

季節(jié)劃分為春季3月至5月、夏季6月至8月、秋季9月至11月、冬季12月至次年2月。

選取了沈陽市、大連市、天津市、北京市、濟南市、石家莊市6個不同海拔高度且人口相對密集、經(jīng)濟較為發(fā)達的典型城市進行重點分析。

2 結果與分析

2.1 環(huán)渤海地區(qū)AOD空間分布特征

2.1.1 年均值的空間分布特征

2000—2015年間環(huán)渤海地區(qū)的AOD年均值空間分布如圖2所示。天津市、山東省西南部和河北省東南部等人口密集、經(jīng)濟發(fā)達的地區(qū)AOD較高；而河北省北部太行山沿線的AOD較低，原因是這里海拔高、植被多且重工業(yè)不發(fā)達。

與環(huán)渤海地區(qū)地勢圖(見圖1)比較發(fā)現(xiàn)，AOD分布與海拔密切相關，高海拔地區(qū)AOD普遍較低，而低海拔地區(qū)AOD普遍較高。當然，局地污染源對AOD也有影響，如遼寧省東部和山東省中部海拔相差不大，但是山東省中部的AOD要比遼寧省東部高。

圖2 2000—2015年環(huán)渤海地區(qū)AOD年均值空間分布Fig.2 Annual AOD spatial distribution around Bohai Sea region during 2000-2015

2.1.2 不同季節(jié)的空間分布特征

圖3分析了2000—2015年間環(huán)渤海地區(qū)不同季節(jié)的AOD多年平均值空間分布特征。由圖3(a)可見，春季AOD相對較低，均在0.9以下，0.6以上的區(qū)域主要集中在環(huán)渤海地區(qū)中部和南部，主要是工業(yè)污染的影響。由圖3(b)可見，夏季污染最嚴重，山東省西部、北部以及北京市、天津市、河北省南部的AOD均處于>1.2～1.5，這是因為夏季高溫高濕的天氣狀況有利于霧霾的形成，導致AOD增高。由圖3(c)可見，秋季AOD的高值區(qū)主要分布在山東省西南部和河北省南部。由圖3(d)可見，冬季AOD的高值區(qū)主要分布在山東省西部和遼寧省東北部，這兩個地區(qū)是海拔相對較高的丘陵或山地，冬季風速較小，不利于氣溶膠的擴散。

有不少學者研究了AOD與AQI以及PM2.5或PM10之間的相關關系[13-19]，但年均值的直接相關程度較低，在考慮季節(jié)后相關程度大有提高。這里分析了2014—2015年的環(huán)渤海地區(qū)6個典型城市的AQI與AOD相關性，發(fā)現(xiàn)春、夏、秋、冬四季的相關系數(shù)(置信度99%)分別為0.655、0.556、0.432、0.495，相關性較高，特別是春、夏季。

圖3 2000—2015年環(huán)渤海地區(qū)AOD各季節(jié)多年平均值空間分布Fig.3 Spatial distribution of seasonal AOD around Bohai Sea region during 2000-2015

圖4 環(huán)渤海地區(qū)AOD年均值年際變化Fig.4 Inter-annual variations of annual AOD around Bohai Sea region during 2000-2015

2.2 環(huán)渤海地區(qū)AOD時間變化特征

2.2.1 年際變化特征

由圖4可以看出，近16年來環(huán)渤海地區(qū)的AOD年均值為0.465～0.682，2012年最高，2001年最低，平均值為0.592，遠高于全球陸上平均水平0.190。其中濟南市、天津市、石家莊市AOD年均值較高，而沈陽市、北京市、大連市的AOD年均值較低。

2.2.2 年內變化特征

由圖5可以看出，環(huán)渤海地區(qū)的AOD年內變化基本呈現(xiàn)較為典型的雙峰型，最高值出現(xiàn)在7月，最低值出現(xiàn)在12月。各典型城市AOD年內變化與整個環(huán)渤海地區(qū)基本相似。

圖5 環(huán)渤海地區(qū)AOD月均值年內變化Fig.5 Variations of monthly AOD around Bohai Sea region

已有研究表明，AOD年內變化與風速的年內變化相關，當風速較大時，大氣流動加強，擴散與稀釋作用使得顆粒物的濃度顯著降低，AOD較低；而風速較小時，大氣流動減弱，顆粒物容易積聚，導致AOD上升[20-21]。這里研究了環(huán)渤海地區(qū)6個典型城市2014—2015年共24個月的風速與AOD之間的關系。結果表明，平均風速最大的大連市AOD最小，而平均風速最小的石家莊市AOD最大。對每個城市來說，風速大的月份AOD小，風速小的月份AOD大。由此可見，風速與AOD呈負相關關系。

3 結 語

(1) 環(huán)渤海地區(qū)的AOD年均值空間分布與海拔密切相關，高海拔地區(qū)AOD普遍較低，而低海拔地區(qū)AOD普遍較高。污染最嚴重的季節(jié)為夏季，AOD與AQI相關性較高。

(2) 近16年來環(huán)渤海地區(qū)的AOD年均值為0.465～0.682，2012年最高，2001年最低，平均值為0.592，遠高于全球陸上平均水平0.190。年內變化呈現(xiàn)雙峰型，最高值出現(xiàn)在7月，最低值出現(xiàn)在12月。AOD年內變化與風速呈負相關。

[1] 毛節(jié)泰,張軍華,王美華.中國大氣氣溶膠研究綜述[J].氣象學報,2002,60(5)：625-634.

[2] 董海鷹,劉毅,管兆勇.MODIS遙感中國近海氣溶膠光學厚度的檢驗分析[J].南京氣象學院學報,2007,30(3):328-337.

[3] 石廣玉,王標,張華.大氣氣溶膠的輻射與氣候效應[J].大氣科學,2008,32(4)：826-840.

[4] SOLOMON S,QIN D H,MANNING M,et al.Climate change 2007:the physical science basis,contribution of working group Ⅰ to the fourth assessment report of the intergovernmental panel on climate change[M].New York:Cambridge University Press,2007.

[5] ZHENG Y F,LIU J J,WU R J,et al.Seasonal statistical characteristics of aerosol optical properties at a site near a dust region in China[J].Journal of Geophysical Research:Atmospheres,2008,113(D16)：4641-4651.

[6] 柳晶.中國地區(qū)氣溶膠光學特性及輻射強迫的衛(wèi)星遙感觀測研究[D].南京：南京信息工程大學,2008.

[7] 鄧學良,鄧偉濤,何冬艷.近年來華東地區(qū)大氣氣溶膠的時空特征[J].大氣科學學報,2010,33(3)：347-354.

[8] RAMACHANDRAN S.Aerosol optical depth and fine mode fraction variations deduced from Moderate Resolution Imaging Spectroradiometer (MODIS) over four urban areas in India[J].Journal of Geophysical Research:Atmospheres,2007,112(D16)：6802-6811.

[9] 車慧正,石廣玉,張小曳.北京地區(qū)大氣氣溶膠光學特性及其直接輻射強迫的研究[J].中國科學院研究生院學報,2007,24(5)：699-704.

[10] 李本綱,冉陽,陶澍.北京市氣溶膠的時間變化與空間分布特征[J].環(huán)境科學學報,2008,28(7)：1425-1429.

[11] 李曉靜,劉玉潔,邱紅,等.利用MODIS資料反演北京及其周邊地區(qū)氣溶膠光學厚度的方法研究[J].氣象學報,2003,61(5)：580-591.

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

[13] 李成才,毛節(jié)泰,劉啟漢,等.利用MODIS光學厚度遙感產(chǎn)品研究北京及周邊地區(qū)的大氣污染[J].大氣科學,2003,27(5)：869-880.

[14] 羅云峰,呂達仁,周秀驥,等.30年來我國大氣氣溶膠光學厚度平均分布特征分析[J].大氣科學,2002,26(6)：721-730.

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

[16] REMER L A,KLEIDMAN R G,LEVY R C,et al.Global aerosol climatology from the MODIS satellite sensors[J].Journal of Geophysical Research:Atmospheres,2008，113(D14):1-17.

[17] 孫文文,江洪,陳洪濱,等.北京地區(qū)大氣氣溶膠光學特性與區(qū)域霧霾事件分析[J].環(huán)境污染與防治,2016,38(2)：1-6,12.

[18] 孫林,商曉青,孫長奎,等.城市地區(qū)氣溶膠光學厚度與空氣污染指數(shù)的相關性分析[J].山東科技大學學報(自然科學版),2008,27(3)：9-13.

[19] 謝志英,劉浩,唐新明.北京市MODIS氣溶膠光學厚度與PM10質量濃度的相關性分析[J].環(huán)境科學學報,2015,35(10)：3292-3299.

[20] CUHADAROGLU B,DEMIRCI E.Influence of some meteorological factors on air pollution in Trabzon city[J].Energy and Buildings,1997,25(3):179-184.

[21] 王淑英,張小玲.北京地區(qū)PM10污染的氣象特征[J].應用氣象學報,2002,13(1)：177-184.

TemporalandspatialvariationcharacteristicsanalysisofaerosolaroundBohaiSearegionbasedonMODISdata

JIXiaolu,LIANLishu,ZHOUTiantian,YUFengshuo.

(GeographyandTourismCollege,QufuNormalUniversity,RizhaoShandong276826)

Temporal and spatial variation characteristics of aerosol optical depth (AOD) around Bohai Sea region were discussed using MODIS Level 2 product MOD04 data of Aqua satellite from National Aeronautics and Space Administration of America during 2000 to 2015. Results showed that:(1) the spatial distribution of annual AOD around Bohai Sea region was related with altitude. That is to say,AOD was low in high altitude area but high in low altitude area. Pollution was the most serious in summer. And AOD was related with air quality index (AQI). (2) During the last 16 years,annual AOD were 0.465-0.682 with an average of 0.592,much higher than that of global level (0.190). The highest annual AOD was appeared in 2012 and the lowest in 2001. Variation of monthly AOD showed double peaks with highest value in July and lowest value in December. Variation of monthly AOD was negatively related with wind speed.

aerosol optical depth; spatial and temporal distribution; variation characteristics; MODIS data; region around Bohai Sea

紀曉璐，女，1990年生，碩士研究生，主要從事氣溶膠時空變化研究。#

。

*山東省自然科學基金資助項目(No.ZR2015DL001)。

10.15985/j.cnki.1001-3865.2017.11.016

2016-10-10)