韋 云 王 迅 王 浩 李 黎 陳國棟 趙 偉

1 蘇州科技大學(xué)地理科學(xué)與測繪工程學(xué)院，江蘇省蘇州市學(xué)府路99號，215009 2 蘇州科技大學(xué)北斗導(dǎo)航與環(huán)境感知研究中心，江蘇省蘇州市學(xué)府路99號，215009

近年來，利用由GNSS技術(shù)獲取的PWV時空變化趨勢來預(yù)測降雨的方法成為國內(nèi)外學(xué)者關(guān)注的重點[1-3]。傳統(tǒng)的獲取PWV的方法是利用由GNSS技術(shù)得到的對流層延遲(ZTD)減去干延遲(ZHD)，得到濕延遲(ZWD)，再利用ZWD推算得到PWV。計算ZHD需要的參數(shù)有測站處地面大氣壓(Ps)、測站緯度(θ)、測站大地高(H)，ZWD與PWV之間的轉(zhuǎn)換系數(shù)也要用到加權(quán)平均溫度(Tm)、地面氣溫(Ts)等參數(shù)，參數(shù)越多，計算過程越復(fù)雜，會導(dǎo)致數(shù)據(jù)量大、計算效率低，且易產(chǎn)生誤差累積等問題。因此，建立ZTD與GNSS-PWV之間的直接轉(zhuǎn)換模型，以計算特定時段的PWV時序信息十分必要[4-6]。

本文基于2017年東南沿海地區(qū)18個GNSS站的觀測數(shù)據(jù)，采用多元函數(shù)線性回歸法建立該地區(qū)GNSS-PWV與ZTD、地面氣溫(Ts)、地面大氣壓(Ps)之間的直接轉(zhuǎn)換模型，并利用模型預(yù)報2018年GNSS-PWV，以各站實測PWV為參考值，驗證本文模型的預(yù)報精度。

1 數(shù)據(jù)來源及其相關(guān)性分析

1.1 數(shù)據(jù)來源

選取由江蘇省氣象局和中國地震局GNSS數(shù)據(jù)產(chǎn)品服務(wù)平臺(http:∥www.cgps.ac.cn)提供的2017～2018年東南沿海地區(qū)18個GNSS站資料，主要包括PWV、ZTD、地面大氣壓(Ps)和地面氣溫(Ts)等，其中常州、黃山、蚌埠、武夷山4個氣象站用于精度檢驗，其余14個測站用于建模，測站空間位置信息見表1，選取的18個GNSS站均勻分布于108°～122°E、19°～33°N區(qū)域，涵蓋了東南沿海大部分地區(qū)。

表1 GNSS測站空間位置信息

中國地震局GNSS數(shù)據(jù)產(chǎn)品服務(wù)平臺采用Saastamoinen模型計算ZHD[7]，其表達式為：

(1)

用ZTD減去ZHD獲得ZWD，ZWD乘以轉(zhuǎn)換系數(shù)Π便可得到GNSS-PWV：

(2)

式中，ρw為液態(tài)水密度，k′2、k3為大氣折射常數(shù)，Rv為水汽氣體常數(shù)，Tm為加權(quán)平均溫度。

1.2 影響因子相關(guān)性驗證

皮爾森(Pearson)相關(guān)系數(shù)R可用來描述2個變量之間的線性相關(guān)程度，R的絕對值越大，表明其相關(guān)性越強，在0.8～1.0之間表示極強相關(guān)，0.6～0.8之間為強相關(guān)，0.4～0.6之間為中等程度相關(guān)，0.2～0.4之間為弱相關(guān)，0.0～0.2則無相關(guān)性[8]。

PWV與ZTD、Ts、Ps之間的相關(guān)程度如圖1所示，可以看出，ZTD與PWV的相關(guān)系數(shù)為0.98，為極強相關(guān)，表明ZTD對PWV值有極大影響；PWV與Ts的相關(guān)系數(shù)為0.78，為強相關(guān)；PWV與Ps的相關(guān)系數(shù)為-0.65，為負強相關(guān)。

圖1 相關(guān)性分析

表2為各變量間相關(guān)性統(tǒng)計結(jié)果，可以看出，各自變量之間也存在一定的線性關(guān)系。若建模選取的樣本自變量之間本身就存在高度共線性關(guān)系，會導(dǎo)致線性回歸模型不穩(wěn)定，難以區(qū)分各自變量對模型結(jié)果的影響，因此需要對自變量進行共線性分析。

表2 相關(guān)性分析統(tǒng)計

本文利用方差膨脹因子(variance inflation factor，VIF)分析各自變量之間的共線性，其表達式為：

VIFi=1/(1-R2)

(3)

式中，R為自變量x與其余自變量作回歸分析的復(fù)相關(guān)系數(shù)。通常認為，當(dāng)0

2 PWV直接轉(zhuǎn)換模型的建立及精度分析

2.1 建模方法

利用多元函數(shù)線性擬合方法建立PWV直接轉(zhuǎn)換模型。設(shè)多元線性回歸模型為：

y=β0+β1x1+…+βmxm+ε，ε～N(0,σ2)

(4)

式中，y為因變量；x1、x2、…、xm為自變量；β0為常數(shù)項；β1、…、βm為回歸系數(shù)，是與自變量無關(guān)的未知參數(shù)；ε為誤差項[10]，是均值為0、方差σ2>0的不可觀測隨機變量。若有n個獨立觀測數(shù)據(jù)，則可列出n個關(guān)于未知參數(shù)x1、x2、…、xm的方程：

yi=β0+β1xi1+…+βmxim+εi

εi～N(0,σ2),i=1,…,n

(5)

則有：

PWV轉(zhuǎn)換模型可表示為：

Y=Xβ+ε，ε～N(0,σ2En)

(6)

式中，En為n階單位矩陣。該線性回歸模型可用最小二乘法原理求解，其中自變量的系數(shù)βi為最小二乘估計值。

2.1.1 多因子模型

基于PWV與ZTD、Ts、Ps的線性關(guān)系，PWV與其他參數(shù)函數(shù)的關(guān)系為：

PWV=β0+β1ZTD+β2Ts+β3Ps

(7)

將2017年14個測站的PWV與ZTD、Ts、Ps代入式(7)，利用最小二乘法可得到多因子PWV的直接轉(zhuǎn)換模型為：

PWV=-13.380 5+0.162 9ZTD+

0.074 1Ts-0.359 7Ps

(8)

2.1.2 雙因子模型

將2017年14個測站的PWV與ZTD、Ts代入式(7)，得到基于ZTD和Ts的雙因子PWV直接轉(zhuǎn)換模型：

PWV=-377.315 7+0.161 0ZTD+

0.415 7Ts

(9)

將PWV與ZTD、Ps代入式(7)，得到基于ZTD和Ps的雙因子PWV直接轉(zhuǎn)換模型：

PWV=3.335 4+0.165 8ZTD-0.3821Ps

(10)

2.1.3 單因子模型

假設(shè)PWV和ZTD之間的線性關(guān)系為：

PWV=β0+β1ZTD

(11)

將2017年14個測站的PWV和ZTD代入式(11)，得到基于ZTD的單因子PWV直接轉(zhuǎn)換模型：

PWV=-423.991 2+0.182 7ZTD

(12)

2.2 精度檢驗

利用2017年14個測站的數(shù)據(jù)建立PWV模型，分別將2017～2018年的數(shù)據(jù)代入所建模型中，計算得到PWV模型預(yù)測值，并與原始數(shù)據(jù)中的PWV實測值(作為真值)進行比較，驗證模型的精度。

本文將PWV實測值減去PWV模型預(yù)測值，計算得到真值與預(yù)測值之間的偏差(bias)，通過統(tǒng)計分析平均bias與RMS值來驗證PWV模型的精度。圖2為各PWV模型的bias變化趨勢，可以看出，單因子PWV模型預(yù)測值與真值的bias基本在15 mm以內(nèi)；雙因子(無Ps)PWV模型預(yù)測值與真值的bias在10 mm以內(nèi)；雙因子(無Ts)PWV模型預(yù)測值與真值的bias在2 mm以內(nèi)；多因子PWV模型預(yù)測值與真值的bias在1 mm以內(nèi)，精度最高。2018年的數(shù)據(jù)未參與建模，PWV模型預(yù)測值及bias的變化趨勢與2017年情況基本一致。

圖2 2017～2018年不同PWV模型預(yù)測值偏差

圖3為未參與建模的常州、黃山、蚌埠和武夷山等4個測站2017～2018年各PWV模型預(yù)測值bias。由圖可知，常州站和蚌埠站基于雙因子(無Ps)模型和單因子模型的預(yù)測值大部分大于真值，而基于雙因子(無Ts)模型的預(yù)測值大部分小于真值，基于多因子模型的預(yù)測值則均勻分布于真值兩側(cè)；黃山站和武夷山站基于雙因子(無Ps)模型和單因子模型的預(yù)測值均小于真值，基于雙因子(無Ts)模型的預(yù)測值均大于真值，而基于多因子模型的預(yù)測值同樣均勻分布于真值兩側(cè)。從bias的變化范圍來看，4個測站的PWV偏差值與參與建模的PWV偏差值分布基本一致，各類模型的bias變化范圍也與圖2非常接近。

圖3 2017～2018年各測站PWV模型預(yù)測值偏差

表3(單位mm)為各PWV模型預(yù)測值的平均bias和RMS統(tǒng)計結(jié)果，可以看出，基于單因子模型的RMS值為4.66 mm，平均bias為3.73 mm；加入Ts的雙因子模型RMS值降至3.94 mm，平均bias為3.13mm，精度有小幅提升；加入Ps的雙因子模型RMS值降至0.50 mm，平均bias為0.40 mm，精度大幅提高；同時加入Ps和Ts的多因子模型RMS值降至0.33 mm，平均bias降為0.24 mm，精度達到最高。

表3 全年P(guān)WV模型精度統(tǒng)計

由式(1)和式(2)可知，ZHD主要通過Ps計算得到，ZWD主要通過Ts及Tm計算得到，而Tm也由Ts計算得出。由于ZHD約占對流層延遲的80%～90%，加入Ps的雙因子模型能提升ZHD的計算精度，而ZWD在ZTD中的占比較小，因此基于ZTD與Ps的雙因子PWV模型的精度明顯比基于ZTD和Ts的雙因子PWV模型的精度高。

與基于ZTD、Ts和Ps的多因子PWV模型相比，基于ZTD和Ps的雙因子PWV模型的精度稍低，但其RMS值有0.50 mm，可滿足絕大部分GNSS氣象學(xué)研究的需求，且其計算參數(shù)較少，轉(zhuǎn)換效率較高。

3 結(jié) 語

本文利用2017～2018年東南沿海地區(qū)18個GNSS站的數(shù)據(jù)，通過分析PWV與ZTD、地面氣溫(Ts)、地面大氣壓(Ps)等參數(shù)的相關(guān)性，基于多元線性擬合方法構(gòu)建東南沿海地區(qū)的多因子PWV直接轉(zhuǎn)換模型，得出以下結(jié)論：

1)PWV與ZTD之間的相關(guān)性最強，相關(guān)系數(shù)達到0.98；PWV與Ps、Ts之間也具有較強的相關(guān)性，相關(guān)系數(shù)分別為-0.65和0.78。

2)基于ZTD的單因子PWV模型平均bias為3.73 mm，RMS為4.66 mm；基于ZTD和Ts的雙因子PWV模型平均bias為3.13 mm，RMS為3.94 mm；基于ZTD和Ps的雙因子PWV模型平均bias為0.40 mm，RMS為0.50 mm；基于ZTD、Ts和Ps的多因子PWV模型平均bias為0.24 mm，RMS為0.33 mm。多因子PWV直接轉(zhuǎn)換模型的精度最高，基于ZTD和Ps的雙因子模型次之，這兩種模型均可獲得精度優(yōu)于1 mm的PWV結(jié)果。

3)東南沿海地區(qū)的PWV直接轉(zhuǎn)換模型不僅可應(yīng)對傳統(tǒng)PWV計算過程數(shù)據(jù)量大、效率低且易產(chǎn)生誤差累積的問題，也可在基本地面氣象數(shù)據(jù)缺失時，基于ZTD直接計算特定時段的PWV時序信息，以用于GNSS氣象學(xué)研究。