徐東卓 尹海權 朱傳寶 孫非非
1)中國地震局第一監測中心,天津市河東區耐火路7號 300180
2)青海省第三地質礦產勘查院,西寧 810029
2017年8月8日,四川省阿壩州九寨溝縣發生 7.0級地震,震中 33.20°N、103.82°E,震源深度約為20km。九寨溝地震是繼南北地震帶中南段和巴顏喀拉地塊周緣多次強震后發生的又一次顯著地震。此次地震震中位于巴顏喀拉地塊的東北部、東昆侖斷裂帶東段塔藏斷裂附近,根據余震精定位發現,余震條帶表現為紡錘形,長軸走向約為340°~350°,與全新世的塔藏斷裂小角度相交,發震斷層可能是東昆侖斷裂東端的一條次級活動斷層。巴顏喀拉地塊是青藏高原內部受印度板塊俯沖推擠亞歐板塊而形變強烈的活動塊體(陳長云等,2013;徐東卓等,2017),其北邊界是以左旋走滑為主的東昆侖斷裂帶,該斷裂帶曾發生過2001年昆侖山口西8.1級地震、1997年瑪尼7.5級地震以及此次九寨溝7.0級地震等強震;東邊界是以擠壓特征為主的龍門山斷裂帶,由3條近乎平行的主斷裂組成,近年曾發生過2008年汶川8.0級、2013年蘆山7.0級等強震;南邊界是以左旋走滑為主的玉樹-鮮水河斷裂帶,曾發生過1981年道孚6.9級、2010年玉樹7.1級等中強震;西邊界是左旋走滑的阿爾金斷裂帶,歷史強震活動較其他邊界帶偏弱,但也發生過包括2008年于田7.3級地震在內的多次顯著地震。據震例統計結果,我國大陸全部8級地震、80%以上的7級地震都發生在塊體邊界區(張國民等,2000),巴顏喀拉地塊更是如此。
有研究認為,汶川地震后,庫侖應力在龍門山斷裂帶南段加速積累,深部孕震環境形成,進而導致了蘆山地震的發生;1997年以來的巴顏喀拉地塊周緣強震的粘彈性觸發研究表明,東昆侖的瑪尼7.5級地震與昆侖山口西8.1級地震、昆侖山口西8.1級地震與玉樹7.1級地震以及昆侖山口西8.1級地震與汶川8.0級地震之間存在不同程度的觸發關系(沈正康等,2003;萬永革等,2007;汪建軍等,2009;程佳等,2011),萬永革等(2007)分析了 1920年海原8.5級地震以來青藏高原東北部7級以上地震的庫侖應力積累過程及其對后續地震的影響,認為85%的后續地震受到了之前地震的觸發作用,并對青藏高原東北部強震之間的觸發關系進行了論述。關于昆侖山口西8.1級地震與汶川8.0級地震之間的聯系,有研究認為,巴顏喀拉地塊作為青藏高原向東擠出的主要地區,具有整體向東運動的趨勢,2次大地震之間存在一定的應力轉換作用(Burchfiel et al,2008;Zhang et al,2009;Luo et al,2010)。那么,巴顏喀拉地塊邊界斷裂帶上近年來的強震與九寨溝7.0級地震之間是否存在某種聯系?九寨溝地震發生前,該區域諸多的形變資料表現出了異常特征,這種異常表現與地震孕育之間的深層次關系如何?都是值得探究的問題。本文利用震中區周緣、巴顏喀拉地塊周緣的跨斷層觀測、大面積水準測量等資料,聯合已有的GNSS分析結果,借助多種處理分析方法,分析九寨溝地震前后的區域地殼形變演化特征,探討深部孕震構造環境的變化機理,并對未來地震危險區域給出一定推測。
為了研究青藏高原東北緣及巴顏喀拉地塊整體地殼形變與九寨溝7.0級地震間的相關性,本文選取了以上區域重要斷裂帶和塊體邊界斷裂上的跨斷層流動場地,其主要位于西秦嶺構造區、六盤山斷裂帶、龍門山斷裂帶和鮮水河斷裂帶,其中,距震中最近的為西秦嶺構造區,場地也多在距震中100km以外。西秦嶺構造區和六盤山斷裂帶的場地始建于20世紀70、80年代,其中,西秦嶺構造區共有20余處,六盤山斷裂帶有10余處,均積累了至少30年的觀測數據,資料完整而連續,監測手段主要為短水準觀測。短水準測量周期每年3~5次,震情緊迫的情況下,還會進行加密觀測。龍門山、鮮水河斷裂帶的跨斷層場地共有16處。其中,侏倭、格簍壩子、虛墟、溝普、龍燈壩、老乾寧和折多塘等7處為基線和水準綜合觀測場地,其余均為水準觀測。它們多始建于20世紀80年代,近年也存在新增建場地,觀測頻次基本為每月1期,震情緊迫時會加密至每半個月1期,個別年份每2個月觀測1期,數據連續且完整。監測儀器運行正常,各項技術指標符合規范要求,觀測結果精度優良。
本文使用的大面積區域水準資料線路西起四川阿壩,東至重慶,北起甘肅臨潭,南到四川西昌,覆蓋了包括九寨溝7.0級地震震中區在內的四川中西部地區(圖1)。水準資料主要利用2期一等復測資料,第1期施測時間為1991~1998年,第2期為2012~2013年。使用的光學水準儀、電子水準儀以及配套的銦瓦水準標尺和尺臺檢驗項目齊全,各項限差均符合規范要求,施測過程嚴格按照《國家一、二等水準測量規范》(中國國家標準化管理委員會,2006)執行。

圖1 研究區跨斷層場地及水準路線分布
本文將收集到的跨斷層形變觀測數據進行了預處理。首先,消除因干擾引起的觀測異常突跳變化,對觀測資料進行連續性檢驗,剔除重復的數據。然后,對日常觀測中因各種原因不能施測而造成的數據缺失采用線性插值法修補,得到統一觀測周期的可靠的觀測數據。對于投入觀測不久的新場地,由于數據量較少,場地穩定性有待檢驗,暫時不作分析。
對先期處理好的數據運用MATLAB進行編程,將其原始曲線繪出。通過對曲線進行分析發現,九寨溝地震前后,西秦嶺、六盤山、龍門山和鮮水河等斷裂帶上的跨斷層流動場地出現了明顯異常變化。其中,西秦嶺構造區上表現出異常的場地有殿溝、盤古川、柳家溝、毛羽溝、巴沙溝、武都東、黃家壩和硤口驛等;六盤山斷裂帶上有三關口、安國鎮、甘溝窯、和尚鋪和六盤山等;龍門山斷裂帶上有寶興和雁門等;鮮水河斷裂帶上有虛墟、老乾寧、格樓壩子、折多塘、恰叫和榆林等(圖 2、3)。
2.1.1 西秦嶺構造區

圖2 西秦嶺、六盤山構造區跨斷層場地異常曲線

圖3 龍門山、鮮水河斷裂帶跨斷層場地異常曲線
殿溝5-2水準數據積累的時間較短,但是九寨溝地震震后的最新數據顯示其大幅下降,變化量為0.41mm,該場地距九寨溝地震震中相對較近,震前異常不明顯,震后出現了明顯異常,應為地震導致場地所跨斷層發生活化調整所致。盤古川12-10水準是新增測段,資料積累時間尚短,震前水準出現大幅上升,幅度達到1.89mm,震后曲線轉折向下,最新數據顯示異常持續,應為受同震作用的影響,短期內斷層活動發生轉折。柳家溝5-7水準對包括汶川地震在內的幾次中強震有不同程度的反應,2015年11月至今已累計大幅下降2.57mm,斷層震前表現為強烈擠壓運動,震后依然下降,這類似于永登5.8級地震后的趨勢,即5年后再次發生了景泰5.9級地震,這值得持續關注。毛羽溝3-2水準曲線顯示,2008年以后斷層活動以逆斷活動為主,2013年岷縣漳縣6.6級地震前3~4個月斷層擠壓運動加速異常,震后出現趨勢性轉折,轉為拉張上升;九寨溝地震震前,由之前的持續振蕩上升轉變為快速下降,下降幅度0.92mm,震后曲線小幅恢復上升。巴沙溝3-2水準2014年以來改變以往趨勢性轉折上升,2015年底至今在高位波動。武都東2-3水準2016年11月由之前的平穩波動,轉而大幅下降,變化量為1.15mm,并在2017年3月轉折上升0.49mm,震后數據暫無明顯變化。黃家壩2-3水準在2003年岷縣5.2級、2004年岷縣5.0級地震前表現出一定的突跳變化,且九寨溝地震前亦出現波動異常。硤口驛1-7水準在九寨溝地震前表現出明顯的快速拉升變化。
距九寨溝地震震中區最近的西秦嶺構造區多表現出震前中短期壓性運動和異常波動、震后快速恢復的活動特點,與地震震源機制解給出的擠壓特征相符,震后隨著區域應力的釋放,斷層活動得到迅速恢復。個別場地的另類異常活動表現,也印證了西秦嶺地區構造體系和地殼形變的復雜性。西秦嶺構造區因青藏高原東北緣物質向東擠出運動,遇到鄂爾多斯地塊和華南地塊的阻礙,在三方互相作用下而表現出構造活動的復雜性,在中強震前,該區域的跨斷層場地多表現出群體性壓性和異常波動的特點。
2.1.2 六盤山斷裂帶
三關口1-3水準在2000年景泰5.9級地震后、2008年汶川8.0級地震前、2013年岷縣漳縣6.6級地震前后均出現轉折突跳、尖點突跳等異常變化,2015年底出現快速下降異常,到九寨溝地震前已累計下降7.66mm,震后繼續下降。安國2-4、1-4水準在汶川8.0級、岷縣漳縣6.6級地震前近似同步地出現過顯著異常變化,九寨溝7.0級地震前也表現出大幅振蕩上升,多次地震前斷層均表現為壓性轉為張性的變化特點。甘溝窯7-5水準在2000年景泰5.9級、2003年民樂-山丹6.1級地震前出現破趨勢的尖點突跳異常,2016年開始以高于以往的速率快速拉升,加速突跳,持續至今。和尚鋪26-118水準較好地回應了幾次中強震,主要表現為尖點突跳異常,整體表現為緩慢抬升的長趨勢性異常,斷層活動表現為拉張特征。六盤山1-4水準在幾次5級以上地震時均出現尖點突跳異常,2010年以來打破以往運動趨勢,表現出類似周期性的異常特征,斷層活動波動強烈。
六盤山斷裂帶受到青藏高原和鄂爾多斯地塊的共同作用,地質背景上表現為壓性運動特征,在九寨溝地震前,長趨勢上整體表現出以壓性轉張性為主的運動方式。造成該異常變化的原因應是以巴顏喀拉地塊的整體性運動為代表的青藏高原東北緣應力場發生變化,造成了其與鄂爾多斯地塊共同作用的六盤山斷裂帶發生應力調整。三關口場地所處的六盤山次級斷裂,九寨溝地震后依然表現出強烈的壓性特征,應與該次級斷裂的特殊性有關。六盤山斷裂帶上的多數場地九寨溝地震后異常仍然未明顯恢復,顯示該區域的應力積累進一步增強,九寨溝地震的發生并未緩解該斷裂帶的壓力,應持續關注區域地震危險性。震前的異常變化特征包含了九寨溝地震的前兆異常,震后卻未恢復,從長期來看,也可能是下一次地震構造事件的前兆。
2.1.3 龍門山斷裂帶
自觀測以來,寶興X1-B1水準曲線表現為下降趨勢,斷層活動整體表現為持續擠壓逆沖特征,九寨溝地震前出現小幅加速下降,隨后轉折上升。2017年以來雁門D-E水準打破以往的周期性振蕩下降變化,表現為波動幅度加大和突跳,這應該與九寨溝地震有關。
九寨溝地震前后,龍門山斷裂帶整體上異常表現相對其他構造帶少,表現為繼承性的逆沖擠壓運動特征。
2.1.4 鮮水河斷裂帶
虛墟A-B基線在蘆山地震時出現大幅變化,2014年以來繼續振蕩上升,2016年下半年到九寨溝地震前一直在高位異常波動。虛墟A-B水準整體表現出持續下降的長趨勢異常,九寨溝地震前出現1.08mm的上升波動,震后略有恢復。虛墟段斷層活動表現為強烈的左旋走滑特征。格簍壩子A1-C基線2015年8月快速拉升,斷層左旋運動加快,到2016年初累計上升幅度達4.27mm,之后直到九寨溝地震前持續高值波動,震后略有恢復。格簍壩子A-C水準震前表現為振蕩緩慢上升的長趨勢異常。斷層整體活動以左旋運動為主兼有拉張特征。恰叫1-2基線2008年至今曲線周期性持續上升,周期幅度明顯加大,2017年有轉折下降跡象。老乾寧5-3、1-3基線在汶川地震后同步改變,趨勢轉為上升,蘆山地震時發生大幅異常變化,之后持續振蕩上升,九寨溝地震時發生快速下降到上升的異常變化。老乾寧5-3、1-3水準同步表現出持續下降的長趨勢變化。斷層活動整體表現為右旋壓性特征。折多塘B-D基線持續振蕩上升,九寨溝地震前快速拉升1mm,震后又快速轉折下降1.25mm。榆林B-C水準2017年初到九寨溝地震前出現大幅上升,震后表現略有恢復,震前斷層活動表現為拉張運動。龍燈壩A-C基線在2013年蘆山地震時發生大幅異常變化,之后周期性振蕩上升,抬升速率加快,持續至今。侏倭A-B基線表現為持續上升趨勢,左旋運動明顯。侏倭A-B水準同樣表現為長趨勢上升變化。斷層活動表現為左旋運動兼具拉張特點。
由于川滇地塊與巴顏喀拉地塊的差異性運動,鮮水河斷裂帶全新世以來以左旋走滑運動為主,兼具擠壓性質。九寨溝地震前,鮮水河斷裂帶表現為左旋運動加速,同時表現為拉張特點。王敏等(2008)、劉冠中等(2013)通過對鮮水河斷裂帶斷層活動進行地表形變監測發現,鮮水河斷裂帶以北的爐霍段、道孚段的走滑量較大,表現為張性的左行走滑,與本文的分析結果一致,而以南的康定段、石棉段走滑量小,運動形態復雜,這可能與地表的形變和深部物質的流動相關(張培震,2008)。震后個別場地未明顯恢復,與場地所處斷裂帶的不同段落和次級斷層的特性有關,也可能代表大區域上的構造應力在繼續起作用,這種應力可能未因九寨溝地震而立即發生轉變。
對于鮮水河斷裂上數據連續完整、同時有2條同樁基線和水準的場地,采用以下3個公式計算斷層現今的三維活動量,即水平走滑量、水平張壓量和垂直升降量


其中,ΔH為水準垂直變化量(以上盤相對于下盤的下降為正);ΔL為基線伸縮量(以伸長為正);ΔS為斷層活動水平走滑量(兩盤做順扭運動為正);ΔR為斷層活動水平張壓量;α1和α2分別為由斷層線方向逆時針轉動至與2條測線方向重合或者平行時的角度。當ΔS為正值時,表示右旋活動;ΔR為正值時,表示張性活動;若兩者均為負值,則表示斷層活動性質相反(薄萬舉等,1998)。
根據計算所得的鮮水河斷裂帶不同段落的三維運動參數,繪制出三維活動參數時間序列(圖4)。由圖4可見,水平走滑量的計算結果顯示斷層在觀測期內為左旋活動,且累積活動量隨時間呈線性增加,2013年7.0級蘆山地震時發生了大幅異常變化。垂直升降量結果顯示,虛墟場地線性特征明顯,表現出以穩定斜率緩慢下降、斷層表現出持續穩定的活動特點;老乾寧場地的垂直升降量整體變化幅度不大,斷層活動表現出弱張性。水平張壓量時序曲線顯示,虛墟和老乾寧場地均表現出近期大幅持續拉升的態勢,其中,虛墟場地在昆侖山口西8.1級地震后出現破趨勢的異常,之后開始由下降轉為上升,上升期間分別對2008年汶川8.0級地震和2013年蘆山7.0級地震作出映震反應,蘆山地震后開始大幅上升,斷層活動表現為強烈拉張運動。老乾寧場地在汶川地震后出現破趨勢的上升變化,在蘆山地震時大幅下降之后轉折大幅拉升,持續至今。跨斷層場地的三維活動參數時序結果顯示,鮮水河斷裂在九寨溝地震前,表現出強烈的左旋走滑兼具拉張的運動特征。基于以上,小結如下。

圖4 鮮水河斷裂帶跨斷層場地三維活動量
(1)由巴顏喀拉地塊周緣以及距震中區域相對較近的跨斷層流動場地原始觀測曲線可見,西秦嶺構造區西段、六盤山斷裂帶、龍門山斷裂帶和鮮水河斷裂帶的跨斷層場地均不同程度地出現了異常表現。西秦嶺構造區在九寨溝地震前多表現出中短期壓性運動、震后快速恢復的活動特點,同時,多個場地震前表現出斷層活動無序混亂、波動較強的復雜性。六盤山斷裂帶受青藏高原和鄂爾多斯地塊的共同作用,地質背景上表現為壓性運動特征,在九寨溝地震前,整體表現出以壓性轉張性為主的運動學方式。由震例分析可知,區域中強震前西秦嶺和六盤山地區距震中300~400km內,甚至更遠地區,會出現跨斷層異常變化。三關口場地所跨斷層在地震前后表現出持續擠壓運動,可能與區域內局部應力增強有關。九寨溝地震前,龍門山斷裂帶表現為明顯的逆沖擠壓運動特征。九寨溝地震前,鮮水河斷裂帶表現為左旋運動加速,同時表現為拉張特點,尤其在蘆山地震后這一特征更為明顯,可見2次7.0級地震之間存在一定相關性。巴顏喀拉地塊整體性的異常運動,造成了九寨溝地震,同時也是鮮水河斷裂震前震后出現異常變化的原因。
(2)由跨斷層場地原始觀測曲線和三維活動參數時序曲線可見,九寨溝地震前,巴顏喀拉地塊周緣的龍門山斷裂帶和鮮水河斷裂帶均出現異常變化,因此,九寨溝地震與巴顏喀拉地塊整體性SE向旋轉運動相關。通過鮮水河斷裂帶北西段表現為強烈左旋運動以及阿爾金斷裂帶和東昆侖斷裂帶的左旋運動可知,應力可能隨著汶川地震、蘆山地震、玉樹地震、于田地震、昆侖山口西地震轉移到了塊體東北部的東昆侖斷裂帶東端區域,且受到NE向的以逆沖擠壓運動為主的龍門山斷裂帶的阻擋影響,應力在該區域不斷集中積累,孕震環境形成。同時,應力沿著巴顏喀拉地塊邊界斷裂順時針運移,鮮水河斷裂帶南東段表現出了明顯的拉張伸展的運動方式。受包括巴顏喀拉地塊在內的高原內部構造應力變化的影響,青藏高原東北緣也存在應力調整,這直接影響到了西秦嶺構造區和六盤山斷裂帶的運動學方式。
2期一等水準資料預處理后,對其進行平差數據處理。主要采用分段動態線性速率模型,平差計算得到該時間段內的地殼垂直形變速率,選取穩定的基巖水準測點作為基準,平差結果顯示,單位權中誤差小于1.0mm,效果較好。平差后對結果進行預處理,刪除個別量級過大或與相鄰點差異過大的突變點。采用多面函數法對結果進行格網化,并計算每個格網點的梯度模值,繪制垂直形變速率等值線圖與梯度圖(孫啟凱等,2017)。
由計算得出的研究區域垂直運動速率結果可知,區域內存在2處明顯地殼形變表現為隆升的區域,其中,北部的甘東南地區隆升速率大于西部的巴顏喀拉地塊東部地區,四川盆地整體表現出相對穩定的下沉特征,九寨溝地震震中區域地殼形變不明顯(圖5)。
(1)巴顏喀拉地塊東部地區存在大范圍差異性隆升表現。康定附近鮮水河斷裂帶和龍門山斷裂帶交匯周緣的地殼垂直隆升速率為2~4mm/a,與王敏(2009)給出的1.8~3.6mm/a的結果基本一致;馬爾康到紅原附近區域隆升速率為2~4mm/a,并分別在這2處區域形成形變梯度帶,而其他地區為1~2mm/a。來自印度板塊向北持續強有力的推擠作用,造成青藏高原部分物質向東運移,遇到NE向橫亙的龍門山斷裂帶的阻擋,巴顏喀拉地塊東部物質不斷積累,造成了地殼隆升,同時,地塊內部構造的差異性和不均衡性,造成了隆升速率的差異。鮮水河斷裂帶與龍門山斷裂帶交匯區周圍,高原內部物質遇到相對穩定的四川盆地阻擋后,隨著巴顏喀拉地塊SE向順時針旋轉流動至該區域集聚,造成隆升。
(2)青藏高原東北部的甘東南地區表現為隆升特征。在舟曲到宕昌附近存在2~4mm/a的地殼垂直隆升,與地質背景一致。高原物質向多個方向流動,其中,NE向的物質運移至該區域受到深大斷裂和剛硬體的攔阻,導致累積隆升。
(3)四川盆地表現為相對穩定的下沉特征。整體來看,區域沉降速率并不太顯著,為0~-1mm/a,成都附近稍高的沉降速率可能與城市過多抽取地下水有關。總體特征與地質背景一致,盆地處于揚子板塊中表現穩定剛硬的古陸核上,相對于逆沖其上的巴顏喀拉地塊東部,其表現為下沉特點。
(4)九寨溝地震震中區附近表現為相對穩定的弱應變特征。九寨溝到松潘附近的形變速率為0~1mm/a,地殼形變積累不顯著。GPS結果也顯示,該區域附近存在弱地殼形變,九寨溝地震前該區域存在一定程度的地殼形變閉鎖,可能是地震前的先兆特征。

圖5 研究區垂直形變速率等值線圖與梯度圖
太平洋板塊相對于亞歐板塊的NW向俯沖與印度板塊北向俯沖推擠亞歐板塊所形成的大規模的左旋走滑作用,是區域整體的地球動力學背景。作為活動塊體邊界的斷裂帶,由于其切割地殼深度大,塊體相互之間差異運動強烈而非連續性更強,更有利于應力的高強度積累而孕育大地震,進而導致絕大部分的強震活動都發生在活動塊體區,尤其是塊體邊界帶上(Nur et al,1989;張國民等,2000;張培震等,2003)。通過對巴顏喀拉塊體東部及鄰區主要7級以上強震的發震斷裂進行統計發現,這些強震主要集中在區內一些大型塊體的邊界活動斷裂帶上,如東昆侖、龍門山、鮮水河等深大斷裂帶。
九寨溝地震前,大面積水準和GNSS觀測資料均顯示,巴顏喀拉地塊東北部震中區域處于地殼形變閉鎖區,應力和能量在該區域深部持續積累,為地震的孕育發生提供了有利環境。通過跨斷層觀測資料可知,鮮水河斷裂帶整體上表現出左旋運動加速、張性運動明顯的特征,其中,鮮水河斷裂帶南東段出現的以張性運動為主、剪切活動較弱的特征,是由于階列式分段的鮮水河斷裂在走滑作用下造成了走滑系的尾端或者階區的伸展運動所致。龍門山斷裂帶表現出以壓性為主的運動學特征,加之阿爾金、東昆侖斷裂帶的左旋運動,可見九寨溝地震前,巴顏喀拉地塊在周緣巨型走滑斷裂帶的作用下,整體出現了順時針加速旋轉的運動特征。GPS速度場結果和地質學研究表明,巴顏喀拉地塊北東地區與高原東部均具有順時針向SE方向旋轉的特征(武艷強等,2013;楊國華等,2015;許志琴等,2016),并伴有運動增強的特點(熊維等,2015;姚鑫等,2017)。九寨溝地震前震中區地殼形變并不顯著,而在距震中區較遠的區域,存在地殼形變異常特征,與汶川地震、蘆山地震前的地殼形變特征類似(徐東卓等,2017)。同時,遠場構造帶上的場地震后恢復程度不等,這與場地所處的地塊屬性、構造體系和斷裂帶的不同段落特性等有關。震后未恢復,可能代表大區域的構造應力在繼續起作用,這種大區域的大規模應力也可能不會由于此次九寨溝7.0級地震立即發生轉變。如蘆山地震后,龍門山斷裂帶南段區域依然存在高應力積累。因不同的動力學環境、構造體系和微板塊屬性等,西秦嶺、六盤山、龍門山和鮮水河構造帶上的場地在九寨溝地震前后出現了不同的反應,甚至同一構造帶上的場地也會由于斷裂分段的復雜性而表現出了不同程度的變化特征。
中強地震前,多數情況下區域地殼會出現形變異常,同時,同一板塊周緣的多個強震之間會存在某種聯系。2000年以后,巴顏喀拉地塊周緣邊界斷裂帶上少見的強震多發的特點,揭示著該地塊進入了構造活動強烈期以及中強震多發幕期。隨著汶川8.0級地震等多次中強震的發生,巴顏喀拉地塊東邊界龍門山斷裂帶與北邊界東昆侖斷裂帶交匯的地塊東北區域成為了中強震缺少區,應力可能隨著汶川地震、蘆山地震、玉樹地震、于田地震、昆侖山口西地震轉移到了塊體東北部的東昆侖斷裂帶東端區域,且受到NE向的以逆沖擠壓運動為主的龍門山斷裂帶攔截影響,不斷在該區域積累,孕震環境形成。Wen等(2007)從地震空區角度分析也認為,作為巴顏喀拉地塊北東邊界的東昆侖斷裂帶東段是地震空區之一,地震危險性較高。從構造地質和巖石圈層理論研究來看,具有粘彈性的下地殼和上地幔弛豫作用可能會使得相距較遠的2次或者多次地震之間建立關聯性。中長期時間尺度上,大震會造成同一板塊深部塑性流動的殼幔物質發生調整運移,向應力易集聚區域集中,這會對相距較遠的斷裂活動產生作用,包括汶川地震、蘆山地震在內的多次強震,對九寨溝地震均產生了一定的促震作用。由震源破裂特征可見,巴顏喀拉地塊南北邊界斷裂帶上的地震多以走滑運動為主,東邊界的汶川地震則表現為強烈的逆沖特征為主,西邊界阿爾金斷裂帶的于田地震則表現出一定的拉張性質。因此,從動力學環境來看,這些強震均為巴顏喀拉地塊整體性運動的產物,伴隨著地塊的整體性東向活動,先后產生了一系列強震,并且強震之間存在一定的關聯。區域性大地震的發生,可引起同一動力學機制下相同板塊周緣邊界斷裂帶的應力調整。據以往震例研究,六盤山、西秦嶺構造區附近出現大量尖點、突跳型跨斷層場地異常,在青藏高原北緣斷裂帶,往往具備發生中強震的可能,并且中強震前西秦嶺、六盤山地區距震中300~400km內,甚至更遠地區,會出現大量跨斷層異常變化(王雙緒等,2012;張希等,2014;Xu et al,2016)。
青藏高原物質持續地向東流動,遇到了剛硬的四川盆地阻擋,隨后沿著龍門山斷裂帶分別向NE和SW運移集聚,NE方向的應力積累導致了九寨溝地震的發生,SW方向上龍門山斷裂帶的南端與鮮水河-安寧河斷裂帶交匯的“三岔口”區域,并未因汶川地震、蘆山地震的發生而使得區域地殼解鎖和能量釋放,而是隨著包括九寨溝地震在內的強震的發生,該區域庫侖應力可能進一步累積,依然存在中強震危險。在地震預測研究中,利用大面積水準資料得出的地殼形變梯度帶來分析7級以上強震發震地點具有一定準確性(張祖勝等,1996)。通過震例分析發現,強震基本上都發生在區域形變速率圖和梯度圖上所顯示的形變曲線密集、形變梯度較高的地區。“三岔口”區域存在地殼垂直形變高梯度帶,對未來7級以上地震有重要指示意義。同時,區域GPS速度場結果顯示,龍門山斷裂帶南段并未因汶川地震、蘆山地震的發生而大范圍解鎖,依然保持遠場加載、近場鎖定狀態(Wu et al,2013)。
基于九寨溝7.0級地震前后震中區周緣以及巴顏喀拉地塊周緣的跨斷層觀測資料、大面積水準測量資料、地震活動性資料以及其他研究結果等,綜合分析了地震前后的區域地殼形變以及巴顏喀拉地塊整體性運動與地震間的關系。
(1)九寨溝地震前后,遠場的西秦嶺構造區、六盤山斷裂帶、龍門山斷裂帶和鮮水河斷裂帶上的跨斷層流動場地均出現不同程度的異常變化,震中區域地殼形變不顯著。西秦嶺構造區在九寨溝地震前多表現出中短期壓性運動和異常波動特征,震后異常快速恢復。六盤山斷裂帶在九寨溝地震前的一段時期,整體表現出以壓性轉張性為主的運動學方式。九寨溝震前,龍門山斷裂帶表現為繼承性的明顯逆沖擠壓運動特征;鮮水河斷裂帶表現為左旋運動加速、同時表現出拉張活動特點。
(2)由巴顏喀拉地塊周緣的斷層異常活動和地震活動性分析可知,巴顏喀拉地塊整體性的SE向順時針加速運動,加快了地塊東北部九寨溝震中區域的應力積累,孕育了此次地震,同時,包括蘆山地震在內的幾次強震對九寨溝地震有一定的促震作用。
(3)九寨溝地震后,區域應力積累得到一定程度的釋放,具備相同動力學背景的 “三岔口”區域應力可能會加速積累。
致謝:感謝兩位審稿專家和編輯老師的辛苦審閱,薄萬舉研究員給予了有益的建議,部分圖件使用GMT軟件(Wessel et al,1995)繪制,一并表示感謝。