周 琳 王慶良 李長(zhǎng)軍 郝 明 宋尚武

1 中國(guó)地震局第二監(jiān)測(cè)中心，西安市西影路316號(hào)，710054 2 防災(zāi)科技學(xué)院學(xué)科與研究生處，北京市東燕郊學(xué)院街465號(hào)，101601

?

基于GPS和水準(zhǔn)資料的拉脊山斷裂帶西段地殼形變研究

周 琳1王慶良1李長(zhǎng)軍1郝 明1宋尚武2

1 中國(guó)地震局第二監(jiān)測(cè)中心，西安市西影路316號(hào)，710054 2 防災(zāi)科技學(xué)院學(xué)科與研究生處，北京市東燕郊學(xué)院街465號(hào)，101601

利用拉脊山地區(qū)1983～1995年水準(zhǔn)資料，結(jié)合1999～2007年、2009～2014年GPS數(shù)據(jù)對(duì)拉脊山斷裂帶西段進(jìn)行三維地殼形變分析。拉脊山斷裂帶西段垂直運(yùn)動(dòng)速率約為1±0.5 mm/a，表明該段處于隆升狀態(tài)；與GPS數(shù)據(jù)顯示出的水平擠壓一致，拉脊山地區(qū)仍處于地殼水平擠壓縮短狀態(tài)，2008年之后斷層兩側(cè)差異運(yùn)動(dòng)更加明顯，地殼縮短速率增加并表現(xiàn)出明顯的左旋走滑分量，與衛(wèi)星影像解譯和野外地質(zhì)地貌調(diào)查結(jié)果一致。拉脊山斷裂西段受擠壓環(huán)境控制下的構(gòu)造變形現(xiàn)今仍在持續(xù)進(jìn)行，同時(shí)受全新世活動(dòng)的日月山斷裂右旋擠壓應(yīng)力影響。該地段是地震孕育的敏感部位，應(yīng)該引起關(guān)注。

拉脊山斷裂帶；GPS；水準(zhǔn)測(cè)量

拉脊山西起日月山，東至民和龔?fù)?，全長(zhǎng)200 多km，寬10～20 km，呈NNW-SSE向并向NE凸出[1]，是青藏高原腹地與黃土高原之間地貌梯級(jí)帶的分界線。拉脊山斷裂以晚更新世活動(dòng)為主，僅局部段活動(dòng)時(shí)代可能為全新世早期。拉脊山南北兩側(cè)有歷史記載的 5 級(jí)左右中等破壞性地震 20余次[2]。地質(zhì)考察未發(fā)現(xiàn)該地區(qū)存在大范圍發(fā)震的區(qū)域性斷裂，這些地震可能表明拉脊山斷裂帶的全新世活動(dòng)已經(jīng)擴(kuò)展到兩側(cè)的盆地，公元138年的金城-隴西6.75級(jí)地震可能與拉脊山斷裂活動(dòng)在盆地內(nèi)的擴(kuò)展相關(guān)[3]。拉脊山南北兩條邊界斷裂作用顯著，巖石地球化學(xué)非常復(fù)雜，剝露歷史存在爭(zhēng)議，一直是中外地質(zhì)學(xué)家的研究熱點(diǎn)。本文在前人研究基礎(chǔ)上，通過(guò)水準(zhǔn)資料和GPS數(shù)據(jù)資料，進(jìn)一步綜合探討該地區(qū)的三維地殼形變特征。

1 資料概況

本文使用覆蓋該區(qū)內(nèi)主要活動(dòng)斷裂的陜甘寧青地震水準(zhǔn)監(jiān)測(cè)網(wǎng)[4]西寧環(huán)1983～1995年水準(zhǔn)數(shù)據(jù)(圖1中實(shí)心圓所示)，測(cè)線橫跨拉脊山斷裂帶西段，拉脊山與日月山相接地區(qū)東側(cè)，總長(zhǎng)近100 km。本文在“中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)”1999、2001、2004、2007、2009、2011和2013年7期GPS區(qū)域站高精度觀測(cè)資料的基礎(chǔ)上[5-6]，結(jié)合中國(guó)地震行業(yè)重大科研專項(xiàng)“中國(guó)綜合地球物理場(chǎng)觀測(cè)”2010～2011年、2012～2014年GPS新測(cè)資料進(jìn)行研究。

實(shí)心三角形代表GPS點(diǎn)，實(shí)心圓代表水準(zhǔn)點(diǎn)，左下角插圖為研究區(qū)域的地質(zhì)構(gòu)造背景圖1 研究區(qū)域數(shù)據(jù)點(diǎn)位分布圖Fig.1 Distribution of data points in study area

2 數(shù)據(jù)處理

選取西寧附近(西寧環(huán)55)作為基點(diǎn)，利用1983和1995年兩期水準(zhǔn)數(shù)據(jù)，計(jì)算每期數(shù)據(jù)的相對(duì)高差值，并與兩期數(shù)據(jù)時(shí)間差值相除得到研究區(qū)域相對(duì)于基點(diǎn)的垂直形變速率，見圖2。

圖2 拉脊山斷裂帶西段垂直形變速率圖Fig.2 Vertical deformation velocity in the west section of Lajishan fault

GPS觀測(cè)數(shù)據(jù)采用GAMIT/GLOBK軟件進(jìn)行精化處理。首先利用GAMIT軟件計(jì)算單日松弛解，解算時(shí)采用SOPAC數(shù)據(jù)中心提供的精密星歷，加入中國(guó)大陸及周邊IGS測(cè)站的觀測(cè)數(shù)據(jù)進(jìn)行處理；之后通過(guò) GLOBK軟件計(jì)算多時(shí)段綜合解，結(jié)合單日松弛解和全球IGS站單日松弛解，選取全球均勻分布的90多個(gè)穩(wěn)定站為基準(zhǔn)點(diǎn)，將區(qū)域站速度場(chǎng)固定到ITRF2008框架下；最后利用ITRF2008框架下歐亞板塊歐拉極，獲取研究區(qū)域相對(duì)于穩(wěn)定歐亞板塊的水平運(yùn)動(dòng)速度場(chǎng)圖像，結(jié)果見圖3。

圖3 拉脊山地區(qū)GPS水平速度場(chǎng)(相對(duì)于歐亞板塊)Fig.3 Horizontal velocity field for the Lajishan fault area with respect to Eurasia block

3 形變分析

受青藏高原北東向擴(kuò)展的影響，拉脊山斷裂帶自西向東其走向由N60°W轉(zhuǎn)為近EW向，性質(zhì)以擠壓逆沖為主，體現(xiàn)為北東方向突出的弧形斷裂，是NWW向的西秦嶺北緣左旋走滑斷裂帶和NNW向的熱水-日月山右旋走滑斷裂帶之間的擠壓構(gòu)造區(qū)和構(gòu)造轉(zhuǎn)換帶，成為北側(cè)西寧-民和盆地和南側(cè)循化-化隆盆地之間重要的邊界斷裂[6]。拉脊山及鄰區(qū)磁組構(gòu)特征研究表明，在第三紀(jì)后期到第四紀(jì)曾發(fā)生過(guò)北東至南西向的構(gòu)造應(yīng)力擠壓，為拉脊山地區(qū)的持續(xù)抬升提供了重要依據(jù)[7]。本文水準(zhǔn)數(shù)據(jù)結(jié)果見圖4，測(cè)線橫跨拉脊山地區(qū)西段，拉脊山與日月山相接地區(qū)東側(cè)。結(jié)果表明，拉脊山地區(qū)西段仍舊處于隆升狀態(tài)，隆升速率為1±0.5 mm/a，隆升速率由拉脊山南緣斷裂到北緣斷裂不斷變大，靠近西寧盆地隆升速率變小。拉脊山地區(qū)地形起伏大，與北側(cè)西寧盆地形成明顯高差[8]。從圖4可以看出，拉脊山兩側(cè)5 km附近相對(duì)隆升速率變小。野外調(diào)查顯示，該地區(qū)地貌形態(tài)受地表流水侵蝕，山體溝谷發(fā)育，是拉脊山北緣斷裂帶和拉脊山南緣斷裂帶仍舊活動(dòng)的有力證據(jù)。

拉脊山地區(qū)弱震活動(dòng)基本沿拉脊山南、北緣弧形斷裂呈帶狀分布，歷史上發(fā)生過(guò)20余次5級(jí)左右的中等破壞性地震。本文研究區(qū)域在拉脊山南緣斷裂帶和北緣斷裂帶的西段，衛(wèi)星影像解譯和野外地質(zhì)地貌調(diào)查發(fā)現(xiàn)，拉脊山斷裂西段左旋走滑活動(dòng)跡象明顯，兼具逆沖分量和左旋走滑分量，以擠壓逆沖為主[6]。受全新世活動(dòng)的日月山斷裂右旋擠壓應(yīng)力的影響，拉脊山斷裂西段應(yīng)引起關(guān)注和重視。從圖5所示的橫跨拉脊山斷裂帶的GPS水平運(yùn)動(dòng)剖面可以看出，1999～2007年垂直分量表現(xiàn)出明顯的縮短分量(圖5(a))，斷層兩側(cè)相對(duì)差異運(yùn)動(dòng)不明顯；2009～2014年垂直分量同樣表現(xiàn)出明顯的縮短分量(圖5(b))，斷層兩側(cè)表現(xiàn)出明顯的差異運(yùn)動(dòng)。西寧(XININ)點(diǎn)所處位置與斷層不遠(yuǎn)不近，反映斷層活動(dòng)的代表性較強(qiáng)。從圖5(b)可以看出，斷層北側(cè)與斷層南側(cè)地殼縮短速率相差較大，約3±0.5 mm/a。1999～2007年GPS水平運(yùn)動(dòng)剖面左旋走滑分量不明顯(圖5(c))；2009～2014年GPS水平運(yùn)動(dòng)剖面表現(xiàn)出明顯的左旋走滑分量(圖5(d))。結(jié)果表明，拉脊山地區(qū)仍舊處于地殼水平擠壓縮短狀態(tài)，2008年之后斷層兩側(cè)地殼縮短速率增加并表現(xiàn)出明顯的左旋走滑分量?？紤]到2008年汶川8.0級(jí)大震發(fā)生之后，又發(fā)生了若干6級(jí)左右以右旋走滑破裂為主的強(qiáng)余震。從機(jī)制上分析，這些強(qiáng)余震的效應(yīng)應(yīng)該是龍門山斷裂西側(cè)向北東向的推擠增強(qiáng)，從而引起拉脊山斷裂水平擠壓增強(qiáng)。

(a)和(c)為1999～2007年橫跨拉脊山斷裂的GPS運(yùn)動(dòng)速度，(b)和(d)為2009～2014年橫跨拉脊山斷裂的GPS運(yùn)動(dòng)速度圖5 橫跨拉脊山斷裂的GPS運(yùn)動(dòng)速度Fig.5 GPS velocity across the Lajishan fault

地質(zhì)考察和衛(wèi)星影像解譯表明，拉脊山斷裂西端是地震孕育的有利部位，應(yīng)關(guān)注該地區(qū)發(fā)生5～6級(jí)中等地震的危險(xiǎn)性[6]。本文根據(jù)橫跨拉脊山斷裂的水準(zhǔn)資料估算出拉脊山垂直運(yùn)動(dòng)速率為1±0.5 mm/a，表明拉脊山斷裂帶西段處于持續(xù)抬升階段。根據(jù)橫跨拉脊山斷裂的GPS水平運(yùn)動(dòng)速度剖面可以看出，拉脊山地區(qū)仍舊處于地殼水平擠壓縮短狀態(tài)，與衛(wèi)星影像解譯和野外地質(zhì)地貌調(diào)查結(jié)果一致；拉脊山斷裂西段左旋走滑活動(dòng)跡象明顯，以擠壓逆沖為主。水準(zhǔn)數(shù)據(jù)結(jié)果表明，拉脊山地區(qū)以抬升為主，與GPS數(shù)據(jù)顯示出的水平擠壓一致。

GPS水平運(yùn)動(dòng)剖面結(jié)果表明，2008年之后斷層兩側(cè)差異運(yùn)動(dòng)更加明顯，地殼縮短速率增加并表現(xiàn)出明顯的左旋走滑分量?？紤]為2008年汶川8.0級(jí)大震之后若干6級(jí)左右以右旋走滑破裂為主的強(qiáng)余震的影響，龍門山斷裂西側(cè)向北東向的推擠增強(qiáng)，從而引起拉脊山斷裂水平擠壓增強(qiáng)。晚更新世晚期以來(lái)拉脊山斷裂主斷裂活動(dòng)不明顯，斷裂活動(dòng)主要以隱伏斷裂和活動(dòng)褶皺的形式向兩側(cè)盆地內(nèi)部擴(kuò)展。在青藏高原北東向擴(kuò)展的背景下，未來(lái)拉脊山斷裂的活動(dòng)將促使兩側(cè)地層發(fā)生更深刻而劇烈的褶皺和斷裂變形。拉脊山斷裂西段地區(qū)受擠壓環(huán)境控制下的構(gòu)造變形現(xiàn)今仍在持續(xù)進(jìn)行，同時(shí)受全新世活動(dòng)的日月山斷裂右旋擠壓應(yīng)力的影響，該地區(qū)是地震孕育的有利部位，其地震危險(xiǎn)性應(yīng)該引起關(guān)注和重視。

4 結(jié) 語(yǔ)

根據(jù)拉脊山斷裂帶西段水準(zhǔn)數(shù)據(jù)，估算出該地區(qū)垂直運(yùn)動(dòng)速率為1±0.5 mm/a，表明該地區(qū)以抬升為主；與GPS數(shù)據(jù)顯示出的水平擠壓一致，表明拉脊山斷裂帶西段處于隆升狀態(tài)。橫跨拉脊山斷裂帶的GPS水平運(yùn)動(dòng)剖面可以看出，拉脊山地區(qū)仍舊處于地殼水平擠壓縮短狀態(tài)，2008年之后斷層兩側(cè)差異運(yùn)動(dòng)更加明顯，地殼縮短速率增加并表現(xiàn)出明顯的左旋走滑分量，與衛(wèi)星影像解譯和野外地質(zhì)地貌調(diào)查結(jié)果一致。拉脊山斷裂西段受擠壓環(huán)境控制下的構(gòu)造變形現(xiàn)今仍在持續(xù)進(jìn)行，同時(shí)受全新世活動(dòng)的日月山斷裂右旋擠壓應(yīng)力的影響，是地震孕育的有利部位，應(yīng)該關(guān)注和重視該地區(qū)的地震危險(xiǎn)性。

[1] 張旺生，馮光勝，高山，等. 拉脊山-化隆變質(zhì)核雜巖構(gòu)造及其隆升機(jī)制探討[J]. 地球科學(xué)，2003,28(4):407-413(Zhang Wangsheng, Feng Guangsheng, Gao Shan, et al. Metamorphic Core Complex Structure and Uplifting Mechanism in Lajishan-Hualong Area[J]. Earth Science, 2006,13(4):119-135)

[2] 袁道陽(yáng)，張培震，雷中生，等. 青海拉脊山斷裂帶新活動(dòng)特征的初步研究[J]. 中國(guó)地震，2005,21(1):93-102(Yuan Daoyang, Zhang Peizhen, Lei Zhongsheng, et al. A Preliminary Study on the New Activity Features of the Lajishan Mountain Fault Zone in Qinghai Province[J]. Earthquake Research in China, 2005,21(1):93-102)

[3] 袁道陽(yáng)，雷中生，劉小鳳，等.138年金城-隴西6.75級(jí)地震的史料考證與發(fā)震構(gòu)造背景探討[J].地震地質(zhì), 2004, 26(1): 52-60(Yuan Daoyang, Lei Zhongsheng, Liu Xiaofeng, et al. Textual Research on the Historical Data of the 138 AD Jincheng-Longxi M6.75 Earthquake and Discussion on Its Seimogenic Structural Background[J]. Seismology and Geology, 2004, 26(1): 52-60)

[4] 郝明. 基于精密水準(zhǔn)數(shù)據(jù)的青藏高原東緣現(xiàn)今地殼垂直運(yùn)動(dòng)與典型地震同震及震后垂直形變研究[D]. 北京：中國(guó)地震局地質(zhì)研究所, 2012(Hao Ming. Present Crustal Vertical Movement of Eastern Tibetan Plateau and Coseismic and Postseimic Vertical Deformation of Two Typical Earthquakes[J]. Beijing: Institute of Geology, CEA, 2012)

[5] 賈鵬，張希，張曉亮，等. 汶川地震前后川滇及其鄰區(qū)GPS水平運(yùn)動(dòng)演化特征[J]. 西北地震學(xué)報(bào)，2012,34(3):264-273(Jia Peng, Zhang Xi, Zhang Xiaoliang et al. Evolution Features of Crustal Horizontal Movement from GPS Observation in Sichuan-Yunnan Region and Its Adjacent Area before and after the Wenchuan Earthquake[J]. Northwestern Seismological Journal, 2012,34(3):264-273)

[6] 張波.西秦嶺北緣斷裂西端與拉脊山斷裂新活動(dòng)特征研究[D]. 蘭州:中國(guó)地震局蘭州地震研究所，2012(Zhang Bo. The Study of New Activies on Western Segment of Northern Margin of Western Qinling Fault and Lajishan Fault[D]. Lanzhou: Lanzhou Institute of Seismology,CEA, 2012)

[7] 周勇，王二七，李齊，等. 拉脊山及鄰區(qū)磁組特征及其地質(zhì)意義[J]. 中國(guó)科學(xué):D輯，2001,31(增1):187-194(Zhou Yong, Wang Erqi, Li Qi, et al. The Magnetic Characteristics and Geological Significance of Lajishan Area[J]. Science in China: Series D, 2001,31(Supp1):187-194)

[8] 李智敏，田勤儉，屠泓為. 拉脊山斷裂帶遙感特征研究[J]. 高原地震，2009, 21(1):26-31(Li Zhimin, Tian Qinjian, Tu Hongwei. Remote Sensing Characteristics of Lajishan Fault[J]. Plateau Earthquake Research, 2009, 21(1):26-31)

About the first author:ZHOU Lin, postgraduate,assistant engineer,majors in geodesy and geodynamics, E-mail: 429zl@163.com.

The Study of Crustal Deformation on Western End of Lajishan Fault Based on GPS and Leveling Data

ZHOULin1WANGQingliang1LIZhangjun1HAOMing1SONGShangwu2

1 Second Monitoring and Application Center, CEA, 316 Xiying Road , Xi’an 710054 ,China Discipline and Graduate Management Office,Institute of Disaster Prevention, 465 Xueyuan Street, Beijing 101601, China

The Lajishan area is an earthquake tectonic window to reflect on tectonic activity and earthquake activity. The Lajishan fault has caused about 20 moderate earthquakes on the both sides of Lajishan mountain in documented history. Lajishan fault is still an active fault zone, and has structural conditions of medium-strong earthquakes in the region. The Lajishan area is at 3 500 m above sea level and is not convenient for study. It remains doubtful whether the Lajishan regional crustal deformation is still in the uplift state. This paper uses 1983-1995 leveling data, combined with 1999-2007 and 2009-2014 GPS data, to analyze the three dimensional crustal deformation of this area. Research shows that the west section of Lajishan fault has a vertical rate of 1±0.5 mm/a ,which proves that the area is still in the uplift state. GPS data also show horizontal compression in this area. The Lajishan area remains a crustal horizontal shortening deformation. The crustal shortening rate of the region has become smaller, and shows an obvious left-lateral strike-slip component after 2008. The results are corroborated by satellite imagery interpretation and field geological survey results. Under the control of extrusion environment, tectonic deformation in the area is still ongoing. This area is also influenced by the Riyueshan right-lateral strike-slip fault. We should be concerned and put attention to seismic risk at the western end of Lajishan fault.

Lajishan fault; GPS; leveling measurement

The Spark Program of Earthquake Techndogy of CEA, No.XH15064Y; Special Fund for Earthquake Research of CEA,No. 201508009.

2015-11-16

項(xiàng)目來(lái)源：中國(guó)地震局地震科技星火計(jì)劃(XH15064Y)；中國(guó)地震局地震行業(yè)科研專項(xiàng)(201508009)。

周琳，碩士，助理工程師，研究方向?yàn)榇蟮販y(cè)量與地球動(dòng)力學(xué)，E-mail:429zl@163.com。

10.14075/j.jgg.2016.12.005

1671-5942(2016)012-1056-04

P315

A