毛 健，朱長青，蘇 笛

（1.南京師范大學(xué)虛擬地理環(huán)境教育部重點實驗室，江蘇南京210046；2.天津師范大學(xué)城市與環(huán)境科學(xué)學(xué)院，天津300387）

全球?qū)α鲗犹祉斞舆t特征研究

毛 健1，2，朱長青1，蘇 笛2

（1.南京師范大學(xué)虛擬地理環(huán)境教育部重點實驗室，江蘇南京210046；2.天津師范大學(xué)城市與環(huán)境科學(xué)學(xué)院，天津300387）

基于IGS提供的對流層天頂延遲數(shù)據(jù)，分析單站對流層天頂延遲特征，據(jù)此構(gòu)建單站對流層天頂延遲模型，并對其進(jìn)行擬合得到相關(guān)模型參數(shù)；比較模型參數(shù)與地理位置的關(guān)系，得出對流層天頂延遲與測量時間及測站地理位置的相關(guān)性，并分析其在全球范圍內(nèi)的連續(xù)分布情況：全球?qū)α鲗犹祉斞舆t變化具有周年性，且與測站位置和高程密切相關(guān)；其較高值多集中在赤道附近，并以赤道為中心向兩側(cè)遞減。該研究為利用全球IGS站的對流層天頂延遲數(shù)據(jù)建立無需氣象參數(shù)的全球?qū)α鲗犹祉斞舆t模型提供了可靠依據(jù)。

對流層天頂延遲；特征研究；GPS氣象學(xué)

0 引言

對流層改正模型能有效降低對流層延遲引起的電磁波傳播誤差，在GPS導(dǎo)航定位中得到了廣泛應(yīng)用，如何建立精度更高的對流層改正模型，是目前GPS領(lǐng)域研究的熱點。對流層延遲改正模型按其計算方法分為兩類：1）需引入氣象參數(shù)的對流層延遲改正模型，主要有Hopfield、Saastmoinen、EGNOS等［1－3］，引入的氣象參數(shù)一般分為實測氣象參數(shù)與標(biāo)準(zhǔn)氣象參數(shù)。由實測氣象參數(shù)計算的對流層延遲改正精度較高，但其量測復(fù)雜，不宜獲得。而標(biāo)準(zhǔn)氣象參數(shù)只需根據(jù)測站點地理位置即可從相關(guān)氣象模型中計算得到，因此，在實際衛(wèi)星導(dǎo)航定位中應(yīng)用較多；但由于計算標(biāo)準(zhǔn)氣象參數(shù)時存在氣象模型計算誤差，從而降低了對流層延遲改正的精度。2）無需氣象參數(shù)的對流層延遲改正模型，如戴吾蛟等建立了無需氣象參數(shù)的香港地區(qū)精密天頂對流層模型［4］，但其具有一定的區(qū)域局限性，無法在全球范圍推廣。

可見，建立一種適應(yīng)全球范圍且與氣象參數(shù)無關(guān)的對流層改正模型，是克服上述兩類模型不足的最佳途徑，而這種模型的建立須了解全球?qū)α鲗犹祉斞舆t（Zenith Tropospheric Delay，ZTD）分布特征。因此，本文基于IGS提供的高精度對流層天頂延遲數(shù)據(jù)，分析全球?qū)α鲗犹祉斞舆t的周年變化及隨地理位置變化的特征，為全球?qū)α鲗有履Ｐ偷慕⑻峁﹨⒖家罁?jù)，對GPS氣象學(xué)研究和觀測網(wǎng)的布設(shè)及其他相關(guān)對流層的研究也具有一定的指導(dǎo)意義。

1 單站對流層天頂延遲改正模型擬合

本文以全球IGS跟蹤站點的對流層天頂延遲數(shù)據(jù)為研究對象，選擇了具有3～7 a的299個IGS站（圖1，見封3）的天頂延遲時間序列。數(shù)據(jù)采樣間隔為5 min，精度為4 mm［5］。首先對單個測點的對流層天頂延遲特征進(jìn)行研究，據(jù)此構(gòu)建單站對流層天頂延遲模型并進(jìn)行擬合，得到各測點的單站模型參數(shù)。

1.1 單站對流層天頂延遲特征

本文以時間序列為橫軸，ZTD值為縱軸，將各IGS站的多年ZTD值繪制成圖。圖2為4個IGS站的6年ZTD值變化情況，從中可以看出天頂延遲的變化規(guī)律：1）中性大氣天頂延遲具有強(qiáng)而平穩(wěn)的周年特性；2）對流層天頂延遲周年變化的相位有很好的一致性，北半球?qū)α鲗犹祉斞舆t均在7月出現(xiàn)最大值、1月達(dá)到最小值，而南半球?qū)α鲗犹祉斞舆t均在1月出現(xiàn)最大值、7月達(dá)到最小值。

1.2 單站對流層天頂延遲模型擬合

根據(jù)單站對流層天頂延遲的周年變化規(guī)律，采用周期函數(shù)與線性函數(shù)構(gòu)建單站天頂大氣延遲的改正模型［6］：

式中：doy為年積日，a為對流層周年變化振幅，b為相位（年積日），c為天頂對流層延遲年均值（m）。

利用299個IGS站的天頂對流層延遲數(shù)據(jù)，采用最小二乘法分別對單站模型進(jìn)行擬合，得到299組a、b、c參數(shù)值，作為分析全球范圍對流層天頂延遲變化特征的主要數(shù)據(jù)。

圖2 ALIC、AUCK、ALBH、BJFS站6年對流層天頂延遲的變化Fig.2 The ZTD variation of ALIC，AUCK，ALBH and BJFS in 6 years

2 全球?qū)α鲗犹祉斞舆t變化特征分析

2.1 對流層天頂延遲隨空間位置的變化特征

利用擬合得到的299組a、b、c參數(shù)值，從經(jīng)度、緯度及高程方面計算相關(guān)系數(shù)和分析圖形分布特征（表1、圖3），了解全球?qū)α鲗犹祉斞舆t變化規(guī)律。

表1 空間位置與a、b、c值的相關(guān)系數(shù)Table 1 Correlation coefficient of geographical position and a，b，c

圖3 空間位置與a，b，c的關(guān)系Fig.3 Relationship between geographical position and a，b，c

從表1及圖3可以看出：經(jīng)度與a、b、c的相關(guān)系數(shù)很低，最大值僅為0.0737；高程與a、b的相關(guān)性也較低，但與c值高度相關(guān)，其相關(guān)系數(shù)達(dá)0.9313；緯度與a具有較高的相關(guān)性，與b、c相關(guān)性較低，a在緯度方向以（0，0）點為中心呈對稱分布（圖3b），c值在緯度方向以縱軸為中心呈對稱分布（圖3h），但由于受高程影響，造成其與緯度的相關(guān)系數(shù)較低?？紤]到高程與c的高度相關(guān)性，且c值與高程呈一定的線性關(guān)系（圖3i），其隨高度的增加而減少，因此利用線性函數(shù)進(jìn)行擬合得出如下關(guān)系式：

式中：h為大地高程，單位為m。

根據(jù)式（2），將各站的高程均歸算到0，得到經(jīng)高程改正后的c值與緯度的關(guān)系（圖4），其具有較高的相關(guān)性。

圖4 經(jīng)高程改正后的緯度與c的關(guān)系Fig.4 Relationship between latitude and c after elevation correction

通過以上分析得出對流層天頂延遲分布特征：1）天頂對流層延遲與經(jīng)度方向并無太大關(guān)系；2）天頂對流層延遲周年變化相位（b）與測站地理位置無關(guān)，且比較集中分布在一定的區(qū)域內(nèi)，可能為一常數(shù)；3）對流層天頂延遲變化特征主要取決于測站的緯度和高程。

2.2 全球?qū)α鲗犹祉斞舆t連續(xù)分布情況

由于采集的是299個離散分布的IGS站ZTD數(shù)據(jù)，無法了解對流層天頂延遲在全球范圍內(nèi)的連續(xù)分布情況。因此，首先對各站的c值（ZTD年均值）進(jìn)行格網(wǎng)化，然后使用球諧函數(shù)展開，得到全球范圍內(nèi)連續(xù)分布的ZTD平均值。具體步驟如下：

（1）格網(wǎng)化。采用surfer8.0軟件中的克立金法對其中275個站（剩余24個站用于檢驗）的c值進(jìn)行格網(wǎng)化，格網(wǎng)寬度5×5，經(jīng)度0°～360°，緯度－90°～90°。

（2）球諧函數(shù)展開及精度評估。地球物理中常用的球諧分析是將地球表面觀測的某個物理量展成球諧函數(shù)的級數(shù)：

式中：Rnm（θ，λ）＝Pnm（cosθ）cos mλ，Snm（θ，λ）＝Pnm（cosθ）· sin mλ，anm、bnm為常數(shù)，Pnm（cosθ）為諦合勒讓德函數(shù)，可由遞推公式給出［7］，θ為地面一點和地心的連線與地軸的夾角，λ為地球的經(jīng)度。

球諧函數(shù)的展開就是求anm、bnm這兩個常數(shù)，而由正交關(guān)系可得：

其中：dσ＝sinθ·dθ·dλ，因此通過積分即可確定anm和bnm的值。

本文使用Fortran語言，將以上公式編程，利用格網(wǎng)化后得到的數(shù)據(jù)，以15階次展開；對展開后獲得的模型與未參與格網(wǎng)化的24個站點進(jìn)行比較，得到的擴(kuò)展不確定度［8］為64.7 mm，即模型精度為±64.7 mm。為了更直觀地了解全球?qū)α鲗犹祉斞舆t分布情況，根據(jù)球諧函數(shù)展開模型，以0.25× 0.25的格網(wǎng)值，繪制出全球ZTD年均值c的分布圖（圖5，見封3）。

圖1中黃色區(qū)域為全球高程較高地區(qū)，將其與圖5相比較可以看出：圖5中區(qū)域1、5為低海拔區(qū)，但其c值小于2 m，而高海拔區(qū)域4的c值大于2.3 m，這都不符合對流層天頂延遲與高程的相關(guān)性規(guī)律，主要是在上述區(qū)域IGS跟蹤站（圖1）稀少所致；而在具有一定IGS跟蹤站數(shù)量的高海拔區(qū)域2、3、6，其ZTD值符合隨高程增加而減少的規(guī)律；并且對流層天頂延遲較高值多集中在赤道附近，且以赤道為中心向兩側(cè)遞減。

3 結(jié)論

GPS反演水汽技術(shù)是近年來衛(wèi)星測量和氣象領(lǐng)域提出的一種探測對流層水汽含量的新方法，能夠精確測定對流層天頂延遲量，為對流層的研究、特別是為構(gòu)建無需氣象參數(shù)的全球?qū)α鲗犹祉斞舆t改正模型提供了新方法。本文利用經(jīng)GPS反演所得的IGS站對流層天頂延遲數(shù)據(jù)，分析了對流層天頂延遲的分布特征，得出以下結(jié)論：1）利用IGS對流層數(shù)據(jù)構(gòu)建無需氣象參數(shù)的全球?qū)α鲗友舆t改正模型切實可行，而且該模型與測量日期和測站的緯度、高程相關(guān)，但如何處理無IGS站分布區(qū)域還需進(jìn)一步研究；2）全球?qū)α鲗犹祉斞舆t變化具有周年性，且與測站的緯度和高程密切相關(guān)，其較高值多集中在赤道附近，且以赤道為中心向兩側(cè)遞減；3）在進(jìn)行GPS氣象學(xué)研究、特別是布設(shè)GPS氣象網(wǎng)時，應(yīng)考慮測站緯度和高程的影響。

［1］ HOPFIELD H S，ANGUS－LEPPAN P V.Improvements in the tropospheric refraction correction for range measurement and discussion［J］.Philosophical Transactions of the Royal Society A：Mathematical，Physical and Engineering Sciences［1364－503X］，1980，294（1410）：341－352.

［2］ SAASTAMOINEN J H.Atmospheric correction for the troposphere and stratosphere in radio ranging of satellites［A］.The Use of Artificial Satellites for Geodesy［C］.1971.247－251.

［3］ PENNA N，DODSON A，CHEN W.Assessment of EGNOStropospheric correction model［J］.Nigel Penna and Others，1972，54：37－55.

［4］ 戴吾蛟，陳招華，匡翠林，等.區(qū)域精密對流層延遲建模［J］.武漢大學(xué)學(xué)報（信息科學(xué)版），2011，36（4）：392－396.

［5］ GODDARD SPACE FLIGHT CENTER.Troposphere Zenith Path Delay Data［DB／OL］.2011.ftp：／／cddis.gsfc.nasa.gov／ gps／products／troposphere／zpd.2011－10－31.

［6］ 曲偉菁.中國地區(qū)GPS中性大氣天頂延遲研究及應(yīng)用［D］.中國科學(xué)院研究生院，2007.

［7］ 傅承義，陳運(yùn)泰，祁貴仲.地球物理學(xué)基礎(chǔ)［M］.北京：科學(xué)出版社，1985.447.

［8］ 倪育才.實用測量不確定度評定［M］.北京：中國計量出版社，2004.

Abstract：It is necessary to do research on the characteristics of the global zenith tropospheric delay for building a more accuracy global zenith tropospheric delay model without meteorological parameters.Based on the precision troposphere data from IGS，this paper first analyzes the characteristics of a single station zenith tropospheric delay，and sets up the single－station zenith tropospheric delay model which is fitted to get the model parameters.Then the relationship between the model parameters and geographical location is analyzed.The analysis shows that the zenith tropospheric delay is closely related to measurement time and station location.Those results provide a reliable basis for building a more accuracy global tropospheric zenith delay model without meteorological parameters，and also gives a direction for GPS meteorology，especially the establishment of the GPS meteorological network.

Key words：Zenith Tropospheric Delay；research on the characteristics；GPS meteorology

Research on the Characteristics of the Global Zenith Tropospheric Delay

MAO Jian1，2，ZHU Chang－qing1，SU Di2
（1.Key Laboratory of Virtual Geographic Environment，Nanjing Normal University，Ministry of Education，Nanjing 210046；2.College of Urban and Environmental Science，Tianjin Normal University，Tianjin 300387，China）

P228

A

1672－0504（2012）04－0107－04

2012－03－15；

2012－05－24

國家自然科學(xué)基金資助項目（41071245）；江蘇高校優(yōu)勢學(xué)科建設(shè)工程資助項目

毛?。?983－），男，博士研究生，助教，從事GPS氣象學(xué)、空間數(shù)據(jù)安全等方面的研究。E－mail：mao＿jian1018＠163.com