王后茂 王詠梅, 王維和 張卓
(1 中國科學院國家空間科學中心,北京 100190;2 國家衛星氣象中心,北京 100081;3 中國科學院大學天文與空間科學學院,北京 100049)
近年來,我國氣溶膠大氣污染事件(如沙塵、灰霾等)的發生越來越頻繁,尤其在我國北方地區,除受沙塵天氣的影響外,受灰霾的影響也較大[1]。灰霾等污染事件伴隨的主要污染物為氣溶膠和痕量氣體(O3、SO2及NO2等),因此對它們進行大范圍的星載監測有助于我國大氣環境污染的預報。
國際上對痕量氣體和氣溶膠的星載監測起步于20世紀70年代。1979年,TOMS(Total Ozone Mapping Spectrometer)探測儀利用大氣紫外后向散射已連續進行近30年的全球大氣臭氧總量反演,并利用吸收性氣溶膠指數進行吸收性氣溶膠監測[2-3]。之后,一系列的星載紫外探測儀器被發射并進行痕量氣體的探測,如歐州空間管理局的GOME-1[4-5]、SCIAMACHY[6]、GOME-2[4],以及美國的OMI[3]和OMPS[7-8]等。此外,這些儀器還進行吸收性氣溶膠指數的反演[9-11],用于吸收性氣溶膠的監測,如灰霾、沙塵、生物燃燒及火山灰等。
由中國科學院國家空間科學中心研制的臭氧總量探測儀(TOU)是我國自主開發研制的用于全球臭氧總量監測的儀器[12],于2008年開始搭載在氣象衛星風云三號A星(FY-3A)發射,之后的FY-3B星和FY-3C星先后搭載TOU發射成功,至今已有三臺TOU儀器在軌運行。TOU的成功發射和運行實現了我國臭氧總量探測從無到有的過程,從根本上改變了我國大氣臭氧探測和研究的現狀,為我國實時地提供全球的臭氧總量分布監測數據。研究表明,TOU與地基測量值之間的平均相對偏差為-0.16%±4.3%,線性相關系數達到0.95,其與美國同類儀器OMI的相對誤差為2.52%[13]。2013年,FY-3/TOU數據被成功地用于紫外吸收性氣溶膠指數的反演,進而進行我國灰霾、沙塵及生物質燃燒等吸收性氣溶膠的監測。
TOU是一臺采用固定光柵、狹縫陣列式Ebert-Fastie單光柵光譜儀[14]。TOU空間掃描瞬時視場角為3.6°,對應的星下點空間分辨率為52.6 km,行掃描像元個數為31個,掃描范圍為±54°,行掃描時間為8.16 s,對應的刈寬為2908 km。TOU探測波段為6個(表1),每個波段的狹縫函數接近于高斯型,帶寬(FWHM)在1.0~1.3 nm[15]。TOU在軌輻射定標和波長定標是分別通過漫反射板和Hg線燈來實現的,其采用三個漫反射板系統監測定標特性的變化,用以減小星上漫反射板衰變對大氣臭氧總量反演精度的影響,漫反射板的輻射定標精度為3%。TOU利用監測296.8 nm Hg線強度相對變化的方法來實現對儀器波長漂移的監測,波長定標精度為0.03 nm[16]。

表1 TOU的光譜設置及特征Table 1 Characteristics of the six TOU bands
FY-3A/FY-3B/FY-3C TOU利用臭氧總量進行全球臭氧分布及變化的監測,利用反演獲得的AAI指數對吸收性氣溶膠(如灰霾、沙塵和生物質燃燒等)的空間分布特征進行監測。本文將對TOU的探測原理、數據產品及應用進行綜合介紹。
不同波長的紫外后向散射強度的差異與臭氧總量之間存在著密切關系,基于這個關系進行TOU臭氧總量的反演。目前,國際上主要有兩種反演方法,分別為TOMS[17-18]和差分吸收DOAS[19-20]算法。其中,TOMS算法有兩個版本,分別為TOMS V7[17]和TOMS V8[18],其被應用于美國NIMBUS-7/TOMS和AURA/OMI 儀器進行臭氧總量的反演;而DOAS算法則被應用于美國AURA/OMI 儀器、歐州空間管理局GOME系列儀器[19]及SCIAMACHY儀器[20]。FY-3/TOU 臭氧總量算法與TOMS V7 基本類似,但略有區別。TOU反演算法在不同的緯度帶采用不同的通道組合來計算臭氧總量初估值,再分別與360 nm 通道組合對輻射定標的影響進行訂正,給出臭氧總量精確值[21]。
利用太陽紫外后向散射反演臭氧和其他痕量氣體時,通常使用通道輻亮度與太陽輻照度的比值即N值

式中,I為通道后向散射亮度,F為相同通道的太陽輻照度。
TOU臭氧總量反演模式包括4個部分:快速正演模式、表面特性計算、臭氧總量初估值計算和臭氧總量精確值計算,具體方法在文獻18中有詳細描述。
吸收性氣溶膠反演基于如下紫外吸收性氣溶膠指數反演公式

式中,AAI為吸收性氣溶膠指數,ALER為地表反照率,I為強度值,λ為波長,下標meas表示衛星測量值,下標calc表示模型計算值。基于輻射傳輸模型,通過改變地表反照率計算獲得波段λ1的多個模擬觀測值,將這些模擬觀測值與λ1的衛星量測值進行比較,找到最接近于所對應的。由于λ1波段的模擬值與觀測值相等,因此式(2)可以簡化為

由于地表同一地物在紫外波段的地表反照率差異較小,可以認為。因此由式(3)可得,基于模型利用求得λ2觀測模擬值,將其與λ2實測值相結合即可求得吸收性氣溶膠指數AAI。
AAI可以較好地進行云、雪或冰等亮背景區域的吸收性氣溶膠的監測。由于灰霾、沙塵等氣溶膠污染事件的發生常伴隨有云,可見光學監測方法無法有效地將云與吸收性氣溶膠進行區分,因而不能很好地進行污染事件的追蹤與監測。而吸收性氣溶膠在紫外波段具有較大的吸收差異,云在紫外波段有強烈的散射作用,吸收差異則非常小。因此,AAI可以很好地將云和吸收性氣溶膠區別開,避免多云和陰天等帶來的監測影響,從而較好地進行云區的灰霾、沙塵等污染事件的監測。
AAI值的大小與大氣中具有吸收作用的氣溶膠含量相關,因此AAI不僅可以定性反映氣溶膠的多少,還可以反映氣溶膠吸收性的強弱,這對氣溶膠類型有很好的指示作用,為氣溶膠污染事件的來源分析提供監測數據。
基于紫外波段的交叉定標方法[22],本研究對TOU L1B數據與歐州空間管理局MetopB/GOME-2(圖1)和美國NPP/OMPS(圖2)對應數據進行在軌交叉比較分析。交叉定標方法主要分為軌道預報、像元時空匹配、觀測幾何匹配、均勻性匹配及光譜轉換與匹配等。其中,通過計算每個像元的時間和地理位置來進行時空匹配,兩臺儀器的像元時間差異小于5 min,像元空間位置差異小于0.5個像元(25 km,TOU空間分辨率為50 km)[20]。比較結果表明,TOU探測數據與國際其他標準儀器的對應數據具有很好的線性相關性(R2一般大于0.96)。如圖1所示,TOU 312 nm、317 nm、322 nm、331 nm及360 nm的輻亮度數據與GOME-2對應波段的數據具有很好的線性相關性(R2>0.98)。此外,交叉定標方程的斜率位于0.93~1.06,且截距均較小,這表明兩臺儀器的探測響應非常一致。
TOU 308 nm、312 nm、317 nm、322 nm、331 nm及360 nm的輻亮度數據與OMPS對應波段的數據具有很好的線性相關(圖2)。此外,交叉定標方程的斜率位于0.86~1.06,方程的截距均較小,這表明兩臺儀器的探測響應非常一致。

圖1 TOU與GOME-2 L1B輻亮度數據交叉定標比較分析Fig. 1 Inter-comparison of radiance between TOU and GOME-2

圖2 TOU與OMPS L1B輻亮度數據交叉比較分析Fig. 2 Inter-comparison of radiance between TOU and OMPS
TOU通過臭氧總量來實現我國對全球大氣臭氧分布的監測。因此,利用在軌測試期間TOU 的觀測數據,進行了全球臭氧總量反演與研究,尤其對南極及青藏高原等地區(圖3)。圖3 為TOU 2008—2016年監測到的南極臭氧洞的年平均變化情況。可以看出,南極臭氧洞在2008—2011年較大,且沒有明顯的變化;2012和2013年,臭氧洞面積縮小很多;2014—2015年,臭氧洞面積擴大,且擴大面積較為明顯;2016年,臭氧濃度增加,臭氧洞面積有所減小。
衛星AAI可以對連續空間分布的霧霾天氣及其移動發展趨勢進行監測,在空間覆蓋方面具有優勢,為霧霾的產生、發展和消退等過程的研究提供實時大范圍的觀測數據。目前,紫外臭氧總量探測儀反演的AAI可以定性反映霧霾的強度分布情況,同時還有助于氣溶膠類型的識別。
TOU是利用331和360 nm通道計算氣溶膠指數,并以此進行灰霾、沙塵及生物質燃燒等吸收性氣溶膠的監測。吸收性氣溶膠指數能夠將云、雪及冰等亮背景與吸收性氣溶膠進行區分,較準確地進行灰霾和沙塵等的發生、發展和影響范圍的遙感監測。
圖4為2015年2月FY-3/TOU對我國灰霾進行連續空間范圍變化監測的結果。由圖可得,TOU可以很好地對我國北方地區產生的灰霾進行大范圍的監測。

圖3 FY-3衛星TOU南極地區臭氧總量監測(2008—2016年)Fig. 3 Ozone measurements for Antarctic from 2008 to 2016 by FY-3/TOU
圖5為2010年4月FY-3/TOU對我國沙塵型氣溶膠進行大范圍連續空間變化監測的結果。由圖可知,TOU可以很好地對沙塵分布等進行大范圍監測。
FY-3A/B/C衛星上搭載的TOU探測儀利用大氣紫外后向散射進行全球臭氧總量的探測,自2008年,3臺TOU已連續進行全球臭氧分布的特征監測,如南極臭氧洞和青藏高原等地區。TOU L1B輻亮度探測值與歐洲空間管理局儀器GOME-2和美國OMPS對應的L1B數據進行比較分析,它們之間具有很好的線性相關;而L2全球臭氧總量數據與美國儀器OMI的相對誤差為2.52%,足以滿足臭氧總量的探測要求。此外,2013年,基于TOU AAI成功實現了我國灰霾、沙塵等吸收性氣溶膠污染事件的監測,并持續為我國大氣環境監測和預報提供實時的監測數據。但TOU的波段僅為6個分離式通道,限制了其進行更多大氣參數的探測。
FY-3/TOU AAI反演選擇的波段為331和360 nm,它們將受到臭氧不同程度的吸收影響,這將導致反演獲得的AAI產生一定的不確定性。此外,由于FY-3/TOU的空間分辨率為50 km×50 km,這不能很好地滿足城市等人口密集區域的環境污染監測的要求。而且,TOU采用的是獨立式多通道進行的探測,而不是采用連續的高光譜進行探測,這導致其不能應用于其他多種痕量氣體(如SO2、NO2)的觀測。因此,在今后的儀器研制中,需要在波段的選取與擴展、空間分辨率的提高等方面進行改進。由于AAI的反演是基于兩個不同波長對氣溶膠吸收的差異,為了提高AAI的敏感性,波段選取應遵循以下原則:

圖4 紫外臭氧總量探測儀在我國華北地區的灰霾動態監測Fig. 4 Smog observation from North China which were completed by using ultraviolet total ozone unit

圖5 2015年2月17日—21日紫外臭氧總量探測儀在我國西北新疆地區的沙塵動態監測Fig. 5 Dust observation from northwest China using the ultraviolet Total Ozone Unit during 17-21 February 2015
1)盡量避開臭氧吸收和Ring效應的影響;
2)盡量避免波段間的地表反射率差異較大;
3)有較強的輻射強度;
4)可以兼顧與可見光氣溶膠光學厚度的結合。
基于以上要求,綜合考慮大氣透過率在紫外波段的變化特性和太陽紫外后向散射輻射強度,在今后的儀器設計中,將采用紫外-可見光連續光譜的方式進行大氣探測,光譜范圍將覆蓋紫外-短波可見光-近紅外波段,以獲得更多的大氣反演參數,包括吸收性氣溶膠指數、氣溶膠類型、氣溶膠光學厚度及痕量氣體,在分辨率的提高上主要分為以下幾個方向:1)光譜分辨率的提高,光譜分辨率將小于0.6 nm,可同時進行大氣氣溶膠相關參數和痕量氣體的觀測;2)空間分辨率的提高,光譜分辨率滿足氣溶膠的探測需求,使得空間分辨率小于7 km;3)光譜分辨率和空間分辨率的同時提高。
Advances in Meteorological Science and Technology2019年1期