黃超，趙錦慧*，何超，劉玉青，李小莉，劉書凱

1. 湖北大學(xué)資源環(huán)境學(xué)院，湖北 武漢 430062；2. 武漢大學(xué)資源與環(huán)境科學(xué)學(xué)院，湖北 武漢 430072；3. 西南林業(yè)大學(xué)國(guó)家高原濕地研究中心，云南 昆明 650224

黑碳（Black Carbon，BC）是指懸浮在空氣中的黑色碳質(zhì)顆粒物，它是含碳物質(zhì)在不充分燃燒條件下形成的無(wú)定型碳質(zhì)（朱厚玲，2003）。黑碳?xì)馊苣z的排放主要分為自然源和人為源，其中自然源包括火山爆發(fā)、森林大火等；人為源主要是化石燃料和生物燃料的燃燒。自然源具有區(qū)域性和偶然性特征，而人為源則是長(zhǎng)期而持續(xù)的（婁淑娟等，2005）。

黑碳?xì)馊苣z作為大氣中最重要的氣溶膠之一，對(duì)太陽(yáng)輻射具有強(qiáng)烈的吸收作用（張華等，2009），能夠影響大氣環(huán)流、全球水分循環(huán)（董俊玲等，2010），破環(huán)區(qū)域大氣的穩(wěn)定性，直接或間接影響區(qū)域和全球的輻射效應(yīng)、氣候變化等。同時(shí)它能降低大氣能見度，攜帶致癌物質(zhì)，對(duì)人體健康造成危害（Kim et al.，2003），因此關(guān)于BC的研究受到越來(lái)越多的關(guān)注。近幾年來(lái)，Sharma et al.（2004）研究了加拿大北極地區(qū)黑碳濃度的長(zhǎng)期規(guī)律，Murphy et al.（2011）使用連續(xù)BC觀測(cè)數(shù)據(jù)分析了美國(guó)黑碳時(shí)間變化，發(fā)現(xiàn)1994—2004年BC濃度降低了超過(guò)25%，這些連續(xù)的BC觀測(cè)為美國(guó)和北極地區(qū)積累了大量的科研數(shù)據(jù)。而在中國(guó)，有關(guān)黑碳?xì)馊苣z的觀測(cè)研究也取得不少進(jìn)展，2009年劉新春等（2011）271-272對(duì)烏魯木齊冬季黑碳?xì)馊苣z進(jìn)行觀測(cè)研究，2010年Verma et al.（2010）研究了廣州市黑碳?xì)馊苣z的時(shí)間變化，2012年姚青等（2012）對(duì)天津城區(qū)秋冬季黑碳?xì)馊苣z進(jìn)行了觀測(cè)分析，上海、南京等城市的BC氣溶膠觀測(cè)研究也相繼被開展，這些研究揭示了不同城市黑碳?xì)馊苣z的變化特征，闡明了黑碳?xì)馊苣z的氣候效應(yīng)、環(huán)境效應(yīng)和生物地球化學(xué)效應(yīng)。

從黑碳?xì)馊苣z的全球分布來(lái)看，北半球的黑碳?xì)馊苣z濃度明顯高于南半球，而中國(guó)東部就是黑碳?xì)馊苣z濃度高值中心區(qū)之一（曹軍驥等，2011）。武漢市是中國(guó)重要排放城市之一，顆粒物減排存在一定壓力。其中，以黑碳排放量較高，對(duì)區(qū)域氣候和空氣質(zhì)量影響很大。因此，應(yīng)重點(diǎn)關(guān)注武漢地區(qū)黑碳?xì)馊苣z的分布及其環(huán)境效應(yīng)。

目前關(guān)于武漢地區(qū)BC的基礎(chǔ)觀測(cè)研究較少，主要集中在BC與有機(jī)碳、顆粒物綜合分析上。2011年張宇堯（2017）利用氣溶膠的復(fù)折射指數(shù)首次反演得到武漢冬季黑碳和有機(jī)碳柱質(zhì)量濃度；Wei et al.（2015）在2014年分析了PM2.5、PM10的污染特征及其與BC的關(guān)系。本研究利用2015年秋冬季黑碳?xì)馊苣z觀測(cè)數(shù)據(jù)和主要污染物數(shù)據(jù)，運(yùn)用 SPSS分析軟件以及軌跡模型等方法，研究黑碳?xì)馊苣z的變化特征及其與污染物的關(guān)系等，為武漢城區(qū)黑碳?xì)馊苣z觀測(cè)提供數(shù)據(jù)資料，有助于推動(dòng)中國(guó)黑碳?xì)馊苣z的系統(tǒng)觀測(cè)研究工作，為探索灰霾監(jiān)測(cè)評(píng)估方法積累一定的實(shí)踐經(jīng)驗(yàn)。

1　采樣和分析

選取武漢城區(qū)漢陽(yáng)公園、解放公園、中山公園、硚口公園、常青公園、黃鶴樓景區(qū)、東湖風(fēng)景區(qū)、沙湖公園、石門峰公園、南干渠游園共 10個(gè)景點(diǎn)作為黑碳?xì)馊苣z的采樣點(diǎn)，采樣點(diǎn)分布如圖1所示。

在公園綠地地勢(shì)開闊的平地上，將儀器架設(shè)在距離地面2 m高度處進(jìn)行采樣。各采樣點(diǎn)下墊面是植被覆蓋率較高的公園綠地，周圍緊鄰交通主干道、居民樓和辦公樓，周圍500 m范圍內(nèi)無(wú)工廠等大氣污染源。采樣時(shí)間為2015年9月—2016年2月。采樣觀測(cè)儀器為美國(guó)瑪基科技公司（Magee Scientific Co，USA）生產(chǎn)的AE-51微型黑碳儀，其工作原理是利用Aethalometer光學(xué)吸收方法對(duì)空氣中的黑碳?xì)馊苣z進(jìn)行實(shí)時(shí)采樣和分析，屬于光學(xué)灰度測(cè)量法。此次采樣以880 nm波長(zhǎng)發(fā)光二極管作為光源，測(cè)量范圍為 0～1 mg·m-3，測(cè)量分辨率為0.001 μg·m-3，測(cè)量精度為±0.100 μg·m-3·min，儀器流量為50 mL·min-1，每隔60 s獲取1組實(shí)時(shí)監(jiān)測(cè)的武漢市各個(gè)綠地公園上空的黑碳質(zhì)量濃度數(shù)據(jù)并自動(dòng)保存，每次持續(xù)觀測(cè)50 min。觀測(cè)期間共獲得了1 min數(shù)據(jù)12202組，保留了11822組，總體數(shù)據(jù)完整率達(dá)到96.89%（表1）。

本文分析了2015年9月—2016年2月為期6個(gè)月的黑碳質(zhì)量濃度觀測(cè)數(shù)據(jù)，同時(shí)從武漢市環(huán)保局網(wǎng)站（http://www.whepb.gov.cn/）搜集同期武漢市可吸入顆粒物（PM10）、細(xì)顆粒物（PM2.5）、二氧化氮（NO2）、二氧化硫（SO2）、一氧化碳（CO）、臭氧（O3）實(shí)時(shí)監(jiān)測(cè)的24 h日均值等污染物數(shù)據(jù)。觀測(cè)期內(nèi)，共獲得12202組1min黑碳?xì)馊苣z數(shù)據(jù)，對(duì)其中因儀器斷電和更換濾膜等原因產(chǎn)生的負(fù)值和異常數(shù)據(jù)進(jìn)行剔除，數(shù)據(jù)資料完整率達(dá)到96.89%（表1），主要污染物數(shù)據(jù)也進(jìn)行異常值剔除處理。

2　結(jié)果與討論

2.1　BC質(zhì)量濃度時(shí)間變化規(guī)律

2.1.1BC質(zhì)量濃度日變化規(guī)律

圖1　采樣點(diǎn)位置圖Fig. 1　The location map of sampling points

根據(jù)各采樣點(diǎn)BC氣溶膠質(zhì)量濃度的連續(xù)觀測(cè)數(shù)據(jù)，可以得到武漢市2015年9月—2016年2月BC氣溶膠質(zhì)量濃度的日變化，如圖2所示。結(jié)果表明：武漢市秋冬季黑碳質(zhì)量濃度日均值為(5.227±3.12) μg·m-3，范圍為 1.205～16.534 μg·m-3，約61%的日均值濃度數(shù)據(jù)分布在2～6 μg·m-3區(qū)間范圍內(nèi)，具有較好的集中分布趨勢(shì)，并且在月份分布上也呈現(xiàn)明顯的分異特征，其中，2016年1月BC質(zhì)量濃度最高，為(7.778±3.12) μg·m-3；2015 年 9月 BC 質(zhì)量濃度最低，為(2.767±3.12) μg·m-3。

表1　數(shù)據(jù)質(zhì)量控制、統(tǒng)計(jì)Table 1　The control and statistics of data quality

圖2　武漢市秋冬季BC濃度時(shí)間序列圖Fig. 2　Time series of BC in Wuhan from Sep 2015 to Feb 2016

2.1.2BC質(zhì)量濃度月變化規(guī)律

根據(jù)BC觀測(cè)值計(jì)算得到武漢市秋冬季各月黑碳平均質(zhì)量濃度，其中1月最高，為8.052 μg·m-3，11月最低，為 2.816 μg·m-3，表現(xiàn)為：1月＞12月＞2月＞10月＞9月＞11月，如圖3所示。

這個(gè)規(guī)律與冬季大氣BC質(zhì)量濃度高于秋季的一般規(guī)律相吻合，原因包括，（1）氣象因素。武漢市屬亞熱帶季風(fēng)氣候，冬季氣溫較低，空氣潮濕，經(jīng)常伴有霧產(chǎn)生，這種條件有利于近地面大氣層保持穩(wěn)定狀態(tài)，使逆溫強(qiáng)度增大，不利于氣體污染物在垂直和水平方向的擴(kuò)散，加重了顆粒物的積聚污染；同時(shí)秋季溫度較高，空氣對(duì)流運(yùn)動(dòng)增強(qiáng)，降水較多，對(duì)顆粒物的濕沉降作用大，因此相較秋季，冬季大氣BC質(zhì)量濃度較高。（2）人為因素。由于冬季是供暖季節(jié)，武漢雖然沒有北方大規(guī)模的采暖供暖，但是部分地區(qū)仍然存在一定形式的燃燒取暖。另外，冬季也是工廠燃煤的高峰期，對(duì)煤炭需求量大，這兩方面直接導(dǎo)致武漢冬季空氣中的顆粒物大大增加，從而造成BC濃度偏高。

圖3　武漢秋冬季BC濃度月變化Fig. 3　Monthly variations of BC concentration in Wuhan in autumn and winter

2.2　BC濃度的空間分布差異

將漢陽(yáng)公園、黃鶴樓景區(qū)、東湖風(fēng)景區(qū)等 10個(gè)城區(qū)公園綠地樣點(diǎn)BC質(zhì)量濃度數(shù)據(jù)進(jìn)行處理，得到武漢秋冬季BC質(zhì)量濃度空間分布，結(jié)果如圖4所示。

圖4　武漢秋冬季BC濃度空間分布Fig. 4　The spatial distribution of BC concentration in Wuhan in autumn and winter

各采樣點(diǎn)在秋冬季和各月份的空間分布上具有明顯的一致性，總體上，武漢秋冬季大氣BC質(zhì)量濃度最高值出現(xiàn)在石門峰公園（其中硚口公園11月BC數(shù)據(jù)部分缺失，不納入比較），最低出現(xiàn)在東湖風(fēng)景區(qū)，BC質(zhì)量濃度大小順序?yàn)椋菏T峰公園＞解放公園＞沙湖公園＞漢陽(yáng)公園＞南干渠游園＞常青公園＞黃鶴樓景區(qū)＞中山公園＞東湖風(fēng)景區(qū)。圖 5所示為各采樣點(diǎn)秋冬季BC質(zhì)量濃度與距主干道路距離的關(guān)系，計(jì)算得到相關(guān)系數(shù)大致都為-0.3，說(shuō)明在局部氣候并且區(qū)域特征較為一致的大背景下，導(dǎo)致空間分布差異的原因可能與距離車流較多的主路的相對(duì)位置有關(guān)，距離越近，BC質(zhì)量濃度越高。同時(shí)在采樣期間，石門峰公園附近青王路改造工程正在施工，運(yùn)輸車輛造成大量揚(yáng)塵飛起。石門峰公園位于三環(huán)線外，小作坊燃煤和室外秸稈焚燒等無(wú)組織排放源較多，因此空氣中揚(yáng)塵、顆粒物較多，BC濃度較高。而東湖風(fēng)景區(qū)位于空氣質(zhì)量較好的武昌區(qū)，環(huán)境監(jiān)管較嚴(yán)，燃煤和秸稈焚燒等現(xiàn)象基本不存在，同時(shí)采樣點(diǎn)附近的沿湖大道一直以來(lái)都被作為傳統(tǒng)的騎行和觀光道路，機(jī)動(dòng)車流量小，排放尾氣少。再加上廣闊水體和較高的植被覆蓋率，其下墊面狀況有利于空氣凈化，因而BC質(zhì)量濃度較低。

圖5　BC濃度與距主干道距離的關(guān)系Fig. 5　The relationship with BC concentration and the distance from the main road

2.3　武漢市BC濃度及與國(guó)內(nèi)外各大城市的比較

大氣成分的本底濃度指能夠反映某一尺度區(qū)域內(nèi)處于均勻混合狀態(tài)的某種成分的大氣濃度，而均勻混合狀態(tài)僅是一種理想狀態(tài)，因此，一般認(rèn)為所謂本底濃度也就是不受局部污染源直接影響的大氣成分的濃度（趙玉成等，2008）。采用最大頻數(shù)法將11822組數(shù)據(jù)進(jìn)行處理，獲得了BC本底濃度。從圖6可以看出，武漢市大氣BC濃度變化符合對(duì)數(shù)正態(tài)分布特征。最大出現(xiàn)頻次所對(duì)應(yīng)的濃度值（即對(duì)數(shù)正態(tài)分布的平均值）反映了該地區(qū)最具代表性大氣狀況下的BC濃度，可以代表該地區(qū)BC濃度的本底濃度。

如圖6所示，最大出現(xiàn)頻次所對(duì)應(yīng)的濃度值為5.45 μg·m-3。查閱文獻(xiàn)發(fā)現(xiàn)，國(guó)內(nèi)外重點(diǎn)針對(duì)秋冬季BC氣溶膠質(zhì)量濃度的研究并不多，大多數(shù)關(guān)于城市BC觀測(cè)研究的時(shí)間跨度比較長(zhǎng)，因此將具有一定代表性的城市黑碳觀測(cè)數(shù)據(jù)進(jìn)行整理，并與武漢秋冬季BC本底濃度進(jìn)行對(duì)比分析，結(jié)果如表2所示。

圖6　BC濃度頻數(shù)分布Fig. 6　The frequency distribution of BC concentration

表2　武漢秋冬季BC本底濃度與其它城市的對(duì)比Table 2　Comparisons in BC background concentration between Wuhan and other cities

可以發(fā)現(xiàn)，武漢秋冬季BC濃度在所有城市里屬于中等水平，高于西寧、南極、瓦里關(guān)等環(huán)境較好的偏遠(yuǎn)地區(qū)和天津、芝加哥等經(jīng)濟(jì)比較發(fā)達(dá)的城市，低于烏魯木齊、西安等煤炭、燃料消耗較多的城市。

2.4　BC與6類污染物的同源性分析

為了證實(shí)BC與SO2、NOx、O3之類的反應(yīng)性氣體是否存在同源性，本文選取PM2.5、PM10、SO2、NO2、O3、CO等 6種大氣污染物作為研究對(duì)象，與BC進(jìn)行對(duì)應(yīng)分析。

2.4.1相關(guān)性分析

將2015年9月—2016年2月武漢秋冬季BC濃度觀測(cè)值與同期武漢環(huán)保局發(fā)布的主要污染物數(shù)據(jù)進(jìn)行相關(guān)性分析（表3），從相關(guān)系數(shù)矩陣可以看出：BC與 PM2.5、CO呈高度正相關(guān)，與 PM10呈中度正相關(guān)，與 O3呈中度負(fù)相關(guān)，在顯著水平0.05上與 SO2、NO2呈正相關(guān)，相關(guān)系數(shù)分別為0.885、0.8、0.658、-0.436、0.262、0.259。

武漢市BC氣溶膠與PM2.5、PM10、CO的高相關(guān)性說(shuō)明，BC與三者具有同源性。主要來(lái)自化石燃料的燃燒（如機(jī)動(dòng)車尾氣、燃煤）、揮發(fā)性有機(jī)物（VOCs）等，污染源為工廠燃煤、汽車尾氣排放、漆料噴涂行業(yè)等。

2.4.2多元線性回歸分析

由于BC與PM2.5、PM10、CO、O3具有較高的相關(guān)性，以多個(gè)自變量的最優(yōu)組合預(yù)測(cè)或估計(jì)因變量，比只用1個(gè)自變量進(jìn)行預(yù)測(cè)或估計(jì)更有效，更符合實(shí)際，因此多元線性回歸的實(shí)用意義更大。本文以BC作為因變量，以PM2.5、PM10、CO、O3作為自變量，運(yùn)用SPSS多元線性回歸分析，探討它們之間的線性關(guān)系，建立回歸方程，得到模型系數(shù)，結(jié)果如表4所示。

方程中的常數(shù)項(xiàng)為-1136.564，偏回歸系數(shù)b1為4.216，b2為 66.316，b3為-7，b4為 40.324，經(jīng) t檢驗(yàn)，b1、b2、b3、b4的概率 P 值分別為 0.226、0.763，0.112、0.211。按照給定的顯著性水平0.10，均具有顯著性意義。同時(shí)VIF值均小于10，各自變量間的共線性不明顯。經(jīng)過(guò)上述檢驗(yàn)，BC與PM10、CO、O3間關(guān)系式可以表達(dá)為：

式中，Y為 BC 的質(zhì)量濃度；X1、X2、X3、X4分別為PM10、CO、O3、PM2.5的質(zhì)量濃度。上式解釋了BC與多個(gè)自變量之間的線性關(guān)系，為同源性分析提供統(tǒng)計(jì)學(xué)依據(jù)。

2.5　BC后向軌跡分析

為了分析大氣污染物輸送、擴(kuò)散軌跡，運(yùn)用HYSPLIT-4模型進(jìn)行后向軌跡分析。

2.5.1秋冬季節(jié)BC的軌跡分析

（1）秋季

如圖7所示，2015年秋季，100 m高度上起源于美國(guó)東海岸、北大西洋附近海域的氣團(tuán)向東經(jīng)過(guò)地中海、西亞、巴基斯坦、印度北部、西藏，到達(dá)中國(guó)青海地區(qū)，與500 m高度來(lái)自加拿大東北部、北太平洋附近海域的氣團(tuán)經(jīng)過(guò)復(fù)雜運(yùn)動(dòng)后一同匯合，而此時(shí)750 m高度上起源于加拿大北方、北冰洋附近海域的氣團(tuán)向東穿越北半球大部分區(qū)域后，也運(yùn)動(dòng)至此。最后3股氣流一起沿東南方向推進(jìn)經(jīng)過(guò)甘肅、陜西、河南等地到達(dá)武漢地區(qū)。

表3　BC與主要污染物之間的相關(guān)系數(shù)矩陣Table 3　The correlation coefficient matrix between BC and the main pollutant

表4　模型系數(shù)統(tǒng)計(jì)Table 4　The statistics of model coefficients

圖7　武漢市秋冬季節(jié)BC的軌跡分析Fig. 7　The back-trajectory analysis of BC for Wuhan City in autumn and winter

（2）冬季

2015年冬季，100、500、750 m高度上的3個(gè)氣團(tuán)均起源于俄羅斯遠(yuǎn)東和北冰洋附近地區(qū)，并且各自大概沿著北緯 45°緯線圈向東運(yùn)動(dòng)，其中在北極地區(qū)進(jìn)行多次不規(guī)則的重復(fù)運(yùn)動(dòng)，最后一起在俄羅斯北部的北冰洋海域匯合，然后一路南下，穿越俄羅斯、蒙古中部、內(nèi)蒙古、山西、河南等地，到達(dá)武漢地區(qū)。在冬季前期，空氣質(zhì)點(diǎn)的高度均在2000 m以下；到了中期，100 m高度的氣流的空氣質(zhì)點(diǎn)先上升到6000 m以上，后下降至2000 m以下，其他的兩股氣流則無(wú)明顯變化；而在后期，3股氣流的空氣質(zhì)點(diǎn)高度均上升至2000 m以上，變化趨勢(shì)基本一致。

由此可以看出，2015年秋冬季武漢地區(qū)污染物主要受來(lái)自西北和北部氣團(tuán)影響，經(jīng)過(guò)山西、河南等傳統(tǒng)污染較重地區(qū)，是武漢地區(qū)大氣污染物的重要來(lái)源。

2.5.2兩次BC 污染事件的氣團(tuán)軌跡分析

黑碳來(lái)源一般包括局地來(lái)源（城市機(jī)動(dòng)車、居民燃煤和烹飪）、周邊源（如農(nóng)村秸稈燃燒）以及本地源，李楊等（2005）在 2003年秋冬季西安市黑碳?xì)馊苣z觀測(cè)研究中發(fā)現(xiàn)：農(nóng)村秸稈燃燒等非局地的周邊源對(duì)西安秋冬季大氣中的BC貢獻(xiàn)份額達(dá)到1/3以上，超過(guò)機(jī)動(dòng)車尾氣和居民活動(dòng)等來(lái)源的貢獻(xiàn)。因此在不討論本地源和各種氣象條件的影響下，探討秋冬季武漢地區(qū)兩次黑碳污染事件的空氣氣團(tuán)72 h后向軌跡，以探清外來(lái)源的大致路徑，解釋污染原因。觀測(cè)發(fā)現(xiàn)，在BC濃度月均值最高的1月，發(fā)生了兩次較大的BC污染事件（2日和17日），BC日均值濃度分別達(dá)到12.588、11.581 μg·m-3。利用HYSPLIT4后向軌跡方法模擬2日和17日在100、500、1000 m 3個(gè)垂直高度層的空氣氣團(tuán)72 h回推軌跡，結(jié)果如圖8所示。

在2日和17日的72 h后向軌跡中，空氣氣團(tuán)均從較高的垂直層向低處集聚，造成武漢局部的低空污染。不同的是：2日的污染源主要來(lái)自東部和南部，途經(jīng)江西、湖南等地，氣流運(yùn)動(dòng)垂直高度??；17日的污染源來(lái)自北部的河南等傳統(tǒng)農(nóng)業(yè)地區(qū)，氣流運(yùn)動(dòng)垂直高度大，污染物擴(kuò)散快，因而其BC濃度比2日小。兩次BC污染事件可能由來(lái)自不同方向的污染源導(dǎo)致，農(nóng)村秸稈焚燒或許是主因。

3　結(jié)論

武漢地區(qū)大氣BC濃度存在明顯的時(shí)間變化規(guī)律。2015年武漢市秋冬季BC氣溶膠質(zhì)量濃度日均值為(5.227±3.12)μg·m-3，范圍為1.205～16.534 μg·m-3，其中約 61%分布在 2～6 μg·m-3，并且日均值的月份分布存在明顯差異，9月最低，1月最高。由于冬季相較秋季溫度低，對(duì)流運(yùn)動(dòng)弱，不利于空氣擴(kuò)散，加上受冬季工廠燃煤和取暖等因素影響，故冬季BC濃度整體上高于秋季。

在 10個(gè)城區(qū)公園綠地空間分布中，石門峰公園的BC濃度最高，東湖風(fēng)景區(qū)最低，導(dǎo)致差異的原因與下墊面狀況有關(guān)，與離車流較多的主路的距離等都是原因之一。武漢秋冬季 BC本底濃度為5.45 μg·m-3，與國(guó)內(nèi)外主要城市相比，高于西寧、南極、瓦里關(guān)等偏遠(yuǎn)地區(qū)和天津、芝加哥等經(jīng)濟(jì)較發(fā)達(dá)城市，但低于烏魯木齊、西安等城市。

圖8　武漢市72 h后向軌跡Fig. 8　The 72 h back-trajectories for Wuhan City

BC與O3呈負(fù)相關(guān)，與PM2.5、PM10、CO呈正相關(guān)，相關(guān)系數(shù)分別為-0.436、0.885、0.658、0.8。將BC與O3、PM2.5、PM10、CO進(jìn)行多元線性回歸分析，得到回歸方程：Y=-1136.564+4.216X1+66.316X2-7X3+40.324X4。

氣塊后向軌跡與天氣形勢(shì)和季節(jié)有關(guān)，秋季武漢市大氣BC氣溶膠主要受來(lái)自西北方向空氣氣團(tuán)影響，而冬季主要受來(lái)自北方的空氣氣團(tuán)影響，并且經(jīng)過(guò)山西、河南等傳統(tǒng)的北方污染通道，加重了武漢地區(qū)污染物積聚。兩次BC污染事件可能都受到周邊源影響，農(nóng)民秸稈焚燒可能是主因。

目前，關(guān)于中國(guó)城市黑碳?xì)馊苣z的觀測(cè)研究不少，大多集中于BC的時(shí)間變化特征及其與氣象條件、主要污染物之間的關(guān)系研究，對(duì)BC污染源的定性研究較多，而定量研究較少，因此可以考慮運(yùn)用軌跡模型、CMAQ等全球大尺度氣候模型解釋區(qū)域BC運(yùn)動(dòng)規(guī)律，對(duì)BC來(lái)源進(jìn)行定量解析。

參考文獻(xiàn)：

BABICH PETER, DAVEY MARK, ALLEN GEORGE, et al. 2000. Method Comparisons for Particulate Nitrate, Elemental Carbon, and PM2.5Mass in Seven U.S Cities [J]. Air & Waste Management Association,50(7): 1095-1105.

KIM H J, LIU X D, KOBAYASHI T, et al. 2003. Ultrafine carbon black particles inhibit human lung fibroblast mediated collagen gel contraction [J]. Cellular and molecular biology, 28(1): 111-121.

MURPHY D M, CHOW J C, LEIBENSPERGER E M, et al. 2011.Decreases in elemental carbon and fine particle mass in the United States [J]. Atmospheric Chemistry and Physics, 11(10): 4679-4686.

SHARMA S, LAVOUE D, CACHIER H, et al. 2004. Long-term trends of the black carbon concentrations in the Canadian Arctic [J]. Journal of Geophysical Research, doi: 10.1029/2003JD004331.

VERMA R L, SAHU L K, KONDO Y, et al. 2010. Temporal variations of black carbon in Guangzhou, China, in summer 2006 [J]. Atmospheric Chemistry and Physics, 10: 6471-6485.

COOKE W F, JENNINGS S G, SPAIN T G. 1997. Black carbon measurements at Mace Head, 1989—1996 [J].Journal of Geophysical Research, 102: 25339-25346.

HITZENBERGER R, TOHNO S. 2001. Comparison of Black Carbon (BC)Aerosols in Two Urban Areas-Concentrations and Size Distributions[J]. Atmosphere Environment, 35(12): 2153-2167.

GONG W, ZHANG T H, ZHU Z M, et al. 2015. Characteristics of PM1.0,PM2.5, and PM10, and Their Relation to Black Carbon in Wuhan,Central China [J]. Atmosphere, 6: 1377-1387.

PARK S S, KIM Y J, FUNG K C. 2002. PM2.5Carbon Measurements in Two Urban Areas: Seoul and Kwangju, Korea [J]. Atmospheric Environment, 36(8): 1287-1297.

曹軍驥, 占長(zhǎng)林. 2011. 黑碳在全球氣候和環(huán)境系統(tǒng)中的作用及其在相關(guān)研究中的意義.地球科學(xué)與環(huán)境學(xué)報(bào), 33(2): 177-183.

國(guó)家統(tǒng)計(jì)局. 2016. 2016武漢統(tǒng)計(jì)年鑒[J]. 北京: 中國(guó)統(tǒng)計(jì)出版社.

李楊, 曹軍驥, 張小曳, 等. 2005. 2003年秋季西安大氣中黑碳?xì)馊苣z的演化特征及其來(lái)源解析[J]. 氣候與環(huán)境研究, 10(2): 229-237.

劉新春, 鐘玉婷, 何清, 等. 2011. 2009年冬季烏魯木齊大氣中黑碳?xì)馊苣z觀測(cè)研究[C]//中國(guó)氣象學(xué)會(huì)年會(huì). 2011年第二十八屆中國(guó)氣象學(xué)會(huì)年會(huì)論文集. 廈門.

婁淑娟, 毛節(jié)泰, 王美華. 2005. 北京地區(qū)不同尺度氣溶膠中黑碳含量的觀測(cè)研究[J]. 環(huán)境科學(xué)學(xué)報(bào), 15(1): 17-22.

楊溯, 張武, 韓晶晶, 等. 2008. 上海市浦東新區(qū)秋冬季黑碳?xì)馊苣z特性[J]. 蘭州大學(xué)學(xué)報(bào)(自然科學(xué)版), 44(4): 66-70.

姚青, 蔡子穎, 韓素芹, 等. 2012. 天津城區(qū)秋冬季黑碳?xì)馊苣z觀測(cè)與分析[J]. 環(huán)境化學(xué), 31(3): 324-329.

張華, 王立志. 2009. 黑碳?xì)馊苣z氣候效應(yīng)的研究進(jìn)展[J]. 氣候變化研究進(jìn)展, 5(6): 311-317.

張宇堯. 2017. 基于氣溶膠光學(xué)特性對(duì)武漢黑碳和有機(jī)碳的研究[J]. 測(cè)繪科學(xué)技術(shù), 5(2): 67-75.

趙玉成, 德力格爾, 馬元倉(cāng), 等. 2008. 2005—2006年秋冬西寧大氣中黑碳?xì)馊苣z的濃度變化特征[J]. 城市環(huán)境與城市生態(tài), 21(2): 26-29.

朱厚玲. 2003. 我國(guó)地區(qū)黑碳?xì)馊苣z時(shí)空分布研究[D]. 北京: 中國(guó)氣象科學(xué)研究院.

董俊玲, 張仁健, 符淙斌. 2010. 中國(guó)地區(qū)氣溶膠氣候效應(yīng)研究進(jìn)展[J].中國(guó)粉體技術(shù), 16(1): 1-4.