孫天紅，高志鈺，蔣玉祥

（1. 西寧不動(dòng)產(chǎn)登記中心，西寧 810001；2. 蘭州交通大學(xué) 測(cè)繪與地理信息學(xué)院，蘭州 730070；3. 西寧市測(cè)繪院，西寧 810001）

0 引言

可降水量（precipitable water vapor，PWV）直接影響著大氣的垂直穩(wěn)定性，其時(shí)空分布情況與實(shí)際降雨量密切相關(guān)，是中尺度和局部尺度降雨強(qiáng)度預(yù)報(bào)的必要參數(shù)之一[1-2]。隨著全球定位系統(tǒng)（global positioning system，GPS）技術(shù)的不斷發(fā)展，利用地基GPS 接收設(shè)備，可以測(cè)量和反演出整層大氣中的PWV，這為天氣預(yù)報(bào)、天氣變化研究、合成孔徑雷達(dá)干涉（interferometric synthetic aperture radar，InSAR）大氣效應(yīng)去除等涉及大氣水汽工作提供了1 種新的PWV 數(shù)據(jù)獲取手段。

在利用地基 GPS 技術(shù)進(jìn)行 PWV 的反演過(guò)程中，天頂濕延遲（zenith wet delay，ZWD）與轉(zhuǎn)換系數(shù)Π直接影響著實(shí)時(shí)PWV 精確反演，即PWV=無(wú)量綱的轉(zhuǎn)換系數(shù)Π是大氣加權(quán)平均溫度Tm的函數(shù)，在PWV 的計(jì)算過(guò)程中，起決定性作用的是加權(quán)平均溫度Tm的準(zhǔn)確計(jì)算[3-4]，獲取適合局地地區(qū)的加權(quán)平均溫度模型，才能更加有效地保證局地地區(qū) PWV 反演的精度。不少學(xué)者就不同地區(qū)大氣加權(quán)平均溫度進(jìn)行分析，建立了加權(quán)平均溫度模型。文獻(xiàn)[5]首次提出利用探空數(shù)據(jù)，建立了適合美國(guó)中緯度地區(qū)的Bevis 經(jīng)驗(yàn)公式，建立了地表溫度Ts與Tm之間的線性關(guān)系。文獻(xiàn)[6]探討了全球Ts與Tm之間的相關(guān)關(guān)系，并建立了與其全球緯度相關(guān)的線性模型（本文稱(chēng)為姚宜斌模型）。文獻(xiàn)[7]對(duì)全國(guó)按氣候分區(qū)、季節(jié)分區(qū)分別建立了加權(quán)平均溫度模型（本文稱(chēng)為龔紹琦模型）。文獻(xiàn)[8]建立了適合于中國(guó)東北部地區(qū)和不同季節(jié)的加權(quán)平均溫度模型。文獻(xiàn)[9]利用湖南地區(qū)探空數(shù)據(jù)，建立了適合于湖南地區(qū)的加權(quán)平均溫度模型。文獻(xiàn)[10]利用甘肅地區(qū)探空數(shù)據(jù)，進(jìn)行了甘肅地區(qū)加權(quán)平均溫度時(shí)空分布特征分析，并建立了適合于甘肅地區(qū)的加權(quán)平均溫度模型。文獻(xiàn)[11]利用廣西地區(qū)探空數(shù)據(jù)，建立了廣西地區(qū)加權(quán)平均溫度模型。文獻(xiàn)[12]建立了適合于新疆地區(qū)的加權(quán)平均溫度模型。此外，文獻(xiàn)[13-17]分別建立了適合于鄭州地區(qū)、徐州地區(qū)、青島地區(qū)、香港地區(qū)的加權(quán)平均溫度模型，分別為局地地區(qū) PWV 的反演提供依據(jù)。通過(guò)大量的閱讀文獻(xiàn)發(fā)現(xiàn)，適用于銀川地區(qū)的加權(quán)平均溫度模型較少，已有的全球、全國(guó)、分區(qū)加權(quán)平均溫度模型，在銀川地區(qū)的適用性也有待進(jìn)一步檢驗(yàn)。

因此，本文針對(duì)適合于銀川地區(qū)的加權(quán)平均溫度模型較少的問(wèn)題，采用了2008—2017 年為期10 a的探空數(shù)據(jù)，建立起適合于銀川地區(qū)的大氣加權(quán)平均溫度模型，以實(shí)現(xiàn)能更準(zhǔn)確、更可靠地反演出銀川地區(qū)的PWV，從為該地區(qū)氣象應(yīng)用提供依據(jù)。

1 加權(quán)平均溫度的計(jì)算方法

水汽大部分分布在地面上空12 km 內(nèi)的區(qū)域，而由氣球攜帶的無(wú)線電探空儀可以升達(dá) 30 km 以上，探空資料提供的變量包括各標(biāo)準(zhǔn)層的氣壓、高度、溫度、露點(diǎn)溫度、風(fēng)向和風(fēng)速等數(shù)據(jù)。加權(quán)平均溫度Tm（單位為 K）可由測(cè)站上空的水汽壓e（單位為hPa）與大氣溫度T（單位為K）沿天頂方向的積分得到，其定義為

式中h為大氣厚度，單位為m。

常見(jiàn)的加權(quán)平均溫度的計(jì)算方法主要有以下幾種：常數(shù)法、近似積分法、數(shù)值積分法以及經(jīng)驗(yàn)公式法。其中，數(shù)值積分法是最常用的1 種方法，其計(jì)算公式為

式中：iT為第i層大氣平均溫度（單位為 K）；Δih為第i層大氣厚度（單位為m）；ei為第i層大氣平均水汽壓（單位為hPa），計(jì)算式為

式中：es為飽和水汽壓（單位為hPa）；t為露點(diǎn)溫度（單位為℃）；r為相對(duì)濕度。

經(jīng)查閱文獻(xiàn)，整理出符合銀川地區(qū)的經(jīng)驗(yàn)回歸公式。文獻(xiàn)[5]根據(jù)美國(guó)13 個(gè)氣象探空站約2 a 的探測(cè)資料，建立起Tm與地面溫度Ts（單位為 K）的經(jīng)驗(yàn)公式，該模型的精度為4.74 K。由于該公式簡(jiǎn)便、實(shí)用，在全球GPS 氣象研究中得到了廣泛使用。該經(jīng)驗(yàn)公式為

文獻(xiàn)[6]利用 2005—2011 年全球大地觀測(cè)系統(tǒng)大氣的Tm數(shù)據(jù)和歐洲中尺度天氣預(yù)報(bào)中心的Ts數(shù)據(jù)，按緯度建立了全球分區(qū)域線性回歸模型，適用于銀川地區(qū)的模型可表示為

文獻(xiàn)[7]利用全國(guó)不同氣候區(qū)、不同季節(jié)數(shù)據(jù)建立全國(guó)模型，普遍適用于我國(guó)陸地上空的Tm的估算，其模型可表示為

在GPS 氣象學(xué)中，數(shù)值積分法是目前公認(rèn)的、最為精確的計(jì)算對(duì)流層加權(quán)平均溫度模型的方法，利用數(shù)值積分法，可計(jì)算得到每個(gè)采樣時(shí)刻的Tm，結(jié)合與Ts之間的線性關(guān)系建立線性回歸公式。這也是獲取加權(quán)平均溫度本地化模型的重要手段，其主要原因在于[18]：①地面溫度Ts較其他地面氣象因子更容易獲取，Ts與Tm之間有顯明的線性關(guān)系；②GAMIT、Bernese 等高精度的全球衛(wèi)星導(dǎo)航系統(tǒng)（global navigation satellite system, GNSS）數(shù)據(jù)處理軟件，采用的加權(quán)平均溫度公式就是Tm與Ts之間的關(guān)系式。

2 加權(quán)平均溫度模型構(gòu)建方法

2.1 數(shù)據(jù)來(lái)源

選取銀川探空站（詳細(xì)信息如表 1 所示）2008—2017 年的探空數(shù)據(jù)作為實(shí)驗(yàn)數(shù)據(jù)。探空數(shù)據(jù)采樣間隔為12 h，其主要以文本的形式存儲(chǔ)，且只能按月進(jìn)行下載。因?yàn)閿?shù)據(jù)量大、下載過(guò)程繁瑣，所以本文通過(guò)Python 編寫(xiě)腳本的方式進(jìn)行批量下載，并按月份分別進(jìn)行存儲(chǔ)。

由于某些時(shí)刻探空氣球未釋放或者數(shù)據(jù)采集不完整，造成某些時(shí)段探空數(shù)據(jù)的缺失，2008—2017 年只采集到了7 193 個(gè)樣本的數(shù)據(jù)。

表1 銀川探空站詳細(xì)信息

2.2 建模方法

首先需對(duì)同一探空站對(duì)應(yīng)的Tm與Ts的相關(guān)性進(jìn)行計(jì)算，其計(jì)算公式為

然后采用一元線性回歸的方式進(jìn)行加權(quán)平均溫度模型的構(gòu)建，其線性方程為

則其誤差方程可表示為

最后根據(jù)最小二乘原理可得回歸系數(shù)計(jì)算式分別為：

3 銀川地區(qū)加權(quán)平均溫度模型的建立及精度檢驗(yàn)

采用數(shù)值積分方法，可以得到各采樣時(shí)刻的加權(quán)平均溫度值，結(jié)合相應(yīng)的地面溫度值，就能夠建立起銀川地區(qū)加權(quán)平均溫度模型。為檢核所建模型的可靠性，本文采用所建模型與 3 種已有模型進(jìn)行精度比對(duì)的方法。

3.1 Tm 與 Ts 之間的關(guān)系

利用MATLAB 軟件編寫(xiě)代碼，對(duì)探空數(shù)據(jù)進(jìn)行讀??；結(jié)合數(shù)值積分法，計(jì)算得到每個(gè)樣本所對(duì)應(yīng)的加權(quán)平均溫度Tm；繪制Tm與Ts之間的時(shí)間演變關(guān)系，如圖1 所示。

圖1 銀川探空站Tm 與Ts 時(shí)間演變關(guān)系

從圖 1 可以看出，Tm普遍低于Ts，Tm與Ts之間的變化趨勢(shì)保持一致，升降幅度也保持同步，峰值和谷值亦對(duì)應(yīng)得，可以判斷Tm與Ts之間具有良好的相關(guān)性。繪制銀川探空站Tm與Ts散點(diǎn)分布及趨勢(shì)線，如圖2 所示。

圖2 銀川探空站Tm 與Ts 散點(diǎn)分布及趨勢(shì)線

由圖2 可以看出，所有散點(diǎn)均分布在趨勢(shì)線附近，且在趨勢(shì)線上下波動(dòng)，這也證明了Tm與Ts之間具有良好的相關(guān)性。

為了進(jìn)一步地探討Tm與Ts之間的相關(guān)性，計(jì)算 2 者之間的相關(guān)系數(shù)R，計(jì)算結(jié)果如表 2所示。2008—2017 年間，Tm與Ts的相關(guān)系數(shù)為 0.804 3，通過(guò)了置信度α=0.01 的檢驗(yàn)。根據(jù)相關(guān)系數(shù)的劃分等級(jí)，當(dāng) 0.7≤R＜1 時(shí)表示強(qiáng)相關(guān)，即銀川地區(qū)加權(quán)平均溫度與地面溫度具有強(qiáng)相關(guān)性，可為后續(xù)加權(quán)平均溫度模型的建立提供依據(jù)。

3.2 銀川地區(qū)加權(quán)平均溫度模型

利用MATLAB 軟件，根據(jù)式（11）、式（12）計(jì)算出銀川地區(qū)Tm與Ts之間的回歸系數(shù)（即銀川地區(qū)加權(quán)平均溫度模型系數(shù)），其結(jié)果如表2 所示。Tm=61.884 9Ts+0.738 5 就是本文所建立的銀川地區(qū)加權(quán)平均溫度模型。

表2 銀川地區(qū)加權(quán)平均溫度模型

3.3 銀川地區(qū)加權(quán)平均溫度模型精度檢驗(yàn)

為了驗(yàn)證建立的銀川地區(qū)加權(quán)平均溫度模型的精度，統(tǒng)計(jì)了由探空資料計(jì)算得到的加權(quán)平均溫度真值，與模型擬合值之間的殘差分布情況，如圖3 所示。由圖3 可以看出：殘差值在-10～10 mm所占比例為92.36 %；殘差值在-20～-15，-15～-10與10～15 mm 所占的比例分別為0.43 %，4.46 %與2.75 %。這些數(shù)據(jù)可以得出，本文所建立的銀川地區(qū)加權(quán)平均溫度模型較為可靠。

圖3 擬合殘差值分布

結(jié)合已有模型，對(duì)本文所建的銀川地區(qū)加權(quán)平均溫度模型進(jìn)行精度檢驗(yàn)，分別計(jì)算平均偏差、平均相對(duì)誤差（mean absolute error, MAE）、均方根誤差（root mean square, RMS），其結(jié)果如表3所示。

表3 銀川地區(qū)加權(quán)平均溫度模型精度檢驗(yàn)

由表3 可以看出，銀川地區(qū)本地化模型所對(duì)應(yīng)的平均偏差、MAE、RMS 均小于Bevis 經(jīng)驗(yàn)公式、姚宜斌模型、龔紹琦模型。其中，姚宜斌模型與龔紹琦模型所對(duì)應(yīng)的值均小于Bevis 經(jīng)驗(yàn)公式。即本文所建立的銀川地區(qū)加權(quán)平均溫度模型，精度明顯高于已有模型，本地化模型適用性更強(qiáng)、精度更高，可以用于該區(qū)域PWV 的反演計(jì)算。

4 結(jié)束語(yǔ)

本文利用2008—2017 年10 a 間銀川地區(qū)的探空站數(shù)據(jù)，計(jì)算得到每個(gè)采樣點(diǎn)的加權(quán)平均溫度值，再結(jié)合探空站地面溫度值，建立了銀川地區(qū)加權(quán)平均溫度模型并對(duì)其精度進(jìn)行檢驗(yàn)，得出以下結(jié)論：

1）加權(quán)平均溫度普遍低于地面溫度，2 者之間的變化趨勢(shì)基本保持一致，其相關(guān)系數(shù)高達(dá)0.804 3，具有較強(qiáng)的線性相關(guān)性；

2）加權(quán)平均溫度真值與所建模型擬合值之間的殘差，在10 mm 以?xún)?nèi)所占比重為92.36 %，這表明所建立的模型與實(shí)際值相吻合；

3）銀川地區(qū)加權(quán)平均溫度模型，所對(duì)應(yīng)的平均偏差、平均相對(duì)誤差、均方差均小于已有模型，其精度更高、適用性更強(qiáng)，完全可以應(yīng)用于銀川地區(qū)大氣可降水量的反演過(guò)程中。

4）區(qū)域化加權(quán)平均溫度模型精度較高，但其僅僅適用于這個(gè)區(qū)域，因此具有一定的局限性。