摘要：為研究2022年6月1日四川蘆山MS6.1地震的孕育和發(fā)生過程，采用CAP方法反演了2013年蘆山MS7.0主震及MS≥5.0余震的震源機制解，并基于應(yīng)力張量方差與b值時空分布特征，探討了蘆山MS6.1地震的力學(xué)機制和震源區(qū)的應(yīng)力狀態(tài)。結(jié)果表明：2022年蘆山MS6.1地震震源機制表現(xiàn)出與2013年蘆山MS7.0主震和5級余震相似的逆沖型破裂特征，壓應(yīng)力軸方位與龍門山斷裂帶南段區(qū)域應(yīng)力場一致。2013年蘆山MS7.0地震后震中及附近的應(yīng)力張量方差和b值長期處于低值狀態(tài)，2022年蘆山MS6.1地震前震中及附近出現(xiàn)了應(yīng)力張量方差和b值的低值異常，表明蘆山余震區(qū)處于較高的應(yīng)力水平。分析認(rèn)為：巴顏喀拉塊體持續(xù)東向運動受到華南塊體的阻擋，震中所在區(qū)域長期受擠壓逆沖作用，從而使蘆山余震區(qū)長期處于應(yīng)力積累的狀態(tài)，蘆山MS6.1地震也是在這種動力學(xué)背景下發(fā)生的。

關(guān)鍵詞：蘆山MS6.1地震；蘆山MS7.0地震；震源機制解；應(yīng)力張量方差；b值；應(yīng)力狀態(tài)

中圖分類號：P315.72文獻標(biāo)識碼：A文章編號：1000-0666（2024）04-0483-10

doi：10.20015/j.cnki.ISSN1000-0666.2024.0061

0引言

據(jù)中國地震臺網(wǎng)測定，2022年6月1日四川省雅安市蘆山縣（30.37°N，102.94°E）發(fā)生MS6.1地震，震源深度17 km。此次地震位于2013年4月20日蘆山MS7.0地震的西北側(cè)，震中僅相距9 km，屬于蘆山MS7.0地震的晚期強余震。

2013年蘆山MS7.0地震發(fā)生后，眾多學(xué)者就蘆山地震序列震源機制及震區(qū)的應(yīng)力狀態(tài)開展了相關(guān)研究（單斌等，2013；羅艷等，2015；梁春濤等，2018）。蘆山MS7.0地震震源機制斷層面走向NE、傾向NW，為NWW-SEE向水平擠壓應(yīng)力作用下發(fā)生的一次逆沖型地震（劉杰等，2013）；蘆山余震區(qū)最大主壓應(yīng)力（S1）方向存在局部空間差異，余震區(qū)南段S1方向總體呈NWW向，而余震區(qū)北段S1方向表現(xiàn)出由NW經(jīng)EW向NE的逆時針旋轉(zhuǎn)（張致偉等，2015a）；單斌等（2013）計算了蘆山MS7.0地震的同震靜態(tài)庫侖應(yīng)力場變化，發(fā)現(xiàn)主震破裂面附近以及破裂面的東北和西南端為庫侖應(yīng)力增強區(qū)，余震也主要發(fā)生在破裂面附近的庫侖應(yīng)力增強區(qū)；易桂喜等（2013）基于蘆山地震前的資料獲得了龍門山南段的b值圖像，顯示震前在天全—蘆山段應(yīng)力水平相對較高；梁春濤等（2018）利用密集監(jiān)測臺陣的走時定位發(fā)現(xiàn)蘆山余震區(qū)東邊的大邑斷裂及一些隱伏斷裂上地震比較活躍，認(rèn)為應(yīng)力在這些斷裂上有所積累；張致偉等（2023）研究認(rèn)為2022年蘆山MS6.1地震序列發(fā)生在2013年蘆山MS7.0地震發(fā)震斷層北西側(cè)的一條傾向南東的反沖斷層上，且MS7.0地震對MS6.1地震有顯著的觸發(fā)作用。

地震的發(fā)生與應(yīng)力狀態(tài)密切相關(guān)，強震震源區(qū)應(yīng)力狀態(tài)的時空變化對了解地震過程和預(yù)測地震具有重要意義。震源機制一致性和b值作為表征區(qū)域應(yīng)力狀態(tài)的參數(shù)，被廣泛應(yīng)用于地震危險性分析和地震預(yù)測研究中（Wiemer et al，2002；Wyss，Stefansson，2006；易桂喜等，2013；曾憲偉等，2020）。在大地震孕育的過程中，震源區(qū)附近的構(gòu)造應(yīng)力場控制作用會增強，特別是臨近發(fā)震的時候，中小地震的震源機制趨于一致，應(yīng)力張量方差顯著降低，反映了震源區(qū)應(yīng)力水平升高；大震后震源區(qū)應(yīng)變能得以充分釋放，中小地震的震源機制趨于紊亂，應(yīng)力張量方差明顯升高，表明震源區(qū)應(yīng)力水平降低（Wiemer et al，2002；付虹等，2011；張致偉等，2015b）。b值本質(zhì)上反映的是不同震級地震之間的比例關(guān)系，代表了介質(zhì)內(nèi)部應(yīng)力水平的高低，其隨著介質(zhì)應(yīng)力水平的增加而減小，介質(zhì)應(yīng)力水平高，在巖石破裂面的邊界上處于高水平的應(yīng)力點所占的比重越大，大破裂的比例也越大，b值越??；介質(zhì)應(yīng)力水平低，破裂面的邊界上就越容易碰上低水平的應(yīng)力點而停止，小破裂所占的比例大，b值越高（Wyss，Wiemer，2000；易桂喜等，2013）。

時隔9年蘆山MS7.0地震余震區(qū)再次發(fā)生MS6.1地震，震源力學(xué)機制與震源區(qū)應(yīng)力狀態(tài)是理解蘆山MS6.1地震孕育和發(fā)生的關(guān)鍵問題之一。本文采用CAP方法反演了2013年以來蘆山MS≥5.0地震震源機制解，基于余震區(qū)應(yīng)力張量方差與b值時空分布特征，探討了蘆山MS6.1地震前震源區(qū)的應(yīng)力狀態(tài)。

地震研究47卷第4期張致偉等：2022年四川蘆山MS6.1地震前應(yīng)力狀態(tài)研究1計算方法

1.1震源機制反演方法

1.2b值計算方法

2地震序列特征

3震源機制和應(yīng)力張量方差特征

3.1蘆山MS≥5.0余震震源機制特征

3.2應(yīng)力張量方差特征

為了求解空間非均勻震源機制解的應(yīng)力場特征，Michael（1991）提出了疊加應(yīng)力場反演方法，它通過在均勻應(yīng)力場上疊加擾動來模擬非均勻的空間應(yīng)力場，獲得張應(yīng)力軸、中等應(yīng)力軸和壓應(yīng)力軸的方向以及反演方差。應(yīng)力張量方差被定義為單個地震的滑動矢量與在假設(shè)應(yīng)力張量作用下產(chǎn)生的理論滑動矢量之間的夾角與其平均數(shù)之差

基于2013年4月20日—2022年5月31日蘆山129次MS≥3.0余震震源機制解（圖4a），將蘆山余震區(qū)劃分為0.01°×0.01°的網(wǎng)格，選取每個網(wǎng)格節(jié)點及其周圍至少8個地震震源機制解，采用Michael（1991）提出的應(yīng)力場反演方法計算了蘆山余震區(qū)的應(yīng)力張量方差（圖4b）。2022年6月1日蘆山MS6.1地震前，余震區(qū)呈現(xiàn)的應(yīng)力張量方差分布特征與早期的研究結(jié)果（張致偉等，2015b）較為一致。蘆山MS7.0地震兩側(cè)的應(yīng)力張量方差存在明顯差異，主震震中及其北側(cè)的應(yīng)力張量方差明顯低于南側(cè)，余震區(qū)北側(cè)的應(yīng)力張量方差小于0.1，說明該區(qū)域的應(yīng)力水平較高，MS6.1地震就發(fā)生在余震區(qū)北側(cè)的應(yīng)力張量方差低值區(qū)。以往研究也有類似的發(fā)現(xiàn)：Wiemer等（2002）分析認(rèn)為1992年蘭德斯MS7.3和大熊MS6.4地震均發(fā)生在應(yīng)力張量方差低值區(qū)；張致偉等（2015b）研究發(fā)現(xiàn)汶川MS≥5.5強余震也發(fā)生在應(yīng)力張量方差低值分布區(qū)或其邊緣附近；李金等（2015）研究認(rèn)為天山中東段地區(qū)的MS≥5.5地震大多發(fā)生在應(yīng)力張量方差相對較低的區(qū)域，新源—和靜MS6.6地震位于低值區(qū)中心。

4余震區(qū)b值時空變化特征

4.1余震區(qū)b值時間變化特征

基于2000年1月1日—2022年12月31日蘆山余震區(qū)的地震，采用最大似然法估算了蘆山余震區(qū)整個時段的b值（圖5）。鑒于震級-頻次關(guān)系式中高、低震級段出現(xiàn)的“擺尾”和“掉頭”現(xiàn)象會對b值計算產(chǎn)生影響（曾憲偉等，2020），本文在計算b值時同時考慮了高震級截斷地震和最小完整性震級的影響。根據(jù)震級-頻次關(guān)系，選取最小完整性震級MC為0.8，高震級的截斷震級為5.5，扣除高、低震級段的影響后，計算獲得的蘆山余震區(qū)的平均b值為0.75±0.06。

以2013年蘆山MS7.0地震序列為研究對象，選取500個地震樣本作為一個時間窗，且步長取為一個時間窗，逐步滑動樣本窗計算每個窗口內(nèi)的b值，繪制了2013年蘆山MS7.0地震以來余震區(qū)的b值時間進程曲線（圖6）。結(jié)果顯示，地震后b值出現(xiàn)短時間上升，2013年7月達到最高值0.9后開始下降，2018年下降到最低值0.68后開始緩慢上升直至發(fā)生2022年蘆山MS6.1地震。馬鴻慶（1978）在研究華北地區(qū)大震前的b值異常變化時，同樣發(fā)現(xiàn)b值隨時間先升高后降低，達到最低值后又回升直至發(fā)震；高雅婧等（2022）研究認(rèn)為臺灣集集地震發(fā)生后十天至5個月內(nèi)b值迅速且連續(xù)升高，6～7個月后b值出現(xiàn)降低趨勢，表明區(qū)域應(yīng)力水平經(jīng)歷了恢復(fù)后又處于積累的狀態(tài)。

2016—2021年蘆山余震區(qū)的b值低于2000年以來的平均b值（0.75）水平，表明蘆山MS7.0地震后余震區(qū)長期處于低b值狀態(tài)。史海霞等（2018）研究發(fā)現(xiàn)汶川震源區(qū)b值從2002年開始至地震前呈現(xiàn)一個長期趨勢性下降；易桂喜等（2013）分析了龍門山斷裂帶南段應(yīng)力狀態(tài)，認(rèn)為2008年汶川MS8.0地震后，天全—蘆山以及寶興北部地區(qū)b值有較明顯的降低；蔣海昆等（2000）通過花崗巖變形過程中的聲發(fā)射實驗，研究認(rèn)為b值的系統(tǒng)降低表征研究區(qū)所處環(huán)境應(yīng)力的進一步增強。分析認(rèn)為，蘆山余震區(qū)長期處于低b值的現(xiàn)象可能源于巴顏喀拉塊體持續(xù)東向運動受到華南塊體的阻擋，余震區(qū)長期受擠壓逆沖作用所致。

4.2余震區(qū)b值空間分布特征

將蘆山余震區(qū)劃分為0.01°×0.01°的網(wǎng)格，以每個網(wǎng)格節(jié)點為中心點，搜索半徑逐步增大，當(dāng)?shù)卣饠?shù)達到200個時停止搜索并計算b值。根據(jù)余震區(qū)b值時間變化過程，考慮到蘆山MS7.0地震對余震區(qū)b值的影響，本文選取b值相對穩(wěn)定的2016年作為起始時間，將研究時段劃分為2016年1月1日至2020年12月31日（第I時段）、2016年1月1日至2022年5月31日（第II時段）和2016年1月1日至2022年12月31日（第III時段）共3個時段分析b值空間掃描圖像（圖7）。

余震區(qū)b值在不同時段的空間圖像顯示，3個時段蘆山MS7.0主震震中及附近b值均處于低值（＜0.75）狀態(tài)，表明主震震中及附近長期處于應(yīng)力積累。2022年蘆山MS6.1地震發(fā)生前，震中附近的b值（圖7b）較第I時段（圖7a）有所下降，地震發(fā)生后短期內(nèi)震源區(qū)仍表現(xiàn)為低b值（圖7c）。易桂喜等（2011）研究發(fā)現(xiàn)龍門山斷裂帶中北段的綿竹—茂縣段與江油—平武段具有異常低b值，這兩個段落也是汶川5級以上強余震集中發(fā)生的區(qū)域；Nanjo等（2012）研究發(fā)現(xiàn)2011年日本東北MS9.0和2004年蘇門答臘MS9.1兩次地震前震源區(qū)均存在異常低b值；馮建剛等（2016）研究發(fā)現(xiàn)2013年岷縣漳縣MS6.6地震發(fā)生在臨潭—宕昌斷裂與西秦嶺北緣斷裂之間的低b值區(qū)邊緣。對比前文給出的應(yīng)力張量方差空間分布，可以看出2022年蘆山MS6.1地震前，震源區(qū)的低b值與應(yīng)力張量方差低值區(qū)吻合較好。由此可見，應(yīng)力張量方差和b值均可反映區(qū)域應(yīng)力狀態(tài)，結(jié)合應(yīng)力張量方差和b值的空間分布可作為識別強震危險地點的方法之一。

5結(jié)論

基于四川地震臺網(wǎng)記錄的觀測資料，采用CAP方法反演了2013年蘆山MS7.0地震及MS≥5.0余震震源機制解，并基于應(yīng)力張量方差及b值時空分布特征，探討了2022年蘆山MS6.1地震的力學(xué)機制及余震區(qū)的應(yīng)力狀態(tài)。獲得的主要認(rèn)識如下：

（1）2022年蘆山MS6.1地震發(fā)生在2013年蘆山MS7.0地震余震區(qū)北西側(cè)約9 km處，屬于蘆山MS7.0地震的晚期強余震。蘆山MS6.1地震震源機制斷層面走向、傾角和滑動角分別為22°、52°和78°，P軸方位120°，傾角6°。表現(xiàn)出與蘆山MS7.0主震和5級余震相似的逆沖型破裂特征，壓應(yīng)力軸方位與龍門山斷裂帶南段區(qū)域應(yīng)力場一致。揭示蘆山MS6.1地震仍是在巴顏喀拉塊體東向擠壓的動力學(xué)背景下發(fā)生的。

（2）2013年蘆山MS7.0地震后，主震震中附近始終表現(xiàn)為較低的應(yīng)力張量方差和b值，其中應(yīng)力張量方差小于0.1，b值小于0.75。2016—2021年期間蘆山余震區(qū)的b值均低于2000年以來的平均b值（0.75）水平，反映出蘆山MS7.0地震震中及附近長期處于應(yīng)力積累狀態(tài)。上述現(xiàn)象可能源于巴顏喀拉塊體持續(xù)東向運動受到華南塊體的阻擋，震中所在區(qū)域長期受擠壓逆沖作用所致。

（3）2022年蘆山MS6.1地震前，震中及附近的應(yīng)力張量方差和b值均呈現(xiàn)低值異常，表明蘆山MS6.1地震前震源區(qū)處于較高的應(yīng)力水平。由此可見，應(yīng)力張量方差和b值均可作為表征區(qū)域應(yīng)力狀態(tài)的參數(shù)，結(jié)合兩者的空間分布可用于識別強震危險地點。

四川地震臺為本研究提供了地震目錄和波形數(shù)據(jù)，CAP波形反演程序來源于美國圣路易斯大學(xué)的朱露培教授，b值計算和繪圖程序來自ZMAP程序包。在此一并表示衷心感謝！

參考文獻：

馮建剛，張輝，楊萍.2016.2013年岷縣漳縣6.6級地震前地震b值異常特征研究［J］.地震，36（1）：32-37.

Feng J G，Zhang H，Yang P.2016.Anormalies of b value Changes before the 2013 Minxian-Zhangxian MS6.6 Earthquake［J］.Earthquake，36（1）：32-37.（in Chinese）

付虹，王紹晉，李麗，等.2011.滇中地區(qū)震源機制一致性參數(shù)時空分布與強震活動［J］.地震研究，34（2）：114-118.

Fu H，Wang S J，Li L，et al.2011.Spatial-temporal distribution of the parameters of the focal mechanism and strong seismicity in central Yunnan［J］.Journal of Seismological Research，34（2）：114-118.（in Chinese）

高雅婧，孫云強，羅綱.2022.1999年集集地震前后臺灣地區(qū)地震b值及應(yīng)力場時空演化特征［J］.地球物理學(xué)報，65（6）：2137-2152.

Gao Y J，Sun Y Q，Luo G.2022.Temporal and spatial evolution characteristics of b value and stress field in Taiwan before and after the 1999 Chi-Chi earthquake［J］.Chinese Journal of Geophysics，65（6）：2137-2152.（in Chinese）

蔣海昆，張流，周永勝.2000.不同圍壓條件下花崗巖變形破壞過程中的聲發(fā)射時序特征［J］.地球物理學(xué)報，43（6）：812-826.

Jiang H K，Zhang L，Zhou Y S.2000.Behavior of acoustic emission time sequence of granite in deformation and failure process under different confining pressures［J］.Chinese Journal of Geophysics，43（6）：812-826.（in Chinese）

闞榮舉，張四昌，晏鳳桐，等.1977.我國西南地區(qū)現(xiàn)代構(gòu)造應(yīng)力場與現(xiàn)代構(gòu)造活動特征的探討［J］.地球物理學(xué)報，20（2）：96-109.

Kan R J，Zhang S C，Yan F T，et al.1977.Present tectonic stress field and its relation to the characteristics of recent tectonic activity in southwestern China［J］.Chinese Journal of Geophysics，20（2）：96-109.（in Chinese）

李金，周龍泉，龍海英，等.2015.天山地震帶（中國境內(nèi)）震源機制一致性參數(shù)的時空特征［J］.地震地質(zhì)，37（3）：792-803.

Li J，Zhou L Q，Long H Y，et al.2015.Spatial-temporal characteristics of the focal mechanism consistency parameter in TianShan（within Chinese territory）seismic zone［J］.Seismology and Geology，37（3）：792-803.（in Chinese）

梁春濤，黃焱羚，王朝亮，等.2018.汶川和蘆山地震之間地震空區(qū)綜合研究進展［J］.地球物理學(xué)報，61（5）：1996-2010.

Liang C T，Huang Y L，Wang C L，et al.2018.Progress in the studies of the seismic gap between the 2008 Wenchuan and 2013 Lushan earthquakes［J］.Chinese Journal of Geophysics，61（5）：1996-2010.（in Chinese）

劉杰，易桂喜，張致偉，等.2013.2013年4月20日四川蘆山M7.0級地震介紹［J］.地球物理學(xué)報，56（4）：1404-1407.

Liu J，Yi G X，Zhang Z W，et al.2013.Introduction to the Lushan，Sichuan M7.0 earthquake on 20 April 2013［J］.Chinese Journal of Geophysics，56（4）：1404-1407.（in Chinese）

呂堅，鄭勇，倪四道，等.2008.2005年11月26日九江-瑞昌MS5.7、MS4.8地震的震源機制解與發(fā)震構(gòu)造研究［J］.地球物理學(xué)報，51（1）：158-164.

Lyu J，Zheng Y，Ni S D，et al.2008.Focal mechanisms and seismogenic structures of the MS5.7 and MS4.8 Jiujiang-Ruichang earthquakes of Nov.26，2005［J］.Chinese Journal of Geophysics，51（1）：158-164.（in Chinese）

羅艷，趙里，曾祥方，等.2015.蘆山地震序列震源機制及其構(gòu)造應(yīng)力場空間變化［J］.中國科學(xué)：地球科學(xué)，45：538-550.

Luo Y，Zhao L，Zeng X F，et al.2015.Focal mechanisms of the Lushan earthquake sequence and spatial variation of the stress field［J］.Scientia Sinica Terrae，45：538-550.（in Chinese）

馬鴻慶.1978.華北地區(qū)幾次大震前的b值異常變化［J］.地球物理學(xué)報，21（2）：126-141.

Ma H Q.1978.Variations of the b-values before several large earthquakes occurred in north China［J］.Chinese Jaurnal of Geohysics，21（2）：126-141.（in Chinese）

單斌，熊熊，鄭勇，等.2013.2013年蘆山地震導(dǎo)致的周邊斷層應(yīng)力變化及其與2008年汶川地震的關(guān)系［J］.中國科學(xué)：地球科學(xué)，43（6）：1002-1009.

Shan B，Xiong X，Zheng Y，et al.2013.Stress changes on major faults caused by 2013 Lushan earthquake，and its relationship with 2008 Wenchuan earthquake［J］.Scientia Sinica Terrae，43（6）：1002-1009.（in Chinese）

史海霞，孟令媛，張雪梅，等.2018.汶川地震前的b值變化［J］.地球物理學(xué)報，61（5）：1874-1882.

Shi H X，Meng L Y，Zhang X M，et al.2018.Decrease in b value prior to the Wenchuan earthquake（MS8.0）［J］.Chinese Journal of Geophysics，61（5）：1874-1882.（in Chinese）

徐錫偉，陳桂華，于貴華，等.2013.蘆山地震發(fā)震構(gòu)造及其與汶川地震關(guān)系討論［J］.地學(xué)前緣，20（3）：11-20.

Xu X W，Chen G H，Yu G H，et al.2013.Seismogenic structure of Lushan earthquake and its relationship with Wenchuan earthquake［J］.Earth Science Frontiers，20（3）：11-20.（in Chinese）

易桂喜，聞學(xué)澤，辛華，等.2011.2008年汶川MS8.0地震前龍門山-岷山構(gòu)造帶的地震活動性參數(shù)與地震視應(yīng)力分布［J］.地球物理學(xué)報，54（6）：1490-1500.

Yi G X，Wen X Z，Xin H，et al.2011.Distributions of seismicity parameters and seismic apparent stresses on the Longmenshan-Minshan tectonic zone before the 2008 MS8.0 Wenchuan earthquake［J］.Chinese Journal of Geophysics，54（6）：1490-1500.（in Chinese）

易桂喜，聞學(xué)澤，辛華，等.2013.龍門山斷裂帶南段應(yīng)力狀態(tài)與強震危險性研究［J］.地球物理學(xué)報，56（4）：1112-1120.

Yi G X，Wen X Z，Xin H，et al.2013.Stress state and major-earthquake risk on the southern segment of the Longmen Shan fault zone［J］.Chinese Journal of Geophysics，56（4）：1112-1120.（in Chinese）

曾憲偉，龍鋒，任家琪，等.2020.2019年6月17日長寧MS6.0地震前后b值時空變化分析［J］.地震，40（3）：1-14.

Zeng X W，Long F，Ren J Q，et al.2020.Spatial and temporal variation of b value before and after the Changning MS6.0 earthquake on June 17，2019［J］.Earthquake，40（3）：1-14.（in Chinese）

張致偉，龍鋒，石富強，等.2023.2022年6月1日四川蘆山MS6.1地震的發(fā)震構(gòu)造與力學(xué)機制探討［J］.地球物理學(xué)報，66（10）：4095-4110.

Zhang Z W，Long F，Shi F Q，et al.2023.Discussion on seismogenic structure and mechanical mechanism of the 1 June 2022 MS6.1 Sichuan Lushan earthquake［J］.Chinese Journal of Geophysics，66（10）：4095-4110.（in Chinese）

張致偉，周龍泉，程萬正，等.2015a.蘆山MW6.6地震序列的震源機制及震源區(qū)應(yīng)力場［J］.地球科學(xué)，40（10）：1710-1722.

Zhang Z W，Zhou L Q，Cheng W Z，et al.2015a.Focal mechanism solutions of Lushan MW6.6 earthquake sequence and stress field for aftershock zone［J］.Earth Science，40（10）：1710-1722.（in Chinese）

張致偉，周龍泉，龍鋒，等.2015b.汶川8.0和蘆山7.0級地震序列應(yīng)力場時空特征［J］.地震地質(zhì)，37（3）：804-817.

Zhang Z W，Zhou L Q，Long F，et al.2015b.Spatial and temporal characteristic of stress field for Wenchuan MS8.0 and Lushan MS7.0 earthquake sequence［J］.Seismology and Geology，37（3）：804-817.（in Chinese）

鄭勇，馬宏生，呂堅，等.2009.汶川地震強余震（MS≥5.6）的震源機制解及其與發(fā)震構(gòu)造的關(guān)系［J］.中國科學(xué)：地球科學(xué)，39（4）：413-426.

Zheng Y，Ma H S，Lü J，et al.2009.Source mechanism of strong aftershocks（MS≥5.6）of the Wenchuan earthquake and the implication for seismotectonics［J］.Scientia Sinica Terrae，39（4）：413-426.（in Chinese）

中國地震局監(jiān)測預(yù)報司.2020.測震分析預(yù)測技術(shù)方法工作手冊［M］.北京：地震出版社.

Department of Monitoring and Forecasting，China Earthquake Administration.2020.Working manual of seismic analysis and prediction technology［M］.Beijing：Seismological Press.（in Chinese）

Hubbard J，Shaw J H.2009.Uplift of the Longmen Shan and Tibetan Plateau，and the 2008 Wenchuan（MW7.9）earthquake［J］.Nature，458：194-197.

Lu Z，Wyss M，Pulpan H.1997.Detail of stress directions in the Alaska Subduction Zone from fault plane solutions［J］.Journal of Geophysical Research，102（B3）：5385-5402.

Michael A J.1987.Use of focal mechanisms to determine stress：A control study［J］.Journal of Geophysical Research，92（B1）：357-368.

Michael A J.1991.Spatial variations in stress within the 1987 Whittier Narrows，California，aftershock sequence：New techniques and results［J］.Journal of Geophysical Research，96（B4）：6303-6319.

Nanjo K Z，Hirata N，Obara K，et al.2012.Decade-scale decrease in b value prior to the M9-class 2011 Tohoku and 2004 Sumatra quakes［J］.Geophys Res Lett，39：L20304.

Shi Y，Bolt B A.1982.The standard error of the magnitude-frequency b-value［J］.Bull Seism Soc Am，72（5）：1677-1687.

Tan Y，Zhu L，Helmberger D V，et al.2006.Locating and modeling regional earthquakes with two stations［J］.J Geophys Res，111：B11.

Wiemer S，Gerstenberger M，Hauksson E.2002.Properties of the aftershock sequence of the 1999 MW7.1 Hector Mine earthquake：Implications for aftershock hazard［J］.Bull Seism Soc Amer，92（4）：1227-1240.

Wyss M，Stefansson R.2006.Nucleation points of recent mainshocks in Southern Iceland，mapped by b-values［J］.Bull Seism Soc Am，96（2）：599-608.

Wyss M，Wiemer S.2000.Change in the probability for earthquakes in southern California due to the Landers magnitude 7.3 earthquake［J］.Science，290（5495）：1334-1338.

Zhang Z W，Liang C T，Long F，et al.2022.Spatiotemporal variations of focal mechanism solutions and stress mechanism solutions and stress field of the 2019 Changning MS6.0 Earthquake sequence［J］.Front Earth Sci，9：797907.

Study on the Stress State before the 2022，Lushan，Sichuan MS6.1 Earthquake

ZHANG Zhiwei1，2，ZENG Xianwei3，WANG Di1，LU Qian1，WANG Weiming1，YANG Peng1，LONG Feng1

（1.Sichuan Earthquake Agency，Chengdu 610041，Sichuan，China）

（2.College of Geophysics，Chengdu University of Technology，Chengdu 610059，Sichuan，China）

（3.Earthquake Agency of Ningxia Hui Autonomous Region，Yinchuan 750001，Ningxia，China）

Abstract

On June 1，2022，an MS6.1 earthquake struck Lushan Country，Yaan City of Sichuan Province.The earthquake is another significant event that occurred in the aftershock zone of the Lushan MS7.0 earthquake on May 20，2013.To figure out the genesis and occurrence of the Lushan MS6.1 earthquake，the CAP inversion method is used to calculate the focal mechanism of the Lushan MS7.0 mainshock and its MS≥5.0 aftershocks.And the mechanical mechanism and the stress state in the source area before the Lushan MS6.1 earthquake are discussed based on the spatial and temporal distribution of the stress tensor variance and b value.The main insights are as follows：The focal mechanism of the Lushan MS6.1 earthquake is featured as a thrust rupture，which is similar to the feature of the focal mechanism of the Lushan MS7.0 earthquake and its MS≥5.0 aftershocks.The direction of the principal compressive stress is consistent with the regional stress field of the southern section of the Longmenshan fault zone.The stress tensor variance and b value around the epicenter of the Lushan MS7.0 earthquake are in a low state in a long period.What is more，shortly before the Lushan MS6.1 earthquake，in its focal area there appeared low-value anomalies of the stress tensor variance and b value，indicating that the Lushan MS7.0 earthquake aftershock zone had been keeping a higher stress level.The continuous eastward movement of the Bayan Har block is blocked by the South China block.Thus，the aftershock zone of the Lushan MS7.0 earthquake is under the action of the horizontal extrusion and thrust，which leads to the long-term stress accumulation in this zone.Hence in this dynamic background，the Lushan MS6.1 earthquake occurred.

Keywords：the Lushan MS6.1 earthquake；the Lushan MS7.0 earthquake；focal mechanism solution；stress tensor variance；b value；stress state