張小紅,郭 斐,郭博峰,呂翠仙
武漢大學測繪學院,武漢 430079
利用高頻GPS進行地表同震位移監(jiān)測及震相識別
張小紅,郭 斐,郭博峰,呂翠仙
武漢大學測繪學院,武漢 430079
結合日本3.11大地震,利用近震區(qū)和遠震區(qū)若干測站的高頻GPS觀測數據,采用動態(tài)精密單點定位技術分析了強震時復雜的地表運動時序,獲得了地表同震位移.在此基礎上提出利用S變換對高頻GPS地震波進行震相識別,有效拾取了P波初至,反演了P波傳播速度并與地震儀數據、USGS提供的參考值進行了對比和分析,為后續(xù)地震應急、地震預警提供了重要的、科學的決策依據.
高頻GPS,同震位移,日本地震,震相識別,P波初至,S變換
地震發(fā)生后,對其做出快速、準確的評估,對于地震應急、震后減災和救災、震后趨勢預測及強震引起的海嘯預警等至關重要.基于短時速報的地表同震位移為研究震源參數、潛在的斷層及地殼形變等事件提供了寶貴的信息.利用地震臺網布設的地震儀能夠在地震發(fā)生的數分鐘內較為準確地確定震源及震級,然而對于一些破壞力極大的強震,地震儀容易產生振幅飽和、速度或加速度積分經常出錯、噪聲放大和扭曲真實信號等現象[1].此外,基于地震儀估計得到的震級較為保守,容易低估地震的大小及其破壞程度.例如2004年12月26日的印度洋海嘯引發(fā)的蘇門答臘地震,利用地震波短時速報的震級(M8.0~8.5)遠小于其真實大?。∕9.2~9.3),從而嚴重影響了快速、準確的地震預警,造成了巨大的人員傷亡[2].近年來,隨著GPS觀測精度及處理方法的不斷改進與提高,尤其是高頻(1Hz)和超高頻(20~50Hz)GPS技術的出現,GPS逐漸成為地震儀的重要補充,為地震學研究提供了一種新的途徑.
利用高頻GPS進行地表位移監(jiān)測及地震分析最早于1994年由日本Kyoto大學災害預防研究所(DPRI)提出,其利用1Hz的GPS觀測數據,采用事后相對定位方式獲得了平面1~2cm的定位精度,并得出利用GPS能夠有效監(jiān)測近震、大幅的地表運動等結論[3].澳大利亞新南威爾士大學(UNSW)的Ge等將GPS天線及地震儀安置在震蕩的卡車頂部,利用10~20Hz的GPS觀測數據,采用RTK方式測試了高頻GPS用于監(jiān)測表面位移的可行性,結果表明由GPS監(jiān)測的位移與加速度計積分得到的位移兩者吻合較好[4-5].此后,Larson等人基于高頻GPS對2002年Denali M7.9級地震,2003年 Tokachi-oki M8.1級地震,2003年加州San Simeon M6.5 級地震,2004 年 Sumatra-Andaman M9.3級巨震,2008年汶川M8.0地震進行了研究,獲取了地震近場和遠場的地表運動時序[6-13].在地震/海嘯預警方面,Blewitt等利用IGS跟蹤站網15min左右、30s間隔的GPS觀測數據,研究了2004年Sumatra地震的震級、機制及其空間破裂尺度,得出其靈敏度(準確度)受衛(wèi)星軌道精度及采樣率影響等結論[2].然而,基于高頻和超高頻的GPS地震波震相識別及地震(海嘯)預警研究仍然較少.
因此,本文首先利用2011年3月11日日本仙臺港大地震當天若干測站1Hz的GPS觀測數據分析了各站的同震位移,并在此基礎上利用測站速度信息,采用S變換提取地震波P波初至時刻(震相識別),為強震預警提供支持和依據.
由于中國正處于工業(yè)化和城市化加速發(fā)展的階段,也正處于經濟增長和環(huán)境保護矛盾十分突出的時期,環(huán)境形勢依然十分嚴峻。一些地區(qū)環(huán)境污染和生態(tài)惡化還相當嚴重,主要污染物排放量超過環(huán)境承載能力,水、空氣、土壤等污染嚴重,固體廢物、汽車尾氣、持久性有機物等污染增加。發(fā)達國家上百年工業(yè)化過程中分階段出現的環(huán)境問題在中國集中出現, 環(huán)境與發(fā)展的矛盾日益突出。資源相對短缺、生態(tài)環(huán)境脆弱、環(huán)境容量不足,逐漸成為中國發(fā)展中的重大問題。
選取此次日本大地震近震及遠震區(qū)的6個IGS跟 蹤 站 (MIZU、USUD、DAEJ、SHAO、PETS、GUAM)1Hz的GPS觀測數據,同時利用距東京15km左右的JA01站1Hz的GPS觀測數據,采用武漢大學研制的TriP軟件進行動態(tài)精密單點定位解算[14],衛(wèi)星軌道及鐘差采用COD中心提供的15min間隔的精密星歷與5s間隔的精密衛(wèi)星鐘差,衛(wèi)星截止高度角設置為5°.各站分布及震中位置(星形)如圖1所示.MIZU站距震中最近,大約為140km;GUAM站距震中最遠,達到2700余千米.
盡管遠震區(qū)(>600km)的高頻GPS能夠監(jiān)測到地震波信號,但其破壞力較小,同震位移容易被觀測噪聲所淹沒,因此本文僅給出近震區(qū)測站的震時水平向運動軌跡(圖3所示)及其震后位移變化如表1所示.在時域內,圖2中已經詳細反映了測站在N、E、U方向的運動過程,因此,圖3側重于從空間域(水平向)描述測站的點位運動過程.
圖1 震中及高頻GPS測站分布示意圖Fig.1 Location of the HR-GPS stations and the epicenter of the M9.0 Tohoku earthquake
基于高頻GPS可以直接獲得地表的真實位移,捕獲強震時復雜的地表運動時序,為研究斷裂擴展機理提供可視化幫助.圖2自上而下依次為所有測站北、東、垂直方向的地表同震位移時間序列,橫軸代表UTC時間,縱軸代表相應的位移量,為便于描述地震波的傳播過程,將所有測站依震中距大小排列顯示其震前、震時及震后的位移變化.從圖中可以看出地震發(fā)生后,距離震中最近的MIZU站最先感知,隨后地震波傳至其他各站.各站各分量在地震到時均有不同程度的震動,其震動大小主要取決于震中距及地震帶的破裂方向.尤其以震中距較小的MIZU 站 (~140km)、JA01 站 (~390km)及USUD站(~450km)最為明顯,對于數千千米之外的DAEJ站、SHAO站、PETS站和GUAM站利用高頻GPS同樣能夠監(jiān)測到其地震波信號.此外,由于破壞力極強的剪切波(S波)主要造成水平向的擠壓和扭動,因此,水平向位移較豎向位移變化更加顯著.
需要注意的是,本法第四十五條規(guī)定,評估專業(yè)人員簽署虛假評估報告,構成犯罪的,終身不得從事評估業(yè)務。根據《刑法》第二百二十九條的規(guī)定,評估專業(yè)人員簽署虛假評估報告,情節(jié)嚴重的,將構成提供虛假證明文件罪,處五年以下有期徒刑或者拘役,并處罰金。本法第四十五條的規(guī)定是針對此類犯罪的特別規(guī)定,也是對第十一條的例外規(guī)定。因為評估專業(yè)人員簽署虛假評估報告,構成虛假證明文件罪的,是最嚴重的評估違法行為,應當給予最嚴厲的懲處,除應當依法追究刑事責任外,還應當接受終身禁業(yè)的懲罰,即終身不得從事評估業(yè)務。
圖2 GPS同震位移時間序列(北/東/垂向)Fig.2 Time series of coseismic displacement using HR-GPS(North/East/Up)
學校對于在日常生活、學習、工作等場景的靈活遷移應用的示范性輸入材料,則更多為中文紙制材料,主要為英語書蟲等文學經典作品改編的中高階英語閱讀材料,生活常識、科普類的閱讀和音像材料幾乎沒有。適用于入門級的英語繪本等閱讀材料,也非常有限。
圖3 近震測站震時水平向運動軌跡(UTC 05∶46∶00—05∶56∶00)Fig.3 Horizontal trajectory of near-field stations(during UTC 05∶46∶00—05∶56∶00)
圖4 MIZU站震前速度波形S變換譜Fig.4 S-transform spectrum of station MIZU waveform before earthquake
表1 近震區(qū)高頻GPS測站同震位移Table 1 Coseismic permanent displacements of near-field stations
連續(xù)時間序列函數h(t)的S變換公式如下:
地震預警的基本思想是利用地震P波傳播速度快于S波與面波的特點,通過分析P波初到時刻、振幅及頻率等信息,估算隨后而至的破壞性震動的大小及范圍,進而確定預警范圍與級別.盡管預警時間非常短暫(數秒至數十秒),但卻具有極其重要的社會效益.因此基于高頻GPS的強震預警及震源參數反演首先需要解決的是初至波的拾取,即P波的震相識別.常用的震相識別方法有:人工經驗法、STA/LTA法、傅里葉變換法、小波分析法、AIC方法、瞬時頻率法等.
S變換(S-Transform)是由美國地球物理學家Stock well等在1996年提出的一種無損可逆的線性時頻分析方法,它綜合了短時傅里葉變換和小波變換的優(yōu)點,具有較高的時頻分辨率,在信號處理領域得到了廣泛應用[15].由于地震波初至前一般表現為隨機噪聲(背景噪聲),初至后則是地震信號與背景噪聲的混合.因此,地震波初動時將形成振幅能量的突變.當信噪比較高時,在時域上很容易識別;當信噪比較低時,初至前后能量差異小,在時域難以識別,但在時頻域內,由于地震波初至信號與背景噪聲占用不同的時頻區(qū)域,使得初至識別仍然可行.
圖3中坐標原點為各站震前初始位置(Start),地震波(尤其是橫向剪切波)到達時測站開始朝東南方向發(fā)生劇烈的抖動,其中MIZU站的最大位移達到東向3m、南向2m左右;隨著地震波的衰減,測站位置發(fā)生回彈至最終點位(End).相對于 MIZU站,JA01、USUD站的位移量級較小,將其運動軌跡放大(Enlarge)進行圖示,得到圖3中左半部分.表1說明強震已經造成這些測站產生了永久性位移.其中MIZU站向南偏移1.133m,向東偏移2.130m,垂直方向也向下沉降0.1m左右.其他測站也呈向東漂移和向下沉降的趨勢,這在一定程度上也反映強震造成日本所在的板塊整體朝東運動,并略有沉降.
為驗證本文所提方法的正確性及有效性,結合USGS網站提供的P波初到時刻作為參考值,對比分析各站的拾取精度,結果如表2所示.比較表2中各站S變換拾取的P波初到時刻與USGS提供的參考值發(fā)現,幾乎所有測站拾取的結果均比USGS提供的結果晚5~8s左右,對于高頻GPS測站(MIZU、JA01),其水平向拾取精度明顯優(yōu)于高程方向.比較 MIZU(高頻GPS)與IWT011(地震儀)這兩個相鄰的測站,發(fā)現其地震波平均到時也相差2~3s.究其原因,這是因為:其一,受GPS本身觀測精度的限制,基于GPS技術獲取的地震波信噪比不及地震儀數據,導致使用S變換不容易區(qū)分背景噪聲與微弱的地震信號(尤其是垂直方向);而作為區(qū)分背景噪聲與地震信號的能量閾值選取極其關鍵,較小的閾值差異可能引起較大的拾取誤差(達到數秒);且不同方向對地震波初至的敏感程度不同,對其簡單的求取平均到時也會引起拾取結果產生一定的偏差.其二,盡管已將MIZU與IWT011的時間歸化至UTC時,但時鐘本身的誤差導致兩者并非嚴格同步,因此基于GPS和地震儀的拾取結果也會有所差異.此外,由于USGS提供的參考值是基于IASP91地球模型計算得到的理論值,它并未考慮橢球改正及嚴密的介質傳播速度等信息,可能會與實際情況存在較大的差異[16].因此,本文采用S變換拾取的地震波初至與USGS提供的參考結果也并非完全一致.但總體而言,基于GPS與地震儀數據在震相識別時具有較好的吻合性.
本文首先利用MIZU和JA01兩個GPS測站的速度信息,采用S變換提取P波初至時刻,并結合鄰近地震儀的數據(K-NET網,IWT011站,距MIZU站約2km)、美國地質調查局(USGS)提供的參考數據進行對比分析,最后利用拾取的地震波初至時刻與震中距反演出P波的傳播速度.
圖7中從左向右、自上而下依次為N、E、U三個方向對應的時頻譜圖和加速度時序圖.雙向豎線表示P波初到時刻(具體時刻如表2所示,精確到0.01s),且P波的歸一化頻率在各方向上具有較好的一致性(0.02~0.03Hz).
圖4顯示,震前速度波形表現為高頻噪聲,其頻率主要分布在0.2~0.5Hz之間,且豎向背景噪聲的能量較水平向明顯偏大,這也與GPS測量垂直方向較水平方向的觀測噪聲較大相符.北、東、垂直方向基底噪聲的平均能量分別為0.001、0.0015和0.002左右,這些數值將作為后續(xù)地震波P波初至拾取的能量閾值.
為了驗證GPS地震波P波拾取的準確性,采用日本K-NET網絡中與MIZU站相近的IWT011站100Hz的地震儀數據對初到時間進行P波拾取.類似地,采用S變換獲得了該站的震前、震時頻譜圖,如圖7所示.
由圖5、6所知,地震波初至時,N、E、U三個方向在時頻域內均發(fā)生了較大的能量突變,且隨著后續(xù)S波、面波的相繼到達,能量逐漸增強至最大,并隨著地震波的衰弱逐漸消退.基于S變換的P波拾取結果如表2所示,受GPS采樣率(1Hz)的限制,P波拾取精度僅能精確到秒級.由N、E、U三個方向拾取的P波初至時刻不盡相同,差異達到數秒,且水平向較垂直方向的拾取時刻偏早.
圖5 MIZU站震時頻譜圖Fig.5 S-transform spectrum of station MIZU
圖6 JA01站震時頻譜圖Fig.6 S-transform spectrum of station JA01
表2 各測站不同拾取類別對應的P波初到時刻Table 2 Arrival time of P-wave of station MIZU,JA01 and IWT011using S-transform and the IASP91earth model(USGS method)
基于S變換去噪后的MIZU、JA01測站速度時頻譜如圖5、6所示,圖中從左向右、自上而下依次代表北、東、垂直三個方向對應的頻譜圖和速度時序圖,橫軸為UTC時,縱軸為歸一化頻率f或速度V值.
精確飼喂是本設備的核心功能。飼喂站包括精確下料系統和食槽。精確下料系統可以依據個體豬的體重和孕期確定對其供食量。該系統主要由料倉、螺旋下料裝置、電機及傳動系統和落料口等部分組成。
經過80余年的學科建設,該校體育教育專業(yè)已形成明確的人才培養(yǎng)理念:學生能在健康和體育教育以及第二個專業(yè)化教學領域擔任中小學教師,還可以從事社區(qū)衛(wèi)生、娛樂等方面的工作,就業(yè)范圍進一步擴大到教練員、科研人員、體育管理人員、體育新聞報道人員、社區(qū)娛樂及體育健身指導人員等。具體培養(yǎng)目標包括:培養(yǎng)合格的健康和體育教育工作者,掌握應對K-12年級以下教學任務的能力;樹立社會責任感,成為健康促進的領導者和他人增進健康的榜樣;通過參與體育運動,掌握團隊合作等日常人際交往技能;涉獵多學科內容,能夠在實際教學中運用應用科學及社會科學等知識。
圖7 IWT011站震時頻譜圖(A為加速度)Fig.7 S-transform Spectrum of Seismograph IWT011
以MIZU站為例,選取其震前900個歷元(15min)的位移序列進行一階差分消除長周期低頻噪聲,得到測站速度波形,采用S變換將其速度波形轉換至時頻域進行頻譜分析,估計背景噪聲的頻率范圍及其能量大小,得到該站速度背景噪聲的S變換譜,如圖4所示,其中橫軸為歷元數,縱軸為歸一化頻率.
其中,t表示時間,f表示頻率,τ為高斯窗函數的平滑因子.對于1維的函數h(t),變量τ和其中一個頻率fi可確定一個S(τ,fi),我們稱S(τ,fi)所確定的值為一道噪聲.因此,當頻率取不同值時,S變換可以獲得若干道噪聲向量.
綜上所述,三伏灸辨證分型聯合西藥治療支氣管哮喘患者,不僅提高了患者治療效果,還改善了FEV1、FVC,值得推廣。
基于上述拾取的P波,根據各站的震中距反演得到MIZU、JA01兩站的P波傳播速度分別為5.06km/s和5.82km/s.根據P波在不同的介質對應的傳播速度一般為5~8km/s,因此,這也在一定程度上驗證了本文基于S變換拾取的震相為P波.
小學生這個年齡段就是對世界充滿了很強地好奇心還有求知欲,并且,孩子們帶著活潑好動的天性、注意力難以集中等一些性格特征。所以,在進行小學數學授課的過程當中,假若老師選擇先進規(guī)范的教學措施,便可以調動起孩子們學習數學的積極性,進一步提升教學質量;相反,課堂授課死板乏味,孩子們便會難以集中注意力聽講。所以,在進行生活化授課過程當中,老師應該主動創(chuàng)設活潑、和諧的數學課堂氣氛,指引孩子們保持良好的心態(tài)走進課堂學習里面去。
本文以日本仙臺大地震為例,首先利用近震區(qū)和遠震區(qū)若干測站的高頻GPS觀測數據采用動態(tài)精密單點定位技術獲取了高頻GPS地震波,分析了強震時復雜的地表運動時序,獲得了地表的同震位移,結果表明:地震到時各站均有不同程度的震動,其震動大小主要取決于震中距及地震帶的破裂方向,即使對于數千千米之外的DAEJ站、SHAO站、PETS站和GUAM站利用高頻GPS同樣能夠監(jiān)測到其地震波信號;由于破壞力極強的剪切波(S波)主要造成水平向的擠壓和扭動,導致水平向位移較垂直向位移變化更加劇烈;強震造成近震區(qū)內的部分測站已產生永久性變形,日本所在的板塊整體朝東運動,并略有下沉.
提出利用S變換進行高頻GPS地震波的震相識別,基于MIZU和JA01兩個GPS測站的速度信息,采用S變換提取了P波初至時刻,并與USGS提供的參考值進行對比分析,同時選擇距 MIZU 2km左右的IWT011站的地震儀數據進行S變換,將其P波拾取結果與USGS提供的參考值進行對比分析.結果表明:受GPS測量精度及采樣率的限制,同一測站N、E、U三個方向拾取的P波初至時刻差異達到數秒,且水平向較垂直方向的拾取時刻早;不同測站基于S變換拾取得到的結果與USGS提供的參考值差異達到5~8s,相鄰GPS測站與地震儀臺站的地震波初至時刻吻合較好,但仍有2~3s的差異.由地震波初到時刻及震中距反演得到的P波傳播速度為5~6km/s,與實際較為符合.
今后,隨著高頻和超高頻GPS技術的發(fā)展,尤其是高程方向定位精度的提高對于改善高頻GPS地震波信噪比、提高震相識別精度,進而為地震應急、預警提供更加科學的決策依據具有十分重要的意義,也是今后急需解決的問題.
致 謝 衷心感謝IGS、USGS、K-NET網絡以及DLR的Thomas Dautermann為本文提供了豐富的數據及產品服務;在此一并表示感謝!
(References)
[1] Boore D M,Stephens C D,Joyner W B.Comments on baseline correction of digital strong-motion data:examples from the 1999Hector Mine,California,earthquake.Bull.Seismol.Soc.Am.,2002,92(4):1543-1560,doi:10.17850120000926.
[2] Blewitt G,Kreemer C,Hammond W C,et al.Rapid determination of earthquake magnitude using GPS for tsunami warning systems.Geophys.Res.Lett.,2006,33:L11309,doi:10.1029/2006GL026145.
[3] Hirahara K,Nakano T,Hoso Y,et al.An experiment for GPS strain seismometer.∥ Proc.Japanese Symposium on GPS.Tokyo,Japan,1994:67-75.
[4] Ge L.GPS seismometer and its signal extraction.∥ Proc.12th Int.Tech.Meeting of the Satellite Division of the U.S.Inst.of Navigation(ION GPS'99).Nashville,Tennessee,1999.
[5] Ge L L,Han S W,Rizos C,et al.GPS seismometers with up to 20Hz sampling rate.Earth.Planets and Space,2000,52(10):881-884.
[6] Larson K M,Bodin P,Gomberg J.Using 1-Hz GPS data to measure deformations caused by the Denali fault earthquake.Science,2003,300(5624):1421-1424.
[7] Bock Y,Prawirodirdjo L,Melbourne T I.Detection of arbitrarily large dynamic ground motions with a dense highrate GPS network.Geophys.Res.Lett.,2004,31:L06604.doi:10.1029/2003GL019150.
[8] Irwan M,Kimata F,Hirahara K,et al.Measuring ground deformations with 1-Hz GPS data:the 2003Tokachi-oki earthquake(preliminary report).Earth,Planets and Space,2004,56(3):389-393.
[9] Miyazaki S,Larson K M,Choi K,et al.Modeling the rupture process of the 2003September 25Tokachi-Oki(Hokkaido)earthquake using 1-Hz GPS data.Geophys.Res.Lett.,2004,31:L21603,doi:10.1029/2004GL021457.
[10] Ji C,Larson K M,Tan Y,et al.Slip history of the 2003San Simeon earthquake constrained by combining 1-Hz GPS,strong motion,and teleseismic data.Geophys.Res.Lett.,2004,31:L17608,doi:10.1029/2004GL020448.
[11] Ohta Y,Meilano I,Sagiya T,et al.Large surface wave of the 2004Sumatra-Andaman earthquake captured by the very long baseline kinematic analysis of 1-Hz GPS data.Earth,Planets and Space,2006,58(2):153-157.
[12] Shi C,Lou Y,Zhang H,et al.Seismic deformation of the Mw8.0Wenchuan earthquake from high-rate GPS observations.Adv.Space Res.,2010,46(2):228-235,doi:10.1016/j.asr.2010.03.006.
[13] 殷海濤,甘衛(wèi)軍,肖根如等.利用高頻GPS技術進行強震地面運動監(jiān)測的研究進展.地球物理學進展,2009,24(6):2012-2019.Yin H T,Gan W J,Xiao G R,et al.Progress on monitoring strong earthquake ground motions using high-rate GPS.Progress in Geophysics (in Chinese),2009,24(6):2012-2019.
[14] 張小紅,劉經南,Forsberg R.基于精密單點定位技術的航空測量應用實踐.武漢大學學報 (信息科學版),2006,31(1):19-22,46.Zhang X H,Liu J N,Forsberg R.Application of precise point positioning in airborne survey.Geomatics and Information Science of Wuhan University (in Chinese),2006,31(1):19-22,46.
[15] Stockwell R G,Mansinha L,Lowe R P.Localization of the complex spectrum:The S transform.IEEE Transactions on Signal Processing,1996,44(4):998-1001.
[16] Kennett B L N,Engdahl E R.Travel times for global earthquake location and phase identification.Geophys.J.Int.,1991,105(2):429-465.
Coseismic displacement monitoring and wave picking with high-frequency GPS
ZHANG Xiao-Hong,GUO Fei,GUO Bo-Feng,LüCui-Xian
School of Geodesy and Geomatics,Wuhan University,Wuhan 430079,China
High-rate kinematic analysis coupled with global positioning system (GPS)has become one of the effective methods to observe the dynamic crustal deformations caused by earthquake.In this article,kinematic precise point positioning(PPP)technology was used for monitoring the complex surface vibration behavior and coseismic displacements with high-frequency GPS data from stations in near-field and teleseismic areas during the 11March 2011moment magnitude 9.0 Tohoku Earthquake.Furthermore,we proposed the S-transform algorithm for the primary wave(P-wave)arrival picking from high-rate GPS (HR-GPS)seismic waves.And its performances were compared with nearby seismograph records and the theoretical value provided by U.S.Geological Survey(USGS).
High-rate GPS data,Coseismic displacement,Tohoku earthquake,Wave picking,Arrival time,S-transform
10.6038/j.issn.0001-5733.2012.06.012
P228
2011-09-09,2012-03-05收修定稿
國家自然科學基金(41074024)、創(chuàng)新研究群體科學基金(41021061)聯合資助.
張小紅,男,1975年生,教授,主要從事GNSS精密定位及GNSS地學應用等方面的教學與研究工作.E-mail:xhzhang@sgg.whu.edu.cn
張小紅,郭斐,郭博峰等.利用高頻GPS進行地表同震位移監(jiān)測及震相識別.地球物理學報,2012,55(6):1912-1918,
10.6038/j.issn.0001-5733.2012.06.012.
Zhang X H,Guo F,Guo B F,et al.Coseismic displacement monitoring and wave picking with high-frequency GPS.Chinese J.Geophys.(in Chinese),2012,55(6):1912-1918,doi:10.6038/j.issn.0001-5733.2012.06.012.
(本文編輯 胡素芳)