焦佳爽, 張永志, 張凱南, 姜永濤, 槐巖珂

(1.長(zhǎng)安大學(xué)地質(zhì)與測(cè)繪工程學(xué)院,陜西 西安 710000; 2.南陽(yáng)師范學(xué)院,河南 南陽(yáng) 473061)

?

云南地區(qū)及鄰域衛(wèi)星重力變化及其與區(qū)域大震的關(guān)系①

焦佳爽1, 張永志1, 張凱南1, 姜永濤2, 槐巖珂1

(1.長(zhǎng)安大學(xué)地質(zhì)與測(cè)繪工程學(xué)院,陜西 西安 710000; 2.南陽(yáng)師范學(xué)院,河南 南陽(yáng) 473061)

摘要:利用GRACE衛(wèi)星重力場(chǎng)數(shù)據(jù),計(jì)算并分析云南地區(qū)及鄰域2003—2014年的衛(wèi)星重力場(chǎng)變化特征,獲得3個(gè)特征點(diǎn)上的月重力變化時(shí)間序列;結(jié)合該地區(qū)2003—2014年發(fā)生的MS5以上地震分布,探討區(qū)域衛(wèi)星重力變化與區(qū)域大震的關(guān)系。結(jié)果表明:云南地區(qū)及鄰域衛(wèi)星重力場(chǎng)變化對(duì)區(qū)域MS5以上地震有一定的反映,地震多發(fā)生在區(qū)域重力變化的高梯度帶上,地震發(fā)生前后,區(qū)域出現(xiàn)與地震相對(duì)應(yīng)的特征性衛(wèi)星重力變化,可能是對(duì)殼幔物質(zhì)調(diào)整的反映;區(qū)域大震前后,不同震區(qū)特征點(diǎn)的重力變化特征存在差異,說(shuō)明分?jǐn)鄬?、分區(qū)域的重力場(chǎng)變化研究對(duì)地震預(yù)測(cè)有重大意義。

關(guān)鍵詞:云南地區(qū)及鄰域; GRACE; 衛(wèi)星重力變化; 地震; 斷裂帶

0引言

云南緊鄰印度板塊向歐亞板塊中國(guó)大陸的側(cè)向俯沖帶,且處于我國(guó)南北地震帶南段,已有研究表明,該區(qū)域地震活動(dòng)與重力場(chǎng)變化特征密切相關(guān)[1-6]。

GRACE(Gravity Recovery And Climate Experiment)雙星系統(tǒng)由美國(guó)宇航局NASA和德國(guó)空間飛行中心DLR聯(lián)合研制,其提供的時(shí)變重力場(chǎng)是探究地球內(nèi)部物質(zhì)遷移的重要依據(jù)。GRACE能探測(cè)到MS9.0以上剪切型地震和MS7.5以上張裂型地震產(chǎn)生的同震重力變化[7]。2006年研究人員利用GRACE數(shù)據(jù)首次探測(cè)到2004年蘇門(mén)答臘MS9.1地震的同震重力變化,其峰值達(dá)±15×10-8ms-2[8]。此后有部分學(xué)者利用GRACE時(shí)變重力場(chǎng)模型,研究了與地震相關(guān)的重力變化。周新等[9]對(duì)GRACE數(shù)據(jù)進(jìn)行了有效處理,成功提取出2010年智利MW8.8地震的同震重力變化信號(hào),表明GRACE具有檢測(cè)出M<9 地震量級(jí)的能力;段虎榮等[10]利用GRACE衛(wèi)星數(shù)據(jù)研究了汶川地震前后的重力場(chǎng)變化,并用地殼膨脹模型理論解釋了同震衛(wèi)星重力的負(fù)變化現(xiàn)象;鄒正波等[11]計(jì)算了中國(guó)大陸及周緣地區(qū)在汶川地震前后兩年的GRACE衛(wèi)星重力變化,得出衛(wèi)星重力變化能揭示汶川地震對(duì)中國(guó)大陸的影響這一結(jié)論;徐海軍等[12]利用GRACE衛(wèi)星數(shù)據(jù)對(duì)印度MS6.8地震區(qū)域重力場(chǎng)變化特征進(jìn)行研究,證實(shí)強(qiáng)震就發(fā)生在重力變化高梯度帶內(nèi);姜永濤等[13]計(jì)算了青藏高原東緣多個(gè)時(shí)間尺度衛(wèi)星重力場(chǎng)變化,分析了3次大于MS7大震與衛(wèi)星重力場(chǎng)變化特征的關(guān)系;對(duì)于日本MW9.0地震,幾位學(xué)者也分別利用GRACE數(shù)據(jù)做了不同方面的研究[14-17]。

本文利用GRACE衛(wèi)星最新重力場(chǎng)模型,計(jì)算并分析云南地區(qū)及鄰域長(zhǎng)期的衛(wèi)星重力場(chǎng)變化特征,探討區(qū)域衛(wèi)星重力變化與區(qū)域大震的關(guān)系。

1區(qū)域概況

云南地區(qū)及鄰域處于印度板塊與歐亞板塊碰撞帶的東緣,地殼運(yùn)動(dòng)劇烈,地震活動(dòng)水平高且與活動(dòng)塊體關(guān)系密切[18]。該地區(qū)在20世紀(jì)共發(fā)生5級(jí)以上地震333次,其中7級(jí)以上大震13次,僅20世紀(jì)70年代以來(lái)發(fā)生的7級(jí)以上地震就多達(dá)8次,是研究地震活動(dòng)性的理想場(chǎng)所[19-20]。云南區(qū)域現(xiàn)今的斷裂構(gòu)造格局由青藏高原的東向擠出環(huán)境所奠定,在這種地殼運(yùn)動(dòng)背景下[21-22]形成了云南地區(qū)特殊的地形地貌。

圖1為云南地區(qū)地形與現(xiàn)今活躍的主要斷裂帶示意圖[23]。圖1標(biāo)注了2003—2014年該地區(qū)發(fā)生的MS5以上地震(http://earthquake.usgs.gov/),這些地震主要發(fā)生在在三塊區(qū)域,即怒江斷裂南段與龍陵—瀾滄斷裂交匯處的西北面;麗江—小金河斷裂帶、紅河斷裂、則木河斷裂和小江斷裂所共同圍成的區(qū)域;大涼山斷裂帶東方向與馬邊斷裂西南方向的附近區(qū)域。

2重力場(chǎng)變化計(jì)算方法

2.1重力場(chǎng)變化計(jì)算原理

在橢球面上,利用GRACE重力場(chǎng)模型計(jì)算時(shí)變重力變化的表達(dá)式為:

(1)

圖1　云南地區(qū)地形(單位:m)Fig.1　Topography in Yunnan region and its adjacent areas (Unit:m)

2.2GRACE去相關(guān)性濾波

GRACE衛(wèi)星觀測(cè)的僅僅是地球的中長(zhǎng)波重力變化,高頻信號(hào)的精確性非常低,重力場(chǎng)模型中的高頻噪音信號(hào)很多,在時(shí)變重力場(chǎng)圖像上表現(xiàn)為南北條帶異常,因此在使用GRACE重力場(chǎng)模型前,須對(duì)位系數(shù)進(jìn)行去相關(guān)和平滑處理。

本文所使用的是ICGEM(http://icgem.gfz-potsdam.de/ICGEM/ICGEM.html)提供的已經(jīng)過(guò)DDK5處理的GFZ Release 05數(shù)據(jù),其有效性已被證實(shí)[24]。DDK5是德國(guó)地學(xué)中心Kusche研究出的一個(gè)卷積矩陣較為簡(jiǎn)單的去相關(guān)平滑方法[25]。該方法考慮在反演月重力場(chǎng)模型過(guò)程中,GRACE軌道數(shù)據(jù)和星間距離測(cè)量數(shù)據(jù)存在的信號(hào)協(xié)方差和誤差協(xié)方差矩陣,通過(guò)“懲罰加權(quán)”來(lái)建立月重力場(chǎng)模型的去相關(guān)濾波矩陣,但由于該矩陣的濾波系數(shù)太多,考慮去相關(guān)濾波矩陣系數(shù)的特點(diǎn)[26],最終簡(jiǎn)化得出了一個(gè)“次卷積”去相關(guān)濾波核。

3重力場(chǎng)計(jì)算及其特征

3.1區(qū)域自由空氣重力異常

圖2為利用最新高階重力場(chǎng)模型EIGEN-6C4[27]計(jì)算得出的云南地區(qū)自由空氣重力異常。通過(guò)圖2和圖1的對(duì)比可以發(fā)現(xiàn),區(qū)域地形與自由空氣重力異常有很好的對(duì)應(yīng)關(guān)系,且在斷裂帶附近多會(huì)出現(xiàn)重力變化異常現(xiàn)象,其中瀾滄江斷裂、怒江斷裂、紅河斷裂、安寧河斷裂、大涼山斷裂和馬邊斷裂的負(fù)重力異常現(xiàn)象尤為顯著。這種現(xiàn)象說(shuō)明云南區(qū)域斷層多位于重力變化異常區(qū),其原因一方面可能是斷裂帶及其毗鄰區(qū)地面起伏對(duì)重力異常結(jié)果產(chǎn)生影響,另一方面,斷層的阻隔作用使得殼幔物質(zhì)運(yùn)動(dòng)受阻,斷裂帶成為一個(gè)地殼密度異常帶。

3.2區(qū)域年差分重力變化

為了盡量避免季節(jié)性因素對(duì)區(qū)域差分重力場(chǎng)的影響,本文以最大限度選取每年的共有月重力場(chǎng)。具體在計(jì)算2003—2014年的差分重力場(chǎng)變化時(shí),采用每年2、4、5、7、11和12月份的GRACE平均重力場(chǎng)模型計(jì)算(其中2012年5月份的數(shù)據(jù)有缺失,用當(dāng)年6月份的數(shù)據(jù)代替;2014年的數(shù)據(jù)只有7個(gè)月份的,選取1、3、4、5、8和9月份),GRACE月重力場(chǎng)模型經(jīng)過(guò)去相關(guān)濾波DDK5處理。由于GRACE雙星系統(tǒng)對(duì)位系數(shù)C20項(xiàng)不敏感,因此在計(jì)算中用人衛(wèi)激光測(cè)距得到的C20項(xiàng)作為替代。計(jì)算結(jié)果見(jiàn)圖3。

圖2　云南地區(qū)自由空氣重力異常(單位:10-5 ms-2)Fig.2　Free air gravity anomalies in Yunnan region and its adjacent areas (Unit:10-5 ms-2)

圖3　區(qū)域年度差分時(shí)變重力場(chǎng)(單位:10-8 ms-2)Fig.3　Regional annual differentiated temporal gravity field (Unit:10-8 ms-2 )

4.3區(qū)域大震與重力變化關(guān)系

從圖3可以看出,保山次級(jí)塊體在2003—2004年呈重力正值變化,且形成了垂直于大盈江斷裂的重力變化梯度帶;2006年重力變化梯度帶由垂直變?yōu)槠叫杏诖笥瓟嗔?并以大盈江斷裂為中心形成了重力正負(fù)變化的高梯度帶,說(shuō)明此期間西北部殼幔物質(zhì)加速南移,使保山次級(jí)塊體南部的騰沖次級(jí)塊體在2007年出現(xiàn)重力正值變化的高值區(qū),2008年8月21日發(fā)生的盈江縣MS5.9地震可能是由這種物質(zhì)高速流動(dòng)而誘發(fā)的斷層破裂。2009—2011年南北殼幔物質(zhì)呈強(qiáng)烈交替運(yùn)動(dòng),在重力變化高梯度帶上發(fā)生了2011年3月10日盈江縣MS5.8地震。2013年保山次級(jí)塊體向東擠壓,保山次級(jí)塊體與滇中次級(jí)塊體交界處形成的重力變化高梯度上發(fā)生了2次震源錯(cuò)動(dòng)方式為正斷層類(lèi)型的MS5以上地震(2013年3月3日洱源縣MS5.4地震和2013年8月31日德欽縣、得榮縣交界處MS5.6地震)。

2003—2004年滇中次級(jí)塊體呈重力正值變化,而2005—2006年重力負(fù)值變化劇烈,2007年又呈重力正值變化,2008—2009年呈連續(xù)的重力負(fù)值變化。這種地殼間反復(fù)的擠壓與松弛運(yùn)動(dòng),最終觸發(fā)了處在重力變化高梯度帶南華—楚雄—建水?dāng)嗔焉系囊Π部h2009年7月9日MS5.7地震。

云南地區(qū)的東部位于滇東塊體的西緣,2004年滇東塊體向西擠壓,在重力正負(fù)變化的高梯度帶上發(fā)生了2004年8月10日魯?shù)镸S5.6地震。重力負(fù)值變化區(qū)在震后向西擴(kuò)張,這可能是地震使原本的阻隔作用減少,殼幔物質(zhì)加速西移,在這一過(guò)程中發(fā)生了2005年8月5日會(huì)澤MS5.0地震。2009—2011年滇中次級(jí)塊體與滇東次級(jí)塊體交界處的重力變化梯度帶一直垂直于小江斷裂的走向,2012年重力變化梯度帶由垂直變?yōu)槠叫?該重力變化的高梯度帶上發(fā)生了2012年9月7日昭通MS5.6地震。此區(qū)域的重力變化梯度帶于2014年再次平行于小江斷裂走向,2014年8月3日在小江斷裂帶上發(fā)生了魯?shù)镸S6.2大震,這種現(xiàn)象可能說(shuō)明臨震前形成的平行于斷裂帶的重力變化梯度帶有利于地震破裂的發(fā)生[28],2014年4月4日發(fā)生的邵通MS5.4地震也有類(lèi)似的重力變化現(xiàn)象。

4三個(gè)特征點(diǎn)的月重力變化序列

為精確了解云南地區(qū)地震高發(fā)區(qū)域的重力場(chǎng)變化情況,本文在這3個(gè)區(qū)域(怒江斷裂南段與龍陵—瀾滄斷裂交匯處的西北面;麗江—小金河斷裂帶、紅河斷裂、則木河斷裂和小江斷裂所共同圍成的區(qū)域;大涼山斷裂帶東方向與馬邊斷裂西南方向的附近區(qū)域)內(nèi)各選擇一個(gè)特征點(diǎn),如圖1中A、B、C所示。每個(gè)特征點(diǎn)的經(jīng)度和緯度分別為其所屬區(qū)域內(nèi)2003—2014年所發(fā)生的MS5以上地震震中位置的平均值,即A(97.9°E,24.8°N)、B(101.5°E,25.8°N)和C(103.8°E,27.6°N)(圖1)。

本文利用ICGEM發(fā)布的2003—2014年共131個(gè)月的GRACE月重力場(chǎng)模型,通過(guò)計(jì)算獲取3個(gè)特征點(diǎn)位跨度為12年的月重力變化時(shí)間序列,并將3個(gè)地震高發(fā)區(qū)2003—2014年發(fā)生的大于MS5地震標(biāo)在對(duì)應(yīng)的圖上(圖4)。分析發(fā)現(xiàn)3個(gè)特征點(diǎn)的重力異常變化存在一定差異。

(1) 特征點(diǎn)A(97.9°E,24.8°N)位于怒江斷裂南段與龍陵—瀾滄斷裂交匯處的西北面。該區(qū)域的盈江縣先后發(fā)生了2008年8月21日MS6.0地震、2011年3月10日MS5.5地震、2014年5月23日MS5.8地震和2014年5月30日MS5.9地震。觀察這4次地震發(fā)生月份A點(diǎn)的重力異常變化值,發(fā)現(xiàn)這些變化值都很小,趨近于0。其中2008年8月21盈江縣MS6.0地震發(fā)生前期重力變化速率較大,在地震前一個(gè)月重力變化達(dá)到增量約11×10-8ms-2,形成了一個(gè)小的峰值,并于地震后的幾個(gè)月又恢復(fù)相對(duì)穩(wěn)定的波動(dòng),這說(shuō)明地震使震前地殼運(yùn)動(dòng)所積累的能量得到了釋放,同時(shí)也印證了地震發(fā)生后的位場(chǎng)恢復(fù)理論[29],這與圖3中結(jié)果相對(duì)應(yīng)。

(2) 特征點(diǎn)B(101.5°E,25.8°N)位于麗江—小金河斷裂帶、紅河斷裂、則木河斷裂和小江斷裂所共同圍成的區(qū)域。該區(qū)域發(fā)生了2003年7月21日大姚MS6.0地震、2003年10月16日大姚MS5.6地震、2008年8月30日川滇交匯處會(huì)理縣MS5.9級(jí)地震、2008年8月31日會(huì)理縣MS5.5地震和2009年7月9日會(huì)理縣MS5.7地震。該地區(qū)發(fā)震月份的重力異常變化具有統(tǒng)一性,即發(fā)震月份的重力異常變化都為正且重力變化量都約為5×10-8ms-2,每一個(gè)月份都處在波峰的位置。其中2003年10月16日大姚MS5.6地震發(fā)生后的一個(gè)月負(fù)重力變化量較大,達(dá)到-9.3×10-8ms-2,這可能是震后其區(qū)域地殼松弛造成的。

(3) 特征點(diǎn)C(103.8°E,27.6°N)位于大涼山斷裂帶東方向與馬邊斷裂西南方向的附近區(qū)域。該區(qū)域發(fā)生了2004年8月10日昭通MS5.1地震、2012年9月7日彝良縣MS5.6與MS5.3地震、2014年4月4日永善縣MS5.4地震和2014年8月3日魯?shù)榭hMS6.2地震。該區(qū)域重力變化無(wú)明顯特征,原因可能是它緊挨川滇菱形塊體東側(cè),川滇菱形塊體東邊界是對(duì)于菱形塊體活動(dòng)最為敏感的條帶,重力場(chǎng)變化極為復(fù)雜,難以從一個(gè)點(diǎn)位分析出其特征。但是可以發(fā)現(xiàn)彝良縣MS5.6與MS5.3地震發(fā)生前一個(gè)月的重力場(chǎng)變化與2008年8月21盈江縣MS6.0地震發(fā)震前的重力變化相似。

圖4　三個(gè)特征點(diǎn)2003—2014年的月重力變化時(shí)間序列(單位:10-8 ms-2)Fig.4　Time series of monthly gravity changes at three points from 2003 to 2014 (Unit:10-8 ms-2 )

5結(jié)論

本文利用GRACE月重力場(chǎng)模型,計(jì)算并分析云南地區(qū)的衛(wèi)星重力場(chǎng)動(dòng)態(tài)演化過(guò)程,最后分析了衛(wèi)星重力變化與區(qū)域大震的關(guān)系,得出主要結(jié)論如下:

(1) 云南地區(qū)的區(qū)域地形與自由空氣重力異常有很好的對(duì)應(yīng)關(guān)系,且斷裂帶附近多出現(xiàn)負(fù)重力異?，F(xiàn)象,其中瀾滄江斷裂、怒江斷裂、紅河斷裂、安寧河斷裂、大涼山斷裂和馬邊斷裂的負(fù)重力異?，F(xiàn)象尤為顯著。

(2) 從云南地區(qū)2003—2014年的逐年差分時(shí)變重力場(chǎng)變化圖像發(fā)現(xiàn),地震多發(fā)生在區(qū)域重力變化的高梯度帶上,且出現(xiàn)了與地震相對(duì)應(yīng)的特征性衛(wèi)星重力變化,反映了與地震相關(guān)的殼幔物質(zhì)的擠壓與松弛運(yùn)動(dòng)。

(3) 3個(gè)地震高發(fā)區(qū)特征點(diǎn)上的重力變化特征與地震的對(duì)應(yīng)關(guān)系各不相同。A點(diǎn)的重力變化值在地震發(fā)生月份都處于低值;B點(diǎn)發(fā)震月份的重力異常變化特征具有統(tǒng)一性,即發(fā)震月份都為重力正變化且量值約為5×10-8ms-2,同時(shí)每一個(gè)發(fā)震月份都處在重力變化波峰位置;C點(diǎn)的重力變化無(wú)明顯特征,原因可能是這一區(qū)域緊挨塊體活動(dòng)最復(fù)雜的川滇菱形塊體東邊界,重力變化極為復(fù)雜,難以從一個(gè)點(diǎn)分析出其變化特征。

參考文獻(xiàn)(References)

[1]褚慶忠,武澤,邵先杰,等.地震重力異常形成機(jī)理及其在地震預(yù)報(bào)中的作用[J].地震工程學(xué)報(bào),2014,36(1):201-206.

CHU Qing-zhong,WU Ze,SHAO Xian-jie,et al.Mechanism of Gravity Anomaly before Earthquake and Its Role in Earthquake Prediction[J].China Earthquake Engineering Journal,2014,36(1):201-206.(in Chinese)

[2]馬麗,吳兵.我國(guó)中強(qiáng)震前后重力場(chǎng)的時(shí)空變化特征[J].西北地震學(xué)報(bào),1988,26(4):96-108.

MA Li,WU Bing.Characteristics of Spatial and Temporal Gravity Changes before and after Moderate Strong Earthquakes in China[J].Northwestern Seismological Journal,1988,26(4):96-108.(in Chinese)

[3]劉仲全.云南及鄰區(qū)MS≥6.0級(jí)地震的重力異常指標(biāo)體系研究[J].地震研究,1995,18(4):378-387.

LIU Zhong-quan.System Studies of Gravity Anomaly Critiria Before Earthquakes withMS≥6.0 in Yunnan and Its Adjacent Areas[J].Journal of Seismological Research,1995,18(4):378-387.(in Chinese)

[4]吳國(guó)華,羅增雄,賴(lài)群.麗江7.0級(jí)地震前后滇西實(shí)驗(yàn)場(chǎng)的重力異常變化特征[J].地震研究,1997,20(1):103-109.

WU Guo-hua,LUO Zeng-xiong,LAI Qun.The Variation Characteristics of Gravity Anomaly in the Earthquake Prediction Test Site in Western Yunnan before and after theMS7.0 Lijiang Earthquake[J].Journal of Seismological Research,1997,20(1):103-109.(in Chinese)

[5]付廣裕,張勇軍,劉少明,等.昆侖山口西8.1級(jí)大震前后重力場(chǎng)動(dòng)態(tài)變化[J].大地測(cè)量與地球動(dòng)力學(xué),2003,23(1):42-46.

FU Guang-yu,ZHANG Yong-jun,LIU Shao-ming,et al.Dynamic Gravity Changes before and after Earthquake ofMS8.1 Occurred in West to Kunlun Mountain Pass[J].Journal of Geodesy and Geodynamics,2003,23(1):42-46.(in Chinese)

[6]祝意青,李輝,朱桂芝,等.青藏塊體東北緣重力場(chǎng)演化與地震活動(dòng)[J].地震學(xué)報(bào),2004,26(增刊1):71-78,175.

ZHU Yi-qing,LI Hui,ZHU Gui-zhi,et al.Gravity Evolution and Earthquake Activities of the Northeastern Edge of Qinghai—Xizang Block[J].Acta Seismologica Sinica,2004,26(Supp 1):71-78,175.(in Chinese)

[7]Sun W,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].Journal of Geophysical Research:Solid Earth (1978—2012),2004:109(B4).

[8]Han S C,Shum C K,Bevis M,et al.Crustal Dilatation Observed by GRACE after the 2004 Sumatra—Andaman Earthquake[J].Science,2006,313(5787):658-662.

[9]周新,孫文科,付廣裕.重力衛(wèi)星GRACE檢測(cè)出2010年智利MW8.8地震的同震重力變化[J].地球物理學(xué)報(bào),2011,54(7):1745-1749.

ZHOU Xin,SUN Wen-ke,FU Guang-yu.Gravity Satellite GRACE Detects Coseismic Gravity Changes Caused by 2010 ChileMW8.8 Earthquake[J].Chinese Journal of Geophysics,2011,54(7):1745-1749.(in Chinese)

[10]段虎榮,張永志,劉鋒,等.利用GRACE衛(wèi)星數(shù)據(jù)研究汶川地震前后重力場(chǎng)的變化[J].地震研究,2009,32(3):295-298,322.

DUAN Hu-rong,ZHANG Yong-zhi,LIU Feng,et al.Study of the Gravity Variance in Chinese Mainland before and after Wenchuan Earthquake with GRACE Gravity Models[J].Journal of Seismological Research,2009,32(3):295-298,322.(in Chinese)

[11]鄒正波,李輝,康開(kāi)軒,等.汶川地震與衛(wèi)星重力變化[J].大地測(cè)量與地球動(dòng)力學(xué),2013,33(1):5-7.

ZOU Zheng-bo,LI Hui,KANG Kai-xuan,et al.Wenchuan Earthquake and Satellite Gravity Variation[J].Journal of Geodesy and Geodynamics,2013,33(1):5-7.(in Chinese)

[12]徐海軍,張永志,段虎榮.印度MS6.8地震區(qū)域重力場(chǎng)變化特征研究[J].大地測(cè)量與地球動(dòng)力學(xué),2012,32(1):10-12,16.

XU Hai-jun,ZHANG Yong-zhi,DUAN Hu-rong.Gravity Variation Detected by GRACE in IndianMS6.8 Earthquake Region[J].Journal of Geodesy and Geodynamics,2012,32(1):10-12,16.(in Chinese)

[13]姜永濤,張永志,王帥,等.基于衛(wèi)星重力的青藏高原東緣多尺度時(shí)變重力場(chǎng)研究[J].地震工程學(xué)報(bào),2014,36(2):353-359.

JIANG Yong-tao,ZHANG Yong-zhi,WANG Shuai,et al.Study on Multi-scale Time-varying Gravity Filed in the Eastern Margin of the Tibetan Plateau Based on Satellite Gravity[J].China Earthquake Engineer Journal,2014,36(2):353-359.(in Chinese)

[14]鄒正波,羅志才,吳海波,等.日本MW9.0地震前GRACE衛(wèi)星重力變化[J].測(cè)繪學(xué)報(bào),2012,41(2):171-176.

ZOU Zheng-bo,LUO Zhi-cai,WU Hai-bo,et al.Gravity Changes Observed by GRACE before the JapanMW9.0 Earthquake[J].Acta Geodaetica et Cartographica Sinica,2012,41(2):171-176.(in Chinese)

[15]張永志,夏朝龍,王衛(wèi)東,等.日本9.0級(jí)地震區(qū)重力梯度的時(shí)空分布[J].大地測(cè)量與地球動(dòng)力學(xué),2013,33(6):1-4.

ZHANG Yong-zhi,XIA Chao-long,WANG Wei-dong,et al.Temporal and Space Contribution of Gravity Gradients in JapanMW9.0 Earthquake Area[J].Journal of Geodesy and Geodynamics,2013,33(6):1-4.(in Chinese)

[16]王武星,顧國(guó)華,陳石.利用 GRACE 觀測(cè)資料分析日本MW9.0地震前區(qū)域重力變化特征[J].地震地質(zhì),2014,36(2):523-535.

WANG Wu-xing,GU Guo-hua,CHEN Shi.Study on Regional Gravity Changes before theMW9.0 Japan Earthquake Detected by GRACE[J].Seismology and Geology,2014,36(2):523-535.(in Chinese)

[17]邢樂(lè)林,李輝,玄松柏,等.GRACE衛(wèi)星觀測(cè)到的日本9.0級(jí)大地震重力前兆信息[J].大地測(cè)量與地球動(dòng)力學(xué),2011,31(2):1-3.

XING Le-lin,LI Hui,XUAN Song-bai,et al.Gravity Caused byMS9.0 Strong Eathquake in Japan Precursor Detected by GRACE[J].Journal of Geodesy and Geodynamics,2011,31(2):1-3.(in Chinese)

[18]皇甫崗,秦嘉政,李忠華,等.云南地震類(lèi)型分區(qū)特征研究[J].地震學(xué)報(bào),2007,29(2):142-150,229.

HUANGFU Gang,QIN Jia-zheng,LI Zhong-hua,et al.Subarea Characteristics of Earthquake Types in Yunnan[J].Acta Seismologica Sinica,2007,29(2):142-150,229.(in Chinese)

[19]皇甫崗,秦嘉政.云南地區(qū)大震活動(dòng)規(guī)律研究[J].地震地質(zhì),2006,28(1):37-47.

HUANGFU Gang,QIN Jia-zheng.Study of the Seismicity of Strong Earthquakes in Yunnan Area[J].Seismology and Geology,2006,28(1):37-47.(in Chinese)

[20]劉漢雄,張志強(qiáng).云南地震活動(dòng)的特征[J].西北地震學(xué)報(bào),1995,17(1):43-46.

LIU Han-xiong,ZHANG Zhi-qiang.Characteristics of Seismicity in Yunnan Province[J].Northwestern Seismological Journal,1995,17(1):43-46.(in Chinese)

[21]馬宗晉,張家聲,汪一鵬.青藏高原三維變形運(yùn)動(dòng)學(xué)的時(shí)段劃分和新構(gòu)造分區(qū)[J].地質(zhì)學(xué)報(bào),1998,72(3):211-227.

MA Zong-jin,ZHANG Jia-sheng,WANG Yi-peng.The 3-d Deformational Movement Episodes and Neotectonic Domains in the Qinghai—Tibet Plateau[J].Acta Geologica Sinica,1998,72(3):211-227.(in Chinese)

[22]張培震,沈正康,王敏,等.青藏高原及周邊現(xiàn)今構(gòu)造變形的運(yùn)動(dòng)學(xué)[J].地震地質(zhì),2004,26(3):367-377.

ZHANG Pei-zhen,SHEN Zheng-kang,WANG Min,et al.Continuous Deformation of the Tibetan Plateau from Global Positioning System Data[J].Seismology and Geology,2004,26(3):367-377.(in Chinese)

[23]鄧起東.中國(guó)活動(dòng)構(gòu)造圖[M].北京:地震出版社,2007.

DENG Qi-dong.Active Tectonics Map of China[M].Beijing:Seismological Prese,2007.(in Chinese)

[24]姜永濤,張永志,王帥,等.GRACE時(shí)變重力場(chǎng)濾波方法研究[J].測(cè)繪通報(bào),2014(11):1.

JIANG Yong-tao,ZHANG Yong-zhi,WANG Shuai,et al.Study on Filtering Methods of GRACE Temporal Gravity Variation[J].Bulletin of Surveying and Mapping,2014(11):1.(in Chinese)

[25]Kusche J,Schmidt R,Petrovic S,et al.Decorrelated GRACE Time-variable Gravity Solutions by GFZ,and Their Validation Using a Hydrological Model[J].Journal of Geodesy,2009,83(10):903-913.

[26]Swenson S,Wahr J.Post-processing Removal of Correlated Errors in GRACE Data[J].Geophysical Research Letters,2006,33(8):L08402.

[27]F?rste C,Bruinsma S,Abrikosov O,et al.EIGEN-6C4-The Latest Combined Global Gravity Field Model Including GOCE Data up to Degree and Order 1949 of GFZ Potsdam and GRGS Toulouse[C].EGU General Assembly Conference Abstracts.2014,16:3707.

[28]申重陽(yáng),李輝,孫少安,等.2008年于田MS7.3地震前重力場(chǎng)動(dòng)態(tài)變化特征分析[J].大地測(cè)量與地球動(dòng)力學(xué),2010,30(4):1-7.

SHEN Chong-yang,LI Hui,SUN Shao-an,et al.Characteristic Analysis of Dynamic Gravity Change before YutianMS7.3 Earthquake 2008[J].Journal of Geodesy and Geodynamics,2010,30(4):1-7.(in Chinese)

[29]祝意青,梁偉鋒,湛飛并,等.中國(guó)大陸重力場(chǎng)動(dòng)態(tài)變化研究[J].地球物理學(xué)報(bào),2012,55(3):804-813.

ZHU Yi-qing,LIANG Wei-feng,ZHAN Fei-bing,et al.Study on Dynamic Change of Gravity Field in China Continent[J].Chinese Journal Of Geophysics,2012,55(3):804-813.(in Chinese)

Satellite Gravity Changes in Yunnan and Adjacent Regions and Its Relation with Regional Strong Earthquakes

JIAO Jia-shuang1, ZHANG Yong-zhi1, ZHANG Kai-nan1, JIANG Yong-tao2, HUAI Yan-ke1

(1.SchoolofGeologyEngineeringandGeomatics,Chang’anUniversity,Xi’an710000,Shaanxi,China;2.NanyangNormalUniversity,Nanyang473061,Henan,China)

Abstract:Yunnan is one of the main seismic activity areas in China.It lies close to the subduction zone where the Indian plate moves towards the Eurasian plate and belongs to the southern section of the north-south seismic belt in China.Historical earthquake activity in this area shows that earthquakes of high frequency and large magnitude occur.The Yunnan region and its adjacent areas are ideal locations for the study of seismicity.Existing research shows that regional seismic activity is closely related to gravity field changes that are mainly caused by geological activities such as crustal deformation and fault material change.Such changes reflect crustal movement and the associated processes involved in large earthquakes.Therefore,gravity changes can be used to predict earthquakes.The GRACE (Gravity Recovery and Climate Experiment) binary star system,jointly developed by NASA and the DLR,provides a time-varying gravity field and thus important information for the exploration of Earth's interior materials.GRACE satellite observations clearly reflect the motion of the Earth’s plates and expand the scope of gravity anomaly seismic observations that can be used for long-term forecasts of earthquakes.In this paper,we use gravity field data from GRACE satellites to calculate and analyze the satellite gravity field variation characteristics from 2003 to 2014 in the Yunnan region and its adjacent areas.We also obtained the monthly gravity change sequence of three feature points,and in this study we analyze their change characteristics by combining the distribution of M≥5 earthquakes from 2003 to 2014 in the study area.Results show that the MS≥5 earthquakes in the Yunnan region and adjacent areas are reflected in satellite images of gravity field changes,and that large earthquakes often occur on the high gradient zone of regional gravity changes.The earthquake zone often shows a corresponding change in satellite gravity characteristics both before and after large earthquakes,which may reflect crust-mantle material adjustment.The corresponding earthquake gravity change characteristics of the three points are different,which may indicate that research into the gravity field changes in various regions and faults is significant for earthquake prediction.

Key words:Yunnan region and its adjacent areas; GRACE; satellite gravity change; earthquake; fault

DOI:10.3969/j.issn.1000-0844.2016.02.0260

中圖分類(lèi)號(hào):P315.72+6

文獻(xiàn)標(biāo)志碼:A

文章編號(hào):1000-0844(2016)02-0260-09

作者簡(jiǎn)介:焦佳爽(1990-), 男,在讀碩士生,主要從事衛(wèi)星重力研究。E-mail:jiaojiashuang@foxmail.com。

基金項(xiàng)目:國(guó)家自然科學(xué)基金項(xiàng)目(41374028,41274083,41304013);陜西省教育廳自然科學(xué)科研項(xiàng)目(12JK0798)

收稿日期:①2015-05-11