利用GPS高頻信號分析瑪多地震產生的動態位移特征

張彩紅 譚 凱 王東振 李 瑜 趙 斌

1 中國地震局地震大地測量重點實驗室，武漢市洪山側路40號，430071 2 中國地震臺網中心，北京市三里河南橫街5號，100045

2021-05-22青海省瑪多縣發生MW7.3地震，震中位于34.613°N 、98.246°E，深度10 km。該地震引發近千次余震，破壞性極強。一個世紀以來，該地震震中距250 km以內已發生2次7級以上地震[1-2]，因此研究該地區因地震產生的位移變化尤為重要。斷層破裂的走向以及特征表明，瑪多地震是一個典型的東西向左旋走滑事件，這主要和青藏高原的南緣喜馬拉雅山一直受到印度洋板塊的俯沖有關，俯沖造成的位移每年高達4～5 cm[3]，誘發了多次地震。隨著GPS軟硬件的日益完善，科研工作者可以在地震后對震源周圍地區受到的影響給出快速響應[4-6]。本文在地震發生后利用Bernese GNSS Software Version 5.2 (簡稱Bernese5.2)[7]對震中附近的11個中國大陸構造環境監測網絡(以下簡稱陸態網)[8]測站1 Hz的GPS觀測數據進行處理分析，并評估地震對周圍地表的破壞性影響，為后續的相關研究提供參考。

1 數據的獲取

本文使用的11個陸態網測站GPS高頻觀測數據由中國地震臺網中心提供，測站全部位于震中距400 km范圍以內，分別為青海德令哈站(DLHA)、甘肅瑪曲站(GSMA)、青海班瑪站(QHBM)、青海都蘭站(QHDL)、青海格爾木站(QHGE)、青海瑪多站(QHMD)、青海瑪沁站(QHMQ)、青海同仁站(QHTR)、青海西寧站(XNIN)、西藏昌都站(XZCD)、四川甘孜站(SCGZ)，具體位置分布如圖1所示。其中QHMD站距離震中最近(37 km)，XZCD站距離震中最遠(396 km)。配置的接收機均為天寶NET8/9，除DLHA采用TRM59900.00 SCIS天線外，其他站均采用TRM59800.00 SCIS天線。由于瑪多縣西部海拔較高、地形復雜，GPS站建造條件受到限制，只設有QHGE一站，無法準確分析瑪多西部的地表形變。

圖1 高頻GPS站分布Fig.1 Distribution of the high rate GPS stations

2 數據解算策略

在解算高頻觀測數據時，采用瑞士伯爾尼大學研制開發的高精度數據處理軟件Bernese5.2。該軟件不僅能夠實現GPS高精度數據處理，服務于地殼形變研究，還能夠實現GPS高頻動態數據處理，獲得cm級的定位結果。為了提高解算效率并保證提取出地震引起的地表動態位移變化，本文解算地震發生時刻之前2 min和震后7 min的高頻觀測數據。由于精密星歷滯后兩周左右，為快速得到周圍地表動態位移變化，在解算過程中采用IGS 快速星歷。采用全球電離層模型進行高階電離層改正，正常解算時需要采用CODE數據中心公布的P1P2和P1C1差分碼偏差參數，但該參數按月公布，因此在解算過程中忽略差分碼偏差改正。采用維也納投影函數計算對流層天頂延遲，Bernese5.2雖然提供了單點定位動態解算模式，但是由于沒有對差分碼進行改正，并且采用了IGS快速星歷，一些共模誤差難以削弱。因此，本文基于基線模式進行單歷元動態定位，以距離震中超過2 000 km的IGS核心站——上海站(SHAO)為參考點，以STAR超長基線模式輻射出11條長度約為2 000 km的基線，從而減弱軌道誤差等共模誤差的影響，提高動態定位精度。解算時采用的具體模型及參數如表1所示。

表1 高頻動態定位采用的模型及參數

3 結果分析與討論

根據上述解算策略，約束遠離震中不受地震影響的GPS測站SHAO，獲得每個GPS觀測站2021-05-22 02：02～02：11 在ITRF2014參考框架下的高頻定位結果。其中東西向和南北向動態位移隨時間變化的序列如圖2所示，由于垂向高頻定位精度不高，本文只分析水平方向位移。圖中垂直于時間軸的紅色虛線為地震發生時刻，為了更清晰地表征各站點震中距的遠近和地震引起的地表動態位移大小，GPS測站從下到上按照震中距從近到遠的順序排列。

圖2 11個GPS測站東西向和南北向地表動態位移時間序列Fig.2 The time series of the eastern and northern dynamic displacements at the eleven GPS stations

由圖2(a)可見，瑪多地震對QHMD站東西向的影響最大，約為25 cm，對南北向的影響較小，約為10 cm。瑪多地震對周圍各站點東西向位移的影響總體上隨著震中距的增加而減小，但對位于震中西北方向QHGE站東西向位移的影響卻很小，這可能與青藏高原地下物質屬性不同從而影響了地震波的衰減速度有關。結合圖1可知，在圖1紅色虛線框內的GPS測站，其東西向位移會受到瑪多地震的影響，與震中距相關性較大、與方位角相關性較小。由圖2(b)可見，與東西向位移特征不同，瑪多地震對周圍地表各GPS站點南北向動態位移的影響隨震中距增大而衰減的現象并不明顯。結合圖1中各GPS站點和震中的相對位置關系可知，此現象與站點的方位角有很大的相關性。瑪多地震對QHMQ站和GSMA站的影響最為顯著，位移均為約27 cm。由圖2可知，無論東西方向還是南北方向，在地震發生之前，并沒有明顯的位移變化，地震使青海瑪多站的位置產生不可逆的變化，說明該地區的地下塑性物質在過程中起主導作用；而地震對其他GPS站點的影響能很快恢復到震前狀態，說明這些地區的地下彈性物質在過程中起主導作用。

4 結 語

本文利用超長基線模式對11個位于瑪多地震震中距400 km范圍內的陸態網高頻GPS觀測站的數據進行解算，得到地震前后共9 min的東西向和南北向動態位移時間序列。結果表明：1)地震使QHMD站發生不可逆的位移變化，說明地下的塑性物質在過程中起主導作用，而其他GPS站的位移在震后很快恢復到震前狀態，說明地下彈性物質在過程中起主導作用。2)瑪多地震對周圍各站位移產生的影響，尤其對南北向位移的影響，與站點方位角有較大的相關性。該研究可為今后瑪多地震形變分析提供參考。