羅 宇，羅林艷，范嘉智，段思汝

(1. 中國(guó)氣象局氣象干部培訓(xùn)學(xué)院湖南分院，湖南 長(zhǎng)沙 410125；2. 湖南省氣象信息中心，湖南 長(zhǎng)沙 410118)

地球水資源極為豐富，總儲(chǔ)量約1.386×1018m3，其中以水汽形式存在的僅占總量的0.001%，卻在各種大氣活動(dòng)中扮演重要角色，是天氣和氣候變化的主要驅(qū)動(dòng)力。理論上，大氣水汽含量可反映一個(gè)地區(qū)在不穩(wěn)定層結(jié)條件下的可降水量，因此高時(shí)空分辨率的水汽資料對(duì)于大氣水循環(huán)監(jiān)測(cè)、災(zāi)害性天氣監(jiān)測(cè)、中小尺度數(shù)值天氣預(yù)報(bào)和強(qiáng)對(duì)流天氣短臨預(yù)報(bào)等極為重要[1]。但水汽時(shí)空分布不均勻，變化迅速，作為常規(guī)探測(cè)手段的無(wú)線電探空站點(diǎn)水平距離平均在200～300 km，且每天一般只進(jìn)行兩次觀測(cè)，無(wú)法獲得高時(shí)空分辨率的大氣水汽信息。20世紀(jì)90年代基于地基全球定位系統(tǒng)(GPS)探測(cè)水汽技術(shù)的出現(xiàn)彌補(bǔ)了常規(guī)大氣水汽探測(cè)手段的不足，可提供快速變化的水汽信息，展現(xiàn)出廣泛的應(yīng)用前景[2-3]。相關(guān)研究表明，GPS可降水量(GPSPWV)探測(cè)具有較高精度。Duan等[3]改進(jìn)了可降水量(PWV)計(jì)算方法，得到GPSPWV探測(cè)精度為1.0～1.5 mm。杜明斌等[4]對(duì)比分析2010年華東區(qū)域地基GPS/MET資料和無(wú)線電探空資料，認(rèn)為二者探測(cè)所得可降水量具有很好的一致性，平均誤差在1 mm以下。羅宇等[5]對(duì)比湖南省2016年4月和7月GPS/MET和無(wú)線電探空水汽資料，認(rèn)為兩者有較強(qiáng)相關(guān)性。湖南省氣象局與湖南省國(guó)土資源廳合作建設(shè)了GPS/MET探測(cè)網(wǎng)，可為災(zāi)害性天氣監(jiān)測(cè)預(yù)警、氣候預(yù)測(cè)提供時(shí)空分辨率更高的水汽初始場(chǎng)。

GPS衛(wèi)星信號(hào)穿過(guò)大氣層，受大氣折射影響產(chǎn)生電離層延遲和中性大氣延遲，其中電離層延遲可以采用雙頻觀測(cè)技術(shù)基本消除。中性大氣延遲可利用GPS數(shù)據(jù)處理軟件解算得到，其由靜力延遲(ZHD)和濕延遲(ZWD)兩部分組成，其中ZHD一般采用基于理想大氣條件推導(dǎo)出的天頂靜力延遲模型計(jì)算得到。但實(shí)際大氣溫濕層千變?nèi)f化，采用模型計(jì)算得到的ZHD可能存在系統(tǒng)誤差，進(jìn)而影響GPSPWV反演精度。

針對(duì)上述問(wèn)題，本文利用湖南省GPS/MET探測(cè)網(wǎng)數(shù)據(jù)及相關(guān)氣象資料，對(duì)天頂靜力模型誤差進(jìn)行分析并改進(jìn)，提升湖南省GPSPWV反演精度，使其更好地應(yīng)用到降水預(yù)報(bào)、防災(zāi)減災(zāi)等氣象業(yè)務(wù)中。

1 數(shù)據(jù)與方法

湖南省地基GPS/MET站網(wǎng)設(shè)計(jì)站點(diǎn)平均間距47.7 km，可對(duì)湖南省大氣水汽狀況進(jìn)行實(shí)時(shí)監(jiān)測(cè)。為檢驗(yàn)天頂靜力延遲模型對(duì)GPSPWV反演精度的影響及改進(jìn)效果，本文選取長(zhǎng)沙、懷化和郴州3個(gè)高空氣象探測(cè)站探空秒級(jí)資料和湘鄉(xiāng)、鳳凰和桂陽(yáng)3個(gè)GPS/MET站數(shù)據(jù)作為對(duì)比分析樣本(長(zhǎng)沙、懷化和郴州未建設(shè)GPS/MET觀測(cè)站，因此選取直線距離最近的GPS /MET站數(shù)據(jù)進(jìn)行對(duì)比)，測(cè)站具體信息如表1所示。資料的時(shí)間年限為2016-04-01—2017-03-31，其中探空秒級(jí)資料每天2次，GPS/MET資料每小時(shí)1次。

表1 高空氣象觀測(cè)站和GPS/MET站基本信息

GPS信號(hào)穿過(guò)大氣時(shí)由于折射產(chǎn)生延遲，大氣折射指數(shù)可表示為

N=k1Pd/T+k2e/T+k3e/T2.

(1)

式中:Pd為干空氣氣壓(hPa);e為水汽壓(hPa);T為絕對(duì)溫度(K);k1，k2和k3為實(shí)驗(yàn)室常數(shù)，取值分別為k1=77.689 0 K/hPa，k2=71.295 2 K/hPa，k3=375 463 K2/hPa。式(1)中k1Pd/T稱(chēng)為干項(xiàng)，對(duì)其積分可得天頂靜力延遲(ZHD)，k2e/T+k3e/T2稱(chēng)為濕項(xiàng)，對(duì)其積分可得天頂濕延遲(ZWD)，二者之和為天頂總延遲(ZTD)。在高精度的GPS數(shù)據(jù)處理中，一般利用Bernese或GAMIT等解算軟件求得ZTD后，由天頂靜力模型計(jì)算得到ZHD，相減得到ZWD。目前常用的普適ZHD模型主要有Black模型[6]、Hopfield模型[7]和Saastamoinen模型[8](以下分別簡(jiǎn)稱(chēng)B模型、H模型和SA模型)，其算式分別為

ZHDB=2.343×(Ts-4.12)×Ps/Ts,

(2)

ZHDH=15.52×Ps×[40.082+

0.148 98×(Ts-273.16)-Hs]/Ts,

(3)

ZHDSA=2.276 8×Ps/

(1-0.002 66×cos 2θ-0.000 28×Hs).

(4)

式中：ZHDB、ZHDH和ZHDSA分別對(duì)應(yīng)B模型、H模型和SA模型ZHD(mm)，θ，Hs，Ps和Ts分別為測(cè)站緯度(°)、高度(km)、氣壓(hPa)和絕對(duì)溫度(K)。由于上述ZHD模型均在理想大氣和靜力平衡假設(shè)下推導(dǎo)，存在因系數(shù)不準(zhǔn)確產(chǎn)生模型誤差可能性。因此，利用長(zhǎng)沙、懷化和郴州探空秒級(jí)資料對(duì)式(1)中干項(xiàng)進(jìn)行積分，得到探空Z(yǔ)HD后，對(duì)式(2)、式(3)和式(4)進(jìn)行線性回歸建模，訂正模型系數(shù)，并利用湘鄉(xiāng)、鳳凰和桂陽(yáng)三站的GPS/MET數(shù)據(jù)對(duì)訂正前后的ZHD模型進(jìn)行檢驗(yàn)。訂正模型為

ZHD=γ0+γ1×f.

(5)

其中f為與氣象要素相關(guān)的函數(shù)，與B模型、H模型和SA模型對(duì)應(yīng)的f分別為

f=(Ts-4.12)×Ps/Ts,

(6)

f=Ps×[40.082+0.148 98×

(Ts-273.16)-Hs]/Ts,

(7)

f=Ps/(1-0.002 66×cos 2θ-0.000 28×Hs).

(8)

在此基礎(chǔ)上反演GPSPW，與探空PWV進(jìn)行對(duì)比分析，評(píng)估ZHD模型對(duì)GPS水汽反演精度的影響。GPSPW和探空PWV的算式分別為[2]

(9)

(10)

(11)

相關(guān)研究顯示Tm與Ts呈現(xiàn)高度線性相關(guān)[2-3]，即Tm=a+bTs，采用基于湖南省2004—2015年21 935次探空觀測(cè)擬合所得的本地化Tm公式，其中a=105.940 5，b=0.597 4。

考慮到國(guó)內(nèi)L波段探空系統(tǒng)濕度資料存在整體偏干[23]及中低空偶發(fā)極度偏干現(xiàn)象[24]，有必要對(duì)探空濕度數(shù)據(jù)進(jìn)行質(zhì)控和訂正。300 hPa高度以下，探空相對(duì)濕度數(shù)據(jù)小于5%且連續(xù)出現(xiàn)厚度達(dá)到200 hPa以上時(shí)認(rèn)為該探空濕度數(shù)據(jù)異常，不納入對(duì)比分析數(shù)據(jù)集。在此基礎(chǔ)上對(duì)大于60%的相對(duì)濕度數(shù)據(jù)進(jìn)行訂正，訂正分段函數(shù)為

(12)

式中：RHob為探空相對(duì)濕度觀測(cè)值；RHthr為相對(duì)濕度訂正閾值，取60%。考慮到對(duì)流層以上大氣對(duì)天頂靜力延遲存在15%左右貢獻(xiàn)[8]，對(duì)探空觀測(cè)高度上界(一般為30 km左右)至100 km高度間大氣層采用中國(guó)標(biāo)準(zhǔn)大氣模型[11](見(jiàn)表2)進(jìn)行續(xù)補(bǔ)后計(jì)算探空Z(yǔ)HD。

表2 中國(guó)標(biāo)準(zhǔn)大氣模型層面

各層大氣溫度和氣壓計(jì)算方法為

(13)

式中：γ為氣溫梯度(K/km);R為理想氣體常數(shù);m為氣體摩爾質(zhì)量;g為重力加速度;hi-1，hi，Pi-1，Pi和Ti-1，Ti分別為第i層大氣下界和上界高度(m)、氣壓(hPa)和溫度(K)。

2 結(jié)果分析

2.1 ZTD模型的改進(jìn)

利用2016-04-01—2017-03-31長(zhǎng)沙、懷化和郴州3站探空秒級(jí)資料計(jì)算得到探空Z(yǔ)HD，與基于式(2)～式(4)計(jì)算得到的ZHD進(jìn)行對(duì)比分析，可見(jiàn)探空Z(yǔ)HD和模型ZHD具有非常顯著的相關(guān)性(相關(guān)系數(shù)均大于0.97)，除郴州站H模型ZHD偏小外，其他模型ZHD普遍較探空Z(yǔ)HD偏大，其中又以B模型平均偏差最大，3站分別平均偏大45.529 0 mm以上。為減小ZHD模型的系統(tǒng)偏差，根據(jù)式(5)～式(8)，利用探空Z(yǔ)HD進(jìn)行線性回歸建模，得到改進(jìn)后的ZHD模型(見(jiàn)表3)。改進(jìn)后的模型ZHD與探空Z(yǔ)HD相比，平均偏差(Bias)均為0 mm，均方誤差(RMSE)由5.230 5～6.615 5 mm降至3.534 9～4.633 8 mm，平均降低32.46%，明顯降低模型誤差。

表3 探空Z(yǔ)HD與改進(jìn)后前后的模型ZHD對(duì)比

2.2 GPS可降水量精度分析

為研究改進(jìn)后的ZHD模型對(duì)GPSPWV反演的影響，選取距離長(zhǎng)沙、懷化和郴州高空氣象探測(cè)站最近的GPS/MET站(湘鄉(xiāng)、鳳凰、桂陽(yáng))數(shù)據(jù)，結(jié)合湖南省本地化Tm模型反演GPSPWV，并與基于式(9)得到的探空PWV進(jìn)行對(duì)比分析。

通過(guò)對(duì)長(zhǎng)沙-湘鄉(xiāng)(410個(gè)樣本)、懷化-鳳凰(390個(gè)樣本)和郴州-桂陽(yáng)(366個(gè)樣本)的數(shù)據(jù)對(duì)比分析可知，GPSPWV與探空PWV存在良好相關(guān)性，兩者相關(guān)系數(shù)均在0.971 4以上。由圖1可知，相較于探空PWV，采用改進(jìn)前的3種常用ZHD模型反演得到的GPSPWV平均分別偏低7.212 3 mm、0.800 3 mm和1.585 0 mm以上，均方誤差在4.357 9 mm以?xún)?nèi)，其中B模型反演GPSPWV系統(tǒng)偏差最大。采用改進(jìn)后ZHD模型反演GPSPWV相較于探空PWV的平均偏差降至±0.5 mm以?xún)?nèi)，降低59.54%以上，均方誤差平均降低4%左右。由此可見(jiàn)改進(jìn)后的3種ZHD模型均可較明顯降低GPSPWV反演系統(tǒng)偏差，并略微降低均方誤差，提升湖南GPSPWV反演精度。

3 結(jié) 論

1)利用普適ZHD模型(B模型、H模型和SA模型)計(jì)算ZHD，與基于湖南長(zhǎng)沙、懷化和郴州3個(gè)高空氣象探測(cè)站實(shí)測(cè)ZHD進(jìn)行對(duì)比，發(fā)現(xiàn)ZHD模型普遍存在系統(tǒng)偏差，其中B模型計(jì)算值較探空實(shí)測(cè)ZHD平均偏大40 mm以上。

2)利用湖南探空實(shí)測(cè)ZHD對(duì)普適ZHD模型進(jìn)行訂正，可明顯減小模型平均偏差和均方誤差，因此有必要利用GPS/MET測(cè)站所在地或附近的氣象探空資料建立訂正模型，改進(jìn)ZHD模型精度。

3)利用普適ZHD模型反演湖南省GPSPWV存在較明顯的系統(tǒng)偏差，其中B模型系統(tǒng)偏差最大，平均偏小7.212 3 mm以上，H模型和SA模型反演精度相當(dāng)；改進(jìn)后的ZHD模型可明顯減小GPSPWV反演的系統(tǒng)偏差，并略微降低均方誤差，其中采用改進(jìn)后的SA模型反演所得GPSPWV與探空PWV相比，平均偏差不超過(guò)0.230 7 mm，可較明顯提升GPSPWV反演準(zhǔn)確性。

(a)長(zhǎng)沙-湘鄉(xiāng) (b)懷化-鳳凰 (c)郴州-桂陽(yáng)圖1 改進(jìn)前后的ZHD模型對(duì)GPSPWV反演精度的影響

[1] RIESCO J,MORA M,DE PABLO F,et al.Regimes of intense precipitation in the Spanish Mediterranean area[J].Atmospheric Research,2014,137:66-79.

[2] BEVIS M,BUSINGER S,HERRING T A,et al.GPS meteorology: Remote sensing of atmosphere water vapor using the global positioning system[J].Journal of Geophysical Research,1992,97(D14): 15787-15801

[3] DUAN J,BEVIS M,FANG P.GPS meteorology: Direct estimation of the absolute value of precipitable water[J].Journal of Applied Meteorology,1996,35(7): 830-838.

[4] 杜明斌,尹球,劉敏,等.地基GPS/MET探測(cè)水汽等相關(guān)參數(shù)精度分析[J].大氣與環(huán)境光學(xué)學(xué)報(bào),2013,8(2): 138-145.

[5] 羅宇,羅林艷,呂冠儒.湖南省GPS/MET與探空水汽資料對(duì)比分析[J].貴州氣象,2016,40(6): 28-31.

[6] BLACK H D.An easily implemented algorithm for the tropospheric range correction[J].Journal of Geophysical Research,1978,83(B4): 1825-1828.

[7] HOPFIELD H S.Tropospheric effect on electromagnetically measured range: Prediction from surface weather data[J].Radio Science,1971,6(3): 357-367.

[8] SAASTAMOINEN J.Contributions to the theory of atmospheric refraction[J].Bulletin Géodésique,1973,107(1): 13-34.

[9] 郝民,龔建東,王瑞文,等.中國(guó)L波段探空濕度觀測(cè)資料的質(zhì)量評(píng)估及偏差訂正[J].氣象學(xué)報(bào),2015,73(1): 187-199.

[10] 唐南軍,劉艷,李剛,等.中低空探空相對(duì)濕度觀測(cè)數(shù)據(jù)的新問(wèn)題——基于中國(guó)L波段探空系統(tǒng)適度觀測(cè)異常偏干現(xiàn)象的初步分析[J].熱帶氣象學(xué)報(bào),2014,30(4): 643-653.