鄒海波，單九生，吳珊珊，劉熙明

(1.江西省氣象科學(xué)研究所，江西南昌 330046;2.中山大學(xué)大氣科學(xué)系∥季風(fēng)與環(huán)境研究中心，廣東廣州 510275;3.江西省氣候中心，江西南昌 330046）

GPS衛(wèi)星發(fā)出的信號(hào)穿越大氣時(shí)，受大氣成分折射的影響，衛(wèi)星信號(hào)會(huì)發(fā)生延遲現(xiàn)象，根據(jù)這種延遲信號(hào)來(lái)測(cè)定大氣中的溫度和水汽含量、監(jiān)測(cè)氣候變化等叫做GPS氣象學(xué) (GPS/Meteorology，簡(jiǎn)寫(xiě)為GPS/MET)。20世紀(jì)90年代以來(lái)，隨著GPS/MET的不斷發(fā)展，GPS水汽反演技術(shù)已經(jīng)被廣泛地應(yīng)用至大氣水汽監(jiān)測(cè)業(yè)務(wù)應(yīng)用和研究中［1－8］，上海、北京、湖北、四川、廣東和江西等地先后利用GAMIT軟件和本地地基GPS觀測(cè)網(wǎng)建立了GPS大氣可降水量監(jiān)測(cè)業(yè)務(wù)應(yīng)用平臺(tái)。文獻(xiàn)［3－4，9］等先后發(fā)現(xiàn)華東、河北和江西等地區(qū)降水產(chǎn)生前，GPS大氣可降水量 (PW)都有一個(gè)急升的過(guò)程。

GPS衛(wèi)星信號(hào)穿越大氣時(shí)的延遲 (總延遲)信號(hào)可分為天頂靜力延遲 (ZHD)(即天頂干延遲)和天頂濕延遲 (ZWD)［10－11］。ZWD近似正比于沿信號(hào)傳播路徑上的PW，根據(jù)這一特點(diǎn)，Bevis等［10］還推導(dǎo)出了PW與ZWD的函數(shù)關(guān)系，在已知ZWD后可通過(guò)這一關(guān)系反算出PW，而ZWD可通過(guò)天頂總延遲 (ZTD)減ZHD獲得，其中ZTD可由GAMIT軟件直接計(jì)算得到，ZHD則要通過(guò)Saastamoinen 模 型［12］、Hopfield 模 型［13］或 Black 模型［14］和地面氣壓計(jì)算得到，而這些模型中的系數(shù)也會(huì)隨地理位置的變化而變化。徐桂榮等［15］先后利用宜昌、武漢、恩施等探空站資料建立了適合湖北地區(qū)的天頂靜力延遲模型，有效地提升了湖北地區(qū)GPS水汽的反演精度。為了提升江西省GPS水汽的反演精度，有必要對(duì)天頂靜力延遲模型進(jìn)行改進(jìn)，為提升短時(shí)臨近預(yù)報(bào)水平和改善數(shù)值模式初始場(chǎng)服務(wù)。

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

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

江西省氣象局與江西省測(cè)繪局已經(jīng)合作建成了62個(gè)布設(shè)在氣象臺(tái)站的地基GPS觀測(cè)站，這些觀測(cè)站自2009年6月開(kāi)始正常觀測(cè)。南昌和贛州兩地不僅有GPS觀測(cè)站，還有探空觀測(cè)站，這不僅有利于建立適合江西的ZHD模型，還有利于分析GPS PW的反演精度。因此，本文選用南昌和贛州2個(gè)GPS觀測(cè)站的數(shù)據(jù)和探空站數(shù)據(jù)作分析樣本，探空數(shù)據(jù) (每天2次)和GPS/MET數(shù)據(jù)的時(shí)間長(zhǎng)度均為2010年1月1日至12月31日。此外，為了提高GAMIT軟件的解算精度，本文還加入了上海、昆明、武漢和臺(tái)北4個(gè)長(zhǎng)基線GPS觀測(cè)站數(shù)據(jù)，安福、德安、崇仁、德興、奉新、高安、吉安、南城、景德鎮(zhèn)和萍鄉(xiāng)等10個(gè)本地GPS觀測(cè)站的數(shù)據(jù)。用于分析臺(tái)風(fēng)的路徑和中心氣壓資料來(lái)源于日本氣象廳提供的最佳路徑資料 (Best Track Data)。用于水汽收支診斷分析的大尺度分析資料來(lái)源NASA提供的MERRA MAI3CPASM資料，其水平分辨率1.25°×1.25°，時(shí)間分辨率為每次3 h。

1.2 處理方法

當(dāng)GPS衛(wèi)星發(fā)出信號(hào)穿越大氣時(shí)，受大氣折射的影響而產(chǎn)生延遲，其中大氣折射數(shù)N與大氣溫濕的關(guān)系可表示為:

其中，T為絕對(duì)溫度 (單位K)，Rd為干空氣氣體常數(shù)，ρ為大氣總密度 (單位g/m3)，e為水汽壓(單位 hPa)，k1、k2、k3為常數(shù)。Thayer［16］在 1974年給出了 k1=77.604，k2=64.79，k3=3.776×105K/hPa。2002 年 Rueger［17］認(rèn)為 (1)式中的 3個(gè)常數(shù)應(yīng)采用可 k1=77.6890，k2=71.2952，k3=375463 K/hPa。2009 年徐桂榮［15］在改善湖北地區(qū)天頂靜力延遲模型時(shí)便采用了Rueger給出的常數(shù)，并得到了理想的結(jié)果。因此，對(duì)于公式 (1)中的3個(gè)常數(shù)本文也選用2002年Rueger給出的數(shù)值。

公式 (1)中右邊第一項(xiàng)稱(chēng)為干項(xiàng)Nd=k1Rdρ，對(duì)干項(xiàng)從地面積分至大氣層頂便可得ZHD。

其中，H為測(cè)站海拔高度，z為測(cè)站上空的垂直高度 (m)。(1)式中右邊第二、三項(xiàng)合稱(chēng)為干項(xiàng)濕項(xiàng)Nw=k2e/T+k3e/T2，對(duì)濕項(xiàng)從地面積分至大氣層頂便可得ZWD，即:

其中，水汽壓e(即露點(diǎn)溫度的飽和水汽壓)通過(guò)1966年世界氣象組織建議采用的飽和水汽壓計(jì)算計(jì)算公式計(jì)算獲得，即:

其中，T1=273.16 K為水的三相溫度，Td是絕對(duì)露點(diǎn)溫度 (K)，即露點(diǎn)溫度的絕對(duì)溫度。

ZHD可以通過(guò) (2)式積分直接得到，但在平流層及其以上仍有空氣存在，即空氣密度ρ不等于0 g/m3，而目前常用的氣象探空資料的最大高度通常為100 hPa(約16 km)，微波輻射計(jì)的探測(cè)高度約為10 km。因此，根據(jù) (2)式用常規(guī)氣象探空資料和微波輻射計(jì)資料積分不能準(zhǔn)確地求得ZHD。(3)式中的水汽壓e是由大氣中的水汽含量決定，大氣中的水汽主要集中于對(duì)流層中低層，隨著高度的增加水汽迅速減小，到了對(duì)流層中上層 (300 hPa，約9.2 km)，大氣中的水汽可以忽略不計(jì)，故利用常規(guī)探空資料或微波輻射計(jì)資料運(yùn)用 (3)式積分可獲得準(zhǔn)確的ZWD。徐桂榮［15］利用宜昌站的加密探空資料分別對(duì) (2)和 (3)式進(jìn)行了積分，其結(jié)果也表明從地面往上至30 km，ZHD隨高度變化一直較大，而ZWD在7 km以下隨高度變化較大，但在7 km以上隨高度變化非常小。因此，通過(guò)GAMIT軟件計(jì)算出ZTD，減去由 (3)式計(jì)算的ZWD并可以獲得可靠的ZHD，本文將以這個(gè)ZHD作為“觀測(cè)值”用以檢驗(yàn)天頂靜力延遲模型的精度。

2 ZHD本地化建模

利用2010年1月1日至2010年12月31日南昌和贛州兩個(gè)探空站的溫度和露點(diǎn)溫度廓線資料(1日2次，08時(shí)與20時(shí))，先運(yùn)用 (4)計(jì)算出這兩站的濕度廓線數(shù)據(jù)，再運(yùn)用 (3)式積分得到這兩站的ZWD。然后，運(yùn)用GAMIT軟件解算出2010年1月1日至12月31日的逐時(shí)ZTD，用這個(gè)ZTD減去 (3)式積分得到的ZWD作為ZHD的觀測(cè) 值 (即 探 空 ZHD)。 Saastamoinen［12］、Hopfield［13］和 Black［14］分別指出 ZHD 模型為:

其中，ZHDS代表 Saastamoinen模型，ZHDH代表Hopfield模型，ZHDB代表Black模型，Ps為地面氣壓 (單位hPa)，H為海拔高度 (單位km)，Ts為地面絕對(duì)溫度 (單位K)，φ為緯度。ZHD模型(5)－(7)是采用文獻(xiàn)［18］給出的折射率公式并在理想大氣和靜力平衡假定條件下推導(dǎo)出來(lái)的，不僅折射率公式本身存在偏差，而且大氣也不完全滿足靜力平衡，這給 (5)－(7)的系數(shù)帶來(lái)了不確定性，因而ZHD模型 (5)－(7)會(huì)有模型偏差。

圖1為探空Z(yǔ)HD與運(yùn)用 (5)－(7)式計(jì)算得到的ZHD之差ΔZHD(探空Z(yǔ)HD減天頂靜力延遲模型計(jì)算的ZHD)分布圖。圖1a和圖1b顯示整個(gè)2010年絕大多數(shù)ΔZHD都大于0 mm(空心圓點(diǎn)多位于0 mm線以上)，且最大差值接近150 mm，這表明ZHDS和ZHDH明顯較探空Z(yǔ)HD值偏小。從圖1c可以看出，探空Z(yǔ)HD與ZHDB的差ΔZHD明顯較前兩個(gè)模型 (圖1a和圖1b)小，其中最大差值接近100 mm，但在1－3月和11－12月ZHDB明顯較探空Z(yǔ)HD值偏大 (空心圓點(diǎn)多位于0 mm線以下)。為了進(jìn)一步分析 Saastamoinen模型、Hopfield模型和Black模型的模型誤差，我們將探空Z(yǔ)HD與各模型計(jì)算的ZHD值作了偏差 (平均值之差)和標(biāo)準(zhǔn)差分析 (見(jiàn)表1)。從表1可以看出探空Z(yǔ)HD與ZHDS的偏差最大，為30.0 mm，按照PW=0.15×ZWD 的公式［15］折算得 PW=4.7 mm，與ZHDH的偏差次之，為29.1 mm，折算成PW為4.4 mm，與ZHDB的偏差最小，為－5.6 mm，折算成PW為－0.9 mm。從表1還可以看出探空Z(yǔ)HD與ZHDS、ZHDH和ZHDB的標(biāo)準(zhǔn)差的差異不大，約為27.8 mm，折算成PW為4.2 mm。

以上分析證明了天頂靜力延遲模型 (5)－(7)在江西地區(qū)有著明顯的模型偏差，為了減小 (5)－(7)式的模型偏差，使其能更好地應(yīng)用至江西地區(qū)，本文將對(duì) (5)－(7)式進(jìn)行改進(jìn)，使 (5)－(7)的模型偏差降到最小。首先將天頂靜力延遲模型 (5)－(7)式中的系數(shù)改為待定系數(shù):

圖1 2010年江西探空Z(yǔ)HD與ZHD模型模擬的ZHD之差Fig.1 The ZHD difference between sounding ZHD and simulated ZHD in Jiangxi in 2010

表1 2010年江西地區(qū)探空Z(yǔ)HD與天頂靜力延遲模型計(jì)算的ZHD偏差和標(biāo)準(zhǔn)差對(duì)比Table 1 The bias and standard deviation of sounding ZHD and simulation ZHD in Jiangxi in 2010 mm

rs0、rs1、rh0、rh1、rb0和rb1為各模型中的系數(shù)，然后可利用探空Z(yǔ)HD，運(yùn)用最小二乘法對(duì) (8)－(10)式進(jìn)行回歸建模，便可求得待定系數(shù)rs0、rs1、rh0、rh1、rb0和rb1，將這些待定系數(shù)引入 (8)－(10)便可以得到適合江西地區(qū)的天頂靜力延遲模型。徐桂榮等［15］用該方法研究發(fā)現(xiàn)適合湖北地區(qū)的天頂靜力延遲模型中的系數(shù)為rs0=0.5576，rs1=0.0017，rh0=40.082，rh1=0.141898，rb0=0.5576和rb1=0.0017，認(rèn)為改善后的Saastamoinen模型、Hopfield模型和Black模型明顯提升了ZHD的計(jì)算精度。本文利用2010年南昌和贛州兩探空站資料計(jì)算的ZHD對(duì) (8)－(10)式進(jìn)行了回歸建模，并得到了適合江西地區(qū)的天頂靜力模型(改善后的天頂靜力延遲模型)為:

圖2為2010年探空Z(yǔ)HD與運(yùn)用改善后的天頂靜力延遲模型 (11)－(13)式模擬ZHD之差ΔZHD的分布圖。與圖1中ΔZHD的分布明顯不同，圖2中的ΔZHD在整個(gè)2010年都均勻地分布在0 mm線的兩側(cè)，相比圖1(改善前)，ΔZHD的值也明顯減小 (特別是圖2a和圖2b)，ΔZHD的最大值約為100 mm(圖1a和圖1b中約為150 mm)。這表明改善后的天頂靜力延遲模型 (11)－(13)有效地提升了ZHD的計(jì)算精度。另外，從表1可以看出，探空Z(yǔ)HD與ZHDS、ZHDH和ZHDB的偏差均為0 mm，標(biāo)準(zhǔn)差約為26.3 mm，折算成對(duì)PW的影響為約為3.9 mm。可見(jiàn)，改善后的天頂靜力延遲模型 (11)－(13)可以很好地?cái)M合探空Z(yǔ)HD，減小了原模型 (5)－(7)的系統(tǒng)偏差，并有效地提升了天頂靜力延遲模型的精度。

圖2 探空Z(yǔ)HD與改善后的ZHD模型模擬的ZHD之差Fig.2 The ZHD difference between sounding ZHD and simulated ZHD by improved ZHD models

3 GPS反演PW的精度分析

以上研究表明，運(yùn)用改善后的天頂靜力延遲模型 (11)－(13)計(jì)算的ZHD，明顯較改善前計(jì)算的ZHD的精度高，那么將改善后的天頂靜力延遲模型 (11)－(13)引入GAMIT軟件中能否解算出精度更高的PW?為了驗(yàn)證這一問(wèn)題，分別將改善前后的天頂靜力延遲模型 (5)－(7)式與 (11)－(13)式引入GAMIT軟件中，對(duì)2010年南昌和贛州兩站的PW進(jìn)行解算 (共6組試驗(yàn))，獲得了逐時(shí)GPS PW值，并與探空PW進(jìn)行對(duì)比分析，以確定采用改善后的天頂靜力延遲模型 (11)－(13)能否有效提升GPS PW的反演精度。

圖3為2010年江西地區(qū)探空PW與采用改善前后的Saastamoinen天頂靜力延遲模型解算的GPS PW之差ΔPW(探空PW－GPS PW)的分布圖。從圖3a可以看出，采用改善前的Saastamoinen模型所解算的GPS PW明顯較探空PW大 (ΔPW值多小于0 mm)，這主要是在總的ZTD不變的情況下，采用ZHDS偏小 (圖1a)，ZWD=ZTD－ZHD將偏大，最終使得PW偏大。與圖3a明顯不同，圖3b中ΔPW比較均勻地分布在0 mm線的兩側(cè)，這表明采用了改善后的Saastamoinen模型所解算的GPS PW更接近觀測(cè)值 (即探空PW)。與Saastamoinen模型類(lèi)似，采用改善前的Hopfield模型所解算的GPS PW明顯較探空PW大 (圖4a)，這主要也是由于改善前的ZHDH偏小所致，而采用改善后的Hopfield模型解算的GPS PW也更接近探空PW(圖4b)。與改善前的Saastamoinen模型和Hopfield模型不同，采用改善前的Black模型所解算的GPS PW明顯較探空PW小 (ΔPW值多大于0 mm)，這主要是由于改善前的ZHDB偏大所致，而采用改善后的Black模型所解算的GPS PW也更接近于探空PW(圖5b)。綜上所述，采用改善后的Saastamoinen模型、Hopfield模型和Black模型所解算的GPS PW都較改善前更接近觀測(cè)值。

表2給出了2010年江西地區(qū)探空PW與采用不同天頂靜力延遲模型解算的GPS PW的偏差、標(biāo)準(zhǔn)差和相關(guān)系數(shù)。從表2可以看出，采用改善前后的天頂靜力延遲模型解算得到的GPS PW與探空PW都具有很好的相關(guān)性，其相關(guān)系數(shù)都為0.97;相對(duì)改善前的天頂靜力延遲模型，采用改善后的天頂靜力延遲模型解算的GPS PW明顯減小了與探空PW的偏差，改善前Hopfield模型和Saastamoinen模型的偏差均為負(fù)值，分別為 －3.97和 －2.89 mm，Black模型的偏差為正值，其值為1.73 mm，這與圖3a、圖4a和圖5a中ΔPW的分布情況吻合。改善后Saastamoinen模型、Hopfield模型和Black模型的偏差均小于等于0.80 mm，其中偏差最小的為改善后的Hopfield模型，其值為0.37 mm，相比改善前的－3.97 mm有顯著的較小;采用改善前后的天頂靜力延遲模型解算得到的GPS PW與探空PW的均方差變化不大，其值在4.31～4.42 mm之間，但相比改善前的均方差，改善后的均方差均有所減小，減小范圍為0.3～0.9 mm，與偏差一樣減小最明顯的是Hopfield模型，改善后較改善前減小了0.9 mm。

圖3 2010年江西探空PW與采用Saastamoinen模型解算的GPS PW之差Fig.3 The PW difference between sounding PW and calculated PW using Saastamoinen model

表2 2010年江西探空PW與采用不同天頂靜力延遲模型解算的GPS PW的偏差、標(biāo)準(zhǔn)差和相關(guān)系數(shù)Table 2 The bias，standard deviation and correlation coefficient of sounding PW and calculated PW using different ZHD model

圖4 2010年江西探空PW與采用Hopfield模型解算的GPS PW之差Fig.4 The PW difference between sounding PW and calculated PW using Hopfield model

圖5 2010年江西探空PW與采用Black模型解算的GPS PW之差Fig.5 The PW difference between sounding PW and calculated PW using Black model

以上分析表明，相比采用改善前的Saastamoinen模型、Hopfield模型和Black模型解算的GPS PW，采用改善后的Saastamoinen模型、Hopfield模型和Black模型解算的GPS PW更接近探空PW(圖3b、圖4b和圖5b)，且與探空PW相比，偏差都明顯減小，標(biāo)準(zhǔn)差也有所較小，相關(guān)系數(shù)不變。其中，采用改善后的Hopfield模型所解算的GPS PW與探空PW的偏差最小，為0.37 mm，這明顯優(yōu)于湖北地區(qū)［16］天頂靜力延遲模型改善后的效果，其改善后的最小偏差為1.2 mm。

4 GPS PW在臺(tái)風(fēng)暴雨中的應(yīng)用分析

大氣中的水汽不僅是臺(tái)風(fēng)產(chǎn)生降水的源泉，還是臺(tái)風(fēng)維持系統(tǒng)所需能量源的主要成分，因此對(duì)高時(shí)空分辨率的GPS PW資料的分析，有利于提升對(duì)臺(tái)風(fēng)降水的預(yù)報(bào)。Lion［19］和丁金才［20］研究都發(fā)現(xiàn)了在臺(tái)風(fēng)影響期間，GPS PW的變化與降水有較好的對(duì)應(yīng)關(guān)系，即臺(tái)風(fēng)降水產(chǎn)生前GPS PW都有一個(gè)急升過(guò)程，GPS PW劇降后降水即將結(jié)束。本節(jié)將利用改善的Hopfield模型解算的GPS PW資料，對(duì)2012年8月臺(tái)風(fēng)“?？?HAIKUI)影響江西期間進(jìn)行分析，以期在今后臺(tái)風(fēng)降水預(yù)報(bào)中使用GPS PW資料提供參考。

4.1 臺(tái)風(fēng)“?？钡母艣r

利用日本氣象廳提供的最佳臺(tái)風(fēng)路徑資料，我們制作了臺(tái)風(fēng)“?？钡穆窂綀D (圖6)。從圖6可以看出，2012年第11號(hào)臺(tái)風(fēng)“?？庇?月1日20時(shí)在日本沖繩縣東偏南方約1850 km的西北太平洋洋面發(fā)展為熱帶低壓，中心氣壓為1004 hPa，隨后逐漸加強(qiáng)西移，3日08時(shí)在沖繩縣東偏南方約1360 km的洋上加強(qiáng)為熱帶風(fēng)暴 (中心氣壓為992 hPa)，并轉(zhuǎn)向西偏北方向移動(dòng)，5日17時(shí)進(jìn)入我國(guó)東海東部海面 (中心氣壓為980 hPa)，5日20時(shí)加強(qiáng)為強(qiáng)熱帶風(fēng)暴，中心氣壓為975 hPa。此后“?？崩^續(xù)加強(qiáng)并先西移再往西北方向移動(dòng)，7日20時(shí)在距離浙江象山東南面大約180 km的東海洋面上加強(qiáng)為臺(tái)風(fēng)，中心氣壓降至965 hPa，此時(shí)“?？钡耐鈬菪h(huán)流開(kāi)始影響江西東北地區(qū) (圖略)，8日03時(shí)在浙江象山登陸，登陸后“?？笨焖贉p弱并繼續(xù)向西北方向移動(dòng)，8日14時(shí)在浙江杭州附近減弱為熱帶風(fēng)暴，中心氣壓升至980 hPa，此后“?？本徛郎p弱并轉(zhuǎn)向西偏北方向移動(dòng)，9日08時(shí)，移至安徽黃山市 (中心氣壓為992 hPa)，此后一直穩(wěn)定維持在安徽西南部至贛皖邊界一帶，其強(qiáng)度減弱緩慢，9日20時(shí)減弱為熱帶低壓 (中心氣壓為994 hPa)，11日08時(shí)“?？钡闹行臍鈮荷?000 hPa，11日14時(shí)后在贛皖邊界中部附近消亡。

圖6 2012年8月11號(hào)臺(tái)風(fēng)“?？钡囊苿?dòng)路徑Fig.6 The moving path of typhoon HAIKUI

受臺(tái)風(fēng)“?？钡挠绊懀?月9－10日江西東北部出現(xiàn)了大暴雨過(guò)程，先后有19個(gè)縣 (市)出現(xiàn)大暴雨，4個(gè)縣 (市)出現(xiàn)特大暴雨，其中日降水量最大出現(xiàn)在8月9日景德鎮(zhèn)的浮梁，為397.2 mm。臺(tái)風(fēng)“?？睅?lái)的暴雨致使昌江、樂(lè)安河水位一度全線超警戒，江西東北部地區(qū)出現(xiàn)洪澇、內(nèi)澇和山洪等災(zāi)害，這對(duì)農(nóng)業(yè)、水利、交通等造成了嚴(yán)重的影響。這次臺(tái)風(fēng)過(guò)程是江西歷史上致災(zāi)最嚴(yán)重的臺(tái)風(fēng)之一。據(jù)統(tǒng)計(jì)，臺(tái)風(fēng)“?？敝率咕暗骆?zhèn)、九江、上饒、鷹潭、撫州和南昌6市192.9萬(wàn)人受災(zāi)，緊急轉(zhuǎn)移安置26.6萬(wàn)人，農(nóng)作物受災(zāi)面積12.7萬(wàn)hm2，江西直接經(jīng)濟(jì)損失達(dá)32.9億元，其中景德鎮(zhèn)的災(zāi)害最為嚴(yán)重。

4.2 GPS PW與降水的關(guān)系

圖7 2012年8月7日12時(shí)至11日23時(shí)GPS PW和逐時(shí)降水量的演變過(guò)程Fig.7 The evolutions of GPS PW and rainfall during 1200 BTC 7 to 2300 BTC，August 11th，2012

圖7給出了臺(tái)風(fēng)“?？庇绊懡髌陂g，江西北部景德鎮(zhèn)、鄱陽(yáng)和南昌三站GPS PW與逐時(shí)降水量的演變過(guò)程。從圖7可以看出，在臺(tái)風(fēng)影響前(8月7日18時(shí)前)景德鎮(zhèn)、鄱陽(yáng)和南昌的GPS PW均小于45 mm，無(wú)降水出現(xiàn);8月7日18時(shí)后，隨著臺(tái)風(fēng)“?？睆恼憬瓥|面洋面不斷往西偏北方向移動(dòng) (圖6)，“?？遍_(kāi)始自東向西影響江西地區(qū)，GPS PW也隨之自東向西開(kāi)始明顯增加，位置偏東的景德鎮(zhèn)于7日19時(shí)開(kāi)始明顯增加，鄱陽(yáng)在7日21時(shí)開(kāi)始明顯增加，而位置偏西的南昌則在7日22時(shí)才開(kāi)始明顯增加，8日13時(shí)景德鎮(zhèn)的GPS PW率先超過(guò)60 mm，并在隨后2 h出現(xiàn)明顯降水 (1小時(shí)降水量超過(guò)1 mm)，8日18時(shí)鄱陽(yáng)的GPS PW超過(guò)60 mm，在隨后2 h也出現(xiàn)明顯降水，而南昌的GPS PW則在9日2時(shí)才超過(guò)60 mm，但同樣在隨后2 h出現(xiàn)了明顯降水;隨著“?？毕虬不漳喜恐邻M皖邊界附近的移動(dòng)，景德鎮(zhèn)、鄱陽(yáng)和南昌三站離臺(tái)風(fēng)中心的距離持續(xù)減小，GPS PW繼續(xù)增加，到10日05時(shí)這三站的GPS PW都到達(dá)了最大值，其中景德鎮(zhèn)為77.6 mm，鄱陽(yáng)為79.7 mm，南昌為76.2 mm，與之對(duì)應(yīng)景德鎮(zhèn)在10日5時(shí)前后出現(xiàn)了超過(guò)32 mm/h的強(qiáng)降水(圖7a)，鄱陽(yáng)湖和南昌也出現(xiàn)了超過(guò)12 mm/h的強(qiáng)降水 (圖7b和圖7c);此后，雖然“海葵”的中心位置維持在贛皖邊界附近，但其強(qiáng)度持續(xù)減弱，GPS PW也隨之持續(xù)減小，10日22時(shí)距離臺(tái)風(fēng)中心最遠(yuǎn)的南昌站的PGS PW率先降至60 mm以下，且降水早在GPS PW降至62 mm時(shí) (10日18時(shí))就結(jié)束了，11日08時(shí)隨著“海葵”的中心氣壓升至1000 hPa，景德鎮(zhèn)、鄱陽(yáng)和南昌三站GPS PW都有一個(gè)劇降過(guò)程，景德鎮(zhèn)和鄱陽(yáng)的GPS PW也在這次劇降過(guò)程中降至60 mm以下，降水也在隨后結(jié)束。從圖7還可以看出，在臺(tái)風(fēng)“海葵”影響江西期間，雖然景德鎮(zhèn)的降水明顯較鄱陽(yáng)的降水強(qiáng)，但鄱陽(yáng)的最大GPS PW卻比景德鎮(zhèn)的大，造成這一現(xiàn)象的主要原因有兩個(gè)，一是在臺(tái)風(fēng)南側(cè)偏西風(fēng)氣流下，景德鎮(zhèn)位于懷玉山脈的迎風(fēng)坡，迎風(fēng)坡的輻合抬升有利于降水的發(fā)生，二是景德鎮(zhèn)的海拔較鄱陽(yáng)高，PW是大氣的整層含水量，在相同環(huán)流場(chǎng)下高海拔地區(qū)的PW不如低海拔地區(qū)高。

4.3 PW的診斷分析

其中P為降水量，E為蒸發(fā)量 (雨天近似為0.0 mm)，q為比濕，vh為水平風(fēng)，t為時(shí)間，重力加速度g=9.8 m/s2。利用分部積分法，(14)式可寫(xiě)為:

(15)式中右邊第一項(xiàng)由水平風(fēng)的散度造成，稱(chēng)為散度項(xiàng)，第二項(xiàng)由比濕的平流造成，稱(chēng)為平流項(xiàng)。根據(jù) (15)式，利用時(shí)間分辨率為3 h每次、高空間分辨率為1.25°×1.25°經(jīng)緯網(wǎng)格的MERRA MAI3CPASM分析資料，對(duì)“?？庇绊懡髌陂g作診斷分析，以便找出引起PW變化的主要原因。

雖然MERRA資料的空間分辨率相對(duì)較低，但該資料在景德鎮(zhèn)東偏北方向約30 km處 (116.875°E，29.375°N)有一網(wǎng)格點(diǎn) (命名為A)，故選用該點(diǎn)的資料近似代替景德鎮(zhèn)的資料作PW的診斷分析，而鄱陽(yáng)和南昌站離MERRA資料的網(wǎng)格點(diǎn)較遠(yuǎn)，本文將不作診斷分析。圖8a為格點(diǎn)A的PW(用比濕資料計(jì)算所得)在臺(tái)風(fēng)“海葵”影響江西期間的演變過(guò)程，對(duì)比圖8a和圖7a可以發(fā)現(xiàn)格點(diǎn)A的PW與景德鎮(zhèn)的GPS PW的演變非常類(lèi)似:8月9日05時(shí)前都為顯著的上升過(guò)程，其中在7日02時(shí)前后PW的增加都有一個(gè)減弱趨勢(shì);9日06時(shí)至10日17時(shí)都為PW的高值時(shí)段，且PW在高值時(shí)段的波動(dòng)也非常相似，只是格點(diǎn)A的PW的最大值比景德鎮(zhèn)的GPS PW最大值稍小;10日17后都為PW的減弱過(guò)程，且在11日05至08時(shí)都有一個(gè)明顯的劇降過(guò)程。這表明用格點(diǎn)A的資料來(lái)近似診斷景德鎮(zhèn)的PW變化的機(jī)理是可行的。

根據(jù) (15)式右端的第一項(xiàng) (散度項(xiàng))和第二項(xiàng) (平流項(xiàng))，利用格點(diǎn)A的資料，計(jì)算出了臺(tái)風(fēng)“海葵”影響江西期間，平流項(xiàng)和散度項(xiàng)對(duì)PW的貢獻(xiàn)圖 (圖8b)。從圖8b可以看出，8日14時(shí)前景德鎮(zhèn)PW的增加主要是由平流項(xiàng) (強(qiáng)的比濕平流)貢獻(xiàn)所致，這主要是由于臺(tái)風(fēng)外圍環(huán)流將濕空氣輸送至景德鎮(zhèn)附近所致;隨著臺(tái)風(fēng)“?？辈粩嗫拷?，平流項(xiàng)的貢獻(xiàn)不斷減小，最后位于0 mm/h線附近波動(dòng)，這主要是由于比濕梯度減小所致，而散度項(xiàng)則取代平流項(xiàng)成為景德鎮(zhèn)PW增加或維持的主要貢獻(xiàn)者，這主要是隨著臺(tái)風(fēng)的靠近，臺(tái)風(fēng)中心的強(qiáng)烈輻合區(qū)控制了景德鎮(zhèn)地區(qū)所致;10日17時(shí)后散度項(xiàng)的貢獻(xiàn)也快速減弱至0 mm/h線附近，這主要是由于臺(tái)風(fēng)減弱所伴隨中心區(qū)輻合減弱所致。以上分析表明在臺(tái)風(fēng)影響初期平流項(xiàng)是PW增加的主要原因，隨后散度項(xiàng)代替平流項(xiàng)成為PW增加或維持的主要原因，而后期平流項(xiàng)和散度項(xiàng)的貢獻(xiàn)都很小，那么后期PW減小是什么原因造成的?從 (15)式還可以看出降水 (P＞0)會(huì)使空氣中的PW減小，而景德鎮(zhèn)在9日8時(shí)至10日18時(shí)都有明顯的降水，且在10日6時(shí)前后的降水強(qiáng)度非常大，很明顯10日05時(shí)后景德鎮(zhèn)PW的減小主要是由降水造成，雖然10日18時(shí)后景德鎮(zhèn)的降水明顯減小，但同樣10日17時(shí)后散度項(xiàng)對(duì)PW的正貢獻(xiàn)也快速減小為0 mm/h，故降水仍會(huì)使PW繼續(xù)減小，11日9是后景德鎮(zhèn)的降水基本結(jié)束，PW也維持在58 mm附近波動(dòng) (圖7a)。

圖8 2012年8月7日14時(shí)至11日23時(shí)MERRA資料計(jì)算的PW、散度項(xiàng)和平流項(xiàng)的演變過(guò)程Fig.8 The evolutions of PW，divergence term and advection term calculated using MERRA data during 1400 BTC，Aug 7th to 2300 BTC，Aug.11th，2012

5 討論與結(jié)論

利用2010年1月1日至12月31日江西省南昌、贛州兩個(gè)探空站每日2次 (00時(shí)和20時(shí))的探空資料計(jì)算出 ZHD作為“觀測(cè)值”，與利用Saastamoinen、Hopfield和Black天頂靜力延遲模型所計(jì)算的ZHD進(jìn)行對(duì)比分析發(fā)現(xiàn)，這三個(gè)天頂靜力延遲模型在江西地區(qū)計(jì)算的ZHD與探空Z(yǔ)HD有明顯的偏差，為了減小這些偏差，我們采用最小二乘法對(duì)這些模型進(jìn)行了改善。改善后的ZHDS、ZHDH和ZHDB不但與探空Z(yǔ)HD的偏差為0 mm，且與改善前的ZHDS、ZHDH和ZHDB相比，還減小了與探空Z(yǔ)HD的標(biāo)準(zhǔn)差。為了驗(yàn)證改善后的Saastamoinen、Hopfield和Black模型能否有效提升江西地區(qū)GPS PW的解算精度，我們利用GAMIT軟件分別采用改善前后的 Saastamoinen、Hopfield和Black模型 (共6組試驗(yàn))對(duì)2010年南昌和贛州的GPS PW進(jìn)行了解算，并與探空PW對(duì)比分析發(fā)現(xiàn)，采用改善后的Saastamoinen、Hopfield和Black模型所解算的PW不僅減小了與探空PW的偏差，還減小了與探空PW的標(biāo)準(zhǔn)差，其中改善后的Hopfield模型解算的GPS PW與探空PW的偏差最小，為0.37 mm。

利用改善后的Hopfield模型解算的GPS PW，對(duì)2012年8月7日至11日臺(tái)風(fēng)“?？庇绊懡鳀|北部期間作了分析，結(jié)果發(fā)現(xiàn)GPS PW與臺(tái)風(fēng)降水密切相關(guān):臺(tái)風(fēng)影響景德鎮(zhèn)、鄱陽(yáng)和南昌前其GPS PW均小于45 mm，都未出現(xiàn)降水;隨著“?？钡牟粩嗫拷暗骆?zhèn)、鄱陽(yáng)和南昌的GPS PW逐漸增加，且在GPS PW超過(guò)60 mm后2小時(shí)均出現(xiàn)了明顯的降水;隨著“?？钡倪M(jìn)一步靠近，10日05時(shí)景德鎮(zhèn)、鄱陽(yáng)和南昌的GPS PW都達(dá)到了最大值 (均超過(guò)了76 mm)，與之對(duì)應(yīng)強(qiáng)降水都出現(xiàn)在10日05時(shí)前后，即GPS PW的高值段;隨著“?？钡臏p弱消亡，景德鎮(zhèn)、鄱陽(yáng)和南昌的GPS PW逐漸減弱，當(dāng)GPS PW降至60 mm以下時(shí)，降水隨之結(jié)束。

本文還利用Trenberth推導(dǎo)的水汽收支方程和MERRA MAI3CPASM分析資料，對(duì)2012年8月7日至11日臺(tái)風(fēng)“?？庇绊懢暗骆?zhèn)期間GPS PW的變化作了診斷分析，其結(jié)果表明:在臺(tái)風(fēng)“?？庇绊懢暗骆?zhèn)初期，比濕平流過(guò)程是景德鎮(zhèn)GPS PW增加的主要原因;隨著“海葵”的靠近，比濕平流的作用逐漸減小，而水平風(fēng)輻合則轉(zhuǎn)變?yōu)镚PS PW增加的主要原因;在臺(tái)風(fēng)“海葵”影響后期，景德鎮(zhèn)GPS PW的減小主要是由10日17時(shí)后水平風(fēng)輻合的快速減弱和降水過(guò)程消耗大氣中的水汽所造成。

［1］俞炳，周媛，聞春華.江西省GPS/MET水汽監(jiān)測(cè)系統(tǒng)設(shè)計(jì)［J］.氣象與減災(zāi)研究，2010，33(2):60－63.

［2］鄒海波，單九生，吳珊珊，等.利用GAMIT對(duì)江西省GPS可降水量的反演應(yīng)用［J］.氣象與減災(zāi)研究，2010，33(3):56 －60.

［3］丁金才，黃炎，葉其欣，等.2002年臺(tái)風(fēng)Ramasun影響華東沿海期間可降水量的GPS觀測(cè)和分析［J］.大氣科學(xué)，2004，28(4):613 －624.

［4］陳小雷，景華，仝美然，等.地基GPS遙測(cè)大氣可降水量在天氣分析診斷中的應(yīng)用［J］.氣象，2007，33(6):19－24.

［5］姚建群，丁金才，王堅(jiān)捍，等.用GPS可降水量資料對(duì)一次大暴雨過(guò)程的分析［J］.氣象，2005，31(4):48－52.

［6］楚艷麗，郭英華，張朝林，等.地基GPS水汽資料在北京“7·10”暴雨過(guò)程研究中的應(yīng)用［J］，氣象，2007，33(12):16－22.

［7］李國(guó)翠，李國(guó)平，連志鸞，等.不同云系降水過(guò)程中GPS可降水量的特征——華北地區(qū)典型個(gè)例分析［J］，高原氣象，2008，27(5):1066 －1073.

［8］張志全，王文旭，黃立人，等.天津市GPS連續(xù)運(yùn)行參考站網(wǎng)系統(tǒng)［J］.城市勘測(cè)，2007，(6):51－55.

［9］廖敏，鄒海波，單九生，等.地基GPS水汽資料在江西一次大暴雨過(guò)程的應(yīng)用分析［J］.氣象與減災(zāi)研究，2012，35(3):33－38

［10］BEVIS M，BUSINGER S，HERRING T，et al.GPS meteorology:Remote sensing of atmospheric water vapor using the global positioning system［J］.J Geophys Res，1992，97(D14):15787 －15801.

［11］DAVIS J L，HERRING T A，SHAPIRO I I，et al.Geodesy of radio interferometry:Effects of atmospheric modeling errors on estimates of baseline length［J］.Radio Sci，1985，20(6):1593 －1607.

［12］SAASTAMOINEN J.Contributions to the theory 0f atmospheric refraction［J］.Bull Geod，1973，107(1):13 －34.

［13］HOPFIELD H S.Tropospheric effect on electromagnetically measured range:Prediction from surface weather data［J］.Radio Sci，1971，6(3):357 －367.

［14］BLACK H D.An easily implemented algorithm for the tropospheric range correction［J］.J Geophys Res，1978，83(B4):1825－1828.

［15］徐桂榮，萬(wàn)蓉，李武階，等.地基GPS反演大氣可降水量方法的改進(jìn)［J］.暴雨災(zāi)害，2009，28(3):203－209.

［16］THAYER G D.An improved equation for the radio refractive index of air［J］.Radio Sci，1974，9(10):803 －807.

［17］RUEGER J M.Refractive index formulae for radio wave［C］//FIG XXII International Congress.Washington D C，USA，2002.

［18］SMITH E K，WEINTRAUB S.The constants in the equation for atmospheric refractive index at radio frequencies［J］.Proc IRE，1953，41(8):1035 －1037.

［19］LIOU Y1 A1，HUANG C Y.GPS observation of PW during the passage of a typhoon［J］.Earth Planet Space，2000，52(10):709 －712.

［20］丁金才，黃炎，葉其欣，等.2002年臺(tái)風(fēng)Ramasun影響華東沿海期間可降水量的GPS觀測(cè)和分析［J］.大氣科學(xué)，2004，28(4):613 －624

［21］TRENBERTH K E，UILLEMOT C J.Evaluation of the global atmospheric moisture budget as seen from analyses［J］.J Climate，1995，8(9):2255 －2272.