劉磊

摘 要：本文利用河南省本地的地基GPS網(wǎng)，以2018年9月20日的數(shù)據(jù)為研究范圍，并與MODIS衛(wèi)星大氣可降水量數(shù)據(jù)集MOD05相比，討論了GPS探空可降水量與MODIS衛(wèi)星數(shù)據(jù)的對(duì)比分析。結(jié)果表明：兩者相關(guān)系數(shù)較好，相關(guān)系數(shù)在0.8左右，相對(duì)誤差平均值在30%以?xún)?nèi)。說(shuō)明MOD05產(chǎn)品與GPS數(shù)據(jù)反演大氣可降水量在一定程度上保持一致性;MOD05產(chǎn)品與GPS數(shù)據(jù)反演大氣可降水量和地面相對(duì)濕度無(wú)明顯相關(guān)。此外，衛(wèi)星反演的可降水量不能說(shuō)明當(dāng)時(shí)地面相對(duì)濕度的狀態(tài)。

關(guān)鍵詞：GPS;大氣可降水量;MOD05;相對(duì)濕度

中圖分類(lèi)號(hào)：P412.2;P457.6 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1003-5168（2019）05-0144-05

Abstract： In this paper， based on the local ground-based GPS network in Henan Province， the comparative analysis of GPS sounding precipitable water and MODIS satellite data was discussed in comparison with MODIS Atmospheric precipitable water data set MOD05. The results showed that the correlation coefficients were good， the correlation coefficients were about 0.8， and the average relative error was less than 30%. It showed that the retrieval of Atmospheric precipitable water from MOD05 and GPS data was consistent to a certain extent， and there was no significant correlation between the retrieval of atmospheric precipitable water from MOD05 and GPS data and the relative humidity of the ground. In addition， the precipitation retrieved by satellite could not explain the relative humidity of the ground at that time.

Keywords： GPS;atmospheric precipitation;MOD05;relative humidity

大氣中的水汽在天氣演變和氣候變化中扮演著重要角色，其影響大氣輻射和太陽(yáng)輻射，從而影響氣候。目前，水汽的三維空間分布主要依靠國(guó)家業(yè)務(wù)探空網(wǎng)來(lái)提取。探空站之間的間距均在200、300km以上，很難捕捉空間尺度小于100km的中尺度天氣系統(tǒng)。而在已有的替代方法中，雙通道微波輻射計(jì)等觀測(cè)手段，不但費(fèi)用昂貴，不能全天候觀測(cè)，而且需經(jīng)常標(biāo)定，難以大范圍密集設(shè)置站網(wǎng)和實(shí)現(xiàn)觀測(cè)業(yè)務(wù)化。此外，氣象衛(wèi)星資料直接反演的水汽精度又較低，誤差約為20%[1]。而利用稠密的地基GPS技術(shù)探測(cè)大氣水汽含量具有費(fèi)用低和可以全天候觀測(cè)等優(yōu)點(diǎn)，因此，其在氣象領(lǐng)域得到廣泛應(yīng)用。

在GPS氣象學(xué)方面，王小亞[2]等通過(guò)天頂干延遲和天頂濕延遲的分離及天頂濕延遲到可降水量的轉(zhuǎn)換，得到了精度為1～2mm的可降水量的計(jì)算結(jié)果。將這些結(jié)果與實(shí)測(cè)探空儀資料計(jì)算的結(jié)果相比較，兩者基本符合。試驗(yàn)的結(jié)果初步驗(yàn)證了地面GPS觀測(cè)為氣象服務(wù)的可行性和可靠性。施闖[3]等分析了國(guó)家高精度GPS網(wǎng)中的粗差觀測(cè)量的特點(diǎn)，提出了采用相關(guān)分析的粗差探測(cè)理論，對(duì)國(guó)家高精度GPS網(wǎng)的各子網(wǎng)和相鄰子網(wǎng)進(jìn)行粗差探測(cè)，并采用穩(wěn)健估計(jì)的平差方法，消除粗差對(duì)平差結(jié)果的影響，大大提高了國(guó)家高精度GPS網(wǎng)的整體平差精度。李杰[4]等詳細(xì)介紹應(yīng)用GAMIT軟件進(jìn)行高精度GPS定位解算。陳小雷[5]等詳細(xì)分析了地基GPS遙測(cè)大氣可降水量在天氣分析診斷中的應(yīng)用。

河南省GPS基準(zhǔn)站點(diǎn)網(wǎng)絡(luò)是河南省氣象局聯(lián)合省地礦廳和測(cè)繪局共同興建的，儀器由地礦和測(cè)繪部門(mén)購(gòu)買(mǎi)，維護(hù)由氣象部門(mén)負(fù)責(zé)，儀器安裝在氣象觀測(cè)站，數(shù)據(jù)為各單位共享。

本文將GPS數(shù)據(jù)反演得到的大氣可降水量與衛(wèi)星數(shù)據(jù)進(jìn)行對(duì)比，分析了解GPS得到的大氣可降水量與衛(wèi)星紅外通道得到的可降水量的相關(guān)性。

2 數(shù)據(jù)與數(shù)據(jù)處理方法

2.1 數(shù)據(jù)

2.1.1 GPS數(shù)據(jù)。河南省氣象局與河南省測(cè)繪局、河南省地礦廳等單位協(xié)作，共建立了38個(gè)GPS基準(zhǔn)站進(jìn)行資料實(shí)時(shí)共享，使陸上基準(zhǔn)站的密度小于100km，河南省的GPS站點(diǎn)分布如圖1所示，不同的站點(diǎn)符號(hào)代表來(lái)自不同的部門(mén)（其中☆為省測(cè)繪局，△為地礦廳）。

2.1.2 GAMIT軟件介紹。GAMIT是由麻省理工學(xué)院（MIT）、美國(guó)加利福尼亞斯克瑞布斯（SCRIPPS）海洋研究所（SIO）等研制的用于大地測(cè)量的GPS分析軟件，之后經(jīng)許多人不斷改進(jìn)而成為應(yīng)用面較為廣泛的高精度GPS數(shù)據(jù)分析和處理軟件，是目前國(guó)際上著名的GPS定位定軌軟件之一。GAMIT集GPS定位定軌于一體，可以解算衛(wèi)星軌道和站坐標(biāo)以及其他一些輔助參數(shù)，如鐘差、大氣延遲改正參數(shù)和整周模糊度參數(shù)等。

GAMIT軟件的Q文件和O文件等經(jīng)計(jì)算得到大氣總延遲ZTD，計(jì)算可降水量PWV經(jīng)編程處理氣象要素，分離出天頂濕延遲，最后計(jì)算出可降水量。

2.1.3 MOD05大氣可降水量產(chǎn)品。研究中使用的第二種遙感數(shù)據(jù)為美國(guó)宇航局（NASA）地球觀測(cè)計(jì)劃（EOS）發(fā)射的中分辨率成像光譜儀（MODIS）數(shù)據(jù)。其共有兩顆星，上午星TERRA于每日地方時(shí)上午10∶30和晚上22∶30過(guò)境，大氣可降水量數(shù)據(jù)集為MOD05L2;下午星AQUA于下午14∶30和凌晨02∶30每日地方時(shí)13∶30左右過(guò)境，大氣可降水量數(shù)據(jù)集為MYD05L2。每顆星每次均包括紅外和近紅外兩種數(shù)據(jù)集。由于近紅外在夜間無(wú)法使用，故每天均有6次數(shù)據(jù)可用，其中近紅外分辨率為1km，可用數(shù)據(jù)每天2次;紅外分辨率5km，可用數(shù)據(jù)每天4次。

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

2.2.1 MOD05數(shù)據(jù)處理。對(duì)于MOD05數(shù)據(jù)，可以使用ENVI軟件直接處理。通過(guò)對(duì)MOD05數(shù)據(jù)進(jìn)行輻射定標(biāo)、投影轉(zhuǎn)換、圖像校正、圖像拼接、去除異常值和無(wú)效值等步驟后，得到近紅外和紅外兩個(gè)波段的降水量數(shù)據(jù)集。如果數(shù)據(jù)量較大，可以使用專(zhuān)門(mén)處理MODIS大氣數(shù)據(jù)的軟件HEG，HEG軟件的可視化界面如圖2所示。軟件可以在MODIS網(wǎng)站的工具頁(yè)面下載。本文使用HEG軟件處理MOD05數(shù)據(jù)。

為了最大限度地利用MOD05大氣可降水量，減少空間不穩(wěn)定性，使異常數(shù)據(jù)得到平滑處理，故可將所有的數(shù)據(jù)加權(quán)平均，公式如下：

式中[，PWV]為過(guò)程平均紅外可降水量;[PWVnir]為近紅外大氣可降水量;[PWVir]為紅外大氣可降水量。根據(jù)式（1）可以得出河南省2018年9月20日上午10：30左右紅外大氣可降水量，其數(shù)值分布在0～38mm。MOD05大氣可降水量產(chǎn)品數(shù)據(jù)見(jiàn)圖3。

2.2.2 氣象數(shù)據(jù)與GPS探空可降水量。本研究的GPS數(shù)據(jù)分析處理按照四個(gè)步驟進(jìn)行。首先，數(shù)據(jù)準(zhǔn)備;其次，用GAMIT軟件處理各時(shí)段GPS的觀測(cè)數(shù)據(jù)，得到各個(gè)時(shí)段的解;再次，用GLOBK（卡爾曼濾波）進(jìn)行多時(shí)段綜合解算，以得到平差的結(jié)果及測(cè)站的平均坐標(biāo)和速度等參數(shù);最后，把平差后得到的測(cè)站坐標(biāo)帶回GAMIT，對(duì)待求測(cè)站進(jìn)行強(qiáng)約束，反演對(duì)流層的天頂延遲值。得到天頂延遲（ZTD）還需要同時(shí)段站點(diǎn)的氣象數(shù)據(jù)，包括氣壓與溫度，最后計(jì)算出數(shù)據(jù)中所用時(shí)段的值[6]。由于氣象數(shù)據(jù)使用的是逐小時(shí)，所有計(jì)算出來(lái)的PWD時(shí)間分辨率為小時(shí)。本文使用的2018年9月20日上午11：00的GPS數(shù)據(jù)反演的大氣可降水量數(shù)值。同時(shí)，從自動(dòng)站數(shù)據(jù)庫(kù)中提取2018年9月20日上午11：00的小時(shí)相對(duì)濕度。GPS數(shù)據(jù)反演的大氣可降水量數(shù)值與同時(shí)刻相對(duì)濕度見(jiàn)圖4。

3 結(jié)果與分析

GPS數(shù)據(jù)反演大氣可降水量與MOD05大氣可降水量產(chǎn)品相比線性相關(guān)較好，相關(guān)系數(shù)在0.8左右，如圖5（a）所示。在數(shù)據(jù)上GPS反演的數(shù)值普遍大于MODIS產(chǎn)品。主要原因是反演的原理不同：MODIS產(chǎn)品是通過(guò)紅外和近紅外波段的亮溫利用分裂窗算法得出;GPS把各個(gè)部分相互作用的延遲量引入解算模型中，考慮誤差來(lái)源并予以消除后對(duì)高精度的大氣延遲量與定位參數(shù)一起求解，然后通過(guò)雙頻技術(shù)訂正電離層延遲到毫米量級(jí)，通過(guò)地面氣壓觀測(cè)量訂正干空氣的延遲到毫米量級(jí)，最后得到毫米量級(jí)的濕空氣延遲，由此提取大氣總的水汽含量。紅外和近紅外受到氣溶膠、云的影響相對(duì)更大。

GPS數(shù)據(jù)反演大氣可降水量與相對(duì)濕度線性一般[見(jiàn)圖5（b）]。相對(duì)濕度觀測(cè)的是地面數(shù)據(jù)，對(duì)整層可降水量影響較小。豫西部分觀察站海報(bào)較高，對(duì)相對(duì)濕度的數(shù)值也有影響。MODIS產(chǎn)品對(duì)相對(duì)濕度的線性相關(guān)較差，相關(guān)系數(shù)在0.4左右，如圖5（c）所示。

4 結(jié)論

本文利用MOD05產(chǎn)品與GPS數(shù)據(jù)反演大氣可降水量數(shù)據(jù)進(jìn)行對(duì)比，結(jié)果表明：兩者相關(guān)系數(shù)較好，相關(guān)系數(shù)在0.8左右，相對(duì)誤差平均值在30%以?xún)?nèi)?？梢?jiàn)，MOD05產(chǎn)品與GPS數(shù)據(jù)反演大氣可降水量在一定程度上保持一致性。MOD05產(chǎn)品在晴空下反演的數(shù)值較GPS數(shù)據(jù)反演大氣可降水量明顯小。由于衛(wèi)星能在各種氣象條件下監(jiān)測(cè)地球大氣，其長(zhǎng)期穩(wěn)定性、分辨率、覆蓋范圍和精度都是前所未有的。其為地球表面上難于進(jìn)行定期測(cè)量的地區(qū)提供了一個(gè)全新的測(cè)定大氣剖面的方法，這對(duì)廣闊的海洋、沙漠和近極地等荒漠地區(qū)上空的大氣可降水量研究具有特定的意義。

MOD05產(chǎn)品和GPS數(shù)據(jù)反演大氣可降水量與地面相對(duì)濕度無(wú)明顯相關(guān)。衛(wèi)星反演的可降水量不能說(shuō)明當(dāng)時(shí)地面相對(duì)濕度的狀態(tài)。

參考文獻(xiàn)：

[1]劉敏，郭鵬，葉其欣，等.上海地區(qū)地基GPS水汽三維層析技術(shù)和初步應(yīng)用[J].天文學(xué)報(bào)，2010（3）：299-308.

[2]王小亞，朱文耀，嚴(yán)豪健，等.地面GPS探測(cè)大氣可降水量的初步結(jié)果[J].大氣科學(xué)，1999（5）：605-612.

[3]施闖，劉經(jīng)南.國(guó)家高精度GPS網(wǎng)整體平差中的粗差分析[J].武漢測(cè)繪科技大學(xué)學(xué)報(bào)，1999（2）：107-111.

[4]李杰，王曉強(qiáng).應(yīng)用GAMIT軟件進(jìn)行高精度GPS定位解算[J].內(nèi)陸地震，2001（4）：355-359.

[5]陳小雷，馬翠平.GPS遙測(cè)大氣可降水量原理與方法[J].河北氣象，2005（4）：37-38.

[6]李成才，毛節(jié)泰.GPS地基遙感大氣水汽總量分析[J].應(yīng)用氣象學(xué)報(bào)，1998（4）：470-477.