珠江三角洲對流層HCHO柱濃度遙感監測及影響因子

咸 龍,葛建團,2*,徐 敏,陳雪萍,王 爽,謝順濤,胡文文

?

珠江三角洲對流層HCHO柱濃度遙感監測及影響因子

咸 龍1,葛建團1,2*,徐 敏3,陳雪萍1,王 爽1,謝順濤1,胡文文4

(1.西北師范大學地理與環境科學學院,甘肅 蘭州 730070；2.嘉峪關市政府,甘肅 嘉峪關 735100；3.蘭州交通大學環境與市政工程學院,甘肅 蘭州 730070；4.合水縣環境保護局,甘肅 慶陽 745400)

基于OMI衛星遙感反演數據,對珠江三角洲地區2009年~2016年對流層甲醛柱濃度時空分布特征及其影響因素進行研究.結果表明,珠江三角洲甲醛柱濃度時間變化特征為:8年來呈波動變化趨勢,年均值為13.11×1015molec/cm2,最低值出現于2012年,最高值出現于2016年;最大降低率為5.8%,最大增長率為6.3%.每年夏季最高,冬季最低,大小依次為夏季>秋季>春季>冬季,8a來96個月甲醛月際變化幅度較大,呈單峰結構,其中每年6月最高;空間變化特征為:甲醛柱濃度值由西北往東南遞減,其中以肇慶東北大部、佛山北部和廣州西部組成高值區分布中心,以佛山中南部、廣州東南半部和江門西北半部組成三級次高級分布區,以惠州、東莞、深圳、中山、珠海和江門等珠江三角洲近海岸地區為一二級低值濃度區;影響因素中氣溫與氣壓等氣象因素對HCHO的生成和分布有著促進作用,植被對HCHO的產生有一定的貢獻,甲醛柱濃度的變化與汽車保有量、地區生產總值等經濟發展要素呈現正相關關系,能源消耗總量與工業廢氣排放總量的增加與甲醛柱濃度增長密切相關,人為因素是甲醛柱濃度變化的主要原因.

OMI；甲醛柱濃度；時空分布；影響因素；珠江三角洲；相關分析

甲醛是具有強烈氣味的刺激性氣體,略重于空氣,是一種揮發性有機化合物,也是室內環境的主要污染物之一[1-2].甲醛(HCHO)主要分布于與人類生產生活密切聯系的對流層中,是對流層大氣中主要的揮發性有機物(VOCs)之一[3],在城市大氣光化學反應中發揮著重要作用,是大氣光化學煙霧的重要組成成分[4-6].甲醛對人體健康具有很大的危害,目前甲醛已被世界衛生組織(WHO)認定為第一類致癌物——對人致癌物質[7-8].大氣中甲醛來源比較廣泛,甲醛在全球范圍內的最大來源為甲烷光化學氧化,但甲醛在年內濃度的變化卻是由非甲烷類有機物光化學氧化而來,夏天53%的甲醛來源于碳氫化合物的氧化,19%來源于甲烷的氧化,16%為VOCs的氧化,由有機過氧化合物和乙醛氧化來源的分別為3%和7%[9-12].除此之外,某些天然生成過程以及植物排放也會對甲醛濃度具有一定貢獻[13-15].而以甲醛為原料的化工廠,化石燃料及生物質的未完全燃燒、工業廢氣與機動車尾氣的排放以及室內裝飾裝修材料的甲醛釋放等人為因素也會對甲醛的濃度變化產生一定影響[16-17].

目前,國內對甲醛的研究主要為分析室內及公共場所等小范圍內大氣中甲醛的濃度及時空分布狀況[1],而應用遙感手段進行長時間序列、大尺度范圍下的甲醛時空分布研究較少, Zhu等[18]與Marais等[19]分別利用OMI衛星以及相關監測數據對大范圍內甲醛時空分布進行了分析,并判斷出生物異戊二烯以及揮發性有機物(VOCs)對HCHO的貢獻,同時也證明了溫度是控制HCHO釋放量的重要氣象要素; Bauwens等[20]利用OMI數據,結合全球IMAGESV2模型、GOME-2衛星數據、各地氣象數據等,利用不同數據方法獲得的HCHO數據,對不同地區HCHO的來源及控制因素進行了判斷與分析.遙感監測可以實現對大尺度范圍下對流層甲醛污染監測的同時也更易獲得長時間序列下甲醛濃度的動態遷移過程,從而對大氣中甲醛氣體的時空分布研究提供便利條件[21-22].珠江三角洲作為我國最具有活力的經濟區之一,NO2和O3等污染物都受到關注與報道[23-25],但對于甲醛的研究未見,本文利用OMI傳感器解譯出的甲醛柱濃度數據,分析珠江三角洲對流層甲醛垂直柱濃度的時空分布特征,并對其影響因子之間的關系進行了研究.

1 研究方法

1.1 研究區概述

珠江三角洲(珠江三角洲),其地理位置為(112°E~115.5°E,21.5°N~24°N)(圖1).該地區海拔位于200m以下,中心地勢較為平坦,氣候以南亞熱帶季風氣候為主[26],終年溫暖濕潤,年均溫在21~23℃之間,最冷月1月均溫13~15℃,7月均溫28℃以上,6~10月,雨熱同期,高溫高濕,年均降水量1500mm以上,土壤肥沃,河道縱橫.夏季南風,冬季北風,植被為亞熱帶常綠闊葉林.

珠江三角洲土地面積為42831.5km2,包括廣州,深圳等九個城市,是我國乃至亞太地區最具有活力的經濟區之一,同時也是世界人口和面積最大的城市群,作為全省交通運輸最繁忙、最發達的地區,現已形成了以廣州為中心,鐵路、公路、水路、民航等多種運輸方式相配合,銜接港澳、溝通全省和全國的較為發達的綜合交通網絡.

圖1 研究區域示意

1.2 數據來源及方法

1.2.1 數據來源 目前用于測量大氣甲醛柱濃度的遙感手段主要有GOME-1(分辨率為320km× 40km,發射于1996年)[27]、GOME-2(分辨率為80×40km,發射于2006年)[28]、SCIAMACHY(分辨率為60×30km,發射于2002年)[29]以及OMI(分辨率為24×13km,發射于2004年)[30],本文將這4種技術所測得的甲醛柱濃度值進行對比,最終選擇精度與分辨率最高的OMI進行數據獲取[27].本文所下載的2009~2016年珠江三角洲對流層L2-v003HCHO的逐日數據,其是源自美國宇航局(NASA)從臭氧檢測儀(簡稱OMI,表1)所監測的數據,采用HDF-EOS格式進行存儲.OMI對流層HCHO柱濃度產品已開展了很多驗證評價工作,與地基Brewer MK3分光計數據產品相關系數可達0.9以上,保證了數據產品的精度與可靠性[5,31].

本文在研究過程中所用歸一化差值植被指數數據集(NDVI) 為美國宇航局(NASA)根據統一算法開發的MODIS 3級陸地產品,本文所用的是MOD13A3,為Terra極地軌道環境遙感衛星全球月尺度的３級產品,其采用正弦曲線投影,空間分辨率為1000m,時間分辨率為月尺度,該數據產品圖像具有現勢性好、覆蓋面廣,多光譜特性和不接受任何費用等特點,加之對氣溶膠、云、水等進行了處理,提高了數據的精度,所以被廣泛應用于反映植被覆蓋度的變化;本文選取時段為2009年1月~2016年12月,研究區域為h28v6條帶.

表1 OMI主要參數

珠江三角洲的經濟發展、能源消耗的各個指標以及汽車保有量等統計數據來源于廣東省統計年鑒,氣象數據(包括氣壓、溫度等)來源于中國氣象科學數據共享平臺.(http://data.cma.cn/).

1.2.2 研究方法 本文選取2009~2016年對流層L2-v003HCHO的逐日數據對研究區域甲醛柱濃度進行研究,運用HDFView軟件提取“Longtitude” (經度)“Latitude”(緯度)、ColumnAmount” (對流層甲醛垂直柱濃度,單位為1015molec/cm2) “AMFCloudFraction”(云量)4個要素,并將提取邊界經緯度范圍擴大0.3°,篩選出AMFCloudFraction<0.2的柱濃度數據,篩選云量是由于云會影響云層下甲醛濃度的精度,篩選可以提高數據精度[32].然后對極大極小濃度值進行剔除后,將所得點數據進行克里金插值、并通過數據掩膜裁剪出研究區域的月、季、年濃度值,進行研究區甲醛柱濃度變化和空間分布特征分析.該數據產品單位為molec/cm2,本文為方便統計,將甲醛柱濃度單位定義為1015molec/cm2.

本文選取的MOD13A3數據,利用MRTSwath對數據進行處理,經過格式轉換和重投影,應用GIS軟件裁剪出珠江三角洲區域,得到投影系統為WGS84的NDVI序列,并對所得圖像進行空間拼接、合成跟裁剪,最終獲得歸一化差值植被指數NDVI[33-35].

本研究采用基于像元的空間分析法,分析了甲醛柱濃度與植被覆蓋度的相關性.

甲醛柱濃度植被覆蓋度的相關系數計算公式[36]如下:

2 結果與討論

2.1 珠江三角洲對流層甲醛柱濃度時間變化特征

2.1.1 甲醛柱濃度年際變化 如圖2所示,8a來珠江三角洲年均甲醛柱濃度值呈現波動變化,年均值由2009年的13.05×1015molec/cm2上升到2016年的13.70×1015molec/cm2,其中2011、2012與2015年均呈現下降趨勢,最大降速為5.8%(2012年),其他年份為上升趨勢,最大增速為6.3%(2016年).8a間甲醛濃度值2012年最低,為12.49×10×1015molec/ cm2;2016年最高,為13.70×1015molec/cm2;2010~ 2012年呈現明顯下降趨勢,這主要是由于:①自2010年起珠江三角洲經濟進入減速換擋時期,根據《廣東省統計年鑒》2010~2012年珠江三角洲生產總值增長率分別為:11.2%、9.9%、8.1%;②根據《廣東省國民經濟和社會發展第十二個五年規劃綱要》,為響應控制大氣污染態勢,建設綠色廣東的號召,珠江三角洲區域積極淘汰落后水泥、鋼鐵產能,關停小火電,推進節能減排工作,大力消減了能源消耗產業的甲醛排放.珠江三角洲甲醛柱濃度的年均值范圍為10.79×1015~14.62×1015molec/cm2,與單源對中國中東部年均甲醛柱濃度范圍6.5×1015~19.0× 1015molec/cm2較接近,與謝順濤等[37]對西部的蘭州市年均甲醛柱濃度范圍7.95×1015~13.32× 1015molec/cm2相比較高,且蘭州市年均最高值與珠江三角洲8a平均值相近.

圖2 2009~2016年甲醛柱濃度年際變化

2.1.2 甲醛柱濃度季節變化 圖3為2009~2016年珠江三角洲不同季節(春季:3~5月、夏季:6~8月、秋季:9~11月、冬季:12~次年2月)甲醛柱濃度值變化.珠江三角洲甲醛柱濃度表現為夏季始終高于其他3季,冬季最低,這可能與甲醛在高溫高濕的條件下有很強散發力有關.8a季節平均值由大到小依次為:夏季(14.69×1015molec/cm2)>秋季(13.42× 1015molec/cm2)>春季(13.31×1015molec/cm2)>冬季(10.09×1015molec/cm2).其中春秋兩季濃度差較小,在8a中交替且幅度變化較大,如2010年、2013年與2016年表現為春季高于秋季,其余年份均秋季高于春季.另外,2014年春夏秋3季甲醛柱濃度差距量最小,與當時三季溫度接近有關.四季中春夏季甲醛柱濃度變化趨勢相似,珠江三角洲夏季甲醛柱濃度在13.85×1015~16.06×1015molec/cm2之間變化,最低值出現于2012年,最高值出現在2016年,其中2010~ 2012年呈現下降趨勢,2012~2016呈現上升趨勢,8a來珠江三角洲夏季甲醛柱濃度增長7.5%,最大增速為2015~2016年15.3%.珠江三角洲春季甲醛柱濃度在12.09×1015~14.49×1015molec/cm2之間,最低值出現在2009年,最高值出現在2016年,8a來春季甲醛柱濃度增幅為19.9%.珠江三角洲秋季甲醛柱濃度最大值出現于2011年,為14.41×1015molec/cm2,最小值出現于2013年為12.20×1015molec/cm2,8年來秋季甲醛柱濃度降幅為19.9%;珠江三角洲冬季甲醛柱濃度呈現小幅上升趨勢,最大值出現于2015年為11.60×1015molec/cm2,最小值出現于2012年為× 1015molec/cm2,8a間冬季甲醛柱濃度上升1.5%.

圖3 2009~2016年甲醛柱濃度季節變化

圖4 2009~2016年甲醛月均變化

2.1.3 甲醛柱濃度月變化 圖4中年內月際變化幅度較大但具有明顯的月變化周期,呈現單峰結構.每年逐月甲醛柱濃度先增大后減小,即甲醛柱濃度自1月不斷上升,于每年6~8月達到峰值,8a均值為15.17×1015molec/cm2,明顯高于年平均值,其中最高月出現于6月.9月后開始下降,最低月出現于12月,年內谷值出現于每年12~次年2月為10.07× 1015molec/cm2,均低于8a平均值.甲醛月均變化圖與氣溫變化呈現出高度一致性,在6~8月太陽輻射較強,氣溫較高,大氣中光化學氧化作用強烈易形成甲醛[38];加之交通、工業排放和生物源燃燒等釋放,夏季偏南風造成污染物向東北的大氣輸送與匯聚,冬季以東北風為主的大氣輸送利于污染物的擴散[39].

圖5采用正弦模型模擬對流層HCHO柱濃度變化趨勢特征,如果人類影響較弱,人為排放源較少,模擬圖則呈現以自然源排放為主的周期性規律,濃度值為正弦曲線,模擬精度較高[40-41].相反,在人類活動頻繁的區域則模擬效果較差,甲醛柱濃度正弦變化規律不明顯.由圖5可見,珠江三角洲96個月的逐月濃度值正弦擬合程度較差,說明該區域甲醛柱濃度受人為影響更大,自然影響較弱.96個月月均甲醛柱濃度值在8.74×1015(2012年1月)~18.37×1015molec/ cm2(2011年6月)之間變化.

圖5 2009~2016年甲醛柱濃度96月月均正弦擬合

2.2 珠江三角洲對流層甲醛柱濃度空間分布特征

2.2.1 甲醛柱濃度年際空間分布 為了便于比較和研究,將珠江三角洲甲醛柱濃度分為4個濃度等級,具體為一級:(11~12)×1015molec/cm2,二級:(12~13) ×1015molec/cm2,三級:(13~14)×1015molec/cm2,四級(14~15)×1015molec/cm2.

由圖6可知,珠江三角洲甲醛柱濃度年際分布呈現出明顯的由西北往東南遞減的空間分布特征,其中以肇慶東北大部、佛山北部和廣州西部組成高值區分布中心,以佛山中南部、廣州東南半部和江門西北半部組成三級次高級分布區,以惠州、東莞、深圳、中山、珠海和江門等珠江三角洲近海岸地區為一二級低值濃度區.三四級高值濃度區大多位于廣佛肇經濟圈,佛肇慶經濟圈較高甲醛柱濃度分布與其經濟水平較高呈現出高度的一致性[42],其中三級濃度區域呈現擴散趨勢,由珠江三角洲西北部往東南部近海岸擴散,其擴散趨勢以2015年~2016年最為明顯,因為2015年9月29日珠江三角洲國家自主創新示范區正式獲得國務院批復,九市協力打造該區域的轉型升級引領區,最終實現我國的開放創新先行區、協同創新示范區、創新創業生態區和國際一流的創新創業中心的目標,珠江三角洲地區2016年生產總值達47897.25億元,比2015年(下同)增長了45.8%,占全省79.1%;人均GDP達84563元,折合13454美元,增長29.6%.

2.2.2 甲醛柱濃度季節空間分布 將甲醛柱濃度季均值分為4個濃度等級,具體為一級:(9~11)× 1015molec/cm2,二級:(11~13)×1015molec/cm2,三級: (13~15)×1015molec/cm2,四級(15~17)×1015molec/cm2.如圖7,與年際空間分布的共同點是甲醛柱濃度西北部高于東南部,且具有明顯的四季變化規律.春秋季,甲醛柱濃度整體呈現二級與三級中值區分布,二級與三級具有明顯的分界線,珠江三角洲中西北部呈現三級分布,珠江三角洲近海岸區域呈現二級分布;夏季是四季中甲醛柱濃度最高的季節,夏季除惠州東部沿海小塊地區呈現二級柱濃度等級外,其他區域均呈現三四級濃度,其中四級濃度高值區主要以廣佛肇經濟圈為主;冬季為四級濃度最低的季節,珠江三角洲中西部為二級濃度分布,其余為一級濃度分布.

圖6 2009~2016年甲醛柱濃度年際空間分布

圖7 2009~2016年甲醛柱濃度季節空間分布

珠江三角洲對流層甲醛柱濃度季節空間分布從整體來看,四季濃度分布皆呈現自東南至西北遞增趨勢,夏季甲醛柱濃度高值分布區域最大,高值濃度區域主要分布于珠江三角洲西北部,冬季甲醛柱濃度低值分布區域最大,這可能是由于夏季氣溫高,冬季氣溫低,夏季植被的生物量最高,冬季植被的生物量較低有關.

2.2.3 甲醛柱濃度月均值空間分布 如圖8,甲醛柱濃度月均值呈現出明顯的月份變化規律.將甲醛柱濃度月均值分為5個濃度等級,具體為一級:(8~10) ×1015molec/cm2,二級:(10~12)×1015molec/cm2,三級: (12~14)×1015molec/cm2,四級(14~16)×1015molec/m2,五級:(16~18)×1015molec/m2.1月、2月與12月僅出現一二級低值區且主要以二級濃度為主,一級濃度分布于惠州北部、惠州東南部沿海地區以及肇慶西部等小范圍內,整體以二級濃度為主;3月甲醛柱濃度整體呈現二級三級分布且主要以三級分布為主,與1月2月相比,一級濃度消失的同時三級濃度出現說明甲醛柱濃度呈現上升趨勢,其中二級分布于珠江三角洲東南部沿海以及惠州廣州北部交界區域,其他區域為三級分布;4月、5月甲醛柱濃度整體呈現三級四級濃度分布,與3月相比,三級濃度分布擴散至珠江三角洲全部區域,4月在肇慶區域出現四級濃度分布且在5月擴散至廣佛肇經濟圈大部分區域;6月為甲醛柱濃度月均值最大值月份,甲醛柱濃度分布整體呈現四級五級分布,與5月相比四級濃度區繼續擴散至珠江三角洲大部分區域,并出現五級濃度區,三級濃度區域僅出現于惠州、深圳、江門與珠海近海岸地區;7月甲醛柱濃度整體變化不大,主要呈現四級分布,其中五級僅出現于肇慶北部與佛山江門交界處,四級濃度分布略微增加,珠江三角洲近海岸地區與6月相比保持三級分布;8月、9月與10月甲醛柱濃度整體呈現三級四級分布,四級濃度分布開始縮小但保持以佛肇慶經濟圈為核心區域,三級濃度區域整體保持不變.11月四級高值濃度區急劇縮小僅出現于肇慶北部,甲醛柱濃度整體呈現三級分布,二級濃度區再次出現并分布于珠江三角洲沿海近海岸地區.

圖8 2009~2016年甲醛柱濃度月均值空間分布

2.3 對流層甲醛柱濃度的影響因素分析

HCHO柱濃度在空間與時間上分布與變化特征不僅與海陸風(北風)、氣溫、氣壓、植被等自然要素有著密切關系,更多受到人類活動影響,如經濟發展、工業排放及能源消耗等因素.本文從氣象要素與人為要素兩方面對甲醛柱濃度影響因素進行分析.

2.3.1 氣象要素 影響甲醛柱濃度的氣象因素較多,除了風的影響以外,本文通過從其他氣象因素中篩選出相關性最高的氣溫與氣壓因素,與珠江三角洲5個主要城市(廣州、深圳、中山、珠海和東莞)進行相關性分析,具體見表2.

表2 甲醛柱濃度與氣象要素的相關性

注:**.在0.01水平(雙側)上顯著相關;*.在0.05水平(雙側)上顯著相關.

由表2可知,珠江三角洲主要城市甲醛柱濃度與溫度呈現顯著正相關,相關系數分別為0.892、0.864、0.839、0.828、0.716,當溫度上升時,促進甲烷與揮發性有機化合物發生強烈氧化反應,影響甲醛光解消除速率,導致甲醛產生速率高于甲醛光解消除速率[43],所以珠江三角洲夏季柱濃度高于其他三季;冬季由于溫度低于其他三季,甲醛氧化反應低于其他三季,同時在冬季南方不集中供暖,煤炭消耗量低,故冬季甲醛柱濃度明顯低于其他三季.

珠江三角洲主要城市甲醛柱濃度與氣壓呈現顯著負相關,相關系數為-0.843、-0.814、-0.818、-0.836、-0.695,在低氣壓時,珠江三角洲地區大氣通常是由熱帶氣旋外圍下沉影響,出現高溫和高日照數天氣,促進光化學反應生成甲醛;當出現高氣壓時,由于北方冷空氣南下,氣溫下降,不利于甲醛生成,且冷空氣南下時的大風有利于大氣污染物的運輸與擴散,因此甲醛濃度較低[44].此外,5市甲醛柱濃度與氣象要素的相關性在空間上呈現一定的差異性.

2.3.2 植被因素 珠江三角洲常年氣候溫暖濕潤,光熱資源充足,具有豐富的植物資源,這些植物可排放出大量的異戊二烯,異戊二烯作為植物排放量最大的VOCs對大氣甲醛的產生具有重要貢獻[45].本文引用植被指數(NDVI)用于反映植物覆蓋度、季節動態變化、長勢等特征,來反映植被生長狀態和空間密度.圖9為甲醛柱濃度平均值空間分布圖與植被覆蓋度相關性分析,與甲醛柱濃度呈現正相關的區域面積占珠江三角洲總面積49.8%,主要以肇慶、江門與惠州等地區為主,其中肇慶地區多山,中低山丘陵占全市總面積的81%,森林覆蓋度高且城市化水平低,與甲醛柱濃度相關性系數多在0.5~1之間變動,且肇慶甲醛柱濃度在空間分布上呈現高值分布,說明影響肇慶地區甲醛柱濃度的影響因素主要為植被排放.

圖9 2009~2016年甲醛柱濃度平均值空間分布與植被覆蓋度相關性

2.3.3 人為要素 經濟發展要素:珠江三角洲的經濟發展水平主要由地區生產總值來衡量,其中工業產業是主要影響甲醛柱濃度的指標之一,同時機動車尾氣的排放作為甲醛的主要來源之一,也是影響其濃度變化的重要因素,截止2016年末珠江三角洲民用汽車保有量已達到129.8萬輛[46].由于廣東省成功申辦2010年亞運會,分別實施了各項污染減排等控制措施,故本文著重選取2011~2016年汽車保有量、地區生產總值、工業總產值及第二產業為主要因素,從而探究經濟發展要素對對流層甲醛柱濃度的影響.圖10表明,4個經濟發展要素與HCHO柱濃度均呈現顯著的正相關關系,且相關系數2均達到0.8以上,說明經濟發展對大氣甲醛柱濃度的變化具有重大貢獻.

能源消耗要素:煤炭、石油等生物質的未完全燃燒,以甲醛為工業生產原料被工廠使用后釋放的有害氣體等都會影響到對流層甲醛柱濃度的變化[47],因此本文選取珠江三角洲2011~2016年工業廢氣排放總量與能源消耗總量兩個指標進行分析,由圖11可知工業廢氣排放總量與甲醛柱濃度呈正相關關系,且相關系數2為0.62,說明珠江三角洲對流層甲醛柱濃度中工業廢氣排放占有較大比重;能源消耗總量與甲醛柱濃度呈現正相關關系,相關系數R2為0.76,說明隨著能源消耗總量的上升,珠江三角洲HCHO柱濃度呈現上升趨勢.

圖10 甲醛柱濃度與經濟發展要素的相關性

圖11 甲醛柱濃度與能源消耗要素的相關性

3 結論

3.1 8a來珠江三角洲甲醛柱濃度年際變化呈現波動變化趨勢,最小值出現于2012年,最大值出現于2016年,其中最大增速為2015年~2016年,達到了6.3%;季節性變化表現為:夏季>秋季>春季>冬季;年內月際變化幅度較大但具有明顯周期性,呈現單峰結構,最低值出現于12月,最高值出現于6月;珠江三角洲96個月的逐月濃度值正弦模擬擬合程度較低,說明該區域甲醛柱濃度受人為影響更大,自然影響較弱.

3.2 8a來珠江三角洲甲醛柱濃度低值區域主要分布于珠海、東莞,惠州、中山,江門及深圳等珠江三角洲近海岸地區;高值區域主要分布于廣佛肇經濟圈,2016年高濃度區域范圍最大;空間分布整體呈現出明顯的由西北往東南遞減的分布特征;季節空間分布表現為夏季甲醛高濃度區域最大,冬季甲醛低濃度區域最大.

3.3 在不同的時空尺度下,溫度與氣壓等要素對甲醛柱濃度的影響有一定的差異性,與溫度的相關性表現為:廣州>深圳>中山>珠海>東莞,與氣壓的相關性絕對值表現為:廣州>珠海>中山>深圳>東莞;在植被量豐富且城市化水平低的地區,甲醛濃度受植被因素影響較大.

3.4 甲醛柱濃度的變化與汽車保養量、地區生產總值、工業總產值及第二產業產值等經濟發展要素呈現正相關關系,能源消耗總量與工業廢氣排放總量的增加也會對甲醛柱濃度增長起促進作用;此外,珠江三角洲不同區域甲醛柱濃度空間分布與當地的經濟水平、煤炭工業等分布情況以及地形差異有一定的關系.

[1] 王 琨,李玉華,趙慶良,等.室內空氣中的甲醛檢測分析及其預測模型 [J]. 中國環境科學, 2004,24(6):19-22.

[2] Li Z, Schwier A N, Sareen N, et al. Reactive processing of formaldehyde and acetaldehyde in aqueous aerosol mimics: surface tension depression and secondary organic products [J]. Atmospheric Chemistry & Physics Discussions, 2011,11(7):19477-19506.

[3] Liu L, Andreani-Aksoyoglu S, Keller J, et al. A photochemical modeling study of ozone and formaldehyde generation and budget in the Po basin [J]. Journal of Geophysical Research Atmospheres, 2007,112(D22):1-16.

[4] Zhu L, Jacob D J, Mickley L J, et al. Anthropogenic emissions of highly reactive volatile organic compounds in eastern Texas inferred from oversampling of satellite (OMI) measurements of HCHO columns [C]//AGU Fall Meeting. AGU Fall Meeting Abstracts, 2014:114004.

[5] Millet D B, Jacob D J, Boersma K F, et al. Spatial distribution of isoprene emissions from North America derived from formaldehyde column measurements by the OMI satellite sensor [J]. Journal of Geophysical Research Atmospheres, 2008,113(D2):

[6] 林 旭,朱 彬,安俊琳,等.南京北郊VOCs對臭氧和二次有機氣溶膠潛在貢獻的研究 [J]. 中國環境科學, 2015,35(4):976-986.

[7] 鐘 禾.甲醛被列為一類致癌物質[J]. 福建質量管理, 2005,(2): 44-45.

[8] 張 瑩,邵 毅,王式功,等.北京市空氣污染物對呼吸系統疾病門診人數的影響 [J]. 中國環境科學, 2014,34(9):2401-2407.

[9] Nicholas L Boeke, Julian D Marshall, Sergio Alvarez, et al. Formaldehyde columns from the Ozone Monitoring Instrument: Urban versus background levels and evaluation using aircraft data and a global model [J]. Journal of Geophysical Research: Atmospheres, 2011,116(D5):D05303.

[10] Stavrakou T, Müller J F, Smedt I D, et al. Evaluating the performance of pyrogenic and biogenic emission inventories against one decade of space-based formaldehyde columns [J]. Atmospheric Chemistry & Physics Discussions, 2009,9(3):1037-1060.

[11] Bryan N Duncan, Yasuko Yoshida, Jennifer R Olson, et al. Application of OMI observations to a space-based indicator of NO, and VOC controls on surface ozone formation [J]. Atmospheric Environment, 2010,44(18):2213-2223.

[12] 彭 偉,王 瑛,高先池,等.華北農村夏季大氣甲醛濃度變化特征[J]. 環境科學研究, 2016,29(8):1119-1127.

[13] Luecken D J, Hutzell W T, Strum M L, et al. Regional sources of atmospheric formaldehyde and acetaldehyde, and implications for atmospheric modeling [J]. Atmospheric Environment, 2012,47(2): 477-490.

[14] Kaiser J, Wolfe G M, Bohn B, et al. Evidence for an unidentified ground-level source of formaldehyde in the Po Valley with potential implications for ozone production [J]. Atmospheric Chemistry and Physics, 2014,15(3):1289-1298.

[15] 張玉潔,龐小兵,牟玉靜.北京市植物排放的異戊二烯對大氣中甲醛的貢獻[J]. 環境科學, 2009,30(4):976-981.

[16] Guven B B, Olaguer E P. Ambient formaldehyde source attribution in Houston during TexAQS II and TRAMP [J]. Atmospheric Environment, 2011,45(25):4272-4280.

[17] Largiuni O, Giacomelli M C, Piccardi G. Concentration of Peroxides and Formaldehyde in Air and Rain and Gas-Rain Partitioning [J]. Journal of Atmospheric Chemistry, 2002,41(1):1-20.

[18] Zhu L, Mickley L J, Jacob D J, et al. Variability of HCHO over the Southeastern United States observed from space: Implications for VOC emissions[C]// AGU Fall Meeting. AGU Fall Meeting Abstracts, 2012.

[19] Marais E A, Jacob D J, Kurosu T P, et al. Isoprene emissions in Africa inferred from OMI observations of formaldehyde columns [J]. Atmospheric Chemistry & Physics Discussions, 2012,12(3):7475- 7520.

[20] Bauwens M, Stavrakou T, Müller J F, et al. Nine years of global hydrocarbon emissions based on source inversion of OMI formaldehyde observations [J]. Atmospheric Chemistry & Physics, 2016,16(15):1-45.

[21] 張彥軍,牛 錚,王 力,等.基于OMI衛星數據的城市對流層NO2變化趨勢研究[J]. 地理與地理信息科學, 2008,(3):96-99.

[22] 吳炳方,張 淼.從遙感觀測數據到數據產品[J]. 地理學報, 2017,72(11):2093-2111.

[23] 劉顯通,鄭騰飛,萬齊林,等.OMI遙感珠江三角洲城市群NO2的時空分布特征及人類活動影響分析[J]. 熱帶氣象學報, 2015,31(2): 193-201.

[24] 陳 皓,王雪松,沈 勁,等.珠江三角洲秋季典型光化學污染過程中的臭氧來源分析[J]. 北京大學學報(自然科學版), 2015,51(4): 620-630.

[25] 楊柳林,曾武濤,張永波,等.珠江三角洲大氣排放源清單與時空分配模型建立 [J]. 中國環境科學, 2015,35(12):3521-3534.

[26] 竺可楨.中國近五千年來氣候變遷的初步研究[J]. 中國科學, 1973,(2):168-189.

[27] Chance K, Palmer P I, Spurr R J D, et al. Satellite observations of formaldehyde over North America from GOME [J]. Geophysical Research Letters, 2000,27(21):3461-3464.

[28] Callies J, Corpaccioli E, Eisinger M, et al. GOME-2-Metop's second-generation sensor for operational ozone monitoring [J]. ESAbulletin, 2000,102:28-36.

[29] 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.

[30] Levelt P F, Oord G H J V D, Dobber M R, et al. The ozone monitoring instrument [J]. IEEE Transactions on Geoscience & Remote Sensing, 2006,44(5):1093-1101.

[31] Baek K H, Kim J H, Park R J. Validation of OMI HCHO data and its analysis over Asia [J]. Science of the Total Environment, 2014,490: 93-105.

[32] Millet D B, Jacob D J, Turquety S, et al. Formaldehyde distribution over North America: Implications for satellite retrievals of formaldehyde columns and isoprene emission [J]. Journal of Geophysical Research Atmospheres, 2006,111(D24):4057-4065.

[33] 葉勤玉,高陽華,楊世琦,等.基于MODIS數據的重慶市植被覆蓋度時空變化分析 [J]. 高原山地氣象研究, 2016, 36(2):53-58.

[34] 王宇航,趙鳴飛,康慕誼,等.黃土高原地區NDVI與氣候因子空間尺度依存性及非平穩性研究 [J]. 地理研究, 2016,35(3):493-503.

[35] 陳艷英,游揚聲,唐云輝,等.重慶市植被指數與高度/坡度的關系研究 [J]. 西南師范大學學報(自然科學版), 2016,41(11):122-129.

[36] 穆少杰,李建龍,陳奕兆,等. 2001~2010年內蒙古植被覆蓋度時空變化特征 [J]. 地理學報, 2012,67(9):1255-1268.

[37] 謝順濤,巨天珍,張惠娥.基于OMI數據的蘭州地區對流層甲醛時空變化研究[J]. 環境科學學報, 2017,37(11):4253-4261.

[38] 陳月霞,曹 琳,楊淑萍.氣溫對大氣主要污染物的影響[J]. 現代農業科技, 2011,(10):277-278.

[39] 何啟超,蔡旭暉,宋 宇,等.珠江三角洲區域大氣擴散和輸送特征診斷分析[J]. 北京大學學報(自然科學版), 2013,49(6):945-954.

[40] 張興贏,張 鵬,張 艷,等.近10a中國對流層NO2的變化趨勢、時空分布特征及其來源解析[J]. 中國科學(D輯:地球科學), 2007, (10):1409-1416.

[41] 李 兵.基于衛星遙感數據的甘肅省對流層NO2時空分布特征研究[D]. 蘭州:西北師范大學, 2016.

[42] 吳 丹,吳仁海.不同地區經濟增長與環境污染關系的VAR模型分析——基于廣州、佛山、肇慶經濟圈的實證研究[J]. 環境科學學報, 2011,31(4):880-888.

[43] 王少麗,王會祥,劉 斌.北京市大氣甲醛濃度研究[J]. 環境科學研究, 2008,(3):27-30.

[44] 黃 俊,廖碧婷,吳 兌,等.廣州近地面臭氧濃度特征及氣象影響分析 [J].環境科學學報,1-17[2017-12-28].https://doi.org/10.13671/ j.hjkxxb.2017.0256.

[45] Guether A B, Monson R K, Flynt M S. 2, 4Dinitrophenylhydrazones emissions from deciduous trees [J]. Atmospheric Environment, 2000, 34:3027-3032.

[46] 田賀忠,郝吉明,陸永琪,等.中國氮氧化物排放清單及分布特征[J]. 中國環境科學, 2001,21(6):14-18.

[47] 曹 旭.甲醛氣體的危害及污染治理[J]. 環境保護與循環經濟, 2014,34(9):48-50+63.

Remote sensing monitoring of tropospheric HCHO column concentration and influential factors over pearl river delta, China.

XIAN Long1, GE Jian-tuan1,2*, XU Min3, CHEN Xue-ping1, WANG Shuang1, XIE Shun-tao1, HU Wen-wen4

(1.College of Geography and Environmental Science, Northwest Normal University, Lanzhou 730070, China；2.Jiayuguan City Government, Jiayuguan 735100, China；3.College of Environmental and Municipal Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China；4.Heshui County Environmental Protection Agency, Qingyang 745400, China)., 2018,38(9)：3221~3231

Satellite remote Sensing data Retrieval based on OMI, the temporal and spatial distribution characteristics of formaldehyde column concentration in the troposphere and its influencing factors in the Pearl River Delta region from 2009 to 2016were studied. The results show that the concentration change of formaldehyde column in the Pearl River Delta was characterized by fluctuations in the past 8years, with an average annual value of 13.11×1015molec/cm2, the lowest value appeared in 2012 and the highest value appeared in 2016, the maximum reduction rate was 5.8%, the maximum growth rate was 6.3%. The highest in summer and the lowest in winter, the order was summer>autumn>spring>winter, and the monthly variation of formaldehyde in 8a to 96months was larger, showing a single peak structure, with the highest in June every year; The spatial variation was characterized by: The concentration of formaldehyde column decreases from northwest to southeast. Among them, the distribution center of the northeastern part of Zhaoqing, the northern part of Foshan and the western part of Guangzhou constitutes the distribution center of high value area. It consists of the third-level sub-level distribution area in the south-central part of Foshan, the southeastern part of Guangzhou and the northwest half of Jiangmen. In Huizhou, Dongguan, Shenzhen, Zhongshan, Zhuhai and Jiangmen, the coastal areas of the Pearl River Delta are the first-level low-value concentration areas. Meteorological factors such as temperature and pressure affect the formation and distribution of HCHO, and vegetation has a certain contribution to the production of HCHO. There is a positive correlation between the change of formaldehyde column concentration and the factors of economic development, such as automobile ownership, regional GDP, etc, the increase of total energy consumption and total industrial waste gas emission was closely related to the increase of formaldehyde column concentration, and the main reason for the change of formaldehyde column concentration was human factors.

OMI；HCHO column concentrationa；temporal and spatial distribution；influencing factors；Pearl River Delta；correlation analysis

X511

A

1000-6923(2018)09-3221-11

咸 龍(1991-)男,甘肅定西人,西北師范大學碩士研究生,主要研究大氣污染物遙感研究.

2018-01-20

國家重點研發計劃項目(2016YFC0500907)甘肅省重點研發計劃項目(17YF1FA120)

* 責任作者, 副研究員, 2995616335@qq.com