王飛，韓業(yè)林，趙穎

太原市大氣顆粒物PM10和PM2.5多時間尺度變化規(guī)律研究

王飛1,2，韓業(yè)林2，趙穎3*

1. 山西大學生命科學學院，山西 太原 030006；2. 山西大學體育學院，山西 太原 030006；3.山西省環(huán)境科學研究院，山西 太原 030024

選取太原市城區(qū)10個監(jiān)測點2014—2016年P(guān)M10和PM2.5日變化數(shù)據(jù)，分析和探討了其時空變化特征，及其與人類經(jīng)濟活動的同步性規(guī)律；采用小波連續(xù)變換的功率譜方法識別顆粒物周期變化特征，采用可視化主成分分析法識別不同時間尺度下顆粒物變化的影響因素。結(jié)果表明，太原市大氣顆粒污染物PM10和PM2.5質(zhì)量濃度的變化存在明顯的時空差異，新興經(jīng)濟發(fā)展區(qū)較傳統(tǒng)老工業(yè)區(qū)污染嚴重，顆粒物污染程度在冬季較為嚴重。小波分析結(jié)果顯示，PM10和PM2.5時間序列的變化周期均以4～8 d的短周期為主（P＜0.05），污染物的質(zhì)量濃度變化與城市經(jīng)濟活動的周波動變化相一致；PM10和PM2.5質(zhì)量濃度最大值出現(xiàn)在周波動的中間時段，最小值出現(xiàn)在周末?？梢暬疨CA結(jié)果揭示，大氣顆粒物PM10和PM2.5季節(jié)性波動均受冬季影響較強；周波動周期內(nèi)均受周三影響最大；一天之內(nèi)PM10和PM2.5質(zhì)量濃度分別受夜晚和早晨影響最大，但白天顆粒物質(zhì)量濃度變化是造成其日變化特征的主要因素。研究結(jié)果有利于從不同時間尺度辨析能源城市大氣顆粒物污染的多變特征，有針對性地開展大氣污染防控，也可為管理部門制定相關(guān)標準和規(guī)范提供科學依據(jù)。

大氣顆粒物；多時間尺度；周期結(jié)構(gòu)；主成分分析；能源城市

隨著經(jīng)濟和城市化的快速發(fā)展，城市大氣污染防治問題已成為大氣環(huán)境領(lǐng)域的重要研究內(nèi)容（Li et al.，2015）。不同粒徑的大氣顆粒污染物（particulate matter，PM）經(jīng)累積導致的健康損害及與顆粒物不同化學組成密切相關(guān)的健康風險已經(jīng)成為關(guān)乎城市人群健康的重要熱點內(nèi)容（Lai，2016）。以氣溶膠形式懸浮在城市空氣中的顆粒物極具潛在毒理風險（Dergham et al.，2015），尤其是可吸入顆粒物的生理毒害以及經(jīng)地球化學循環(huán)通過干濕沉降在區(qū)域土壤、水體、動植物中的富集（Eeftens et al.，2015），使得城市大氣污染防治成為城市環(huán)境治理的重中之重。城市區(qū)域大氣污染以PM10和PM2.5為代表，其污染程度的加重，一方面使得大氣的能見度降低，霧霾等天氣出現(xiàn)頻繁，嚴重影響人們的正常生產(chǎn)生活及交通秩序；另一方面，過高的顆粒物質(zhì)量濃度提高了呼吸系統(tǒng)和循環(huán)系統(tǒng)疾病的發(fā)病率，帶來嚴重健康風險（王蕊等，2015），并且伴隨著靜穩(wěn)氣候條件持續(xù)時間的延長而顯著增加（Lai，2016）。許多研究或從不同季節(jié)研究大氣顆粒物污染來源、分布及潛在健康風險（王占山等，2015；Liu et al.，2014），或從不同月份探究了顆粒物質(zhì)量濃度變化情況（陶燕等，2014），或顆粒物在周循環(huán)變化的“周末效應”（雷瑜等，2015）以及顆粒物日質(zhì)量濃度變化（李名升等，2016；李名升等，2013）。這些研究從不同的時間尺度揭示了短期大氣顆粒物質(zhì)量濃度變化規(guī)律，但對顆粒物質(zhì)量濃度變化的周期劃分較為簡單，研究方法較為單一，在不同時間尺度上揭示人類經(jīng)濟活動和大氣顆粒物質(zhì)量濃度變化的密切關(guān)系仍不夠深入，缺乏系統(tǒng)性。

山西省太原市是全國典型的能源與重工業(yè)基地，其經(jīng)濟發(fā)展短時期內(nèi)還無法完全擺脫對煤炭的依賴，在可預見的時期內(nèi)煤炭仍將會在山西能源消費結(jié)構(gòu)中仍占有較高比例。由此造成的日益突出的大氣污染問題，政府雖然加強了工業(yè)污染的監(jiān)督和治理，積極調(diào)整能源產(chǎn)業(yè)結(jié)構(gòu)，加大城市綠化等一系列措施來改善城市大氣質(zhì)量，但大氣污染總體狀況仍不容樂觀（趙穎等，2015）。本文太原市以為例，重點分析城區(qū)范圍內(nèi)大氣顆粒污染物PM10和PM2.5質(zhì)量濃度的多尺度時空變化特征，探討其污染變化與人類日常生產(chǎn)生活的同步性規(guī)律及影響，并基于可視化的主成分分析（Principal Component Analysis，PCA）確定各監(jiān)測點PM10和PM2.5不同時間尺度下變化的特征，以期為預測評價太原市大氣顆粒污染的狀況及為城市的大氣污染防治提供理論基礎(chǔ)和科學依據(jù)。

1 材料與方法

1.1 研究區(qū)域概況

研究區(qū)位于山西省的省會太原市（111°30'～113°09'E，37°27'～38°25'N），平均海拔約 800 m，西、北、東三面環(huán)山，中、南部為河谷平原，整個地形北高南低。中、南部為河谷平原，整個地形北高南低．研究區(qū)域?qū)贉貛Ъ撅L性氣候，室外大氣壓年均920 hPa，年均降雨量456 mm，空氣濕度年均49%，室外近地面風速年均0.3 m?s-1，年平均氣溫9.5 ℃，1月最冷，平均氣溫 6.8 ℃；7月最熱，平均氣溫23.5 ℃。

1.2 監(jiān)測數(shù)據(jù)的采集及獲取

太原市城區(qū)各監(jiān)測點分布如圖1所示。各站點采樣口距地面高度為 10.5 m（塢城采樣點為 13.5 m），樣品分析方法依據(jù)中國《環(huán)境空氣質(zhì)量標準》（GB3905—2012）（環(huán)境保護部，2012）。本研究數(shù)據(jù)包括2014年2月16日—2016年2月16日共計733 d的PM10和PM2.5數(shù)據(jù)，顆粒物監(jiān)測設(shè)備24 h連續(xù)自動監(jiān)測、自動儲存?zhèn)鬏敂?shù)據(jù)，監(jiān)測數(shù)據(jù)有多種顯示方式，如按半小時、1 h、8 h、12 h、24 h顯示，本研究數(shù)據(jù)顯示方式為 8 h，即每日不定時間隔8 h顯示顆粒物質(zhì)量濃度，每天3個有效監(jiān)測數(shù)據(jù)，數(shù)據(jù)記錄總數(shù)為21867條（有效監(jiān)測天數(shù)見表1）。研究數(shù)據(jù)來源于自山西省環(huán)境監(jiān)測站。

表1 各監(jiān)測點研究期間的缺失值統(tǒng)計Table 1 Summary statistics of missing values of ten monitoring sites

各監(jiān)測點有效監(jiān)測天數(shù)和缺失天數(shù)統(tǒng)計如表 1所示，研究期間內(nèi)各監(jiān)測點有效天數(shù)基本滿額，僅部分站點缺失 1～7 d的數(shù)據(jù)，缺失比例僅為0.14%～0.96%，對于各站點缺失值的處理采取簡單刪除個案的方法。

1.3 數(shù)據(jù)整理與分析

圖1 研究區(qū)域大氣監(jiān)測站點分布圖Fig. 1 Location of study area and distribution of monitoring sites

數(shù)據(jù)整理：首先對于部分站點缺失值進行簡單刪除個案。在匯總統(tǒng)計多時間尺度數(shù)據(jù)庫的構(gòu)建中，以各監(jiān)測站點為基準匯總統(tǒng)計，并對無效值進行簡單刪除處理，對同一時間段僅存?zhèn)€別站點數(shù)據(jù)的情況亦進行簡單刪除，基本保證多時間尺度下所有監(jiān)測點數(shù)據(jù)的完備性。分別按照季節(jié)尺度、周尺度和日尺度進行分析研究，其中，日時間尺度劃分基于 8 h采集頻率，定義 01:00—09:00為早（Morning），10:00—18:00為白天（Daytime），19:00—24:00 為夜（Night）。

數(shù)據(jù)分析：運用經(jīng)典描述性統(tǒng)計分析方法對PM10和PM2.5數(shù)據(jù)的基本特征進行統(tǒng)計分析；運用HeatMap分析各監(jiān)測點PM10和PM2.5質(zhì)量濃度的區(qū)域性和季節(jié)性變化特征；采用Morlat連續(xù)小波變換的功率譜和總體方差分析揭示PM10和PM2.5質(zhì)量濃度變化的周期結(jié)構(gòu)及其顯著性；基于可視化主成分分析PCA確定多時間尺度下的PM10和PM2.5的變化特征。

本文中，太原市區(qū)域內(nèi)各監(jiān)控點PM10和PM2.5濃度變化的 HeatMap熱圖譜分析采用 R語言“pheatmap”包pheatmap函數(shù)，熱圖中各監(jiān)測點間的的聚類劃分是基于各監(jiān)測點PM10和PM2.5質(zhì)量濃度的歐式距離來確定的，由 pheatmap函數(shù)自動完成；周期結(jié)構(gòu)及顯著性的小波功率譜和小波方差運用Matlab2012b軟件編程完成；太原市各監(jiān)測點大氣顆粒污染物參數(shù)的PCA分析采用R語言（R 3.2.3版本）prcomp（）函數(shù)完成，PCA排序圖由“bioplot”包 bioplot（）函數(shù)完成。可視化 PCA圖設(shè)置突出變量格式，其解讀規(guī)則為：變量箭頭長短代表對前兩軸的貢獻大小，即變量箭頭越長，影響越大；變量間的相關(guān)性可由變量間相互夾角大小衡量，夾角越小，相關(guān)性越強。

2 結(jié)果分析與討論

2.1 大氣顆粒物質(zhì)量濃度季節(jié)性變化的基本特征

表2所示為太原市10個監(jiān)測點PM10和PM2.5日質(zhì)量濃度變化。統(tǒng)計結(jié)果顯示，PM10日質(zhì)量濃度的最大值均出現(xiàn)在冬季（其值分別為 758.0 μg?m-3和 696.0 μg?m-3），PM2.5日質(zhì)量濃度的最大值出現(xiàn)的季節(jié)不同，2014年出現(xiàn)在秋季（其值為 682.0 μg?m-3），2015 年出現(xiàn)在冬季（其值為 572.0 μg?m-3）；PM10平均質(zhì)量濃度的最高值所處的季節(jié)相同，均為冬季[2014 年值為(96.1±64.3) μg?m-3，2015 年值為(91.5±95.2) μg?m-3]，均高于國家二級年平均質(zhì)量濃度標準（70 μg?m-3），PM2.5平均質(zhì)量濃度的最高值所處的季節(jié)同樣均為冬季[2014年值為(147.4±81.7)μg?m-3，2015 年值為(135.3±91.7) μg?m-3]。2014 年，PM10和PM2.5在夏季的平均值均相對較低，其余3個季節(jié)（春、秋、冬）均高于 70 μg?m-3；而 2015年，PM10和PM2.5在冬季的質(zhì)量濃度相對較高，其余3個季節(jié)（春、夏、秋）和年均質(zhì)量濃度均低于70 μg?m-3，表明太原市大氣顆粒污染整體情況雖比較嚴重，但2015年較2014年，大氣顆粒污染物造成的環(huán)境污染情況有所好轉(zhuǎn)。圖2所示為太原城區(qū)各監(jiān)測PM10和PM2.5月質(zhì)量濃度變化的熱圖譜分析結(jié)果，依據(jù)污染物24 h平均值劃分為5級，分別為優(yōu)（0～35 μg?m-3）、良（35～75 μg?m-3）、輕度污染（75～115 μg?m-3）、中度污染（115～150 μg?m-3）、重度污染（150～250 μg?m-3），嚴重污染（250 μg?m-3及以上）。由圖可知，尖草坪站點、桃園站點和太原站點具有相似性，在地理分布上屬于太原市區(qū)北中部，是典型的傳統(tǒng)工業(yè)區(qū)（Group1）；小店站點、晉源站點、塢城站點具有相似性，在地理位置上處于太原南部，是近年來的經(jīng)濟發(fā)展的核心地帶（Group2）；位于太原市北部的上蘭站點、澗河站點、南寨站點和金勝站點相對經(jīng)濟發(fā)展較慢，處于過渡區(qū)（Group3）。整個監(jiān)測期間，各監(jiān)測點PM10和PM2.5污染程度在時間分布上存在明顯差異。在2014年2—5月（冬春季），2014年 10—12月（秋冬季），以及2015年11月—2016年1月（秋冬季），各監(jiān)測點PM10和PM2.5的污染情況相對嚴重；從空間上看，人類活動頻繁的新興經(jīng)濟區(qū)（塢城、晉源、小店）較其他監(jiān)測站點顆粒物污染嚴重。

表2和圖2的分析結(jié)果均表明，大氣顆粒污染物PM10和PM2.5在冬季均表現(xiàn)出污染程度相對嚴重的現(xiàn)象。造成冬季顆粒物污染嚴重的原因可能有兩方面，一是與太原市冬季供暖有關(guān)（太原市供暖法定時間：11月1日—次年3月31日），太原市的取暖主要靠化石燃料的燃燒（朱紅霞等，2014），是產(chǎn)生顆粒污染物的重要人為源，產(chǎn)生的顆粒物對大氣顆粒物污染貢獻值相對較大；二是與冬季的氣象因素有關(guān)，冬季樹葉凋零，植物對顆粒物的吸附作用大幅度降低，而且太原冬季降水少，氣候干燥多風，更易導致大氣顆粒污染問題。

表2 太原市大氣顆粒污染物質(zhì)量濃度季節(jié)性變化結(jié)果統(tǒng)計Table 2 Seasonal changes of particle matter concentrations μg?m-3

圖2 太原市各監(jiān)測點大氣顆粒污染物PM10（a）和PM2.5（b）月質(zhì)量濃度熱圖譜Fig. 2 Heat map of monthly changes in each monitoring site for PM10(a) and PM2.5(b)

圖2還揭示了在2014年期間，冬季污染程度比較嚴重，春、秋季節(jié)次之，而進入 2015年后，春、秋季的污染程度相對于 2014年同一時段有所減弱，冬季雖然污染程度比較嚴重，但相對于2014年亦有所好轉(zhuǎn)，造成這種現(xiàn)象的原因可能與近年來山西省加強工業(yè)污染的監(jiān)督和治理，積極調(diào)整能源產(chǎn)業(yè)結(jié)構(gòu)，加大城市綠化等一系列措施有關(guān)（馬麗梅等，2014）。

2.2 大氣顆粒物質(zhì)量濃度變化的周期結(jié)構(gòu)分析

2.2.1 顆粒物質(zhì)量濃度變化的小波功率譜和小波方差分析

采用小波分析方法對太原市大氣顆粒物質(zhì)量濃度變化的周期結(jié)構(gòu)進行分析，小波分析對于長時間序列周期結(jié)構(gòu)的識別具有明顯優(yōu)勢，在小波變換中，小波功率直接由其變換數(shù)值絕對值平方獲得（Wang et al.，2012）。小波變換數(shù)值絕對值的平方顯示了不同尺度或周期結(jié)構(gòu)下大氣顆粒污染物質(zhì)量濃度變化功率的相對強弱。

在監(jiān)測期間，太原市城區(qū)大氣顆粒物（PM10和 PM2.5）日質(zhì)量濃度變化的小波功率見圖 3。圖3a和圖3c所示為小波變化的實際振幅及其變化量級。小波功率譜圖揭示了城區(qū)顆粒物質(zhì)量濃度在不同時間尺度下的最高能量變化特征，PM10和PM2.5的時間序列均顯示相對較高能量的變化周期，以4～8 d的周期為主。精確的小波總體方差變化均揭示了 PM10和 PM2.5序列約 4～8 d的周期結(jié)構(gòu)大于95%的統(tǒng)計置信水平（圖3b和圖3d）。此外，小波總體方差還揭示了PM10和PM2.5時間序列在8～16 d的周期結(jié)構(gòu)均具有接近95%置信水平的特征。

通常，諸如小波模、小波實部分析中得出的周期結(jié)構(gòu)并無嚴格的物理意義。本文中Morlet小波功率譜及小波方差分析得出的周期結(jié)構(gòu)所代表的物理意義是：該周期結(jié)構(gòu)在所研究的時間段內(nèi)具有顯著的能量變化，代表了該時間段內(nèi)大氣顆粒污染物（PM10和PM2.5）的實際波動幅度。城市大氣顆粒污染物（PM10和PM2.5）變化出現(xiàn)4～8 d的周期性特征可能與城市人類生產(chǎn)生活的周期性有關(guān)。眾所周知，人類生產(chǎn)生活活動具有周變化特征（weekly cycles）（Sanchez-Lorenzo et al.，2012；Yang et al.，2016），即周末城市人類生產(chǎn)活動較工作日有所減少，在工業(yè)城市尤為如此，而這種周變化是導致大氣等因子呈現(xiàn)周期變化的重要因素，同時這種結(jié)構(gòu)還是城市熱源和空氣污染的重要影響因素（Earl et al.，2016；Javed et al.，2015）。

2.2.2 大氣顆粒污染物周波動結(jié)構(gòu)分析

為了進一步驗證大氣顆粒污染物（PM10和PM2.5）與城市人類生產(chǎn)生活的周期性關(guān)系，以周為時間尺度分別繪制PM10和PM2.5的熱圖譜（圖4）。如圖所示，PM10和PM2.5出現(xiàn)了相似的聚類劃分，即星期日被單獨聚為一類，周三和周四被劃分為了一類，周一、周二被劃分為一類，可見污染物在一周的變化中有明顯的特征。此外，在2014年3—5月（春），2014年10月—2015年1月，以及2015年 12月這段時間，PM10和 PM2.5污染程度在一周內(nèi)均維持較高值，其原因可能是這些月份均處于大氣顆粒污染物高發(fā)季和易發(fā)季，本身污染情況比較嚴重，不易區(qū)分一周內(nèi)的污染變化。除這幾個月之外的時間里，可以看出同一周時間內(nèi)，周六和周日的污染情況在大部分時候比其他幾天污染程度輕。聚類結(jié)果揭示的顆粒物周內(nèi)變化特征和一周內(nèi)周末污染情況相對較輕微的現(xiàn)象，從某種程度上驗證大氣顆粒污染變化與人們的日常生產(chǎn)生活有著密切的關(guān)系。

圖3 太原市大氣顆粒污染物PM10和PM2.5質(zhì)量濃度的小波功率譜和小波方差Fig. 3 The wavelet power spectrum PM10(a) and PM2.5(c)功率譜圖中，高分辨率區(qū)域是小波影響錐，其中黑色實線輪廓代表白噪背景下 95%的顯著性置信水平；小波總體方差中，藍色實線代表小波總體方差，藍色虛線是與功率譜一致背景下95%的顯著性置信水平The relative low-resolution region is the cone of influence, where zero padding has reduced the variance. Black contour is the 5% significance level, using a white-noise background spectrum. The global wavelet variance (solid line) for PM10(b) and PM2.5(d), and the dashed line is the confidence for the global wavelet variance, assuming the same confidence level and background spectrum as in power spectrum

圖4 太原市各監(jiān)測點大氣顆粒污染物PM10（a）和PM2.5（b）周變化熱圖譜Fig. 4 Heat map of weekly changes in each monitoring site for PM10(a) and PM2.5(b)

2.3 大氣顆粒物不同時間尺度的主成分分析（PCA）

在不同時間尺度下，應用 PCA對太原市各監(jiān)測點大氣顆粒污染物監(jiān)測樣本數(shù)據(jù)進行 PCA分析之前，先進行消除趨勢的對應分析（Detrended Correspondence Analysis，DCA）（Wang et al.，2014）。DCA結(jié)果顯示，季節(jié)時間尺度、周時間尺度和日時間尺度下其排序前2個軸中最大值均小于2個單位（月尺度變化未能符合主成分分析前提），表明各監(jiān)測點大氣顆粒污染物PM10和PM2.5監(jiān)測樣本滿足季節(jié)時間尺度、周時間尺度和日時間尺度下 PCA分析的前提條件。

由表3可知，太原市大氣顆粒污染物監(jiān)測樣本在不同時間尺度PCA的前2軸方差累計百分比均在70%左右，即前2個主成分已經(jīng)涵蓋了各影響變量所涵蓋的絕大部分信息。以季節(jié)作為劃分依據(jù)時，PM10前 2個軸累計捕捉方差為 67.8%，PM2.5前2個軸累計捕捉方差為74.6%；以星期作為劃分依據(jù)時，PM10前2個軸累計捕捉方差為87.4%，PM2.5前2個軸累計捕捉方差為84.9%；以天作為劃分依據(jù)時，PM10前2個軸累計捕捉方差為88.9%，PM2.5前2個軸累計捕捉方差為90.2%。

表3 太原市大氣顆粒污染物在不同時段主成分分析（PCA）前三軸的統(tǒng)計特征Table 3 The statistics of the previous three axis of PCA for PM10and PM2.5in different time scales

圖5所示為太原市各監(jiān)測點大氣顆粒物質(zhì)量濃度不同時間尺度下的 PCA分析。傳統(tǒng)工業(yè)區(qū)、新興經(jīng)濟發(fā)展區(qū)和過渡區(qū)分別用 Group1、Group2和Group3表示，其劃分依據(jù)聚類結(jié)果（圖 2）。圖 5中各組線圈包含了各站點96%以上的數(shù)據(jù)信息。圖5（a）和5（b）是以季節(jié)為劃分依據(jù)的大氣顆粒污染物參數(shù)的PCA分析，圖5（c）和5（d）是以星期為劃分依據(jù)的大氣顆粒污染物參數(shù)的PCA分析，圖5（e）和5（f）是以星期為劃分依據(jù)的大氣顆粒污染物參數(shù)的 PCA分析。以季節(jié)作為時間尺度對PM10和PM2.5進行PCA分析時，發(fā)現(xiàn)影響PM10和PM2.5季節(jié)性變化的指標既有共性也存在差異性，冬季指標箭頭最長，而春、夏、冬季指標箭頭相對較短且長度相當，表明影響太原市大氣 PM10和PM2.5的季節(jié)性變化受冬季影響最大。PCA結(jié)果還顯示，不同站點PM10和PM2.5的季節(jié)性變化均存在差別。此外，PM10和PM2.5受季節(jié)性影響也存在差異，秋季PM10與PC1軸有較強的相關(guān)性（圖5a），春季PM2.5與PC1軸有較強的相關(guān)性（圖5b）。綜上，季節(jié)尺度 PCA結(jié)果揭示了冬季是影響顆粒物季節(jié)變化的主要因素。這與前人對大氣顆粒物季節(jié)性變化的結(jié)果相一致（Lai，2016；王占山等，2015）。

以周為時間尺度對PM10和PM2.5進行PCA分析，結(jié)果揭示了影響PM10和PM2.5質(zhì)量濃度變化的周時間指標具有相似性。結(jié)果均顯示，周三（Wednesday）指標箭頭最長（圖5c和圖5d），表明一周之內(nèi)太原市城區(qū)顆粒物污染受周三影響最大，周五和周六次之，各變量之間相關(guān)性揭示周二和周四在PM10和PM2.5中具有一致的聚類特征。綜上，以周為時間尺度的 PCA結(jié)果揭示了周三是影響PM10和PM2.5周變化特征的主要因素。Cerveny et al.（1998）指出區(qū)域污染和人類活動有較強關(guān)系，表現(xiàn)出周循環(huán)（weekly cycle）的顯著特征。Shi et al.（2014）同樣指出北京地區(qū)綜合大氣污染指數(shù)具有周波動特征。世界其他城市大氣顆粒物的質(zhì)量濃度變化也被觀測到具有周波動特征（Delene et al.，2002；Gong et al.，2007；Jin et al.，2005），這種周波動特征與周期性的人類經(jīng)濟活動，尤其是明顯的交通特征有密切關(guān)系。

以天為時間尺度對PM10和PM2.5進行PCA分析，結(jié)果揭示了影響PM10和PM2.5質(zhì)量濃度變化的日時間區(qū)段差異性相對較大。PCA結(jié)果表明，PM10污染受夜晚影響最大（圖5e），PM2.5污染受早晨影響最大（圖5f），但白天（Daytime）太原各站點PM10和PM2.5在PCA分析中均與PC1軸有最大相關(guān)性，表明白天大氣顆粒物排放是造成夜間 PM10污染和早晨 PM2.5污染最明顯的主要原因。一天之內(nèi)不同時間區(qū)間下的 PCA結(jié)果揭示了白天顆粒物質(zhì)量濃度變化是造成 PM10和 PM2.5日變化特征的主要因素。王占山等（2013）研究表明白天時段顆粒物排放處于平緩的低谷區(qū)，大部分站點在8:00或21:00時出現(xiàn)峰值（王占山等，2013），且具有空間差異性。本研究表明，日尺度內(nèi)不同時間段PM10和PM2.5污染空間差異不大，但隨著時間尺度增大，空間變異性也增大，如圖5a和圖5b顯示，季節(jié)性變化中過渡區(qū) Group3（藍色圈）相比新興經(jīng)濟發(fā)展區(qū)（Group2）（綠色圈）范圍更大，顯然其受各時間參數(shù)的影響相對較弱。

圖5 監(jiān)測點大氣顆粒物PM10（a，c，e）和PM2.5（b，d，f）季節(jié)尺度、周尺度和日尺度PCA分析Fig. 5 Analysis of PCA for PM10(a, c, e) and PM2.5(b, d, f) on seasonally, weekly and daily scales

3 結(jié)論

（1）太原市大氣顆粒污染物PM10和PM2.5質(zhì)量濃度的變化存在明顯的季節(jié)差異。冬季顆粒物污染程度相對嚴重；2015年顆粒物污染相對于2014年同期有明顯好轉(zhuǎn)。

（2）太原市大氣顆粒污染物PM10和PM2.5質(zhì)量濃度的變化存在明顯的空間差異，新型經(jīng)濟發(fā)展區(qū)較傳統(tǒng)老工業(yè)區(qū)污染嚴重。太原市大氣顆粒污染物的濃度與城市生產(chǎn)生活的周波動變化相一致，存在周時間尺度的變化規(guī)律，PM10和PM2.5質(zhì)量濃度最大值出現(xiàn)在一周的中間時段（周三和周四），最小值出現(xiàn)在周末（周日）。

（3）PCA結(jié)果表明：太原市大氣顆粒物 PM10和 PM2.5季節(jié)性變化受冬季影響最大；一周之內(nèi)受周三影響最大；一天之內(nèi)PM10和PM2.5質(zhì)量濃度分別受夜晚和早晨影響最大，但白天顆粒物質(zhì)量濃度變化是造成其日變化特征的主要因素。

CERVENY R S, BALLING R C. 1998. Weekly cycles of air pollutants,precipitation and tropical cyclones in the coastal NW Atlantic region[J]. Nature, 394(6693): 561-563.

DERGHAM M, LEPERS C, VERDIN A, et al. 2015. Temporal-spatial variations of the physicochemical characteristics of air pollution Particulate Matter (PM2.5～0.3) and toxicological effects in human bronchial epithelial cells (BEAS-2B) [J]. Environmental Research,137: 256-267.

DELENE D J, OGREN J A. 2002. Variability of aerosol optical properties at four North American surface monitoring sites [J]. Journal of the Atmospheric Sciences, 59(6): 1135-1150.

EARL N, SIMMONDS I, TAPPER N. 2016. Weekly cycles in peak time temperatures and urban heat island intensity [J]. Environmental Research Letters, 11(7): 074003.

EEFTENS M, PHULERIA H C, MEIER R, et al. 2015. Spatial and temporal variability of ultrafine particles, NO2, PM2.5, PM2.5absorbance, PM10and PMcoarse in Swiss study areas [J]. Atmospheric Environment, 111:60-70.

GONG D Y, HO C H, CHEN D, et al. 2007. Weekly cycle of aerosolmeteorology interaction over China [J]. Journal of Geophysical Research: Atmospheres, 112(D22202): 1-9.

JAVED W, WEXLER A S, MURTAZA G., et al. 2015. Spatial, temporal and size distribution of particulate matter and its chemical constituents in Faisalabad, Pakistan [J]. Atmósfera, 28(2): 99-116.

JIN M, SHEPHERD J M, KING M D. 2005. Urban aerosols and their variations with clouds and rainfall: A case study for New York and Houston [J]. Journal of Geophysical Research: Atmospheres,110(D10S20): 1-12.

LAI L W. 2016. Public health risks of prolonged fine particle events associated with stagnation and air quality index based on fine particle matter with a diameter ＜2.5 μm in the Kaoping region of Taiwan [J].International Journal of Biometeorology, 60(12): 1907-1917.

LI H, WANG J, WANG Q G, et al. 2015. Chemical fractionation of arsenic and heavy metals in fine particle matter and its implications for risk assessment: A case study in Nanjing, China [J]. Atmospheric Environment, 103: 339-346.

LIU J, CUI S. 2014. Meteorological influences on seasonal variation of fine particulate matter in cities over Southern Ontario, Canada [J].Advances in Meteorology, (1): 1-15.

SANCHEZ-LORENZO A, LAUX P, FRANSSEN H J H, et al. 2012.Assessing large-scale weekly cycles in meteorological variables: a review [J]. Atmospheric Chemistry and Physics, 12(13): 5755-5771.

SHI M, WU H, ZHANG S, et al. 2014. Weekly cycle of magnetic characteristics of the daily PM2.5and PM2.5～10in Beijing, China [J].Atmospheric Environment, 98: 357-367.

WANG F, WANG X, ZHAO Y, et al. 2012. Nutrient response to periodic hydrological fluctuations in a recharging lake: a case study of Lake Baiyangdian [J]. Fresenius Environmental Bulletin, 21(5a): 1254-1262.

WANG F, WANG X, ZHAO Y, et al. 2014. Temporal variations of NDVI and correlations between NDVI and hydro-climatological variables at Lake Baiyangdian, China [J]. International journal of biometeorology,58(7): 1531-1543.

YANG X, LI Z, LIU L, et al. 2016. Distinct weekly cycles of thunderstorms and a potential connection with aerosol type in China [J]. Geophysical Research Letters, 43(16): 8760-8768.

環(huán)境保護部, 國家質(zhì)量監(jiān)督檢驗檢疫管理總局. 2012. 環(huán)境空氣質(zhì)量標準(GB 3905—2012)[S]. 北京: 中國環(huán)境科學出版社.

雷瑜, 張小玲, 唐宜西, 等. 2015. 北京城區(qū)PM2.5及主要污染氣體“周末效應”和“假日效應”研究[J]. 環(huán)境科學學報, 35(5): 1520-1528.

李名升, 任曉霞, 于洋, 等. 2016. 中國大陸城市 PM2.5污染時空分布規(guī)律[J].中國環(huán)境科學, 36(3): 641-650.

李名升, 張建輝, 張殷俊, 等. 2013. 近10年中國大氣PM10污染時空格局演變[J].地理學報, 68(11): 1502-1510.

馬麗梅, 張曉. 2014. 區(qū)域大氣污染空間效應及產(chǎn)業(yè)結(jié)構(gòu)影響[J]. 中國人口資源與環(huán)境, 24(7): 157-164.

陶燕, 劉亞夢, 米生權(quán), 等. 2014. 大氣細顆粒物的污染特征及對人體健康的影響[J]. 環(huán)境科學學報, 34(3): 592-597.

王蕊, 蘇玉紅, 卓少杰, 等. 2015. 我國 10城市冬季大氣顆粒物中多環(huán)芳烴污染及呼吸暴露風險評價[J].生態(tài)毒理學報, 10(4): 96-104.

王占山, 李云婷, 陳添, 等. 2015. 2013年北京市PM2.5的時空分布[J]. 地理學報, 70(1): 110-120.

朱紅霞, 趙淑莉. 2014. 中國典型城市主要大氣污染物的濃度水平及分布的比較研究[J]. 生態(tài)環(huán)境學報, 23(5): 791-796.

趙穎, 王飛, 溫彥平. 2015. 煤炭工業(yè)城市近地面臭氧變化特征及其與氣象條件的關(guān)聯(lián)[J]. 環(huán)境化學, 34(5): 842-852.

Abstract: Daily time-series data of PM10and PM2.5on ten locations in Taiyuan City were employed to analyze spatiotemporal characteristics, as well as the synchronous variations with the anthropogenic activities. Main influence variables on different time-scales were identified based on PCA. The results indicated that the obvious spatiotemporal differences were observed on PM10and PM2.5changes. The relatively higher concentrations of PM2.5and PM10in the newly downtown zones, comparing with those in the traditional heavy energy-consumption zones. The particle matter concentrations in winter showed relatively higher values. The wavelet analysis revealed that the period structures of both PM10and PM2.5time series were 4～8 days (P＜0.05), which indicated the weekly cycle fluctuations synchronized with human activities. The higher values of PM10and PM2.5concentrations occurred in midweek, but lower values on the weekend. The PCA results revealed the main influence variables of winter on seasonal scale,Wednesday on weekly scale, and Night and Morning on daily scale for PM10and PM2.5,separately. The variations of PM10and PM2.5concentrations on daily scale can be explained by the changes occurred on Daytime scale because of the approximate horizontal variable showed in PCA result. The research demonstrated the variations of PM10and PM2.5on different time-scale in a heavy energy-consumption city, and which also provided a new perspective for effective prevention and control of atmospheric pollution, as well as the scientific basis of relative standards and specifications for government.

Key words: atmospheric particle matter; different time scales; periodic structure; PCA; energy-consumption city.

Spatial-temporal Variations of PM10and PM2.5on Different Time-scales in Taiyuan

WANG Fei1, HAN Yelin2, ZHAO Ying3*

1. School of Life Science, Shanxi University, Taiyuan 030006, China; 2. School of Physical Education, Shanxi University, Taiyuan 030006, China;3. Shanxi Academy of Environmental Sciences, Taiyuan 030024, China

10.16258/j.cnki.1674-5906.2017.09.010

X513

A

1674-5906（2017）09-1521-08

王飛, 韓業(yè)林, 趙穎. 2017. 太原市大氣顆粒物 PM10和 PM2.5多時間尺度變化規(guī)律研究[J]. 生態(tài)環(huán)境學報, 26(9):1521-1528.

WANG Fei, HAN Yelin, ZHAO Ying. 2017. Spatial-temporal variations of PM10and PM2.5on different time-scales in Taiyuan [J].Ecology and Environmental Sciences, 26(9): 1521-1528.

國家自然科學基金項目（41401020；41601202）

王飛（1980年生），男，副教授，博士，主要研究方向為環(huán)境資源與利用。E-mail: nemo@sxu.edu.cn*通信作者：趙穎（1982年生），女，高級工程師，博士，主要研究方向為生態(tài)環(huán)境保護。E-mail: shadowying210@163.com

2016-10-01