北京市平原區(qū)裸露地風(fēng)蝕揚(yáng)塵排放量

劉奧博,吳其重,陳雅婷,趙天成,程 曉 (北京師范大學(xué)全球變化與地球系統(tǒng)科學(xué)研究院,北京 100048)

目前,北京市以可吸入顆粒物(PM10)、細(xì)顆粒物(PM2.5)為代表的區(qū)域性大氣環(huán)境問題依然嚴(yán)峻,危及人民群眾的身體健康[1-2].準(zhǔn)確完整的大氣污染源排放清單為區(qū)域大氣污染治理工作提供了數(shù)據(jù)支撐和理論依據(jù),是控制和改善區(qū)域環(huán)境空氣質(zhì)量的重要基礎(chǔ)[3-5].研究表明,揚(yáng)塵排放是我國(guó)北方城市大氣顆粒物中最主要的貢獻(xiàn)源之一[6],郊區(qū)地表風(fēng)蝕起塵對(duì)城區(qū)空氣質(zhì)量影響巨大[7].據(jù)環(huán)保部門發(fā)布的大氣顆粒物源解析結(jié)果,北京市、天津市、石家莊市的PM2.5來(lái)源中,揚(yáng)塵源占本地排放的比例分別為14.3%、30%、22.5%,對(duì)PM10的貢獻(xiàn)可高達(dá)40%[8-10].

典型的揚(yáng)塵污染源包括交通揚(yáng)塵、料堆揚(yáng)塵、施工揚(yáng)塵和裸露面揚(yáng)塵等[11].裸露面揚(yáng)塵是揚(yáng)塵排放的重要一類,是自然裸露地面受風(fēng)侵蝕排放的顆粒物,極端過程如沙塵暴,可以顯著影響大氣環(huán)境質(zhì)量.土壤風(fēng)蝕是中國(guó)北方干旱與半干旱區(qū)土地退化或荒漠化的主要過程之一,也是空氣中懸浮顆粒物的主要貢獻(xiàn)源[12].顆粒物風(fēng)蝕揚(yáng)塵過程的主要影響因素包括氣候條件、土壤性質(zhì)、地形和植被覆蓋等[13].北京雖屬半濕潤(rùn)區(qū),但在冬春季節(jié)植物休眠地表裸露,與干旱大風(fēng)天氣同期,風(fēng)蝕現(xiàn)象較為嚴(yán)重.

裸露面風(fēng)蝕揚(yáng)塵排放的常規(guī)測(cè)算以定點(diǎn)儀器監(jiān)測(cè)和野外取樣分析為主,雖測(cè)量結(jié)果較為準(zhǔn)確,但在空間上的可擴(kuò)展性較差,難以全面展現(xiàn)北京市平原區(qū)裸露地的揚(yáng)塵污染源狀況和分布特征,不利于大氣環(huán)境質(zhì)量及其污染來(lái)源的宏觀分析.遙感技術(shù)能夠動(dòng)態(tài)地反映研究區(qū)域內(nèi)土地覆蓋的變化過程,因而被廣泛用于揚(yáng)塵源的識(shí)別[14].故本研究以NASA的Landsat系列遙感衛(wèi)星影像為數(shù)據(jù)源,采用全自動(dòng)裸土提取算法,高效準(zhǔn)確地提取出北京平原區(qū)的裸土圖斑,以此結(jié)合氣候觀測(cè)資料和揚(yáng)塵排放模型,開展了北京市平原區(qū)裸露地風(fēng)蝕揚(yáng)塵排放清單的編制工作,對(duì)裸露地風(fēng)蝕揚(yáng)塵排放量進(jìn)行了估算,為北京市環(huán)境統(tǒng)計(jì)和污染檢測(cè)提供了一種新的技術(shù)手段,并為揚(yáng)塵污染控制的政策制定提供了輔助信息和理論依據(jù).

1 裸露地提取

1.1 研究區(qū)域和數(shù)據(jù)源

本文的研究區(qū)域?yàn)楸本┦衅皆瓍^(qū)(圖1),選用的數(shù)據(jù)源為美國(guó)航空航天局(NASA)發(fā)布的Landsat衛(wèi)星遙感影像.為保證裸露地提取精度,研究區(qū)范圍內(nèi)積雪和云霧覆蓋需少于15%.經(jīng)篩選,合格的影像數(shù)據(jù)分布如表1所示.

1.2 構(gòu)建增強(qiáng)型裸土指數(shù)

根據(jù)《揚(yáng)塵源顆粒物排放清單編制技術(shù)指南(試行)》[15](簡(jiǎn)稱《指南》),裸露地包括無(wú)植被覆蓋的農(nóng)田、未硬化或未綠化的空地、干涸的河谷、裸露山體、灘涂等.本文研究區(qū)內(nèi)主要的地表覆蓋物類型為水體、植被、建筑物和裸露地.其中,水體和植被光譜特征明顯,較容易區(qū)分,而裸露地和建筑物存在“異物同譜”的現(xiàn)象,難以準(zhǔn)確區(qū)分.

圖1 研究區(qū)域北京市平原Fig.1 The study domain Beijing plain area

表1 Landsat系列遙感數(shù)據(jù)源Table 1 Landsat remote sensing data series

本文借用吳志杰等[16]提出的增強(qiáng)型裸土指數(shù)(EBSI)進(jìn)行裸露地信息提取,其理論依據(jù)是裸地在裸土指數(shù)(BSI)[17]和修正歸一化水體指數(shù)(MNDWI)[18]圖像上的亮度反差最大,通過差值可以最大程度地增強(qiáng)裸露地信息.其模型表達(dá)式為:

式中: BSI和MNDWI分別為裸土指數(shù)和修正歸一化水體指數(shù)圖像經(jīng)拉伸的灰度值.BSI和MNDWI的計(jì)算公式為:

式中: TM1～5分別代表TM/ETM+影像的1～5波段,對(duì)應(yīng)Landsat8OLI陸地成像儀的2～6波段.在EBSI影像上,水體和植被的亮度值小于0,建筑用地的亮度值小于裸露地且接近于0.1.

1.3 實(shí)現(xiàn)動(dòng)態(tài)閾值分割

對(duì)遙感影像進(jìn)行增強(qiáng)型裸土指數(shù)計(jì)算后,為實(shí)現(xiàn)裸露地信息的全自動(dòng)提取,還需實(shí)現(xiàn)動(dòng)態(tài)閾值分割方法.在EBSI特征影像上,裸露地和建筑物、植被和水體的亮度值相近,在直方圖上表征為典型的雙峰分布,而常用的動(dòng)態(tài)閾值分割算法如迭代法[19]、最大類間方差法[20]和直方圖凹面分析法[21]確定的分割閾值都位于雙峰間的谷底,不適用于裸露地提取.受植被覆蓋影響,冬春季EBSI影像均值小于夏季,需要進(jìn)行一定的修正,經(jīng)試驗(yàn)發(fā)現(xiàn)全局最佳閾值與EBSI和土壤調(diào)節(jié)植被指數(shù)SAVI)[22]近似滿足如下關(guān)系:

式中: MeanEBSI和MeanSAVI分別為EBSI和SAVI影像的均值,β為比例因子,取0.2.

式中: TM3、TM4分別代表TM/ETM+影像的3、4波段,對(duì)應(yīng)OLI的4、5波段;L為土壤調(diào)節(jié)因子,取0.5時(shí)可將土壤亮度差異減到最小.

1.4 分類精度評(píng)定

表2 遙感方法裸露地提取精度驗(yàn)證Table 2 Validation of remote sensing method of bare soil extraction

對(duì)EBSI影像進(jìn)行二值分割(裸露地和非裸露地)后,采用隨機(jī)抽樣方式選擇5期影像共1000個(gè)驗(yàn)證像元,目視解譯的參照影像為同時(shí)期的TM4(R)、3(G)、2(B)假彩色合成圖像.總分類精度可達(dá)88.9%,Kappa系數(shù)為0.77(表2).

2 揚(yáng)塵排放計(jì)算

2.1 揚(yáng)塵排放模型

本研究針對(duì)風(fēng)蝕揚(yáng)塵產(chǎn)生的可吸入顆粒物(PM10)和細(xì)顆粒物(PM2.5)在不同區(qū)域或不同時(shí)期內(nèi)的排放系數(shù)進(jìn)行計(jì)算.本文采用《指南》推薦的風(fēng)蝕方程[15,23],其計(jì)算公式如下:

式中:裸露地?fù)P塵中PMi(空氣動(dòng)力學(xué)粒徑在0～iμm間的顆粒物)的排放系數(shù),t/(km2·a); Iwe為土壤風(fēng)蝕指數(shù),t/(hm2·a); ki為PMi在土壤揚(yáng)塵中的百分含量; f為地面粗糙因子,對(duì)于光滑地面取值為1,粗糙地面為0.5; L為無(wú)屏障寬度因子,指沒有明顯阻擋物的最大范圍,當(dāng)無(wú)屏蔽寬度＜300m時(shí),L＝0.7,當(dāng)無(wú)屏蔽寬度為300～600m時(shí),L＝0.85,當(dāng)無(wú)屏蔽寬度≥600m時(shí),L=1.0; V為植被覆蓋因子,是指裸露土壤面積占總計(jì)算面積的比例,地面完全裸露時(shí)為1; C為氣候因子,表征氣候因素對(duì)土壤揚(yáng)塵的影響; u為年平均風(fēng)速,m/s; Pe為桑氏威特( Thornthwaite )降水-蒸發(fā)指數(shù); P為年降水量,mm; E*為年潛在蒸發(fā)量,mm;Ta為年平均溫度,℃.

2.2 確定模型參數(shù)

由式(6)～式(9)可知,計(jì)算揚(yáng)塵排放系數(shù)的基礎(chǔ)參數(shù)為土壤風(fēng)蝕指數(shù)Iwe,其定義為每公頃裸露土地每年潛在的風(fēng)蝕揚(yáng)塵的噸數(shù),反映了土壤的可風(fēng)蝕性.本文參考彭應(yīng)登等[24]測(cè)得的北京市各區(qū)裸露農(nóng)田的風(fēng)蝕指數(shù)作為區(qū)域內(nèi)的土壤風(fēng)蝕指數(shù)參照值.百分含量ki采用推薦值,PM10為0.3,PM2.5為0.05;地面粗糙因子f取0.5,無(wú)屏障寬度因子L取1,植被覆蓋因子V取0.9;模型內(nèi)的氣候參數(shù)包括溫度、降水和風(fēng)速,均采用中國(guó)氣象局公開的北京市15個(gè)氣象站點(diǎn)的氣候統(tǒng)計(jì)資料,數(shù)據(jù)獲取時(shí)間為1980～2010年,站點(diǎn)分布如表3所示.

表3 北京市氣象站點(diǎn)分布Table 3 Distribution of meteorological stations in Beijing

2.3 年排放系數(shù)計(jì)算

由于各區(qū)域土質(zhì)類型和氣候因子的不同,自然裸露地面的起塵能力有所差異,反映為揚(yáng)塵排放潛力的強(qiáng)弱.故將各氣象站點(diǎn)觀測(cè)的年均風(fēng)速、降水和溫度代入式(7)～式(9),計(jì)算得到風(fēng)蝕揚(yáng)塵的氣候因子(C).然后將2.2節(jié)確定的模型參數(shù)代入式(6), 最終計(jì)算得到各區(qū)的裸露地風(fēng)蝕揚(yáng)塵排放系數(shù)(Qei).

3 結(jié)果與討論

3.1 裸露地面積的時(shí)間變化

北京市平原區(qū)主要的土地利用類型為城鎮(zhèn)建筑用地、水體、農(nóng)用地和自然綠地等.在城區(qū),裸露地主要為建筑施工裸地,而在郊區(qū)和廣大農(nóng)村地區(qū),裸露地主要為未利用的空地和作物收獲后的裸露農(nóng)田.據(jù)北京市統(tǒng)計(jì)年鑒,北京市農(nóng)田面積在3000km2以上,故裸露地面積存在顯著的季節(jié)變化.1～3月裸露地面積最大,可高達(dá)4500km2;夏季草地和農(nóng)田被植被覆蓋時(shí),裸地面積最小,8月平均裸地面積僅為500km2.

因遙感數(shù)據(jù)源無(wú)法保證影像數(shù)據(jù)的均勻覆蓋,簡(jiǎn)單年平均計(jì)算無(wú)法準(zhǔn)確表示裸露地面積的年際變化.故本文采用滑動(dòng)平均的方法,計(jì)算裸露地面積的5a滑動(dòng)平均值,并同樣處理氣候參數(shù)計(jì)算揚(yáng)塵的PM10排放系數(shù),結(jié)果如圖2.

據(jù)唐新明等[25]的研究,1989～2012年間,北京市生態(tài)用地的面積減少了443km2,與之相應(yīng)的是建設(shè)用地面積增加了437km2.而本研究中1987～2016年間裸露地面積共減少了約600km2(圖2),故推測(cè)裸露地面積減小的主因是城市發(fā)展、城區(qū)擴(kuò)張,自然裸露地如農(nóng)田、草地被開發(fā)為建城區(qū).而揚(yáng)塵速率的年際變化波動(dòng)較為明顯,裸露地?fù)P塵源PM10排放系數(shù)的高值可達(dá)7～9t/(km2·a),最高值出現(xiàn)在1999～2003年間,恰與北京沙塵天氣頻發(fā)的時(shí)期相吻合.這是由于沙塵暴和揚(yáng)塵在起沙機(jī)制方面具有相似之處,且沙塵暴之后常伴有沙塵附著在裸露地表面,受擾動(dòng)后極易形成二次揚(yáng)塵,推測(cè)沙塵天氣對(duì)裸露地?fù)P塵源的揚(yáng)塵排放系數(shù)有著較為顯著的影響.同時(shí),我們也注意到,從2011年開始,PM10揚(yáng)塵排放系數(shù)存在上升趨勢(shì).

圖2 裸露地面積和PM10排放系數(shù)的5a滑動(dòng)平均變化Fig.2 Annual variations of bare soil area and PM10 emission factor by 5-years Moving Average

3.2 裸露地的空間分布

以2016年12月的分類結(jié)果影像為例(圖3).北京市平原區(qū)裸露地分布廣泛,且裸露地主要集中在郊區(qū)和農(nóng)村地區(qū),首都核心區(qū)和功能拓展區(qū)裸露地分布較少.圖中分類結(jié)果未經(jīng)濾波或掩膜處理.首都各功能區(qū)的裸露地面積分布如表4所示,2016年冬季北京市平原區(qū)裸露地面積近3600km2,約占全市總面積的21.86%、平原區(qū)面積的56.59%.裸露地主要分布在城市發(fā)展新區(qū)和生態(tài)涵養(yǎng)發(fā)展區(qū),分別占60%和34%,其中,大興區(qū)和順義區(qū)的裸露地面積均超過550km2,通州區(qū)和延慶縣次之,門頭溝區(qū)、懷柔區(qū)和房山區(qū)分布著大面積的林地,所以裸露地面積較小.首都核心區(qū)僅存在零星的裸露地分布,且多為易受侵蝕的屋頂,城市功能拓展區(qū)的裸露地主要分布在五環(huán)路和六環(huán)路之間,包括建筑施工裸地、農(nóng)用地和未利用的空地.

圖3 2016年12月裸露地分布Fig.3 Spatial distribution of bare soil in December 2016

3.3 顆粒物排放特征

北京市四大功能區(qū)16個(gè)區(qū)的裸露地風(fēng)蝕揚(yáng)塵顆粒物排放量計(jì)算結(jié)果如表4所示.因?yàn)椴煌筋w粒物排放量在模型中的差異僅反映為百分比含量值,故PM10與PM2.5的排放量存在著對(duì)應(yīng)關(guān)系.北京市裸露地風(fēng)蝕揚(yáng)塵在空間上分布不均勻,各區(qū)縣的風(fēng)蝕程度和風(fēng)蝕強(qiáng)度存在很大差異,其中大興區(qū)、通州區(qū)和延慶縣的土壤風(fēng)蝕現(xiàn)象相對(duì)較為嚴(yán)重,大興區(qū)的裸露地風(fēng)蝕揚(yáng)塵排放量可達(dá)1913t/a,占全市總排放量的25%.分析認(rèn)為大興區(qū)的土壤類型主要為砂土,土壤風(fēng)蝕指數(shù)大,易發(fā)生風(fēng)蝕,且大興區(qū)恰處于風(fēng)廊帶,大風(fēng)天氣頻發(fā),故土壤的風(fēng)蝕作用較強(qiáng).

表4 北京市平原區(qū)各區(qū)裸露地?fù)P塵排放量Table 4 Bare Soil Dust emissions over different districts in Beijing plain area

北京市各區(qū)縣的PM10排放系數(shù)范圍為(2.1±1.3) t/(km2·a),而聶磊等[26]研究獲得北京市郊區(qū)農(nóng)田風(fēng)蝕揚(yáng)塵產(chǎn)生的PM10的平均排放因子為0.82t/(km2·a),與本文的研究結(jié)果相近.排放總量上,以PM10排放量為例,全市總排放量約為7600t/a,這一數(shù)字大于北京市平原區(qū)民用散煤燃燒[27]顆粒物的貢獻(xiàn),但小于北京市建筑施工裸地的揚(yáng)塵排放量[28]和交通揚(yáng)塵量[29].

徐媛倩等[5]同樣采用土壤風(fēng)蝕模型對(duì)鄭州市的裸露面風(fēng)蝕揚(yáng)塵進(jìn)行研究,但其選擇了2013年4月和8月兩景衛(wèi)星影像的裸露區(qū)域作為永久性裸露地面(未考慮季節(jié)性裸露的農(nóng)田),故提取的裸露地面積僅為208km2.但因?yàn)猷嵵菔械臍夂蛞蛩叵噍^于北京更易于土壤風(fēng)蝕作用的形成(溫度更高,降水量更少,風(fēng)速更大),所以其估算的PM10年排放量可達(dá)3581t/a.王瑋等[30]對(duì)北京市中心城區(qū)(即核心區(qū)和功能拓展區(qū))的裸露面揚(yáng)塵排放進(jìn)行研究,估算得到的PM10排放量為420.7t/a,而本文相應(yīng)的PM10排放量為506t/a,結(jié)果相近.葉芝祥等[31]采用監(jiān)測(cè)儀器和FDM模型對(duì)成都市某處未施工裸地在2014年4月～5月的PM10排放因子進(jìn)行測(cè)算,結(jié)果為2.58×10-4t/(km2·a),比本文估算的排放因子大2個(gè)數(shù)量級(jí).

據(jù)分析,受氣候因素影響,土壤風(fēng)蝕作用的季節(jié)性差異顯著,所以采用年均氣候參數(shù)估算的揚(yáng)塵排放量比起實(shí)驗(yàn)觀測(cè)可能存在系統(tǒng)性的低估.

3.4 不確定性分析

在裸露地?fù)P塵排放估算過程中,計(jì)算結(jié)果存在一定的不確定性,主要因?yàn)?

(1)受算法和衛(wèi)星影像分辨率限制,裸露地提取結(jié)果存在一定的誤差,且被識(shí)別為裸露地的像元中可能含有植被或建筑物信息,模型中植被覆蓋因子取值有待評(píng)估.

(2)揚(yáng)塵排放系數(shù)的計(jì)算過程中由于缺乏實(shí)測(cè)數(shù)據(jù)或驗(yàn)證數(shù)據(jù)帶來(lái)的誤差.具體表現(xiàn)為PMi百分比含量采用的是推薦值,而非動(dòng)力學(xué)粒徑譜儀的實(shí)測(cè)值.氣候因子計(jì)算采用的是氣象站點(diǎn)的數(shù)據(jù),誤差相對(duì)較小,但缺乏空間連續(xù)性;而空間插值方法也會(huì)帶來(lái)不確定性.

(3)土壤風(fēng)蝕指數(shù)Iwe的測(cè)定僅限于一時(shí)一地,其測(cè)量值可能隨時(shí)間、季節(jié)變化,對(duì)揚(yáng)塵排放估算結(jié)果造成較大影響.

(4)模型本身也存在一定的待改進(jìn)空間.因?yàn)槭軞夂驐l件,如風(fēng)速和降水的影響,風(fēng)蝕作用主要發(fā)生在冬春季[12],而模型中采用的氣候參數(shù)均為年均值,弱化了自然過程原本顯著的季節(jié)變化,這將導(dǎo)致裸露地?fù)P塵源PM10、PM2.5排放量的低估.

3.5 季節(jié)、月份計(jì)算的揚(yáng)塵排放量

針對(duì)模型中氣候參數(shù)采用年均值將導(dǎo)致顆粒物排放低估的問題,本研究進(jìn)一步對(duì)風(fēng)蝕模型進(jìn)行改進(jìn):①采用月均值和季度均值取代式(7)～(9)中相應(yīng)的年均值;②修訂式(9)中的時(shí)間項(xiàng),用對(duì)應(yīng)的天數(shù)替代;③考慮到冬季低溫導(dǎo)致E*趨近于0,故對(duì)(9)的Ta項(xiàng)進(jìn)行截?cái)嗵幚?設(shè)置最小溫度閾值為-2℃;④采用2010～2016年各區(qū)縣裸露地面積分布的平均值,乘以對(duì)應(yīng)逐月、季度的揚(yáng)塵排放系數(shù)并累計(jì),最終計(jì)算得到各月、各季度的裸露地?fù)P塵源PM10排放量.

結(jié)果顯示,分季節(jié)計(jì)算,北京市平原區(qū)裸露地?fù)P塵源的PM10年排放量可達(dá)39294t,逐月計(jì)算的PM10年排放量可達(dá)55175t,均遠(yuǎn)大于表5的結(jié)果7591.7t.這表明現(xiàn)有研究及模型可能忽視了裸露地?fù)P塵源揚(yáng)塵排放的季節(jié)變化,采用平均的年值氣候因子可能導(dǎo)致了揚(yáng)塵排放量的系統(tǒng)性低估.

本文完成了北京市平原區(qū)裸露地風(fēng)蝕揚(yáng)塵源顆粒物排放清單的初步編制工作,但土壤揚(yáng)塵源顆粒物排放清單的編制工作還需進(jìn)一步加強(qiáng).尤其是需要針對(duì)風(fēng)蝕作用的季節(jié)性變化做進(jìn)一步的詳細(xì)研究,并利用空氣質(zhì)量模型分析土壤揚(yáng)塵源對(duì)北京市大氣環(huán)境空氣質(zhì)量的影響,以探索適合北京市實(shí)際情況的揚(yáng)塵防治策略.

4 結(jié)論

4.1 北京市平原區(qū)的年均裸露地面積約為3587km2,不同季節(jié)存在較大差異,1～3月裸露地面積最大,約為4500km2,8月裸露地面積最小,約為500k

4.2 1987～2016年間,北京市平原區(qū)的裸露地面積減少了約600km2,據(jù)推測(cè)主要是因城市發(fā)展和城區(qū)擴(kuò)張?jiān)斐?

4.3 區(qū)域分布上,裸露地主要分布在城市發(fā)展新區(qū)和生態(tài)涵養(yǎng)發(fā)展區(qū),分別占60%和34%,其中,大興區(qū)和順義區(qū)的裸露地面積均超過550km2,通州區(qū)和延慶縣次之,首都核心區(qū)僅存在零星的裸露地分布.

4.4 以氣候因子的年均值計(jì)算,北京市平原區(qū)的土壤揚(yáng)塵PM10排放量為7591t/a,PM2.5排放量為1265t/a.

4.5 對(duì)風(fēng)蝕模型進(jìn)行改進(jìn),分別用逐月和季度氣候參數(shù)進(jìn)行計(jì)算,PM10排放量分別為55175t和39294t.這一結(jié)果表明,采用平均的年值氣候因子可能導(dǎo)致了揚(yáng)塵排放量的系統(tǒng)性低估.土壤風(fēng)蝕揚(yáng)塵對(duì)大氣顆粒物的貢獻(xiàn)可能大于預(yù)期.

[1] Yang F, Tan J, Zhao Q, et al. Characteristics of PM2.5speciation in representative megacities and across China [J]. Atmospheric Chemistry & Physics, 2011,11(11):1025-1051.

[2] Hallquist M, Wenger J C, Baltensperger U, et al. The formation,properties and impact of secondary organic aerosol: current and emerging issues [J]. Atmospheric Chemistry and Physics, 2009,9(14):5155-5236.

[3] 潘月云,李 楠,鄭君瑜,等.廣東省人為源大氣污染物排放清單及特征研究 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2015,35(9):2655-2669.

[4] 韓力慧,張 鵬,張海亮,等.北京市大氣細(xì)顆粒物污染與來(lái)源解析研究 [J]. 2016,36(11):3203-3210.

[5] 徐媛倩,姜 楠,燕啟社,等.鄭州市裸露地面風(fēng)蝕揚(yáng)塵排放清單研究 [J]. 環(huán)境污染與防治, 2016,38(4):22-27.

[6] 胡 敏,唐 倩,彭劍飛,等.我國(guó)大氣顆粒物來(lái)源及特征分析[J]. 環(huán)境與可持續(xù)發(fā)展, 2011,36(5):15-19.

[7] 陳 莉,李 濤,韓婷婷,等.WEPS模型下天津郊區(qū)風(fēng)蝕塵對(duì)城區(qū)空氣質(zhì)量的影響 [J]. 中國(guó)環(huán)境科學(xué), 2012,32(8):1353-1360.

[8] 北京市環(huán)境保護(hù)局.北京市PM2.5來(lái)源解析正式發(fā)布[EB/OL].http://www.bjepb.gov.cn/bjhrb/xxgk/jgzn/jgsz/jjgjgszjzz/xcjyc/xw fb/607219/index.html.

[9] 天津市環(huán)境保護(hù)局.天津發(fā)布顆粒物源解析結(jié)果 [EB/OL].http://www.tjhb.gov.cn/news/news_headtitle/201410/t20141009_570.html.

[10] 中國(guó)環(huán)境監(jiān)測(cè)總站.石家莊市大氣污染源解析結(jié)果 [EB/OL].http://www.cnemc.cn/publish/totalWebSite/news/news_42659.html.

[11] 王社扣,王體健,石 睿,等.南京市不同類型揚(yáng)塵源排放清單估計(jì) [J]. 中國(guó)科學(xué)院大學(xué)學(xué)報(bào), 2014,31(3):351-359.

[12] 宣 捷.中國(guó)北方地面起塵總量分布 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2000,20(4):426-430.

[13] Gillette D A, Passi R. Modeling dust emission caused by wind erosion [J]. Journal of Geophysical Research Atmospheres, 1988,93(D11):14233-14242.

[14] Holdt J R V, Eckardt F D, Wiggs G F S. Landsat identifies aeolian dust emission dynamics at the landform scale [J]. Remote Sensing of Environment, 2017,198:229-243.

[15] 環(huán)境保護(hù)部.揚(yáng)塵源顆粒物排放清單編制技術(shù)指南(試行)[EB/OL]. http://www.zhb.gov.cn/gkml/hbb/bgg/201501/t20150107_293955.htm.

[16] 吳志杰,趙書河.基于TM圖像的“增強(qiáng)的指數(shù)型建筑用地指數(shù)”研究 [J]. 國(guó)土資源遙感, 2012,24(2):50-55.

[17] Rikimaru A. Landsat TM data processing guide for forest canopy density mapping and monitoring model [C]//ITTO workshop on utilization of remote sensing in site assessment and planning for rehabilitation of logged-over forest, Bangkok, Thailand, 1996:1-8.

[18] 徐涵秋.利用改進(jìn)的歸一化差異水體指數(shù)(MNDWI)提取水體信息的研究 [J]. 遙感學(xué)報(bào), 2005,9(5):589-595.

[19] 范九倫,趙 鳳.灰度圖像的二維Otsu曲線閾值分割法 [J]. 電子學(xué)報(bào), 2007,35(4):751-755.

[20] 孫 璐,陳洪海.最大類間方差法在圖像分割中的應(yīng)用 [J]. 煤炭技術(shù), 2008,27(7):144-145.

[21] 張文娟,潘曉嵐.基于灰度直方圖的閾值分割算法分析與比較[J]. 科技資訊, 2006,(14):12-13.

[22] Huete A R. A soil-adjusted vegetation index (SAVI). Remote sensing of environment [J]. Remote Sensing of Environment,1988,25(3):295-309.

[23] Fryrear D W, Chen W N, Lester C. Revised wind erosion equation[J]. Annals of Arid Zone, 2001,40(3):265-279.

[24] 彭應(yīng)登,周 莉,田 剛.北京市裸露農(nóng)田可風(fēng)蝕性分析 [C]//北京綠色奧運(yùn)環(huán)境保護(hù)技術(shù)與發(fā)展研討會(huì), 2005:346-351.

[25] 唐新明,劉 浩,李 京,等.北京地區(qū)霾/顆粒物污染與土地利用/覆蓋的時(shí)空關(guān)聯(lián)分析 [J]. 中國(guó)環(huán)境科學(xué), 2015,35(9):2561-2569.

[26] 聶 磊,邵 霞,樊守彬,等.保護(hù)性耕作對(duì)減輕農(nóng)田風(fēng)蝕揚(yáng)塵PM10排放的作用 [C]//全國(guó)保護(hù)性耕地與農(nóng)機(jī)農(nóng)藝結(jié)合技術(shù)交流研討會(huì), 2009:181-183.

[27] 趙文慧,姜 磊,張立坤,等.北京平原區(qū)平房冬季燃煤量及污染物排放估算 [J]. 中國(guó)環(huán)境科學(xué), 2017,37(3):859-867.

[28] 徐 謙,李令軍,趙文慧,等.北京市建筑施工裸地的空間分布及揚(yáng)塵效應(yīng) [J]. 中國(guó)環(huán)境監(jiān)測(cè), 2015,(5):78-85.

[29] 樊守彬,張東旭,田靈娣.AP-42道路交通揚(yáng)塵排放模型評(píng)估及其在北京市的應(yīng)用 [J]. 環(huán)境工程學(xué)報(bào), 2016,10(5):2501-2506.

[30] 中國(guó)環(huán)境科學(xué)研究院.北京市可吸入顆粒物污染源信息平臺(tái)構(gòu)建和示范研究 [R]. 科技基礎(chǔ)性工作和社會(huì)公益研究專項(xiàng),2007,199-201.

[31] 葉芝祥,楊 松,楊懷金,等.成都市裸地?fù)P塵源排放因子研究[C]//葉芝祥.2015年中國(guó)環(huán)境科學(xué)學(xué)會(huì)學(xué)術(shù)年會(huì)論文集(第二卷), 2015:1612-1615.