魏夜香,張霄羽,張 紅

中國(guó)二氧化硫的時(shí)空分布及主要排放來(lái)源研究

魏夜香,張霄羽*,張 紅

(山西大學(xué)環(huán)境與資源學(xué)院,山西 太原 030006)

隨著城市化進(jìn)程的加快和經(jīng)濟(jì)的快速發(fā)展,中國(guó)已成為世界第三大二氧化硫(SO2)排放國(guó),人類(lèi)健康和社會(huì)經(jīng)濟(jì)的可持續(xù)發(fā)展面臨著嚴(yán)峻的挑戰(zhàn).基于OMSO2e產(chǎn)品分析了2005～2020年中國(guó)區(qū)域大氣SO2柱濃度的時(shí)空格局變化,并利用MEIC排放數(shù)據(jù)和社會(huì)經(jīng)濟(jì)數(shù)據(jù)探討影響中國(guó)10大經(jīng)濟(jì)分區(qū)SO2變化的主要因素.結(jié)果表明:16年間SO2呈現(xiàn)波動(dòng)下降的趨勢(shì),新疆和青藏高原地區(qū)略微增長(zhǎng),其余地區(qū)均呈現(xiàn)顯著降低的趨勢(shì).SO2年際變化較大,整體呈現(xiàn)東高西低的空間分布格局,北部沿海和黃河中游地區(qū)年均值和變異系數(shù)最高.SO2與工業(yè)源、居民源、電力源、交通源及第二產(chǎn)業(yè)占GRP比重呈現(xiàn)正相關(guān),與人口密度、工業(yè)廢氣治理設(shè)施和工業(yè)廢氣治理費(fèi)用呈現(xiàn)負(fù)相關(guān).排放源的變化直接影響了SO2的變化.本研究為各地區(qū)制定相關(guān)的減排政策和生態(tài)環(huán)境保護(hù)提供了理論依據(jù)和參考.

二氧化硫柱濃度；時(shí)空分布；排放源；隨機(jī)森林；影響因素

隨著氣候變暖的加劇,氣候變化已成為21世紀(jì)各國(guó)探討的熱門(mén)話(huà)題.目前我們不僅面臨著氣候變化的嚴(yán)峻挑戰(zhàn),同時(shí)也需要應(yīng)對(duì)空氣污染造成的嚴(yán)重危害[1-2].二氧化硫(SO2)作為大氣污染物的一種,是一種無(wú)色且有刺激性氣味的有毒氣體,對(duì)人體健康和大氣環(huán)境造成巨大的影響.根據(jù)綠色和平環(huán)境信托組織(GET,2019)的報(bào)告,中國(guó)是全球第三大SO2排放國(guó),貢獻(xiàn)了約8%的全球人為SO2排放量,僅次于印度和俄羅斯[3].根據(jù)2021年中國(guó)環(huán)境狀況公報(bào),98.2%的城市SO2濃度達(dá)到國(guó)家一級(jí)標(biāo)準(zhǔn),酸雨區(qū)面積占國(guó)土總面積的3.8%,主要分布在長(zhǎng)江以南-云貴高原以東地區(qū)[4].

空氣污染物與溫室氣體具有同根、同源、同時(shí)的特性,化石燃料的燃燒、交通運(yùn)輸和居民生活等均是環(huán)境污染物與溫室氣體的主要來(lái)源[1-2].污染氣體的治理對(duì)改善大氣環(huán)境、推動(dòng)雙碳目標(biāo)具有重要的意義.針對(duì)大氣污染問(wèn)題,政府頒布并實(shí)施一系列的法律政策,如《大氣污染物綜合排放標(biāo)準(zhǔn)》、《大氣污染防治行動(dòng)計(jì)劃》和《京津冀及周邊地區(qū)大氣污染治理行動(dòng)計(jì)劃》等,用于控制SO2的排放,這些嚴(yán)格的排放標(biāo)準(zhǔn)和政策的實(shí)施對(duì)SO2排放的減少發(fā)揮了舉足輕重的作用.研究SO2濃度的時(shí)空分布及影響因素對(duì)政策實(shí)施的有效性及人類(lèi)健康有至關(guān)重要的意義.

認(rèn)識(shí)SO2的時(shí)空分布主要包括地面監(jiān)測(cè)和衛(wèi)星遙感兩種手段,地面監(jiān)測(cè)站可以提供近地表在線(xiàn)微量氣體濃度的監(jiān)測(cè),具有準(zhǔn)確度高,實(shí)時(shí)檢測(cè),數(shù)據(jù)可靠等特點(diǎn),但受到人為干擾較大且站點(diǎn)分布主要集中在城市區(qū)域.而衛(wèi)星遙感具有連續(xù)性強(qiáng),覆蓋范圍廣等優(yōu)勢(shì),成為探測(cè)污染氣體的重要手段[5].臭氧檢測(cè)儀(OMI)是研究長(zhǎng)時(shí)間大氣SO2變化的主要傳感器,可以實(shí)現(xiàn)大范圍,長(zhǎng)時(shí)間對(duì)SO2濃度的全球觀測(cè)[6].

目前關(guān)于SO2濃度的研究集中于探究其時(shí)空分布格局,驅(qū)動(dòng)因子和前后向軌跡[7-10].時(shí)空變化的研究從全球到區(qū)域均有涉及,以京津冀為研究區(qū)域較多.2021年zhai等人研究了二氧化硫在京津冀地區(qū)的分布和變化,指出SO2呈現(xiàn)下降趨勢(shì),主要?dú)w功于政府實(shí)施的管控措施[11].同時(shí),一些學(xué)者對(duì)二氧化硫的影響因素進(jìn)行了深入研究,研究表明SO2濃度的大小,一方面受到本地污染源排放的影響,另一方面與當(dāng)?shù)氐奶鞖?、地形等緊密相關(guān).例如:劉濤濤等人分析山西省二氧化硫的影響因素,得出人類(lèi)活動(dòng)對(duì)SO2的貢獻(xiàn)度較大,但氣象要素也發(fā)揮著舉足輕重的作用[12].魏毅[13]研究烏魯木齊的二氧化硫污染特征時(shí),指出燃煤、地形和氣象等因素均影響了SO2的濃度,特殊的盆地地形使大氣污染物不易擴(kuò)散和稀釋.關(guān)于SO2驅(qū)動(dòng)因素的研究大多從自然和人為兩方面進(jìn)行分析,考慮排放清單對(duì)大氣中SO2的影響研究較少.

本研究基于2005～2020年OMSO2e的數(shù)據(jù),在探討不同區(qū)域SO2的時(shí)空分布特征的基礎(chǔ)上,研究了MEIC排放數(shù)據(jù)(民用源、工業(yè)源、電力排放源和機(jī)動(dòng)車(chē)排放源)和社會(huì)經(jīng)濟(jì)(人口密度、GDP、夜間燈光)對(duì)SO2分布的影響,并借助隨機(jī)森林模型分析每個(gè)地區(qū)SO2濃度的主要影響因素,為各地區(qū)有效控制大氣污染,制定相關(guān)減排政策提供理論依據(jù)和參考.

1 數(shù)據(jù)和方法

1.1 研究區(qū)概況

根據(jù)經(jīng)濟(jì)發(fā)展、氣候和地形特點(diǎn),本文將中國(guó)大陸劃分為十個(gè)地區(qū)(圖1).分別為:東部沿海地區(qū)(EC)、珠江中上游地區(qū)(MUPR)、長(zhǎng)江中上游地區(qū)(MUYR)、黃河中游地區(qū)(MYR)、北部沿海地區(qū)(NC)、東北地區(qū)(NE)、青藏高原地區(qū)(QTP)、東南沿海地區(qū)(SC)、黃河上游地區(qū)(UYR)和新疆地區(qū)(XJ)[14-15].

圖1 中國(guó)十大區(qū)域分布示意

(審圖號(hào):GS(2020)4630號(hào))下同

1.2 數(shù)據(jù)來(lái)源

1.2.1 衛(wèi)星遙感數(shù)據(jù) 研究中SO2濃度數(shù)據(jù)來(lái)源于搭載在Aura衛(wèi)星的OMI傳感器反演得到的OMSO2e產(chǎn)品,其空間分辨率為0.25°,時(shí)間分辨率為1d,單位為多布森單位,簡(jiǎn)稱(chēng)DU(Dobson units).Aura衛(wèi)星于2004年7月15日發(fā)射升空,其軌道高度為705km,軌道周期為98.8min,在我國(guó)的過(guò)境時(shí)間為13:40～13:50,是一顆近極地、太陽(yáng)同步軌道衛(wèi)星. OMI可以監(jiān)測(cè)大氣中的O3、NO2、SO2、氣溶膠等氣體(https://daac.gsfc.nasa.gov/),已經(jīng)被廣泛應(yīng)用于污染氣體監(jiān)測(cè),空氣質(zhì)量預(yù)報(bào)和氣候變化研究等領(lǐng)域[16-17].本研究下載了2005年1月1日至2020年12月31日的OMSO2e產(chǎn)品,并計(jì)算了月及年SO2平均柱濃度.

1.2.2 地面站點(diǎn)數(shù)據(jù) 二氧化硫?qū)崪y(cè)數(shù)據(jù)來(lái)自中國(guó)科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心的城市空氣質(zhì)量監(jiān)測(cè)數(shù)據(jù)(https://www.resdc.cn/).在全國(guó)300多個(gè)城市,共設(shè)有1400多個(gè)監(jiān)測(cè)站點(diǎn),提供了PM2.5、O3、NO2、SO2、CO的質(zhì)量濃度[18-19].由于環(huán)境監(jiān)測(cè)站可獲取的時(shí)間較短,本文選取2015年1月1日到2020年12月31日每小時(shí)近地面的SO2濃度,監(jiān)測(cè)范圍覆蓋331個(gè)城市(含直轄市、地級(jí)市、地區(qū)、自治州和盟),作為對(duì)衛(wèi)星產(chǎn)品的驗(yàn)證數(shù)據(jù).

1.2.3 排放清單數(shù)據(jù) 中國(guó)多分辨率排放清單模型(MEIC)由清華大學(xué)基于云計(jì)算平臺(tái)開(kāi)發(fā)的中國(guó)大氣污染物和溫室氣體人為源排放清單模型,該清單主要包括溫室氣體和10種大氣污染物(SO2、NO、CO、NMVOC、NH3、CO2、PM2.5、PM10、BC和OC),并劃分了工業(yè)、居民、交通、農(nóng)業(yè)、電力五部門(mén)(http://meicmodel.org.cn)[18,20].本文利用2005～2020年中國(guó)大陸地區(qū)SO2各部門(mén)的年排放數(shù)據(jù),空間分辨率為0.25o′0.25o,分析SO2時(shí)空變化的潛在原因.

1.2.4 社會(huì)經(jīng)濟(jì)數(shù)據(jù) 社會(huì)經(jīng)濟(jì)數(shù)據(jù)包括:夜間燈光、GDP、人口密度、第二產(chǎn)業(yè)占GRP(地區(qū)生產(chǎn)總值)的比重、工業(yè)廢氣治理費(fèi)用和工業(yè)廢氣治理設(shè)施數(shù)(表1).由于SO2與各因子數(shù)據(jù)集的空間分辨率不同,借助ArcGIS重采樣功能,將其與SO2統(tǒng)一為相同的分辨率,以便進(jìn)一步的計(jì)算分析.

表1 社會(huì)經(jīng)濟(jì)數(shù)據(jù)

1.3 研究方法

1.3.1 相關(guān)分析 皮爾遜相關(guān)分析用來(lái)比較兩個(gè)變量之間的相關(guān)程度的強(qiáng)弱,其值位于-1到1之間.當(dāng)皮爾遜相關(guān)系數(shù)值大于零時(shí)表示正相關(guān),小于零時(shí)表示負(fù)相關(guān).其絕對(duì)值越接近1,表明相關(guān)性越強(qiáng);相反,兩個(gè)變量之間的相關(guān)性越弱[24].本文利用皮爾遜相關(guān)來(lái)分析SO2柱濃度和各影響因子之間的關(guān)系.其公式如下:

1.3.2 趨勢(shì)分析 趨勢(shì)分析法是對(duì)隨時(shí)間變化的變量進(jìn)行線(xiàn)性回歸分析,從而來(lái)預(yù)測(cè)其變化趨勢(shì)的方法.常用于大氣污染物的時(shí)間動(dòng)態(tài)變化[25],Tiao指出,當(dāng)p值>0.05或者在95%的置信區(qū)間時(shí),趨勢(shì)被認(rèn)為是顯著的[26].本文采用趨勢(shì)分析法探究在2005—2020年中國(guó)各區(qū)域SO2的時(shí)間變化特征.其計(jì)算方法如下:

式中:slope為多年每個(gè)柵格像元的傾向率,SO2i為第年的SO2的柵格像元值,為時(shí)間序列.若slope>0時(shí),表示該柵格多年SO2呈現(xiàn)上升趨勢(shì);當(dāng)Slope=0,表示該柵格多年SO2基本不變;當(dāng)Slope<0時(shí),表示該柵格多年SO2呈現(xiàn)下降趨勢(shì).

1.3.3 隨機(jī)森林 隨機(jī)森林模型(Random forest, RF)是由Breiman提出的一種基于決策樹(shù)的機(jī)器學(xué)習(xí)算法[27],既可用于構(gòu)建分類(lèi)模型,也可用于構(gòu)建回歸模型.該模型相比于傳統(tǒng)的回歸方法,具有高精度,能夠高效處理高維度數(shù)據(jù),不必?fù)?dān)心過(guò)擬合以及實(shí)現(xiàn)簡(jiǎn)單等優(yōu)點(diǎn).RF由多棵樹(shù)集成而成,具有二分叉特性,能夠提供特征重要性(feature importance)這一獨(dú)特排序?qū)傩?因此可以通過(guò)此屬性來(lái)判斷解釋變量對(duì)因變量的擬合程度[28].本研究通過(guò)建立SO2與影響因子的關(guān)系,獲取影響因子對(duì)SO2的重要性排序.

本文將影響因素的多年傾向率作為驅(qū)動(dòng)因子, SO2多年傾向率作為因變量,構(gòu)建隨機(jī)森林回歸模型．將原始數(shù)據(jù)按3:7的比例劃分訓(xùn)練集和測(cè)試集,通過(guò)調(diào)試最終選定的參數(shù)如下:隨機(jī)森林中樹(shù)的數(shù)目(n_tree)為1000,每個(gè)節(jié)點(diǎn)可能分裂的變量數(shù)(mtry)為4,探討近地表SO2與其影響因素之間的關(guān)系,進(jìn)而預(yù)測(cè)變量的相對(duì)重要性.

圖2 衛(wèi)星遙感與地面站點(diǎn)驗(yàn)證

(a)皮爾遜相關(guān)性空間分析圖(b) 石家莊市擬合圖

2 結(jié)果與討論

2.1 OMI衛(wèi)星數(shù)據(jù)的驗(yàn)證

研究統(tǒng)計(jì)了全國(guó)331個(gè)城市地面觀測(cè)的SO2,并與衛(wèi)星反演的SO2進(jìn)行了比較.圖2a給出了2015年1月—2020年12月中國(guó)OMI反演的SO2柱濃度與地面觀測(cè)站點(diǎn)月均值的相關(guān)系數(shù)空間分布圖[29].其中通過(guò)顯著性檢驗(yàn)的城市有96個(gè),相關(guān)性大于0.5的城市有102個(gè),主要分布在黃河中游、北部沿海和黃河上游的部分地區(qū),在西南部和新疆的極個(gè)別城市呈現(xiàn)負(fù)相關(guān)趨勢(shì).與此同時(shí),以濃度相對(duì)較高的河北省石家莊市為例,利用回歸分析法分析了行星邊界層SO2與地面監(jiān)測(cè)站SO2濃度的關(guān)系,圖2b為地面站點(diǎn)觀測(cè)的SO2濃度和衛(wèi)星遙感反演得到的SO2柱濃度的散點(diǎn)擬合圖,可以看出,擬合優(yōu)度(R2)達(dá)到60%,通過(guò)了95%的置信度檢驗(yàn),說(shuō)明衛(wèi)星觀測(cè)數(shù)據(jù)集基本能夠反映當(dāng)?shù)氐奈廴厩闆r[30].

圖3 各區(qū)域SO2濃度的時(shí)間變化趨勢(shì)(a)及中國(guó)SO2濃度傾向率的空間分布(b)

東部沿海地區(qū)(EC)、珠江中上游地區(qū)(MUPR)、長(zhǎng)江中上游地區(qū)(MUYR)、黃河中游地區(qū)(MYR)、北部沿海地區(qū)(NC)、東北地區(qū)(NE)、青藏高原地區(qū)(QTP)、東南沿海地區(qū)(SC)、黃河上游地區(qū)(UYR)和新疆(XJ)

2.2 中國(guó)二氧化硫的時(shí)空變化

2.2.1 時(shí)間變化特征 本文利用SO2隨年份變化的速率表征其多年時(shí)間變化.圖3a顯示了各區(qū)域SO2隨時(shí)間的變化趨勢(shì)擬合圖,可以看出,除新疆和青藏高原地區(qū)SO2呈現(xiàn)略微增長(zhǎng)的趨勢(shì)外,其他地區(qū)呈現(xiàn)顯著降低的趨勢(shì).其中北部沿海地區(qū)的降幅最大,為-0.04DU/a,這種現(xiàn)象的產(chǎn)生與國(guó)家對(duì)該地區(qū)大氣污染治理的重視程度有關(guān),相應(yīng)地區(qū)大部分制定了較為嚴(yán)格的污染排放控制政策,如《京津冀及周邊地區(qū)大氣污染治理行動(dòng)計(jì)劃》等[31];對(duì)西北地區(qū)來(lái)說(shuō),一些省份為了加快能源和地方經(jīng)濟(jì)發(fā)展,免除了SO2排放控制的規(guī)定[32].從而東部和南部SO2的降低速率大于西北部.

由SO2傾向率可以看出(圖3b),在新疆烏魯木齊附近和重慶西南等地SO2濃度降幅較為明顯,烏魯木齊實(shí)行“煤改氣”的政策,遏制了該地區(qū)大氣污染惡化的狀況[33],降幅最大達(dá)到-0.09DU/a;重慶等地SO2濃度降低主要受工業(yè)布局的影響,工業(yè)企業(yè)搬遷和政策的實(shí)施,使該地區(qū)SO2降幅較為明顯,降幅最大達(dá)到-0.06DU/a.總的來(lái)說(shuō),近16年二氧化硫污染狀況得到了明顯的改善,主要與我國(guó)推出的一系列節(jié)能減排政策有關(guān).

2.2.2 空間分布特征 多年SO2平均濃度呈現(xiàn)東高西低的特征(圖4a),尤其在北部沿海(0.76DU)和黃河中游(0.57DU)地區(qū)濃度達(dá)到了最大,超過(guò)全國(guó)均值(0.34DU).由于該地區(qū)是我國(guó)主要的工業(yè)生產(chǎn)基地和燃煤產(chǎn)區(qū),尤其冬季采暖期以燃煤為主,且降水較少,導(dǎo)致SO2濃度較高[34].新疆烏魯木齊SO2濃度也較高,該地區(qū)的能源消耗結(jié)構(gòu)主要以燃煤為主,同時(shí)地形多為盆地,一定程度上減弱了SO2的擴(kuò)散[13].整體來(lái)說(shuō),我國(guó)SO2的空間分布與地區(qū)能源消耗和工業(yè)活動(dòng)有關(guān).

圖4 中國(guó)SO2多年均值(a)和變異系數(shù)圖(b)

圖4b可以看出,北部沿海和黃河中游地區(qū)的變異系數(shù)較高,說(shuō)明該地區(qū)SO2濃度年際波動(dòng)大,分別為0.266和0.208,高于全國(guó)尺度的變異系數(shù)(0.077),該地區(qū)受人為干擾的影響較大,因此變異系數(shù)較大,年際波動(dòng)大.在青藏高原和西北地區(qū)SO2的變異系數(shù)較低,由于該地區(qū)自然景觀較為完整,生態(tài)環(huán)境良好,工業(yè)活動(dòng)強(qiáng)度小,SO2在這些地區(qū)受人為因素和自然因素影響的程度相對(duì)較低,具有較低的年際波動(dòng).總體而言,除西北某些小部分變異系數(shù)較低,其他地區(qū)變異系數(shù)較大,且東南部的變異系數(shù)高于西北部,說(shuō)明東南部年際變化高于西北部.

2.3 SO2的主要來(lái)源分析

大氣中SO2的主要來(lái)源包括含硫燃料的燃燒,有色金屬的冶煉,化工產(chǎn)業(yè)和交通出行等生產(chǎn)活動(dòng)[35].本文選取MEIC公布的人為污染源排放清單和社會(huì)經(jīng)濟(jì)數(shù)據(jù),探討了SO2的主要來(lái)源,圖5給出了SO2與各影響因子的相關(guān)性分布圖.

結(jié)果表明,SO2柱濃度與工業(yè)源、民用源、電力源、交通排放源和第二產(chǎn)業(yè)占GRP比重呈現(xiàn)正相關(guān).燃煤產(chǎn)生大量的SO2,我國(guó)是世界上最大的煤炭生產(chǎn)和消費(fèi)國(guó),排放的SO2約90%來(lái)自于燃煤.除此而外,隨著城市化進(jìn)程的加快,機(jī)動(dòng)車(chē)生產(chǎn)和使用量急劇增長(zhǎng),排放的氣體對(duì)大氣環(huán)境污染日趨嚴(yán)重,2022年《中國(guó)移動(dòng)源環(huán)境管理年報(bào)》非道路移動(dòng)源二氧化硫排放達(dá)到16.8萬(wàn)t,因此工業(yè)污染、生活污染和機(jī)動(dòng)車(chē)產(chǎn)生的SO2是中國(guó)大氣中SO2的主要來(lái)源.其中在黃河中游、長(zhǎng)江中上游、東部沿海和北部沿海地區(qū)的相關(guān)性較強(qiáng),相關(guān)系數(shù)達(dá)到了0.8以上,這些地區(qū)經(jīng)濟(jì)發(fā)達(dá),人口較多,工業(yè)的發(fā)展主要依靠化石燃料和煤炭的燃燒,排放大量的SO2造成嚴(yán)重的空氣污染.

人口密度、工業(yè)廢氣治理設(shè)施和工業(yè)廢氣治理費(fèi)用和SO2呈現(xiàn)負(fù)相關(guān).工業(yè)廢氣治理設(shè)施和工業(yè)廢氣治理費(fèi)用的逐年增長(zhǎng),推動(dòng)了SO2空氣濃度的降低,尤其中部和東部地區(qū)明顯,負(fù)相關(guān)達(dá)到了0.8,說(shuō)明這些區(qū)域工業(yè)污染物處理技術(shù)先進(jìn),國(guó)家投資充足,污染治理成效顯著.

由圖5可知,工業(yè)排放源、民用源、電力源、機(jī)動(dòng)車(chē)排放源和第二產(chǎn)業(yè)占GRP比重等因素影響著空氣中的SO2濃度.為進(jìn)一步探討各研究區(qū)SO2變化的主要驅(qū)動(dòng)因子,我們借助隨機(jī)森林模型評(píng)估各變量變化引起SO2變化的重要程度(圖6).

在黃河上游、黃河中游、北部沿海和新疆地區(qū),工業(yè)源的變化是影響SO2變化的主導(dǎo)因子.我國(guó)主要的煤炭開(kāi)采和煤炭深加工基地位于這些地區(qū),高硫煤的儲(chǔ)量占煤炭總儲(chǔ)量的20%以上,因此在煤炭的利用過(guò)程中,會(huì)排放大量的SO2.在東北、東部沿海和青藏高原地區(qū),由于以煤炭為主的重工業(yè)水平較低,而居民爐灶和商業(yè)飲食業(yè)燃煤排放的SO2對(duì)城市大氣環(huán)境質(zhì)量影響占據(jù)主要位置,因此居民源的變化直接影響空氣中SO2濃度的變化.交通源的變化是東南沿海和長(zhǎng)江中上游地區(qū)SO2變化的主要驅(qū)動(dòng)因子,該地區(qū)經(jīng)濟(jì)的快速發(fā)展吸引了大量人口,因此對(duì)交通的需求量大,汽車(chē)的保有量較大,汽車(chē)SO2的排放量大.在珠江中上游地區(qū),電力源的變化對(duì)SO2變化影響較大,該地區(qū)耗電量較大,電力行業(yè)的主體是火力發(fā)電,主要?jiǎng)恿κ敲禾?貴州作為我國(guó)南方最大的產(chǎn)煤省區(qū),2019年煤電占全省總裝機(jī)容量的51%,廣西煤電占總裝機(jī)的39%,因此燃煤電廠是排放SO2等氣態(tài)酸物質(zhì)的主要污染源.綜上,排放源的變化是引起二氧化硫的變化的主導(dǎo)因子.

3 結(jié)論

3.1 二氧化硫柱濃度多年整體呈現(xiàn)下降趨勢(shì).在新疆(0.0008DU/a)和青藏高原(0.0004DU/a)地區(qū)SO2呈現(xiàn)略微增長(zhǎng)的趨勢(shì),在其他地區(qū)呈現(xiàn)顯著降低的趨勢(shì),其中北部沿海地區(qū)的降幅最大,為-0.04DU/a,超過(guò)全國(guó)平均傾向率(-0.005DU/a).

3.2 二氧化硫柱濃度在空間分布上呈現(xiàn)東高西低的特征,北部沿海和黃河中游地區(qū)濃度最大,分別為0.76DU和0.57DU,超過(guò)全國(guó)均值(0.34DU),新疆地區(qū)二氧化硫濃度最低,為0.24DU.

3.3 工業(yè)源、居民源、電力源和機(jī)動(dòng)車(chē)排放源是大氣中SO2的主要來(lái)源,工業(yè)廢氣治理設(shè)施和工業(yè)廢氣治理費(fèi)用的逐年增長(zhǎng),推動(dòng)SO2濃度的降低.

3.4 排放源的變化主要影響了二氧化硫的變化.其中黃河上游、黃河中游、北部沿海和新疆地區(qū),工業(yè)源的變化是影響二氧化硫變化的主導(dǎo)因子;東北、東部沿海和青藏高原地區(qū),居民源的變化主要影響了二氧化硫濃度的變化;交通源的變化是東南沿海和長(zhǎng)江中上游地區(qū)二氧化硫變化的主要驅(qū)動(dòng)因子;在珠江中上游地區(qū)電力源的變化對(duì)二氧化硫變化影響較大.

[1] Qian H Q, Xu S D, Cao J, et al. Air pollution reduction and climate co-benefits in China’s industries [J]. Nature Sustainability, 2021,4(5).

[2] 蔡博峰.中國(guó)城市二氧化碳排放空間特征及與二氧化硫協(xié)同治理分析 [J]. 中國(guó)能源, 2012,34(7):33-37,19. Cai B F.Spatial characteristics of urban carbon dioxide emissions in China and analysis of synergistic treatment with sulfur dioxide [J]. Energy of China, 2012,34(7):33-37,19.

[3] Greenpeace Environment Trust (GET), 2019. Global SO2emission hotspot database: Ranking the world’s worst sources of SO2pollution. Greenpeace Environment Trust. https://www.greenpeace.org.au/wp/ wp-content/uploads/2019/08/Global-Hotspot-and-Emission-Sources-for-SO2_August-2019_AU_final.pdf.

[4] 中國(guó)環(huán)境公報(bào) [R]. 2022. www.cnemc.cn/jcbg/zghjzkgb/202206/ P020220624327755600688.pdf.China Environmental Bulletin [R]. 2022. www.cnemc.cn/jcbg/ zghjzkgb/202206/P020220624327755600688.pdf.

[5] 燕 麗,賀晉瑜,楊曉玥,等.基于衛(wèi)星遙感數(shù)據(jù)和污染控制政策分析河南省2001～2018年SO2濃度變化特征 [J]. 環(huán)境科學(xué)研究, 2021, 34(3):531-537.Yan L, He J Y, Yang X Y, et al. Analysis of SO2concentration in Henan Province from 2001to 2018based on satellite remote sensing data and pollution control policy [J]. Research of Environmental Sciences, 2021,34(3):531-537.

[6] Karplus V J, Zhang S. Quantifying coal power plant responses to tighter SO2emissions standards in China [J]. Energy Weekly News, 2018,115(27),7004–7009.

[7] 張連華,馬鵬飛,周春艷,等.2005～2018年全國(guó)SO2柱濃度時(shí)空分布特征分析 [J]. 遙感信息, 2020,35(5):8. Zhang L H, Ma P F, Zhou C Y, et al. Spatial-temporal distribution characteristics of SO2column concentration in China from 2005 to 2018 [J]. Remote Sensing Information, 2020,35(5):8.

[8] Ranaarif S, Yuwono A S. Analysis of the distribution of sulfur dioxide (SO2) pollutant in Bali Island for the 2011～2020 period [J]. IOP Conference Series: Earth and Environmental Science, 2021,871(1).

[9] Yang X, Wang S J, Zhang W Z, et al. The impact of anthropogenic emissions and meteorological conditions on the spatial variation of ambient SO2concentrations: A panel study of 113 Chinese cities [J]. Science of the Total Environment, 2017:584-585.

[10] Wang T, Wang P, Theys N, et al. Spatial and temporal changes in SO2regimes over China in the recent decade and the driving mechanism, [M]. Atmos. Chem. Phys., 18,18063–18078.

[11] Zhai H R, Yao J Q, Qi J W, et al. Spatial-temporal characteristics and influencing factors of NO2and SO2in the Beijing-Tianjin-Hebei region based on satellite remote sensing [J]. IOP Conference Series: Earth and Environmental Science, 2021,687(1).

[12] 劉濤濤,王勇輝,劉西剛.山西省SO2時(shí)空變化特征及影響要素分析 [J]. 環(huán)境工程, 2019,37(12):153-160. Liu T T, Wang Y H, Liu X G. Analysis of spatiotemporal variation characteristics and influencing factors of SO2in Shanxi Province [J]. Environmental Engineering, 2019,37(12):153-160.

[13] 魏 毅,萬(wàn)旭榮,李朝陽(yáng),等.烏魯木齊市二氧化硫污染特征與防治對(duì)策 [J]. 四川環(huán)境, 2011,30(6):76-80.Wei Y, Wan X R, Li C Y, et al. Characteristics and prevention and control countermeasures of sulfur dioxide pollution in Urumqi [J]. Sichuan Environment, 2011,30(6):76-80.

[14] Liu P F, Song H Q, Wang T H , et al. Effects of meteorological conditions and anthropogenic precursors on ground-level ozone concentrations in Chinese cities [J]. Environmental Pollution, 2020, 262(prepublish).

[15] Li X Y, Song H Q, Zhai S Y, et al. Particulate matter pollution in Chinese cities: Areal-temporal variations and their relationships with meteorological conditions (2015～2017) [J]. Environmental Pollution, 2018,246.

[16] Zhang X Y, Wang Z, Cheng M M, et al. Long-term ambient SO2concentration and its exposure risk across China inferred from OMI observations from 2005 to 2018 [J]. Atmospheric Research, 2021,247.

[17] Jiang L, He S X, Zhou H F, et al. Coordination between sulfur dioxide pollution control and rapid economic growth in China: Evidence from satellite observations and spatial economic models [J]. Structural Change and Economic Dynamics, 2021,57:279-291.

[18] Si Y D, Wang H M, Cai K, et al. Long-term (2006～2015) variations and relations of multiple atmospheric pollutants based on multi- remote sensing data over the North China Plain [J]. Environmental Pollution, 2019,255:311-323.

[19] Shen N, Zhao X, Li L, et al. Spatial and temporal variation characteristics of atmospheric NO2and SO2in the Beijing-Tianjin- Hebei region before and after the COVID-19 outbreak. Air Qual Atmos Health, 2021,14:1175–1188.

[20] Li M, Liu H, Geng G N, et al. Corrigendum to anthropogenic emission inventories in China: A review [J]. National Science Review, 2018, 5(4):603.

[21] Wu Y, Shi K, Chen Z, et al. Developing improved time-series DMSP- OLS-like data (1992～2019) in China by integrating DMSP-OLS and SNPP-VIIRS [J]. IEEE Transactions on Geoscience and Remote Sensing, 2021,10.1109/TGRS.2021.3135333.

[22] 徐新良.中國(guó)GDP空間分布公里網(wǎng)格數(shù)據(jù)集.資源環(huán)境科學(xué)數(shù)據(jù)注冊(cè)與出版系統(tǒng)(http://www.resdc.cn/DOI),2017.DOI:10.12078/ 2017121102). Xu X L. Kilometer grid dataset of spatial distribution of China's GDP. Resource and environmental science data registration and publication system.(http://www.resdc.cn/DOI),2017.DOI:10.12078/2017121102).

[23] Stevens F R, Gaughan A E, Linard C, et al. Disaggregating census data for population mapping using random forests with remotely-sensed and ancillary data. [J]. PloS one, 2015,10(2):e0107042.

[24] Huang R., Ju T, Dong H, et al. Analysis of atmospheric SO2in Sichuan-Chongqing region based on OMI data [J]. Environ. Monit. Assess., 2021,193,849.

[25] 柯碧欽,何 超,楊 璐,等.華北地區(qū)地表臭氧時(shí)空分布特征及驅(qū)動(dòng)因子 [J]. 中國(guó)環(huán)境科學(xué), 2022,42(4):1562-1574.Ke B Q, He C, Yang L, et al. Spatial-temporal distribution characteristics and driving factors of surface ozone in North China [J]. China Environmental Science, 2022,42(4):1562-1574.

[26] Tiao G C, Reinsel G C, Xu D M, et al. Effects of autocorrelation and temporal sampling schemes on estimates of trend and spatial correlation [J]. Journal of Geophysical Research: Atmospheres, 1990, 95(D12).

[27] Breiman L. Random forests [J]. Machine Learning, 2001,45:5–32.

[28] 夏曉圣,陳菁菁,王佳佳,等.基于隨機(jī)森林模型的中國(guó)PM2.5濃度影響因素分析 [J]. 環(huán)境科學(xué), 2020,41(5):2057-2065. Xia X S, Chen J J, Wang J J, et al. Analysis of influencing factors of PM2.5concentration in China based on random forest model [J]. Environmental Science, 2020,41(5):2057-2065.

[29] 戚佩霓.多源衛(wèi)星遙感分析中國(guó)區(qū)域污染氣體特征 [D]. 南京信息工程大學(xué), 2021. Qi P N. Multi-source satellite remote sensing analysis of regional polluting gas characteristics in China [D]. Nanjing University of Information Science and Technology, 2021.

[30] 衛(wèi) 瑋,王黎娟,靳澤輝,等.基于OMI數(shù)據(jù)汾渭平原大氣SO2時(shí)空分布特征分析 [J]. 生態(tài)環(huán)境學(xué)報(bào), 2018,27(12):2276-2283. Wei W, Wang L J, Jin Z H, et al. Spatial and temporal distribution characteristics of atmospheric SO2in Fenwei Plain based on OMI data [J]. Ecology and Environmental Sciences, 2018,27(12):2276-2283.

[31] Miao Z, Liu S C, Chen X D. Driving factors and spatio-temporal features underlying industrial SO2emissions in “2+26” in North China and extended cities [J]. Chinese Journal of Population, Resources and Environment, 2020,18(4).

[32] Ling Z, Huang T, Zhao Y, et al. OMI-measured increasing SO2emissions due to energy industry expansion and relocation in northwestern China, Atmos. Chem. Phys., 17,9115–9131.

[33] 冒 瑩.烏魯木齊市煤改氣工程對(duì)采暖期SO2和NO2濃度變化的影響分析 [J]. 干旱環(huán)境監(jiān)測(cè), 2014,28(3):109-112. Mao Y. Analysis of the influence of coal-to-gas conversion project on the change of SO2and NO2concentration during heating period in Urumqi [J]. Arid Environmental Monitoring, 2014,28(3):109-112.

[34] 肖鐘湧,趙伯維,陳雅文,等.2005～2016年中國(guó)大氣邊界層SO2的時(shí)空變化趨勢(shì) [J]. 中國(guó)環(huán)境科學(xué), 2018,38(10):3621-3627. Xiao Z Y, Zhao B W, Chen Y W, et al. Spatial-temporal variation trend of atmospheric boundary layer SO2in China in 2005～2016 [J]. China Environmental Science, 2018,38(10):3621-3627.

[35] 史 家.基于OMI數(shù)據(jù)的中國(guó)二氧化硫時(shí)空分布及影響因素研究 [D]. 蘭州:西北師范大學(xué), 2020. Shi J. Spatial and temporal distribution and influencing factors of sulfur dioxide in China based on OMI data [D]. Lanzhou: Northwest Normal University, 2020.970166572@qq.com.

Spatial and temporal distribution of sulfur dioxide and main emission sources in China.

WEI Ye-xiang, ZHANG Xiao-yu*, ZHANG Hong

(College of Environment and Resources, Shanxi University, Taiyuan 030006, China)., 2023,43(11)：5678～5686

s：With the rapid development of urbanization and economy, China has become the world's third largest emitter of sulfur dioxide (SO2), posing severe challenges to human health and the sustainable development of social and economic. Based on OMSO2e products, this paper analyzed the spatial and temporal changes of atmospheric SO2column concentrations in China from 2005 to 2020 and used MEIC emission data and socio-economic data to explore the main factors affecting SO2changes in China's 10 economic regions. The results show that: (1) Over the past 16 years, SO2has decreased in a fluctuation way, with a slight increase in Xinjiang and the Qinghai-Tibet Plateau, while other regions showed a significant decreasing trend. (2) SO2exhibited large inter-annual variation and displayed a spatial distribution pattern of high levels in the east and low levels in the west. The northern coastal areas (NC) and the middle reaches of the Yellow River (MYR) had the highest annual average concentrations and coefficients of variation. (3) SO2was positively correlated with the proportion of industrial sources, residential sources, power sources, transportation sources and secondary industries in GRP while it was negatively correlated with population density, industrial emission gas treatment facilities and associated expenses. (4) Changes in emission sources were the main factors influencing variations in SO2.This paper provides a theoretical basis and reference for formulating relevant emission reduction policies and ecological environmental protection in China.

sulfur dioxide column concentration (SO2)；spatial-temporal distribution；sources of emissions；Random forest；influencing factors

X51

A

1000-6923(2023)11-5678-09

魏夜香(1996-),女,山西忻州人,山西大學(xué)碩士研究生,主要從事大氣環(huán)境遙感相關(guān)研究.970166572@qq.com.

魏夜香,張霄羽,張 紅.中國(guó)二氧化硫的時(shí)空分布及主要排放來(lái)源研究 [J]. 中國(guó)環(huán)境科學(xué), 2023,43(11):5678-5686.

Wei Y X, Zhang X Y, Zhang H.Spatial and temporal distribution of sulfur dioxide and main emission sources in China [J]. China Environmental Science, 2023,43(11):5678-5686.

2023-04-25

國(guó)家自然科學(xué)基金-山西煤基低碳聯(lián)合基金(U1910207)

* 責(zé)任作者, 教授, zhang_xyhz@sxu.edu.cn