翟　墨，原　帥

（佳木斯氣象衛(wèi)星地面站，黑龍江佳木斯154007）

基于EOS/MODIS的黑龍江省大氣可降水量反演方法研究

翟墨，原帥

（佳木斯氣象衛(wèi)星地面站，黑龍江佳木斯154007）

本文以大氣水汽反演為中心，將雙通道和三通道反演算法應(yīng)用于實(shí)際的黑龍江省MODIS影像中，通過對多個通道多組數(shù)據(jù)反演結(jié)果的分析和驗(yàn)證，證實(shí)利用該算法反演大氣水汽是切實(shí)可行的，所獲結(jié)果反映了實(shí)際水汽分布的大致規(guī)律。

大氣水汽含量；透過率；MODIS；衛(wèi)星遙感；反演

1　引言

大氣水汽含量在大氣過程的各種時間空間尺度中飾演著重要的角色。精準(zhǔn)地確定大氣水汽的含量具有十分重要的意義，不僅能夠準(zhǔn)確預(yù)測降水和災(zāi)害性天氣，同時對于全球氣候變化的研究、水循環(huán)、地-氣系統(tǒng)中物質(zhì)和能量的傳輸也有非常重要的作用，因此監(jiān)測水汽是大氣遙感的主要目的之一。而且，由于輻射傳輸在衛(wèi)星遙感反演的過程中會受到水汽的影響，因此為了獲取精確的地表參數(shù)，水汽數(shù)據(jù)的另外一個作用就是可以作為大氣校正的輸入數(shù)據(jù)。據(jù)估計(jì)，利用通用的分裂窗算法反演陸地表面溫度時，1 K精度的地表溫度需要0.6 g/cm2精度大氣水汽數(shù)據(jù)的支持［1］。

縱觀半個世紀(jì)，對于全球網(wǎng)絡(luò)的水汽探測一直是利用傳統(tǒng)的氣球探空方式進(jìn)行，由于觀測儀器多變性而且缺少即時有效的訂正與校準(zhǔn)，因此依然無法在氣候模型中使用這些探空數(shù)據(jù)［2］。另外，通常大氣水汽含量的計(jì)算需要探空資料，而其觀測的水汽數(shù)據(jù)是以點(diǎn)狀分布的，需要對其實(shí)施空間插值擴(kuò)展，才能在接下來的氣候模型中運(yùn)用，但是“以點(diǎn)帶面”又會帶來較大誤差導(dǎo)致模型的模擬結(jié)果缺乏可靠性和準(zhǔn)確性。因此，缺乏長期穩(wěn)定且精準(zhǔn)度較高的水汽數(shù)據(jù)，勢必會嚴(yán)重阻礙水汽在氣候系統(tǒng)中作用與影響的深入研究。

從目前國內(nèi)外的研究現(xiàn)狀來看，隨著對于大氣水汽的相關(guān)研究方法的發(fā)展與改進(jìn)，反演算法的精度的要求也愈來愈高，但迄今為止大多數(shù)大氣水汽的研究是利用AVHRR數(shù)據(jù)［3］和ATSR數(shù)據(jù)進(jìn)行的。而由于MODIS擁有比AVHRR更低的內(nèi)在噪音和相對更高的精密度，且在衛(wèi)星影像的光譜功能和波長上，MODIS的29、31和32通道與AVHRR的4和5通道有極大的相似性。為此，可以利用MODIS反演大氣水汽含量。

2　反演原理

近紅外通道的大氣輻射傳輸方程可表示成［4］：

式中，Lsensor（λ）是衛(wèi)星傳感器接收到的總輻射；Lsun（λ）是大氣上界的入射太陽輻射；τ（λ）是大氣總透過率，即從太陽到地球表面，再從地球表面到傳感器的大氣路徑上的透過率；

ρ（λ）是下墊面反射率；Lpath（λ）是大氣的路徑輻射。由于氣溶膠在近紅外區(qū)域的光學(xué)厚度極小，因此，大氣路徑輻射Lpath（λ）僅為地表反射輻射的百分之幾［5］，可以忽略不計(jì)。

定義星上的反射率為：

MODIS第i（i＝2，5，17，18，19）通道的大氣輻射傳輸方程可用以下公式表示：

上式中，λi為MODIS第i通道的中心波長。

將簡化后的（1）式兩邊同時除以Lsun（λ）ρ（λ）得：

上式中τ（λ）為大氣透過率，對于不同波長的ρ（λ）下墊面反射率基本上不同，因此不可能通過單一通道的輻射值計(jì)算出水汽的透過率。Gao和Kaufman通過實(shí)驗(yàn)證實(shí)，波長在0.85-1.25的范圍內(nèi)，地物的反射率隨波長呈線性的變化［6］，因此地面的反射率不隨波長變化。若地物反射率為常數(shù)，即ρ（λ）=ρ，則由（4）式可得：

即兩通道大氣透過率的比值為兩者的星上反射率之比。當(dāng)大氣路徑輻射Lpath（λ）忽略不計(jì)時，水汽吸收通道上的水汽透過率就可通過吸收通道（即λ1）的反射率與大氣窗通道（即λ2）的反射率相比獲得，且τ（λ2）≈1，此方法稱為二通道比值法，其中以大氣窗0.865 μm（即2通道）為例，具體公式如下：

若地面反射率與波長滿足線性關(guān)系，與雙通道不同，水汽吸收通道上水汽的透射率是通過吸收通道與兩個大氣窗通道之和相比獲得，此法稱為三通道比值法，具體公式如下：

其中：C1=0.8，C2=0.2；λ為MODIS吸收通道17、18、19波長；ρ（λ）為波段反射率。

MODIS近紅外940 nm（17，18，19通道）波長附近具有對水汽的強(qiáng)吸收特性，CO2等氣體與其相比吸收很弱，可以忽略不計(jì)。同時，在通道2和通道5附近窄波段的大氣透過率實(shí)際值大于0.99，接近于1。因此它們可以近似看成大氣窗波段，可通過水汽吸收通道和大氣窗通道的比值來獲取大氣水汽的含量。另外，近紅外波段水汽反演算法只適用于對白天的數(shù)據(jù)進(jìn)行觀測，可獲取全球范圍內(nèi)晴空條件下的陸地、陸地上空和海洋上空云層以上的大氣水汽含量。若在晴空的洋面區(qū)域，只能對太陽耀斑區(qū)的水汽進(jìn)行反演。

整層大氣水汽含量表達(dá)式：

對于符合性地表，α=0.02，β=0.651［7］。

3　基于MODIS影像數(shù)據(jù)的大氣水汽反演應(yīng)用實(shí)例

MODIS的全稱是中分辨率成像光譜儀（Moderate Resolution Imaging Spectrum-radiometer），是搭載在美國EOS系列衛(wèi)星中TERRA與AQUA上的重要傳感器。TERRA和AQUA都是太陽同步近極地軌道，主要任務(wù)都是對地觀測。MODIS是現(xiàn)今全球新一代“圖譜合一”的光譜傳感器，包含36個光譜通道，波譜寬度極廣，在0.4-14 μm的范圍內(nèi)。其中1-19和26通道為可見光和近紅外通道，其他16個通道為熱紅外通道。MODIS數(shù)據(jù)空間分辨率分別為250、500和1000 m，還具有很高的信噪比，掃描寬度為2330 km。每1-2 d將提供覆蓋全球的白天可見光和白天/夜間紅外圖像，觀測過程中數(shù)據(jù)傳輸速度可達(dá)6.1 Mbps，可同時獲得來自大氣、陸地和海洋表面信息。

本文所用MODIS數(shù)據(jù)來源于國家衛(wèi)星氣象中心的衛(wèi)星遙感數(shù)據(jù)服務(wù)網(wǎng)（http：//satellite.cma.gov. cn），且均為AQUA產(chǎn)品，選擇晴空大氣下且質(zhì)量較好的圖像數(shù)據(jù)，分別是2014年6-9月的7幅MOD02影像（6月1日、6月13日、6月29日、7月27日、8月8日、9月12日和9月30日）。實(shí)測資料為黑龍江省4個探空站（哈爾濱、齊齊哈爾、嫩江、伊春）與7幅MODIS遙感影像相對應(yīng)時間的探空實(shí)測數(shù)據(jù)，來源于中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)（http：//cdc.cma. gov.cn）。

3.1數(shù)據(jù)預(yù)處理

圖像預(yù)處理是遙感圖像處理前的必要環(huán)節(jié)，其中包括幾何校正與圖像裁剪等。MODIS水汽產(chǎn)品自帶的經(jīng)緯度文件未經(jīng)過幾何精度校正，首先需要將經(jīng)緯度文件與MODIS影像相匹配，利用ENVI軟件的波段運(yùn)算功能對經(jīng)緯度信息進(jìn)行還原，然后進(jìn)行配準(zhǔn)。

對水汽產(chǎn)品進(jìn)行幾何校正后，以黑龍江省的矢量邊界作為Mask，從MODIS影像中裁剪出研究子區(qū)部分作為下一步的研究數(shù)據(jù)（以2014年6月1日數(shù)據(jù)為例）。

3.2反演算法實(shí)現(xiàn)

在這里將運(yùn)用兩通道與三通道比值法對上一步裁剪出的研究子區(qū)數(shù)據(jù)進(jìn)行大氣水汽反演。依據(jù)公式（6）、（7）進(jìn)行兩通道與三通道比值法反演，運(yùn)用ENVI提供的band math運(yùn)算方法完成算法的實(shí)現(xiàn)。

3.3結(jié)果分析

3.3.1兩通道比值法與三通道比值法反演結(jié)果對比

分別利用兩通道比值法和三通道比值法對黑龍江大氣水汽含量進(jìn)行反演計(jì)算，將反演結(jié)果與黑龍江省四個探空站（哈爾濱、齊齊哈爾、嫩江、伊春）相同時間的探空實(shí)測數(shù)據(jù)進(jìn)行統(tǒng)計(jì)對比分析，結(jié)果表明兩種反演方法結(jié)果與探空實(shí)測數(shù)據(jù)均有顯著的正相關(guān)性，而且兩通道比值法反演結(jié)果與實(shí)測數(shù)據(jù)的相關(guān)系數(shù)（0.944）大于三通道比值法（0.811），說明對于本文研究子區(qū)兩通道比值法反演大氣水汽含量效果優(yōu)于三通道比值法（圖1）。

圖1?。╝）兩通道比值法與（b）三通道比值法大氣水汽含量反演結(jié)果相關(guān)性對比分析圖

3.3.2黑龍江水汽含量的變化特點(diǎn)

3.3.2.1時間變化特征

圖2　黑龍江省不同時次兩通道比值法的大氣水汽含量反演結(jié)果（全省平均值）

由圖2可見，黑龍江省大氣水汽含量從6月初至8月上旬呈上升趨勢，9月則明顯減少，其中黑龍江全省平均水汽含量最大值出現(xiàn)在8月上旬（最大值為3.73 g/cm2），最小值出現(xiàn)在9月下旬（最小值為0.64 g/cm2）。上述水汽含量變化趨勢比較明顯，說明夏季雨季，濕度較大；秋季氣候干燥，濕度較小。

3.3.2.2空間變化特征

以2014年6月1日的MODIS反演圖像為例，利用ENVI軟件提供的感興趣區(qū)（ROI：Region of Interest）對上步得到的反演結(jié)果進(jìn)行分區(qū)域比較分析。

region1的主要地表類型是陸地，包含植被種類繁多，有草地、耕地和林地等，代表的是陸地上空大氣水汽含量區(qū)；region2的地表類型主要是水體，代表了水體上空大氣水汽含量區(qū)；region3代表全省的大氣水汽含量分布狀況。

表1　感興趣區(qū)（ROI）水汽含量統(tǒng)計(jì)表（g/cm2）

為了便于下一步的結(jié)果分析，分別針對三個特征區(qū)的最大值、最小值和平均值進(jìn)行了統(tǒng)計(jì)，具體結(jié)果見表1。

從表1可以看出，水汽的最大值4.21 g/cm2出現(xiàn)在水體上空，對于陸地表面最大值1.37 g/cm2，結(jié)果符合自然條件下水汽的實(shí)際分布規(guī)律。region2的水汽平均值3.69 g/cm2同樣大于region1的平均值0.95 g/cm2，此結(jié)果也符合水體上空水汽總量大于陸地上空水汽總量的客觀規(guī)律。對于平均值而言，region3的值在region1和region2之間，同樣也符合實(shí)際情況。

4　結(jié)論

本文利用兩通道比值法和三通道比值法對黑龍江省MODIS影像數(shù)據(jù)進(jìn)行水汽含量的反演，從目前的研究結(jié)果來看，反演結(jié)果與實(shí)際探空資料均有極顯著的相關(guān)性，且符合水汽實(shí)際分布規(guī)律，說明本文所述利用MODIS影像數(shù)據(jù)反演大氣水汽的方法是可行的。兩通道比值法反演結(jié)果的相關(guān)系數(shù)大于三通道比值法反演結(jié)果，因此認(rèn)為，利用MODIS影像數(shù)據(jù)對黑龍江大氣水汽含量反演時兩通道比值法較合適。

但是本文所采用的反演方法存在一定的局限性，對數(shù)據(jù)要求晴空大氣而且僅能反演云層以上的水汽含量，并且反演得到的水汽圖與實(shí)際的水汽圖還是存在一定的誤差，因?yàn)闆]有考慮到當(dāng)氣溶膠達(dá)到一定光學(xué)厚度時，對近紅外波段水汽反演所產(chǎn)生的影響。

［1］Li Z L，Su Z，et al.A new Approach for Retrieving Precipitable Water from ATSR2 Split-window ChannelData over Land Area［J］.International Journal of Remote Sensing，2003，24:5059-5117.

［2］Brian JS.Atmospheric physics enlightening water vapor［J］.Nature，2000，406；247-248.

［3］毛克彪，覃志豪等.針對MODIS數(shù)據(jù)的大氣水汽含量反演及31和32波段透過率計(jì)算［J］.國土資源遙感，2005（01）:26-29.

［4］Liu Y J，Yang Z D.The Principle and Arithmetic of MODIS［M］.Beijing:Science Press.2001.

［5］GAO B C，Yoram J Kaufman.The MODIS Near-IR Water Vapor Algorithm.Algorithm Theoretical Basis Document［Z］.ATBD-MYD-03，NASA Goddard Space Flight Center，1988.

［6］GAO B C，Alexander F H Geotz.Column Atmospheric Water Vapor and Vegetation Liquid Water Retrievals from Airborne Imaging Spectrometer Data［J］. Journal of Geophysical Research，1990，（95）:3549-3564.

［7］Kaufman Y J.Gao B C.Remote Sensing of Water Vapor in the Near IR from EOS/MODIS［J］.IEEE Transaction on Geo science and Remote Sensing. 1992.5（30）:871-884.

1002-252X（2016）03-0020-03

2016-6-1

翟墨（1988-），男，黑龍江省哈爾濱市人，東北農(nóng)業(yè)大學(xué)，本科生，助理工程師.