1961—2012年?yáng)|北地區(qū)霾的時(shí)空變化特征分析

時(shí)一文，溫 娜，*，姚俊英

(1.南京信息工程大學(xué) 大氣科學(xué)學(xué)院，南京 210044；2.黑龍江省氣象服務(wù)中心，哈爾濱 150036)

1961—2012年?yáng)|北地區(qū)霾的時(shí)空變化特征分析

時(shí)一文1，溫 娜1，*，姚俊英2

(1.南京信息工程大學(xué) 大氣科學(xué)學(xué)院，南京 210044；2.黑龍江省氣象服務(wù)中心，哈爾濱 150036)

利用1961—2012年?yáng)|北地區(qū)霾日數(shù)統(tǒng)計(jì)資料，對(duì)東北地區(qū)霾的時(shí)空分布特征進(jìn)行了詳細(xì)分析，并利用EOF分解對(duì)霾日數(shù)進(jìn)行了分型討論。結(jié)果表明：東北地區(qū)1961—2012年霾日數(shù)雖然整體上呈下降趨勢(shì)，但存在明顯的年代際變化特征，霾日數(shù)在20世紀(jì)60—70年代較多，80年代之后下降；霾日數(shù)的月、季分布總體是春季霾日數(shù)最多，冬季次之，隨后是夏季，秋季最少；EOF分析的第一模態(tài)主要表達(dá)了東北地區(qū)霾日數(shù)為80年代之前顯著上升，80年代開(kāi)始突然下降；1961—2012年大部分區(qū)域霾日數(shù)沒(méi)有增加，且其中有約50%的區(qū)域呈現(xiàn)減少的趨勢(shì)；霾日數(shù)的變化和分布與工農(nóng)業(yè)發(fā)展?fàn)顩r密切相關(guān)，呈現(xiàn)的減少趨勢(shì)是由于獨(dú)特的歷史造成的。

東北地區(qū)；霾；時(shí)空特征；EOF分解

0 引 言

霾是一種對(duì)視程造成障礙的天氣現(xiàn)象, 霾發(fā)生時(shí)，有大量極細(xì)微的干塵粒等均勻地浮游在空氣中, 并伴隨著空氣混濁，使水平能見(jiàn)度<10 km[1]。近年來(lái)，不少學(xué)者對(duì)霾的變化趨勢(shì)進(jìn)行過(guò)研究[2-11]，認(rèn)為全國(guó)霾日數(shù)不僅有明顯增加的趨勢(shì)，而且存在季節(jié)變化特征，同時(shí)認(rèn)為這種霾天氣的頻繁發(fā)生與經(jīng)濟(jì)活動(dòng)密切相關(guān)。在區(qū)域分布特征的研究方面[12-13]，認(rèn)為霾的高發(fā)區(qū)主要分布在珠三角地區(qū)、華北地區(qū)以及長(zhǎng)三角地區(qū)等地。胡亞旦等[14]研究認(rèn)為在距離相近的區(qū)域中，大中城市的霾天氣比鄉(xiāng)村明顯偏多，并且深度討論了風(fēng)、降水和季風(fēng)環(huán)流等大氣環(huán)境對(duì)霾的影響。此外，研究者普遍的共識(shí)是：霾大多來(lái)源于下墊面的排放，在一定時(shí)段內(nèi)，致使霾的顆粒物得以輸送、堆積、稀釋擴(kuò)散的天氣系統(tǒng)背景、大氣層結(jié)條件和有著明顯的季、月變化，這種變化大于地面顆粒物排放的變化。

鑒于在眾多關(guān)于霾的研究中，還沒(méi)有專(zhuān)門(mén)針對(duì)東北地區(qū)霾的研究。因此本文將針對(duì)1961—2012年?yáng)|北地區(qū)霾日數(shù)的時(shí)空變化特征進(jìn)行系統(tǒng)分析，并力求對(duì)其成因進(jìn)行初步探索。

1 資料和方法

1.1 資料

文中東北地區(qū)包括黑龍江省、遼寧省、吉林省以及內(nèi)蒙古東部(指赤峰市、通遼市、興安盟和呼倫貝爾市)。1961—2012年?yáng)|北地區(qū)87個(gè)國(guó)家地面觀測(cè)站的霾觀測(cè)資料，來(lái)自中國(guó)氣象局國(guó)家氣象信息中心，該資料已經(jīng)通過(guò)了國(guó)家氣象信息中心資料室的質(zhì)量檢驗(yàn)，質(zhì)量和完整性很高。

1.2 霾日數(shù)的統(tǒng)計(jì)

測(cè)站處有霾發(fā)生的當(dāng)天記錄成1個(gè)霾日，月、季、年霾日數(shù)依據(jù)日值記錄統(tǒng)計(jì)。

1.3 經(jīng)驗(yàn)正交函數(shù)分解(EOF分析)

經(jīng)驗(yàn)正交函數(shù)分解(EOF分析)[15]的基本原理是將包含p個(gè)空間點(diǎn)的場(chǎng)隨時(shí)間變化的變量場(chǎng)進(jìn)行分解。假設(shè)樣本容量為n的資料，則時(shí)空?qǐng)鲋忻恳粋€(gè)空間點(diǎn)i和每一個(gè)時(shí)間點(diǎn)j的數(shù)值(xij)mn可看作由p個(gè)空間函數(shù)vik和p個(gè)時(shí)間函數(shù)tkj(k=1,2,…,p)合成的組合函數(shù)。

EOF按如下矩陣形式展開(kāi)：

(1)

式中Lp k為空間載荷向量陣(LV)；Fk n為時(shí)間系數(shù)陣(PC)，Kk為xp n的相關(guān)系數(shù)陣Rpp的特征值。由于時(shí)間系數(shù)和空間向量的正交性質(zhì)，可以用前幾個(gè)方差較大的分量的時(shí)間函數(shù)和空間函數(shù)來(lái)反映場(chǎng)的主要特征，用其乘積作為原氣象要素場(chǎng)的估計(jì)。

1.4 滑動(dòng)t-檢驗(yàn)法

滑動(dòng)t-檢驗(yàn)法(Movingt-test technique)[16]是用于檢驗(yàn)兩個(gè)隨機(jī)樣本平均值的顯著性差異的方法，把要研究的隨機(jī)變量x分成兩個(gè)子集x1與x2。xi的平均值、方差和樣本長(zhǎng)度(i=1,2)分別用ui、和ni代表，其中ni的長(zhǎng)度需要人為地定義。

假設(shè)原假設(shè)H0:u1-u2=0,統(tǒng)計(jì)量定義為：

(2)

1.5 空間分布圖的制作方法

采用GIS進(jìn)行空間分布圖的制作，采用自然斷點(diǎn)分級(jí)法進(jìn)行分級(jí)。

(3)

式中A為數(shù)組(數(shù)組長(zhǎng)度為N)；meani-j為每個(gè)等級(jí)中的平均值。該方法可用 GIS 軟件實(shí)現(xiàn)。

2 東北地區(qū)霾的時(shí)空分布特征

2.1 霾的時(shí)間變化特征

2.1.1 霾的月、季變化特征

東北地區(qū)霾的月、季變化曲線見(jiàn)圖1。從月變化曲線可見(jiàn)一年中5月霾日數(shù)最多，4月和1月次多，9月最少，6月次少；從季節(jié)變化曲線可見(jiàn)，東北地區(qū)春季霾日數(shù)最多，冬季次之，隨后是夏季，秋季最少。這與全國(guó)霾日數(shù)季節(jié)變化特點(diǎn)(冬季多，夏季少，春季秋季居中)[2，17]稍有不同。分析原因可能與東北地區(qū)冬季處于供暖期，采暖設(shè)施的使用量增加，燃燒煤的煙塵排放至空中，且冬季上層空氣逆溫情況普遍，這種層結(jié)狀況不利于污染物的排放，加之冬季沒(méi)有雨水，有利于低層大氣中氣溶膠粒子的聚集，形成霾天氣。而春季的霾日數(shù)略多于冬季可能因?yàn)槎究諝庵械臒焿m容易積累至春季，而且到了春季雖然逆溫情況有所減弱，但大風(fēng)天氣增多，容易使地面灰塵等揚(yáng)起送到空中，導(dǎo)致沙塵等氣溶膠增加，致使霾日數(shù)增加。夏秋兩季雨水較多，對(duì)空氣中的氣溶膠和顆粒物起到?jīng)_刷作用，因此夏季開(kāi)始霾日數(shù)減少。

圖1 東北地區(qū)霾的月、季變化曲線Fig.1 Month season change curve of haze days of Northeast of China

2.1.2 霾日數(shù)的年變化特征

東北地區(qū)1961—2012年霾的年變化曲線見(jiàn)圖2。由圖2可見(jiàn)，整個(gè)東北地區(qū)1961—2012年霾日數(shù)整體上呈下降趨勢(shì)，平均下降16.8 d/10 a。東北地區(qū)總霾日數(shù)的年變化大致分為3個(gè)階段。第一個(gè)階段為20世紀(jì)60年代至70年代末的明顯上升期，60年代平均80 d/a，到了70年代升至120 d/a，在此期間，平均增加4～5 d/a。到了20世紀(jì)80年代初霾日數(shù)突然呈階梯式下降，然后維持在低位徘徊，為40 d/a左右的平臺(tái)期，到90年代時(shí)雖然略有升高但仍以低位震蕩為主。這個(gè)結(jié)果與丁一匯等[17]的研究結(jié)論比較一致。

圖2 東北地區(qū)霾的年變化曲線Fig.2 Annual change curve of haze days of Northeast of China

東北地區(qū)年霾日數(shù)的滑動(dòng)t檢驗(yàn)趨勢(shì)線見(jiàn)圖3。由圖3可見(jiàn)東北地區(qū)霾日數(shù)在1978年超過(guò)顯著性檢驗(yàn)的臨界值之后，在1981年達(dá)到明顯的低谷，可見(jiàn)在20世紀(jì)80年代初東北霾日數(shù)出現(xiàn)明顯的下降突變。

圖3 東北地區(qū)年總霾日數(shù)滑動(dòng)t突變檢驗(yàn)趨勢(shì)線Fig.3 Moving t-test method curve of annual haze days Northeast of China

2.2 霾的空間變化特征

2.2.1 霾日數(shù)的空間分布特征

東北地區(qū)1961—2012年總霾日數(shù)的空間分布見(jiàn)圖4。由圖4可見(jiàn)，東北大部分地區(qū)霾日數(shù)不足60 d，平均發(fā)生1次/a。其中有約1/2的區(qū)域總霾日數(shù)不足25 d，平均發(fā)生1次/(2～5)a。霾日數(shù)的高值區(qū)主要分布在遼寧省西部，以及黑龍江省齊齊哈爾市周邊。而遼寧省內(nèi)主要有兩個(gè)明顯的高值區(qū)，一個(gè)是以彰武為中心的遼寧省西部地區(qū)，另一個(gè)位于遼寧省南部的綏中。最大高值區(qū)達(dá)到400 d以上，比低值區(qū)高出幾十倍，這可能與遼寧省的工業(yè)發(fā)展程度明顯高于其它地區(qū)有關(guān)[17]。

圖4 東北地區(qū)1961—2012年總霾日數(shù)的空間分布Fig.4 Annual haze days spatial distribution of Northeast of China from 1961—2012

從東北地區(qū)1961—2012年霾日數(shù)變化氣候傾向率(圖5)可見(jiàn)：東北地區(qū)1961—2012年大部分區(qū)域霾日數(shù)沒(méi)有增加，且其中有約1/2的區(qū)域呈現(xiàn)減少趨勢(shì)，僅存在遼寧省小塊區(qū)域的增加現(xiàn)象，這與全國(guó)霾日數(shù)明顯增加的特征并不一致。這不僅與東北地區(qū)的經(jīng)濟(jì)和工業(yè)水平相對(duì)全國(guó)較為滯后有關(guān)，還由于新中國(guó)成立初期開(kāi)始至20世紀(jì)80年代之前，東北地區(qū)作為老工業(yè)基地進(jìn)行高度的工業(yè)活動(dòng)，以及大范圍的土地開(kāi)發(fā)，致使霾頻繁發(fā)生、持續(xù)增長(zhǎng)[18]。而80年代之后開(kāi)始強(qiáng)調(diào)工業(yè)結(jié)構(gòu)調(diào)整和環(huán)境治理改善，工業(yè)活動(dòng)顯著減少，環(huán)境治理初見(jiàn)成效，使得霾日數(shù)呈現(xiàn)減少趨勢(shì)[5]。由此可見(jiàn)，由于東北地區(qū)獨(dú)特的歷史和地理因素，霾的整體變化情況主要受到工農(nóng)業(yè)發(fā)展的影響。

圖5 東北地區(qū)1961—2012年霾日數(shù)氣候傾向率空間分布圖Fig.5 Spatial distribution of annual haze days climate tendency trend of Northeast of China from 1961—2012

2.2.2 四季霾日數(shù)的空間分布特征

從東北地區(qū)四季霾日數(shù)空間分布圖(圖6)可見(jiàn)，春季高值區(qū)主要為東北地區(qū)南部的遼寧省的工業(yè)集中區(qū)，以及黑龍江省西部的齊齊哈爾市，最大值出現(xiàn)在遼寧省南部；與春季霾區(qū)域分布明顯不同，夏季高值區(qū)主要分布在東北地區(qū)北部，這個(gè)高值區(qū)面積廣且數(shù)值較高，最大值出現(xiàn)在黑龍江省西部；秋季和冬季的霾日數(shù)高值區(qū)的分布相似，僅有遼寧省西北部的工業(yè)發(fā)展集中區(qū)一個(gè)高值區(qū)，但冬季高值區(qū)的范圍和數(shù)值都略大于秋季?？v觀四季分布圖，在東北地區(qū)北部由漠河延伸到齊齊哈爾的區(qū)域，存在明顯的異?，F(xiàn)象：夏季總霾日數(shù)約為40 d，顯著高于其他季節(jié)該區(qū)域的霾日數(shù)，也明顯大于周?chē)鷧^(qū)域，這與夏季霾日數(shù)少于春季的總體分布特征有所不同，其原因有待于進(jìn)一步研究。整體來(lái)看，霾的四季分布情況除夏季外，高值區(qū)均集中在東北地區(qū)南部的工業(yè)發(fā)展集中區(qū)。

3 基于EOF分解的時(shí)空變化特征分析

通過(guò)對(duì)東北87站 1961—2012年逐年霾日數(shù)進(jìn)行EOF分析，得到第一、二、三模態(tài)的特征向量場(chǎng)，發(fā)現(xiàn)其分別存在綏中、彰武和莊河3個(gè)明顯的高值區(qū)，分別提取這3個(gè)站的年變化折線圖發(fā)現(xiàn)其與東北87個(gè)站的霾日數(shù)EOF分析第一、二、三模態(tài)的時(shí)間序列一致。說(shuō)明這3個(gè)模態(tài)分別反映了綏中、彰武和莊河3個(gè)測(cè)站的變化情況。因此去掉綏中、彰武和莊河3個(gè)測(cè)站的資料后，用EOF方法進(jìn)行分析，得到各個(gè)EOF模態(tài)，表1顯示了前5個(gè)主成分的特征值和解釋方差。前3個(gè)EOF模態(tài)的解釋方差之和約占52%，因而本文主要研究前3個(gè)主成分所對(duì)應(yīng)的模態(tài)，其特征場(chǎng)見(jiàn)圖7。前3個(gè)時(shí)間系數(shù)所對(duì)應(yīng)的演變曲線見(jiàn)圖8。

第一模態(tài)的空間場(chǎng)表現(xiàn)為全區(qū)一致的形式(圖7(a))，且解釋方差明顯高于其余的模態(tài)，因此是東北地區(qū)霾日數(shù)變化的最主要模態(tài)。由圖8(a)可見(jiàn)，第一模態(tài)的時(shí)間系數(shù)在60年代初開(kāi)始顯著增加，一直延續(xù)到70年代末，峰值達(dá)到24。而后開(kāi)始急劇下降，隨后一致維持在-5左右的平穩(wěn)狀態(tài)。在2003年經(jīng)歷了一次上升，隨后迅速回落至-5。結(jié)合第一模態(tài)空間場(chǎng)(圖7(a))來(lái)看，東北地區(qū)的霾日變化的第一模態(tài)為80年代之前顯著上升，80年代開(kāi)始突然下降，這也與圖2的年總霾日的變化相似。圖7(a)還顯示，大面積的正值高值區(qū)集中在黑龍江省西部，其余的幾個(gè)高值區(qū)比較分散。

(a)春季(3～5月)

(b)夏季(6～8月)

(c)秋季(9～11月)

(d)冬季(12～2月)

Table 1 Variance contribution and accumulated variance contribution of the first five EOF modes

序號(hào)特征值方差貢獻(xiàn)率/%累計(jì)方差貢獻(xiàn)率/%12555.190920.269370.2693721228.080690.129470.3988431144.329220.120640.519484736.785950.077670.597155514.773320.054270.65142

第二模態(tài)的空間場(chǎng)是南北相反的變化情況，南負(fù)北正，正的高值中心位于黑龍江省西部。其時(shí)間序列大體上呈上升趨勢(shì)，1975年之前為負(fù)值，1975年之后開(kāi)始為正值。即1975年之前東北地區(qū)大體上北部減少南部增多，1975年之后北部增多南部減少。第三模態(tài)的空間場(chǎng)(圖7(c))整體為南正北負(fù)，結(jié)合時(shí)間序列(圖8(c))可見(jiàn)，1970—1980年以及2000年前后，經(jīng)歷了東北地區(qū)南部增加、北部減少的過(guò)程。

圖7 EOF特征向量場(chǎng)Fig.7 The first three EOF modes

圖8 1961—2012東北地區(qū)霾日數(shù)EOF展開(kāi)的時(shí)間系數(shù)變化Fig.8 Variations of time coefficients of the first three EOF modes

4 結(jié) 論

1)東北地區(qū)1961—2012年霾日數(shù)雖然整體上呈下降趨勢(shì)，但存在明顯的年代際變化特征。20世紀(jì)60年代至70年代末為明顯上升期，20世紀(jì)80年代初霾日數(shù)出現(xiàn)明顯的下降突變。通過(guò)對(duì)1961—2012年逐年?yáng)|北霾日數(shù)進(jìn)行EOF分析也得到相似的結(jié)論。這種變化與霾的發(fā)生易受到人類(lèi)活動(dòng)尤其是工農(nóng)業(yè)活動(dòng)的影響有關(guān)， 20世紀(jì)60—70年代東北地區(qū)老工業(yè)基地建設(shè)和集中土地開(kāi)發(fā)的歷史背景導(dǎo)致霾出現(xiàn)明顯上升；80年代東北地區(qū)開(kāi)始產(chǎn)業(yè)結(jié)構(gòu)調(diào)整與東北老工業(yè)基地改造，生態(tài)環(huán)境得到改善，霾日數(shù)也隨之下降。

2)東北地區(qū)一年中春季霾日數(shù)最多，冬季次之，隨后是夏季，秋季最少。霾的季節(jié)差異可能受到大氣環(huán)流條件和污染源排放的影響，而霾的出現(xiàn)對(duì)于污染物排放的響應(yīng)有一定的滯后性。霾的四季分布情況除夏季外，高值區(qū)均集中在東北地區(qū)南部的工業(yè)發(fā)展集中區(qū)。春季高值區(qū)主要為東北地區(qū)南部的遼寧省的工業(yè)集中區(qū)，以及黑龍江省西部的齊齊哈爾市;夏季高值區(qū)主要分布在東北地區(qū)北部；秋季和冬季僅有遼寧省西北部的工業(yè)發(fā)展集中區(qū)一個(gè)高值區(qū)，但冬季高值區(qū)的范圍和數(shù)值都略大于秋季。

3)東北大部分地區(qū)霾天氣平均1次/a，高值區(qū)主要分布在遼寧省西南部，以及黑龍江省齊齊哈爾市。遼寧省內(nèi)主要有兩個(gè)明顯的高值區(qū)，一個(gè)是以彰武為中心的遼寧省西部地區(qū)，另一個(gè)是位于遼寧省南部的綏中。霾的重污染區(qū)大都在工業(yè)城市以及人口較為密集的城市。

4)從空間分布來(lái)看， 1961—2012年?yáng)|北大部分區(qū)域霾日數(shù)沒(méi)有增加，且其中有約1/2的區(qū)域呈現(xiàn)減少趨勢(shì)。再一次說(shuō)明霾日數(shù)的變化和分布與工農(nóng)業(yè)發(fā)展?fàn)顩r相關(guān)。

[1]中國(guó)氣象局．地面氣象觀測(cè)規(guī)范[S]．北京：氣象出版社，2003:23．

[2]高 歌．1961—2005年中國(guó)霾日氣候特征及變化分析[J]．地理學(xué)報(bào)，2008，63(7)：762-768．

[3]吳 兌．關(guān)于霾與霧的區(qū)別和灰霾天氣預(yù)警的討論[J]．氣象，2005，31(4)：3-7．

[4]吳 兌．再論都市霾與霧的區(qū)別[J]．氣象，2006，32(4)：9-14．

[5]吳 兌,吳曉京,李 菲,等．1951—2005 年中國(guó)大陸霾的時(shí)空變化[J]．氣象學(xué)報(bào)， 2010，68(5)：680-688．

[6]王業(yè)宏,盛春巖,楊曉霞,等．山東省霾日時(shí)空變化特征及其與氣候要素的關(guān)系[J]．氣候變化研究進(jìn)展，2009，5(1)：24-28．

[7]廖玉芳,吳賢云,潘志祥,等.1961—2006年湖南省霾現(xiàn)象的變化特征[J].氣候變化研究進(jìn)展，2007，3(5)：260-265．

[8]童堯青,銀 燕,錢(qián) 凌,等．南京地區(qū)霾天氣特征分析[J]．中國(guó)環(huán)境科學(xué)，2007，27(5)：584-588．

[9]顏嬌瓏．基于衛(wèi)星遙感資料對(duì)中國(guó)地區(qū)霾時(shí)空特征分析[D]．蘭州：蘭州大學(xué)，2014．

[10]宋連春，高 榮,李 瑩,等．1961—2012年中國(guó)冬半年霾日數(shù)的變化特征及氣候成因分析[J]．氣候變化研究進(jìn)展，2013，9(5)：313-318．

[11]潘本鋒，汪 巍,李 亮，等．我國(guó)大中型城市秋冬季節(jié)霧霾天氣污染特征與成因分析[J]．環(huán)境與持續(xù)發(fā)展，2013，(1):33-36．

[12]Zhang X Y,Wang Y Q,Niu T．et al．Atmospheric aerosol compositions in China: Spatial/temporal variability,chemical signature,regional haze distribution and comparisons with global aerosols[J]．Atmospheric Chemistry and Physics，2011，(11)：26 571-26 615．

[13]孫 彧,馬振峰,牛 濤,等．最近40年中國(guó)霧日數(shù)和霾日數(shù)的氣候變化特征[J]．氣候與環(huán)境研究，2013，18 (3)：397-406．

[14]胡亞旦,周自江．中國(guó)霾天氣的氣候特征分析[J]．氣象，2009，35(7)：73-78．

[15]黃嘉佑．氣象統(tǒng)計(jì)分析與預(yù)報(bào)方法[M]．2版．北京：氣象出版社，2000：135-145．

[16]符淙斌，王 強(qiáng)．氣候突變的定義和檢測(cè)方法[J]．大氣科學(xué)，1992，16(4)：482-493．

[17]丁一匯,柳艷菊．近 50 年我國(guó)霧和霾的長(zhǎng)期變化特征及其與大氣濕度的關(guān)系[J]．中國(guó)科學(xué)：地球科學(xué)， 2014，44(1)：37-48．

[18]薛繼坤．東北老工業(yè)基地形成演變過(guò)程的歷史分析[J]．長(zhǎng)春大學(xué)學(xué)報(bào)，2005，15(5): 8-11．

Spatial and temporal variations of the haze from 1961—2012 in Northeast of China

SHI Yi-Wen1，WEN Na1，*，YAO Jun-Ying2

(1.College of Atmospheric Science，Nanjing University of Information Science and Technology，Nanjing 210044,China；2.Heilongjiang Meteorological Service Center, Haerbin 150036，China)

Based on haze the data days of chinese northeast during 1961—2012, the characteristics and changes of temporal and spatial distributions of haze days are analysed, as well as the method of vector EOF expansion is used to carry on a classification discussion. The results show that although the haze days in northeast of China during 1961—2012 exhibits declining trends as a whole, but shows obvious decadal variation characteristics; There are more haze days during 1960s—1970s century, and then reduce significantly 1980s; Haze day shows the character of more in winter half year and less in summer half year. The haze days are the most in spring, then followed by winter, summer, autumn; The first mode of EOF analysis expresses that the haze days in the northeast experiences a marked increase tendency in the days before 1980s, and then begins to fall suddenly; In most area of the northeast of China, the trend is not increasing in 1961—2012, and over 50% shows the tendency of decrease; The industrial and agricultural development are possible reasons of the variations, and the declining trend is due to the unique history background.

Northeast of China; haze; spatial and temporal variations; EOF expansion

10.13524/j.2095-008x.2015.04.057

2015-07-30

中國(guó)科技部軟科學(xué)研究計(jì)劃項(xiàng)目(2012GXS4B071)

時(shí)一文(1993-)，女，黑龍江哈爾濱人，本科，學(xué)士，研究方向：大氣科學(xué)研究，E-mail：ssussi831@163.com；*通訊作者：溫 娜(1979-)，女，河南南陽(yáng)人，講師，博士，研究方向：大氣科學(xué),E-mail:wenna@nuist.edu.cn。

P427.122

A

2095-008X(2015)04-0016-07