杜君平，朱玉霞，劉　銳,，謝　濤，姚　新

1.中科宇圖(北京)資源環(huán)境科學(xué)研究院，北京 100101

2.中科宇圖天下科技有限公司，北京 100101

臭氧是大氣中重要的微量成分之一，約有90％的臭氧分布在10～50 km高度的平流層大氣中，其余10％分布在對(duì)流層[1]。臭氧對(duì)太陽(yáng)紫外輻射(波長(zhǎng)為0.2～0.29 μm)有強(qiáng)烈的吸收作用：一方面可以阻擋強(qiáng)紫外線(xiàn)到達(dá)地面，保護(hù)地球上的生命；另一方面可以使平流層增溫，對(duì)流層降溫，對(duì)大氣的溫度場(chǎng)和大氣環(huán)流具有重要的作用。此外，臭氧在紅外波段9.6μm附近有一個(gè)很強(qiáng)的吸收帶，大量吸收地球紅外輻射，使低層大氣增溫，是一種重要的溫室氣體。在近地表，臭氧本身還是一種污染氣體，高濃度的臭氧可以對(duì)人和動(dòng)物的呼吸系統(tǒng)造成傷害，并損害植被的健康。由于臭氧總量及其時(shí)空變化對(duì)氣候環(huán)境、生態(tài)環(huán)境以及人類(lèi)生存環(huán)境都會(huì)產(chǎn)生重要的影響，全球臭氧監(jiān)測(cè)日益成為人們關(guān)注和研究的熱點(diǎn)之一。

世界氣象組織(WMO)于20世紀(jì)70年代啟動(dòng)了“全球臭氧觀(guān)測(cè)系統(tǒng)”(GO3OS)，并于1989年開(kāi)始組建全球大氣監(jiān)測(cè)網(wǎng)(GAW)，建立了包括溫室氣體、臭氧等的資料中心[2]。這些地基觀(guān)測(cè)數(shù)據(jù)有力地促進(jìn)了人們對(duì)臭氧的研究和認(rèn)識(shí)，為臭氧的模型反演提供了重要的地面驗(yàn)證數(shù)據(jù)，但其局限于單點(diǎn)觀(guān)測(cè)，無(wú)法獲得臭氧的空間分布。

隨著遙感空間技術(shù)的發(fā)展，衛(wèi)星遙感為獲取全球或區(qū)域尺度的臭氧監(jiān)測(cè)數(shù)據(jù)提供了可能。自1978年美國(guó)Nimbus-7衛(wèi)星發(fā)射以來(lái)，臭氧總量繪圖譜儀(TOMS)和太陽(yáng)后向散射紫外光譜儀(SBUV)在逐日臭氧分布監(jiān)測(cè)，尤其是極地臭氧變化監(jiān)測(cè)中發(fā)揮了巨大作用[3]：美國(guó)NOAA系列衛(wèi)星(NOAA-9、NOAA-11、NOAA-14等)自1985年開(kāi)始搭載SUBV/2，進(jìn)行可操作性的臭氧監(jiān)測(cè)；1996年7月美國(guó)發(fā)射的Earth Probe衛(wèi)星搭載了TOMS，用以接替1993年5月失效的TOMS/Nimbus-7繼續(xù)臭氧總量的監(jiān)測(cè)；隨著2006年11月TOMS-EP的失效，2007年5月長(zhǎng)達(dá)30 a的TOMS觀(guān)測(cè)計(jì)劃宣布結(jié)束。1995年4月21日，歐洲ERS-2衛(wèi)星發(fā)射升空，星上搭載的全球臭氧監(jiān)測(cè)試驗(yàn)儀(GOME-1)可以對(duì)臭氧、SO2等痕量氣體進(jìn)行全球監(jiān)測(cè)[4]，該儀器于2011年7月停止運(yùn)行。掃描成像大氣吸收光譜儀(SCIAMCHY)是歐洲空間局(ESA)2002年3月1日發(fā)射的大型環(huán)境衛(wèi)星ENVISAT-1上搭載的十大載荷之一，較GOME-1在光譜范圍和空間分辨率上都有較大提高，在全球痕量氣體監(jiān)測(cè)方面發(fā)揮了重要作用[5]，但其在2012年4月失效。

目前用于臭氧總量監(jiān)測(cè)的在軌衛(wèi)星傳感器主要有：美國(guó)EOS-AURA衛(wèi)星搭載的臭氧監(jiān)測(cè)儀(OMI)，2007年4月運(yùn)行[6]；歐洲 MetOP系列衛(wèi)星搭載的GOME-2，其中MetOP-A和MetOP-B分別于2006年10月和2012年9月成功發(fā)射；中國(guó)FY-3系列衛(wèi)星搭載的紫外臭氧總量探測(cè)儀(TOU)，其中FY-3 A和 FY-3 B分別于2008年5月和2011年11月成功發(fā)射；美國(guó)Suomi NPP衛(wèi)星搭載的臭氧繪圖儀和廓線(xiàn)儀套件(OMPS)，包括下視天底繪圖儀、天底廓線(xiàn)儀和臨邊廓線(xiàn)儀，于2011年10月發(fā)射升空。這些傳感器可以監(jiān)測(cè)全球臭氧變化，為臭氧研究提供更加豐富精確的數(shù)據(jù)源。

為了深入了解最近幾年中國(guó)上空臭氧總量的變化，該研究基于OMI臭氧格網(wǎng)資料，提取了2005—2011年中國(guó)上空的臭氧總量，并對(duì)其空間分布特征、時(shí)間變化規(guī)律及緯向差異進(jìn)行了分析和研究。

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

臭氧總量數(shù)據(jù)來(lái)源于搭載在EOS-Aura衛(wèi)星上的OMI傳感器。Aura衛(wèi)星于2004年7月15日在美國(guó)成功發(fā)射，是一顆由多國(guó)共同研制的近極地、太陽(yáng)同步科學(xué)探測(cè)衛(wèi)星，軌道高度705 km，約在當(dāng)?shù)貢r(shí)間13:45過(guò)境。OMI是繼GOME和SCIAMACHY之后的新一代大氣成分探測(cè)傳感器，由荷蘭航空局和芬蘭氣象所提供，波長(zhǎng)范圍270～500 nm，平均波譜分辨率0.5 nm，星下點(diǎn)空間分辨率13 km×2 km，掃描幅寬2 600 km，覆蓋全球只需1 d。

該研究使用的是OMI Level 3全球臭氧總量格網(wǎng)產(chǎn)品OMITO3e(單位DU)，空間分辨率達(dá)到0.25°(緯度)×0.25°(經(jīng)度)，存儲(chǔ)格式為HDF-EOS 5。它采用的是OMI-TOMS(TOMS V8)算法反演臭氧總量，該算法僅使用2個(gè)波長(zhǎng)(317.5、331.2 nm，在太陽(yáng)天頂角度較高時(shí)使用331.2、360 nm)，較長(zhǎng)的波長(zhǎng)用于提取地表反射率，較短的波長(zhǎng)由于可以被臭氧強(qiáng)烈吸收而用于提取臭氧柱的總量，數(shù)據(jù)平均精度優(yōu)于1％，并且相對(duì)其他數(shù)據(jù)產(chǎn)品對(duì)太陽(yáng)天頂角度和緯度無(wú)明顯的依賴(lài)性[7]。

由于OMITO3e僅是從2004年10月開(kāi)始工作，因此研究只對(duì)2005—2011年臭氧總量進(jìn)行分析?；贓NVI/IDL批量讀取HDF原始數(shù)據(jù)，轉(zhuǎn)換成通用的TIFF格式(0.25°×0.25°)，提取2005年1月1日—2011年12月31日中國(guó)上空的臭氧總量信息，在逐日數(shù)據(jù)的基礎(chǔ)上合成得到逐月和逐年的臭氧總量，并計(jì)算其緯向偏差，進(jìn)而分析中國(guó)臭氧總量的時(shí)空變化。某個(gè)網(wǎng)格(lat_i，lon_j)的臭氧總量緯向偏差Dlat_i，lon_j，通過(guò)式(1)計(jì)算得到。

(1)

2　中國(guó)臭氧總量時(shí)空分布特征

2.1臭氧總量空間分布特征

通過(guò)計(jì)算發(fā)現(xiàn)，2005—2011年，中國(guó)各年年均臭氧總量空間分布、緯向差異分布類(lèi)似，最終選取多年平均狀況進(jìn)行分析，并以2011年為代表分析臭氧總量的逐月空間分布特征，詳見(jiàn)圖1、圖2和圖3。

圖1　中國(guó)2005—2011年平均臭氧總量空間分布

圖2　中國(guó)2005—2011年平均臭氧總量緯向偏差分布

圖3　中國(guó)2011年1—12月臭氧總量逐月空間分布

從圖1可以看出，2005—2011年中國(guó)年均臭氧總量從西南到東北隨著緯度的增加而增加，呈現(xiàn)出明顯的緯度地帶性；峰值出現(xiàn)在中國(guó)東北地區(qū)，最大值約為372DU；青藏高原(78°E～103°E，28°N～38°E)、中國(guó)東南部大氣臭氧含量較低，最小值出現(xiàn)在青藏高原南部、海南島和臺(tái)灣島局部，約256DU。

從圖2可見(jiàn)，隨著緯度的增加，臭氧總量從西南到東北緯向偏差不斷增加；低緯地區(qū)緯向偏差為負(fù)，即該區(qū)域臭氧總量小于緯圈的平均值；中緯地區(qū)緯向偏差大致以105°E經(jīng)線(xiàn)為界，東部為正、西部為負(fù)，間接說(shuō)明該區(qū)域中國(guó)東部臭氧總量大于中國(guó)西部臭氧總量，這與之前的一些研究結(jié)果一致[8-9]。研究表明，不同緯度帶臭氧總量的差異，主要與太陽(yáng)輻射的緯度差異有關(guān)；而同一緯度臭氧總量分布不均勻，主要是由于行星波擾動(dòng)對(duì)臭氧的動(dòng)力輸送作用[10-11]導(dǎo)致的，該擾動(dòng)作用隨著緯度的增加而增強(qiáng)，且在冬季較強(qiáng)、夏季較弱。由平流層行星波作用形成的東亞臭氧高值中心(約在120°E～150°E, 50°N～70°N范圍)影響了中國(guó)中緯度地區(qū)臭氧總量的空間分布格局，是中國(guó)東北部臭氧高值區(qū)出現(xiàn)的直接原因[10,12]。

在青藏高原(78° E～103°E，28°N～38°N)、天山山脈地區(qū)出現(xiàn)明顯的緯向偏差低值區(qū)，尤其在青藏高原緯向偏差最低可至-25 DU(圖 2)。

進(jìn)一步分析2011年逐月臭氧總量變化(圖3)，從1—4月，中低緯度月均臭氧總量空間分布基本與緯圈平行；而從5月開(kāi)始青藏高原臭氧低谷開(kāi)始出現(xiàn)，并不斷發(fā)展，到10月之后逐漸消退。周秀驥等認(rèn)為，動(dòng)力和熱力過(guò)程是形成青藏高原夏季臭氧低值中心的主要作用，而化學(xué)過(guò)程作用較?。?月，南亞高壓從低緯地區(qū)經(jīng)中印半島向西北方向移動(dòng)，帶來(lái)低緯地區(qū)低濃度的臭氧，在青藏高原開(kāi)始形成臭氧低值中心；夏季，青藏高原上空盛行上升氣流運(yùn)動(dòng)，對(duì)流層低濃度臭氧向平流層輸送及低空污染物在平流層引起的物理化學(xué)作用，維持并加強(qiáng)了青藏高原的臭氧低值中心[13-14]。鄒捍等在研究了大尺度山地(青藏高原、落基山脈和安第斯山脈)上空臭氧總量的分布和季節(jié)變化規(guī)律后指出，由于大尺度山地對(duì)大氣的動(dòng)力和熱力作用，相比于同緯度無(wú)山地區(qū)，山地上空臭氧總量在各個(gè)季節(jié)存在著明顯虧損[15]。而從圖 2可知，該山地效應(yīng)在天山山脈和臺(tái)灣山脈也有明顯的體現(xiàn)。

從圖 1和圖 2可以看出，與在高原和山地出現(xiàn)的臭氧低值中心現(xiàn)象相反，中國(guó)四川盆地(103°E～108°E，28°N～33°N)多年平均臭氧總量高于周?chē)貐^(qū)，存在一個(gè)小的臭氧高峰。進(jìn)一步分析2011年臭氧總量的逐月分布發(fā)現(xiàn)，四川盆地臭氧峰值常年存在，但在春季尤為明顯，隨著時(shí)間的推移逐漸減弱(圖3)。這可能是由于西風(fēng)氣流在青藏高原周?chē)a(chǎn)生繞流和爬坡，于東部邊緣較低地形處匯合和下沉，帶來(lái)了高緯度地區(qū)的高濃度臭氧，而盆地地形的下沉和匯合作用更加明顯[9]；同時(shí)，盆地地形不利于氣流的交換，容易累積污染物，而且四川盆地濕度較大，有利于對(duì)流層中發(fā)生光化學(xué)反應(yīng)，產(chǎn)生臭氧[16]，兩者綜合作用導(dǎo)致四川盆地上空的臭氧總量常年高于周?chē)貐^(qū)。由此可見(jiàn)，地形對(duì)臭氧的空間分布具有較大的影響。

由圖 3可知，2011年1—12月，中國(guó)東北部常年為臭氧總量高值區(qū)，而西南部常年為臭氧總量低值區(qū)。臭氧總量在春冬季節(jié)空間變異程度較大，其中1月的月均臭氧總量標(biāo)準(zhǔn)差最高，約為56DU；臭氧總量在夏秋季節(jié)空間變異程度較低，其中8月月均臭氧總量標(biāo)準(zhǔn)差最低，僅11DU。結(jié)合圖4可知，臭氧總量在其他年份也有類(lèi)似的空間變異特征：總體而言，2005—2011年臭氧總量在2月空間離散程度最大，8月空間離散程度最小。

2.2臭氧總量時(shí)間變化特征

圖 4為2005—2011年中國(guó)月均臭氧總量在空間上的均值、最大值、最小值、標(biāo)準(zhǔn)差的逐月變化。從圖4可以看出，2005—2011年，中國(guó)上空月均臭氧總量(對(duì)應(yīng)圖中均值)大致在273～346DU之間呈正弦曲線(xiàn)變化，具有明顯的季節(jié)特征，其中最大值多出現(xiàn)在3月，而最小值多出現(xiàn)在9、10月。

圖4　2005—2011年中國(guó)月均臭氧總量空間統(tǒng)計(jì)特征

為了進(jìn)一步分析中國(guó)上空臭氧總量的季節(jié)變化，計(jì)算了2005—2011年平均的中國(guó)月均臭氧總量(圖 5)。

圖 5　2005—2011年平均的中國(guó)逐月臭氧總量變化

從圖5可以看出，就2005—2011年中國(guó)地區(qū)平均狀況而言，中國(guó)上空臭氧總量在3月達(dá)到峰值(平均約330DU)，之后開(kāi)始下降，并在10月達(dá)到最低值(平均約277DU)，之后臭氧總量開(kāi)始上升，即：從春末到秋初，中國(guó)上空臭氧總量逐漸消耗；從秋末到春初，中國(guó)上空臭氧總量逐漸恢復(fù)，且在春初達(dá)到峰值。臭氧總量的季節(jié)變化主要與太陽(yáng)輻射、大氣環(huán)流的季節(jié)性變化有關(guān)。

沈凡卉等利用TOMS和OMI數(shù)據(jù)分析了1979—2007年中國(guó)臭氧總量變化，指出近30 a中國(guó)臭氧總量年均值為298.61DU，廣州地區(qū)臭氧總量呈現(xiàn)明顯下降趨勢(shì)[17]。鄭向東和韋小麗將香河、昆明、瓦里關(guān)山、龍鳳山地基觀(guān)測(cè)與衛(wèi)星臭氧總量長(zhǎng)期觀(guān)測(cè)進(jìn)行比較，分析得出：1979—2007年期間，4個(gè)站點(diǎn)地區(qū)的臭氧總量在1993年顯著降低，后于1995—1996年逐漸回升，其中瓦里關(guān)站在2001年回升[18]。

圖6為2005—2011年中國(guó)年均臭氧總量變化統(tǒng)計(jì)結(jié)果。

圖6　2005—2011年中國(guó)年均臭氧總量變化

由圖6可以看出：2005—2011年，中國(guó)7 a平均臭氧總量為300.75DU，較30 a均值線(xiàn)高出2DU左右，與上述4個(gè)基站類(lèi)似，中國(guó)臭氧總量近期有回升趨勢(shì)；年均臭氧總量極大值出現(xiàn)在2010年，為310DU；極小值出現(xiàn)在2008年，為296DU；其他年份臭氧總量在298～302DU范圍內(nèi)變化，平均值為299.22DU，略高于30 a均值線(xiàn)。此外，2010年臭氧總量明顯高于其他年份，具體原因有待進(jìn)一步深入研究。

3　結(jié)論

該研究基于OMI數(shù)據(jù)，提取和分析了2005—2011年中國(guó)上空臭氧總量，并對(duì)其時(shí)空分布特征進(jìn)行了分析：

1)中國(guó)上空臭氧總量總體上呈現(xiàn)北高南低、東高西低的特征。在動(dòng)力和化學(xué)過(guò)程作用下，青藏高原常年存在臭氧低值中心，且在夏季尤為明顯；由于地形影響，四川盆地較周?chē)貐^(qū)常年存在一個(gè)較小的臭氧高峰；中國(guó)東北部常年為臭氧高值區(qū)。此外，臭氧總量空間分布在春冬季節(jié)變異程度較大，夏秋季節(jié)離散程度較低。

2)受太陽(yáng)輻射、大氣環(huán)流季節(jié)性變化等因素的影響，中國(guó)臭氧總量的時(shí)間變化呈現(xiàn)明顯的季節(jié)特征。從春末到秋初為臭氧消耗過(guò)程，反之為臭氧恢復(fù)過(guò)程，最大值多出現(xiàn)在3月，最小值多出現(xiàn)在9、10月。

3)2005—2011年，中國(guó)多年平均臭氧總量為301DU，但在2010年存在一個(gè)明顯的高峰，具體原因仍需進(jìn)一步深入研究。

致謝：該研究使用的OMI數(shù)據(jù)產(chǎn)品來(lái)自美國(guó)戈達(dá)德地球科學(xué)數(shù)據(jù)和信息服務(wù)中心(Goddard Earth Sciences Data and Information Services Centre)，在此表示感謝。

參考文獻(xiàn)：

[1]盛裴軒，毛節(jié)泰，李建國(guó)，等. 大氣物理學(xué)[M]. 北京: 北京大學(xué)出版社,2008.

[2]高梅，接連淑，張文華. 氣象科研數(shù)據(jù)共享系統(tǒng)建設(shè)[J]. 應(yīng)用氣象學(xué)報(bào),2004,15(Z1):17-25.

[3]Bracher A, Lamsal L, Weber M, et al. Global satellite validation of SCIAMACHY O3columns with GOME WFDOAS[J]. Atmospheric Chemistry and Physics, 2005,5(2):2 357-2 368.

[4]Burrows J P, Weber M, Buchwitz M, et al. The global ozone monitoring experiment (GOME): Mission concept and first scientific results[J]. Journal of the Atmospheric Sciences,1999,56(2):151-175.

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

[6]Levelt P F, Hilsenrath E, Leppelmeier G W, et al. Science objectives of the ozone monitoring instrument[J]. Geoscience and Remote Sensing, IEEE Transactions on,2006,44(5):1 199-1 208.

[7]Balis D, Kroon M, Koukouli M E, et al. Validation of Ozone Monitoring Instrument total ozone column measurements using Brewer and Dobson spectrophotometer ground-based observations[J]. Journal of Geophysical Research,2007,112(D24): D24S-D46S.

[8]王躍啟，江洪，肖鐘湧，等. 基于OMI數(shù)據(jù)的中國(guó)臭氧總量時(shí)空動(dòng)態(tài)信息提取[J]. 環(huán)境科學(xué)與技術(shù), 2009,32(6):177-180.

[9]韋惠紅，鄭有飛. 我國(guó)臭氧總量的時(shí)空分布特征[J]. 南京氣象學(xué)院學(xué)報(bào),2006,29(3):390-395.

[10]楊景梅，邱金桓. 利用衛(wèi)星資料分析我國(guó)北方東西部臭氧分布差異[J]. 應(yīng)用氣象學(xué)報(bào),2009, 20(1):1-7.

[11]鄭彬，陳月娟，施春華. 平流層臭氧緯向分布季節(jié)變化和行星波的關(guān)系[J]. 高原氣象,2006,25(3): 366-374.

[12]鄭光，吳統(tǒng)文. 北半球臭氧總量與平流層環(huán)流關(guān)系的分析[J]. 高原氣象, 1991,10(3):277-286.

[13]周秀驥，羅超，李維亮，等. 中國(guó)地區(qū)臭氧總量變化與青藏高原低值中心[J]. 科學(xué)通報(bào),1995,40(15):1 396-1 398.

[14]周秀驥，李維亮，陳隆勛，等. 青藏高原地區(qū)大氣臭氧變化的研究[J]. 氣象學(xué)報(bào), 2004,62(5):513-527.

[15]鄒捍，郜永祺. 大尺度山地上空臭氧虧損及其變化趨勢(shì)[J]. 山地研究, 1997, 15(4): 209-213.

[16]肖鐘湧，江洪. 亞洲地區(qū)OMI和SCIAMACHY臭氧柱總量觀(guān)測(cè)結(jié)果比較[J]. 中國(guó)環(huán)境科學(xué), 2011, 31(4): 529-539.

[17]沈凡卉，王體健，沈毅，等. 中國(guó)近30 a 臭氧氣候場(chǎng)特征[J]. 大氣科學(xué)學(xué)報(bào),2011,34(3):288-296.

[18]鄭向東，韋小麗. 中國(guó)4個(gè)地點(diǎn)地基與衛(wèi)星臭氧總量長(zhǎng)期觀(guān)測(cè)比較[J]. 應(yīng)用氣象學(xué)報(bào),2010,21(1):1-10.