近20年來京津冀地區(qū)土壤風(fēng)蝕揚塵顆粒物排放變化

趙志平,齊永青,漢瑞英,肖能文,李俊生,*

1 中國環(huán)境科學(xué)研究院,環(huán)境基準(zhǔn)與風(fēng)險評估國家重點實驗室, 北京 100012 2 中國科學(xué)院遺傳與發(fā)育生物學(xué)研究所農(nóng)業(yè)資源研究中心, 石家莊 050022

土壤風(fēng)蝕是風(fēng)力作用導(dǎo)致表土物質(zhì)脫離原空間位置的過程[1—2],其主要危害包括沙漠化[3—4]和各種風(fēng)沙災(zāi)害[5—6]。土壤揚塵源是裸露地面的顆粒物在自然或人力的作用下形成的揚塵[7],其中土壤風(fēng)蝕導(dǎo)致的揚塵源是細顆粒物的主要來源[8]。裸露地面的顆粒物進入大氣中造成大氣污染,導(dǎo)致PM10和PM2.5超標(biāo)[9]。目前灰霾污染已是影響我國城市空氣質(zhì)量的主要因素[10],成為影響人們生活幸福感的重要因素之一[11]。雖然我國空氣質(zhì)量在逐年改善,大氣顆粒物仍是影響我國很多城市空氣質(zhì)量的首要污染物,其中土壤揚塵是造成城市顆粒物污染嚴(yán)重的重要因素之一[12]。

隨著大氣污染問題日益突出,近年來國內(nèi)很多地區(qū)開展了土壤揚塵源顆粒物排放清單編制工作[13—18]。2015年全國風(fēng)蝕揚塵總懸浮顆粒物(TSP)、PM10和PM2.5的年排放量分別為2.27×107t、6.77×106t和1.17×106t[19]。京津冀地區(qū)是土壤風(fēng)蝕揚塵高發(fā)區(qū),其中北京市西北和東南地區(qū)土壤風(fēng)蝕揚塵PM2.5排放因子較高[8],天津市北辰區(qū)PM10和PM2.5的最大排放源是工地[20],河北省土壤揚塵源PM2.5排放量為39699 t,農(nóng)田揚塵排放占85.3%[21],此外2013年石家莊市土壤揚塵TSP、PM10和PM2.5的排放量分別為12569 t、3771 t和947 t[22]。未來氣候變化情景下,河北壩上地區(qū)土壤風(fēng)蝕揚塵年排放速率分別增高25%、54%、35%和54%[23]。綜上所述,目前關(guān)于京津冀地區(qū)土壤揚塵顆粒物排放量研究只有零星報道,且缺乏多年變化狀況研究。整體上近20年來京津冀地區(qū)土壤揚塵顆粒物排放量變化及其主要影響因素仍然未知。

京津冀地區(qū)是我國大氣污染嚴(yán)重區(qū)域,并且容易受到來自內(nèi)蒙的沙塵天氣影響[24—25]。根據(jù)生態(tài)環(huán)境部發(fā)布的空氣質(zhì)量狀況,2019年京津冀及周邊地區(qū)優(yōu)良天數(shù)只有53.1%,低于全國82.0%的平均水平[26]。因此,土壤揚塵顆粒物排放量變化研究對于改善京津冀地區(qū)空氣質(zhì)量具有重要意義。本研究收集2000—2019年京津冀地區(qū)氣候、土壤、植被覆蓋等數(shù)據(jù),估算近20年來京津冀地區(qū)土壤揚塵顆粒物排放量及其變化,揭示其年際變化的影響因素,以期為生態(tài)環(huán)境管理部門提出相應(yīng)措施、降低京津冀地區(qū)土壤揚塵顆粒物排放量、協(xié)同改善空氣質(zhì)量提供科學(xué)支撐。

1數(shù)據(jù)和方法1.1研究區(qū)概況

京津冀位于我國北方地區(qū),東臨大海、西倚太行山,北為燕山,南部為華北平原,地理上界于北緯36°05′—42°40′,東經(jīng)113°27′—119°50′之間,總面積約21.8×104km2,行政區(qū)劃包括北京市、天津市和河北省[27]。氣候上該區(qū)屬于暖溫帶半濕潤大陸性季風(fēng)氣候,四季分明。年平均溫度8—12.5℃,年降水量600—1000 mm。主要土壤質(zhì)地類型有壤土、壤砂土、砂土、輕粘土、砂質(zhì)粘壤土等。該區(qū)北部植被結(jié)構(gòu)復(fù)雜,種類繁多,南部平原區(qū)大多數(shù)以耕地和城市人工植被為主[28]。該區(qū)人類活動強度日益增大[29],使得灰霾天氣頻發(fā)[26],困擾著京津冀地區(qū)環(huán)境保護工作。

1.2 數(shù)據(jù)

本研究用到的氣象數(shù)據(jù)來自于國家氣象科學(xué)數(shù)據(jù)中心(網(wǎng)址http://data.cma.cn/),包括2000—2019年京津冀及周邊地區(qū)氣象站點年平均溫度、年降水量、年平均風(fēng)速,其中年平均溫度和年降水量數(shù)據(jù)利用Anusplin軟件插值成1 km空間柵格數(shù)據(jù),年平均風(fēng)速數(shù)據(jù)利用反距離加權(quán)插值成1 km空間柵格數(shù)據(jù)。土壤質(zhì)地數(shù)據(jù)[30]來源于黑河計劃數(shù)據(jù)管理中心(網(wǎng)址http://westdc.westgis.ac.cn)(圖1)。2000—2019年歸一化植被指數(shù)(NDVI)數(shù)據(jù)來源于MODIS MOD13A2產(chǎn)品,空間分辨率1 km,通過對每16天NDVI求均值得到NDVI年均值。2000、2005、2010和2015年土地利用狀況數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心(網(wǎng)址http://www.resdc.cn/),空間分辨率1 km。

圖1 京津冀地區(qū)土壤質(zhì)地類型Fig.1 Soil texture types in the Beijing-Tianjin-Hebei region

1.3 方法

本文研究對象為土壤揚塵源排放量,以期為生態(tài)環(huán)境管理部門提出相應(yīng)措施、協(xié)同改善空氣質(zhì)量提供科學(xué)參考,因此本研究采用生態(tài)環(huán)境部2015年發(fā)布的《揚塵源顆粒物排放清單編制技術(shù)指南(試行)》[7]中土壤揚塵源排放量計算公式：

WSi=ESi×AS

(1)

式中,WSi為土壤揚塵中PMi(空氣動力學(xué)粒徑在0—iμm間的顆粒物,下同)總排放量,t/a；AS為土壤揚塵源的面積,m2；ESi為土壤揚塵源的PMi排放系數(shù),t m-2a-1,計算公式如下：

ESi=Di×C×(1-η)×10-4

(2)

式中,Di為PMi的起塵因子,t 10-4m-2a-1；C為氣候因子,表征氣象因素對土壤揚塵的影響；η為污染控制技術(shù)對揚塵的去除效率,%,本研究設(shè)定η=0。Di計算公式如下：

Di=ki×Iwe×f×L×V

(3)

式中,ki為PMi在土壤揚塵中的百分含量,推薦值TSP為1、PM10為0.30、PM2.5為0.05；Iwe為土壤風(fēng)蝕指數(shù),推薦值見表1；f為地面粗糙因子,取值為0.5；L為無屏蔽寬度因子,本研究取L=1.0。V為植被覆蓋因子,是指裸露土壤面積占總計算面積的比例,本研究利用NDVI數(shù)據(jù)來代替,該數(shù)據(jù)來源于MODIS MOD13A2產(chǎn)品。

氣候因子C計算公式如下：

C=0.504×u3/PE2

(4)

式中,u為年平均風(fēng)速,m/s；PE為桑氏威特降水—蒸發(fā)指數(shù),計算公式如下：

PE=1.099×P/(0.5949+0.1189×T)

(5)

式中,P為年降水量,mm；T為年平均溫度,℃。

表1 土壤風(fēng)蝕指數(shù)參考值/(t hm-2 a-1)

圖2 2000—2019年京津冀地區(qū)風(fēng)速變化 Fig.2 Changes of wind speed in the Beijing-Tianjin-Hebei region from 2000 to 2019

2 結(jié)果

2.1 氣候因子變化

2.1.1風(fēng)速

2000—2019年京津冀地區(qū)年平均風(fēng)速均值為2.54 m/s(圖2)?？傮w上2000—2019年京津冀地區(qū)年平均風(fēng)速具有略微上升趨勢,但上升速率非常小為0.52×10-2m/s a-1,且達到一定的顯著水平(P<0.1)。

2.1.2溫度和降水

2000—2019年京津冀地區(qū)年平均溫度均值為9.77℃,其中2008年平均溫度最高為12.81℃,2010年平均溫度最低為6.92℃,近20年京津冀地區(qū)年平均溫度呈現(xiàn)升高趨勢(圖3),升高速率約為0.04℃/a,但未達到顯著水平(P>0.05)。年降水量也呈現(xiàn)增加趨勢,增加速率約為4.12 mm/a,且達到一定顯著水平(P<0.1)。2000—2019年京津冀地區(qū)年降水量均值為502.23 mm,其中2016年降水量最多為632.50 mm,2002年最低為396.40 mm。

2.1.3氣候因子和桑氏威特蒸發(fā)指數(shù)

氣候因子是土壤揚塵源顆粒物排放系數(shù)計算的重要因子。2000—2019年京津冀地區(qū)氣候因子呈現(xiàn)下降趨勢(圖4),但未達到顯著水平(P>0.05)。其中2002年最大為1.51×10-4,2012最小為0.46×10-4。桑氏威特蒸發(fā)指數(shù)由年平均溫度和年降水量計算得到,是構(gòu)成氣候因子的重要參數(shù)。2000—2019年京津冀地區(qū)桑氏威特蒸發(fā)指數(shù)呈現(xiàn)上升趨勢,但未達到顯著水平(P>0.05)。

2.2 植被覆蓋變化

2000—2019年京津冀地區(qū)植被NDVI均值為0.38,2000年最低為0.34,2017年最高為0.40。2000—2019年京津冀地區(qū)植被NDVI年均值呈現(xiàn)極顯著增加趨勢(P<0.01),增加速率為0.21×10-2/a(圖5),與NDVI均值相比約為5.53%/10a,雖然反映了近20年來該地區(qū)植被生長狀況整體呈現(xiàn)變好趨勢,但變化幅度較小。

圖3 2000—2019年京津冀地區(qū)年平均溫度和年降水量變化Fig.3 Annual temperature and precipitation in Beijing-Tianjin-Hebei region during 2000—2019

圖4 2000—2019年京津冀地區(qū)氣候因子和桑氏威特指數(shù)變化Fig.4 Changes in climate factor and Thornthwaite index in the Beijing-Tianjin-Hebei region from 2000 to 2019

圖5 2000—2019年京津冀地區(qū)植被NDVI變化 Fig.5 Changes of vegetation NDVI in the Beijing-Tianjin-Hebei region from 2000 to 2019

2.3 土壤揚塵源排放系數(shù)變化

空間上2019年京津冀地區(qū)土壤揚塵源TSP排放系數(shù)均值約為2.30 t km-2a-1,其中最大為23.95 t km-2a-1,最小為0.03 t km-2a-1。TSP排放系數(shù)最高地區(qū)主要分布在河北省西南部地區(qū),包括保定市西南部、石家莊市西部、邢臺市中西部和邯鄲市西北部地區(qū)(圖6)。TSP排放系數(shù)較高地區(qū)主要分布在天津市南部、滄州市東北部以及張家口市局部地區(qū)。TSP排放系數(shù)較低地區(qū)主要分布在京津冀北部地區(qū)。2019年京津冀地區(qū)土壤揚塵源PM10和PM2.5排放系數(shù)分布格局與TSP一致。

圖6 2019年京津冀地區(qū)土壤揚塵源總懸浮顆粒物(TSP)排放系數(shù)Fig.6 TSP emission coefficient of soil dust source in the Beijing-Tianjin-Hebei region in 2019

2000—2019年京津冀地區(qū)土壤揚塵源TSP排放系數(shù)均值為1.79 t km-2a-1,其中PM10占8.99%,PM2.5占0.25%。近20年TSP排放系數(shù)具有下降趨勢(圖7),但未達到顯著水平(P>0.05)。其中2002年最高為3.27 t km-2a-1,2003年最低為1.01 t km-2a-1。相關(guān)分析顯示,近20年土壤揚塵源TSP排放系數(shù)變化與氣候因子呈現(xiàn)極顯著正相關(guān)(P<0.01),與植被NDVI指數(shù)呈現(xiàn)負相關(guān),但未達到顯著水平(P>0.05),表明京津冀地區(qū)土壤揚塵源TSP排放系數(shù)變化與氣候因子有關(guān),而氣候因子與溫度和降水具有極顯著相關(guān)關(guān)系(P<0.01)。2000—2019年京津冀地區(qū)土壤揚塵源PM10和PM2.5排放系數(shù)變化過程與TSP一致。

2.4 土壤揚塵源排放量空間格局

2000—2019年京津冀地區(qū)土壤揚塵源TSP、PM10和PM2.5排放系數(shù)均值分別為1.79、0.16、0.45×10-2t km-2a-1(表2)。滄州市土壤揚塵源TSP、PM10和PM2.5排放系數(shù)均值最高為3.17、0.29、0.80×10-2t km-2a-1,其次是天津市和石家莊市,分別為2.96、0.27、0.75×10-2、2.62、0.24、0.66×10-2t km-2a-1。承德市土壤揚塵源TSP、PM10和PM2.5排放系數(shù)均值最低為0.60、0.05、0.15×10-2t km-2a-1,其次是秦皇島市0.71、0.06、0.18×10-2t km-2a-1。

2000—2019年京津冀地區(qū)土壤揚塵源TSP、PM10和PM2.5排放總量均值分別為3.81×105t/a、0.34×105t/a、0.96×103t/a。其中張家口市土壤揚塵源TSP排放量均值最大為0.73×105t/a,占京津冀地區(qū)總量的19.18%,其次是保定市和滄州市分別為0.49×105t/a和0.44×105t/a,分別占12.98%和11.63%。秦皇島市土壤揚塵源TSP排放量均值最小為0.06×105t/a,占京津冀地區(qū)總量的1.44%,其次是廊坊市和唐山市分別為0.12×105t/a和0.15×105t/a,分別占京津冀地區(qū)總量的3.27%和4.00%。

2000—2019年承德市、秦皇島市、邢臺市和邯鄲市土壤揚塵源TSP、PM10和PM2.5排放系數(shù)具有上升趨勢,但未達到顯著水平(P>0.05)。其余地區(qū)均具有下降趨勢,其中北京市、衡水市土壤揚塵源TSP、PM10和PM2.5排放系數(shù)下降趨勢達到極顯著水平(P<0.01),天津市、廊坊市和保定市下降趨勢達到顯著水平(P<0.05)。上述土壤揚塵源TSP、PM10和PM2.5排放系數(shù)顯著下降區(qū)域的排放量僅占京津冀總量的35.63%,京津冀大部分區(qū)域土壤揚塵源顆粒物排放狀況未得到有效遏制。

本研究利用2000、2005、2010和2015年京津冀地區(qū)4期土地利用空間狀況數(shù)據(jù)來統(tǒng)計2000—2019年京津冀地區(qū)土壤揚塵源TSP排放量,結(jié)果顯示耕地排放量最大為2.29×105t/a(圖8),占京津冀地區(qū)總量的59.83%,是抑制土壤揚塵源顆粒物排放的重點關(guān)注對象,其次為草地為0.60×105t/a,占15.66%。林地和居民地土壤揚塵源TSP排放量分別為0.42×105t/a和0.36×105t/a,分別占11.04%和9.50%。水體排放量為0.14×105t/a,占3.64%。裸地排放量最低為0.13×104t/a,占京津冀地區(qū)總量的0.33%。

表2 2000—2019年京津冀各市土壤揚塵源TSP、PM10和PM2.5排放量

圖7 2000—2019年京津冀地區(qū)土壤揚塵源TSP、PM10和PM2.5排放系數(shù)變化Fig.7 Changes in the emission coefficient of TSP,PM10 and PM2.5 of soil dust sources in the Beijing-Tianjin-Hebei region from 2000 to 2019

圖8 2000—2019年京津冀地區(qū)各土地利用類型土壤揚塵源TSP、PM10和PM2.5排放量Fig.8 Emissions of TSP, PM10 and PM2.5 from soil dust sources of various land use types in the Beijing-Tianjin-Hebei region from 2000 to 2019

3 討論

2021年3月初以來,京津冀地區(qū)遭遇多場沙塵暴襲擊,是揚塵污染嚴(yán)重區(qū)。本研究參考生態(tài)環(huán)境部2015年發(fā)布的《揚塵源顆粒物排放清單編制技術(shù)指南(試行)》[7],估算京津冀地區(qū)年度土壤揚塵顆粒物排放。該方法中用到的變量均為年值,因此不能用于季節(jié)土壤揚塵顆粒物排放計算,否則會出現(xiàn)較大誤差。與以往京津冀地區(qū)研究結(jié)果比較發(fā)現(xiàn),本研究2000—2019年北京市土壤揚塵PM2.5排放系數(shù)均值為0.45×10-2t km-2a-1,其中2017年為0.44×10-2t km-2a-1,以往學(xué)者認為2017年北京市土壤揚塵PM2.5排放系數(shù)為0.18 t km-2a-1,與本研究相比存在一定低估[8]。此外,還有研究認為2015年河北省各地級市中承德市、廊坊市、唐山市、秦皇島市土壤揚塵顆粒物排放系數(shù)較低,滄州市、邢臺市和邯鄲市土壤揚塵顆粒物排放系數(shù)較高[21]。本研究結(jié)果與之相比土壤揚塵顆粒物排放系數(shù)空間格局大體一致,即滄州市、石家莊市土壤揚塵顆粒物排放系數(shù)較高,承德市、秦皇島市、唐山市土壤揚塵顆粒物排放系數(shù)較低。2013年石家莊市土壤揚塵TSP排放為12569 t、PM10為3771 t、PM2.5為947 t[22],轉(zhuǎn)換后平均排放系數(shù)分別為0.87、0.26、0.07 t km-2a-1,分別為本研究(2.62、0.24、0.07×10-1t km-2a-1)的0.33倍、1.08倍和10倍,相比本研究TSP存在一定低估、PM10接近、PM2.5存在一定高估。總體上本研究中京津冀地區(qū)土壤揚塵顆粒物排放量估算相對合理。

本研究發(fā)現(xiàn)2000—2019年土壤揚塵源TSP排放系數(shù)變化與氣候因子呈現(xiàn)極顯著正相關(guān)(圖9,P<0.01),因此京津冀地區(qū)土壤揚塵源TSP排放系數(shù)變化主要與氣候因子有關(guān)。2000—2019年京津冀地區(qū)植被NDVI呈現(xiàn)極顯著上升趨勢,反映了近20年來該地區(qū)植被生長狀況整體呈現(xiàn)變好,但由于上升幅度較小,僅為5.53%/10a,與土壤揚塵源TSP排放系數(shù)變化呈現(xiàn)一定顯著水平負相關(guān)(圖9,P<0.1),因此植被覆蓋是土壤揚塵源TSP排放系數(shù)變化的次要因素。進一步分析發(fā)現(xiàn)氣候因子中年降水量與土壤揚塵源TSP排放系數(shù)變化具有極顯著負相關(guān)關(guān)系(圖9,P<0.01),而年平均風(fēng)速和年平均溫度與土壤揚塵源TSP排放系數(shù)變化相關(guān)性不顯著(P>0.05),表明京津冀地區(qū)土壤揚塵源TSP排放系數(shù)變化主要受年降水量影響。

圖9 2000—2019年京津冀地區(qū)土壤揚塵源總懸浮顆粒物(TSP)排放系數(shù)與氣候因子、覆蓋因子、年降水量相關(guān)性Fig.9 Correlations of Total Suspended Particulate (TSP) emission coefficient and climate factor, coverage factor, and annual precipitation in the Beijing-Tianjin-Hebei region from 2000 to 2019

本研究假設(shè)京津冀地區(qū)土壤揚塵源可能影響的高度為100 m,同時利用2019年各市土壤揚塵源排放系數(shù)計算得到土壤揚塵源PM10和PM2.5在空氣中含量,并與2019年京津冀各地區(qū)發(fā)布的環(huán)境質(zhì)量公報結(jié)果進行對比。結(jié)果顯示土壤揚塵源排放顆粒對PM10影響較大,對PM2.5影響較小(表3)。2019年邢臺市PM10觀測值115 μg/m3,其中土壤揚塵源排放顆粒物含量占比最高為12.66%。石家莊市和天津市PM10觀測值分別為118 μg/m3和76 μg/m3,其中土壤揚塵源排放顆粒物含量占比也較高分別為11.09%和10.30%。滄州市和邯鄲市PM10觀測值分別為89 μg/m3和124 μg/m3,其中土壤揚塵源排放顆粒物含量占比分別為8.63%和8.02%。上述地區(qū)應(yīng)關(guān)注土壤揚塵源顆粒物排放對空氣質(zhì)量影響。秦皇島市和承德市土壤揚塵源排放顆粒物含量占PM10觀測值比例分別為2.55%和2.99%,可能與當(dāng)?shù)刂脖桓采w狀況較好有關(guān)。

表3 2019年京津冀各市土壤揚塵源排放顆粒物含量占比

本研究估算了2000—2019年京津冀地區(qū)土壤揚塵源排放,從長時間尺度上闡明該地區(qū)土壤揚塵顆粒物排放空間格局和變化過程,發(fā)現(xiàn)京津冀大部分區(qū)域土壤揚塵源顆粒物排放狀況未得到有效遏制。人為翻動表土?xí)?dǎo)致嚴(yán)重的土壤風(fēng)蝕發(fā)生,人類擾動會導(dǎo)致土壤風(fēng)蝕加劇。本研究中耕地揚塵源TSP排放量最大占京津冀地區(qū)總量的59.83%,是抑制土壤揚塵源顆粒物排放的重點關(guān)注對象,其次為草地占15.66%。該結(jié)果可為京津冀地區(qū)各市相關(guān)管理部門制定降低大氣顆粒物含量、減少土壤揚塵顆粒物排放措施的依據(jù)。同時相關(guān)部門可以引入空氣質(zhì)量模型分析土壤揚塵源對大氣環(huán)境空氣質(zhì)量的影響,以探索適合京津冀地區(qū)實際情況的土壤揚塵防治策略。主要包括京津冀地區(qū)沙塵抑制生態(tài)屏障協(xié)同建設(shè),統(tǒng)籌協(xié)調(diào)退耕還林還草和自然植被恢復(fù)工作,以及沙塵天氣動態(tài)監(jiān)測和預(yù)警系統(tǒng)建設(shè)[31]。由于氣候因子和土壤質(zhì)地數(shù)據(jù)來源和分辨率的差異,本研究結(jié)果與前人的結(jié)果存在一定差異,下一步工作應(yīng)該加強土壤揚塵顆粒物排放研究的細致對比分析工作,同時搜集實測數(shù)據(jù),結(jié)合其他源的排放清單,同顆粒物源解析結(jié)果進行比較驗證。

4 結(jié)論

(1)2000—2019年京津冀地區(qū)土壤揚塵源TSP排放系數(shù)均值為1.79 t km-2a-1,其中PM10占8.99%,PM2.5占0.25%。近20年土壤揚塵源TSP排放系數(shù)具有下降趨勢,PM10和PM2.5排放系數(shù)變化過程與TSP一致。上述變化主要受氣候因子影響,其次受植被覆蓋度影響。進一步分析發(fā)現(xiàn)近20年來京津冀地區(qū)土壤揚塵源TSP排放系數(shù)變化主要受年降水量影響。

(2)2000—2019年滄州市土壤揚塵源TSP、PM10和PM2.5排放系數(shù)均值最高,其次是天津市和石家莊市。張家口市土壤揚塵源TSP排放量均值最大占京津冀地區(qū)總量的19.18%,其次是保定市和滄州市分別占12.98%和11.63%。北京市、衡水市、天津市、廊坊市和保定市土壤揚塵源TSP、PM10和PM2.5排放系數(shù)顯著下降,然而上述區(qū)域土壤揚塵排放量僅占京津冀總量的35.63%,京津冀大部分區(qū)域土壤揚塵源顆粒物排放狀況未得到有效遏制。耕地土壤揚塵排放量最大占京津冀地區(qū)總量的59.83%,是抑制土壤揚塵源顆粒物排放的重點關(guān)注對象,其次為草地占15.66%。林地和居民地分別占11.04%和9.50%。

(3)2019年邢臺市土壤揚塵源PM10排放占觀測值比例最高為12.66%,石家莊市和天津市占比也較高分別為11.09%和10.30%,滄州市和邯鄲市占比分別為8.63%和8.02%。上述地區(qū)環(huán)境管理部門均應(yīng)關(guān)注土壤揚塵源顆粒物排放對空氣質(zhì)量的影響。秦皇島市和承德市土壤揚塵源排放顆粒物含量占PM10觀測值比例分別為2.55%和2.99%,可能與當(dāng)?shù)刂脖桓采w狀況較好有關(guān)。