GRACE和地面重力測(cè)量監(jiān)測(cè)到的中國(guó)大陸長(zhǎng)期重力變化

邢樂林，李 輝，玄松柏，汪 健

1 中國(guó)地震局地震研究所，武漢 430071

2 武漢大學(xué)測(cè)繪學(xué)院，武漢 430079

GRACE和地面重力測(cè)量監(jiān)測(cè)到的中國(guó)大陸長(zhǎng)期重力變化

邢樂林1，2，李 輝1，玄松柏1，2，汪 健1

1 中國(guó)地震局地震研究所，武漢 430071

2 武漢大學(xué)測(cè)繪學(xué)院，武漢 430079

自2002年以來，GRACE衛(wèi)星探測(cè)計(jì)劃可提供高精度的時(shí)變地球重力場(chǎng)，用以探測(cè)地球系統(tǒng)的物質(zhì)分布.自1998年中國(guó)大陸重力監(jiān)測(cè)網(wǎng)建立以來，利用FG5絕對(duì)重力儀和LCR－G型相對(duì)重力儀每2年對(duì)該網(wǎng)進(jìn)行重復(fù)測(cè)量獲取重力場(chǎng)時(shí)變信息.基于此，本文利用GRACE和地面重力測(cè)量獲得了中國(guó)大陸重力場(chǎng)的長(zhǎng)期年變率，利用位錯(cuò)理論根據(jù)USGS發(fā)布的斷層模型計(jì)算了2008年汶川Ms8.0級(jí)地震的同震重力變化并進(jìn)行了300km高斯濾波.GRACE衛(wèi)星重力和地面重力結(jié)果均表明華北地區(qū)地下水流失嚴(yán)重，在絕對(duì)重力基準(zhǔn)站上，GRACE衛(wèi)星重力與絕對(duì)重力變化率較為一致，汶川區(qū)域的地面重力變化結(jié)果可視為大地震前兆信息.

GRACE衛(wèi)星，地面重力測(cè)量，重力變化，地震，位錯(cuò)

1 引 言

地震的孕育和發(fā)展伴隨著構(gòu)造活動(dòng)、質(zhì)量遷移和密度變化等物理過程，從而引起地球重力場(chǎng)的非潮汐變化.自20世紀(jì)60年代拉科斯特重力儀出現(xiàn)后，許多國(guó)家開展了重力測(cè)量工作，并把觀測(cè)與地殼運(yùn)動(dòng)有關(guān)的區(qū)域重力場(chǎng)變化作為地震預(yù)測(cè)的一種手段［1］.研究表明，地震前后重力變化過程具有明顯的上升－下降－恢復(fù)特征，強(qiáng)震的孕育具有10年或更長(zhǎng)的時(shí)間尺度的應(yīng)力積累［2－4］.我國(guó)的流動(dòng)重力測(cè)量開始于邢臺(tái)地震后的20世紀(jì)60年代，在80和90年代取得較快的發(fā)展，特別90年代末期以來，隨著中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)、中國(guó)數(shù)字地震觀測(cè)網(wǎng)絡(luò)工程建設(shè)運(yùn)行以來，初步形成了覆蓋全國(guó)的基于絕對(duì)重力控制的相對(duì)重力聯(lián)測(cè)網(wǎng).網(wǎng)絡(luò)工程自1998年以來，每2～3年開展一期基準(zhǔn)站絕對(duì)重力測(cè)量和基本網(wǎng)相對(duì)重力聯(lián)測(cè)（1998，2000，2002和2005年），數(shù)字地震觀測(cè)網(wǎng)絡(luò)于2007年完成一期觀測(cè)［5－6］.

2002年3月美德聯(lián)合發(fā)射的GRACE（Gravity Recovery and Climate Experiment）衛(wèi)星為研究全球重力場(chǎng)及時(shí)間變化提供了精細(xì)信息.截至目前，GRACE提供了大量高精度觀測(cè)數(shù)據(jù)，為了解地球深部構(gòu)造及全球物質(zhì)分布、特別是季節(jié)性變化提供了重要參考依據(jù)［7］.過去由于GRACE衛(wèi)星重力資料的周期比較短，主要利用GRACE研究大尺度的季節(jié)信號(hào)變化，隨著資料的累計(jì)和數(shù)據(jù)處理方法的改進(jìn)，關(guān)注GRACE衛(wèi)星重力的長(zhǎng)期變化成為新的熱點(diǎn).

本文基于大空間尺度的GRACE衛(wèi)星重力和覆蓋中國(guó)大陸的地面重力測(cè)量數(shù)據(jù)，確定中國(guó)大陸重力場(chǎng)長(zhǎng)期變化特征，同時(shí)結(jié)合部分基準(zhǔn)站的長(zhǎng)期絕對(duì)重力時(shí)間序列，探索區(qū)域重力場(chǎng)變化機(jī)制，提取與汶川大地震相關(guān)的地震前兆信息.

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

GRACE重力衛(wèi)星計(jì)劃由美國(guó)NASA和德國(guó)DLR聯(lián)合開發(fā)，旨在獲取地球重力場(chǎng)中長(zhǎng)波部分及全球重力場(chǎng)的時(shí)變特征［8］，截至目前基于GRACE重力衛(wèi)星觀測(cè)數(shù)據(jù)已解出了9年多的月地球重力位模型.研究表明，利用GRACE位模型推算的地球表面質(zhì)量變化結(jié)果可以獲得高于750km空間分辨率的1.5cm等效水質(zhì)量變化［9］.因此諸如地下水儲(chǔ)量、冰川消融、海平面等變化，同震、震后形變及重力變化以及冰后回彈等地球物理效應(yīng)能夠被GRACE所檢測(cè).自2004年以來，大多數(shù)GRACE文獻(xiàn)主要集中在陸地水文或海洋變化［10－15］.除了與大地震相關(guān)的部分文獻(xiàn)外，涉及到的重力變化方面的較少［16－17］.隨著資料的累積增加以及數(shù)據(jù)處理技術(shù)的提高，科學(xué)家們開始關(guān)注長(zhǎng)期重力變化方面的研究.

本文使用的是CSR發(fā)布的GRACE level－2（Release－04）60階次的正則化月重力場(chǎng)模型，該版本與之前的版本比較，做了一些重要改進(jìn)，因此GRACE重力場(chǎng)中去除了大氣和海洋的影響.在研究區(qū)域內(nèi)，在2°×2°分辨率的格網(wǎng)點(diǎn)上根據(jù)每個(gè)月（2002年3月至2010年9月，共99個(gè)）重力場(chǎng)模型按公式（1）計(jì)算經(jīng)300km高斯平滑后的重力異常序列Δg［18］為

式（1）中，Re＝6378km為地球平均半徑，θ為余緯，λ為經(jīng)度，和為每個(gè)月重力場(chǎng)模型球諧系數(shù)與99個(gè)月重力場(chǎng)模型球諧系數(shù)均值之差，C20利用LAGEOS－1／2SLR數(shù)據(jù)處理結(jié)果，（cosθ）為完全正則化勒讓德函數(shù)，Wl為高斯平滑因子.

設(shè)年變率為B，振幅為Ci和Di的周期項(xiàng)ωi，則每個(gè)格網(wǎng)點(diǎn)的重力變化可通過式（2）表達(dá)［19］：

式（2）中，t為月重力場(chǎng)模型時(shí)間，i＝1和i＝2分別代表年周期項(xiàng)和半年周期項(xiàng)，且年周期項(xiàng)的振幅為代表與S2半日波相關(guān)的161天周期項(xiàng).

根據(jù)（1）、（2）式，由99個(gè)月重力場(chǎng)模型計(jì)算全球及中國(guó)大陸及周邊區(qū)域的長(zhǎng)期重力年變率，結(jié)果如圖1和圖2所示.

圖1 GRACE監(jiān)測(cè)到的全球重力年變率Fig.1 Secular gravity changes of global from GRACE

圖2 GRACE監(jiān)測(cè)到的中國(guó)大陸及鄰區(qū)重力年變率Fig.2 Secular gravity changes in Chinese mainland and its vicinity from GRACE

從圖1可看出，GRACE監(jiān)測(cè)到的較明顯的全球重力年變率主要體現(xiàn)了一些著名的區(qū)域，包括亞馬遜、剛果、密西西比、葉尼塞等大河流的水儲(chǔ)量變化；格陵蘭島、南極冰川和阿拉斯加冰川消融引起的重力效應(yīng)［20］；以及北加拿大冰后回彈效應(yīng).從圖2不難看出，中國(guó)大陸及其周邊區(qū)域含有5個(gè)比較明顯的重力變化特征：恒河流域上游和喜馬拉雅冰川消融區(qū)域，三峽大壩蓄水區(qū)域，三江源保護(hù)區(qū)，天山冰川消融區(qū)域和華北地區(qū).其中華北地區(qū)重力年變化率呈負(fù)變化，特別是京津冀地區(qū)變化尤為明顯，約為－0.3μGal／a的下降趨勢(shì)，由于華北地區(qū)重力觀測(cè)歷史資料較多，下文將進(jìn)一步進(jìn)行分析.

3 流動(dòng)重力數(shù)據(jù)處理

圖3 重力測(cè)站分布圖，籃點(diǎn)代表相對(duì)聯(lián)測(cè)點(diǎn)，紅星代表絕對(duì)基準(zhǔn)點(diǎn)Fig.3 The distribution of gravity stations.Blue circles represent relative stations and red stars represent absolute stations

為了獲得高精度的絕對(duì)重力變化，需要在重力網(wǎng)的絕對(duì)重力基準(zhǔn)站上進(jìn)行至少1個(gè)點(diǎn)次的絕對(duì)重力測(cè)量，同時(shí)使用相對(duì)重力儀進(jìn)行全網(wǎng)重力段差測(cè)量［21］.目前覆蓋我國(guó)大陸的流動(dòng)重力網(wǎng)有兩個(gè)：1998年開始運(yùn)行的中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)，由23個(gè)絕對(duì)重力基準(zhǔn)點(diǎn)和361個(gè)相對(duì)聯(lián)測(cè)點(diǎn)組成，已經(jīng)完成4期（1998、2000、2002、2005年）測(cè)量任務(wù)；2006年運(yùn)行的中國(guó)數(shù)字地震重力網(wǎng)，由10個(gè)絕對(duì)重力基準(zhǔn)點(diǎn)和453個(gè)聯(lián)測(cè)點(diǎn)組成，于2007年完成1期測(cè)量任務(wù).流動(dòng)重力網(wǎng)中絕對(duì)重力基準(zhǔn)點(diǎn)使用FG5絕對(duì)重力儀進(jìn)行測(cè)量，聯(lián)測(cè)點(diǎn)使用LCR－G型相對(duì)重力儀，形成絕對(duì)重力控制的具有時(shí)空基準(zhǔn)的重力網(wǎng)，如圖3所示.

3.1 絕對(duì)重力數(shù)據(jù)處理

FG5絕對(duì)重力儀的標(biāo)稱精度優(yōu)于5μGal，且不同儀器之間的較差約為1～2μGal，不存在明顯的系統(tǒng)偏差［22－23］.全國(guó)重力網(wǎng)中的絕對(duì)重力基準(zhǔn)站使用FG5絕對(duì)重力儀進(jìn)行絕對(duì)重力測(cè)量，每次觀測(cè)至少25h，每小時(shí)測(cè)量1組，每組100次下落，每臺(tái)儀器在每個(gè)測(cè)站的有效落體數(shù)不少于2400次下落，絕對(duì)重力值精度要求優(yōu)于5μGal，由于絕對(duì)重力基準(zhǔn)站地基穩(wěn)定，大多數(shù)測(cè)站的組內(nèi)精度約為7～15 μGal，點(diǎn)值精度優(yōu)于2μGal.在進(jìn)行絕對(duì)重力測(cè)量時(shí)，需要利用2臺(tái)LCR－G型相對(duì)重力儀進(jìn)行重力垂直梯度測(cè)量，為了保證測(cè)量結(jié)果的精度，在每個(gè)測(cè)站上地面與130cm高度處至少進(jìn)行5個(gè)高—低—高或低—高—低的往返測(cè)量.

為了得到高精度的地面絕對(duì)重力值，絕對(duì)重力測(cè)量數(shù)據(jù)處理軟件g允許對(duì)測(cè)試塊加速度進(jìn)行改正，主要包括大氣壓力、極移、基準(zhǔn)高度、潮汐等改正.

3.1.1 大氣壓力

在每個(gè)測(cè)站上基于在觀測(cè)期間測(cè)到的大氣壓力應(yīng)用改正使觀測(cè)的重力值被歸算到標(biāo)稱壓力，公式為

式（3）中A為大氣壓導(dǎo)納因子.該值通常在0.30～0.42之間，推薦值為0.3μGal／mBar，C（p）為以μGal為單位的大氣壓改正，P（o）為大氣壓觀測(cè)值，P（n）為標(biāo)稱大氣壓.

3.1.2 極移

此項(xiàng)改正補(bǔ)償由于重力測(cè)站到地球旋轉(zhuǎn)軸距離的改變而引起的離心加速度變化.改正量通常由測(cè)站上最接近觀測(cè)時(shí)間的極點(diǎn)位置計(jì)算得到.公式為：

式（4）中δg為以μGal為單位的極移改正，ω為地球旋轉(zhuǎn)角速度（rad／s），a為參考橢球體長(zhǎng)半軸（m），φ為大地緯度（rad），λ為大地經(jīng)度（rad），x，y為在IERS系統(tǒng)的極坐標(biāo)（rad）.

3.1.3 基準(zhǔn)高度

在落體倉(cāng)內(nèi)，重力值實(shí)際確定在落體倉(cāng)的頂部.對(duì)于每次下落的重力觀測(cè)值，轉(zhuǎn)換到用戶需要的特定高度.典型重力基準(zhǔn)高度是0cm、100cm或130cm.利用測(cè)量參考高度加出廠高度與基準(zhǔn)高度的差乘以測(cè)站重力梯度對(duì)重力值進(jìn)行調(diào)整，從而進(jìn)行轉(zhuǎn)換計(jì)算.

3.1.4 潮汐改正

計(jì)算日月引力對(duì)每次下落觀測(cè)的影響.重力因子通常取1.16.g提供Berger和ETGTAB2種模型的潮汐改正，與Berger模型相比，ETGTAB更為完善，通常采用ETGTAB模型進(jìn)行固體潮、海潮負(fù)荷改正.

3.1.5 統(tǒng)計(jì)

g自動(dòng)拒絕的限差默認(rèn)值是3，即任何相對(duì)平均值大于3σ的下落被拒絕.g對(duì)每組和最終工程的重力值都計(jì)算總的不確定度，公式為

3.2 相對(duì)重力數(shù)據(jù)處理

3.2.1 聯(lián)測(cè)誤差方程

相對(duì)重力網(wǎng)的平差以單程重力段差作為觀測(cè)值來建立平差模型，誤差方程為

式（6）中，vij為重力段差Δgij＝gi－gj的誤差，gi為測(cè)站i讀數(shù)的預(yù)處理重力值，ˉgi為測(cè)站i的平差值，D為儀器的漂移改正數(shù)，ti為i測(cè)站的觀測(cè)時(shí)刻.

3.2.2 基準(zhǔn)控制方程

設(shè)有起始重力點(diǎn)的絕對(duì)基準(zhǔn)平差，每個(gè)點(diǎn)的誤差方程為

式（7）中，gAi表示測(cè)點(diǎn)i的絕對(duì)重力值，vgi是絕對(duì)重力值的誤差改正數(shù).

3.2.3 觀測(cè)值權(quán)的確定

設(shè)某臺(tái)儀器有2個(gè)相鄰的重力段差為

即

設(shè)gi的方差為Dgr＝，gAi的方差為DgA＝，則基于絕對(duì)重力控制的相對(duì)聯(lián)測(cè)的協(xié)方差矩陣為

聯(lián)合（6）～（11）式分別對(duì)1998和2007年進(jìn)行重力網(wǎng)平差數(shù)據(jù)處理，各期點(diǎn)值平均精度優(yōu)于15μGal，2期差分即為基于地面重力測(cè)量獲得的中國(guó)大陸長(zhǎng)期重力變化，如圖4所示.

4 討論與分析

從圖4看出，中國(guó)大陸長(zhǎng)期重力變化呈現(xiàn)比較錯(cuò)亂的形態(tài)，各種形態(tài)應(yīng)該為不同地球物理信號(hào)的響應(yīng).

新疆西部：呈較明顯的正變化，與GRACE結(jié)果相反，主要體現(xiàn)了不同條件下的牛頓引力效應(yīng)［24］：天山區(qū)域雪融化的不斷累積在該區(qū)域?qū)е碌孛嬷亓χ翟黾?，該區(qū)域不呈現(xiàn)量級(jí)相同大小相反的圖像，因此可以不考慮地震孕育因素.

華北地區(qū)：呈較明顯的負(fù)變化，與GRACE結(jié)果相同，表明華北區(qū)域地下物質(zhì)大量流失或地面隆升，為了對(duì)引起重力變化的因素進(jìn)行提取，利用泰安和北京2個(gè)絕對(duì)重力基準(zhǔn)站的長(zhǎng)期觀測(cè)數(shù)據(jù)結(jié)合GRACE進(jìn)行比較和分析.

圖5表明泰安和北京2個(gè)基準(zhǔn)站的絕對(duì)重力值測(cè)定精度大部分優(yōu)于2.0μGal，地面絕對(duì)重力年變率均為負(fù)值，分別為－0.13μGal／a和－0.94μGal／a，GRACE衛(wèi)星重力年變率也均為負(fù)值，分別為－0.11μGal／a和－0.34μGal／a，而中國(guó)大陸地殼運(yùn)動(dòng)觀測(cè)網(wǎng)GPS的結(jié)果分別為－0.15cm／a和－0.03cm／a.假設(shè)測(cè)站地殼形變上升過程中，地表物質(zhì)同時(shí)上升，在區(qū)域半徑遠(yuǎn)大于發(fā)生形變的前提下，可利用Bouguer梯度來估算重力變化與形變過程中高程變化的關(guān)系，即

其中Gr＝－3.086μGal／cm，δ＝2.67g·cm－3為地球表層巖石的密度.GPS結(jié)果表明地面沉降，引起的重力效應(yīng)分別約為0.3μGal／a和0.06μGal／a，則地下物質(zhì)遷移引起的重力變化分別為－0.43μGal／a和－1.0μGal／a.據(jù)河北省水利部門公布的數(shù)據(jù)稱自1978年至今，地下水超采量達(dá)到1200億噸，約45億噸／年，區(qū)域陸地水量以大約2.4cm／a的速率下降［25］.因此，GRACE和利用消除高程因素影響的絕對(duì)重力長(zhǎng)期變化結(jié)果均表明，華北地區(qū)地下水流失嚴(yán)重.

龍門山斷裂帶區(qū)域：強(qiáng)地震給人類經(jīng)濟(jì)和社會(huì)發(fā)展造成了巨大的損失，為此，地震預(yù)測(cè)和抗震救災(zāi)成為當(dāng)今國(guó)家和各級(jí)政府部門與廣大地球物理工作者深切關(guān)注的科學(xué)問題.中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)和數(shù)字地震觀測(cè)網(wǎng)絡(luò)重力觀測(cè)數(shù)據(jù)的科學(xué)目標(biāo)是以地震預(yù)測(cè)預(yù)報(bào)為主，地震孕育過程中的地殼變形和介質(zhì)變化會(huì)引起震區(qū)周圍重力場(chǎng)變化，而強(qiáng)地震的應(yīng)力與能量累積具有10年及以上的時(shí)間尺度，因此利用長(zhǎng)期重力變化結(jié)果有可能為預(yù)測(cè)強(qiáng)震提供科學(xué)依據(jù).2008年汶川MS8.0地震發(fā)生在龍門山斷裂帶上，本文根據(jù)位錯(cuò)理論［26］按照USGS發(fā)布的斷層位錯(cuò)模型［27］計(jì)算汶川區(qū)域的重力變化，并對(duì)結(jié)果進(jìn)行300km空間高斯平滑，如圖6所示.

Sun給出了變形地球表面和空間固定點(diǎn)的同震重力變化計(jì)算公式，二者除了在數(shù)量級(jí)上有所不同外，符號(hào)基本相反［28］.2004年蘇門答臘大地震產(chǎn)生的重力變化被GRACE衛(wèi)星所檢測(cè)，其振幅為15 μGal［17］.從圖6b知，適用于GRACE衛(wèi)星觀測(cè)結(jié)果的理論重力變化振幅為2μGal，而GRACE對(duì)高頻重力變化不敏感，因此汶川同震重力變化很難被GRACE檢測(cè)到.分析其原因，可能包括以下幾個(gè)方面：汶川大震較Sumatra－Adaman大震的震級(jí)小且發(fā)震模式不同；大陸地震與海洋地震的差異；GRACE的空間分辨率和精度指標(biāo)尚未達(dá)到.

比較圖4與圖6a，在汶川區(qū)域的重力變化圖像相似，但符號(hào)相反.多個(gè)研究結(jié)果表明［1，4，6］，大地震發(fā)生的數(shù)年尺度內(nèi)重力呈正變化，地震發(fā)生后重力會(huì)產(chǎn)生較明顯的反向變化.從圖4不難看出，汶川恰在龍門山斷裂帶上，孕震區(qū)域所處的近北西向正負(fù)重力變化過度的梯度帶，形成較為明顯正負(fù)對(duì)稱重力變化，與理論同震重力變化較為接近.為了進(jìn)一步分析該區(qū)域的長(zhǎng)期重力變化，利用向上廷拓技術(shù)將重力變化廷拓至40km高度處，其結(jié)果反應(yīng)的是深度為20km地殼以下的質(zhì)量變化效應(yīng)［29］，扣除該變化后提取地殼構(gòu)造運(yùn)動(dòng)與質(zhì)量遷移引起的重力變化，如圖7所示，可以清晰地看到負(fù)正對(duì)稱但與圖6a符號(hào)相反的重力變化.汶川區(qū)域正處于較大規(guī)模的重力變化梯度帶，在重力變化正負(fù)過渡部位，能夠反映出地下物質(zhì)最大差異運(yùn)動(dòng).

圖6 同震重力變化：（a）地球表面，（b）空間固定點(diǎn)結(jié)果Fig.6 Co－seismic gravity changes：（a）is at the surface point and（b）at a space－fixed point

成都基準(zhǔn)點(diǎn)是距離汶川大地震震中最近的絕對(duì)重力測(cè)站.地震發(fā)生前后，利用FG5絕對(duì)重力儀進(jìn)行了多次絕對(duì)重力測(cè)量.這些觀測(cè)結(jié)果一方面為相對(duì)重力測(cè)量提供高精度絕對(duì)起始基準(zhǔn)，另一方面也為研究重力場(chǎng)的變化規(guī)律以及應(yīng)用于地球動(dòng)力學(xué)、環(huán)境監(jiān)測(cè)等研究積累了第一手觀測(cè)資料.從圖8看出，絕對(duì)重力變化呈增大趨勢(shì)直至2007年下半年（4.8μGal／a），特別是2002—2003年，絕對(duì)重力變化十分明顯（14.6μGal），但同震重力變化不明顯.震后絕對(duì)重力變化呈下降趨勢(shì)，至2010年9月與震前最后一次結(jié)果相比下降了約4μGal，但仍未恢復(fù).

5 結(jié)論與展望

GRACE和地面重力長(zhǎng)期觀測(cè)結(jié)果揭示了中國(guó)大陸重力場(chǎng)變化特征，在地基穩(wěn)定的基準(zhǔn)站上GRACE和地面絕對(duì)重力的長(zhǎng)期變化特征較為一致.GRACE全球重力場(chǎng)變化主要反映了空間質(zhì)量的分布，主要有大河流域的水儲(chǔ)量、冰山消融以及冰后回彈等地球物理現(xiàn)象.我國(guó)大陸陸地水量變化引起的重力變化，主要有華北地區(qū)、三峽大壩、黃河上游以及三江源，其中華北地區(qū)地下水抽取嚴(yán)重，引起的重力效應(yīng)約為－0.5μGal／a.天山冰川消融引起的地面重力變化較大，與衛(wèi)星結(jié)果相反.目前GRACE的精度難以探測(cè)到汶川Ms8.0地震引起的同震重力變化，地面重力長(zhǎng)期觀測(cè)結(jié)果與地球表面理論同震重力變化一致性較好，說明該地震孕育發(fā)生的前兆或影響范圍較大，震中西南部呈現(xiàn)正重力變化，東北部負(fù)重力變化，與發(fā)震斷層方向一致的較大規(guī)模重力變化梯度帶，成都絕對(duì)重力時(shí)序顯示的重力變化呈增大—加速增大—減速—發(fā)震—恢復(fù)的過程，與唐山地震孕育不同階段的重力場(chǎng)變化特征類似，為強(qiáng)震的中長(zhǎng)期預(yù)測(cè)提供了位置和時(shí)間的基本參考依據(jù).

2010年開始實(shí)施的中國(guó)大陸構(gòu)造環(huán)境監(jiān)測(cè)網(wǎng)絡(luò)是國(guó)家高技術(shù)產(chǎn)業(yè)發(fā)展計(jì)劃的12個(gè)國(guó)家重大科技基礎(chǔ)設(shè)施實(shí)施建設(shè)之一，其中高精度流動(dòng)重力網(wǎng)由608個(gè)相對(duì)重力聯(lián)測(cè)點(diǎn)和100個(gè)絕對(duì)重力基準(zhǔn)點(diǎn)構(gòu)成，將提供高質(zhì)量的中國(guó)大陸重力場(chǎng)動(dòng)態(tài)變化圖，同時(shí)在所有基準(zhǔn)站上進(jìn)行連續(xù)GPS觀測(cè)，具有空間高分辨率的多手段觀測(cè)網(wǎng)為研究中國(guó)大陸地殼運(yùn)動(dòng)與質(zhì)量遷移引起的重力變化及其與強(qiáng)震的關(guān)系提供了重要機(jī)遇.

（References）

［1］ Zhu Y Q，Zhan F B，Zhou J C，et al.Gravity measurements and their variations before the 2008Wenchuan earthquake.Bulletin of Seismological Society of America，2010，100（5B）：2851－2824.

［2］ 陳運(yùn)泰，顧浩鼎，盧造勛.1975年海城地震與1976年唐山地震前后的重力變化.地震學(xué)報(bào)，1980，2（1）：21－31.Chen Y T，Gu H D，Lu Z X.Variations of gravity before and after the Haicheng earthquake，1975and the Tangshan earthquake，1976.Acta Seismologica Sinica（in Chinese），1980，2（1）：21－31.

［3］ Li R H，F(xiàn)u Z Z.Local gravity variations before and after the Tangshan earthquake（M＝7.8）and the dilatation process.Tectonophysics，1983，97（1－4）：159－169.

［4］ 申重陽(yáng)，李輝，孫少安等.重力場(chǎng)動(dòng)態(tài)變化與汶川Ms8.0地震孕育過程.地球物理學(xué)報(bào)，2009，52（10）：2547－2557.Shen C Y，Li H，Su S A，et al.Dynamic variations of gravity and the preparation process of the Wenchuan Ms8.0earthquake.Chinese J.Geophys.（in Chinese），2009，52（10）：2547－2557.

［5］ 李輝，申重陽(yáng)，孫少安等.中國(guó)大陸近期重力場(chǎng)動(dòng)態(tài)變化圖像.大地測(cè)量與地球動(dòng)力學(xué)，2009，29（3）：1－10.Li H，Shen C Y，Sun S A，et al.Dynamic gravity change in recent years in China continent.Journal of Geodesy and Geodynamics（in Chinese），2009，29（3）：1－10.

［6］ 祝意青，王慶良，徐云馬.我國(guó)流動(dòng)重力監(jiān)測(cè)預(yù)報(bào)發(fā)展的思考.國(guó)際地震動(dòng)態(tài)，2008，（9）：19－25 Zhu Y Q，Wang Q L，Xu Y M.Thoughts on the development of earthquake monitoring and prediction in mobile gravity.Recent Developments in World Seismology（in Chinese），2008，9：19－25.

［7］ Cazenave A，Chen J L.Time－variable gravity from space and present－day mass redistribution in the Earth system.Earth and Planetary Science Letters，2010，298（3－4）：263－274.

［8］ Tapley B D，Bettadpur S，Watkins M，et al.The gravity recovery and climate experiment：mission overview and early results.Geophys.Res.Lett.，2004，31（9）：L09607.

［9］ Wahr J，Swenson S，Zlotnicki V，et al.Time－variable gravity from GRACE：First results.Geophys.Res.Lett.，2004，31（11）：L11501.

［10］ Tapley B D，Bettdpur S，Ries J C，et al.GRACE measurements of mass variability in the Earth system.Science，2004，305（5683）：503－505.

［11］ Chen J L，Wilson C R，F(xiàn)amiglietti J S，et al.Spatial sensitivity of the Gravity Recovery and Climate Experiment（GRACE）time－variable gravity observations.J.Geophys.Res.，2005，110：B08408.

［12］ Chen J L，Tapley B D，Wilson C R.Alaskan mountain glacial melting observed by satellite gravimetry.Earth and Planetary Science Letters，2006，248（1－2）：368－378.

［13］ Chen J L，Wilson C R，Tapley B D.Satellite gravity measurements confirm accelerated melting of Greenland ice sheet.Science，2006，313（5795）：1958－1960.

［14］ Chen J L，Wilson C R，F(xiàn)amiglietti J S，et al.Attenuation effect on seasonal basin－scale water storage changes from GRACE time－variable gravity.Journal of Geodesy，2007，81（4）：237－245.

［15］ Chen J L，Wilson C R，Tapley B D.The 2009exceptional Amazon flood and interannual terrestrial water storage change observed by GRACE.Water Resour.Res.，2010，46（12）：W12526.

［16］ Chen J L，Wilson C R，Tapley B D，et al.GRACE detects coseismic and postseismic deformation from the Sumatra－Andaman earthquake.Geophys.Res.Lett.，2007，34：L13302.

［17］ Han S C，Shum C K，Bevis M，et al.Crustal dilatation observed by GRACE after the 2004Sumatra－Andaman earthquake.Science，2006，313（5787）：658－662.

［18］ Jekeli C.Alternative methods to smooth the Earth′s gravity field.Rep.327，Dep.of Geod，Departhment of Geodetic Science and Surveying，The Ohio State University，Columbus，1981.

［19］ Steffen H，Gitlein O，Denker H，et al.Present rate of uplift in Fennoscandia from GRACE and absolute gravimetry.Tectonophysics，2009，474（1－2）：69－77.

［20］ Sun W，Miura S，Sato T，et al.Gravity measurements in southeastern Alaska reveal negative gravity rate of change caused by glacial isostatic adjustment.J.Geophys.Res.，2010，115：B12406.

［21］ Hwang C，Cheng T C，Cheng C C，et al.Land subsidence using absolute and relative gravimetry：a case study in central Taiwan.Survey Review，2010，42（315）：27－39.

［22］ Xing L L，Li H，Li J C，et al.Comparison of absolute gravity measurements obtained with FG5／232and FG5／214 instruments.Geo－spatial Information Science，2009，12（4）：307－310.

［23］ 邢樂林，申重陽(yáng)，李輝等.歐洲絕對(duì)重力儀比對(duì)觀測(cè)（ECGS′07）.地球物理學(xué)進(jìn)展，2009，24（6）：2054－2057.Xing L L，Shen C Y，Li H，et al.European comparison of absolute gravimeters（ECGS′07）.Progress in Geophys.（in Chinese），2009，24（6）：2054－2057.

［24］ Boy J P，Hinderer J.Study of the seasonal gravity signal in superconducting gravimeter data.Journal of Geodynamics，2006，41（1－3）：227－233.

［25］ Zhong M，Duan J B，Xu H Z，et al.Trend of China land water storage redistribution at medi－and large－spatial scales in recent five years by satellite gravity observations.Chinese Science Bulletin，2009，54（5）：816－821.

［26］ Sun W K，Okubo S.Surface potential and gravity changes due to internal dislocations in a spherical Earth，II.Application to a finite fault.Geophys.J.Int.，1998，132（1）：79－88.

［27］ Ji C，Hayes G.Preliminary result of the May 12，2008 Mw7.9Eastern Sichuan，China Earthquake.http：／／earthquake.usgs.gov／eqcenter／eqinthenews／2008／us2008ryan／finite＿fault.php.

［28］ Sun W K，Okubo S.Coseismic deformations detectable by satellite gravity missions：a case study of Alaska（1964，2002）and Hokkaido（2003）earthquakes in the spectral domain.J.Geophys.Res.，2004，109：B04405.

［29］ Meng X H，Guo L H，Chen Z X，et al.A method for gravity anomaly separation based on preferential continuation and its application.Applied Geophysics，2009，6（3）：217－225.

Long－term gravity changes in Chinese mainland from GRACE and terrestrial gravity measurements

XING Le－Lin1，2，LI Hui1，XUAN Song－Bai1，2，WANG Jian1
1 Instiute of Seismology，China Earthquake Administration，Wuhan 430071，China
2 School of Geodesy and Geomatics，Wuhan University，Wuhan 430079，China

Since 2002，the GRACE satellite mission has been providing precise survey data of the Earth′s time－variable gravity field，and has greatly improved the understanding of mass distribution on and near the surface of the earth.A gravity network in Chinese mainland is established to determine gravity variations using LCR－G relative gravimeter and FG5absolute gravimeter about every 2years since 1998.Based on these observations，the secular trend of gravity changes in Chinese mainland is obtained and analyzed，and co－seismic gravity changes caused by the 2008Wenchuan Ms8.0earthquake are computed using the dislocation model of USGS and treated with 300km Gaussian filter.Both GRACE and terrestrial observation results show that in the north part of China，the features have strongest resemblance，indicating the presence of large－scale changes in ground water storage，at the absolute gravity datum stations，the features of the gravity variations observed by GRACE and the absolute gravity measurements agree reasonably well，the terrestrial gravity changes around Wenchuan area could be viewed as aprecursor of the great earthquake.

GRACE satellites，Terrestrial gravity measurements，Gravity change，Earthquake，Dislocation

10.6038／j.issn.0001－5733.2012.05.013

P228

2011－09－26，2012－03－28收修定稿

國(guó)家自然科學(xué)基金項(xiàng)目（41004030）資助.

邢樂林，男，助理研究員，主要從事地震重力觀測(cè)與解釋研究.E－mail：xinglelin＠163.com

邢樂林，李輝，玄松柏等.GRACE和地面重力測(cè)量監(jiān)測(cè)到的中國(guó)大陸長(zhǎng)期重力變化.地球物理學(xué)報(bào)，2012，55（5）：1557－1564，

10.6038／j.issn.0001－5733.2012.05.013.

Xing L L，Li H，Xuan S B，et al.Long－term gravity changes in Chinese mainland from GRACE and terrestrial gravity measurements.Chinese J.Geophys.（in Chinese），2012，55（5）：1557－1564，doi：10.6038／j.issn.0001－5733.2012.05.013.

（本文編輯 胡素芳）