基于GAMIT系統(tǒng)的GPS/MET數(shù)據(jù)反演及其在閩南地區(qū)獲取PWV的應(yīng)用

杜言霞 ，梁 鶯 ，吳勇凱 ，陳州川 ，林添水

（1.泉州市氣象局，福建 泉州362000；2.海峽氣象開(kāi)放研究實(shí)驗(yàn)室，福建 廈門361012；3.福建省大氣探測(cè)技術(shù)保障中心，福建 福州350008；4.安溪縣氣象局，福建 安溪362400）

隨著全球?qū)Ш叫l(wèi)星和低軌衛(wèi)星技術(shù)的進(jìn)步，GPS氣象學(xué)也得到了迅速發(fā)展。GPS技術(shù)具有自動(dòng)、精確、快速、成本低、時(shí)空分辨率高和不易受天氣變化影響等優(yōu)勢(shì)，是一種極具潛力的探測(cè)手段[1]，因此得到很多國(guó)家和地區(qū)的關(guān)注。研究表明，利用地基GPS觀測(cè)數(shù)據(jù)和氣象數(shù)據(jù)反演出的大氣可降水量具有較高的精度，在暴雨預(yù)測(cè)、數(shù)值天氣預(yù)報(bào)和其它極端天氣預(yù)測(cè)中都起到了重要的作用。

GPS氣象學(xué)（GPS/MET）分為地基GPS氣象學(xué)和空基GPS氣象學(xué)。20世紀(jì)80年代后期美國(guó)的Askne和Nordius最早提出用地基GPS網(wǎng)探測(cè)地球大氣，他們推導(dǎo)了PWV和大氣濕延遲之間的關(guān)系[2]。Bevis和Businger進(jìn)一步分析和討論了這個(gè)轉(zhuǎn)換關(guān)系，從而提出地基GPS氣象學(xué)[3]。

GPS/MET作為一種大氣遙感探測(cè)新技術(shù)，國(guó)外發(fā)達(dá)國(guó)家在技術(shù)研究及業(yè)務(wù)應(yīng)用上，已經(jīng)取得顯著成果，世界上許多國(guó)家都建立了GPS觀測(cè)網(wǎng)。最早是美國(guó)在1996年就開(kāi)始組建覆蓋全國(guó)的GPS/MET網(wǎng)，德國(guó)在地基和空基GPS氣象學(xué)方面也做了大量的工作，德國(guó)氣象局和地學(xué)研究中心（CFZ&DWD）利用組建的GPS站網(wǎng)數(shù)據(jù)，每小時(shí)提供一次全國(guó)GPS/PWV分布圖。近年來(lái)，我國(guó)氣象部門也積極開(kāi)展GPS/MET的研究和業(yè)務(wù)化試驗(yàn)，取得了一定的成績(jī)，把建設(shè)地基GPS探測(cè)大氣和電離層參數(shù)的系統(tǒng)列入了現(xiàn)代化氣象探測(cè)系統(tǒng)的一項(xiàng)重要內(nèi)容，并制定了“中國(guó)氣象局導(dǎo)航衛(wèi)星接收站網(wǎng)建設(shè)規(guī)劃”。截至目前，我國(guó)已組建覆蓋全國(guó)的實(shí)時(shí)GPS/MET觀測(cè)網(wǎng)，提供1 h間隔的全國(guó)范圍的可降水量圖。

1 GAMIT軟件簡(jiǎn)介

GAMIT軟件系統(tǒng)是由美國(guó)麻省理工學(xué)院（MIT）開(kāi)發(fā)的，由于它具有開(kāi)源性以及版本更新速度快等諸多優(yōu)點(diǎn)，因此得到了較為廣泛的應(yīng)用。它主要由以下幾個(gè)程序模塊構(gòu)成：?jiǎn)尾钭詣?dòng)修復(fù)周跳模塊、雙差自動(dòng)修復(fù)周跳模塊、數(shù)據(jù)文件的格式轉(zhuǎn)換、軌道積分模塊、觀測(cè)值模型、人工交互式修復(fù)周跳模塊、參數(shù)估計(jì)、用于創(chuàng)建SOLVE所需的M文件和一些輔助程序。并且在數(shù)據(jù)編輯完成之后，可利用SOLVE模塊求解出相關(guān)的未知參數(shù)[4]。

GAMIT軟件的主要功能包括：估算衛(wèi)星狀態(tài)矢量、太陽(yáng)光壓參數(shù)、衛(wèi)星Y方向偏差、地球自轉(zhuǎn)參數(shù)、臺(tái)站坐標(biāo)和對(duì)流層天頂延遲等[5]。

GAMIT采用雙差做基本觀測(cè)量，可以消除站鐘和星鐘之間的主要誤差；軌道誤差對(duì)測(cè)站相對(duì)位置不敏感，有利于精密定位；對(duì)于臺(tái)站初始坐標(biāo)的精度要求也不高[6]。但是減少了有效觀測(cè)，各站有效觀測(cè)因其在網(wǎng)內(nèi)位置的變化而變化，且不能用于單點(diǎn)定位。

2 GAMIT數(shù)據(jù)處理過(guò)程

GAMIT軟件數(shù)據(jù)處理過(guò)程分為3個(gè)步驟：GPS數(shù)據(jù)準(zhǔn)備、GAMIT參數(shù)配置和GAMIT數(shù)據(jù)解算。

2.1 GPS數(shù)據(jù)準(zhǔn)備

數(shù)據(jù)處理過(guò)程主要用到3種數(shù)據(jù)文件，分別為觀測(cè)文件、星歷文件和導(dǎo)航文件，這些文件可以從IGS 網(wǎng)站（ftp：//garner/ucsd.edu/pub）下載，但下載后需轉(zhuǎn)換為RINEX格式的文件。若下載的數(shù)據(jù)文件在WINDOWS系統(tǒng)下進(jìn)行處理時(shí)，可通過(guò)dos2unix命令進(jìn)行格式轉(zhuǎn)換[7]。

在單天目錄里用sh_link_rinex命令進(jìn)行相應(yīng)文件的鏈接，也可以直接把原始觀測(cè)文件、星歷文件和導(dǎo)航文件拷貝到相應(yīng)的單天目錄。然后進(jìn)入tables目錄，建立與/gamit/tables目錄的連接（%links.table J2000 2001）。

GPS衛(wèi)星星歷的質(zhì)量對(duì)GPS衛(wèi)星定位的精度起著至關(guān)重要的作用，利用GAMIT軟件解算GPS數(shù)據(jù)，對(duì)于不同的研究采用不同的星歷。如果反演的大氣可降水量用于天氣預(yù)報(bào)，那么應(yīng)采用精密預(yù)報(bào)星歷進(jìn)行解算，如果用于長(zhǎng)期性的氣象及氣候?qū)W研究，則可采用各種精密星歷（表1）。

表1 星歷一覽表

2.2 軟件參數(shù)配置

為使GAMIT的內(nèi)置參數(shù)能與GPS數(shù)據(jù)匹配，需重新配置 sittbl.、sestbl.、station.info、session.info 等文件，配置完成后一起放入表目錄[8]。

（1）臺(tái)站坐標(biāo)文件。獲得測(cè)站概略坐標(biāo)的方法主要有：利用原始觀測(cè)RINEX文件獲得、利用命令gapr_to_l生成L文件、利用命令sh_rx2apr生成單測(cè)站的L文件，得到的L文件需進(jìn)行手工編輯。

在GAMIT運(yùn)行過(guò)程中，初始坐標(biāo)的誤差必須在10 m以內(nèi)，否則可能會(huì)導(dǎo)致最后迭代不收斂。初始坐標(biāo)的誤差在1 m以內(nèi)，才能保證最終的臺(tái)站坐標(biāo)精度好于1 mm。但有時(shí)因?yàn)榻邮諜C(jī)性能不穩(wěn)定，或者是周圍觀測(cè)條件不佳，會(huì)影響某些時(shí)段的觀測(cè)，導(dǎo)致單點(diǎn)定位精度不好或者定位時(shí)不收斂，那么需要重新選擇其他時(shí)段的RINEX觀測(cè)文件。為了獲得較精準(zhǔn)的初始坐標(biāo)，通?？梢远噙x取幾天的RINEX觀測(cè)文件，進(jìn)行單點(diǎn)定位，對(duì)多次定位結(jié)果取均值。

（2）臺(tái)站信息文件。station.info文件包含了所有需要解算的臺(tái)站信息，包括臺(tái)站的天線相位中心對(duì)于地球標(biāo)準(zhǔn)橢球面在X、Y、Z方向的偏差以及接收機(jī)型號(hào)、天線型號(hào)等信息。這些信息都可從每個(gè)臺(tái)站的觀測(cè)文件的頭文件中獲取，其中所謂天線相位中心在X、Y、Z方向的偏差，實(shí)際只有天線相位中心的高度，因?yàn)樘炀€相位中心在Y和X方向的偏差一般為零。可以參考/templates/station.info對(duì)文件按照規(guī)定的格式進(jìn)行手工編輯，并根據(jù)研究需要引入相應(yīng)測(cè)站。另外，工程文件的名字也在這里定義。

值得注意的是，有些臺(tái)站的信息因?yàn)楦鼡Q接收機(jī)或遷址等原因，會(huì)進(jìn)行信息更新，在station.info文件中同一個(gè)臺(tái)站會(huì)有幾行信息，必須按時(shí)間先后順序排列，將最新的信息放在最后一行，GAMIT默認(rèn)讀取最后一行信息。

（3）運(yùn)行控制文件。sestbl.文件是數(shù)據(jù)處理方案的核心控制文件，很多參數(shù)都需要在該文件里配置。具體要根據(jù)GPS網(wǎng)類型、所要解決的問(wèn)題等采取相應(yīng)的配置方案。主要對(duì)sestbl.文件中的幾項(xiàng)進(jìn)行設(shè)置：①choice of experiment。RELAX：定軌、定位、解ERP；BASELINE：僅僅定位；ORBIT：僅僅定軌。一般為解算水汽獲取大氣總延遲的時(shí)候，僅僅需要選擇BASELINE，不需要對(duì)衛(wèi)星進(jìn)行定軌，廣播星歷的精度已經(jīng)足夠。

②choice of observable。LC_HELP：用 LC 觀測(cè)解模糊度；LC_RANGE：用于LC_HELP，但更強(qiáng)調(diào)偽距的作用；LC_ONLY：用 LC觀測(cè)不解模糊度；L1_ONLY：僅使用L1解模糊度；L2_ONLY：僅使用L2 解模糊度；L1，L2_INDEPEND：分別使用 L1，L2解模糊度。

③type of analysis。QUICK：每一個(gè)缺口加周跳標(biāo)志，對(duì)每一個(gè)周跳標(biāo)志（在連續(xù)資料中或在缺口處）求解bias，這樣解得參數(shù)多了，解的精度也就低了，但對(duì)周跳不敏感；0-ITER：ARC（optional），MODEL，AUTCLN （optional），SOLVE，SCANDD；1-ITER：兩個(gè)0-ITER序列，但第一個(gè)是QUICK解；2-ITER：在1-ITER基礎(chǔ)上再加一序列，用于確定軌道；對(duì)于僅僅解算基線，獲取大氣延遲量，只需要選取0-ITER即可。

（4）臺(tái)站約束文件。通過(guò)設(shè)置sittbl.文件可以完成測(cè)站約束，另外大氣模型等參數(shù)也在這個(gè)文件里設(shè)置。首先要檢查是否包含所有需要解算的臺(tái)站的信息，設(shè)置的約束值的大小要依據(jù)L文件的精度，如果L文件中初始臺(tái)站坐標(biāo)比較準(zhǔn)確，約束值可以設(shè)置比較小；如果L文件誤差大，約束值就要比較大，一般為5~10 m，否則在最后解算時(shí)會(huì)導(dǎo)致迭代不收斂。

（5）時(shí)段信息文件?？梢允止ぞ庉媠ession.info文件，也可以用sh_makexp命令創(chuàng)建，采樣間隔和衛(wèi)星號(hào)等參數(shù)都在這個(gè)文件里設(shè)置。

2.3 數(shù)據(jù)處理

GAMIT數(shù)據(jù)處理過(guò)程主要分為3個(gè)階段：首先利用GAMIT軟件處理下載的觀測(cè)數(shù)據(jù)，然后利用卡爾曼濾波進(jìn)行綜合解算，最后求出對(duì)流層天頂延遲[9]。

用戶可以根據(jù)研究情況需要，對(duì)約束條件和參數(shù)進(jìn)行設(shè)置。衛(wèi)星截止高度角、數(shù)據(jù)采樣率、星歷等參數(shù)也可以根據(jù)研究需要進(jìn)行選擇。GPS數(shù)據(jù)處理流程如圖1所示。

首先，利用GAMIT軟件解算出載波相位的理論值，然后根據(jù)理論值和觀測(cè)值求出相位殘差，最終求出坐標(biāo)以及各個(gè)大氣參數(shù)等。

以外文期刊為例，嘗試構(gòu)建一個(gè)基于學(xué)科與用戶行為的外文期刊的資源導(dǎo)航，該模型如何搭建？需要從資源保障和績(jī)效評(píng)估兩方面構(gòu)建各種指標(biāo)體系。

ARC模塊的作用是得到星歷表T文件；MODEL模塊的作用是將理論值與觀測(cè)值的偏差存入C文件，并對(duì)其進(jìn)行編輯和估算；AUTCLN模塊由兩個(gè)模塊組成，分別為單差修復(fù)周跳模塊（SINCLN）和雙差修復(fù)周跳模塊（DBCLN），在GPS數(shù)據(jù)處理中，正確探測(cè)與修復(fù)周跳是一個(gè)非常重要的問(wèn)題[10]；CFMRG模塊的作用是生成觀測(cè)方程M文件；SOLVE模塊的作用是利用最小二乘方法解算未知參數(shù)。H文件是GLOBK和GAMIT的數(shù)據(jù)交換中介。

圖1 GPS數(shù)據(jù)處理流程圖

GAMIT求解天頂延遲量的數(shù)據(jù)處理過(guò)程主要分為7個(gè)步驟[11]，具體流程見(jiàn)圖2。

（1）建立測(cè)站先驗(yàn)坐標(biāo)文件L-file。值得注意的是，當(dāng)用于水汽反演研究時(shí)，要求初始坐標(biāo)必須為精確的坐標(biāo)。

（2）執(zhí)行makexp命令。這一步主要是給定數(shù)據(jù)處理的時(shí)間，包括開(kāi)始時(shí)刻、采樣頻率和計(jì)算歷元，查看所需的文件是否齊全，建立所有準(zhǔn)備文件的輸出及一些模塊的輸入文件。這一步結(jié)束以后，會(huì)生成一個(gè)session.info文件，其中包含了數(shù)據(jù)處理相關(guān)的時(shí)間信息和衛(wèi)星信息。

（3）執(zhí)行sh_sp3fit命令。這一步是使用偽觀測(cè)（從觀測(cè)文件所得的位置和速度）獲得精確的衛(wèi)星初始狀態(tài)，得到T文件和Y文件。

（4）通過(guò)執(zhí)行makej命令，得到J文件。

（5）執(zhí)行makex命令。通過(guò)執(zhí)行這條命令可以生成兩個(gè)文件，分別為X文件和K文件。X文件主要是對(duì)觀測(cè)文件分析后按GAMIT需要的格式重新寫(xiě)過(guò)，K文件是利用臺(tái)站的初始坐標(biāo)、衛(wèi)星廣播星歷和偽距信息計(jì)算出的接收機(jī)鐘差。

（7）執(zhí)行csh b*.bat。批處理工作由ARC（軌道積分）、MODEL（組成觀測(cè)方程）、AUTCLN（自動(dòng)修復(fù)周跳）和SOLVE（最小二乘法求解參數(shù)）等模塊組成。

圖2 GAMIT求解天頂延遲量的數(shù)據(jù)處理過(guò)程

3 實(shí)例分析

3.1 GPS站的選取

在局域GPS觀測(cè)網(wǎng)（小網(wǎng)）反演大氣可降水量時(shí)，如果想提高反演結(jié)果的精度，則應(yīng)該適當(dāng)?shù)奶砑右恍┹o助站，組成一個(gè)監(jiān)測(cè)網(wǎng)進(jìn)行解算，這樣就能減小監(jiān)測(cè)網(wǎng)內(nèi)水汽的相關(guān)性[12]。本文以閩南地區(qū)的廈門站（站號(hào)：59134，下同）為基準(zhǔn)站，另外選取了該地區(qū)的其他5個(gè)站點(diǎn)聯(lián)合進(jìn)行解算，分別為：莆田（58946）、平和（59125）、平潭（58944）、南靖（59124）、福清（58942）。

3.2 GAMIT軟件數(shù)據(jù)處理參數(shù)配置

本文采用GAMIT軟件的版本號(hào)為10.34版。選取6個(gè)GPS永久跟蹤站2017年5月BT的觀測(cè)數(shù)據(jù)、以及精密星歷和導(dǎo)航文件等數(shù)據(jù)。數(shù)據(jù)處理是在慣性坐標(biāo)系J2000中進(jìn)行的，數(shù)據(jù)處理觀測(cè)量為無(wú)電離層延遲LC觀測(cè)值的雙差組合觀測(cè)量[13]，處理時(shí)對(duì)廈門跟蹤站的坐標(biāo)進(jìn)行強(qiáng)約束。資料處理的參數(shù)約束和方案見(jiàn)表2。

表2 資料處理中的參數(shù)約束和方案

相關(guān)的文件配置好以后，就可以進(jìn)行后面的數(shù)據(jù)解算過(guò)程。主要的數(shù)據(jù)處理過(guò)程如下：

（1）利用廈門及其它GPS基準(zhǔn)站組成的觀測(cè)網(wǎng)來(lái)綜合解算廈門地區(qū)的對(duì)流層天頂總延遲（Zenith Total Delay，ZTD）；

（2）利用相對(duì)精度較高Sastamonien模型求出廈門地區(qū)的天頂干延遲量（Zenith Hydrostatic Delay，ZHD），由ZTD減去ZHD，得到對(duì)流層濕延遲（Zenith Wet Delay，ZWD）[14]；

（3）利用廈門地區(qū)2016年的實(shí)際氣象資料統(tǒng)計(jì)擬合出當(dāng)?shù)氐募訖?quán)平均溫度模型，從而得到加權(quán)平均溫度Tm和反演算法中的轉(zhuǎn)換因子Π[15]；

Π稱為無(wú)量綱比例因子，即水汽轉(zhuǎn)換因子。k1、k2為常數(shù)，ρw表示濕大氣的摩爾質(zhì)量，Tm為大氣加權(quán)平均溫度（K），Rv為水汽常數(shù)。

（4）濕延遲量乘上由加權(quán)平均溫度求出的轉(zhuǎn)換因子，即PWV=Π·ZWD，便得出廈門地區(qū)的大氣水汽含量。

3.3 GAMIT解算結(jié)果分析

判斷結(jié)果是否正確主要有2個(gè)判斷標(biāo)準(zhǔn)：（1）是否具有可進(jìn)行合理估計(jì)的足夠數(shù)據(jù)；（2）數(shù)據(jù)的噪聲水平是否與選用的模型匹配[16]。第一條標(biāo)準(zhǔn)的首要指標(biāo)是基線分量誤差的大小，如果比預(yù)先設(shè)定的站點(diǎn)坐標(biāo)和軌道參數(shù)的先驗(yàn)約束還要大，則會(huì)發(fā)現(xiàn)Q-文件或autcln.sum文件中大量數(shù)據(jù)被autcln丟棄。對(duì)于第二條標(biāo)準(zhǔn)，主要的指標(biāo)是解決nrms（標(biāo)準(zhǔn)化的均方根誤差），也就是chi-square x2的自由度。在Q文件中可以查看解算出的結(jié)果是否符合標(biāo)準(zhǔn)，本例的nrms為0.3，符合要求結(jié)果低于0.5的標(biāo)準(zhǔn)。6個(gè)測(cè)站的基線邊長(zhǎng)度及精度見(jiàn)表3。

表3 基線邊解算精度統(tǒng)計(jì)

本文以解算廈門2017年4月BT水汽結(jié)果為例。圖3為利用GAMIT軟件解算出的對(duì)流層天頂總延遲（Zenith Total Delay，ZTD）時(shí)間序列圖。圖4為利用薩氏模型計(jì)算出的對(duì)流層干延遲分量（Zenith Hydrostatic Delay，ZTD）時(shí)間序列圖。圖5為利用天頂總延遲與干延遲分量相減間接求出的對(duì)流層濕延遲（Zenith Wet Delay，ZTD）時(shí)間序列圖。

利用經(jīng)驗(yàn)?zāi)Ｐ陀?jì)算濕延遲需要很多的地面氣象參數(shù)，給很多實(shí)際應(yīng)用造成不便，而且會(huì)增加計(jì)算濕延遲的誤差。因此，地基GPS氣象學(xué)中通常是從天頂總延遲中減去干延遲分量從而得到天頂濕延遲，而不選用經(jīng)驗(yàn)?zāi)Ｐ陀?jì)算天頂濕延遲[17]。

3.4 廈門地區(qū)加權(quán)平均溫度的局地算式

圖3 廈門站ZTD的時(shí)間序列圖

圖4 廈門站ZHD的時(shí)間序列圖

圖5 廈門站ZWD的時(shí)間序列圖

大氣加權(quán)平均溫度Tm的精度制約著水汽轉(zhuǎn)換因子的精度，因此要想反演出高精度的水汽結(jié)果，就必須保證加權(quán)平均溫度的精度[18]。Tm的獲取途徑有4 個(gè)[19]：（1）常數(shù)法，取 Tm=281 K，但是這種方法的精度較低；（2）利用當(dāng)?shù)靥娇諗?shù)據(jù)采樣積分公式計(jì)算；（3）利用數(shù)值預(yù)報(bào)得到加權(quán)平均溫度的值；（4）利用當(dāng)?shù)貧庀髷?shù)據(jù)算出適合當(dāng)?shù)禺?dāng)時(shí)段的回歸經(jīng)驗(yàn)公式。目前，由于氣象探空數(shù)據(jù)的時(shí)空分辨率較低，而且高空的數(shù)據(jù)不易獲得，常數(shù)法的精度又太低，數(shù)值預(yù)報(bào)的應(yīng)用還不夠成熟，因此一般都采用回歸經(jīng)驗(yàn)公式的方法。本文選取廈門2016年的氣象數(shù)據(jù)資料，經(jīng)過(guò)統(tǒng)計(jì)分析，求出當(dāng)?shù)氐募訖?quán)平均溫度模型。圖6為廈門地區(qū)地面氣溫與加權(quán)平均溫度的散點(diǎn)圖及趨勢(shì)線。

圖6 廈門地區(qū)地面氣溫與加權(quán)平均溫度的散點(diǎn)圖及擬合線

從兩者地面氣溫和加權(quán)平均溫度的趨勢(shì)圖可看出，所有的散點(diǎn)基本上分布在一條直線附近，說(shuō)明加權(quán)平均溫度與地面氣溫之間是線性相關(guān)的，因此可擬合出Tm和Ts的線性關(guān)系：

通過(guò)對(duì)加權(quán)平均溫度和地面氣溫值進(jìn)行統(tǒng)計(jì)分析并擬合，求出廈地區(qū)加權(quán)平均溫度的線性回歸關(guān)系為Tm=106.7+0.583 1Ts，均方根誤差RMS=3.045 K。

3.5 水汽結(jié)果求解與分析

局地加權(quán)平均溫度求出之后就可以確定轉(zhuǎn)換系數(shù)Π，將天頂濕延遲ZWD乘上一個(gè)轉(zhuǎn)換系數(shù)Π可以轉(zhuǎn)換得到大氣可降水量PWV[20]。圖7為廈門站2017年4月BTGPS反演的PWV與由探空資料反演求得的PWV的序列比較圖，圖8為2種方法求得的PWV之差的時(shí)間序列圖。

從圖7可以看出，GPS反演的水汽含量較好的反應(yīng)出水汽的積累和釋放過(guò)程，水汽總量大致是以迅速增加、相對(duì)穩(wěn)定、迅速降低的規(guī)律呈周期變化；并通過(guò)與探空資料反演的PWV比較可以看出，兩種方法得到的大氣可降水量具有良好的相關(guān)性。

圖7 GPS與探空反演的PWV比較

圖8 GPS與探空反演的PWV之差

從圖8可以看出，GPS反演的PWV與由探空資料反演的PWV相比，兩者偏差最小為0.3 mm，最大達(dá)到8 mm，其偏差的均方差為4.3 mm，均在允許的誤差范圍之內(nèi)。從而說(shuō)明，利用GAMIT軟件可以很好的解算PWV。

值得一提的是，GPS反演的PWV與探空PWV的相關(guān)性受大氣水汽含量的影響，當(dāng)大氣水汽含量較大時(shí)，GPS反演的PWV與探空PWV的平均偏差會(huì)增大，兩者的相關(guān)系數(shù)也會(huì)變小，這可能與GPS水汽反演方法有關(guān)，而與GAMIT軟件求解的天頂延遲無(wú)關(guān)[21]。

4 結(jié)論

基于閩南地區(qū)廈門及其他5個(gè)GPS站點(diǎn)的觀測(cè)資料和從IGS下載的精密星歷數(shù)據(jù)，通過(guò)GAMIT軟件計(jì)算出對(duì)流層濕延遲，并利用廈門地區(qū)的實(shí)際氣象資料，通過(guò)回歸的方法得到廈門地區(qū)加權(quán)平均溫度局地算式，進(jìn)而求出了該地區(qū)的大氣可降水量，并將最終反演結(jié)果與通過(guò)探空方法計(jì)算出的大氣可降水量進(jìn)行對(duì)比分析。

（1）利用GAMIT軟件可以很好的解算PWV。

（2）大氣加權(quán)平均溫度進(jìn)行本地化訂正有利于提高GPS反演PWV的準(zhǔn)確度。

大氣水汽含量的多少對(duì)GPS反演PWV準(zhǔn)確度的影響是一個(gè)值得重視的研究問(wèn)題。