康開軒 李 輝 劉少明 郝洪濤 鄒正波

1 中國地震局地震研究所（地震大地測量重點實驗室），武漢市洪山側路40號，430071

2 中國地震局地殼應力研究所科技創(chuàng)新基地，武漢市洪山側路40號，430071

地震前后重力變化過程具有明顯的上升－下降－恢復特征，而強震的孕育具有10a或更長時間尺度的應力積累［1］。2015－04－25尼泊爾地震屬于印度板塊俯沖、碰撞歐亞大陸的動力作用結果，其孕震尺度的空間范圍必然涉及印度板塊與青藏高原的較大區(qū)域。本文搜集了尼泊爾Ms8.1地震前15a間我國西藏及周邊地區(qū)的8期流動重力觀測資料（圖1），獲取該區(qū)域長期重力變化分布特征，并初步探討此次尼泊爾Ms8.1地震孕育的重力場較長期或背景性效應。

圖1 研究區(qū)構造及重力測點分布Fig.1 Location of the relative and absolute gravity stations and faults geometry

1 數(shù)據(jù)處理與結果

本文綜合使用地殼重力網(wǎng)5期、數(shù)字地震重力網(wǎng)1期和陸態(tài)重力網(wǎng)2期的觀測資料，獲取西藏及周邊地區(qū)1998～2013每1～2a的相對重力變化。首先，為綜合使用不同觀測網(wǎng)絡的資料，對歷史數(shù)據(jù)的測點編號依據(jù)現(xiàn)行統(tǒng)一的編碼規(guī)范進行清理；其次，利用WDD99 固體潮模型、基于Ferrell格林函數(shù)的海潮負荷模型及基于導納原理的氣壓和溫度改正模型，去除原始觀測中的主要重力潮汐信號及由環(huán)境因素引起的區(qū)域重力非潮汐效應；最后，每期資料采用絕對重力基準控制下的相對重力聯(lián)測的經(jīng)典間接平差原理，建立基于重力段差觀測值的平差函數(shù)模型和隨機模型，利用最小二乘約束獲得測網(wǎng)中各測點的平差最優(yōu)解。相對重力聯(lián)測網(wǎng)平差誤差方程為：

式中，vij為Δgij＝gi－gj的誤差，gi為i測點上讀數(shù)的預處理重力值，為測點i的平差值，D為儀器的漂移改正數(shù)，ti為測點i的觀測時刻。

絕對重力觀測結果以5μGal精度定權，相對重力聯(lián)測觀測值的權重則按儀器型號的不同，在解算過程中迭代計算確定。相對重力儀彈簧零漂改正在網(wǎng)平差計算中作為未知參數(shù)通過整體平差解算，以獲取合理的儀器漂移系數(shù)及漂移改正量。重力儀格值系數(shù)先作為未知參數(shù)在相關期的絕對重力測量控制下整體解算，獲取合適的格值系數(shù)改正值，再對原始數(shù)據(jù)作相應的格值系數(shù)改正［2］。

測區(qū)中各相對聯(lián)測點的相對重力變化由每期平差結果減去多期結果平均值獲得，重力年變率由多期平差結果加權最小二乘回歸擬合獲得，權重采用各期數(shù)據(jù)平差后的測點重力值精度，年變率的精度由1－σ的標準差估算，相當于67%的置信度。解算結果表明，測區(qū)內重力年變率分布于－17～20μGal／a，平均年變率約為3.7μGal／a；年變率精度優(yōu)于10μGal／a，分布于0.4～9.8μGal／a，大部分集中在3μGal／a附近，少數(shù)復測次數(shù)較少的測點年變率精度較低。拉薩絕對重力點年變率為－2.66±0.75μGal／a，與文獻［3－5］的結果在給定誤差范圍內是一致的。圖2給出了精度優(yōu)于5μGal／a的重力年變率分布，可以看出，青藏高原西南邊緣喜馬拉雅地區(qū)及北部可可西里山－阿爾金山地區(qū)呈現(xiàn)明顯的正重力變化特征，最大重力變化約11.21μGal／a。圖3給出了獅泉河、于田等典型測點的時間序列，可以看出，這些測點的時間序列均表現(xiàn)出明顯的線性變化特征，表明測點附近區(qū)域的長期重力變化在15a時間尺度上以線性上升趨勢為主，反映了該區(qū)域逐年累積增加的重力背景效應。

圖2 1998～2013我國西藏及周邊區(qū)域的重力變化率Fig.2 The gravity annual rate in Tibet and its vicinity in the period of 1998－2013

圖3 典型測點的時間序列Fig.3 The time series at typical gravity stations

2 討論與分析

2.1 陸地水儲量變化引起的區(qū)域重力效應

2.1.1 冰川儲量變化

近50a來中國西部82.2%的冰川處于退縮和消融狀態(tài)，冰川面積減少了4.5%，但退縮的幅度和比例不同，存在明顯的區(qū)域差異：冰川末端退縮、物質平衡變化和面積變化的幅度在青藏高原邊緣和周邊地區(qū)要大于高原內部［6］。由圖1 可見，我國西藏及周邊區(qū)域的流動重力測點大多集中在高原內部，只有少數(shù)測點離冰川分布區(qū)較近。孫文科［4］認為，就長期重力變化而言，由冰川融化和高原降水引起的重力變化率應該是較小的，并且這些參數(shù)中的季節(jié)性變化可以通過長期觀測取平均消除，它們對重力變化率的貢獻同全球水平衡有關，應是一個微小的量。

2.1.2 高原湖泊水量變化

青藏高原高山地區(qū)冰川補給占主導的湖泊集中分布于藏南喜馬拉雅地區(qū)和藏東南地區(qū)，這些冰湖近40a來面積和水位總體上呈迅速增加態(tài)勢［7］。Zhang［8］對青藏高原111個湖泊2003～2009的高程變化研究表明，湖泊亞區(qū)Ⅰ、Ⅱ、Ⅲ、Ⅳ及整個青藏高原的湖泊水位變化速率分別為0.12、0.26、0.19、－0.11和0.2m／a（圖4）。近30a來我國喜馬拉雅山區(qū)冰湖總面積增大了29%，區(qū)域性的水量變化與湖泊面積變化相對應，藏北的湖泊水量以增加為主，高山地區(qū)冰川補給占主導的湖泊水位上升非常迅速［8－9］。湖泊水量的逐年增加是區(qū)域重力場變化的主要因素之一。比較圖2、圖4可知，研究區(qū)內重力變化率為正值的測點大多分布在湖泊水位上升的區(qū)域，這表明地面重力觀測的重力變化信號中包含了測點周圍湖泊水量逐年增加引起的重力變化。為研究與地震相關的地殼構造運動及深部物質運移的重力信號，必須扣除陸地水儲量變化引起的重力效應。精細扣除這部分水負荷效應需要測點附近河流、湖泊、地下水和常年降水量等水文氣象資料，由于目前缺乏上述資料，本研究將利用Zhang［8］給出的青藏高原主要湖泊區(qū)水位變化率，基于地表負荷響應的理論計算測區(qū)內湖泊水量變化引起的重力效應，并估算測區(qū)內各重力測點的水負荷重力效應改正量。

圖4 青藏高原74個湖泊水位變化［8］Fig.4 Water level change of 74lakes of the Tibetan plateau［8］

陸地水儲量變化對重力場觀測的影響可分為3部分：1）陸地水質量變化引起的直接牛頓引力效應；2）陸地水質量負荷作用下的彈性地球形變效應；3）由于地球形變使地球內部質量重新分布產(chǎn)生的附加效應［10］。根據(jù)地表負荷響應的理論計算方法，P點t時刻的重力變化δg可由下列公式表示［11］：

其中，gN為直接牛頓引力效應；gE與格林函數(shù)相關，代表地球彈性變化引起的重力響應；地表質量負荷變化δg（P，t）與湖泊水位變化l（p，t）及水密度ρw＝1 000kg／m3直接相關，即δg（P，t）＝ρwl（p，t）；me為地球質量；g為地表平均絕對重力值；Pn（cosθ）為n階勒讓德多項式；hn、kn為負荷Love數(shù)；θ為觀測點P與負荷單元p之間的角距離。

圖5 扣除高原湖泊水負荷重力效應后的重力年變率分布Fig.5 The gravity annual rate with the hydrological loading effect removed

圖5為扣除湖泊水負荷重力效應的各測點重力年變率分布。比較圖2、圖5可以看出，經(jīng)湖泊水負荷改正后的重力年變率變化較大的測點均為靠近高原湖泊的測點。

2.2 尼泊爾地震震前我國西藏及周邊區(qū)域的重力長期變化

巨大地震的孕育和發(fā)生過程中在震源區(qū)附近會產(chǎn)生物質遷移和質量再分布，重力場也會隨之發(fā)生變化。青藏高原巖石圈活動的重力效應主要表現(xiàn)為兩種機制：地表隆升和地下深部物質運移［3，5，12］。基于平均球形地球模型的垂直重力梯度約－3.08μGal／cm，結合近期多位學者［13－14］給出的青藏高原地區(qū)現(xiàn)今地殼垂直運動結果，用自由空氣改正估算由地表垂直位移引起的重力變化約為－0.6～－1.5μGal／a。同地面觀測到的重力變化率（圖5）相比可知，由地殼垂直位移引起的重力變化量級很小。由此可以推斷，在扣除陸地水負荷重力效應后，本文基于地面重力觀測得到的西藏及周邊地區(qū)的重力變化主要反映了該地區(qū)深部殼幔物質遷移和質量再分布。

圖6為利用陸地水負荷改正后各測點的重力年變率經(jīng)50km 平滑半徑的高斯濾波后的重力長期變化空間分布。圖6表明，我國西藏及周邊區(qū)域的重力長期變化在空間分布上具有顯著的不均勻性和重力變化分區(qū)現(xiàn)象，這與青藏高原復雜的變形構造和動力學系統(tǒng)密切相關。徐紀人［15］關于青藏高原構造運動特征與地震震源斷層活動類型的研究表明，在青藏高原周緣海拔較低的地區(qū)，地震主要以逆斷層型或者走滑斷層型活動發(fā)生，顯示出在印度板塊向歐亞大陸擠壓碰撞的作用下青藏高原不斷隆升的構造運動特征；而在青藏高原中部高海拔地區(qū)構造運動趨勢相反，促使不斷隆升的青藏高原基本保持著重力均衡狀態(tài)。圖6給出的15a時間尺度長期重力變化分布表明，青藏高原周緣地區(qū)，尤其是青藏高原南段的喜馬拉雅活動構造帶和北段的阿爾金－昆侖斷裂帶呈明顯的正重力變化趨勢，變化范圍約為2～12 μGal／a；高原中部的藏北地塊及青南地塊為大范圍的負重力變化區(qū)，變化率集中分布在－2～－8 μGal／a。地面重力觀測給出的青藏高原周緣地區(qū)和高原中部地區(qū)截然相反的重力長期背景變化趨勢，可能與徐紀人［15］指出的青藏高原周緣與中部地區(qū)相反趨勢的構造運動有關。占偉等［16］的研究表明，喜馬拉雅主邊界斷裂存在大范圍擠壓應變積累，尼泊爾地震震源區(qū)處于近南北向擠壓應變累積高值區(qū)的過渡區(qū)域。一般來說，區(qū)域物質運動質量的總體增加，其所攜帶的能量亦增加，有利于地震孕育積累。地震源于地殼深處，是介質應力應變能量積累到一定程度而造成的剪切破裂。大震的孕育往往伴隨著能量的積累，比如地殼內部物質不斷積聚變形運動，地表重力變化上升或正重力變化應是有利于能量積累的可能標志之一［17］。圖6給出的喜馬拉雅活動構造帶顯著的正重力變化趨勢，反映了尼泊爾Ms8.1地震孕育的中長期時間尺度（15a）信息。大范圍的正重力變化趨勢可能和印度板塊與歐亞板塊存在的持續(xù)擠壓變形引起地下物質重新分布與調整有關，反映了大震孕育過程中地殼變形和介質變化引起的震區(qū)周圍應力與能量的累積。

圖6 1998～2013我國西藏及周邊區(qū)域的重力長期變化空間分布Fig.6 The spatial distribution of the gravity annual rate in Tibet and its vicinity in the period of 1998－2013

3 結 語

流動重力典型測點的時間序列表現(xiàn)出明顯的線性變化特征，表明測點附近區(qū)域的長期重力變化在15a時間尺度上以線性趨勢為主，反映了區(qū)域重力逐年累積增加的背景效應。我國西藏及周邊區(qū)域的重力長期變化在空間分布上具有顯著的不均勻性和重力變化分區(qū)現(xiàn)象，這與青藏高原復雜的變形構造和動力學系統(tǒng)密切相關。喜馬拉雅活動構造帶在15a時間尺度上呈明顯的正重力變化趨勢，可能與印度板塊與歐亞板塊存在的持續(xù)擠壓變形引起的地下物質重新分布與調整有關，反映了大震孕育過程中地殼變形和介質變化引起的震區(qū)周圍應力與能量的累積。

本文中陸地水儲量變化只考慮了冰川消融和高原湖泊水位變化引起的重力效應，更為嚴密的陸地水模型仍需考慮高原凍土、土壤含水層、地下水儲量和年季降水變化等因素。除此之外，青藏高原冰后期回彈效應和地表剝蝕等因素也會引起重力場的變化，在后續(xù)研究中應予以詳細探討。

［1］李瑞浩.唐山地震前后區(qū)域重力場變化機制［J］.地震學報，1997，19（4）：399－407（Li Ruihao.The Regional Gravity Variation Mechanism before and after Tangshan Earthquake［J］.Acta Seismologica Sinica，1997，19（4）：399－407）

［2］李輝，劉冬至，劉紹府.地震重力監(jiān)測網(wǎng)統(tǒng)一平差模型的建立［J］.地殼形變與地震，1991，11（增刊）：68－74（Li Hui，Liu Dongzhi，Liu Shaofu.Integtated Adjustment Models for the Seismic－Gravity Network［J］.Crustal Deformation and Earthquake，1991，11（Supp）：68－74）

［3］Sun W K，Wang Q，Li H，et al.Gravity and GPS Measurements Reveal Mass Loss Beneath the Tibetan Plateau：Geodetic Evidence of Increasing Crustal Thickness［J］.Geophysical Research Letters，2009，36（2）：206－218

［4］孫文科，長谷川崇，張新林，等.高斯濾波在處理GRACE數(shù)據(jù)中的模擬研究：西藏拉薩的重力變化率［J］.中國科學：地球科學，2011，41（9）：1 327－1 333（Sun Wenke，Hasegawa T，Zhang Xinlin，et al.Effects of Gaussian Filter in Processing GRACE Data：Gravity Rate of Change at Lhasa，Southern Tibet［J］.Sci China Earth Sci，2011，54：1 378－1 385）

［5］邢樂林，孫文科，李輝，等.用拉薩點大地測量資料檢測青藏高原地殼的增厚［J］.測繪學報，2011，40（1）：41－58（Xing Lelin，Sun Wenke，Li Hui，et al.Present－Day Crust Thickness Increasing Beneath the Qinghai－Tibetan Plateau by Using Geodetic Data at Lhasa Station［J］.Acta Geodaetica et Cartographica Sinica，2011，40（1）：41－58）

［6］段建平，王麗麗，任賈文，等.近百年來中國冰川變化及其對氣候變化的敏感性研究進展［J］.地理科學進展，2009，28（2）：231－237（Duan Jianping，Wang Lili，Ren Jiawen，et al.Progress in Glacier Variations in China and Its Sensitivity to Climatic Change during the Past Century［J］.Progress in Geography，2009，28（2）：231－237）

［7］李治國.近50a氣候變化背景下青藏高原冰川和湖泊變化［J］.自然資源學報，2012，27（8）：1 431－1 443（Li Zhiguo.Glaciers and Lakes Changes on the Qinghai－Tibet Plateau under Climate Change in the Past 50 Years［J］.Journal of Natural Researches，2012，27（8）：1 431－1 443）

［8］Zhang G Q，Xie H J，Kang S C.Monitoring Lake Level Change on the Tibetan Plateau Using ICES at Altimetry Data（2003－2009）［J］.Remote Sensing of Environment，2011，115（7）：1 733－1 742

［9］Matsuo K，Heki K，Time－Variable Ice Loss in Asian High Mountains from Satellite Gravimetry［J］.Earth and Planetary Science Letters，2010，290：30－36

［10］Pfeffer J，Boucher M，Hinderer J，et al.Local and Global Hydrological Contributions to Time－Variable Gravity in Southwest Niger［J］.Geophs J Int，2011，184：661－672

［11］Farrell W E.Deformation of the Earth by Surface Loads［J］.Rev Geophys，1972，10（3）：761－797

［12］滕吉文，阮小敏，張永謙.青藏高原地殼與上地幔成層速度結構與深部層間物質的運移軌跡［J］.巖石學報，2012，28（12）：4 077－4 100（Teng Jiwen，Ruan Xiaomin，Zhang Yongqian，et al.The Stratificational Velocity Structure of Crust and Covering Strata of Upper Mantle and the Orbit of Deep Interaquifer Substance Locus of Movement for Tibetan Plateau［J］.Acta Petrologica Sinica，2012，28（12）：4 077－4 100）

［13］王敏.GPS 觀測結果的精化分析與中國大陸現(xiàn)今地殼形變場研究［D］.北京：中國地震局地質研究所，2009（Wang Min.Analysis of GPS Data with High Precision and Study on Present－Day Crustal Deformation in China［D］.Beijing：Institute of Geology，CEA，2009）

［14］王偉，楊少敏，趙斌，等.中國大陸現(xiàn)今地殼運動速度場［J］.大地測量與地球動力學，2012，32（6）：29－32（Wang Wei，Yang Shaomin，Zhaobin，et al.Present－Day Crustal Movement Velocity Field in Chinese Mainland［J］.Journal of Geodesy and Geodynamics，2012，32（6）：29－32）

［15］徐紀人，趙志新.青藏高原及其周邊地區(qū)區(qū)域應力場與構造運動特征［J］.中國地質，2006，33（2）：275－285（Xu Jiren，Zhao Zhixin.Characteristics of the Regional Stress Field and Tectonic Movement on the Qinghai－Tibet Plateau and in Its Surrounding Areas［J］.Geology in China，2006，33（2）：275－285）

［16］占偉，武艷強，梁洪寶，等.GPS 觀測結果反映的尼泊爾Mw7.8 地震孕震特征［J］.地球物理學報，2015，58（5）：1 818－1 826（Zhan Wei，Wu Yanqiang，Liang Hongbao，et al.Characteristics of the Seismogenic Model for the 2015 Nepal Mw 7.8Earthquake Derived from GPS Data［J］.Chinese J Geophys，2015，58（5）：1 818－1 826）

［17］申重陽，李輝，孫少安，等.2008年于田Ms7.3地震前重力場動態(tài)變化特征分析［J］.大地測量與地球動力學，2010，30（4）：1－7（Shen Chongyang，Li Hui，Sun Shaoan，et al.Characteristic Analysis of Dynamic Gravity Change before Yutian Ms 7.3Earthquake，2008［J］.Journal of Geodesy and Geodynamics，2010，30（4）：1－7）