劉嚴(yán)萍 王 勇 聞德保

1)中南大學(xué)土木工程學(xué)院，長(zhǎng)沙 410075

2)天津城建大學(xué)地質(zhì)與測(cè)繪學(xué)院，天津 300384

3)長(zhǎng)沙理工大學(xué)交通運(yùn)輸工程學(xué)院，長(zhǎng)沙 410004

基于經(jīng)驗(yàn)?zāi)B(tài)分解的GPS可降水量變化研究*

劉嚴(yán)萍1)王 勇2)聞德保3)

1)中南大學(xué)土木工程學(xué)院，長(zhǎng)沙 410075

2)天津城建大學(xué)地質(zhì)與測(cè)繪學(xué)院，天津 300384

3)長(zhǎng)沙理工大學(xué)交通運(yùn)輸工程學(xué)院，長(zhǎng)沙 410004

選取中國(guó)地殼運(yùn)動(dòng)監(jiān)測(cè)網(wǎng)2000～2004年數(shù)據(jù)，使用高精度定位定軌軟件GAMIT處理，得到各測(cè)站可降水量時(shí)間序列，繼而采用經(jīng)驗(yàn)?zāi)B(tài)分解并提取各站可降水量變化趨勢(shì)。對(duì)處于不同氣候類型區(qū)域的GPS站可降水量進(jìn)行分析，結(jié)果顯示，處于熱帶季風(fēng)氣候、亞熱帶季風(fēng)氣候、溫帶大陸性氣候和青藏高原高寒氣候等區(qū)域的GPS站的可降水量變化趨勢(shì)為先升后降;處于溫帶季風(fēng)氣候區(qū)域的GPS站點(diǎn)的可降水量變化為先降后升。

GPS;經(jīng)驗(yàn)?zāi)B(tài)分解;可降水量;中國(guó)地殼運(yùn)動(dòng)監(jiān)測(cè)網(wǎng)絡(luò)

大氣中水汽含量的變化，對(duì)水平尺度100 km左右、生命史僅有幾小時(shí)的中小尺度災(zāi)害性天氣的監(jiān)測(cè)具有重要意義［1］。1992年 Bevis等［2］提出地基GPS的大氣水汽含量估算原理。國(guó)內(nèi)外在該領(lǐng)域出現(xiàn)一些代表性成果［3－9］，為GPS技術(shù)在氣象學(xué)領(lǐng)域的應(yīng)用奠定了基礎(chǔ)。中國(guó)大陸構(gòu)造環(huán)境監(jiān)測(cè)網(wǎng)絡(luò)于2012-03正式運(yùn)行，涵蓋260個(gè)連續(xù)觀測(cè)站，氣象學(xué)應(yīng)用是其重要應(yīng)用領(lǐng)域。隨著多個(gè)省級(jí)、市級(jí)GPS連續(xù)觀測(cè)網(wǎng)的建成運(yùn)行，歷史觀測(cè)數(shù)據(jù)逐步積累，從而使得長(zhǎng)時(shí)序的GPS數(shù)據(jù)在氣候領(lǐng)域的研究呈現(xiàn)廣闊研究空間。本文將針對(duì)2000～2004年的中國(guó)地殼運(yùn)動(dòng)監(jiān)測(cè)網(wǎng)絡(luò)數(shù)據(jù)，開(kāi)展短時(shí)氣候變化研究，以期為國(guó)內(nèi)GPS氣象學(xué)的氣候變化研究提供參考。

1 GPS對(duì)流層延遲解算與可靠性驗(yàn)證

1.1 對(duì)流層延遲解算處理

選取中國(guó)地殼運(yùn)動(dòng)監(jiān)測(cè)網(wǎng)絡(luò)2000～2004年的數(shù)據(jù)，使用GAMIT軟件處理。星歷采用IGS精密星歷，采樣間隔30 s，UTC時(shí)間00:00～24:00，截止高度角10°，模式為Relax，每小時(shí)估算一個(gè)對(duì)流層延遲值。對(duì)于中國(guó)地殼運(yùn)動(dòng)監(jiān)測(cè)網(wǎng)絡(luò)而言，測(cè)站間距大于500 km，對(duì)流層延遲為絕對(duì)值。GAMIT處理GPS數(shù)據(jù)后，提取5 a各測(cè)站的對(duì)流層延遲序列。

1.2 GPS結(jié)果可靠性驗(yàn)證

IGS提供的GPS站對(duì)流層延遲值具有較高精度。中國(guó)地殼運(yùn)動(dòng)監(jiān)測(cè)網(wǎng)絡(luò)有若干GPS站，本文以LHAS、WUHN兩站為驗(yàn)證樣本，與對(duì)應(yīng)的IGS提供的對(duì)流層延遲值比較，見(jiàn)圖1。圖中實(shí)線為IGS對(duì)流層延遲，三角形點(diǎn)線為作者利用中國(guó)地殼運(yùn)動(dòng)監(jiān)測(cè)網(wǎng)絡(luò)解算的對(duì)流層延遲。

由圖1可見(jiàn)，WUHN站對(duì)流層延遲與IGS對(duì)流層延遲非常吻合，LHAS站對(duì)流層延遲與IGS結(jié)果基本一致，證明本文解算的對(duì)流層延遲結(jié)果準(zhǔn)確可靠。

圖1 作者解算對(duì)流層延遲與IGS對(duì)流層延遲比較Fig.1 Comparison of tropospheric delay solved by authors with that by IGS

2 基于經(jīng)驗(yàn)?zāi)B(tài)分解的GPS可降水量趨勢(shì)項(xiàng)提取

經(jīng)驗(yàn)?zāi)B(tài)分解(empirical mode decomposition，簡(jiǎn)稱EMD)是Huang等［10］于1998年提出的，可對(duì)非平穩(wěn)數(shù)據(jù)進(jìn)行平穩(wěn)化處理，將復(fù)雜信號(hào)分解為有限本征模函數(shù)(IMF)和趨勢(shì)項(xiàng)，各IMF分量含原信號(hào)不同時(shí)間尺度的局部特征。經(jīng)驗(yàn)?zāi)B(tài)分解的理論及實(shí)現(xiàn)過(guò)程見(jiàn)文獻(xiàn)［11－12］。

在提取中國(guó)地殼運(yùn)動(dòng)監(jiān)測(cè)網(wǎng)絡(luò)各GPS站2000～2004年對(duì)流層延遲序列的基礎(chǔ)上，結(jié)合同期氣壓、溫度觀測(cè)值，經(jīng)Black模型處理，得到相應(yīng)站的可降水量序列［9］。部分 GPS站點(diǎn)因無(wú)相應(yīng)的氣象觀測(cè)數(shù)據(jù)，未得可降水量序列。對(duì)能獲可降水量序列的14個(gè)站，使用經(jīng)驗(yàn)?zāi)B(tài)分解并提取趨勢(shì)項(xiàng)，用于分析不同氣候類型的GPS可降水量變化情況。

3 不同氣候類型的GPS可降水量變化

QION處于熱帶季風(fēng)氣候區(qū)域。由圖2可見(jiàn)，該站的GPS可降水量變化趨勢(shì)為2000～2002年緩升，2002～2004年降(從42 mm降至39 mm)。

圖2 GPS可降水量變化趨勢(shì)(熱帶季風(fēng)氣候類型)Fig.2 The trend of GPS precipitable water vapor(tropical monsoon climate)

中國(guó)地殼運(yùn)動(dòng)監(jiān)測(cè)網(wǎng)絡(luò)的 WUHN，SHAO、KMIN、GUAN、XIAM位于亞熱帶季風(fēng)氣候類型區(qū)域［13］。由圖3可見(jiàn)，處于該氣候區(qū)域的大部分GPS站，除 KMIN站外，可降水量總體變化趨勢(shì)為降。GUAN和XIAM兩站可降水量變化為平緩上升，后再持續(xù)下降。SHAO在2002-06之前升，此后降。WUHN站可降水量變化為降-緩升-降。KMIN站點(diǎn)海拔明顯高于其他站，因而該站可降水量變化趨勢(shì)為先降后升。

圖3 GPS可降水量變化趨勢(shì)(亞熱帶季風(fēng)氣候類型)Fig.3 The trend of GPS precipitable water vapor(sub-tropical monsoon climate)

溫帶季風(fēng)氣候范圍為北緯35°～55°。從圖4可見(jiàn)，該氣候類型的4個(gè)站的可降水量變化趨勢(shì)基本一致，變化趨勢(shì)為先降后升。

URUM站點(diǎn)屬于溫帶大陸性氣候。從圖5可以看出，URUM站可降水量變化趨勢(shì)為先升后降。

圖4 GPS可降水量變化趨勢(shì)(溫帶季風(fēng)氣候類型)Fig.4 The trend of GPS precipitable water vapor(temperate monsoon climate)

圖5 GPS可降水量變化趨勢(shì)(溫帶大陸性氣候類型)Fig.5 The trend of GPS precipitable water vapor(temperate continental climate)

處于青藏高原高寒氣候區(qū)域的GPS站的可降水量明顯低于其他區(qū)域，原因可能是其溫度低、海拔高。由圖6可知，XNIN與LHAS兩站的可降水量變化趨勢(shì)為上升，而DLHA站的變化趨勢(shì)為緩升-下降。

4 結(jié)論

圖6 GPS可降水量變化趨勢(shì)(青藏高原高寒氣候類型)Fig.6 The trend of GPS precipitable water vapor(Alpine climate)

選取中國(guó)地殼運(yùn)動(dòng)監(jiān)測(cè)網(wǎng)絡(luò)的GPS數(shù)據(jù)，結(jié)合氣象數(shù)據(jù)，處理并獲得各站GPS可降水量序列，繼而用經(jīng)驗(yàn)?zāi)B(tài)分解得到各站點(diǎn)的可降水量變化趨勢(shì)，按照中國(guó)氣候類型對(duì)各站點(diǎn)進(jìn)行劃分并比較。除溫帶季風(fēng)氣候類型的GPS站可降水量變化為先降后升外，其他氣候類型區(qū)域的大部分GPS站的可降水量變化均呈現(xiàn)先升后降的趨勢(shì)特點(diǎn)。本文僅基于2000～2004年數(shù)據(jù)進(jìn)行研究，下一步可望針對(duì)更多的中國(guó)陸態(tài)環(huán)境監(jiān)測(cè)網(wǎng)絡(luò)GPS站、更長(zhǎng)期的觀測(cè)數(shù)據(jù)，進(jìn)行該方法的研究和驗(yàn)證。

1 王勇，柳林濤，許厚澤，等.利用GPS技術(shù)反演中國(guó)大陸水汽變化［J］.武漢大學(xué)學(xué)報(bào):信息科學(xué)版，2007，32(2):152－155.(Wang Yong，Liu Lintao，Xu Houze，et al.Retrieving change of precipitable water vapor in chinese mainland by GPS technique［J］.Geomatics and Information Science of Wuhan University，2007，32(2):152 －155)

2 Bevis M，Businger S，Herring T A，et al.GPS meteorology:remote sensing of atmospheric water vapor using the global positioning system［J］.J Geophys Res，1992，97:15 787 －15 801.

3 袁招洪，顧松山，丁金才.數(shù)值模式預(yù)報(bào)延遲量與GPS測(cè)量的比較研究［J］.南京氣象學(xué)院學(xué)報(bào)，2006，29(5):581－290.(Yuan Zhaohong，Gu Songshan，Ding Jincai.Comparison of MM5 model-simulated and GPS-observed zenith delay［J］.Journal of Nanjing Institute of Meteorology，2006，29(5):581 －290)

4 李國(guó)翠，李國(guó)平，陳小雷.強(qiáng)降雪天氣中GPS可降水量與地面空氣濕度的綜合分析［J］.高原氣象，2011，30(6):1 626－1 632.(Li Guocui，Li Guoping，Chen Xiaolei.Comprehensive analyses on GPS precipitable water vapor and surface humidity parameters during severe snowfall weather［J］.Plateau Meteorology，2011，30(6):1 626 －1 632)

5 李國(guó)平，陳嬌娜，郝麗萍.基于GPS-PWV的不同云系降水個(gè)例的綜合分析［J］.武漢大學(xué)學(xué)報(bào):信息科學(xué)版，2011，36(4):384 － 388.(Li Guoping，Chen Jiaona，Hao Liping.Case study of the rainfall processes in different cloud systems based on GPS-PWV data in chengdu plain［J］.Geomatics and Information Science of Wuhan University，2011，36(4):384 －388)

6 宋淑麗.地基GPS網(wǎng)對(duì)水汽三維分布的監(jiān)測(cè)及其在氣象學(xué)中的應(yīng)用［D］.上海:中科院上海天文臺(tái)，2004.(Song Shuli.Sensing three dimensional water vapor structure with ground-based GPS network and the application in meteorology［D］.Shnghai:Shanghai Astronomical Observatory，Chinese A-cademy of Sciences，2004)

7 Duan J，Bevis M，F(xiàn)ang P，et al.GPS Meteorology:direct estimation of the absolute value of precipitable water［J］.J Appl Meteorol，1996，35:830 －838.

8 李征航，徐曉華，羅佳，等.利用GPS觀測(cè)反演三峽地區(qū)對(duì)流層濕延遲的分布及變化［J］.武漢大學(xué)學(xué)報(bào):信息科學(xué)版，2003，28(4):393 －396.(Li Zhenghang，Xu Xiaohua，Luo Jia，et al.Inversion of the distribution and variation of ZWD over the Three Gorge area with GPS observation［J］.Geomatics and Information Science of Wuhan University 2003，28(4):393－396)

9 王勇，劉嚴(yán)萍.地基GPS氣象學(xué)原理與應(yīng)用研究［M］.北京:測(cè)繪出版社，2012.(Wang Yong，Liu Yanping.Theory and application of ground-based GPS meteorology［M］.Beijing:Surveying and Mapping Press，2012)

10 Norden E H，Shen Z，Long S R，et al.The empirical mode decomposition and the Hilbert spectrum for nonlinear and nonstationary time series analysis［D］.Royal Society A:Mathematical，Physical and Engineering Sciences，1998.

11 王勇，楊晶，張立輝，等.基于經(jīng)驗(yàn)?zāi)B(tài)分解與神經(jīng)網(wǎng)絡(luò)的信號(hào)預(yù)測(cè)［J］.大地測(cè)量與地球動(dòng)力學(xué)，2011(6):121－124.(Wang Yong，Yang Jing，Zhang Lihui，et al.Analysis of signal prediciton based on EMD and ANN［J］.Journal of Geodesy and Geodynamics，2011(6):121 －124)

12 朱金龍，邱曉輝.正交多項(xiàng)式擬合在EMD算法端點(diǎn)問(wèn)題中的應(yīng)用［J］.計(jì)算機(jī)工程與應(yīng)用，2006，23(2):72 －74.(Zhu Jinlong，Qiu Xiaohui.Application with orthogonal polynomial algorithm in the empirical mode decomposition endpoint problem［J］.Computer Engineering and Application，2006，23(2):72 －74)

13 鄒友峰，王勇，聞德保，等.利用不同氣候類型的GPS可降水量比較研究［J］.武漢大學(xué)學(xué)報(bào):信息科學(xué)報(bào)，2012，37(5):573 － 576.(Zou Youfeng，Wang Yong，Wen Debao，et al.GPS precipitable water vapor of different climate types［J］.Geomatics and Information Science of Wuhan University，2012，37(5):573 －576)

STUDY ON CHANGE OF PRECIPITABLE WATER VAPOR WITH GPS BASED ON EMPIRICAL MODE DECOMPOSITION

Liu Yanping1)，Wang Yong2)，Wen Debao3)
1)School of Civil Engineering，Central South University，Changsha 410075
2)School of Geology and Geomatics，Tianjin Chengjian Univesity，Tianjin 300384
3)School of Traffic and Transportation Engineering，Changsha University of Science＆ Technology，Changsha410004

A method of empirical mode decomposition was used to processe precipitable water vapor time series during 2000－2004 from GPS data of Crustal Movement Observation Network of China(CMONOC)and change trend of precipitable water vapor at every GPS station was obtained.The result of classification and comparison for the change trend indicates that the change trend in tropical monsoon climate zone，monsoon climate of subtropical zone，temperate continental climate zone and alpine climate zone is first up and then down;the change trend in temperate monsoon climate zone is first down and then up.

GPS;empirical mode decomposition;precipitable water vapor;CMONOC

P228.42

A

1671-5942(2014)05-0156-04

2013-05-17

湖南省研究生科研創(chuàng)新項(xiàng)目(CX2012B061);大地測(cè)量與地球動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室開(kāi)放基金項(xiàng)目(SKLGED2013-5-5-E);河北省教育廳重點(diǎn)項(xiàng)目(ZH2012060);國(guó)家自然科學(xué)基金項(xiàng)目(41174001)。

劉嚴(yán)萍，女，1979年生，講師，博士生，主要從事GPS應(yīng)用研究。E－mail:liuxiawy@126.com。

王勇，男，1978年生，博士，教授，主要從事GPS氣象學(xué)、InSAR大氣校正研究。E－mail:wangyongjz@126.com。