潘正洋，何建坤，盧雙疆

（1．中國(guó)科學(xué)院青藏高原研究所 大陸碰撞與高原隆升重點(diǎn)實(shí)驗(yàn)室，北京 100101；2．中國(guó)科學(xué)院大學(xué)，北京 100049）

0 引 言

全球定位系統(tǒng)（GPS）利用24顆繞地衛(wèi)星進(jìn)行實(shí)時(shí)精確的測(cè)定地面點(diǎn)位，以地面GPS固定站點(diǎn)之間的位置及其基線(xiàn)長(zhǎng)度隨時(shí)間變化來(lái)測(cè)定局部地區(qū)以及全球地殼的相對(duì)運(yùn)動(dòng)。GPS觀(guān)測(cè)信息豐富，全球覆蓋面廣，具有高精度、全天候、全球性、實(shí)時(shí)性的特點(diǎn)，為得到高時(shí)空分辨率地殼形變的可靠結(jié)果提供了保證［1－2］。高精度GPS的相對(duì)定位已經(jīng)達(dá)到上千千米長(zhǎng)基線(xiàn)的觀(guān)測(cè)中誤差可達(dá)毫米級(jí)精度的水平，這樣的精度足以檢測(cè)出地殼運(yùn)動(dòng)與構(gòu)造變形的細(xì)微變化［3－7］。

近年來(lái)，國(guó)內(nèi)外學(xué)者先后在帕米爾地區(qū)展開(kāi)了一系列的GPS觀(guān)測(cè)及其研究。王琪等針對(duì)天山及塔里木地區(qū)，在1992～1999年間多次利用GPS測(cè)量直接觀(guān)測(cè)到跨天山西段（76°E）每年約20mm的地殼快速縮短率及天山東段（87°E）每年約4mm的匯聚變形率［8－9］。楊少敏等利用在1992～2006年間天山境內(nèi)外地區(qū)近400個(gè)測(cè)站的GPS觀(guān)測(cè)，獲取了天山現(xiàn)今地殼運(yùn)動(dòng)速度場(chǎng)［10－11］。牛之俊等利用天山境內(nèi)外地區(qū)在1992～2005年382個(gè)GPS測(cè)站的原始觀(guān)測(cè)資料計(jì)算了天山地區(qū)現(xiàn)今地殼運(yùn)動(dòng)速率，并得出天山的匯聚速率由西向東逐步減小的結(jié)論［12］。Reigber等于1994～1998年間沿著塔吉克盆地北邊和東邊、北帕米爾高原、塔里木盆地到費(fèi)爾干納盆地及其南天山到哈薩克地臺(tái)，布置了92個(gè)GPS觀(guān)測(cè)臺(tái)站（CATS），得到北天山匯聚速率為每年（12±3）mm，北北西向到東東南向縮短速率為每年（23±3）mm［13］。Banerjee等在喜馬拉雅西北地區(qū)主沖斷帶MFT布設(shè)了26個(gè)GPS臺(tái)站，測(cè)得喀喇昆侖西北部的滑動(dòng)速率為每年（11±4）mm，并且得出青藏高原南部區(qū)域呈東西向伸展，西北部喜馬拉雅區(qū)域相對(duì)于南加帕爾巴特向西運(yùn)動(dòng)［14］。Mohadjer等在30°N～44°N和60°E～76°E區(qū)域布置了一系列的GPS站點(diǎn)，得到阿萊—南天山的匯聚速率為每年（12±2）mm，恰曼—伽德茲—庫(kù)納斷裂帶滑動(dòng)速率為每年（18±1）mm，達(dá)瓦孜—卡拉庫(kù)爾斷裂帶為每年（11±2）mm，塔拉斯—費(fèi)爾干納斷裂帶匯聚速率小于每年2mm［15］。Zubovich等于1994～2003年間在中亞地區(qū)布置了144個(gè)連續(xù)GPS觀(guān)測(cè)點(diǎn)和400個(gè)左右流動(dòng)GPS觀(guān)測(cè)點(diǎn)，得到塔里木盆地相對(duì)于歐亞板塊的匯聚速率為每年（20±2）mm；GPS觀(guān)測(cè)數(shù)據(jù)結(jié)合地質(zhì)數(shù)據(jù)顯示，帕米爾北緣俯沖到阿萊谷地至少需要速率為每年10mm，最大達(dá)每年15mm，塔拉斯—費(fèi)爾干納斷裂帶現(xiàn)今滑動(dòng)的最大速率為每年2mm［16］。

這些GPS觀(guān)測(cè)點(diǎn)主要分布在帕米爾高原西南側(cè)及其東北側(cè)天山附近，在喀喇昆侖斷裂地區(qū)及其阿里盆地地區(qū)的GPS觀(guān)測(cè)點(diǎn)非常稀少。筆者通過(guò)對(duì)帕米爾地區(qū)GPS測(cè)點(diǎn)進(jìn)行合理設(shè)置及補(bǔ)充，分別在南北向及其東西向布置了2條剖面進(jìn)行觀(guān)測(cè)，認(rèn)識(shí)青藏高原的形成過(guò)程和大陸動(dòng)力學(xué)理論，探討巖石圈變形的行為方式。

1 研究區(qū)概況

新生代以來(lái)，印度板塊與歐亞大陸碰撞和持續(xù)匯聚作用在中亞地區(qū)造成了強(qiáng)烈的陸內(nèi)形變，引起帕米爾高原大幅度向北推移旋轉(zhuǎn)以及西昆侖山、天山等古老造山帶的復(fù)活，使帕米爾及其鄰區(qū)成為大陸內(nèi)部構(gòu)造運(yùn)動(dòng)十分強(qiáng)烈的地區(qū)［17－18］。

在整個(gè)青藏高原及其附近地區(qū)的GPS速率圖上，可以明顯看出帕米爾地區(qū)的GPS觀(guān)測(cè)相對(duì)不足（圖1）。現(xiàn)有這些GPS觀(guān)測(cè)資料難以約束帕米爾高原一些大型斷裂的運(yùn)動(dòng)學(xué)特征，為了更加緊密具體地得到帕米爾地區(qū)一些主要斷裂的滑動(dòng)速率，更好地約束整個(gè)青藏高原及其帕米爾活動(dòng)區(qū)域的應(yīng)變場(chǎng)［11，19］，筆者在前人工作基礎(chǔ)上，于2011年7月在該地區(qū)進(jìn)行了首次觀(guān)測(cè)，并取得了一些結(jié)果（圖2）。

2 GPS觀(guān)測(cè)剖面布設(shè)與首期觀(guān)測(cè)

本期GPS測(cè)點(diǎn)位于帕米爾高原區(qū)域（36°N～42°N、68°E～78°E），南邊剖面從塔吉克盆地東北部到塔里木盆地西側(cè)，西邊剖面從費(fèi)爾干納盆地到帕米爾構(gòu)造結(jié)內(nèi)弧。野外選取無(wú)遮擋的出露基巖布設(shè)GPS測(cè)點(diǎn)，采用Trimble NetRs接收機(jī)和Zephry天線(xiàn)的GPS觀(guān)測(cè)設(shè)備，設(shè)置采樣率為30s，每個(gè)點(diǎn)位進(jìn)行了不低于36h的連續(xù)觀(guān)測(cè)。

3 GPS數(shù)據(jù)處理

本研究采用GAMIT／GLOBK軟件（10．40版）對(duì)首期GPS觀(guān)測(cè)數(shù)據(jù)進(jìn)行處理。GAMIT／GLOBK軟件是一種綜合性GPS精密分析軟件，由美國(guó)麻省理工學(xué)院（MIT）、哈佛大學(xué)天體物理研究中心（CFA）和斯克里普斯海洋研究所（SIO）共同開(kāi)發(fā)和維護(hù)，可用來(lái)解算地面觀(guān)測(cè)站的相對(duì)位置和衛(wèi)星軌道及其他參數(shù)的數(shù)據(jù)處理。

圖1 青藏高原GPS速率分布Fig．1 Distribution of GPS Rates in Tibetan Plateau

圖2 帕米爾地區(qū)主要GPS測(cè)點(diǎn)分布Fig．2 Distribution of GPS Sites in Pamirs Plateau

首期GPS觀(guān)測(cè)數(shù)據(jù)處理主要分為2步進(jìn)行：①利用GAMIT軟件進(jìn)行單天解解算，估計(jì)站位置、衛(wèi)星軌道參數(shù)等；②利用GLOBK軟件獲得時(shí)間序列檢驗(yàn)坐標(biāo)的一致性?；€(xiàn)處理時(shí)，選擇周?chē)腎GS（International GNSS Service）站與自己剖面測(cè)量的數(shù)據(jù)一起解算。解算時(shí)，主要考慮如下因素及其相關(guān)設(shè)置［20－24］：①電離層折射影響用LC觀(guān)測(cè)值消除；②使用8個(gè)IGS站（BJFS，BRAO，HYDE，IRKT，POL2，KUNM，SELE，TASH）數(shù)據(jù)作為參考站約束；③使用來(lái)自IGS的精密衛(wèi)星軌道和地球定向參數(shù)；④同時(shí)估計(jì)衛(wèi)星軌道誤差和測(cè)站坐標(biāo)，即松弛IGS軌道；⑤截止高度角為15°，歷元間隔為30 s；⑥對(duì)流層折射根據(jù)標(biāo)準(zhǔn)大氣模型用薩斯坦莫寧（Saastamoinen）模型改正；⑦衛(wèi)星鐘差的模型改正（用廣播星歷中的鐘差參數(shù)）；⑧接收機(jī)鐘差的模型改正（用根據(jù)偽距觀(guān)測(cè)值計(jì)算出的鐘差）；⑨衛(wèi)星和接收機(jī)天線(xiàn)相位中心改正，衛(wèi)星與接收機(jī)天線(xiàn)L1、L2相位中心偏差采用GAMIT軟件的設(shè)定值；⑩測(cè)站位置的潮汐改正以及光壓模型BERNE（與SIO的全球解一致）。IGS跟蹤站（NER測(cè)站）坐標(biāo)約束為0．005、0．005、0．010m；GPS監(jiān)測(cè)站（NER測(cè)站）坐標(biāo)約束為9．999、9．999、9．999m。

在完成基線(xiàn)計(jì)算之后，進(jìn)行GLOBK平差計(jì)算。計(jì)算中，采用ITRF08框架作為網(wǎng)平差的基礎(chǔ)，主要約束為：①I(mǎi)GS核心站約束，共采用44個(gè)測(cè)站，每站坐標(biāo)的約束為0．003、0．003、0．010m；②軌道約束，與SIO的全球解H文件聯(lián)合求解，放松軌道至20m。

利用GAMIT解算時(shí)，首先需要建立工程目錄（project），然后準(zhǔn)備 GPS數(shù)據(jù)文件（rinex文件），導(dǎo)航文件（n文件）和星歷文件（sp3文件），接下來(lái)更新表文件（tables），并建立表文件目錄與該工程的連接，最后編輯控制文件包括處理控制文件（process．defaults）、站處理控制文件（site．defaults）、時(shí)段控制表（sestbl．）、站控制表（sittbl．）、站信息文件（station．info）和站坐標(biāo)文件（lfile．）及其準(zhǔn)備參考框架文件（itrf08．a(chǎn)pr），準(zhǔn)備工作完成之后，即可運(yùn)行GAMIT的批處理命令進(jìn)行解算。在解算過(guò)程中，GAMIT通過(guò)makexp和makex準(zhǔn)備數(shù)據(jù)，用fixdrv模塊生成批處理文件，用arc模塊進(jìn)行衛(wèi)星軌道積分，用model模塊計(jì)算觀(guān)測(cè)值殘差和偏導(dǎo)數(shù)，用autcln模塊自動(dòng)編輯數(shù)據(jù)，最后用solve模塊進(jìn)行最小二乘分析計(jì)算參數(shù)。

利用GLOBK解算時(shí)，首先須將GAMIT的松弛解H文件由ASCII格式轉(zhuǎn)換成二進(jìn)制格式，然后準(zhǔn)備好先驗(yàn)坐標(biāo)文件、地球極移文件、計(jì)算控制文件等，利用glred和glorg以批處理形式進(jìn)行解算。

4 結(jié)果分析

GAMIT計(jì)算得到的單天解標(biāo)準(zhǔn)化均方差（normalized rms，NRMS）是衡量單天解質(zhì)量的重要指標(biāo)之一。據(jù)國(guó)內(nèi)外GPS數(shù)據(jù)處理經(jīng)驗(yàn)，NRMS值一般應(yīng)小于0．3mm，若標(biāo)準(zhǔn)化均方差太大，則說(shuō)明處理過(guò)程中周跳可能沒(méi)有完全得到修正。本研究數(shù)據(jù)時(shí)間范圍為2011年年積日為第199天至第228天，檢查q文件后，處理30d數(shù)據(jù)的標(biāo)準(zhǔn)化均方差基本符合解算要求［25－28］（表1）。

表1 單天解標(biāo)準(zhǔn)化均方根誤差值Tab．1 NRMS of Day Solutions

在GAMIT得到單天解后，可以得到衛(wèi)星與測(cè)站的單程均方根殘差，該指數(shù)可以在單天解總結(jié)文件（sh＿gamit＿ddd．summary）中查看，并且在用sh＿make＿sky＿gifs或者批處理命令進(jìn)行批處理生成每個(gè)測(cè)站的天空?qǐng)D時(shí)，對(duì)均方根殘差有一個(gè)直觀(guān)的可視化呈現(xiàn)（圖3）。其中，0°對(duì)應(yīng)子午線(xiàn)。

另外，在評(píng)價(jià)GAMIT解算過(guò)程中，可以從GAMIT預(yù)處理結(jié)算中LC相位殘差得到站點(diǎn)的矯正模型，并且得到一些高度角與殘差的關(guān)系，以KHUF測(cè)站為例，均方殘差為4．88mm，基本符合測(cè)站的要求（圖4、5）。

圖3 KHUF測(cè)站點(diǎn)在第210天的天空?qǐng)DFig．3 Sky Plots of Day 210at Site KHUF

圖4 KHUF測(cè)站的高度角與均方根殘差的關(guān)系Fig．4 Relationship Between Elevation Angle and RootMean Square Residual at Site KHUF

在用GLOBK結(jié)算時(shí)，也可以對(duì)站點(diǎn)的信息進(jìn)行統(tǒng)計(jì)分析，檢驗(yàn)GAMIT處理結(jié)果，用gdl文件及其GMT進(jìn)行畫(huà)圖處理。圖6左邊縱軸表示GPS站點(diǎn)名稱(chēng)，右邊縱軸表示站點(diǎn)個(gè)數(shù)，方形代表單天解的個(gè)數(shù)（即每個(gè)站點(diǎn)的觀(guān)測(cè)時(shí)段，一天對(duì)應(yīng)一個(gè)觀(guān)測(cè)文件）。

利用GLOBK軟件對(duì)觀(guān)測(cè)期間各測(cè)點(diǎn)的點(diǎn)位重復(fù)性進(jìn)行計(jì)算，并得到每個(gè)測(cè)點(diǎn)的時(shí)間序列，結(jié)果顯示測(cè)點(diǎn)的點(diǎn)位重復(fù)性北向分量精度在3mm以?xún)?nèi)，東向分量在3mm以?xún)?nèi)，垂向分量在6mm以?xún)?nèi)。選取觀(guān)測(cè)時(shí)間最長(zhǎng)的DASH站點(diǎn)為代表，從中可以看出，該站點(diǎn)水平方向北向分量精度為1．04mm，東向分量為1．01mm，垂向分量為4．41mm。這體現(xiàn)了該站點(diǎn)在水平方向具有很好的重復(fù)精度（圖7）。

圖5 KHUF測(cè)點(diǎn)模型Fig．5 Model of Site KHUF

圖6 站點(diǎn)文件信息統(tǒng)計(jì)分布Fig．6 Distribution of File Information of Sites

圖7 測(cè)點(diǎn)DASH的時(shí)間序列分布Fig．7 Distribution of Temporal Series of Site DASH

多天解利用GLOBK軟件對(duì)觀(guān)測(cè)期間所有測(cè)點(diǎn)進(jìn)行聯(lián)合解算，并將結(jié)果歸算到ITRF08參考框架下［29－34］。多天整體解利用 GAMIT 軟件30d（年積日為第199天至第228天）的處理結(jié)果H文件，結(jié)合SIO的全球解H文件，并且將30d的單天解進(jìn)行合并，生成合并解GLX文件，利用全球stab＿site文件及其core＿site文件聯(lián)合itrf08．a(chǎn)pr文件，將其結(jié)果歸算到ITRF08參考框架中，再用GLOBK命令作綜合處理，通過(guò)平差計(jì)算得到最終點(diǎn)位和基線(xiàn)結(jié)果，結(jié)果歸算到2011年8月16日，點(diǎn)位精度水平分量為2～3mm，高程分量為3～7mm（表2）。

表2 測(cè)站坐標(biāo)及其精度Tab．2 Coordinates and Precisions of Sites

5 結(jié) 語(yǔ)

2011年7月至8月，對(duì)帕米爾高原及其鄰近地區(qū)進(jìn)行了GPS觀(guān)測(cè)剖面的布設(shè)和首期GPS觀(guān)測(cè)，這是2條真正意義上的GPS觀(guān)測(cè)剖面。筆者利用GAMIT／GLOBK軟件對(duì)首期觀(guān)測(cè)數(shù)據(jù)進(jìn)行處理，引入周邊8個(gè)IGS站點(diǎn)參與聯(lián)合解算，并將結(jié)果歸算到ITRF08下，解算得到的準(zhǔn)觀(guān)測(cè)值均滿(mǎn)足所須精度要求，其中所有站點(diǎn)的標(biāo)準(zhǔn)化均方差均小于0．3mm。在觀(guān)測(cè)區(qū)域單天解的基礎(chǔ)上，利用GLOBK合并全球解，并且進(jìn)行全區(qū)域的綜合平差，得到GPS站點(diǎn)的精確坐標(biāo)，結(jié)果顯示所有站點(diǎn)的點(diǎn)位水平分量為1～2mm，高程分量為3～7mm。

有必要在未來(lái)繼續(xù)進(jìn)行2～3期GPS觀(guān)測(cè)，完善現(xiàn)有的GPS觀(guān)測(cè)資料，幫助認(rèn)識(shí)帕米爾高原的運(yùn)動(dòng)學(xué)特征和青藏高原隆升機(jī)制。

QOCA（Quasi－observation Combination Analysis）學(xué)習(xí)班的董大南老師、中國(guó)地震局地震預(yù)測(cè)研究所武艷強(qiáng)副研究員和中國(guó)地震局地殼應(yīng)力研究所田云鋒副研究員在GAMIT／GLOBK軟件處理中提供了指導(dǎo)，圖件都由GMT完成，在此一并致謝。

［1］ 黨亞民．GPS和地球動(dòng)力學(xué)進(jìn)展［J］．測(cè)繪科學(xué)，2004，29（2）：77－79．DANG Ya－min．Progress in GPS and Its Geodynamics［J］．Science of Surveying and Mapping，2004，29（2）：77－79．

［2］ 宋淑麗，朱文耀，廖新浩．GPS應(yīng)用于地球動(dòng)力學(xué)研究的進(jìn)展［J］．天文學(xué)進(jìn)展，2003，21（2）：95－112．SONG Shu－li，ZHU Wen－yao，LIAO Xin－h(huán)ao．Recent Progress of GPS in Geodynamics［J］．Progress in Astronomy，2003，21（2）：95－112．

［3］ 王 敏，沈正康，董大南．非構(gòu)造形變對(duì)GPS連續(xù)站位置時(shí)間序列的影響和修正［J］．地球物理學(xué)報(bào)，2005，48（5）：1045－1052．WANG Min，SHEN Zheng－kang，DONG Da－nan．Effects of Non－tectonic Crustal Deformation on Continuous GPS Position Time Series and Correction to Them［J］．Chinese Journal of Geophsysics，2005，48（5）：1045－1052．

［4］ 王 敏，張祖勝，許明元，等．2000國(guó)家GPS大地控制網(wǎng)的數(shù)據(jù)處理和精度評(píng)估［J］．地球物理學(xué)報(bào)，2005，48（4）：817－823．WANG Min，ZHANG Zu－sheng，XU Ming－yuan，et al．Data Processing and Accuracy Analysis of National 2000GPS Geodetic Control Network［J］．Chinese Journal of Geophsysics，2005，48（4）：817－823．

［5］ 張培震，王 琪，馬宗晉．中國(guó)大陸現(xiàn)今構(gòu)造運(yùn)動(dòng)的GPS速度場(chǎng)與活動(dòng)地塊［J］．地學(xué)前緣，2002，9（2）：430－441．ZHANG Pei－zhen，WANG Qi，MA Zong－jin．GPS Velocity Field and Active Crustal Blocks of Contemporary Tectonic Deformation in Continental China［J］．Earth Science Frontiers，2002，9（2）：430－441．

［6］ 張培震，王 琪，馬宗晉．青藏高原現(xiàn)今構(gòu)造變形特征與 GPS速度場(chǎng)［J］．地學(xué)前緣，2002，9（2）：442－450．ZHANG Pei－zhen，WANG Qi，MA Zong－jin．GPS Velocity Field and Active Crustal Deformation in and Around the Qinghai－Tibet Plateau［J］．Earth Science Frontiers，2002，9（2）：442－450．

［7］ 張永志，胡 斌，陳宏賓．應(yīng)用GPS觀(guān)測(cè)青藏高原東北緣應(yīng)力場(chǎng)變化［J］．長(zhǎng)安大學(xué)學(xué)報(bào)：地球科學(xué)版，2003，25（4）：80－84．ZHANG Yong－zhi，HU Bin，CHEN Hong－bin．Variation of Stress Computed from GPS Measurement in Northeastern Margin of Qinghai－Tibet Plateau，China［J］．Journal of Chang'an University：Earth Science Edition，2003，25（4）：80－84．

［8］ 王 琪，張培震，馬宗晉．中國(guó)大陸現(xiàn)今構(gòu)造變形GPS觀(guān)測(cè)數(shù)據(jù)與速度場(chǎng)［J］．地學(xué)前緣，2002，9（2）：415－429．WANG Qi，ZHANG Pei－zhen，MA Zong－jin．GPS Database and Velocity Field of Contemporary Tectonic Deformation in Continental China［J］．Earth Science Frontiers，2002，9（2）：415－429．

［9］ 王 琪，丁國(guó)瑜，喬學(xué)軍，等．天山現(xiàn)今地殼快速縮短與南北地塊的相對(duì)運(yùn)動(dòng)［J］．科學(xué)通報(bào)，2000，45（14）：1543－1547．WANG Qi，DING Guo－yu，QIAO Xue－jun，et al．Rapid Crustal Shortening and Relative Motion of Crustal Blocks in the Tianshan［J］．Chinese Science Bulletin，2000，45（14）：1543－1547．

［10］ 楊少敏，李 杰，王 琪．GPS研究天山現(xiàn)今變形與斷層活動(dòng)［J］．中國(guó)科學(xué)：D輯，2008，38（7）：872－880．YANG Shao－min，LI Jie，WANG Qi．The Deformation Pattern and Fault Rate in the Tianshan Mountains Inferred from GPS Observations［J］．Science in China：Series D，2008，38（7）：872－880．

［11］ 楊少敏，王 琪，游新兆．中國(guó)現(xiàn)今地殼運(yùn)動(dòng)GPS速度場(chǎng)的連續(xù)變形分析［J］．地震學(xué)報(bào)，2005，27（2）：128－138．YANG Shao－min，WANG Qi，YOU Xin－zhao．Numerical Analysis of Contemporary Horizontal Tectonic Deformation Fields in China from GPS Data［J］．Acta Seismologica Sinica，2005，27（2）：128－138．

［12］ 牛之俊，游新兆，楊少敏．利用GPS分析天山現(xiàn)今地殼形變特征［J］．大地測(cè)量與地球動(dòng)力學(xué)，2007，27（2）：1－9．NIU Zhi－jun，YOU Xin－zhao，YANG Shao－min．Analysis of Contemporary Crustal Deformation Characteristics with GPS Data of Tianshan Mountain［J］．Journal of Geodesy and Geodynamics，2007，27（2）：1－9．

［13］ REIGBER C，MICHEL G W，GALAS R，et al．New Space Geodetic Constraints on the Distribution of Deformation in Central Asia［J］．Earth and Planetary Science Letters，2001，191（1／2）：157－165．

［14］ BANERJEE P，BURGMANN R．Convergence Across the Northwest Himalaya from GPS Measurements［J］．Geophysical Research Letters，2002，29（13）：1652．

［15］ MOHADJER S，BENDICK R，ISCHUK A，et al．Partitioning of India－Eurasia Convergence in the Pamir－Hindu Kush from GPS Measurements［J］．Geophysical Research Letters，2010，37（4）：L04305．

［16］ ZUBOVICH A V，WANG X Q，SCHERBA Y G，et al．GPS Velocity Field for the Tian Shan and Surrounding Regions［J］．Tectonics，2010，29（6）：TC6014．

［17］ 陳 杰，李 濤，李文巧，等．帕米爾構(gòu)造結(jié)及鄰區(qū)的晚新生代構(gòu)造與現(xiàn)今變形［J］．地震地質(zhì)，2011，33（2）：241－259．CHEN Jie，LI Tao，LI Wen－qiao，et al．Late Cenozoic and Present Tectonic Deformation in the Pamir Salient，Northwestern China［J］．Seismology and Geology，2011，33（2）：241－259．

［18］ FU B H，NINOMIYA Y，GUO J M．Slip Partitioning in the Northeast Pamir－Tian Shan Convergence Zone［J］．Tectonophysics，2010，483（3）：344－364．

［19］ 孫文斌，和躍時(shí)，常 征，等．帕米爾—興都庫(kù)什地區(qū)板塊俯沖及其應(yīng)力狀態(tài)［J］．地震地質(zhì)，2009，31（2）：207－217．SUN Wen－bin，HE Yue－shi，CHANG Zheng，et al．The Plate Subduction and Stress State in the Pamir－Hindu Kush Region［J］．Seismology and Geology，2009，31（2）：207－217．

［20］ 王 敏，張培震，沈正康，等．全球定位系統(tǒng)（GPS）測(cè)定的印尼蘇門(mén)達(dá)臘巨震的遠(yuǎn)場(chǎng)同震地表位移［J］．科學(xué)通報(bào)，2006，51（3）：365－368．WANG Min，ZHANG Pei－zhen，SHEN Zheng－kang，et al．Far－field Coseismic Displacements Associated with the Great Sumatra Earthquakes of December 26，2004and March 29，2005Constrained by Global Positioning System［J］．Chinese Science Bulletin，2006，51（3）：365－368．

［21］ 王 敏，沈正康，甘衛(wèi)軍．GPS連續(xù)監(jiān)測(cè)鮮水河斷裂形變場(chǎng)動(dòng)態(tài)演化［J］．中國(guó)科學(xué)：D 輯，2008，38（5）：575－581．WANG Min，SHEN Zheng－kang，GAN Wei－jun．GPS Monitoring of Temporal Deformation of the Xianshuihe Fault［J］．Science in China：Sereis D，2008，38（5）：575－581．

［22］ 王閻昭，王恩寧，沈正康，等．基于GPS資料約束反演川滇地區(qū)主要斷裂現(xiàn)今活動(dòng)速率［J］．中國(guó)科學(xué)：D輯，2008，38（5）：582－597．WANG Yan－zhao，WANG En－ning，SHEN Zheng－kang，et al．GPS－constrained Inversion of Present－day Slip Rates Along Major Faults of the Sichuan－Yunnan Region［J］．Science in China：Series D，2008，38（5）：582－597．

［23］ 張培震，王 敏，甘衛(wèi)軍，等．GPS觀(guān)測(cè)的活動(dòng)斷裂滑動(dòng)速率及其對(duì)現(xiàn)今大陸動(dòng)力作用的制約［J］．地學(xué)前緣，2003，10（增1）：81－92．ZHANG Pei－zhen，WANG Min，GAN Wei－jun，et al．Slip Rates Along Major Active Faults from GPS Measurements and Constraints on Contemporary Continental Tectonics［J］．Earth Science Frontiers，2003，10（S1）：81－92．

［24］ 張曉亮，張 希，王雙緒，等．GPS資料揭示的汶川8．0級(jí)地震前地殼動(dòng)態(tài)演化特征［J］．大地測(cè)量與地球動(dòng)力學(xué)，2011，31（4）：1－5．ZHANG Xiao－liang，ZHANG Xi，WANG Shuang－xu，et al．Characteristics of Crustal Dynamic Evolution Revealed by GPS Observation Before Wenchuan 8．0 Earthquake［J］．Journal of Geodesy and Geodynamics，2011，31（4）：1－5．

［25］ 孟國(guó)杰，趙承坤，顧國(guó)華，等．中國(guó)地殼運(yùn)動(dòng)觀(guān)測(cè)網(wǎng)絡(luò)首期觀(guān)測(cè)數(shù)據(jù)處理［J］．地震，2001，21（1）：63－68．MENG Guo－jie，ZHAO Cheng－kun，GU Guo－h(huán)ua，et al．Data Procession of the First Epoch Measurement of the Crustal Movement Observation Network in China［J］．Earthquake，2001，21（1）：63－68．

［26］ 牛之俊，王 敏，孫漢榮，等．中國(guó)大陸現(xiàn)今地殼運(yùn)動(dòng)速度場(chǎng)的最新觀(guān)測(cè)結(jié)果［J］．科學(xué)通報(bào)，2005，50（8）：839－840．NIU Zhi－jun，WANG Min，SUN Han－rong，et al．Contemporary Velocity Field of Crustal Movement of Chinese Mainland from Global Positioning System Measurements［J］．Chinese Science Bulletin，2005，50（8）：839－840．

［27］ 田云鋒，沈正康．GPS坐標(biāo)時(shí)間序列中非構(gòu)造噪聲的剔除方法研究進(jìn)展［J］．地震學(xué)報(bào)，2009，31（1）：68－81．TIAN Yun－feng，SHEN Zheng－kang．Progress on Reduction of Non－tectonic Noise in GPS Position Time Series［J］．Acta Seismologica Sinica，2009，31（1）：68－81．

［28］ 田云鋒，沈正康，李 鵬．連續(xù)GPS觀(guān)測(cè)中的相關(guān)噪聲分析［J］．地震學(xué)報(bào)，2010，32（6）：696－704．TIAN Yun－feng，SHEN Zheng－kang，LI Peng．Analysis on Correlated Noise in Continuous GPS Observations［J］．Acta Seismologica Sinica，2010，32（6）：696－704．

［29］ DONG D，HERRING T A，KING R W．Estimating Regional Deformation from a Combination of Space and Terrestrial Geodetic Data［J］．Journal of Geodesy，1998，72（4）：200－214．

［30］ GAN W J，ZHANG P Z，SHEN Z K，et al．Presentday Crustal Motion Within the Tibetan Plateau Inferred from GPS Measurements［J］．Journal of Geophysical Research，2007，112（B8）：416．

［31］ KING R W，SHEN F，BURCHFIEL B C，et al．Geodetic Measurement of Crustal Motion in Southwest China［J］．Geology，1997，25（2）：179－182．

［32］ SHEN Z K，KING R W，AGNEW D C，et al．A Unified Analysis of Crustal Motion in Southern California，1970－2004：The SCEC Crustal Motion Map［J］．Journal of Geophysical Research，2011，116（B11）：402．

［33］ SHEN Z K，ZHAO C K，YIN A，et al．Contemporary Crustal Deformation in East Asia Constrained by Global Positioning System Measurements［J］．Journal of Geophysical Research，2000，105（B3）：5721－5734．

［34］ SHEN Z K，LU J，WANG M，et al．Contemporary Crustal Deformation Around the Southeast Borderland of the Tibetan Plateau［J］．Journal of Geophysical Research，2005，110（B11）：409．