徐媛倩 姜 楠 燕啟社 張瑞芹# 陳良富 李莘莘

(1.鄭州大學(xué)化學(xué)與分子工程學(xué)院，河南 鄭州 450001；2.中國科學(xué)院遙感與數(shù)字地球研究所，遙感科學(xué)國家重點(diǎn)實(shí)驗(yàn)室，北京 100101)

根據(jù)中國環(huán)境保護(hù)部發(fā)布的74個(gè)城市空氣質(zhì)量狀況，鄭州市2013年、2014年均位列空氣質(zhì)量相對(duì)較差的前十名城市名單中。2013年鄭州市環(huán)境空氣質(zhì)量優(yōu)良天數(shù)為127 d，僅占總天數(shù)的34.79%，嚴(yán)重污染天數(shù)16 d，重度污染天數(shù)56 d，PM2.5年均值為102 μg/m3。根據(jù)鄭州市PM2.5源解析結(jié)果，鄭州市PM2.5主要污染源為煤炭燃燒、機(jī)動(dòng)車污染和揚(yáng)塵污染3大類，其中揚(yáng)塵污染對(duì)大氣環(huán)境PM2.5的貢獻(xiàn)率為26%[1]。PM10受揚(yáng)塵污染影響更大，根據(jù)BI等[2]的研究，濟(jì)南、烏魯木齊、太原、天津、安陽、銀川等6個(gè)城市的城市揚(yáng)塵對(duì)大氣環(huán)境PM10的貢獻(xiàn)率分別為30%、30%、32%、34%、37%、59%。

建立揚(yáng)塵污染源排放清單是區(qū)域制霾必不可少的工作內(nèi)容之一，揚(yáng)塵污染源清單不僅為地方揚(yáng)塵控制提供科學(xué)依據(jù)，而且為區(qū)域空氣質(zhì)量評(píng)估提供數(shù)據(jù)支撐[3-4]。揚(yáng)塵污染源排放清單包括土壤揚(yáng)塵排放清單、道路揚(yáng)塵排放清單、建筑揚(yáng)塵排放清單和堆場揚(yáng)塵排放清單4個(gè)部分[5-7]。土壤揚(yáng)塵是指直接來源于裸露地面的揚(yáng)塵[8]，裸露地面起塵機(jī)制主要有滾動(dòng)、躍遷和懸浮3種形式，起塵過程既受風(fēng)場的動(dòng)力學(xué)影響，也受顆粒物粒徑影響[9]。土壤揚(yáng)塵對(duì)城市揚(yáng)塵污染具有相當(dāng)高的貢獻(xiàn)，由于其起動(dòng)風(fēng)速小，傳輸距離遠(yuǎn)，一定程度上加劇了大氣顆粒物污染[10-11]；方小珍等[12]對(duì)寧波市城市揚(yáng)塵源解析進(jìn)行研究，結(jié)果表明土壤揚(yáng)塵對(duì)城市揚(yáng)塵的貢獻(xiàn)率為34.9%，是貢獻(xiàn)率最大的揚(yáng)塵源。

土壤揚(yáng)塵包括裸露地面風(fēng)蝕揚(yáng)塵和農(nóng)耕地?fù)P塵兩類，其中裸露地面風(fēng)蝕揚(yáng)塵對(duì)土壤揚(yáng)塵的貢獻(xiàn)率為88%，是主要的土壤揚(yáng)塵貢獻(xiàn)源[13]871，故本研究以裸露地面風(fēng)蝕揚(yáng)塵為研究對(duì)象，以2013年為基準(zhǔn)年，參照環(huán)境保護(hù)部發(fā)布的《揚(yáng)塵源顆粒物排放清單編制技術(shù)指南(試行)》(以下簡稱《指南》)和文獻(xiàn)[14]、[15]的研究成果，結(jié)合遙感信息，利用排放系數(shù)法開展鄭州市裸露地面風(fēng)蝕揚(yáng)塵排放清單的首次編制工作，估算裸露地面風(fēng)蝕揚(yáng)塵的排放量，并基于地理信息系統(tǒng)(GIS)對(duì)顆粒物排放量進(jìn)行空間分布[16]，研究裸露地面風(fēng)蝕揚(yáng)塵的分布特征，為鄭州市揚(yáng)塵污染控制提供理論依據(jù)。

1 材料與方法

將鄭州市7個(gè)地區(qū)(鞏義、滎陽、登封、市區(qū)、中牟、新鄭、新密)劃分成n個(gè)1 km×1 km的網(wǎng)格，利用式(1)計(jì)算裸露地面風(fēng)蝕揚(yáng)塵的年排放量。

由式(1)可知，裸露地面風(fēng)蝕揚(yáng)塵年排放量的估算需要裸露地面面積和排放系數(shù)兩個(gè)參數(shù)，其中裸露地面的面積根據(jù)環(huán)境一號(hào)衛(wèi)星遙感數(shù)據(jù)和Google Earth提供的衛(wèi)星影像獲取，風(fēng)蝕揚(yáng)塵的排放系數(shù)通過計(jì)算得到。

1.1 裸露地面面積的獲取

本研究將環(huán)境一號(hào)衛(wèi)星影遙感像作為控制影像，采用目視識(shí)別為主，機(jī)助分類為輔的方法，勾繪裸露地面界線，繪制鄭州市裸露地面的分布圖[17]，得到裸露地面分布、面積等信息，具體方法如下：由于風(fēng)蝕揚(yáng)塵排放具有明顯的季節(jié)性差異，本研究選取質(zhì)量較好的環(huán)境一號(hào)衛(wèi)星遙感影像，以4月代表春季地表土地利用類型，8月代表夏季地表土地利用類型，用ENVI軟件對(duì)其進(jìn)行波段合成，得到研究區(qū)域真彩圖(見圖1)。由圖1可見，鄭州市8月與4月的植被覆蓋情況稍有差異，但部分區(qū)域兩個(gè)月份都是裸露地面，本研究將兩個(gè)季節(jié)都為裸露地面的區(qū)域認(rèn)定為永久性裸露地面，并對(duì)其產(chǎn)生的風(fēng)蝕揚(yáng)塵進(jìn)行研究。對(duì)衛(wèi)星數(shù)據(jù)進(jìn)行波段計(jì)算，得到歸一化植被指數(shù)(NDVI)分布情況，為計(jì)算方便，將NDVI計(jì)算結(jié)果放大1 000倍，使NDVI的取值在-200～700(見圖2)，根據(jù)NDVI取值初步判定土地利用類型。當(dāng)NDVI取值為負(fù)時(shí)，表示地面覆蓋為云、水、雪等，NDVI取值為0時(shí)表示地面有巖石或裸土等存在，當(dāng)NDVI取值為正時(shí)，表示地面有不同程度的植被覆蓋，NDVI取值在兩個(gè)季節(jié)均為0～150時(shí)可認(rèn)定為永久性揚(yáng)塵源[13]880。通過對(duì)真彩圖進(jìn)行非監(jiān)督分類，得到進(jìn)一步的土地利用類型分類，在ArcGIS軟件中以真彩圖為控制影像，以真彩圖信息為主，結(jié)合NDVI取值和非監(jiān)督分類結(jié)果對(duì)4月和8月均存在的裸露地面進(jìn)行勾繪，并參照Google Earth中的高分辨率衛(wèi)星影像對(duì)勾繪結(jié)果進(jìn)行調(diào)整，得到穩(wěn)定存在的裸露地面分布數(shù)據(jù)，再在ArcGIS軟件中對(duì)裸露地面分布圖層和1 km×1 km 網(wǎng)格圖層進(jìn)行處理，最終得到每個(gè)網(wǎng)格中裸露地面的面積，以便做后續(xù)的空間分配。經(jīng)分析，鄭州市裸露地面主要分布在北部、西部，北部為黃河灘涂，西部為嵩山地區(qū)裸露山體，其他區(qū)域的裸露地面主要是中東部建筑施工、土方開挖等形成的空地。

圖1 鄭州市真彩圖Fig.1 True color map of Zhengzhou

圖2 鄭州市NDVI取值分布Fig.2 NDVI distribution of Zhengzhou

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

本研究針對(duì)風(fēng)蝕揚(yáng)塵中3種顆粒物(PM2.5、PM10、總懸浮顆粒物(TSP))在3種類型裸露地面(灘涂、裸露山體、空地)的排放系數(shù)進(jìn)行計(jì)算，計(jì)算公式如下：

E=D×C×(1-η)×100

(2)

式中：E為裸露地面風(fēng)蝕揚(yáng)塵中顆粒物的排放系數(shù)，t/(km2·a)；D為顆粒物的起塵因子，t/(hm2·a)；C為氣候因子；η為污染控制技術(shù)對(duì)風(fēng)蝕揚(yáng)塵的去除率。

起塵因子、氣候因子計(jì)算的經(jīng)驗(yàn)公式分別見式(3)、式(4)：

D=K×Iwe××L×V

(3)

C=0.504×u3/PE2

(4)

式中：K為顆粒物在風(fēng)蝕揚(yáng)塵中的比例；Iwe為土壤風(fēng)蝕指數(shù)，t/(hm2·a)；為地面粗糙因子；L為無屏蔽寬度因子；V為植被覆蓋因子；u為年平均風(fēng)速，m/s；PE為桑氏威特降水指數(shù)(蒸發(fā)指數(shù))，PE計(jì)算的經(jīng)驗(yàn)公式如下：

PE=100×(p/E*)

(5)

E*=(0.594 9+(0.118 9×Ta))×365

(6)

式中：p為年降水量，mm；E*為年潛在蒸發(fā)量，mm；Ta為年平均溫度，℃。

由式(2)可知，排放系數(shù)的決定因素是起塵因子、氣候因子和污染控制技術(shù)對(duì)風(fēng)蝕揚(yáng)塵的去除率，由于本研究中只考慮裸露地面的風(fēng)蝕揚(yáng)塵排放量，設(shè)裸露地面污染控制技術(shù)對(duì)風(fēng)蝕揚(yáng)塵的去除率為0，此時(shí)風(fēng)蝕揚(yáng)塵的排放系數(shù)僅取決于起塵因子和氣候因子兩個(gè)因素。

起塵因子的計(jì)算需要顆粒物在風(fēng)蝕揚(yáng)塵中的比例、土壤風(fēng)蝕指數(shù)、地面粗糙因子、無屏蔽寬度因子和植被覆蓋因子5個(gè)參數(shù)。其中，顆粒物在風(fēng)蝕揚(yáng)塵中的比例和地面粗糙因子采用《指南》推薦值；植被覆蓋因子指實(shí)際裸露地面面積與總計(jì)算面積的比值，本研究以勾繪出的裸露地面面積為總計(jì)算面積，考慮到勾繪出的裸露地面上有少量植物存在，故將植被覆蓋因子取值為0.9；無屏蔽寬度因子指沒有明顯阻擋物的最大范圍，當(dāng)無屏蔽寬度<300 m時(shí)，L=0.70，當(dāng)無屏蔽寬度在300～600 m時(shí)，L=0.85，當(dāng)無屏蔽寬度≥600 m時(shí)，L=1.00。參考文獻(xiàn)[18]、[19]并結(jié)合實(shí)地調(diào)查，認(rèn)為黃河灘涂土質(zhì)為砂土，山體土質(zhì)為黏壤土，空地土質(zhì)為壤質(zhì)砂土，根據(jù)《指南》推薦方法計(jì)算3種顆粒物在不同裸露地面的土壤風(fēng)蝕指數(shù)。起塵因子各計(jì)算參數(shù)的取值見表1。

氣候因子計(jì)算需要年平均溫度、年降水量和年平均風(fēng)速3個(gè)參數(shù)，根據(jù)鄭州市7個(gè)地區(qū)2013年的氣象自動(dòng)站觀測數(shù)據(jù)，得到各地區(qū)的氣象因子計(jì)算參數(shù)，具體見表2。

2 結(jié)果與討論

2.1 裸露地面的分布及面積

鄭州市主要土地利用類型為水體、自然綠地、農(nóng)用地、城鎮(zhèn)建筑用地和交通用地等，其中鄭州市北邊邊界線以南為黃河，是水體的主要分布區(qū)域，自然綠地主要分布在西部嵩山地區(qū)，農(nóng)用地主要分布在鄭州市中部和東部的非城鎮(zhèn)地區(qū)[20]。圖3為本研究勾繪出的鄭州市裸露地面分布情況。利用ArcGIS軟件對(duì)圖3中裸露地面的面積進(jìn)行計(jì)算，得到鄭州市裸露地面面積約208 km2，占鄭州市總面積的2.8%。

表1 起塵因子計(jì)算參數(shù)

表2 7個(gè)地區(qū)2013年氣象參數(shù)

圖3 鄭州市裸露地面分布Fig.3 Distributing of barrens in Zhengzhou

2.2 排放系數(shù)

排放系數(shù)由起塵因子和氣候因子共同決定，反映了單位面積裸露地面的起塵能力。由于土質(zhì)類型的差異，不同類型裸露地面的起塵因子存在一定的差別。各類型裸露地面的起塵因子計(jì)算結(jié)果如圖4所示。從圖4可以看出，3種類型裸露地面的起塵因子排序?yàn)闉┩?空地>裸露山體。

鄭州市7個(gè)地區(qū)的氣候因子計(jì)算結(jié)果見圖5。

圖4 不同裸露地面的起塵因子Fig.4 Particulate emission factor of barrens with different prefecture

從圖5可以看出，7個(gè)地區(qū)中新鄭的氣候因子最大，市區(qū)次之，新密氣候因子最小。氣候因子對(duì)顆粒物排放系數(shù)的影響主要來源于年平均風(fēng)速，這是因?yàn)轱L(fēng)速越大揚(yáng)起的顆粒物越多，粒徑范圍越廣，當(dāng)風(fēng)速為2.0～3.5 m/s時(shí)，粒徑在60 μm以下的顆粒物可被揚(yáng)起；當(dāng)風(fēng)速為1.0～2.0 m/s時(shí)，粒徑在20 μm以下的顆粒物可被揚(yáng)起。

圖5 鄭州市7個(gè)地區(qū)的氣候因子Fig.5 Climate factor of 7 district in Zhengzhou

利用氣候因子和起塵因子計(jì)算鄭州市各地區(qū)不同類型裸露地面的顆粒物排放系數(shù)，結(jié)果見圖6。鄭州市裸露地面風(fēng)蝕揚(yáng)塵中PM2.5、PM10、TSP的排放系數(shù)計(jì)算結(jié)果見圖7。

圖6 鄭州市不同地區(qū)裸露地面的排放系數(shù)Fig.6 Emission factor of barrens with different prefecture in different districts of Zhengzhou

由圖6可見，鄭州市各地區(qū)顆粒物排放系數(shù)存在較大差異，3種裸露地面類型下均為新鄭、市區(qū)的顆粒物排放系數(shù)最大，登封次之，滎陽、新密相對(duì)最小。由圖7可見，鄭州市裸露地面風(fēng)蝕揚(yáng)塵中PM2.5、PM10、TSP的排放系數(shù)分別為(3.36±2.24)、(20.16±13.44)、(67.21±44.81) t/(km2·a)，與2002年江小珂等[21]對(duì)北京市的研究結(jié)果相近(風(fēng)蝕揚(yáng)塵中PM2.5、PM10排放系數(shù)分別為3.31、22.7 t/(km2·a))，說明我國北方城市在城市發(fā)展過程中存在著水平相當(dāng)?shù)穆懵兜孛骘L(fēng)蝕揚(yáng)塵顆粒物污染問題。

圖7 3種顆粒物的排放系數(shù)Fig.7 Emission factor of 3 types particulate matters

2.3 顆粒物排放特征

鄭州市7個(gè)地區(qū)的裸露地面風(fēng)蝕揚(yáng)塵顆粒物排放量計(jì)算結(jié)果見圖8。由圖8可見，7個(gè)地區(qū)顆粒物排放量存在較大差異，具體表現(xiàn)為市區(qū)>中牟>登封>鞏義>滎陽>新鄭>新密。經(jīng)計(jì)算，2013年鄭州市裸露地面風(fēng)蝕揚(yáng)塵源PM2.5、PM10、TSP排放量分別為597、3 581、11 937 t。

圖8 鄭州市7個(gè)地區(qū)風(fēng)蝕揚(yáng)塵中顆粒物的排放量Fig.8 Emission load of particulate matters from wind erosion dust in 7 district of Zhengzhou

以1 km×1 km的網(wǎng)格為研究單位，通過對(duì)上述各階段研究成果進(jìn)行總結(jié)分析，得到鄭州市裸露地面風(fēng)蝕揚(yáng)塵源顆粒物排放量的空間分布，結(jié)果見圖9。鑒于TSP分布情況與PM2.5、PM10總體相似，且PM2.5、PM10為顆粒物研究的重點(diǎn)，本研究僅列出鄭州市PM2.5、PM10的排放量分布情況。

從圖9可以看出，鄭州市裸露地面PM2.5和PM10的高排區(qū)主要集中在研究區(qū)北部邊緣(黃河流域的灘涂區(qū)域)、鞏義和登封兩市的山地區(qū)域及市區(qū)北部和中牟東北部地區(qū)。鞏義和登封兩地的裸露地面以裸露山體為主，該區(qū)域植被稀疏，海拔較高，易受自然風(fēng)的影響而起塵?？刂圃搮^(qū)域的風(fēng)蝕揚(yáng)塵起塵量，需從山體綠化做起，減少裸露山體的面積；市區(qū)和中牟裸露地面PM2.5、PM10的高排區(qū)主要分布在市區(qū)北邊的黃河灘涂區(qū)域，該區(qū)域地表開闊、土質(zhì)疏松、植被稀疏，易受自然風(fēng)的影響而起塵。控制該區(qū)域的風(fēng)蝕揚(yáng)塵起塵量，要對(duì)黃河灘涂區(qū)域進(jìn)行加固和綠化，以減少灘涂裸露地面的暴露；其他地區(qū)裸露地面PM2.5、PM10的排放除分布在灘涂、裸露山體區(qū)域外，還分布在空地區(qū)域，該區(qū)域土質(zhì)疏松，地表無植被防護(hù)，自然風(fēng)吹后易起塵。對(duì)于建筑施工造成的人為空地，建筑工地及土方開挖過程應(yīng)著重加強(qiáng)防塵措施，減少施工過程中裸露地面的暴露；對(duì)于自然形成的空地，應(yīng)加大監(jiān)管和綠化力度，盡量避免自然空地的形成。

圖9 2013年鄭州市裸露地面PM2.5、PM10排放量分布Fig.9 Distributing characteristic of PM2.5，PM10 from barrens in Zhengzhou,2013

2.4 不確定性分析

在顆粒物排放量估算過程中，排放因子、活動(dòng)水平等因素給計(jì)算結(jié)果帶來一定不確定性。具體涉及的不確定性因素主要來自以下兩個(gè)方面：(1)在目視識(shí)別裸露地面獲取污染源活動(dòng)水平過程中，人為識(shí)別存在一定的誤差，衛(wèi)星影像分辨率的限制進(jìn)一步增大了目視識(shí)別的誤差；(2)排放因子計(jì)算過程中由于缺少實(shí)測數(shù)據(jù)或數(shù)據(jù)來自間接推估造成的誤差，具體表現(xiàn)在起塵因子計(jì)算過程中的誤差和氣候因子空間分布的誤差。起塵因子計(jì)算過程中顆粒物在風(fēng)蝕揚(yáng)塵中的比例采用《指南》推薦值，土壤風(fēng)蝕指數(shù)為文獻(xiàn)中的推估數(shù)值，上述參數(shù)的不確定性在一定程度上影響排放系數(shù)的計(jì)算結(jié)果。氣候因子計(jì)算采用鄭州市7個(gè)地區(qū)自動(dòng)氣象站數(shù)據(jù)，誤差相對(duì)較小，但是在空間分配時(shí)一個(gè)行政地區(qū)采用一個(gè)監(jiān)測點(diǎn)的數(shù)據(jù)，無空間連續(xù)性，這對(duì)排放量計(jì)算結(jié)果和空間分配都會(huì)造成一定的誤差。

3 結(jié)論與建議

(1) 鄭州市2013年裸露地面風(fēng)蝕揚(yáng)塵源PM2.5、PM10、TSP排放系數(shù)分別為(3.36±2.24)、(20.16±13.44)、(67.21±44.81) t/(km2·a)；PM2.5、PM10、TSP年排放量分別為597、3 581、11 937 t；裸露地面風(fēng)蝕揚(yáng)塵源顆粒物的排放主要分布在鄭州市北部黃河灘涂，西部嵩山裸露山體、中東部建筑施工、土方開挖等形成的空地，減少鄭州市風(fēng)蝕揚(yáng)塵起塵量的主要措施是加大綠化和管控力度。

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

(致謝：本研究是在中國科學(xué)院遙感與數(shù)字地球研究所,遙感科學(xué)國家重點(diǎn)實(shí)驗(yàn)室各位老師和同學(xué)的幫助下完成，在此表示感謝！)

[1] GENG Ningbo,WANG Jia,XU Yifei,et al.PM2.5in an industrial district of Zhengzhou,China:chemical composition and source apportionment[J].Particuology,2013,11(1):99-109.

[2] BI Xiaohui,FENG Yinchang,WU Jianhui,et al.Source apportionment of PM10in six cities of northern China[J].Atmospheric Environment,2007,41(5):903-912.

[3] 趙斌，馬建中．天津市大氣污染源排放清單的建立[J].環(huán)境科學(xué)學(xué)報(bào)，2008，28(2)：368-375．

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

[5] PARK S H,GONG S L,GONG W,et al.Relative impact of windblown dust versus anthropogenic fugitive dust in PM2.5on air quality in North America [J].Journal of Geophysical Research,2010,115(D16):751-763.

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

[7] MARTINS H,MIRANDA A L,BORREGO C,et al.Air pollution modeling and its application ⅩⅪ [M].Berlin: Springer Netherlands,2012.

[8] HJ/T 393—2007，防治城市揚(yáng)塵污染技術(shù)規(guī)范[S].

[9] 韓旸，白志鵬，姬亞芹，等．裸土風(fēng)蝕型開放源起塵機(jī)制研究進(jìn)展[J].環(huán)境污染與防治，2008，30(2)：77-82．

[10] 張利文，白志鵬，郭光煥，等．開放源對(duì)環(huán)境空氣質(zhì)量影響的評(píng)估技術(shù)與實(shí)例[J].環(huán)境科學(xué)研究，2006，19(3)：18-23．

[11] TAO Minghui,CHEN Liangfu,SU Lin,et al.Satellite observation of regional haze pollution over the North China Plain[J].Journal of Geophysical Research,2012,117(D12):48-50.

[12] 方小珍，孫列，畢曉輝，等．寧波城市揚(yáng)塵化學(xué)組成特征及其來源解析[J].環(huán)境污染與防治，2014，36(1)：55-59．

[13] KIM D,CHIN M,BIAN Huisheng,et al.The effect of the dynamic surface bareness on dust source function,emission,and distribution[J].Journal of Geophysical Research,2013,118(2).

[14] WOODRUFF N,SIDDOWAY F.A wind erosion equation[J].Soil Science Society of America Journal,1965,29(5):602-608.

[15] XUAN Jie.Dust emission factors for environment of Northern China[J].Atmospheric Environment,1999,33(11):1767-1776.

[16] 鄭君瑜，王水勝，黃志炯，等．區(qū)域高分辨率大氣排放源清單建立的技術(shù)方法與應(yīng)用[M].北京：科學(xué)出版社，2014．

[17] 趙庚星，李玉環(huán)，徐春達(dá)．遙感和GIS支持的土地利用動(dòng)態(tài)監(jiān)測研究——以黃河三角洲墾利縣為例[J].應(yīng)用生態(tài)學(xué)報(bào)，2000，11(4)：573-576．

[18] 鐘國敏，張學(xué)雷，段金龍，等．鄭州市土壤多樣性和土地利用多樣性研究及關(guān)聯(lián)分析[J].土壤通報(bào)，2013，44(3)：513-521．

[19] 環(huán)境保護(hù)部.揚(yáng)塵源顆粒物排放清單編制技術(shù)指南(試行)[R].北京：環(huán)境保護(hù)部，2015.

[20] 段金龍，張學(xué)雷，王玉龍．鄭州市土地資源動(dòng)態(tài)變化的遙感監(jiān)測與分析[J].河南科技大學(xué)學(xué)報(bào)(自然科學(xué)版)，2013，34(1)：98-100．

[21] 江小珂，唐孝炎．北京市大氣污染控制對(duì)策研究[Z].北京：北京市環(huán)境保護(hù)局，2002.