GPS/PWV在臺(tái)風(fēng)“瑪莉亞”襲閩期間的變化特征研究

(福建省測(cè)繪院，福建 福州 350003)

可降水量(PWV)反映大氣中水汽轉(zhuǎn)化為降水的潛在動(dòng)能，對(duì)研究水汽的時(shí)空分布與變化特性具有重要價(jià)值。常規(guī)PWV監(jiān)測(cè)手段主要利用探空氣球每天進(jìn)行早晚兩次觀測(cè)，雖然精度較高，但由于時(shí)空分辨率過(guò)低和成本過(guò)高，很難得到大范圍的推廣普及。自從文獻(xiàn)[1]開創(chuàng)了GPS氣象學(xué)以來(lái)，利用高時(shí)空分辨率的地基GPS監(jiān)測(cè)PWV的方法獲得越來(lái)越廣泛的應(yīng)用，為氣象學(xué)研究和天氣預(yù)報(bào)提供了基礎(chǔ)資料。例如，文獻(xiàn)[2]率先在國(guó)內(nèi)開展GPS氣象學(xué)研究，獲得了可靠的PWV結(jié)果；文獻(xiàn)[3]研究發(fā)現(xiàn)，GPS/PWV能夠用于改進(jìn)數(shù)值預(yù)報(bào)模式。眾多研究表明，GPS/PWV的變化與降水之間具有很強(qiáng)的對(duì)應(yīng)關(guān)系[4-6],與霧霾之間具有顯著的相關(guān)特性[7-9]，為空間環(huán)境學(xué)的研究提供了良好的契機(jī)[10]。

臺(tái)風(fēng)是一種強(qiáng)烈的熱帶氣旋，在福建夏季多發(fā)，具有較強(qiáng)的破壞性，然而，由于臺(tái)風(fēng)能夠帶來(lái)豐富的水汽，臺(tái)風(fēng)期間對(duì)PWV的監(jiān)測(cè)有助于深入認(rèn)識(shí)大氣中水汽的傳輸過(guò)程。利用地基GPS資料已被證實(shí)能夠探測(cè)到臺(tái)風(fēng)天氣高精度的PWV[11-12],而對(duì)GPS/PWV資料分析有助于直觀認(rèn)識(shí)臺(tái)風(fēng)期間水汽的分布及動(dòng)態(tài)變化過(guò)程[13-14]。本文利用FJCORS基準(zhǔn)站GPS資料反演PWV，在檢驗(yàn)反演的精度基礎(chǔ)上，分析2018年第8號(hào)臺(tái)風(fēng)“瑪莉亞”登陸福建期間的GPS/PWV時(shí)空變化及其與降水量的關(guān)系，討論臺(tái)風(fēng)中水汽的交換和輸送過(guò)程，從而進(jìn)一步檢驗(yàn)GPS探測(cè)水汽的性能。

1 GPS/PWV資料分析

1.1 臺(tái)風(fēng)回顧

2018年第8號(hào)臺(tái)風(fēng)“瑪莉亞”于北京時(shí)間7月11日09：10(世界協(xié)調(diào)時(shí)間UTC 01:10)在福建連江黃岐半島登陸，沿西北西方向深入內(nèi)陸，UTC 04:00經(jīng)過(guò)閩清縣，UTC 08:00經(jīng)過(guò)順昌縣，UTC 12:00左右進(jìn)入江西省境內(nèi)。圖1展示了“瑪莉亞”的移動(dòng)路徑，來(lái)源于中央氣象臺(tái)·臺(tái)風(fēng)網(wǎng)，詳見網(wǎng)址：http:∥typhoon.nmc.cn/web.html。“瑪莉亞”給福建全省帶來(lái)強(qiáng)降水，其中，東北部和中北部出現(xiàn)大雨或暴雨，局地特大暴雨，累計(jì)降水量超過(guò)250 mm。據(jù)不完全統(tǒng)計(jì)，“瑪莉亞”共造成福建23.8萬(wàn)人受災(zāi)，700余間房屋不同程度損壞，直接經(jīng)濟(jì)損失2.4億元，破壞性較為嚴(yán)重。

1.2 資料來(lái)源及處理方法

本文利用FJCORS基準(zhǔn)站反演GPS/PWV，并引入國(guó)內(nèi)及周邊12個(gè)IGS(International GNSS Service)跟蹤站資料進(jìn)行同步解算，以削弱區(qū)域網(wǎng)PWV反演結(jié)果的相關(guān)性。其中，IGS站包括：BJFS、CHAN、CUSV、 DAEJ、 HYDE、 LHAZ、PIMO、POL2、SHAO、TWTF、ULAB、URUM。數(shù)據(jù)下載地址為ftp:∥garner.ucsd.edu/。探空資料來(lái)自美國(guó)懷俄明州立大學(xué)官方網(wǎng)站，網(wǎng)址為http:∥weather.uwyo.edu/。降水資料來(lái)自中國(guó)氣象數(shù)據(jù)網(wǎng)，網(wǎng)址為http:∥data.cma.cn/。

GPS可降水量源于GPS定位的對(duì)流層延遲誤差，對(duì)GPS定位進(jìn)行逆過(guò)程反演便得到了GPS可降水量。本文采用高精度解算軟件GAMIT[15]處理GPS觀測(cè)資料，從中解算出天頂方向?qū)α鲗涌傃舆t,然后扣除穩(wěn)定的干空氣產(chǎn)生的流體靜力學(xué)延遲，得到天頂方向濕延遲，最后通過(guò)濕延遲公式轉(zhuǎn)換得到GPS/PWV。其中，天頂對(duì)流層總延遲的計(jì)算，是在利用Saastamoinen模型[16]進(jìn)行改正的基礎(chǔ)上，通過(guò)附加天頂延遲參數(shù)進(jìn)一步修正。主要計(jì)算公式為

ZTD=ZHD+ZWD

(1)

PWV=ZWD·κ/ρ

(2)

式中，ZTD為天頂對(duì)流層總延遲；ZHD、ZWD為天頂流體靜力學(xué)延遲和天頂濕延遲；ρ為液態(tài)水的密度；κ為一個(gè)常數(shù)，其計(jì)算公式詳見文獻(xiàn)[1]。

1.3 精度檢驗(yàn)

采用2018年1—8月GPS/PWV資料與高精度探空PWV資料進(jìn)行對(duì)比，以檢驗(yàn)水汽反演的精度。其中，GPS/PWV為逐小時(shí)值，探空PWV為每天UTC 0:00和UTC 12:00兩次觀測(cè)值，將同時(shí)刻數(shù)據(jù)提取進(jìn)行分析，以福州和廈門站為例繪制了散點(diǎn)對(duì)比圖，如圖2所示。圖中顯示，散點(diǎn)均勻地分布在直線y=x兩側(cè)，進(jìn)一步的統(tǒng)計(jì)結(jié)果表明，兩站的相關(guān)系數(shù)均達(dá)到0.98，呈高度相關(guān)。在精度方面，福州站平均偏差為1.10 mm，標(biāo)準(zhǔn)差為3.72 mm；廈門站平均偏差為0.67 mm，標(biāo)準(zhǔn)差為3.37 mm。如果忽略FJCORS基準(zhǔn)站和探空觀測(cè)站不同址(福州間距為5 km，高差為40 m；廈門間距為13 km，高差為80 m)而引起的微弱差異，則上述結(jié)果表明，GPS/PWV和探空PWV之間具有很好的一致性，而且精度相當(dāng)，能夠反映較為真實(shí)的水汽含量。

2 臺(tái)風(fēng)期間GPS/PWV空間動(dòng)態(tài)演變

在精度檢驗(yàn)的基礎(chǔ)上，本文采用2018年7月11日臺(tái)風(fēng)“瑪莉亞”登陸期間的GPS/PWV資料分析水汽傳輸過(guò)程。由于臺(tái)風(fēng)“瑪莉亞”登陸期間部分區(qū)域斷電，有完整數(shù)據(jù)資料的基準(zhǔn)站實(shí)際只有39個(gè)，利用這些站的逐小時(shí)資料內(nèi)插得到空間上連續(xù)分布的PWV，如圖3所示。從圖3可以清晰地看到，臺(tái)風(fēng)登陸前(臺(tái)風(fēng)于UTC 01:10登陸)，福建東北部和中北部區(qū)域?yàn)镻WV高值區(qū)，主要在60 mm以上，局地達(dá)到70 mm；而南部區(qū)域?yàn)镻WV低值區(qū)，主要在50 mm以下。臺(tái)風(fēng)登陸后，東北部和中北部PWV高值區(qū)不斷擴(kuò)大，同時(shí)，南部PWV低值區(qū)逐漸縮小；到UTC 06:00，西部已經(jīng)全變?yōu)镻WV高值區(qū)，并開始超過(guò)東部高值區(qū)；到UTC 09:00，西部PWV高值區(qū)已經(jīng)明顯超過(guò)東部高值區(qū)；UTC 12:00以后，西部PWV高值區(qū)逐漸向西北方向移動(dòng)，與此同時(shí)，東部PWV低值區(qū)逐漸擴(kuò)大，到UTC 18:00，東部PWV值主要在50 mm以下。上述PWV的變化梯度大致沿著連江-延平-將樂(lè)-建寧一帶展開，與臺(tái)風(fēng)“瑪莉亞”的移動(dòng)路徑(如圖1所示)具有較好的對(duì)應(yīng)關(guān)系，表明PWV的變化能夠較好地反映臺(tái)風(fēng)“瑪莉亞”登陸前后水汽的空間動(dòng)態(tài)傳輸過(guò)程。

3 臺(tái)風(fēng)期間GPS/PWV時(shí)變及與降水的對(duì)比

為進(jìn)一步研究臺(tái)風(fēng)“瑪莉亞”登陸前后水汽的時(shí)變特征，本文選取臺(tái)風(fēng)路徑上的連江、延平、將樂(lè)、建寧4個(gè)基準(zhǔn)站(站址分布如圖3所示)的GPS/PWV結(jié)果進(jìn)行分析。如圖4所示，其中時(shí)間以年積日表示，采用UTC時(shí)間，GPS可降水量為逐小時(shí)值，豎直虛線表示“瑪莉亞”登陸連江黃岐半島的時(shí)刻。

從圖4明顯看出，臺(tái)風(fēng)登陸前，各站PWV值總體下降，最低值為35～40 mm；隨著臺(tái)風(fēng)登陸，PWV值均迅速上升，達(dá)到峰值70 mm左右,其中，最大升速超過(guò)10 mm/h；在維持?jǐn)?shù)小時(shí)峰值后，PWV值快速下降，最大降速超過(guò)5 mm/h；臺(tái)風(fēng)經(jīng)過(guò)各站1天后，PWV降至谷值40 mm左右,然后呈現(xiàn)隨機(jī)波動(dòng)。從基準(zhǔn)站PWV值升降的時(shí)間順序來(lái)看，依次為連江、延平、將樂(lè)、建寧，與臺(tái)風(fēng)路徑非常吻合。

圖5顯示GPS可降水量與實(shí)際降水的對(duì)比,其中橫坐標(biāo)以年積日表示時(shí)間，左側(cè)縱軸表示逐小時(shí)GPS可降水量，右側(cè)縱軸表示逐小時(shí)降水量。從圖5可以看出，在臺(tái)風(fēng)登陸時(shí)，PWV的迅速上升伴隨著降水的急劇增加，PWV峰值前后出現(xiàn)降水的峰值，PWV的迅速下降伴隨著降水的急劇減少。上述結(jié)果較好地展示了臺(tái)風(fēng)期間水汽的交換和傳輸過(guò)程，表明PWV和降水之間具有極強(qiáng)的耦合性。

4 結(jié) 論

本文利用FJCORS基準(zhǔn)站GPS資料反演可降水量，并檢驗(yàn)其精度，在此基礎(chǔ)上分析臺(tái)風(fēng)“瑪莉亞”登陸期間水汽的變化特征，主要得到以下結(jié)論：

(1) GPS可降水量與探空觀測(cè)資料具有很好的一致性，其平均偏差不超過(guò)1.10 mm，標(biāo)準(zhǔn)偏差不超過(guò)3.72 mm，表明GPS可降水量具有較高的精度和良好的穩(wěn)定性，GPS具備良好的水汽探測(cè)性能。

(2) GPS/PWV的變化與臺(tái)風(fēng)過(guò)程具有很好的對(duì)應(yīng)關(guān)系，空間上的變化梯度與臺(tái)風(fēng)路徑一致，時(shí)間上的變化與臺(tái)風(fēng)過(guò)程吻合，同時(shí)與降水具有明顯的關(guān)聯(lián)性，表明GPS可降水量能夠較好地反映臺(tái)風(fēng)期間水汽的時(shí)空動(dòng)態(tài)傳輸過(guò)程。

致謝：感謝國(guó)家氣象信息中心共享氣象觀測(cè)資料，感謝美國(guó)懷俄明州立大學(xué)共享探空觀測(cè)資料，感謝美國(guó)麻省理工學(xué)院(MIT)授權(quán)使用GAMIT軟件包，感謝福建省測(cè)繪院應(yīng)急保障分院FJCORS中心在本文數(shù)據(jù)處理分析過(guò)程中給予的建議和幫助。