王坦 李瑜 張銳 師宏波 王閱兵

摘要：回顧了GPS技術(shù)應用于我國地震監(jiān)測中的發(fā)展歷程，重點從GPS數(shù)據(jù)處理與時間序列、中國大陸構(gòu)造變形速度場與應變率場、中國大陸活動地塊運動定量化和地震研究等方面，介紹了GPS在我國地震監(jiān)測中的應用現(xiàn)狀，并針對當前面臨的問題與挑戰(zhàn)，進行了討論，從地震預測研究、基礎研究、GPS臺網(wǎng)監(jiān)測布局和GPS解算精度與時效4個方面，對GPS在我國地震監(jiān)測中的發(fā)展進行了展望。

關(guān)鍵詞：GPS;地震監(jiān)測;陸態(tài)網(wǎng)絡;應用現(xiàn)狀;發(fā)展展望

中圖分類號：P315.725?? 文獻標識碼：A?? 文章編號：1000-0666（2021）02-0192-16

0 引言

地殼運動與地震有著密切的關(guān)系，因此，對地殼運動的監(jiān)測一直是地震預測預報研究的重要依據(jù)和手段。但受限于觀測技術(shù)和精度，1966年邢臺地震后至1998年，我國主要是采用傳統(tǒng)大地測量方法開展大地形變測量（王琪等，2020），包括水準測量、激光測距和三角測量等。傳統(tǒng)大地測量效率低，觀測周期長，長期以來只能在某些點和小范圍進行監(jiān)測，難以對大范圍、特別是斷層帶的形變場實施動態(tài)監(jiān)測。雖然也曾有研究人員通過大規(guī)模人工觀測獲取中國大陸垂直形變圖像（張祖勝等，1989;薄萬舉等，2007），但在時空分辨率和精度上并不能滿足構(gòu)造變形研究和地震預測的需要，嚴重制約了地震監(jiān)測預報的發(fā)展。

20世紀90年代初出現(xiàn)的GPS對地觀測技術(shù)為大地測量和地殼運動監(jiān)測帶來了一場革命。GPS高精度、大范圍和準實時的地殼運動觀測結(jié)果，使得多種空間尺度、多種時間分辨率和寬頻率域的現(xiàn)今地殼運動觀測成為可能（賴錫安等，2004）。隨著人們對導航定位性能要求的不斷提升，以及航天、通信及衛(wèi)星導航技術(shù)本身的發(fā)展，原始GPS在性能上的局限性已越來越明顯，新一代GNSS，如現(xiàn)代化的GPS、不斷改造的GLONASS以及新建的歐盟Galileo、我國的北斗衛(wèi)星導航系統(tǒng)（BDS），都有了新的改進，處于迅速發(fā)展之中?，F(xiàn)今的導航定位技術(shù)無論在空間與時間的覆蓋性上還是在精度上，都取得了革命性的進步，應用領域得到了極大拓寬，這些給我們帶來前所未有的機遇的同時，也向我們提出了新的挑戰(zhàn)。本文主要對GPS技術(shù)應用于我國地震監(jiān)測的30年歷史和取得的成就進行了回顧和評述，并在此基礎上，結(jié)合國外的最新研究進展，對當前面臨的困境進行總結(jié)，并提出了未來的發(fā)展方向。

1 GPS觀測臺網(wǎng)

1.1 早期觀測

我國最早的GPS地殼運動觀測網(wǎng)是中國地震局1988年建設的滇西地震預報實驗場GPS觀測網(wǎng)（王琪，Seeber，1998）。在20世紀90年代初，我國地震部門陸續(xù)在川滇、河西走廊（黃立人，馬青，2003）、青藏高原（游新兆，王啟梁，1994;蔡宏翔等，1997;陳俊勇等，2001;張培震等，2002a）、新疆天山（王琪等，2003）、華北（李延興，趙承坤，1998）、福建東南沿海（劉序儼等，1999）等地區(qū)零星地開展了小范圍、定期觀測的GPS地殼運動監(jiān)測研究工作（1988—1997年）。同一時期，中國科學院上海天文臺建設的國家攀登計劃項目“現(xiàn)今地殼運動與動力學”建成了覆蓋中國大陸22個點位的GPS觀測網(wǎng)，該網(wǎng)分別于1992、1994和1996年進行了3期聯(lián)測，大部分測點的速度精度都優(yōu)于3 mm/a（朱文耀等，1998;李延興等，2000;孟國杰，2001）。國家測繪局完成了國家A、B級網(wǎng)共800多個點位的布設和觀測，東部較密，西部稀疏，主要服務于國民經(jīng)濟建設。總參測繪局完成了國家GPS一、二級網(wǎng)500多個點位的布設和觀測，均勻地分布于全國各地。

以上這些早期觀測通常都采用流動觀測模式進行（王小亞等，2002），每次觀測時長在1到數(shù)天，主要關(guān)注于一些局部區(qū)域的形變特征，難以提供全國尺度連續(xù)的地殼運動特征。

1.2 中國地殼運動觀測網(wǎng)絡

1998年，“九五”期間，由中國地震局牽頭啟動實施了國家重大科學工程“中國地殼運動觀測網(wǎng)絡工程項目”（簡稱“網(wǎng)絡工程”）。整個工程由基準網(wǎng)、基本網(wǎng)、區(qū)域網(wǎng)和數(shù)據(jù)傳輸與分析處理系統(tǒng)4大部分組成?；鶞示W(wǎng)是網(wǎng)絡工程的基本框架，由25個GPS連續(xù)觀測站構(gòu)成;基本網(wǎng)由56個定期復測的GPS站組成;區(qū)域網(wǎng)由1 000個不定期復測的GPS站組成;數(shù)據(jù)系統(tǒng)由1個數(shù)據(jù)中心和3個數(shù)據(jù)共享子系統(tǒng)構(gòu)成。網(wǎng)絡工程在中國大陸建立了由25個連續(xù)GPS基準站，56個基本站和1 000個流動GPS區(qū)域站組成的GPS監(jiān)測網(wǎng)絡。（張祖勝，2001;牛之俊等，2002;陳俊勇，1997）。基準站相鄰站間距離平均約為700 km，主要功能是監(jiān)測中國大陸一級塊體?；鶞收鞠噜徴鹃gGPS基線長度年變化率測定的精度優(yōu)于2 mm/a（王琪等，2003;游新兆等，2001）。基本站作為基準站的補充，主要用于一級塊體本身及塊體間的地殼變動的監(jiān)測。它與基準站一起均勻布設，平均站距約為350 km，相鄰站間GPS基線每期測定精度：水平分量優(yōu)于5 mm，垂直分量優(yōu)于15 mm。區(qū)域站中約700個站集中分布在主要構(gòu)造帶和地震帶上，用于監(jiān)測其活動狀況，為地震預測預報服務;約300個站均勻分布在全國，作為基準站和基本站的補充，相鄰站間GPS基線每期測定精度：水平分量優(yōu)于5 mm，垂直分量優(yōu)于15 mm。網(wǎng)絡工程以監(jiān)測地殼運動服務于地震預報為主要目標，并兼顧大地測量和國防建設的需要，以GPS技術(shù)為主，輔之以VLBI、SLR、重力和水準觀測，其規(guī)模、觀測密度和精度等方面都將中國大陸GPS地殼運動觀測水平推向了一個新高度。

1.3 中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡

作為網(wǎng)絡工程的延伸，“十五”期間建成的中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡（簡稱“陸態(tài)網(wǎng)絡”）（甘衛(wèi)軍等，2007，2012;游新兆，馬海建，2012），使得我國大陸連續(xù)運行的GPS基準站由網(wǎng)絡工程的25個擴展到260個，區(qū)域流動GPS觀測站點由1 000個擴充至2 000個。一期工程的基本站在二期工程中將不再單獨存在，其原有的55個定期觀測的GPS基本站，除部分被升級改造為連續(xù)觀測基準站外，其余站點歸入2 000個區(qū)域站之中。

基準站布設方面，在中國大陸的20個主要活動塊體上各布設3個基準站，共計60站;在我國總長約32 000 km的活動斷層上平均每200 km（大致為一個8級地震的斷裂破裂長度）布設1個基準站，共計160站;為實現(xiàn)對活動強烈的主要斷裂的有效監(jiān)測，在地震活動最為強烈的中國大陸西部地區(qū)布設垂直于斷裂的密集線狀觀測臺陣和以監(jiān)測慢地震為目標的斷裂兩側(cè)的成對臺陣，共計32站;在中國大陸周邊海域的島礁上布設相應臺站，共計4站;在監(jiān)測塊體涉及的中國大陸周邊國家布設相應臺站，共計4站?；鶞收菊鹃g距縮小至160～350 km，水平分量測定精度優(yōu)于2 mm/a。

區(qū)域站在網(wǎng)絡工程的基礎上新增1 000個，采用按需要不定期復測的方法以補充基準站的不足。新建的區(qū)域站中，有400個區(qū)域站組成50個跨越主要活動斷裂帶的密集線狀觀測點陣，用于研究斷層結(jié)構(gòu)和形變特征;300個用于25個地震危險區(qū)的加密，點間距30～70 km;另300個在全國均勻分布。全國的大部分地區(qū)區(qū)域站的平均間距縮小至100 km左右，有可能監(jiān)測到5級地震產(chǎn)生的變形。區(qū)域站坐標變化量測定精度：水平分量優(yōu)于5 mm，垂直分量優(yōu)于10 mm。

陸態(tài)網(wǎng)絡與美國PBO和日本GEONET地球觀測網(wǎng)絡一同成為世界上性能指標最先進的三大地殼運動觀測網(wǎng)絡，成為國際地球科學研究與發(fā)展的基礎平臺（甘衛(wèi)軍等，2012）。陸態(tài)網(wǎng)絡在建站技術(shù)、觀測技術(shù)、集成技術(shù)及其綜合應用方面，實現(xiàn)了與世界先進水平同步發(fā)展。在布網(wǎng)設計上，采取基準站與區(qū)域站有機結(jié)合的方式，形成自成特色的站網(wǎng)體系，臺站的空間分布達到與國外同類網(wǎng)絡相同的水平，時間分辨率在重點地區(qū)接近世界先進水平。陸態(tài)網(wǎng)絡所采取的整體規(guī)劃、統(tǒng)一布局、重點布設、統(tǒng)一標準的建設模式，在充分發(fā)揮整體性優(yōu)勢的同時，更具針對性和實用性，并保證了高質(zhì)量數(shù)據(jù)的產(chǎn)出，提高了國家投資效益。陸態(tài)網(wǎng)絡的建成實現(xiàn)了對中國大陸主要活動塊體和活動斷層的覆蓋，為獲取中國大陸地殼運動細部特征提供基礎，使我國綜合利用空間對地觀測技術(shù)，監(jiān)測大陸巖石圈、大氣圈動態(tài)變化，特別是與地震有關(guān)的變形監(jiān)測的水平達到國際先進（甘衛(wèi)軍等，2012）。與此同時，中國地震局自2009年開始，針對地震重點危險區(qū)和重點監(jiān)視區(qū)，先后實施了華北強震強化跟蹤監(jiān)測項目、綜合地球物理場觀測計劃（青藏高原東緣地區(qū)和鄂爾多斯地塊周緣地區(qū)）和大華北地區(qū)項目等，不斷完善GPS觀測。

1.4 GPS數(shù)據(jù)資源共享與發(fā)展

中國地震局自2014年開始以陸態(tài)網(wǎng)絡為紐帶，逐步開展了同氣象、測繪和國土等部門的全國GPS數(shù)據(jù)資源共享工作。截至2020年，已形成含陸態(tài)網(wǎng)絡260個共計超過2 000個連續(xù)GPS觀測站的全國共享GPS觀測網(wǎng)，進一步擴大了監(jiān)測范圍，提升了監(jiān)測能力，提高了資源利用率。東部GPS流動觀測逐步以共享站代替，重點加強西部地區(qū)GPS流動觀測，實現(xiàn)了GPS觀測網(wǎng)絡的大發(fā)展。至此，中國大陸GPS觀測實現(xiàn)了從全國大范圍空間覆蓋以流動觀測為主向固定連續(xù)站發(fā)展的進步。連續(xù)觀測不僅實現(xiàn)了高時間分辨信息的獲取，同時由于GPS包含的非構(gòu)造不規(guī)則的周期信息，連續(xù)觀測還實現(xiàn)了震間長期高精度地殼運動速率的獲取。

這些日益豐富的GPS觀測與研究，開拓了中國大陸現(xiàn)今地殼運動-變形和動力學研究新領域，直接支持了大陸動力學、地震科學、地震大地測量學（周碩愚等，2017）和地震預測的創(chuàng)新，同時也為中國大陸的地殼運動和地球動力學研究提供了十分重要的基礎資料和約束條件。

2 GPS在我國地震監(jiān)測中的應用現(xiàn)狀

2.1 GPS數(shù)據(jù)處理與時間序列

2.1.1 GPS數(shù)據(jù)處理軟件

目前，用于地殼運動研究的GPS數(shù)據(jù)處理主要采用的是精密單點定位軟件（如GIPSY軟件、Panda軟件）、雙差定位軟件（如GAMIT軟件）和同時具備這兩種功能的軟件（如Bernese軟件）（顧國華等，2001，2005，2007）。

精密單點定位（PPP，Precise Point Positioning）技術(shù)由美國噴氣推進實驗室（JPL）的Zumberge等（1997）提出。它是通過外部高精度衛(wèi)星軌道數(shù)據(jù)和鐘差產(chǎn)品，對單個測站的非差導航衛(wèi)星觀測偽距值及載波相位觀測值進行處理，獲得分米級、厘米級甚至毫米級定位結(jié)果的一種導航數(shù)據(jù)處理技術(shù)。Zumberge等（1997）還基于平方根濾波方法研制出精密單點定位解算軟件GIPSY并進行相關(guān)實驗，單天解的精度為：水平方向1 cm，高程方向2 cm。經(jīng)過20多年的快速發(fā)展，國內(nèi)外學者研發(fā)出了許多成熟的PPP軟件。瑞士伯爾尼大學天文研究所在其研制的BERNESE5.0中添加了PPP.PCF模塊，可進行GPS/GLONASS PPP 解算（Dach et al，2007）。武漢大學劉經(jīng)南等基于信息平方根濾波法（SRIF）研制出PPP解算軟件 PANDA（Liu，Ge，2003;趙齊樂等，2005）。

雙差定位即在兩測站之間以及兩顆衛(wèi)星之間分別做差，是較早出現(xiàn)的、成熟的高精度定位方法，廣泛應用于高精度GPS測量數(shù)據(jù)處理軟件。因為要做差，就需要選擇基準站和流動站，選擇參考衛(wèi)星和非參考衛(wèi)星。雙差定位有松組合模型（LCM，Loosely Combined Model）和緊組合模型（TCM，Tightly Combined Model）兩種?，F(xiàn)在主流的方法是松組合模型，在每個衛(wèi)星導航系統(tǒng)內(nèi)分別獨立選擇各自的參考衛(wèi)星，只在系統(tǒng)內(nèi)組建雙差觀測方程，再聯(lián)合不同系統(tǒng)的雙差方程進行組合定位解算。緊組合模型，即在多系統(tǒng)組合定位時，不同的GNSS系統(tǒng)在組雙差觀測方程時，只選取某一顆衛(wèi)星作為參考衛(wèi)星，其它衛(wèi)星均與其做差，這樣可以形成系統(tǒng)內(nèi)和系統(tǒng)間的雙差觀測值（張小紅等，2016）。在松組合模型中，各個系統(tǒng)的觀測值相互獨立，實現(xiàn)簡單，己被廣泛采用。但是，松組合模型沒有考慮不同GNSS之間的相關(guān)性，沒有充分利用多系統(tǒng)觀測數(shù)據(jù)，不能真正意義上實現(xiàn)GNSS高精度相對定位中GNSS系統(tǒng)間的兼容和互操作性。相比之下，緊組合模型充分考慮了不同GNSS系統(tǒng)間的融合處理，減少了參考衛(wèi)星的個數(shù)，增加了觀測值的個數(shù)，提高了定位模型強度以及相對定位的精度、可靠性和可用性。

2.1.2 GPS坐標時間序列

GPS在地震監(jiān)測研究中的基礎應用主要是通過高精度、自動化、快速甚至實時的GPS觀測數(shù)據(jù)處理，得到地殼形變的時間序列和空間分布信息。通過長時間的GPS坐標時間序列處理，進行數(shù)據(jù)補缺、粗差剔除、噪聲分析、功率譜分析、共模誤差分析、階躍修復和震后馳豫形變等手段，可獲取前兆信息，進而研究地殼運動機理或模型。

為了對GPS坐標時間序列進行分析，需要建立相應的時間序列模型。Davis等（2012）提出的對數(shù)模型目前被廣泛應用于GPS坐標時間序列分析領域，其數(shù)學模型如下（Blewitt，Lavallée，2002;Langbein，2008;Montillet，Yu，2015;黃立人，2006;賀小星，2013;侯曉偉，2018;許家琪，2019）：

y（ti）=a+bti+csin（2πti）+dcos（2πti）+esin（4πti）+fcos（4πti）+∑njj=1giH（ti-Thj）+∑nhj=1hjH（ti-Thj）ti+vi（1）

式中：a為截距;b為線性速率;csin（2πti）+dcos（2πti）+esin（4πti）+fcos（4πti）為周年半周年變化; ∑njj=1giH（ti-Thj）為同震階躍; ∑nhj=1hjH（ti-Thj）ti為非同震階躍（儀器變更）;vi為觀測噪聲。該公式包括了測站的線性速率b（該值即為后面比較的速度場的值）、周年變化半周年變化、同震階躍和非同震階躍。

基于上述對GPS時間序列的分析，扣除粗差，并對GPS時間序列中的階躍（地震同震、儀器更換等原因造成的）進行探測和剔除，最終得到相對“干凈”的GPS站點位置時間序列（圖1）。

2.1.3 GPS基線時間序列

GPS站點數(shù)據(jù)解算還可獲得基線時間序列（圖2）。由于GPS基線時間序列能在一定程度上削弱系統(tǒng)誤差和共模噪聲的影響，可以明顯提高形變場觀測精度，因此比單個測站坐標時間序列能更靈敏地反映周圍構(gòu)造環(huán)境的變化，連續(xù)觀測點日誤差僅為約1 mm，在斷層兩側(cè)成對設立GPS可以用來有效監(jiān)測斷層內(nèi)、特別是脆韌轉(zhuǎn)換帶內(nèi)的細微形變變化，捕捉斷裂帶上與構(gòu)造運動相關(guān)的精確時空信息（王敏等，2008;張風霜等，2011）。王敏等（2008）利用跨鮮水河斷裂的兩對GPS連續(xù)觀測點資料，獲得了跨鮮水河斷裂高精度形變場（誤差約1 mm）的動態(tài)演化軌跡。該軌跡清晰地反映了鮮水河斷裂乾寧段和道孚段形變場的明顯不同，揭示當前斷裂南段轉(zhuǎn)換層強度可能高于北段，對應于更長的發(fā)震周期。

2.2 中國大陸構(gòu)造變形速度場與應變率場

自1988年中德合作滇西地震預報實驗場開展GPS地殼形變監(jiān)測以來，獲取中國大陸構(gòu)造變形的完整圖像始終是研究人員們最基礎、最重要的科學目標。中國大陸地殼運動與變形監(jiān)測經(jīng)歷了監(jiān)測區(qū)域由小到大、站點由稀疏到相對密集的一個發(fā)展過程，與之相對應的是地殼運動速度場的不斷完善和精化（王敏，沈正康，2020）。

2.2.1 水平速度場

早期的GPS測站十分稀疏，區(qū)域控制范圍極其有限，人們對中國大陸整體的構(gòu)造形變特征認知不夠。在前人工作基礎上，王琪等（2001，2002，2006）和Wang（2001a）首次利用國內(nèi)及周邊國家地區(qū)的229個GPS觀測站資料進行統(tǒng)一解算，獲得了更為精細的中國大陸地區(qū)地殼運動水平速度場，在國際上產(chǎn)生重要影響。馬宗晉等（2001）利用1998—2001年基準站和基本站觀測約束了中國大陸一級構(gòu)造塊體的運動特征。王敏等（2003）通過分析中國地殼運動觀測網(wǎng)絡GPS數(shù)據(jù)特別是1999年與2001年區(qū)域網(wǎng)數(shù)據(jù)，初步得到了中國大陸地殼運動速度場，并用統(tǒng)計分析的方法從高密度臺站速度場中區(qū)分出9個獨立活動塊體和2個廣泛形變帶，求出活動塊體剛體運動歐拉極和相鄰塊體間相互運動速率。牛之俊等（2005）利用3期及其他一些區(qū)域流動觀測結(jié)果共1 199個觀測站點，通過GAMIT軟件得到中國大陸現(xiàn)今地殼運動速度場。該速度場清晰地揭示了不同地塊的變形方式差異，最明顯的特征是，青藏高原和川滇地區(qū)繞東喜馬拉雅構(gòu)造結(jié)順時針旋轉(zhuǎn)，GPS觀測的速度矢量在西藏東部向北運動，到川西一帶轉(zhuǎn)為向SEE向，到云南地區(qū)逐漸轉(zhuǎn)為向SSE、SW方向運動。此外中國西部速度矢量由南向北逐漸遞減，表明其青藏高原的地殼縮短吸收了絕大部分印度和歐亞板塊之間的相對運動，剩余部分則被天山及以北的地殼縮短所吸收。Gan等（2007）用青藏高原及周邊726個GPS站，對青藏高原現(xiàn)今地殼運動給出了精準的描述和推斷，認為整個青藏高原的剛性旋轉(zhuǎn)可以被描述為相對歐亞框架下的（24.38±0.42）°N、（102.37± 0.42）°E，（0.7096± 0.0206）°/Ma歐拉矢量，剛性旋轉(zhuǎn)組件吸收了至少50%的來自印度板塊向東北的推力，并主導了北部的向東擠壓。

由于2008年陸態(tài)網(wǎng)絡的實施，GPS站點密度得到極大提升，我國的GNSS觀測資料日趨豐富，中國大陸構(gòu)造變形圖像由最初的幾十個點發(fā)展至目前的幾千個點。李強等（2012a，b）以網(wǎng)絡工程1 056個區(qū)域站1998—2007年觀測資料及陸態(tài)網(wǎng)絡2 056個區(qū)域站2009和2011年兩期觀測資料為基礎，剔除其中位移異常站，實際選取總計1 979個測站，利用GAMIT、GIPSY、Bernese、Panda獨立計算等權(quán)合并，盡可能增強觀測結(jié)果的可靠性，繪制出了空間分辨率高，構(gòu)造變形特征更為豐富、翔實的中國大陸地區(qū)地殼運動水平速度場。Zheng等（2017）利用1991—2015年2 574個GPS速度場測量數(shù)據(jù)，給出了印度—歐亞大陸最完整、準確和最新的速度場。Xu和Stamps（2019）等利用1996—2017年陸態(tài)網(wǎng)絡GPS數(shù)據(jù)，產(chǎn)出了可靠的速度場，并據(jù)此估計了中國大陸的地震活動率。Wang等（2020）對近25年超過4 000個GPS觀測點的數(shù)據(jù)資料進行了嚴格處理，得到了當前最新、最精細的中國大陸地殼運動速度場，更進一步地印證了印度與歐亞板塊相對運動被青藏高原變形所吸收，變形發(fā)生在整個青藏高原，而不僅在少數(shù)斷裂帶上，呈現(xiàn)近南北向縮短，東西向拉長，繞喜馬拉雅東構(gòu)造結(jié)和秦嶺構(gòu)造結(jié)發(fā)生明顯的順時針旋轉(zhuǎn)特征。

2.2.2 垂向速度場

在垂向速度場方面，甘衛(wèi)軍課題組利用青藏高原及其周邊十多年的流動和連續(xù)GPS資料，對青藏高原現(xiàn)今地殼運動，尤其是垂向地殼運動進行研究，并分析探討青藏高原現(xiàn)今隆升狀況與長期地形的相關(guān)性及其所反映的地球動力學含義，揭示了青藏高原相對其北部穩(wěn)定塊體的整體隆升趨勢（Liang et al，2013;梁詩明，2015）。Hao等（2016）通過對青藏高原東南部連續(xù)GPS站和GRACE數(shù)據(jù)聯(lián)合分析，確定了青藏高原東南部的垂直地殼運動速率，揭示了青藏高原東南部正在經(jīng)歷隆升，而川南—云南南部相對于北部地區(qū)正在下沉的現(xiàn)象，這也是地殼均勻伸展假設所不能解釋的。Zhan等（2017）利用2010—2016年235個GPS連續(xù)站的數(shù)據(jù)分析了大空間尺度的中國大陸地殼垂向運動，又以云南和天津兩個典型地區(qū)為重點研究區(qū)域進行了重點分析，對中國大陸垂向運動以及周期性特征和區(qū)域性差異給出了定量描述。

隨著網(wǎng)絡工程、陸態(tài)網(wǎng)絡的建設，我國的GPS觀測資料日趨豐富，中國大陸構(gòu)造變形圖像由最初的幾十個點發(fā)展至目前的幾千個點，為定量、精細研究塊體運動、斷層活動提供了數(shù)據(jù)基礎，GPS觀測資料已全面、系統(tǒng)、精確地展示出中國大陸地區(qū)地殼構(gòu)造變形的方式及幅度。

2.2.3 應變率場

GPS速度場反映的是構(gòu)造活動的運動學特征，應變率場給出的是變形樣式和強度，在一定程度上反映了構(gòu)造應力場特征。在GPS點位分布不均勻的地區(qū)，利用GPS速度場計算的應變率場能更好地描述構(gòu)造變形特征（王敏，沈正康，2020）。另外，不同的應變矢量，可以反映不同方向的變形分布（武艷強等，2020）。應變率場研究對于地震危險性評估具有重要意義，利用GPS資料獲取地殼運動與應變積累一直是國內(nèi)外重視的地震中長期預測的技術(shù)途徑（王靜，2019）。計算應變率場的方法有很多（Shen et al，1996;李延興等，2001;楊少敏等，2012;石耀霖等，2006;沈正康等，2003;黃立人，王敏，2003;江在森等，2003，2010;武艷強等，2009，2011），使用相同或類似的GPS資料求應變，不同研究者給出的結(jié)果會有顯著差別（石耀霖，朱守彪，2006）。不同應變率計算方法有其各自的特點，在具體使用中，需要根據(jù)區(qū)域大小、地質(zhì)資料、GPS站點數(shù)量及空間分布選擇合適的方法。

根據(jù)動態(tài)的應變率場圖像再計算即可得到主應變、面膨脹和剪應變，這對研究區(qū)域應變積累過程具有重要意義。通過剪應變可以比較直觀、清楚地顯示出，中國大陸內(nèi)部不同構(gòu)造單元對周邊板塊和板內(nèi)深部構(gòu)造動力共同作用的變形響應。中國大陸地區(qū)最大剪應變率以南北地震帶為界，呈西強東弱態(tài)勢，中國大陸內(nèi)部剪應變率最大的三個區(qū)域依次是喜馬拉雅、昆侖山斷裂帶、阿爾金斷裂帶，藏南、羌塘地塊中西部部分地區(qū)、昆侖山中部至鮮水河—安寧河—小金斷裂帶一帶、東天山等屬最大剪應變率高值區(qū)（李延興等，2003;江在森等，2001，2003;沈正康等，2003;王敏，沈正康，2020）。另外，安寧河、小江斷裂帶、鄂爾多斯東北角、祁連山和天山的剪應變值也相對較高。

通過GPS推算的現(xiàn)今地殼應變場的空間變化與大陸強震具有較好的對應關(guān)系（張培震等，2003）。1995年以來發(fā)生的MS≥6地震的震中一般都分布在剪應變高值區(qū)，尤其是2001年昆侖山口西MS8.1地震，正好發(fā)生在面膨脹應變率和最大剪應變率高值區(qū)邊緣的梯度帶上（Lin et al，2002）。隨著觀測精度的進一步提高和站點密度的增加，GPS必將在地殼構(gòu)造變形和地球動力學研究以及地震預測中發(fā)揮越來越大的作用。

2.3 中國大陸活動地塊運動定量化

中國大陸巖石圈新生代和現(xiàn)代構(gòu)造變形的最顯著特征是巨大的晚第四紀活動斷裂十分發(fā)育，將中國大陸切割成為不同級別的活動塊體（鄧起東等，2002）?；顒拥貕K邊界構(gòu)造活動強烈，內(nèi)部相對穩(wěn)定，絕大多數(shù)強烈地震都發(fā)生在地塊邊界的活動構(gòu)造帶上。根據(jù)晚第四紀構(gòu)造變形、地震活動和地球物理場的差異，可以將中國大陸分成若干活動地塊，GPS所揭示的中國大陸現(xiàn)今地殼運動場清晰地表現(xiàn)出了分塊特征，不同的活動地塊具有不同的水平運動和變形方式（黃立人，王敏，2003）。

丁國瑜和盧演儔（1986）根據(jù)地震資料、活動構(gòu)造資料首次定量研究了我國大陸內(nèi)部各個塊體的相對運動和總的變形情況，獲得了百余個觀測段上斷層滑動平均速率矢量的水平投影值，從地質(zhì)學角度揭示了斷層的長期滑動速率。張培震等（2002b，2003）結(jié)合測震、活動構(gòu)造及GPS資料，將中國大陸及其鄰區(qū)劃分為22個活動地塊，并定量描述了每個地塊運動的速率大小和方位。李延興等（2003）利用1 598個GPS觀測站資料，從地殼運動與應變的角度，定量計算了各地塊的運動參數(shù)和應變參數(shù)。

隨著GPS觀測資料的不斷積累，尤其是號稱“世界屋脊”的青藏高原地區(qū)GPS觀測資料的增多，為研究其構(gòu)造變形模式提供了有力的數(shù)據(jù)支持。

國外學者應用已發(fā)表的GPS速度場數(shù)據(jù)，從塊體運動的角度進一步量化了青藏高原以及周邊滑動斷裂帶的運動速率?；赯hang等（2004）發(fā)表的青藏高原地區(qū)GPS速度場數(shù)據(jù)，Meade（2007）將其與地質(zhì)觀測數(shù)據(jù)整合到內(nèi)部一致的三維區(qū)域尺度模型中，對彈性上地殼的彈性行為進行解釋，繪制出了青藏高原及周邊地區(qū)地殼運動各地塊邊界滑動速率圖;Thatcher（2007）同樣選取了其中349個GPS站的速度，利用準剛性塊模型，將青藏高原及其周圍地區(qū)劃分成11個塊體。

GPS的出現(xiàn)與發(fā)展，為定量化研究現(xiàn)代板塊內(nèi)部活動構(gòu)造提供了前所未有的機遇，其與地震地質(zhì)資料的結(jié)合發(fā)展，豐富和完善了現(xiàn)代板塊內(nèi)部活動構(gòu)造理論，為定量研究地震危險性分析提供了依據(jù)（Shen et al，2005;王閻昭等，2008;楊少敏等，2008;Zhao et al，2017a;Wang et al，2017;Zheng et al，2017）。

2.4 地震研究

GPS可以準確測定地震引起的永久的靜態(tài)同震形變場，為研究地震的同震破裂分布、斷層幾何結(jié)構(gòu)提供重要約束，彌補遠場地震波形對發(fā)震斷層幾何結(jié)構(gòu)不敏感的缺點。此外，高精度、高頻率的GPS監(jiān)測擴充了強地面運動觀測的頻帶范圍，不存在振幅飽和、曲線漂移等缺點，已成為地震儀的有效補充。GPS觀測也被廣泛用于震后形變監(jiān)測和震后機理研究等方面，而以震后形變?yōu)榧s束，可以探討斷層的余滑性質(zhì)、下地殼/上地幔的流變性，為斷層的運動學研究和區(qū)域地球動力學研究提供依據(jù)，是對巖石力學實驗、冰后回彈等手段探測結(jié)果的檢驗和補充（刁法啟，2011;劉剛等，2020;趙斌等，2020）。

2.4.1 同震形變場

中國境內(nèi)的GPS臺網(wǎng)已經(jīng)成功觀測到中國大陸及其周邊地區(qū)多次地震的同震形變場，如2004年印尼蘇門答臘發(fā)生MW9.3地震（王敏等，2006）、2008年汶川MS8.0地震（Wang，2011;伍吉倉等，2020）（圖3）、2011年日本MW9.0宮城地震（王敏等，2011）、2015年尼泊爾MS8.1地震（李瑜等，2015;趙斌等，2015）等大地震，以及2010年青海玉樹MS7.1地震（Meng et al，2013）、2013年四川蘆山MS7.0地震（Jiang et al，2014）、2014年新疆于田MS7.3地震（梁洪寶等，2018）、2014年云南魯?shù)镸S6.5地震（魏文薪等，2018）、2015年新疆皮山MW6.4地震（He et al，2016）、2016年新疆阿克陶MW6.6地震（Li et al，2019）、2017年四川九寨溝MS7.0地震（Zhao et al，2018;王閱兵等，2018;Liu et al，2019）等。值得一提的是GPS在汶川地震同震三維形變場的測定、同震滑動分布模型的反演中發(fā)揮了至關(guān)重要的作用，大大促進了對汶川地震發(fā)震機理的認識（Wang，2011b;顧國華，王武星，2011;殷海濤等，2010;趙靜等，2018）。

2.4.2 高頻GPS應用

GPS與強震儀同址觀測，可利用加速度和位移的關(guān)系，來彌補兩種觀測手段的不足，即通過高頻GPS來矯正強震記錄，以強震記錄來提高高頻GPS的采樣率和精度，達到聯(lián)合反演地震波的影響過程，這在近年應用較為廣泛。Larson等（2003）最早成功利用1 Hz GPS觀測數(shù)據(jù)進行研究，復原了2002年德納利（Denali）MW7.9地震的遠場地震波，并與強震儀記錄比較，認為1 Hz GPS觀測資料可以聯(lián)合地震觀測資料反演研究地震破裂過程。此后，外國學者Ji等（2004）、Bock等（2004）、Miyazaki等（2004）、Emore等（2007）、Bilich等（2008）、Davis和Smalley（2009）和國內(nèi)學者殷海濤等（2010）、Shi等（2010）、方榮新（2010）、張小紅等（2016）、申文豪等（2019）也先后利用高頻GPS數(shù)據(jù)進行了相關(guān)研究，推動了高頻GPS在地震研究中的廣泛應用。Geng等 （2019，2020）將并置的高頻GNSS、加速度計和陀螺儀3種儀器進行融合，以期發(fā)揮每種儀器的優(yōu)勢。融合3種觀測數(shù)據(jù)的地震大地測量學方法，在保證觀測穩(wěn)定性的條件下，將觀測近場形變精度提高至毫米量級，這對于大地震震源機制的研究，乃至地震預警與烈度速報工程都具有重要意義。

2.4.3 GPS震后形變監(jiān)測

隨著GPS觀測臺網(wǎng)的建設和發(fā)展，國內(nèi)也積累了非常豐富的震后形變觀測資料，為震后形變監(jiān)測和震后機理研究提供了豐富的數(shù)據(jù)基礎，可模擬震后應變應力時空發(fā)展過程，為地震預測提供參考。2001年昆侖山口西MW7.8地震是網(wǎng)絡工程建成后發(fā)生的首次大地震（趙斌等，2020）。通過網(wǎng)絡工程GPS測站和震后加密觀測，國內(nèi)學者首次使用GPS技術(shù)大面積高密度監(jiān)測地震同震位移和震后斷層蠕動（任金衛(wèi)，王敏，2005），觀測到了昆侖山斷裂帶南、北兩側(cè)震后形變的不均一性（Ryder et al，2010;賀鵬超等，2018;Liu et al，2019）。2008年汶川MS8.0地震后，GPS觀測站的數(shù)量較之前有很大的提升，觀測到明顯的地震后變形，垂直于斷層方向的形變在2013年后衰減到幾乎為零，直到2013年5月的觀測中，明顯的右旋運動并沒有明顯衰減（Xu et al，2014）。2015年尼泊爾MS8.1地震后，Zhao等（2017b）利用GPS觀測資料，首次提出藏南地區(qū)幔黏彈性松弛在震后早期也在起作用，并分析得到印度—藏南黏滯系數(shù)存在明顯的橫向差異，與青藏高原東緣具有相似性（Huang et al，2014;趙斌等，2020）。

2.5 廣泛的影響

GPS觀測技術(shù)從20世紀90年代初發(fā)展至今，其應用在中國大地上已遍地開花，不論是在軍事測繪保障、交通運輸導航、災害天氣預報方面，還是在地學科學研究、地震預測預報等領域，都發(fā)揮了巨大的作用，并產(chǎn)生了廣泛而深刻的影響（姜衛(wèi)平，2017;陳增強等，2000;王小亞，2002;尹繼堯等，2007;魏明耀等，2010;徐韶光等，2011）。

GPS的出現(xiàn)，在時空尺度上極大地提升了大地測量監(jiān)測地球表層運動與變形的能力。全國尺度的GPS觀測網(wǎng)絡（網(wǎng)絡工程、陸態(tài)網(wǎng)絡等）的建成并投入使用，實現(xiàn)了對板塊及其邊界帶運動的監(jiān)測，建立了空間大地測量板塊運動模型;區(qū)域性和板內(nèi)重要斷裂帶運動的監(jiān)測和研究，提高了對板緣地震和板內(nèi)地震構(gòu)造環(huán)境的認識，推動了大陸動力學研究的進展。時空監(jiān)測能力的提高，大大增強了動力學研究的高精度地殼運動觀測約束，使位錯反演時所依據(jù)的地球模型更接近真實的地球，斷層模型更加復雜多樣，適應了不同時空尺度大地測量反演的需要，推動了位錯理論的發(fā)展與應用，使我國在相關(guān)領域的研究躍入國際先進行列。

3 面臨的問題

3.1 地震監(jiān)測預報難題尚未解決

GPS觀測因建立了全球統(tǒng)一參考基準而保證了大空間尺度觀測結(jié)果具有可比性及高精度特性，對地震同震形變、震后形變以及震間形變等各個方面的研究起到非常重要的作用（周碩愚等，2006，2007，2015;吳云等，2003;顧國華等，2003，2012;薄萬舉，2013;張曉亮等，2006;冒愛泉等，2008）。然而，地震監(jiān)測預報遇到的困難遠遠超出人們的預期，是世界性的科學難題，目前尚處于初期的科學探索階段（陳運泰，2009）。相比于傳統(tǒng)大地測量技術(shù)，GPS能夠獲得多尺度地殼運動和形變動態(tài)信息，但與地震預測對孕育過程信息的需求相比仍存在較大的差距（江在森等，2020）。從國內(nèi)到國外，GPS在地震預測方面的應用都面臨相似的困境。在地震預測研究方面，國內(nèi)外的研究實例均表明基于震間期GPS資料反演的閉鎖程度/滑動虧損分布與大地震發(fā)生地點具有對應關(guān)系，當GPS觀測資料的空間分辨率滿足要求時，此類工作對于強震地點和震級預測具有指示意義。這些回溯性研究為GPS應用于中長期地震預測指明了方向。

3.2 基礎研究缺乏突破

地震工作者基于震前GPS資料針對強震孕育晚期的研判一直在進行著有益的探索（武艷強等，2020）：如基于強震前GPS形變場空間分布與強震發(fā)生地點關(guān)系的探索（剪切變形分布與走滑型強震的關(guān)系、正應變分布與傾滑型強震的關(guān)系等）;大地震同震及震后影響的顯著弱響應區(qū)與后續(xù)強震的關(guān)系探索;強震前地殼變形模式與同震應變釋放模式的非對稱特征探索等。但目前此類研究多為現(xiàn)象總結(jié)和定性機理分析，尚缺乏有效的數(shù)理模型支持（周碩愚等，2006）。

基于GPS時序資料進行震源異常的探索一直被廣大地震工作者所關(guān)注，多位學者研究認為基于GPS時序結(jié)果，在2011年日本“3·11”地震前觀測到了幾年尺度的偏離原有趨勢的異?，F(xiàn)象（回溯性研究）（Ozawa et al，2012;Yokota，Koketso，2015）。另外，基于GPS時序資料的強震短臨異常探索也在積極探索中，比如GPS方向角指數(shù)方法、變形參數(shù)時序變化與強震關(guān)系等。目前此類研究多為零散的現(xiàn)象總結(jié)，尚缺乏合理有效的定性或定量數(shù)理模型支持。

綜合而言，基于GPS資料的強震地點預測取得了一些進展（江在森，武艷強，2012），具有一定的理論模型支撐，亟需加大科研投入，發(fā)展適合于中國的板內(nèi)地震孕育模型、復雜斷層系統(tǒng)強震孕育模型、多活動地塊組合的強震孕育模型。另一方面，基于GPS資料對于強震孕育晚期判斷及短臨預測方面，積累了一些觀測現(xiàn)象，除了繼續(xù)積累震例外，亟需加大科研支撐，盡快發(fā)展起有效的解析或數(shù)值模型。

3.3 站點密度有待提高

現(xiàn)今活動斷層滑動速率和閉鎖深度的判定是進行地震危險性地點預測的基礎。在我國西部大多數(shù)活動斷裂及其附近，GPS站點的間距在80～150 km，而斷層滑動引起的地表位移變化大的地方一般在斷層兩邊100 km的范圍內(nèi)。因此，以現(xiàn)有的GPS站點密度無法可靠地反演斷層參數(shù)，判斷斷層活動的性質(zhì)，因此也無法對地震危險性做出可靠的判斷。在幾次大地震之前，連續(xù)GPS站由于距離震中較遠沒有發(fā)現(xiàn)異常，而區(qū)域站發(fā)現(xiàn)了明顯的異常。目前，國際上應用于地震監(jiān)測的最具代表性的地殼形變站網(wǎng)是美國（南加州地區(qū)）和日本的GPS觀測網(wǎng)。其中美國（南加州地區(qū)）采用“重點地區(qū)密集觀測”策略，觀測站沿重點斷裂帶周邊布設，間距達2～5 km;日本采用“全域網(wǎng)格式均勻分布”策略，在陸地范圍內(nèi)，觀測站平均間距約為17 km，同美國、日本的GPS觀測網(wǎng)相比，國內(nèi)在站間距方面還有很大差距。

3.4 GPS解算精度和時效有待提高

GPS觀測極大地提高了地殼運動水平方向監(jiān)測能力，但受限于GPS系統(tǒng)設計與當前數(shù)據(jù)解算水平，垂向運動精度、數(shù)據(jù)解算時效和非構(gòu)造因素剔除等諸方面仍相當薄弱。

4 發(fā)展方向

4.1 強化地震預測研究目標導向

（1）強化中長期地震預測，完善GPS在發(fā)震地點和震級預測的方法和應用。目前，GPS在地震中長期預測中的應用都獲得廣泛的認同（江在森等，2020）。研究人員一致認為，GPS的觀測與彈性位錯理論吻合很好，成功地解釋了地震的孕育和破裂過程。GPS研究獲得的斷裂滑動速率、閉鎖程度、地震矩積累程度等可應用于地震的中長期預測，能較好地確定中長期地震的震級和地點。

（2）加強GPS用于短臨預測的實踐探索與研究。GPS在地震中長期預測中的應用已經(jīng)基本獲得大部分認可，但是在地震短臨預測中的應用還存在較大的分歧。地震前兆是有的，問題在于是否能觀測、記錄到，即使觀測到了，又是否能識別（馬瑾，2016）。有些專家認為GPS由于時空分布不足以及精度的問題，開展短臨預測還存在困難，但有些專家對此呈樂觀態(tài)度。GPS應用于地震短臨預報缺乏具有信服力的證據(jù)，但值得大力探索。目前，GPS短臨預測研究還相對較少（Chen et al，2014;馮蔚，2016;Feng et al，2016;Bedford et al，2020），需要加強。短臨預測主要依靠GPS連續(xù)觀測資料，這對GPS連續(xù)觀測數(shù)據(jù)的處理和應用方法提出了新的挑戰(zhàn)，需要探索、發(fā)展GPS在短臨預測跟蹤中的數(shù)據(jù)處理方法和應用方法的研究。

（3）強化斷層滑動行為的精細化研究，分析地震危險性?；诿芗腉PS觀測資料（包括速度場、基線或應變率），以震間GPS運動速度場為約束，研究主要活動斷裂滑動速率、應變和地震能量積累，并根據(jù)地震目錄統(tǒng)計斷裂帶的地震矩釋放情況，從而分析評估區(qū)域地震危險性。重點加強與地震地質(zhì)（古地震考察）、地震精定位等相關(guān)研究的交叉融合。

（4）提高斷層滑動瞬態(tài)變化過程的監(jiān)測能力?；诿芗倪B續(xù)GPS觀測站，構(gòu)建跨斷層基線及觀測網(wǎng)，識別并提取出斷層滑動行為隨時間的變化特征，進而研究斷層瞬態(tài)變形過程與地震活動的關(guān)系。

（5）GPS分析研究要從地表到地下，從運動學到動力學轉(zhuǎn)換。GPS觀測到的是地表的運動，要分析預測地下的地震活動，需要加強深部應變、應力變化分析，甚至要從運動學研究轉(zhuǎn)換到動力學研究。

（6）加強GPS與InSAR、地震波等多元數(shù)據(jù)、多學科融合應用。各學科的數(shù)據(jù)都有自己的優(yōu)勢和劣勢，GPS可以觀測大尺度的長期構(gòu)造變形以及近場形變數(shù)據(jù)，InSAR可獲取近場形變的面場變化，并對高程變化敏感（Gong et al，2019），加強GPS與InSAR資料的融合、提高垂向位移信號的信噪比應該是未來努力的方向，而統(tǒng)計地震資料獲得的地震平靜期和空區(qū)信息對預測很有幫助。綜上認為，將GPS與InSAR、地震波等多元數(shù)據(jù)融合應用，可以降低地震分析預測的不確定性。

4.2 加強基礎研究和基礎性工作

要從明確定位、重點工作領域、規(guī)劃布局和推進舉措等方面加強GPS用于地震監(jiān)測預報的基礎研究工作。一是明確基礎研究工作定位。要以支撐地震監(jiān)測預報為導向，開展應用目標明確的基礎研究工作。二是明確基礎性工作重點領域。圍繞地震成因理論、地震觀測技術(shù)和強震潛在危險性等關(guān)鍵科學問題，注重長、中、短臨和不同空間尺度相結(jié)合，強調(diào)地質(zhì)學、地球物理學和大地測量等多學科的綜合，持續(xù)推進GPS地震學基礎性工作，為基礎研究和技術(shù)創(chuàng)新夯實工作基礎。三是加強謀劃布局。開展GPS地震應用工作研討，梳理凝練未來基礎研究和基礎性工作的重點領域方向和科學選題，研究提出未來項目、平臺、設施布局建議。四是明確推進舉措。出臺進一步加強GPS地震學科研基礎工作的指導意見，探索建立跨學科、面向基礎研究的指導機構(gòu)，強化基礎性工作人才隊伍、實驗室建設和項目支持，推進科技數(shù)據(jù)資源共享、共用。

4.3 統(tǒng)籌數(shù)據(jù)資源，優(yōu)化GPS監(jiān)測布局

（1）充分利用連續(xù)站資源，適當減少東部地區(qū)的流動觀測。近幾年GPS流動監(jiān)測得到顯著加強，主要構(gòu)造帶和重點地震監(jiān)測區(qū)域基本實現(xiàn)了每年一期的復測。但從目前數(shù)據(jù)處理情況看，對于東部地區(qū)地殼運動與變形較小的地區(qū)，以當前實際3 mm左右的GPS流動觀測精度來看，這樣密集的重復觀測意義不大，特別是華北、東北、華南地區(qū)采用大量的土層觀測墩，其穩(wěn)定性存在顯著的沉降與傾斜，難以真實反映地殼構(gòu)造運動的實際狀態(tài)。在全國GPS資源共享基礎上，進行共享連續(xù)站點數(shù)據(jù)質(zhì)量篩查，挑選出能夠用于地殼運動監(jiān)測的站點，可代替部分流動觀測。部分區(qū)域如華北和華南的連續(xù)站點密度比陸態(tài)網(wǎng)絡流動站點密度還要大，因此，充分利用連續(xù)站觀測資料，適當減少這些區(qū)域流動觀測，不僅省時省力，而且可獲取連續(xù)觀測數(shù)據(jù)，能較好地分析該區(qū)域形變。應對歷史上發(fā)生過強震的重大斷裂（如郯廬斷裂）加強數(shù)據(jù)分析，必要時增補流動觀測。

（2）加強西部地區(qū)觀測密度?，F(xiàn)有的GPS觀測網(wǎng)在布局上存在東西部布網(wǎng)觀測的不均衡性，東部相對密集，而地殼構(gòu)造運動強烈的廣大西部地區(qū)相對稀疏?，F(xiàn)有100～300 km的基準站密度，80～150 km流動站的密度，要實現(xiàn)對地震發(fā)生的塊體邊界帶和主要斷層的密集連續(xù)監(jiān)測，尤其是西部地區(qū)還遠不夠，有必要持續(xù)加大西部流動GPS觀測或連續(xù)觀測站的密度。建議進一步統(tǒng)籌、優(yōu)化現(xiàn)有資源配置，通過加密建站在短期內(nèi)實現(xiàn)站點間距50～80 km，遠期內(nèi)（5～10年）西部地區(qū)站點間距達到30～50 km，針對重大斷裂帶最好達到10 km間距。若空間分布密度在短時間內(nèi)不能實現(xiàn)，可以在西部地區(qū)適當增加觀測頻次，至少每年觀測一期，重點危險區(qū)可一年進行2～3次復測。

（3）構(gòu)建活動塊體邊界帶的綜合觀測體系。大陸強震主要發(fā)生在活動塊體邊界帶上，單純依靠GPS技術(shù)不可能捕捉到全部地震孕育釋放的信息，應建立GPS、InSAR、定點形變、重力、測震、流體、地磁等綜合觀測體系，甚至進行多學科同址觀測，以期發(fā)揮多種技術(shù)優(yōu)勢，增加約束信息，進行相互比對驗證。定點形變儀、地震儀、重力儀的觀測資料的應用，還可以拓寬GPS監(jiān)測的時間尺度，獲取三維動態(tài)形變場，提高地殼形變分辨能力，對開展斷層構(gòu)造運動、地震能量積累、地震破裂、震后形變、斷層蠕滑等典型動力學過程跟蹤研究具有重要意義。

（4）加強地震重點危險區(qū)細部觀測，獲取精細變形特征。根據(jù)中國大陸地震重點危險區(qū)的研判結(jié)果，地震危險信度和緊迫程度高或相對較高的危險區(qū)，連續(xù)GPS站間距可適當加密至25 km左右。對西部十年尺度危險區(qū)，要獲取斷層閉鎖特征、應變積累速率，可以采用流動站加密觀測為主，輔以新建臨時固定觀測站。喜馬拉雅東構(gòu)造結(jié)，由于臨近南北地震帶，其動態(tài)作用與川滇地區(qū)強震的孕育、發(fā)生密切相關(guān)，需要重點強化觀測。喜馬拉雅西段存在M9強震風險，會對整個青藏高原及中國大陸的動力學環(huán)境產(chǎn)生重大影響，也需要強化對該區(qū)的連續(xù)監(jiān)測能力。

（5）在“十四五”期間，推進GPS臺站加密建設。利用“一帶一路”和“川滇國家地震科學實驗場”等重大項目的推進實施，實現(xiàn)全國和主要是重點區(qū)域GPS連續(xù)臺站觀測密度的顯著提升。

4.4 重視基礎工作，加強GPS解算精度和時效

首先，要開發(fā)自主產(chǎn)權(quán)的GPS數(shù)據(jù)處理軟件。目前的GPS數(shù)據(jù)處理主要依賴于國外的GAMIT、Bernese、GIPSY等軟件。從長遠發(fā)展看，我國需要自主研發(fā)一套能夠進行GPS數(shù)據(jù)處理的軟件，以提升國內(nèi)GPS數(shù)據(jù)處理分析的水平;同時，通過數(shù)據(jù)處理軟件的研發(fā)，培養(yǎng)一批GPS數(shù)據(jù)處理分析的專業(yè)技術(shù)人才，探索提高GPS數(shù)據(jù)解算精度、時效和剔除非構(gòu)造影響因素的方法。

其次，加強對全國垂直形變場研究。加強對GPS垂向形變分量數(shù)據(jù)的研究，結(jié)合區(qū)域水準復測資料，分析GPS應用于地殼垂向運動監(jiān)測的可行性;同時還要分析不同區(qū)域地殼垂向活動與地震活動之間的關(guān)系。

此外，還要加強GPS從業(yè)人員的素質(zhì)培養(yǎng)。GPS觀測資料解算過程相對復雜，需要從業(yè)人員具備較扎實的專業(yè)背景、長期的技術(shù)儲備和穩(wěn)定的工作團隊。為此，需要建立相對穩(wěn)定的GPS數(shù)據(jù)處理、分析預報的從業(yè)人員隊伍，制定教育培訓計劃，定期開展GPS相關(guān)業(yè)務的培訓，同時加強從業(yè)隊伍的工作團隊建設，在實際工作中通過團隊合作提高人員素質(zhì)。

5 結(jié)語

GPS對地觀測技術(shù)為大地測量和地殼運動監(jiān)測帶來了一場革命。全國尺度上短時間獲取地殼運動的定量結(jié)果，不僅對地球幾何形變的定量研究、地震動態(tài)形變過程監(jiān)測、地球動力學和斷裂帶微小地殼形變觀測研究等方面有著重要的推動作用，同時為地震危險區(qū)的預測、地震前兆的觀測和長、中、短、臨的地震預測提供了更好的方法和手段，使地震預報從經(jīng)驗預報向物理預報推進成為可能。

雖然當前GPS在地震預報上的應用主要還是提供中長期預報信息，但也在一些地震的預測和研判中發(fā)揮了巨大的作用，如2016年新疆阿克陶6.7級地震（39.27°N，74.04°E，震源深度10 km）和新疆呼圖壁6.2級地震（43.83°N，86.40°E，震源深度6 km）。相信隨著全國GPS數(shù)據(jù)資源的共享、十四五規(guī)劃、川滇試驗場等國內(nèi)一系列GPS項目的推進，GPS多手段的聯(lián)合，數(shù)據(jù)處理技術(shù)、基礎研究和針對地震預測基礎研究的突破（顧國華，2008），地震短臨預測水平將會得到大幅度提高。

參考文獻：

薄萬舉，章思亞，劉宗堅，等.2007.大地形變資料用于地震預測的回顧與思考[J].地震，27（4）：68-76.

薄萬舉.2013.GPS展示的中國大陸主要相對變形特征及強震活動研究[J].地球物理學進展，28（2）：599-606.

蔡宏翔，宋成驊，劉經(jīng)南.1997.青藏高原1993和1995年地殼運動與形變的GPS監(jiān)測結(jié)果分析[J].中國科學：地球科學，27（3）：233-238.

陳俊勇，王澤民，龐尚益，等.2001.論珠穆朗瑪峰地區(qū)地殼運動[J].中國科學：地球科學，31（4）：265-271.

陳俊勇.1997.中國地殼運動觀測網(wǎng)絡簡介[J].測繪通報，（2）：8-11，15.

陳運泰.2009.地震預測：回顧與展望[J].中國科學：地球科學，39（12）：1633-1658.

陳增強，胡瑞華，劉立嵩.2000.中國地殼運動觀測網(wǎng)絡在地球科學研究中的應用前景[J].地球科學進展，15（4）：421-425.

鄧起東，張培震，冉勇康，等.2002.中國活動構(gòu)造基本特征[J].中國科學：地球科學，32（12）：1020-1030.

刁法啟.2011.基于GPS觀測的同震、震后形變研究[D].北京：中國科學院研究生院.

丁國瑜，盧演儔.1986.對我國現(xiàn)代板內(nèi)運動狀況的初步探討[J].科學通報，31（18）：54-57.

方榮新.2010.高采樣率GPS數(shù)據(jù)非差精密處理方法及其在地震學中的應用研究[D].武漢：武漢大學.

馮蔚.2016.利用 GPS 資料分析震前地表位移階段性特征[D].北京：中國地震局地震預測研究所.

甘衛(wèi)軍，李強，張銳，等.2012.中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡的建設與應用[J].工程研究——跨學科視野中的工程，4（4）：16-23.

甘衛(wèi)軍，張銳，張勇，等.2007.中國地殼運動觀測網(wǎng)絡的建設及應用[J].國際地震動態(tài)，（7）：43-52.

顧國華，牛紅葉，孟國杰，等.2001.中國地殼運動觀測網(wǎng)絡GPS數(shù)據(jù)處理進展[J].全球定位系統(tǒng)，26（1）：5-13.

顧國華，王武星.2011.區(qū)域網(wǎng)GPS觀測得到的汶川大地震前后的地殼垂直運動[J].地震，31（3）：1-8.

顧國華，張晶，王武星.2003.GPS空間大地測量在地震預測研究中的最新進展[J].全球定位系統(tǒng)，28（1）：9-13.

顧國華，張晶，王武星.2005.中國地殼運動觀測網(wǎng)絡GPS觀測得到的地殼運動結(jié)果[J].測繪通報，（7）：2-8.

顧國華.2007.GNSS（GPS）觀測研究地殼運動的新進展[J].國際地震動態(tài)，（7）：9-15.

顧國華.2008.GNSS科學發(fā)展與前景[J].全球定位系統(tǒng)，33（4）：4-9，13.

顧國華.2012.地殼形變與地震前兆探索回顧和展望[J].地震，32（2）：22-30.

國家重大科學工程“中國地殼運動觀測網(wǎng)絡”項目組.2008.GPS測定的2008年汶川MS8.0級地震的同震位移場[J].中國科學：地球科學，38（10）：1195-1206.

賀鵬超，王敏，王琪，等.2018.基于2001年MW7.8可可西里地震震后形變模擬研究藏北地區(qū)巖石圈流變學結(jié)構(gòu)[J].地球物理學報，61（2）：531-544.

賀小星.2013.GPS臺站時間序列分析及其地殼形變應用[D].上海：東華理工大學.

侯曉偉.2018.高精度GPS用于地殼形變研究分析[D].西安：長安大學.

黃立人，馬青.2003.祁連山—河西走廊地區(qū)的現(xiàn)今水平形變[J].大地測量與地球動力學，23（4）：9-13.

黃立人，王敏.2003.中國大陸構(gòu)造塊體的現(xiàn)今活動和變形[J].地震地質(zhì)，25（1）：23-32.

黃立人.2006.GPS 基準站坐標分量時間序列的噪聲特性分析[J].大地測量與地球動力學，26（2）：31-33，38.

江在森，劉經(jīng)南.2010.應用最小二乘配置建立地殼運動速度場與應變場的方法[J].地球物理學報，3（5）：1109，1116-1117.

江在森，馬宗晉，張希，等.2001.青藏塊體東北緣水平應變場與構(gòu)造變形分析[J].地震地質(zhì)，23（3）：337-346.

江在森，馬宗晉，張希，等.2003.GPS初步結(jié)果揭示的中國大陸水平應變場與構(gòu)造變形[J].地球物理學報，46（3）：352-358.

江在森，武艷強，鄒鎮(zhèn)宇，等.2020.GNSS在中國大陸的地震預測應用研究進展與展望[J].中國地震，36（4）：693-707.

江在森，武艷強.2012.地殼形變與強震地點預測問題與認識[J].地震，32（2）：8-21.

姜衛(wèi)平.2017.衛(wèi)星導航定位基準站網(wǎng)的發(fā)展現(xiàn)狀、機遇與挑戰(zhàn)[J].測繪學報，46（10）：1379-1388.

賴錫安，黃立人，徐菊生.2004.中國大陸現(xiàn)今地殼運動[M].北京：地震出版社.

李強，韓宇飛，寧百齊，等.2012a.基于陸態(tài)網(wǎng)絡GPS數(shù)據(jù)的電離層空間天氣監(jiān)測與研究[J].地球物理學報，55（7）：2193-2202.

李強，游新兆，楊少敏，等.2012b.中國大陸構(gòu)造變形高精度大密度GPS監(jiān)測——現(xiàn)今速度場[J].中國科學：地球科學，42（5）：629-632.

李延興，黃城，朱文耀，等.2000.中國大陸水平運動特征與動力學分析[J].地震，20（S1）：29-34.

李延興，黃珹，胡新康，等.2001.板內(nèi)塊體的剛性彈塑性運動模型與中國大陸主要塊體的應變狀態(tài)[J].地震學報，23（6）：565-572.

李延興，楊國華，李智，等.2003.中國大陸活動地塊的運動與應變狀態(tài)[J].中國科學：地球科學，33（S1）：65-81.

李延興，趙承坤.1998.華北地區(qū)GPS監(jiān)測網(wǎng)建設，地殼水平運動與應力場及地震活動性的關(guān)系[J].中國地震，14（2）：116-125.

李瑜，劉靜，梁宏，等.2015.全球定位系統(tǒng)測定的尼泊爾MW7.8地震同震位移[J].科學通報，60（36）：3606-3616.

梁洪寶，武艷強，陳長云，2018.2014年新疆于田MS7.3地震同震位移及位錯反演研究[J].地球物理學報，61（12）：4817-4826.

梁詩明.2015.基于GPS觀測的青藏高原現(xiàn)今三維地殼運動研究[J].國際地震動態(tài)，（7）：40-42.

劉剛，喬學軍，王琪.2020.強震變形的地震大地測量監(jiān)測與研究[J].中國地震，36（4）：718-728.

劉序儼，郭逢英，林繼華，等.1999.福建沿海地殼運動與GPS測量結(jié)果初步分析[J].福建地震，19（1）：40-47.

馬瑾.2016.從“是否存在有助于預報的地震先兆”說起[J].科學通報，61（4）：409-414.

馬宗晉，陳鑫連，葉叔華，等.2001.中國大陸區(qū)現(xiàn)今地殼運動的GPS研究[J].科學通報，46（13）：1118-1120.

冒愛泉，張金華，胡伍生，等.2008.GPS地殼形變和前兆信息與地震預測[J].現(xiàn)代測繪，31（5）：3-6.

孟國杰.2001.基于全球定位系統(tǒng)的現(xiàn)今地殼運動與形變研究[D].北京：中國地震局地球物理研究所.

牛之俊，馬宗晉，陳鑫連，2002.中國地殼運動觀測網(wǎng)絡[J].大地測量與地球動力學，22（3）：88-93.

牛之俊，王敏，孫漢榮，等.2005.中國大陸現(xiàn)今地殼運動速度場的最新觀測結(jié)果[J].科學通報，50（8）：839-840.

任金衛(wèi)，王敏.2005.GPS觀測的2001年昆侖山口西MS8.1地震地殼變形[J].第四紀研究，25（1）：34-44.

申文豪，李永生，焦其松，等.2019.聯(lián)合強震記錄和InSAR/GPS結(jié)果的四川九寨溝7.0級地震震源滑動分布反演及其地震學應用[J].地球物理學報，62（1）：115-129.

沈正康，王敏，甘衛(wèi)軍，等.2003.中國大陸現(xiàn)今構(gòu)造應變率場及其動力學成因研究[J].地學前緣，10（S1）：93-100.

石耀霖，朱守彪.2006.用GPS位移資料計算應變方法的討論[J].大地測量與地球動力學，26（1）：1-8.

王靜.2019.GPS應變率場計算方法研究進展[J].地震，39（2）：124-136.

王敏，李強，王凡，等.2011.全球定位系統(tǒng)測定的2011年日本宮城MW9.0級地震遠場同震位移[J].科學通報，56（20）：13-16.

王敏，沈正康，甘衛(wèi)軍，等.2008.GPS連續(xù)監(jiān)測鮮水河斷裂形變場動態(tài)演化[J].中國科學：地球科學，38（5）：575-581.

王敏，沈正康，牛之俊，等.2003.現(xiàn)今中國大陸地殼運動與活動塊體模型[J].中國科學：地球科學，33（S1）：21-32，209.

王敏，沈正康.2020.中國大陸現(xiàn)今構(gòu)造變形：三十年的GPS觀測與研究[J].中國地震，36（4）：660-683.

王敏，張培震，沈正康，等.2006.全球定位系統(tǒng)（GPS）測定的印尼蘇門達臘巨震的遠場同震地表位移[J].科學通報，51（3）：365-368.

王琪，Seeber.1998.紅河斷裂的GPS監(jiān)測與現(xiàn)代構(gòu)造應力場[J].大地測量與地球動力學，18（2）：49-56.

王琪，牛之俊，石俊成.2003.中國地殼運動觀測網(wǎng)絡基本網(wǎng)觀測精度研究[J].大地測量與地球動力學，23（3）：11-15.

王琪，喬學軍，游新兆.2020.中國地震大地測量——半個世紀的歷程與科學貢獻[J].中國地震，36（4）：647-659.

王琪，張培震，馬宗晉.2002.中國大陸現(xiàn)今構(gòu)造變形GPS觀測數(shù)據(jù)與速度場[J].地學前緣，9（2）：415-429.

王琪，張培震，牛之俊，等.2001.中國大陸現(xiàn)今地殼運動和構(gòu)造變形[J].中國科學：地球科學，31（7）：529-536.

王琪，張培震.2006.用GPS研究中國大陸現(xiàn)今的地殼變形[J].科學觀察，1（6）：43.

王琪.2003.中國大陸現(xiàn)今地殼運動研究[J].地震學報，25（5）：541-547.

王小亞，朱文耀，符養(yǎng)，等.2002.GPS監(jiān)測的中國及其周邊現(xiàn)時地殼形變[J].地球物理學報，45（2）：198-209.

王小亞.2002.GPS在地球物理方面的應用[D].上海：中國科學院研究生院（上海天文臺）.

王閻昭，王恩寧，沈正康，等.2008.基于 GPS 資料約束反演川滇地區(qū)主要斷裂現(xiàn)今活動速率[J].中國科學：地球科學，38（5）：582-597.

王閱兵，甘衛(wèi)軍，陳為濤，等.2018.GNSS觀測的九寨溝7.0級地震同震位移初步結(jié)果[J].地球物理學報，61（1）：161-170.

魏明耀，張浪靜，郭強.2010.GPS技術(shù)在震后地質(zhì)災害監(jiān)測中的應用[J].測繪科學，35（S1）：192-194.

魏文薪，江在森，邵德勝，等.2018.2014年魯?shù)?.5級地震GPS同震位移及反演分析[J].地球物理學報，61（4）：1258-1265.

吳云，孫建中，喬學軍，等.2003.GPS在現(xiàn)今地殼運動與地震監(jiān)測中的初步應用[J].武漢大學學報（信息科學版），28（S1）：79-82.

伍吉倉，宋鑫友，胡鳳鳴，等.2020.聯(lián)合GNSS和InSAR觀測位移反演2008年汶川大地震斷層位錯模型參數(shù)[J].中國地震，36（4）：767-779.

武艷強，江在森，楊國華，等.2009.利用最小二乘配置在球面上整體解算 GPS 應變場的方法及應用[J].地球物理學報，52（7）：1707-1714.

武艷強，江在森，楊國華，等.2011.汶川地震前 GPS 資料反映的應變率場演化特征[J].大地測量與地球動力學，31（5）：20-25，29.

武艷強，江在森，朱爽，等.2020.中國大陸西部GNSS變形特征及其與M≥7.0強震孕育的關(guān)系[J].中國地震，36（4）：756-766.

徐韶光，熊永良，黃丁發(fā)，等.2011.長距離動態(tài)GPS數(shù)據(jù)處理在地震監(jiān)測中的應用研究[J].測繪科學，36（2）：88-90.

許家琪.2019.GNSS坐標時間序列噪聲模型建立影響因素研究[D].南昌：東華理工大學.

楊少敏，蘭啟貴，聶兆生，等.2012.用多種數(shù)據(jù)構(gòu)建2008年汶川特大地震同震位移場[J].地球物理學報，55（8）：2575-2588.

楊少敏，李杰，王琪.2008.GPS 研究天山現(xiàn)今變形與斷層活動[J].中國科學：地球科學，38（7）：872-880.

殷海濤，張培震，甘衛(wèi)軍，等.2010.高頻GPS測定的汶川MS8.0級地震震時近場地表變形過程[J].科學通報，55（26）：2621-2626.

尹繼堯，徐平，吳培稚，等.2007.GPS數(shù)據(jù)用于地震短臨預測的數(shù)據(jù)處理初步探索[J].地震地磁觀測與研究，28（5）：50-56.

游新兆，杜瑞林，王琪，等.2001.中國大陸地殼現(xiàn)今運動的GPS測量結(jié)果與初步分析[J].大地測量與地球動力學，21（3）：1-8.

游新兆，馬海建.2012.國家重大科技基礎設施項目——“中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡”建設成果展示[J].科技成果管理與研究，（6）：78-82.

游新兆，王啟梁.1994.青藏高原1993年GPS觀測成果的精度分析[J].地殼形變與地震，14（3）：27-33.

游新兆，楊少敏，王琪.2012.中國大陸地殼運動觀測研究進展[J].地震，32（2）：31-39.

張風霜，武艷強，孫東穎.2011.中國大陸GPS連續(xù)觀測站基線時間序列分析[J].地震，31（4）：86-97.

張培震，鄧起東，張國民，等.2003.中國大陸的強震活動與活動地塊[J].中國科學：地球科學，33（S1）：12-20.

張培震，王琪，馬宗晉.2002a.青藏高原現(xiàn)今構(gòu)造變形特征與GPS速度場[J].地學前緣，9（2）：442-450.

張培震，王琪，馬宗晉.2002b.中國大陸現(xiàn)今構(gòu)造運動的GPS速度場與活動地塊[J].地學前緣，9（2）：430-441.

張小紅，郭斐，郭博峰，等.2012.利用高頻 GPS 進行地表同震位移監(jiān)測及震相識別[J].地球物理學報，55（6）：1912-1918.

張小紅，吳明魁，劉萬科.2016.BeiDou B2/Galileo E5b短基線緊組合相對定位模型及性能評估[J].測繪學報，45（S2）：1-11.

張曉亮，江在森，王敏，等.2006.利用GPS連續(xù)站資料研究地殼運動與地震的關(guān)系[J].大地測量與地球動力學，26（4）：63-68.

張祖勝，耿世昌，陳德民，等.1989.現(xiàn)代地殼垂直形變速率[M]//馬杏垣.中國巖石圈動力學地圖集.北京：中國地圖出版社.

張祖勝.2001.中國地殼運動觀測網(wǎng)絡[J].中國基礎科學，（3）：45-49.

趙斌，杜瑞林，張銳，等.2015.GPS測定的尼泊爾MW7.9和MW7.3級地震同震形變場[J].科學通報，60（S2）：2758-2766.

趙斌，王敏，胡巖，等.2020.中國及鄰域強震震后變形監(jiān)測及巖石流變性質(zhì)研究[J].中國地震，36（4）：806-816.

趙靜，劉杰，任金衛(wèi)，等.2018.汶川8.0級地震前GPS與跨斷層資料反映的運動與變形演化特征[J].中國地震，34（3）：31-44.

趙齊樂，劉經(jīng)南，葛茂榮，等.2005.用PANDA對GPS和CHAMP衛(wèi)星精密定軌[J].大地測量與地球動力學，25（2）：113-116，122.

周碩愚，吳云，江在森，等.2017.地震大地測量學[M].武漢：武漢大學出版社.

周碩愚，吳云，施順英，等.2006.GPS對地震預報的推進和問題研究[J].大地測量與地球動力學，26（3）：111-117.

周碩愚，吳云，施順英，等.2007.現(xiàn)今地殼運動動力學基本狀態(tài)與地震可預報性研究[J].大地測量與地球動力學，27（4）：92-99.

周碩愚，吳云.2015.由震源到動力學系統(tǒng)──地震模式百年演化[J].大地測量與地球動力學，35（6）：911-918.

朱文耀，程宗頤，王小亞，等.1998.中國大陸地殼運動的背景場[J].天文學進展，16（4）：346-348.

Bedford J R，Moreno M，Deng Z，et al.2020.Months-long thousand-kilometre-scale wobbling before great subduction earthquakes[J].Nature，580（7805）：628-635.

Bilich A，Cassidy J F，Larson K M.2008.GPS Seismology：Application to the 2002 MW7.9 Denali fault earthquake[J].Bulletin of the Seismological Society of America，98（2）：593-606.

Blewitt G，Lavallée D.2002.Effect of annual signals on geodetic velocity[J].Journal of Geophysical Research：Solid Earth，107（B7）：ETG 9-1-ETG 9-11.

Bock Y，Prawirodirdjo L，Melbourne T I.2004.Detection of arbitrarily large dynamic ground motions with a dense high-rate GPS network[J].Geophysical Research Letters，31（6）：177-182.

Chen C H，Wen S，Liu J Y，et al.2014.Surface displacements in Japan before the 11 March 2011 M9.0 Tohoku-Oki earthquake[J].Journal of Asian Earth Sciences，80（2）：165-171.

Dach R，Hugentobler U，F(xiàn)ridez P.2007.Bernese GPS Software Version 5.0[M].Bern：Astronomical Institute，University of Bern.

Davis J L，Wernicke B P，Tamisiea M E.2012.On seasonal signals in geodetic time series[J].Journal of Geophysical Research：Solid Earth，117（B1）：B0143.

Davis J P，Smalley J R.2009.Love wave dispersion in central North America determined using absolute displacement seismograms from high-rate GPS[J].Journal of Geophysical Research：Solid Earth，114（B11）：292-310.

Emore G L，Haase J S，Choi K，et al.2007.Recovering seismic displacements through combined use of 1-Hz GPS and strong-motion accelerometers[J].Bulletin of the Seismological Society of America，97（2）：357-378.

Feng W，Ren J，Jiang Z.2016.GPS station short-term dynamic characteristics of micro displacement before Menyuan M6.4 earthquake[J].Geodesy and Geodynamics，7（4）：237-244.

Gan W，Zhang P，Shen Z K，et al.2007.Present-day crustal motion within the Tibetan Plateau inferred from GPS measurements[J].journal of geophysical research solid earth，112（B8）：B08416.

Geng J，Wen Q，Chen Q，et al.2019.Six-degree-of-freedom broadband seismogeodesy by combining collocated high-rate GNSS，accelerometers，and gyroscopes[J].Geophysical Research Letters，46（2）：708-716.

Geng J，Wen Q，Chen Q.2020.Six-degree-of-freedom seismogeodesy by combining collocated high-rate GNSS，accelerometers and gyroscopes[C]//EGU General Assembly Conference Abstracts.

Gong W，Zhang Y，Li T，et al.2019.Multi-sensor geodetic observations and modeling of the 2017 MW6.3 Jinghe earthquake[J].Remote Sensing，11（18）：2157.

Hao M，F(xiàn)reymueller J T，Wang Q，et al.2016.Vertical crustal movement around the southeastern Tibetan Plateau constrained by GPS and GRACE data[J].Earth and Planetary Science Letters，437（3）：1-8.

He P，Wang Q，Ding K，et al.2016.Source model of the 2015 MW6.4 Pishan earthquake constrained by interferometric synthetic aperture radar and GPS：Insight into blind rupture in the western Kunlun Shan[J].Geophysical research letters，43（4）：1511-1519.

Huang M H，Roland Bürgmann，F(xiàn)reed A M.2014.Probing the lithospheric rheology across the eastern margin of the Tibetan Plateau[J].Earth and Planetary Science Letters，396（6）：88-96.

Ji C，Larson K M，Tan Y，et al.2004.Slip history of the 2003 San Simeon earthquake constrained by combining 1-Hz GPS，strong motion，and teleseismic data[J].Geophysical research letters，31（17）：169-188.

Jiang Z，Wang M，Wang Y，et al.2014.GPS constrained coseismic source and slip distribution of the 2013 MW6.6 Lushan，China，earthquake and its tectonic implications[J].Geophysical Research Letters，41（2）：407-413.

Liu J N，Ge M R.2003.PANDA software and its preliminary result of positioning and orbit determination[J].Wuhan University Journal of Natural Sciences，8（2）：603.

Langbein J.2008.Noise in GPS displacement measurements from Southern California and Southern Nevada[J].Journal of Geophysical Research：Solid Earth，113（B5）：B05405.

Larson K M，Bodin P，Gomberg J.2003.Using 1-Hz GPS data to measure deformations caused by the Denali fault earthquake[J].Science，300（5624）：1421-1424.

Li J，Liu G，Qiao X，et al.2018.Rupture characteristics of the 25 November 2016 Aketao earthquake（MW6.6）in eastern Pamir revealed by GPS and teleseismic data[J].Pure Appl Geophys，175（2）：573-585.

Liang S，Gan W，Shen C，et al.2013.Three-dimensional velocity field of present-day crustal motion of the Tibetan Plateau derived from GPS measurements[J].Journal of Geophysical Research Solid Earth，118（10）：5722-5732.

Lin A，F(xiàn)u B，Guo J，et al.2002.Co-seismic strike-slip and rupture length produced by the 2001 MS8.1 Central Kunlun earthquake[J].Science，296（5575）：2015-2017.

Liu G，Xiong W，Wang Q，et al.2019.Source characteristics of the 2017 MS7.0 Jiuzhaigou，China，earthquake and implications for recent seismicity in eastern Tibet[J].Journal of Geophysical Research：Solid Earth，124（5）：4895-4915.

Meade B J.2007.Present-day kinematics at the India-Asia collision zone[J].Geology，35（1）：81-84.

Miyazaki S，Larson K M，Choi K，et al.2004.Modeling the rupture process of the 2003 September 25 Tokachi-Oki（Hokkaido）earthquake using 1-Hz GPS data[J].Geophysical Research Letters，31（21）：L21603.

Montillet J P，Yu K.2015.Modeling geodetic processes with levy-stable distribution and FARIMA[J].Mathematical Geosciences.，47（6）：627-646.

Meng G，Ren J，Su X，et al 2013.Coseismic deformation of the 2010 MW6.9 Yushu earthquake derived from GPS Data[J].Seismological Research Letter，84（1）：57-64.

Ozawa S，Nishimura T，Munekane H，et al.2012.Preceding，coseismic，and postseismic slips of the 2011 Tohoku earthquake，Japan[J].Journal of Geophysical Research：Solid Earth，117（B7）：B07404.

Ryder I，Bürgmann R，Sun J.2010.Tandem afterslip on connected fault planes following the 2008 Nima-Gaize（Tibet）earthquake[J].Journal of Geophysical Research：Solid Earth，115（B3）：B03404.

Shen Z K，Jackson D D，Ge B X.1996.Crustal deformation across and beyond the Los Angeles basin from geodetic measurements[J].Journal of Geophysical Research：Solid Earth，101（B12）：27957-27980.

Shen Z K，Lü J，Wang M，et al.2005.Contemporary crustal deformation around the southeast borderland of the Tibetan Plateau[J].Journal of Geophysical Research：Solid Earth，110（B11）：B11409.

Shi C，Lou Y，Zhang H，et al.2010.Seismic deformation of the MW8.0 Wenchuan earthquake from high-rate GPS observations[J].Advances in Space Research，46（2）：228-235.

Thatcher W.2007.Microplate model for the present-day deformation of Tibet[J].Journal of Geophysical Research Solid Earth，112（B1）：534-535.

Wang M，Shen Z K.2020.Present-day crustal deformation of continental China derived from GPS and its tectonic implications[J].Journal of Geophysical Research：Solid Earth，125（2）：e2019JB018774.

Wang Q.2001a.Present-day crustal deformation in China constrained by global positioning system measurements[J].Science，294（5542）：574-577.

Wang Q.2011b.Rupture of deep faults in the 2008 Wenchuan earthquake and uplift of the Longmen Shan[J].Nature Geoscience，4（9）：634–640.

Wang Y，Wang M，Shen Z K.2017.Block-like versus distributed crustal deformation around the northeastern Tibetan plateau[J].Journal of Asian Earth Sciences，140（11）：31-47.

Xu C J，F(xiàn)an Q B，Wang Q，et al.2014.Postseismic deformation after 2008 Wenchuan earthquake[J].Survey Review，46（339）：432-436.

Xu R，Stamps D S.2019.A geodetic strain rate and tectonic velocity model for China[J].Geochemistry，Geophysics，Geosystems，20（3）：1280-1297.

Yokota Y，Koketsu K.2015.A very long-term transient event preceding the 2011 Tohoku earthquake[J].Nature communications，6（1）：1-5.

Zhan W，Li F，Hao W，et al.2017.Regional characteristics and influencing factors of seasonal vertical crustal motions in Yunnan，China[J].Geophysical Journal International，210（3）：1295-1304.

Zhang P Z，Shen Z K，Wang M，et al.2004.Continuous deformation of the Tibetan Plateau from global positioning system data[J].Geology，32（9）：809-812.

Zhao B，Bürgmann R，Wang D，et al.2017b.Dominant controls of downdip afterslip and viscous relaxation on the postseismic displacements following the MW7.9 Gorkha，Nepal，earthquake[J].Journal of Geophysical Research：Solid Earth，122（10）：8376-8401.

Zhao B，Zhang C，Wang D，et al.2017a.Contemporary kinematics of the Ordos block，North China and its adjacent rift systems constrained by dense GPS observations[J].Journal of Asian Earth Sciences，135（3）：257-267.

Zhao D，Qu C，Shan X，et al.2018.InSAR and GPS derived coseismic deformation and fault model of the 2017 MS7.0 Jiuzhaigou earthquake in the Northeast Bayanhar block[J].Tectonophysics，726（4）：86-99.

Zheng G，Wang H，Wright T J，et al.2017.Crustal deformation in the India-Eurasia collision zone from 25 years of GPS measurements[J].Journal of Geophysical Research：Solid Earth，122（11）：9290-9312.

Zumberge J F，Heflin M B，Jefferson D C，et al.1997.Precise point positioning for the efficient and robust analysis of GPS data from large networks[J].Journal of Geophysical Research，102（B3）：5005.

GPS in Earthquake Monitoring in China：Current Situation and Prospect

WANG Tan1，2，LI Yu2，ZHANG Rui2，SHI Hongbo2，WANG Yuebing2

（1.Capital Normal University，Beijing 100048，China）

（2.China Earthquake Networks Center，Beijing 100045，China）

Abstract

In this paper we firstly review the development of application of GPS technology to earthquake monitoring in China，and then describe mainly GPS application to the following aspects：data processing，time series，velocity field and strain-rate field of the tectonic deformation of the Chinese mainland，quantification of the movement of active blocks in the Chinese mainland，and seismological research.Finally，we discuss the problems and challenges we are currently facing in GPS in earthquake monitoring in China，and put forward prospects from four aspects：earthquake prediction research，basic research，GPS network monitoring layout，and GPS solution accuracy and timeliness.

Keywords：GPS;earthquake monitoring;CMONOC;current situation of GPS application;development prospect

收稿日期：2021-02-08.

基金項目：中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡運維資助.

第一作者簡介：王坦（1991-），工程師，主要從事地殼形變觀測技術(shù)研究和地震站網(wǎng)運管信息化工作.E-mail：twang@seis.ac.cn