相對濕度對天頂總延遲量的影響

張　偉，李星光

(1.合肥市測繪設計研究院，合肥　230071；2.上海司南衛星導航股份有限公司，上海　201103)

?

相對濕度對天頂總延遲量的影響

張偉1，李星光2*

(1.合肥市測繪設計研究院，合肥230071；2.上海司南衛星導航股份有限公司，上海201103)

摘要：針對相對濕度參數對GPS水汽反演精度的影響的存在性，以香港衛星參考站網數據為例，提出了利用GAMIT軟件實驗分析相對濕度對天頂總延遲量的方法。采用單一控制變量法，分別分析了相對濕度參數設置及其本身對天頂總延遲量解算的影響。結果表明：在相對濕度存在且無誤差時，相對濕度參數設置對ZTD解算影響很小。同時，20%的相對濕度誤差對ZTD解算影響也不大。而當相對濕度量不存在時，默認的相對濕度設置(50%)引起的ZTD誤差較大，在一定程度上會影響GPS水汽反演精度。因此，在有相對濕度量數據時，都應該被使用。

關鍵詞：GAMIT軟件；相對濕度；天頂總延遲；GPS水汽

0引言

美國的全球定位系統(global positioning system，GPS)作為全球衛星導航系統(global navigation satellite system，GNSS)最先建成的導航定位系統，已經在各行各業中獲得了廣泛的應用，其中地基GPS氣象學是近些年研究的熱點問題，利用GPS反演大氣可降水量來分析天氣特征變化的優勢逐漸顯現，尤其在是暴雨(雪)、冰雹、臺風等惡劣天氣狀況下。GPS水汽探測手段主要是在利用高精度GPS處理軟件實現GPS原始數據的處理并在提取高精度的對流層天頂總延遲量(zenith total delay，ZTD)的基礎之上，通過剝離ZTD中的干延遲分量獲取精確可用的濕延遲分量，最終通過轉化關系實現空間水汽含量的提取[1]。因此，解算并提取高精度的ZTD成為首要解決的問題。

研究表明實測的地面氣象數據作為建立對流層延遲模型的主要物理量參與解算能夠有效的提高解算的精度[2]。同時，國內外大量學者研究證實，天頂干延遲(zenith hydrostatic delay，ZHD)是與露點溫度、水汽壓密切相關[3-4]，且地面溫度露點差和水汽壓等反映近地層空氣濕度量被證實與降水(雪)出現的時間和強度存在一定關系[5]。文獻[6]利用最小二乘曲線擬合GPS觀測結果、大氣可降水量(precipitable water vapor,PWV)與相對濕度時間序列能幫助短期暴雨監測預報。而對于相對濕度在ZTD解算中的影響分析還很少見，基于此，本文利用GAMIT軟件解算GPS原始數據處理研究分析相對濕度量對ZTD解算的影響程度，以圖得到一些有益結論。

1數據來源

香港地區衛星參考站網(SatRef)數據只需要注冊申請并通過香港大地測量組許可便能免費下載使用。SatRef由12個連續運行的參考站組成，分布相對較密集(圖1)，且每個GPS站點都配置了Meteo Met*/Met*氣象儀對測站表面的氣象數據(氣壓、溫度、相對濕度)進行采集，有利于進行GPS遙感水汽計算分析。同時，考慮到相對濕度易受天氣變化的影響，選取了2013年第264～266天遭遇臺風“天兔”影響時間段內的數據進行分析。

圖1　香港衛星參考站網分布

2GAMIT軟件參數設置及約束條件

利用GAMIT軟件對香港衛星定位參考站網觀測數據進行處理，其具體固定設置參見表1。同時為了消除區域GPS網ZTD解算的相關性[7]，獲得更精確的PWV估值，同時顧及國際GNSS服務(international GNSS service，IGS)站數量及空間分布對GPS網解算精度的影響[8]，引入LHAZ、SHAO、TWTF、WUHN 4個IGS國際跟蹤站點組成長基線參與解算。

GPS/PWV精度與定位精度密切相關，適合的參數設置能在保證基線解算精度的同時獲取高精度的氣象參數，具體設置綜合參考文獻[8-13]。

GAMIT軟件中相對濕度控制項在tables文件夾下的sestbl.表文件中，采用GPS測站接收機自主交換格式(receiver independent exchange format，RINEX)氣象文件時，其具體設置如下[12]：

Met obs source=RNX UFL GPT 50

表1　GAMIT主要固定解算設置及約束方案

表示為先檢查RINEX氣象文件，然后檢查u文件中的測站，最后返回GPT模型值，最后一個數值為假設的相對濕度值，50為默認設置(%)。當最后一個數值設置為-1(負值)時，表示僅使用RINEX氣象文件中的相對濕度值。若RINEX氣象文件中的相對濕度量不符合GAMIT軟件解算時規范要求(HR∈(0，100))或部分缺失，軟件會自動將相對濕度值更新成設置的值或離該濕度量最近的上一個相對濕度量。通常對于沿海地區相對濕度值較大，且出現100概率很高，為了排除這種情況可能產生的影響，可先對相對濕度量進行檢查并將值為100的更改為99.99。

3實驗分析

以下進行兩組實驗分別驗證相對濕度參數設置和相對濕度本身對天頂總延遲量解算的影響。限于圖幅要求，在同一解算后只選取了HKSC站點用于分析。

3.1相對濕度參數設置對ZTD解算的影響

無GPS站點RINEX氣象數據參與解算，只能比較不同假定相對濕度值設置對ZTD解算的影響。表2列出了不同相對濕度參數設置中的歸一化均方根(normalized root mean square，NRMS)可以看出都滿足基線解算的精度評定指標(<0.3[13])，且不同相對濕度值設定對基線解算精度影響很小。圖2中rhX分別為對應的相對濕度量為X的ZTD值(下同)。可以看出ZTD互差很小，在mm級范圍內，故而對PWV反演精度影響更小，幾乎可以忽略不計。

表2　不同相對濕度參數設置中的NRMS

圖2　ZTD差值時序分布

有GPS站點RINEX氣象文件參與解算時，不同相對濕度值設定對基線解算精度影響很小(表3)。這里假定相對濕度值設置為-1時解算得到的ZTD值為真值，記為rhz。圖3為ZTD解算誤差分布，可以看出兩種相對濕度參數設置引起的ZTD誤差不大，總體在2 mm以內，且誤差曲線除個別時刻外幾乎一致，表明在RINEX氣象文件符合規范時，相對濕度參數設置不能引起ZTD的較大偏差。同時也能明顯看出第266天(48～72 h)ZTD誤差遠小于前2 d，究其原因可能是，第264天、265天前后臺風“天兔”襲擾香港帶來的天氣變化引起的濕度變化對ZTD解算的影響。由此可以推斷ZTD與天氣變化引起的濕度變化是存在某種聯系的。

表3　不同相對濕度參數設置中的NRMS

圖3　ZTD差值時序分布

3.2相對濕度對ZTD解算的影響

在實際解算過程中往往會出現相對濕度量缺失或存在誤差的情況，為了比較這兩種情況下相對濕度對ZTD解算的影響，這里假定原RINEX氣象為真值，同時采用單一控制變量法，人為將HKSC站點的相對濕度量引入20%的誤差，并將獲取的ZTD值記為rh20；將相對濕度量缺失并在解算中采用默認的假設相對濕度值(50%)的情況下得到的ZTD值記為rh50。圖4給出了假定情況下的ZTD誤差分布，可以看出，當相對濕度值全部采用默認值時，ZTD誤差相對較大，必然會對PWV解算精度產生影響。而相對濕度量存在20%的誤差時引起的ZTD誤差總體偏小(ZTD誤差均值0.065 mm，均方根偏差1.07 mm)，因此對PWV精度產生的影響不大，這對于獲取高精度PWV來說是有利的。圖5給出了相對濕度誤差與ZTD誤差的分析，不難發現，相對濕度采用默認值時，ZTD誤差與相對濕度誤差間的對應關系并不明顯，即伴隨著相對濕度增加并未出現ZTD誤差的增大或減小。

圖4　ZTD誤差時序分布

圖5　ZTD誤差與相對濕度誤差對應關系

盡管在多數情況下，相對濕度量設置并不會對ZTD解算精度產生大的影響，但在無RINEX氣象數據參與解算時，氣壓、溫度、相對濕度對ZTD產生的綜合誤差影響是一直存在的，如圖6給出的有實測氣象文件參與解算與無實測氣象文件參與解算的差值分布。因此，在有實測RINEX氣象數據時，是應該被使用的。當無RINEX氣象數據時，可采用氣溫、氣壓插值手段進行獲取，并按照RINEX氣象文件的格式進行編寫后使用。

圖6　HKSC站ZTD誤差分布

4結束語

利用GAMIT軟件解算香港衛星參考站網數據，分析相對濕度設置與相對濕度對天頂總延遲量解算的影響，得到的結論為

(1)在正常天氣狀況下，有無實測RINEX氣象文件參與解算時，相對濕度設置對ZTD解算產生的影響很小，在天氣狀況發生變化時(臺風、降水等)，相對濕度設置對ZTD解算的影響稍大，但總體上并不影響GPS水汽反演的精度。

(2)當相對濕度存在20%的誤差時，由此引起的ZTD誤差總體較小，而相對濕度缺失時采用默認相對濕度(50%)引起的ZTD誤差總體較大且不容忽視，可見相對濕度對ZTD解算的是存在一定影響的。因此，當RINEX文件中相對濕度值缺失時，合理的選擇相對濕度參數設置能在一定程度上提高ZTD解算精度。

參考文獻

[1]李國平，黃丁發，郭杰，等.地基GPS氣象學[M].北京：科學出版社，2010：2-5.

[2]丁曉光.對流層延遲改正在數據處理中的應用與研究[D].西安：長安大學，2009：8.

[3]谷曉平.GPS水汽反演及降雨預報方法研究[D].北京：中國農業大學，2004：47.

[4]劉旭春，張正祿，張鵬，等.GPS反演大氣綜合水汽含量的影響因素分析[J].測繪科學，2007，32(2)：21-23.

[5]李國翠，李國平，陳小雷.強降雪天氣中GPS可降水量與地面空氣濕度的綜合分析[J].高原氣象，2011，30(6)：1626-1632.

[6]王留朋，白征東，過靜珺，等.利用地基GPSPWV序列和相對濕度RH序列研究暴雨過程[J].測繪科學，2007，32(3)：142-143.

[7]DUAN Jingping，BEVISA M，FANG Peng，et al.GPS meteorology：direct estimation of the absolute value of precipitable water[J].Journal of Applied Meteorology，1996，35(6)：830-838.

[8]徐平，尹繼堯，吳培稚，等.北京市GPS網解算時IGS站的選取實驗[J].大地測量與地球動力學，2006，26(1)：49-54.

[9]高旺，高成發，潘樹國，等.基于廣播星歷的GAMIT基線解算方法及精度分析[J].測繪工程，2014，23(8)：54-57.

[10]黃祖登，梅連輝.高精度GPS數據處理中對流層參數估計研究[J].測繪地理信息，2013，38(5)：12-15.

[11]李星光，鄭南山.參數設置對高精度GPS數據解算的影響探討[J].測繪科學，2015，40(1)：33-37.

[12]HERRING T A，KING R W，FLOYD M A，et al.GAMIT referece manual：release 10.6[EB/OL].(2015-06-16)[2015-07-03].http：//www-gpsg.mit.edu/～simon/gtgk/GAMIT_Ref.pdf.

[13]王晶.基于GAMIT的高精度基線解算方法及其在汶川地震分析中的應用[D].武漢：武漢大學，2008：28.

Influence of Relative Humidity on Zenith Total Delay Solution

ZHANGWei1，LIXingguang2

(1.Hefei Surveying and Mapping Institute，Hefei 230061，China；2.ComNav Technology Ltd.，Shanghai 201103，China)

Abstract：For analysing the influence of relative humidity on the retrieval of precipitable water from GPS solution，GAMIT software was used for solving the Hong Kong satellite reference station network data，and analysing the influence of relative humidity and its different settings on zenith total delay solution.The results show that：Relative humidity parameter setting have little effect on ZTD solver as the relative humidity was present and no error.Also 20% relative humidity error have little effect on ZTD solution.When the relative humidity was not existed，the default setting of relative humidity (50%) would cause a large ZTD error，and affect the accuracy of GPS Water Vapor retrieval to a certain extent.So，if we have the relative humidity data，they should be used.

Key words：GAMIT；relative humidity；zenith total delay；GPS water vapor

中圖分類號：P228

文獻標識碼：A

文章編號：2095-4999(2016)-01-0113-04

通訊作者簡介：李星光(1988—)，安徽淮南人，碩士，目前主要研究GNSS遙感應用。

作者簡介：第一張偉(1990—)，安徽合肥人，碩士，主要從事GNSS數據處理研究。

收稿日期：2015-05-18

引文格式：張偉，李星光.相對濕度對天頂總延遲量的影響[J].導航定位學報，2016，4 (1)：113-116.(ZHANG Wei，LI Xingguang.Influence of Relative Humidity on Zenith Total Delay Solution[J].Journal of Navigation and Positioning，2016，4(1)：113-116.)DOI：10.16547/j.cnki.10-1096.20160122.