利用全國陸態(tài)網絡站點反演天頂對流層延遲分布特征

王東振，趙 斌，譚 凱，張彩紅

(中國地震局地震研究所 地震大地測量重點實驗室，湖北 武漢 430071)

利用全國陸態(tài)網絡站點反演天頂對流層延遲分布特征

王東振，趙 斌，譚 凱，張彩紅

(中國地震局地震研究所 地震大地測量重點實驗室，湖北 武漢 430071)

針對對流層延遲研究范圍小、時間跨度短的問題，文中利用GAMIT軟件獲得2015年全國270個態(tài)網絡站點的對流層延遲值，通過GAMIT中的sh_metutil命令調用Metutil程序，計算并提取天頂對流層總延遲量、干延遲分量、濕延遲分量、溫度及氣壓值。對干、濕延遲的年平均量進行分析，獲得全國對流層延遲分布特征；對影響干延遲量的因素進行分析，得到其主要影響因素是氣壓和高程因素；最后分析濕延遲與降水量的關系，得到兩者的分布特征基本一致。

GAMIT；陸態(tài)網絡；對流層延遲；分布特征

“十一五”期間國家重大科技基礎設施建設項目之一的“中國大陸構造環(huán)境監(jiān)測網絡”(以下簡稱“陸態(tài)網絡”)是“九五”國家重大科學工程“中國地殼運動觀測網絡”(以下簡稱“網絡工程”)的延續(xù)。截止2011年已完成260個GPS連續(xù)觀測站的建設，利用這些站點可獲得大量高質量的觀測數(shù)據(jù)。利用這些觀測數(shù)據(jù)可以求得連續(xù)的天頂對流層延遲(ZTD)，同時利用ZTD可以直接反映當?shù)亟邓?，這對用于氣象預報、氣候監(jiān)測有著重要的意義[1-2]。國內外學者利用大量GPS觀測資料來獲得實時或近實時的對流層延遲和大氣水汽含量做了大量工作。曲建光等人提出無需其他轉換,利用對流層延遲可直接推算水汽含量的方法，且結果與水汽計量監(jiān)測結果相當[2]；唐超華等人利用天津地區(qū)CORS站數(shù)據(jù)進行了小區(qū)域天頂對流層延遲年周期變化研究[3]；而對于全國大范圍、長時間的天頂對流層延遲解算研究甚少。本文將利用全國陸態(tài)網絡站點2015年的觀測數(shù)據(jù)，解算得到各個站點天頂對流層總延遲、干延遲以及濕延遲量，分析全國不同地區(qū)對流層干延遲與氣壓、溫度及高程之間的關系，最后對比分析其濕延遲的反演結果與實際降水量的關系。

1 對流層延遲解算理論

由于GPS信號穿過對流層時，傳播速度會發(fā)生改變，傳播路徑也會發(fā)生彎曲，這種影響就叫做對流層延遲效應。所以任意角度ε方向的對流層延遲可表示為[4-5]

ΔL(ε)=ZHDmfh(h)+ZWDmfw(h).

(1)

式中:ZHD和ZWD分別為干、濕分量；mfh(h)和mfw(h)分別為干、濕映射函數(shù)。

鑒于上式中的干、濕映射函數(shù)是基于大氣球對稱假設而建立的，這與實際狀況存在一定的差距，為此引入大氣梯度：

(2)

式中：ε為高度角；φ為方位角；GNS,GEW分別為南北、東西梯度。

由于改正模型只是根據(jù)多年的數(shù)據(jù)積累而建立的，并不能代表實際的對流層延遲，所以在沒有水汽計量觀測的情況下，提高對流層模擬折射精度的方法主要是參數(shù)估計法，常用的有單參數(shù)法、多參數(shù)法、隨機過程法和分段線性(PWL)方法，其實質都是一階高斯馬爾科夫過程的近似[6]。

GAMIT采用PWL方法，PWL用步長為kΔt的離散隨機過程表示對流層延遲隨時間的變化：

(3)

假設在歷元I～I+k之間，天頂對流層延遲隨時間線性變化：

(I≤j

(4)

適當選擇k值使求解的測站天頂方向的對流層折射參數(shù)ρ(I)，ρ(I+k)等的個數(shù)較少，從而可以用最小二乘的方法估計參數(shù)。

2 數(shù)據(jù)準備與處理策略

本文所采用的數(shù)據(jù)資料主要是陸態(tài)網絡連續(xù)站2015年的全網觀測資料，圖1為站點分布。

圖1 陸態(tài)網站點分布

數(shù)據(jù)處理是利用GAMIT軟件獲得測站和衛(wèi)星軌道的單日區(qū)域松弛解，該單日解給出了區(qū)域測站坐標、極移、天頂對流層延遲、衛(wèi)星相位中心偏差和衛(wèi)星參數(shù)的松弛解和方差-協(xié)方差矩陣。解算過程中用到的改正模型見表1。

3 結果與分析

3.1 全國范圍內的對流層延遲

在數(shù)據(jù)解算中每1 h設置一個對流層天頂延遲參數(shù)，衛(wèi)星截止高度角為15°，首先每天取一個平均值作為一天的天頂對流層延遲，然后對2015年的對流層延遲值取平均，得到全國2015年對流層天頂延遲分布，見圖2。

圖2為2015年全國陸態(tài)網絡站點對流層天頂延遲年平均量，從圖2中可以看出，全國整體對流層天頂延遲量在1.3～2.6 m之間，可分為三個層次；其中對流層天頂延遲較大值集中在東南沿海地區(qū)，大約2.4 m；其次是東北、內蒙、新疆及云南地區(qū)，大約為2.0 m；最小值集中在青藏高原地區(qū)。

表1 改正模型

圖2 2015年對流層天頂延遲年平均量

3.2 因素分析

天頂對流層總延遲包括濕延遲和干延遲，所以分別對2015年干、濕延遲的年平均量進行了計算，結果見圖3。

由數(shù)值上可以看出，干延遲占總延遲量的80%～90%，且干延遲量的分布特征與總延遲量的分布特征基本一致。眾所周知，對流層干延遲分量主要與大氣壓、溫度以及高程有關，所以本文計算了各站點年平均氣溫、大氣壓以及站點高程，結果如圖4所示。

由圖4中(b)(c)可以看出，氣壓的大小與高程的大小成反比關系，且兩者分布特征與干延遲量的分布特征基本一致；由圖4(a)和圖3(a)比較可以看出，東北、華北以及華南地區(qū)并沒有因為溫度的較大差異而導致這些地區(qū)干延遲量的明顯差異，說明對于這些地區(qū)決定其對流層干延遲大小的主要因素是氣壓與高程。

圖3 2015年干、濕對流層延遲年平均量

圖4 各站點氣溫、氣壓和高程分布

圖5 濕延遲與降水量分布圖

濕延遲量的影響因素主要是大氣中的水汽含量，由于水汽含量很難測定[7-11]，所以本文采用氣象局提供的年平均降水量分布來驗證反演對流層濕延遲的可靠性，結果如圖5所示。

由圖5可以看出，濕延遲量分為4個層次，與國家氣候中心發(fā)布的全國年平均降水量分布圖基本一致。濕延遲量最大的地區(qū)主要是南部沿海地區(qū)，其次是華南、四川、云南等地，第三是東北、山東、河南、四川北部等地，最后新疆、西藏、內蒙、甘肅、青海等地為第四層次。

4 結 論

通過本文的研究可以看出，天頂對流層延遲存在很強的地域性特征，其主要是受高程和氣壓因素的影響；天頂對流層延遲可分為干延遲和濕延遲，天頂對流層干延遲量占總延遲量的；利用GAMIT軟件對天頂對流層濕延遲進行反演，其結果與年平均降水量結果相吻合，說明該方法可以為氣象中短期預報等方面提供參考。

[1] 王勇, 聞德保, 劉嚴萍,等. 霧霾天氣對GPS天頂對流層延遲與可降水量影響研究[J]. 大地測量與地球動力學, 2014, 34(2):120-123.

[2] 曲建光, 劉基余, 韓中元. 利用天頂對流層延遲數(shù)據(jù)直接推算水汽含量的研究[J]. 武漢大學學報(信息科學版), 2005, 30(7):625-628.

[3] 唐超華, 李陶, 施闖,等. 天津地區(qū)CORS站天頂對流層延遲年周期變化研究[J]. 大地測量與地球動力學, 2009, 29(2):106-110.

[4] HAASE J S, VEDEL H, GE M, et al. GPS zenith tropospheric delay (ZTD) variability in the Mediterranean[J]. Physics & Chemistry of the Earth Part A Solid Earth & Geodesy, 2001, 26(6-8):439-443.

[5] 葉世榕, 張雙成, 劉經南. 精密單點定位方法估計對流層延遲精度分析[J]. 武漢大學學報(信息科學版), 2008, 33(8):788-791.

[6] 殷海濤, 黃丁發(fā), 熊永良,等. GPS信號對流層延遲改正新模型研究[J]. 武漢大學學報(信息科學版), 2007, 32(5):454-457.

[7] 王勇, 劉嚴萍, 柳林濤,等. 無氣象要素的GPS對流層延遲推算可降水量的研究[J]. 測繪科學, 2007, 32(3):122-124.

[8] VEDEL H, HUANG X, HAASE J, et al. Impact of GPS Zenith Tropospheric Delay data on precipitation forecasts in Mediterranean France and Spain[J]. Geophysical Research Letters, 2004, 31(31):231-254.

[9] 劉志成,于江. CORS在浮動車交通信息采集系統(tǒng)中的應用[J]. 交通科技與經濟,2013,14(6):59-61.

[10] HAASE J, VEDEL H, CALAIS E, et al. Accuracy and Variability of GPS Tropospheric Delay Measurements of Water Vapor in the Western Mediterranean[J]. Journal of Applied Meteorology, 2003, 42(11):1547-1568.

[11] 張隨甲,周曉忠. GPS-RTK誤差來源及減弱誤差的實用方法[J]. 交通科技與經濟,2014,15(2):121-124.

[12] 王勇,聞德保,劉嚴萍,等. 霧霾天氣對GPS天頂對流層延遲與可降水量影響研究[J]. 大地測量與地球動力學,2014(2):120-123，127.

[13] 王勇,焦健,曾琪明,等. 基于不同地形的GPS對流層延遲插值方法研究[J]. 大地測量與地球動力學,2010(3):132-136.

[14] 楊晶,顧慧,劉嚴萍,等. 基于小波分析的GPS對流層延遲變化研究[J]. 測繪通報,2015(5):60-62.

[15] 王富有,史先領. 實時氣象參數(shù)改正GPS定位中的對流層延遲[J]. 能源技術與管理,2006(6):132-134.

The distribution of zenith troposphere delay inverted in the stations of “Crustal Movement Observation Network in China”

WANG Dongzhen,ZHAO Bin,TAN Kai,ZHANG Caihong

(Key Lab of Earthquabe Geodesy,Institute of Seismology,China Earthquabe Administration,Wuhai 430071,China)

Aiming at the problem that the existing tropospheric delay studies have a small range and a short time span, the tropospheric delay value of 270 land-based network sites in China is researched by GAMIT software in this paper. First, the metutil program is called by the sh_metutil command in GAMIT to calculate and extract the zenith total delay, tropospheric delay, dry wet delay component and temperature component and pressure. Then, its annual mean amount of dry and wet delay is analyzed, and the characteristics of tropospheric delay distribution are obtained. Furthermore, the factors affecting the dry delay are analyzed, and the results show that the main influencing factors are barometric pressure and altitude factor. Finally, the relation between the wet delay and precipitation is analyzed, and the conclusion is that the distribution characteristics of them are basically the same.

GAMIT; Crustal Movement Observation Network in China;troposphere delay；the characteristics of distribution

2017-01-01

中國地震局地震研究所資助項目(IS201526244)

王東振(1989-)，男，實習研究員.

著錄：王東振，趙斌，譚凱，等.利用全國陸態(tài)網絡站點反演天頂對流層延遲分布特征[J].測繪工程，2018，27(2)：32-36.

10.19349/j.cnki.issn1006-7949.2018.02.006

P228.4

A

1006-7949(2018)02-0032-05

李銘娜]