2005~2017年全球大氣邊界層SO2時空變化

康重陽,趙 軍,宋國富

2005~2017年全球大氣邊界層SO2時空變化

康重陽1,2,趙 軍1*,宋國富3

(1.西北師范大學(xué)地理與環(huán)境學(xué)院,甘肅 蘭州 730070；2.隴東學(xué)院歷史與地理學(xué)院,甘肅 慶陽 745000；3.西藏自治區(qū)生態(tài)環(huán)境遙感監(jiān)測中心,西藏 拉薩 850000)

使用2005~2017年遙感數(shù)據(jù)研究了全球大氣邊界層SO2時空分布特征及變化趨勢.結(jié)果表明:空間分布上SO2呈現(xiàn)空間異質(zhì)性,大氣高SO2柱量值集中在以火山噴發(fā)為代表的自然源區(qū)域和以工業(yè)排放為代表的人為源及附近區(qū)域;一、二、三級SO2柱量值在全球范圍內(nèi)整體呈現(xiàn)緯度地帶性分布特征,北半球受人為影響較為明顯,分界線在陸地區(qū)域向南凸出,海洋上向北凹陷,而南半球受人為影響較小,分界線呈現(xiàn)與緯線平行趨勢;2005~2017年全球大氣邊界層SO2單元柵格年均值整體呈現(xiàn)先增后減趨勢,火山噴發(fā)導(dǎo)致2008、2009、2011年夏季出現(xiàn)了明顯的波動,其余季節(jié)無顯著變化;全球范圍內(nèi)大氣邊界層SO2年內(nèi)變化,伴隨太陽直射點南北移動,2005~2014年一級SO2柱量值延緯向?qū)ΨQ軸的緯度,除10~12月份外,其余月份與對應(yīng)月15日太陽赤緯基本吻合.

全球大氣邊界層；OMI；OMPS；SO2

二氧化硫(SO2)是一種污染氣體,濕沉降危及生態(tài)系統(tǒng)、水生動物和植物健康、建筑的壽命[1],干沉降影響大氣環(huán)境、人體健康,及全球大氣環(huán)流、季風(fēng)、區(qū)域和局地氣候變化[2],甚至?xí)?dǎo)致平流層水汽含量增加[3].煤炭、石化燃燒為主的人為源,貢獻(xiàn)了全球接近70%的SO2排放量[1],其余來自于火山噴發(fā)、山火燃燒等為主的自然源.

20世紀(jì)發(fā)生在全球范圍內(nèi)的八大公害事件,其中4件與大氣中的SO2有直接關(guān)系,造成眾多生命的終結(jié)和巨大的經(jīng)濟(jì)損失.大氣中的SO2存在時間較短,Lee等[4]研究發(fā)現(xiàn),夏季白天SO2的壽命最短,約為19h,冬季白天的壽命約為56h.此外,大氣中 SO2還具有易氧化、易沉降的特點,導(dǎo)致其監(jiān)測比較困難.SO2的監(jiān)測經(jīng)過百年的發(fā)展[5-6],從早期的化學(xué)、電化學(xué)法逐漸發(fā)展為基于光譜技術(shù)的地基監(jiān)測、機(jī)載觀測、遙感監(jiān)測.化學(xué)、電化學(xué)監(jiān)測出現(xiàn)在20世紀(jì)40~70年代歐美發(fā)達(dá)國家,易受操作、采集方式等因素限制;地基監(jiān)測是基于“自下而上”的監(jiān)測方式,具有實時、高精度的優(yōu)點,但受制于操作復(fù)雜、維護(hù)成本高、測量范圍小等約束,難以實現(xiàn)大范圍SO2監(jiān)測;機(jī)載觀測機(jī)動性較高,但成本高昂,多應(yīng)用于數(shù)據(jù)驗證方面;1983年Krueger等[7]利用TOMS數(shù)據(jù),第一次從遙感衛(wèi)星數(shù)據(jù)提取火山噴發(fā)排放的SO2,大氣SO2監(jiān)測從地表轉(zhuǎn)向太空,實現(xiàn)“自上而下”的監(jiān)測,尤其是在紫外高光譜技術(shù)得到廣泛應(yīng)用之后,SO2遙感監(jiān)測技術(shù)取得實質(zhì)性的進(jìn)展.遙感監(jiān)測大氣SO2的傳感器主要有TOMS(Total Ozone Mapping Spectrometer)[7]、GOME (Global Ozone Monitoring Experiment)[8-9]、SCIAMACHY(The Scanning Imaging Absorption Spectrometer for Atmospheric Chartography )[10-11]、OMI(Ozone Monitor Instrument)[12-13]和OMPS(Ozone Mapping and Profiler Suite)[14-15].

近年來,國內(nèi)外學(xué)者應(yīng)用GOME、SCIAMACHY、GOME-2、OMI和OMPS遙感數(shù)據(jù),在大氣SO2監(jiān)測方面展開大量研究,研究主要集中在反演算法的開發(fā)及其數(shù)據(jù)準(zhǔn)確性的驗證[16-21]、區(qū)域SO2監(jiān)測及排放量隨時間變化趨勢[22-35]、SO2排放源頭目錄生成[15,36]、火山噴發(fā)SO2排放量估算[37-41]等方面, OMI和OMPS憑借其較高的空間分辨率、時間分辨率和光譜分辨率的優(yōu)勢,應(yīng)用更為廣泛. NASA科研工作者就OMI和OMPS監(jiān)測SO2準(zhǔn)確性方面進(jìn)行了相關(guān)研究,表明二者在重點監(jiān)測區(qū)域的相關(guān)系數(shù)高達(dá)0.8以上[15],這有利于使用兩種不同傳感器進(jìn)行長時間序列變化趨勢方的研究.先前的研究主要集中在區(qū)域尺度,在全球時空變化方面研究較少.本研究基于2005~2014年OMI SO2L2V003、2015~2017年OMPS NMSO2L2V1.1遙感數(shù)據(jù)中的PBL(The planetary boundary layer) SO2數(shù)據(jù),進(jìn)行2005~2017年全球大氣邊界層 SO2空間分布特征,年際、年內(nèi)變化趨勢研究,并定量的分析了全球PBL SO2年內(nèi)變化與太陽赤緯之間的對應(yīng)關(guān)系,有助于在全球氣候變化背景下,進(jìn)一步認(rèn)識的大氣邊界層SO2的空間分布與變化.

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

1.1 數(shù)據(jù)來源

臭氧監(jiān)測儀(OMI),是繼GOME、GOME-2和SCIAMACHY之后的新一代高光譜大氣成分探測傳感器,基于近紫外/可見光波段近地觀測的CCD光譜儀,搭載于NASA地球觀測系統(tǒng)(EOS)Aura衛(wèi)星之上,該衛(wèi)星由荷蘭、和芬蘭聯(lián)合研發(fā),是一款太陽同步衛(wèi)星,赤道過境時間為13:45,于2004年7月15日發(fā)射升空,同年8月9日開始采集數(shù)據(jù),星下點空間分辨率為13×24km2,每日產(chǎn)生14~15景影像,能夠?qū)崿F(xiàn)日全球覆蓋,每景影像為38MB.其主要目的是研究大氣痕量成分及其對氣候變化的影響,OMI已經(jīng)被廣泛應(yīng)用于大中尺度的污染氣體監(jiān)測等領(lǐng)域[14,21,24-28].

臭氧成像探測儀(OMPS),作為OMI傳感器的繼承者,搭載于美國新一代對地觀測衛(wèi)星Suomi NPP之上,2011年8月28日發(fā)射,2012年1月開始采集數(shù)據(jù)[14-15],運(yùn)行高度為839km,赤道過境時間為13:30,比前代傳感器OMI早10~15min,全球范圍內(nèi)每天產(chǎn)生14~15景影像,每景影像約為1.14MB.主要目的是用于監(jiān)測全球范圍內(nèi)平流層O3分布及變化,延續(xù)了1970年以來對全球平流層O3時空分布與變化特征長期監(jiān)測,同時也監(jiān)測如SO2、NO2、氣溶膠光學(xué)厚度、云量等地球物理量.

以上遙感數(shù)據(jù)都來自于NASA GES DISC數(shù)據(jù)庫,該數(shù)據(jù)庫收集了關(guān)于大氣成分、水和能量循環(huán)、氣候變化相關(guān)的遙感數(shù)據(jù)資料.

1.2 數(shù)據(jù)及反演算法

PBL SO2數(shù)據(jù)反映的是剖面高度小于2km大氣中的SO2總柱量值[15].

SO2反演算法:OMI SO2反演算法經(jīng)歷了DOAS[16](Differential Optical Absorption Spectroscopy)、LF[17](Linear Fit)、BRD[18](Band Residual Difference),到新一代通用的PCA[19](Principal components analysis)算法,監(jiān)測的準(zhǔn)確性不斷得到提升,OMPS NMSO2也延用了PCA算法.PCA[19]算法由Li等2013年提出,該算法使用赤道太平洋無SO2排放區(qū)域的輻射值,主成分能夠捕獲接近99.9999%輻射量的變化.前3個主成分為臭氧吸收、表面反射、旋轉(zhuǎn)拉曼散射效應(yīng),其余主成分來源于測量噪聲和探測器偽成像.PCA算法的瞬時視場角的標(biāo)準(zhǔn)偏差,在赤道太平洋30°N~30°S,無SO2排放區(qū)域為0.5DU,是前一代反演算法BRD標(biāo)準(zhǔn)偏差的一半,在高緯度區(qū)隨著O3傾斜柱量值的增加,標(biāo)準(zhǔn)偏差上升0.7~0.9DU,這依然比BRD算法小一倍.其次使用雅克比行列式代替了先前算法中的空氣質(zhì)量系數(shù)(AMF),相同條件下提高了數(shù)據(jù)的準(zhǔn)確性.本文所使用的遙感數(shù)據(jù)均由PCA算法反演產(chǎn)生,消除了因反演算法不同導(dǎo)致長時間序列研究的差異問題.

1.3 數(shù)據(jù)處理方法

參考OMI、OMPS官方文檔,使用VISAN3.7及python編寫程序,完成OMSO2、NMSO2數(shù)據(jù)提取,處理過程如下:

(1)逐日數(shù)據(jù)剔除及提取:剔除云參數(shù)(Cloud Fraction)>0.3、太陽高度角>70°、SO2柱量值<0DU的數(shù)據(jù),提取全球范圍內(nèi)經(jīng)、緯度及對應(yīng)PBL SO2柱量值.

(2)將產(chǎn)生的數(shù)據(jù)轉(zhuǎn)成點矢量文件,以0.125°網(wǎng)格做普通克里金插值,并進(jìn)行投影轉(zhuǎn)換.

(3)在時間尺度上,獲取大氣邊界層SO2年均柱量值、季均柱量值和月均柱量值,并在空間上完成SO2年均空間分布圖、季節(jié)變化、月均變化趨勢圖.

2 結(jié)果與討論

以往研究以等間距作為標(biāo)準(zhǔn)進(jìn)行劃分,有0.2DU、0.5DU.通過前期研究發(fā)現(xiàn)[37],自然斷點法更能夠準(zhǔn)確反映SO2的空間分布,本文對全球SO2柱量值的劃分如表1所示.

表1 全球SO2柱量值七級劃分標(biāo)準(zhǔn)

自然斷點法是一種基于自然分組,將屬性值相近的值進(jìn)行合并,是類內(nèi)差異最小,類間差異最大,方差和最小的一種統(tǒng)計學(xué)分類方法.該方法可以將分類間隔加以識別,對相似值進(jìn)行七級分組,并能使各個類之間的差異最大化.要素被劃分為多個類,對于這些類,會在數(shù)據(jù)值得差異相對較大的位置設(shè)置分類邊界.

2.1 全球大氣邊界層SO2空間分布

全球范圍內(nèi)大氣邊界層SO2在空間上呈現(xiàn)異質(zhì)性,六、七級高柱量值呈現(xiàn)聚集狀分布,一、二、三級柱量值區(qū)域呈現(xiàn)緯度地帶性,南、北半球表現(xiàn)出不均勻?qū)ΨQ性.

圖1 2005~2017年全球范圍內(nèi)大氣邊界層SO2平均柱量值空間分布

圖中南美洲異常值區(qū)域(South Atlantic Anomaly),以下統(tǒng)一簡稱為SAA

火山噴發(fā)向大氣中排放大量的SO2氣體,大部分火山活動強(qiáng)烈,對SO2濃度的影響高至平流層,而少部分火山比較溫和噴發(fā)高度低,主要引起大氣邊界層SO2濃度增加.如圖1,西太平洋菲律賓海東北馬里納亞群島的安納塔漢(Anatahan)火山、澳大利亞東北部瓦努阿圖共和國的瓦努阿圖群島(Vanuatu)火山、非洲剛果民主共和國東部的尼拉貢戈(Nyiragongo)火山、美國夏威夷的冒納凱阿(Mauna Kea)火山、厄瓜多爾的阿蘇爾(Cerro Azul)火山、巴布亞新幾內(nèi)亞獨立國(Papua,New Guinea)境內(nèi)火山等,都引起大氣邊界層SO2濃度的增加(詳細(xì)區(qū)域請參照NASA全球SO2監(jiān)測首頁:https://so2.gsfc. nasa.gov/).石油、煤炭等化學(xué)能源的大量使用,向大氣集中排放大量SO2氣體,形成人為高SO2排放地區(qū).如:俄羅斯北部克拉斯諾亞爾斯克邊疆區(qū)的諾里爾斯克,是俄羅斯主要的有色金屬工業(yè)基地,周邊存在大型火電站,煤炭消耗量巨大;墨西哥合眾國的中部地區(qū),以其首都墨西哥城為中心,形成的六級SO2柱量值區(qū)域,受地理條件、地形特點的影響,燃料難以充分燃燒,SO2擴(kuò)散困難,產(chǎn)生堆積現(xiàn)象[42];中國高柱量值區(qū)域涵蓋北京、天津、河北、山東全境,遼寧大部、山西中部和南部、河南中部和北部以及安徽省和江蘇省北部,發(fā)達(dá)的鋼鐵、石化、火電傳統(tǒng)高耗能產(chǎn)業(yè),排放大量的SO2氣體;非洲南非共和國的東北部,電廠、工業(yè)排放、擴(kuò)散,造成該地區(qū)出現(xiàn)人為高SO2柱量值[43].

45°N~45°S區(qū)域內(nèi),空間上整體呈現(xiàn)緯度地帶性分布,熱帶地區(qū)最低,伴隨緯度增加,SO2柱量值增加,表現(xiàn)出不均勻?qū)ΨQ性.一級柱量值區(qū)域SO2年均柱量值小于0.512DU,包含了熱帶(23.5°N~ 23.5°S)、部分北溫帶地區(qū)(23.5°N~30°N)和部分南溫帶地區(qū)(23.5°S~30°S),是全球大氣邊界層SO2柱量值最低的區(qū)域;二級柱量值區(qū)域包含了北溫帶30°N~45°N和南溫帶30°S~45°S區(qū)域,SO2柱量值處于0.512~ 0.615DU,南半球一、二級SO2柱量值分界線與緯線平行,平行效果優(yōu)于北半球.

圖中SAA區(qū)域的SO2明顯高于周邊地區(qū),是由于地球的磁場捕獲高能帶電粒子,該粒子會對低地球軌道衛(wèi)星探測器造成影響,使得該區(qū)域的輻照度高于正常值,產(chǎn)生了高SO2柱量值現(xiàn)象[15],所以本文在討論全球SO2空間及時空變化時不考慮該區(qū)域的變化,如何剔除SAA區(qū)域的地磁影響,有待進(jìn)一步工作的開展.

2.2 2005~2017年全球SO2年際變化特征

OMI傳感器從2007年開始出現(xiàn)數(shù)據(jù)壞道,并逐年增加,雖然NASA工作人員嘗試修復(fù)壞道數(shù)據(jù),但大范圍上較難準(zhǔn)確的反映大氣邊界層SO2柱量值,所以本研究選用2005年1月1日至2014年12月31日OMISO2遙感數(shù)據(jù),及2015年1月1日~2017年12月31日OMPS遙感數(shù)據(jù),以進(jìn)行近13年全球大氣邊界層SO2年際變化研究.

由于OMI、OMPS傳感器自身的局限性,采集數(shù)據(jù)受制于太陽照射,太陽直射點及地球在公轉(zhuǎn)軌道的位置不同,采集數(shù)據(jù)的區(qū)域有所差別,再加上火山噴發(fā)排放SO2對大氣邊界層的影響具有不確定性,所以在進(jìn)行年際變化比較時,本文采用單元柵格年均值,即該區(qū)域內(nèi)年均SO2總柱量值除以該區(qū)域的總柵格數(shù).如圖2所示,2005~2017年全球SO2單元柵格年均值整體呈現(xiàn)先增后減趨勢,最大值出現(xiàn)在2011年,單元柵格年均值為0.613DU,最小值在2016年,單元柵格年均值為0.302DU.2005~2007年單元柵格年均值基本保持不變,約為0.55DU,從2008年開始緩慢增加,到2011年出現(xiàn)最大值,隨后2012~ 2013年有所下降,2014年出現(xiàn)第二峰值,約為0.588DU, 2014~2015年出現(xiàn)較大幅度降低,較2014年降低幅度高達(dá)45%,這種較大幅度的變化很大一部分可能是由于傳感器的不同造成,如何消除不同傳感器之間的差異,進(jìn)行長時間序列的時空變化特征研究,有待新的估算算法的開發(fā).

圖2 2005~2017年全球范圍內(nèi)大氣邊界層SO2單元柵格年平均柱量值變化趨勢

圖中2014~2015年之間的巨大降幅主要源于傳感器的差異,圖4同

2.3 全球大氣邊界層SO2季節(jié)變化

為研究不同季節(jié)大氣邊界層SO2時空分布特征,分別取12月、次年1、2月為冬季,3、4、5月為春季,6、7、8月為夏季,9、10、11月為秋季,六、七級柱量值區(qū)域主要集中在火山噴發(fā)影響區(qū)域,如非洲中部的尼拉貢戈火山、美國的夏威夷火山等,和嚴(yán)重的人為排放區(qū)域,如南非共和國、俄羅斯的諾里爾斯克、中國的華北平原、墨西哥合眾國的墨西哥城等.一、二、三級區(qū)域分界線,在北半球人為排放較為嚴(yán)重的區(qū)域向南凸出,海洋上向北凹陷,南半球基本不受人為排放影響,表現(xiàn)出與緯線平行跡象.北半球一、二、三級柱量值分界線隨季節(jié)發(fā)生變化,夏季分界線向北達(dá)到最大,隨后隨太陽直射點的南移向南變化,與冬季一、二級分界線到達(dá)北回歸線附近,到達(dá)年內(nèi)最南端,南半球相反.

2.4 2005~2017年SO2季均值年際變化趨勢

通過圖2發(fā)現(xiàn),2005~2017年全球SO2單元柵格年均值整體呈現(xiàn)先增加后降低過程,為了解2005~2017年不同季節(jié)的變化趨勢,本文通過計算2005~2017 年不同季節(jié)單元柵格季節(jié)均值,來估算不同季節(jié)SO2的年際變化狀況.如圖4,2005~2017年SO2單元柵格季均值均出現(xiàn)不同程度的變化,變化幅度最大的季節(jié)是夏季,其余季節(jié)SO2單元柵格季均值在2005~2017年間無明顯變化,夏季在2008、2009、2011年出現(xiàn)突變,主要受火山噴發(fā)排放大量SO2在大氣邊界層隨高空氣流擴(kuò)散,出現(xiàn)大面積SO2柱量值增加.

圖3 2005~2017年全球范圍內(nèi)大氣邊界層SO2平均柱量值季節(jié)變化

Fig.3 Global seasonal variation of averaged PBL SO2 columns between 2005 and 2017

2.5 2005~2017年全球SO2年內(nèi)變化特征

太陽直射點的南北移動,引起全球范圍內(nèi)眾多氣象要素的周期性變化,如氣壓、溫度、降水等,而這些氣象要素的周期性變化直接影響到大氣邊界層SO2的干沉降、濕沉降、光化學(xué)反應(yīng)等,大氣邊界層SO2柱量值也表現(xiàn)出周期性變化的特征.

全球范圍內(nèi)大氣邊界層SO2年內(nèi)變化,伴隨太陽直射點南北移動,一級SO2柱量值區(qū)域周期性變化.1、2、12月份,一級SO2柱量值區(qū)域包含了整個熱帶及大部分的南溫帶,這三個月南半球處于夏季,高溫多雨不利于SO2在大氣中單獨存在.伴隨太陽直射點由南回歸線向赤道轉(zhuǎn)移,一級SO2柱量值區(qū)域開始向北半球轉(zhuǎn)移,北半球氣溫回暖,采暖結(jié)束,二、三級SO2柱量值區(qū)域緯度范圍增大,至5月份,除俄羅斯大部分區(qū)域外,二級SO2柱量值區(qū)域包含了整個北溫帶及北寒帶的部分地區(qū),南半球氣溫降低,三、四級柱量值區(qū)域面積明顯擴(kuò)大;6、7、8月,北半球處于夏季,SO2的干、濕沉降速度快,大氣逆溫層強(qiáng)度降低、高度升高,SO2擴(kuò)散增強(qiáng),正常情況下,三級及以上SO2柱量值區(qū)域面積急劇減少,2009年6、7月份美國阿拉斯加火山、俄羅斯東部的勘察加火山噴發(fā),火山噴發(fā)排放的大量SO2,在大氣環(huán)流和行星風(fēng)系的影響下,一部分呈緯度地帶性向東延伸至加拿大境內(nèi),最遠(yuǎn)延伸至格陵蘭島,另一部分在西風(fēng)帶內(nèi)受極地西風(fēng)擴(kuò)散,在北半球中緯、中高緯、高緯度地區(qū)出現(xiàn)大面積的五、六、七級SO2柱量值區(qū)域;9月份開始,太陽直射點轉(zhuǎn)入南半球,北半球進(jìn)入秋季,氣溫降低,降水減少,北半球SO2柱量值升高,而南半球進(jìn)入春季,大氣中SO2的生存時間逐漸縮短,南半球大部分處于三級及以下SO2柱量值區(qū)域,并且伴隨南極進(jìn)入極晝,南極大氣邊界層出現(xiàn)SO2空洞.

2.6 全球大氣邊界層SO2年內(nèi)變化與太陽赤緯關(guān)系

太陽輻射作為地球最主要的能量來源和基本動力,推動地球表層幾乎全部的自然物理過程,使地理環(huán)境得以形成和有序發(fā)展[44].通過上文發(fā)現(xiàn)地球大氣邊界層SO2在空間分布上呈現(xiàn)緯度地帶性,通過圖5發(fā)現(xiàn)全球范圍內(nèi)大氣邊界層SO2年內(nèi)變化表現(xiàn)出一定運(yùn)動的特征,為定量研究其特征,本文選取太平洋無SO2排放區(qū)域(經(jīng)度范圍150°W~170°W、緯度范圍70°S~70°N),測量一級SO2柱量值延緯向?qū)ΨQ軸的緯度,發(fā)現(xiàn)除10~12月份外,其余月份對稱軸緯度與對應(yīng)月15日太陽赤緯基本吻合,10~12月份,相差較大與南半球大面積二級柱量值區(qū)域有關(guān),導(dǎo)致一級柱量值區(qū)域被壓縮.如表3所示,基于OMI遙感數(shù)據(jù),2005年~2014年一級SO2柱量值區(qū)域緯向?qū)ΨQ軸與每月15日太陽赤緯之間存在較好的對應(yīng)關(guān)系,且存在相近的運(yùn)動趨勢.

太陽赤緯的計算:太陽直射點的變化決定著太陽高度的變化,直接影響到地區(qū)的太陽輻射量.本文采用黃玫生太陽直射點的計算方法[45],計算每月15日的太陽直射點緯度位置.

表2 2005~2014年一級柱量值緯向?qū)ΨQ軸與每月15日太陽赤緯對應(yīng)關(guān)系

月份一級柱量值對稱軸緯度位置 200920102011201220132014 118.19°S17.69°S17.9°S15.93°S17.67°S18.3°S 211.44°S13.43°S7.6°S9.56°S11.67°S13.56°S 30.63°S4.3°S5.31°S0°2.38°S4.56°S 44.375°N2.9°N3.75°N8.88°N3.75°N2.69°N 514.31°N12.5°N13.06°N16.19°N13.81°N12.9°N 65.75°N17.85°N18.81°N19°N17°N17.38°N 711.25°N15.4°N18.06°N16.69°N16.31°N16.38°N 816.06°N10.3°N9.37°N13.6°N10.68°N7.5°N 98.94°N6.3°N6.88°N6.88°N5°N1.88°N 102.25°N2.6°S14.88°S1.38°N0.2°N2°N 1115.38°S16.88°S18.31°S18.43°S11.5°S7.31°S 1220.19°S19.43°S20.31°S18.85°S20.79°S15.8°S

緯向?qū)ΨQ軸緯度=

A

-(

A

-

B

)/2 (2)

式中:為研究區(qū)域北半球一、二級柱量值分界線最大緯度,為研究區(qū)域南半球一、二級柱量值分界線最大緯度.

3 結(jié)論

3.1 全球范圍內(nèi)大氣邊界層SO2柱量值空間上呈現(xiàn)空間異質(zhì)性,六、七級高柱量值呈現(xiàn)聚集狀分布,主要集中在以火山噴發(fā)為代表的自然源和以工業(yè)排放為代表的人為源及附近區(qū)域; SO2柱量值在熱帶地區(qū)最低,伴隨緯度增加而增長,南、北半球表現(xiàn)出不均勻?qū)ΨQ性.

3.2 2005~2017年全球大氣邊界層SO2單元柵格年均值整體呈現(xiàn)先增后減趨勢,最大值出現(xiàn)在2011年,最小值出現(xiàn)在2016年,受傳感器影響,2014年之后出現(xiàn)較大幅度的下降,很大一部分由傳感器不同引起;除夏季外,單元柵格季均SO2柱量值在2005~2017年之間沒有出現(xiàn)較大幅度的變化,火山噴發(fā)導(dǎo)致2008、2009、2011年夏季出現(xiàn)了明顯的突變,火山噴發(fā)排放大量SO2在大氣邊界層隨高空氣流擴(kuò)散,出現(xiàn)大面積SO2柱量值增加.

3.3 全球范圍內(nèi)大氣邊界層SO2整體呈現(xiàn)緯度地帶性分布特征,北半球受人為影響較為明顯,一、二、三級分界線在陸地區(qū)域向南凸出,海洋上向北凹陷,而南半球受人為影響較小,分界線呈現(xiàn)與緯線平行趨勢.

3.4 全球范圍內(nèi)大氣邊界層SO2年內(nèi)變化,受太陽直射點南北移動的影響.2005~2014年月均一級柱量值緯向?qū)ΨQ軸緯度與每月15日太陽赤緯之間存在良好的對應(yīng)關(guān)系,太陽輻射是造成大氣邊界層SO2柱量值年內(nèi)變化的最主要因素.

[1] Smith S J, Aardenne J V, Klimont Z, et al. Anthropogenic sulfur dioxide emissions: 1850~2005 [J]. Atmospheric Chemistry & Physics Discussions, 2011,10(6):16111-16151.

[2] Xu Q. Abrupt change of the mid-summer climate in central East China by the influence of atmospheric pollution [J]. Atmospheric Environment, 2001,35(30):5029–5070.

[3] Notholt J, Luo B P, Fueglistaler S, et al. Influence of tropospheric SO2emissions on particle formation and the stratospheric humidity [J]. Geophysical Research Letters, 2005,32(7):99-119.

[4] Lee C, Martin R V, Van Donkelaar A, et al. SO2emissions and lifetimes: Estimates from inverse modeling using in situ and global, space-based (SCIAMACHY and OMI) observations [J]. Journal of Geophysical Research Atmospheres, 2011,116(D6):D06307.

[5] 胡利芬,郭燕芬,馮 皓.二氧化硫監(jiān)測方法進(jìn)展[J]. 裝備環(huán)境工程, 2015,12(2):64-69. Hu Li-fen, Guo Yan-fen, Feng Hao. Progress in monitoring methods of sulfur dioxide [J]. Equipment Environmental Engineering, 2015, 12(2):64-69.

[6] 王炳華,趙 明.美國環(huán)境監(jiān)測一百年歷史回顧及其借鑒[J]. 環(huán)境監(jiān)測管理與技術(shù), 2000,13(6):13-17.Wang Bing-hua, Zhao Ming. The Review of USA environmental monitoring in one century and it’s reference to us [J]. Environmental Monitoring Management and Technology, 2000, 13(6):13-17.

[7] Krueger A J. Sighting of El Chichón sulfur dioxide clouds with the nimbus 7 total ozone mapping spectrometer [J]. Science, 1983, 220(7607):1377-1389.

[8] Burrows J P, Weber M, Buchwitz M, et al. The global ozone monitoring experiment (GOME): Mission concept and first scientific results [J]. Journal of Atmospheric Sciences, 1997,56(717):2370- 2352.

[9] Eisinger M, Burrows J P. Tropospheric sulfur dioxide observed by the ERS-2GOME instrument [J]. Geophysical Research Letters, 1998, 25(22):7177–7180.

[10] Bovensmann H, Burrows J P, Buchwitz M, et al. SCIAMACHY: Mission objectives and measurement modes [J]. Journal of Atmospheric Sciences, 1997,56(2):125-150.

[11] Afe O T, Richter A, Sierk B, et al. BrO emission from volcanoes: A survey using GOME and SCIAMACHY measurements [J]. Geophysical Research Letters, 2007,31(27):357-370.

[12] Xia Y, Zhao Y, Nielsen C P. Benefits of China's efforts in gaseous pollutant control indicated by the bottom-up emissions and satellite observations 2000~2017 [J]. Atmospheric Environment, 2016,136:73- 53.

[13] Levelt P F, Hilsenrath E, Leppelmeier G W, et al. Science objectives of the ozone monitoring instrument [J]. IEEE Transactions on Geoscience & Remote Sensing, 2006,77(5):1199-1208.

[14] Yang K, Dickerson R R, Carn S A, et al. First observations of SO2from the satellite Suomi NPP OMPS: Widespread air pollution events over China [J]. Geophysical Research Letters, 2013,70(18):7957- 7962.

[15] Zhang Y, Li C, Krotkov N A, et al. Continuation of long-term global SO2pollution monitoring from OMI to OMPS [J]. Atmospheric Measurement Techniques, 2017,10(7):1-21.

[16] Wagner T, Beirle S, Deutschmann T, et al. Satellite monitoring of different vegetation types by differential optical absorption spectroscopy (DOAS) in the red spectral range [J]. Atmospheric Chemistry & Physics, 2007,7(1):69-79.

[17] Yang K, Krotkov N A, Krueger A J, et al. Retrieval of large volcanic SO2columns from the Aura Ozone Monitoring Instrument: Comparison and limitations [J]. Journal of Geophysical Research Atmospheres, 2007,112(D27):1-17.

[18] Krotkov N A, Carn S A, Krueger A J, et al. Band residual difference algorithm for retrieval of SO2from the aura ozone monitoring instrument (OMI) [J]. IEEE Transactions on Geoscience & Remote Sensing, 2006,77(5):1259-1266.

[19] Li C, Joiner J, Krotkov N A, et al. A fast and sensitive new satellite SO2retrieval algorithm based on principal component analysis: Application to the ozone monitoring instrument [J]. Geophysical Research Letters, 2013,70(23):6317-6318.

[20] 閆歡歡,李曉靜,張興贏,等.大氣SO2柱總量遙感反演算法比較分析及驗證 [J]. 物理學(xué)報, 2016,65(8):177-160. Yan Huan-Huan, Li Xiao-Jing, Zhang Xing-Ying, et al. Comparison and validation of band residual difference algorithm and principal component analysis algorithm for retrievals of atmospheric SO2columns from satellite observations [J]. Acta Physica Sinica, 2016, 65(8):177-160.

[21] Krotkov N A, Mcclure B, Dickerson R R, et al. Validation of SO2retrievals from the Ozone Monitoring Instrument over NE China [J]. Journal of Geophysical Research Atmospheres, 2008,113(D16):16-70.

[22] Carn S A, Krueger A J, Krotkov N A, et al. Sulfur dioxide emissions from Peruvian copper smelters detected by the ozone monitoring instrument [J]. Geophysical Research Letters, 2007,37(9):1093-1101.

[23] Calkins C, Ge C, Wang J, et al. Effects of meteorological conditions on sulfur dioxide air pollution in the North China plain during winters of 2006~2015 [J]. Atmospheric Environment, 2016,177(10):296-309.

[24] 趙 軍,張斌才,樊潔平,等.基于OMI數(shù)據(jù)的蘭州及附近地區(qū)大氣邊界層SO2量值變化初步分析[J]. 遙感技術(shù)與應(yīng)用, 2011,26(6): 808-813. Zhao Jun, Zhang Bincai, Fan jieping, et al. The analysis of SO2value changes in the atmospheric boundary layer over Lanzhou and surrounding areas based on the OMI product [J]. Remote Sensing Technology and Application, 2011,26(6):808-813.

[25] 沈艷潔.基于OMI衛(wèi)星遙感數(shù)據(jù)的能源金三角空氣質(zhì)量初探 [D]. 蘭州:蘭州大學(xué), 2016. Shen Yanjie. A primary study for air quality of the Energy Golden Triangle, Northwest China based on the data of satellite remote sensing [D]. Lanzhou: Lanzhou University, 2016.

[26] 張興贏, Jos van Geffen,張 鵬,等.衛(wèi)星遙感監(jiān)測中國區(qū)域?qū)α鲗佣趸虻臅r空分布及其變化趨勢研究[C]//中國氣象學(xué)會2008年年會大氣環(huán)境監(jiān)測、預(yù)報與污染物控制分會場論文集, 2008,s6:607.Zhang Xing-ying, Jos van Geffen, Zhang Peng, et al.Study of temporal and spatial distribution of tropospheric SO2in china and Its changing trend by remote sensing [C]//Chinese Meteorological Society. Proceedings of the Atmospheric Environment Monitoring, Forecasting and Pollutant Control Sub-committee, 2008 Annual Meeting, 2008, s6:607.

[27] 康重陽,趙 軍,宋國富,等.基于OMI數(shù)據(jù)中國大氣邊界層SO2空間格局[J]. 中國環(huán)境科學(xué), 2018,38(2):735-773. Kang Chong-yang, Zhao Jun, Song Guo-fu, et al. Study on the spatial pattern of SO2in China's atmospheric boundary layer based on OMI data [J]. China Environmental Science, 2018,38(2):735-773.

[28] 肖鐘湧,趙伯維,陳雅文,等.2005~2016年中國大氣邊界層SO2的時空變化趨勢[J]. 中國環(huán)境科學(xué), 2018,38(10):23-29. Xiao Zhong-yong, Zhao Bo-wei, Chen Ya-wen, et al.Trend of temporal and spatial variation of planetary boundary layer SO2over China from 2005 to 2016 [J]. China Environmental Science, 2018, 38(2):735-773.

[29] Khokhar M F, Beirle S, Platt U , et al. 7-year temporal trend of anthropogenic SO2emissions over China identified from GOME observations[C]//Cospar Scientific Assembly. 36th COSPAR Scientific Assembly, 2006.

[30] Dickerson R R, Doddridge B G, Hains J C, et al. The SO2plume over the eastern US: Remote and in situ observations [C]//Egs-agu-eug Joint Assembly. EGS - AGU - EUG Joint Assembly, 2003,5:07797.

[31] Zhang X, Geffen J V, Liao H, et al. Spatiotemporal variations of tropospheric SO2over China by SCIAMACHY observations during 2004–2009 [J]. Atmospheric Environment, 2012,60:238-246.

[32] Mclinden C A, Fioletov V, Boersma K F, et al. Air quality over the Canadian oil sands: A first assessment using satellite observations [J]. Geophysical Research Letters, 2012,39(4):L04804- L04812.

[33] Huanhuan Y, Liangfu C, Lin S , et al. SO2columns over China: Temporal and spatial variations using OMI and GOME-2observations [J]. IOP Conference Series: Earth and Environmental Science, 2014, 17:012027.

[34] Levelt P F , Oord G H , Hilsenrath E , et al. Science goals of EOS-Aura's Ozone Monitoring Instrument (OMI) [C]//Agu Fall Meeting. AGU Fall Meeting Abstracts, 2003,6:02526-02527.

[35] Livingston J M, Redemann J, Russell P B, et al. Comparison of aerosol optical depths from the Ozone Monitoring Instrument (OMI) on Aura with results from airborne sunphotometry, other space and ground measurements during MILAGRO/INTEX-B [J]. Atmospheric Chemistry & Physics, 2009,9(2):537-570.

[36] 陳 輝,王 橋,厲 青,等.大氣環(huán)境熱點網(wǎng)格遙感篩選方法研究[J]. 中國環(huán)境科學(xué), 2018,38(7):2461-2470.Chen Hui, Wang Qiao, Wang Zhong-ting et al. Research on the method for remote sensing screening of atmospheric environment hotspot grid [J]. China Environmental Science, 2018,38(7):2461- 2470.

[37] Hayer C S, Wadge G, Edmonds M, et al. Sensitivity of OMI SO2measurements to variable eruptive behaviour at Soufrière Hills Volcano, Montserrat [J]. Journal of Volcanology and Geothermal Research, 2016,312:1-10.

[38] Sheldrake T E, Aspinall W P, Odbert H M, et al. Understanding causality and uncertainty in volcanic observations: An example of forecasting eruptive activity on Soufrière Hills Volcano, Montserrat [J]. Journal of Volcanology and Geothermal Research, 2017,341:287-300.

[39] Campion R. New lava lake at Nyamuragira volcano revealed by combined ASTER and OMI SO2measurements [J]. Geophysical Research Letters, 2015,41(21):7485-7492.

[40] Carn S A, Krotkov N A, Krueger A J, et al. The first daily measurements of passive volcanic degassing from space [C]//Agu Fall Meeting. AGU Fall Meeting Abstracts, 2006,V43E-05.

[41] Khokhar M F, Frankenberg C, Van Roozendael M, et al. Satellite observations of atmospheric SO2from volcanic eruptions during the time-period of 1996~2002 [J]. Advances in Space Research, 2005, 36(5):879-887.

[42] 方 鳴.墨西哥市大氣污染監(jiān)測與治理[J]. 全球科技經(jīng)濟(jì)瞭望, 1998,(9):49-50. Fang Ming. Monitoring and Governance of Air Pollution in Mexico City [J]. Global Science and Technology Economics Outlook, 1998, (9):49-50.

[43] Josipovic M, Annegarn H J, Kneen M A, et al. Concentrations, distributions and critical level exceedance assessment of SO2, NO2and O3in South Africa [J]. Environmental Monitoring & Assessment, 2010,171(14):181-196.

[44] 伍光和.自然地理學(xué)(第三版) [M]. 北京:高等教育出版社, 2000.Wu Guanghe. Physical Geography [M]. 3rd ed. Beijing: Higher Education Press, 2000.

[45] 黃玫生.利用公式求算太陽直射點和極晝極夜的緯度位置[J]. 中學(xué)地理教學(xué)參考, 1988,(5):11-12.Huang Meisheng. Using the formula to calculate the latitude position of the direct sunlight and the extreme polar night [J]. Middle School Geography Teaching Reference, 1988,(5):11-12.

Temporal and spatial distribution characteristics of global atmospheric boundary layer SO2based on remote sensing data from 2005 to 2017.

KANG Chong-yang1,2, ZHAO Jun1*, SONG Guo-fu3

(1.College of Geography and Environmental Science, Northwest Normal University, Lanzhou 730070, China；2.History and Geography Department, Long Dong College, Qingyang, 745000, China；3.Tibet Autonomous Region Remote Sensing Monitoring Center, Lhasa 850000, China)., 2019,39(10)：4033~4042

As the main precursor of the smog and acid rain in the atmosphere, SO2plays an important role to the near-surface O3pollution. Meanwhile, SO2also affects the global and local atmospheric circulations because its strong absorption on the ultraviolet influences the solar radiation dramatically. Therefore, it is beneficial to study the temporal and spatial distributions of SO2and their variations. In this paper, the spatial and temporal distributing features of SO2at global atmospheric boundary layer were studied by the remote sensing data from 2005 to 2017. A spatial heterogeneity was observed from the analyzing results. Specifically, the high columns of SO2mainly concentrated on the volcanic eruption areas as a natural source and on the industrial emission areas as an anthropogenic source. The first, second and third grades of SO2columns distributed zonally in the whole world. It was evidently affected by the human emissions in the Northern Hemisphere that the boundary line of SO2column grades protruded southward on the land and sunk northward on the ocean. In the Southern Hemisphere, on the contrary, the boundary line paralleled to the latitude. From 2005 to 2017, the annual average of SO2unit grid values at the global atmospheric boundary layer was increasing before 2011 and then decreasing. During this period of time, the volcanic eruptions led to obvious variations in the summer of 2008, 2009 and 2011. The alteration of SO2at the atmospheric boundary layer was associated with the movement of the direct solar radiation point within each year. From 2005 to 2014, the latitudinal axis of symmetry of the first grade SO2columns matched with the 15thday’s solar declination of each month, except October, November, and December.

global atmospheric boundary layer；OMI；OMPS；SO2

X511

A

1000-6923(2019)10-4033-10

康重陽(1983-),男,甘肅寧縣人,碩士,主要從事環(huán)境遙感與GIS應(yīng)用研究.發(fā)表論文2篇.

2019-03-10

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