馮楚豪, 嚴(yán)月天, 馮萬(wàn)鵬*, 王雨晴, 陳大慶, 吳傳勇
1 中山大學(xué)地球科學(xué)與工程學(xué)院, 廣東珠海 519000 2 南方海洋科學(xué)與工程廣東省實(shí)驗(yàn)室, 廣東珠海 519000 3 廣東省地震局, 廣州 510000
在距今約50~55 Ma前,印度板塊與歐亞板塊相互碰撞,造就了宏偉的喜馬拉雅造山帶,并快速促進(jìn)了青藏高原地區(qū)的持續(xù)抬升,最終形成現(xiàn)今世界第三極——青藏高原(Molnar and Tapponnier, 1975).在這持續(xù)的擠壓造山過(guò)程中,青藏高原內(nèi)部還同時(shí)廣泛發(fā)育了一系列啟動(dòng)于中新世中期至第四紀(jì)的近南北展布的裂谷,這些伸展構(gòu)造不僅調(diào)節(jié)了印度板塊與歐亞板塊的碰撞變形,同時(shí)也與青藏高原的隆升歷史密切相關(guān).因此,識(shí)別裂谷帶區(qū)域現(xiàn)今變形特征對(duì)理解青藏高原的生長(zhǎng)與變形模式有著重要的意義(張進(jìn)江和丁林, 2003).
對(duì)于上述裂谷成因,不同學(xué)者提出了多種可能解釋,目前主要可以概括為兩種主流模式(Wang et al., 2014),一種為與重力垮塌相關(guān)的模式(Coleman and Hodges, 1995),另一種為與區(qū)域應(yīng)力場(chǎng)及邊界條件相關(guān)的模式(Cogan et al., 1998),目前尚未達(dá)到廣泛統(tǒng)一.過(guò)去的30年間,地震學(xué)和空間對(duì)地觀測(cè)技術(shù)記錄了大量的發(fā)生在青藏高原內(nèi)部的正斷層地震,這些資料為我們定量化斷層結(jié)構(gòu)和形變過(guò)程奠定了堅(jiān)實(shí)的數(shù)據(jù)基礎(chǔ).2020年6月藏北于田地區(qū)再次發(fā)生了中強(qiáng)正斷層地震,為嘗試探索當(dāng)?shù)亓压壬L(zhǎng)規(guī)律提供了新的機(jī)遇.
據(jù)中國(guó)地震臺(tái)網(wǎng)測(cè)定,2020年6月26日5時(shí)05分(UTC+8)新疆于田地區(qū)發(fā)生了6.4級(jí)的地震,震源深度約10 km,震中位于82.33°E,35.73°N.美國(guó)地質(zhì)調(diào)查局(USGS)的結(jié)果顯示,于田地震可能地震斷層節(jié)面的走向、傾角和滑動(dòng)角分別為(24°,42°,-108°)和(227°,50°,-74°),震級(jí)大小為MW6.3.此次地震發(fā)生的位置距離居民區(qū)較遠(yuǎn),于田縣無(wú)明顯震感,地震未造成嚴(yán)重?fù)p害.
2020年于田地震發(fā)生于青藏高原巴顏喀拉塊體西緣,黑石北湖西側(cè)的正斷層系統(tǒng)附近(圖1a).GPS速度場(chǎng)顯示該區(qū)域現(xiàn)今處于東西拉張構(gòu)造環(huán)境(葛偉鵬等, 2015).2000—2020年間,當(dāng)?shù)毓灿?000余起地震事件被地震臺(tái)網(wǎng)所記錄到,顯示出該區(qū)域較強(qiáng)的地震活動(dòng)性.2008—2020年間,于田地震周邊地區(qū)相繼發(fā)生了4次M>6級(jí)的地震事件,其中3次為正斷事件(圖1),包括2008年MW7.2地震,2012年MW6.2 地震和2020年MW6.3地震(He et al., 2020; 康帥等, 2021),另1次為2014年2月的MW6.9走滑事件(Li et al., 2020;羅毅等, 2021).該地區(qū)北部邊界為阿爾金斷裂西南段—康西瓦斷裂,南部邊界為郭扎錯(cuò)斷裂(圖1),兩斷裂同為左旋走滑斷裂系統(tǒng),滑動(dòng)速率分別為5~10 mm·a-1和10 mm·a-1(Wright et al., 2004).穿過(guò)震區(qū)地形剖面AA′(圖1b)顯示,研究區(qū)內(nèi)海拔高度均在4000 m以上,地勢(shì)起伏多變,山峰與溝谷交錯(cuò)分布,暗示了當(dāng)?shù)卦诘刭|(zhì)歷史中經(jīng)歷了強(qiáng)烈的構(gòu)造變形過(guò)程.相對(duì)于藏南廣泛發(fā)育的南北走向裂谷帶以及相生的正斷層地震事件,位于藏北的于田地區(qū)正斷事件相對(duì)集中,很可能與上述兩個(gè)走滑系統(tǒng)相關(guān).因此,開(kāi)展該地區(qū)正斷層地震的詳細(xì)研究,將有助于進(jìn)一步理解青藏高原的擴(kuò)展模式.
圖1 2020于田地震構(gòu)造背景場(chǎng)和SAR數(shù)據(jù)空間覆蓋(a) 紅星代表于田地震震中;綠線代表剖面圖位置;紅箭頭為GPS速度場(chǎng)(Zhao et al., 2015);彩色點(diǎn)為于田地震精定位后的余震(冉慧敏等, 2020).(b)橫跨2008年于田MW7.3地震及2020年于田MW6.3地震的地形剖面圖.Fig.1 Tectonic background of the 2020 MW6.3 Yutian earthquake and SAR spatial coverage(a) The red star is the epicenter of the Yutian earthquake. The green line represents a local profile A—A′, crossing the footwall and hanging walls of the earthquake fault. The red arrow represents GPS velocity (Zhao et al., 2015) in the regions; The colored dots represent the relocated aftershocks (Ran et al., 2020) of the Yutian earthquake. (b) Surface topography along Profile AA′ and coseismic displacements of the 2020 mainshock. The profile AA′ goes through both the 2008 MW7.3 and 2020 MW6.3 Yutian earthquakes areas.
InSAR地表形變資料已經(jīng)成為研究斷層運(yùn)動(dòng)狀態(tài),約束斷層幾何參數(shù)的重要手段(張文婷等,2020;季靈運(yùn)等,2021).本文中,我們將利用InSAR技術(shù)獲得此地震的同震和震后地表變形序列,并借助地球物理反演方法得出于田地震的斷層參數(shù)和滑動(dòng)分布,并以此為依據(jù),整合前人的研究結(jié)果,進(jìn)一步探究了藏北裂谷生長(zhǎng)模式與特征.
此次研究中,我們收集了沿Sentinel-1(S1) T165降軌和T158升軌軌道TOPS模式的雷達(dá)圖像(表1),包括震后6個(gè)月內(nèi)分別沿T165和T158軌道的21副和20副雷達(dá)數(shù)據(jù).同時(shí),為了增加同震測(cè)量的可靠性,我們還收集了一副ALOS2衛(wèi)星L波段的同震像對(duì)(表1).因ALOS2衛(wèi)星對(duì)當(dāng)?shù)爻上耦l率不足,只找到時(shí)間間隔約1年的同震像對(duì),并包含了約3個(gè)月的震后過(guò)程.
在早期NRCAN/CCRS開(kāi)發(fā)的InSAR自動(dòng)化處理環(huán)境pSAR基礎(chǔ)上(Feng et al., 2016),我們使用開(kāi)源GMTSAR軟件替代早期商業(yè)軟件(Sandwell et al., 2011; Xu et al., 2018),對(duì)TOPS-SAR數(shù)據(jù)開(kāi)展D-InSAR自動(dòng)化處理.處理中,30 m格網(wǎng)大小的SRTM DEM數(shù)據(jù)(Farr et al., 2007)被用于消除干涉相位中的地形影響.有關(guān)研究顯示該地形數(shù)據(jù)在垂直精度方面略優(yōu)于其他地形數(shù)據(jù)集(Li et al., 2021).TOPS數(shù)據(jù)處理過(guò)程中,方位向與距離向的多視比參數(shù)分別設(shè)置為2和8,并使用了高斯平滑和Goldstein方法對(duì)干涉圖進(jìn)行降噪處理.解纏中,選擇相干性0.15作為解纏閥值,該值以下像素做填零處理.考慮大氣在InSAR處理中存在的潛在影響,我們利用基于ERA5 外部大氣數(shù)據(jù)的InSAR大氣改正方法(Wang et al., 2021),進(jìn)行了大氣噪聲抑制操作,部分提升了衛(wèi)星視線(Line of Sight, LOS)向形變場(chǎng)精度.
受雷達(dá)成像幾何和斷層展布的影響,T158軌道的震后干涉圖像的震后形變不顯著.因此,本文中只處理了T158一對(duì)同震圖像用于震源參數(shù)反演(表1),震后形變分析部分主要依托于T165軌道的形變序列(圖2).在震后InSAR形變分析中,除基于外部大氣改正操作外,我們嘗試運(yùn)用干涉網(wǎng)絡(luò)分析的技術(shù),處理中當(dāng)某一時(shí)刻SAR數(shù)據(jù)與等時(shí)長(zhǎng)的前后數(shù)據(jù)進(jìn)行干涉,可以通過(guò)相減直接獲取該時(shí)刻大氣的近似(Tymofyeyeva and Fialko, 2015).然后,采用基于光滑最小二乘形變分離的SBAS方法最終得到震后形變序列(Schmidt and Bürgmann, 2003).震后形變趨勢(shì)的具體分析將在后文進(jìn)一步展開(kāi).
圖2 T165軌道干涉像對(duì)時(shí)空基線圖Fig.2 Temporal-spatial baseline plot of the processed interferometric pairs for T165
表1 用于地震斷層建模的Sentinel-1和ALOS2干涉像對(duì)的基本信息Table 1 The selected interferometric pairs from Sentinel-1 and ALOS2 SAR data used for the earthquake fault modelling
由于植被覆蓋稀疏,于田地震的同震形變場(chǎng)相位清晰連續(xù),自西向東分別展現(xiàn)兩個(gè)顯著形變中心(圖3).其中,S1降軌T165同震形變場(chǎng)顯示,東西兩盤形變趨勢(shì)相反,其中西盤表現(xiàn)為負(fù)形變,表示遠(yuǎn)離衛(wèi)星方向運(yùn)動(dòng),最大LOS位移達(dá)到200 mm,東盤符號(hào)為正,標(biāo)識(shí)朝向衛(wèi)星運(yùn)動(dòng),最大 LOS向位移達(dá)到80 mm(圖3d).相較而言,升軌T158的同震形變場(chǎng)在東盤表現(xiàn)出與降軌T165相反的趨勢(shì),同時(shí)表現(xiàn)為遠(yuǎn)離衛(wèi)星運(yùn)動(dòng),但幅度較小,LOS向最大位移僅為10 mm(圖3a).可見(jiàn),升降軌道D-InSAR結(jié)果同時(shí)顯示出于田地震的形變模式表現(xiàn)為西側(cè)形變區(qū)的變形幅度大于東側(cè).斷層的上盤形變一般較下盤更劇烈,我們可以定性判定西側(cè)形變區(qū)對(duì)應(yīng)此次發(fā)震斷層上盤,故而發(fā)震斷層可能西傾.同時(shí),考慮地表三維形變分量在InSAR升降軌形變貢獻(xiàn)的差異性,它們?cè)跂|側(cè)形變中心區(qū)表現(xiàn)出的極性相反的特征,說(shuō)明于田地震引起的形變具有一定的水平分量.另外,在干涉圖中可見(jiàn)一定范圍的失相干區(qū)域,可能與該地區(qū)的地表積雪有關(guān).其中ALOS2干涉圖(圖3g)與S1 T158同屬于升軌數(shù)據(jù),同震條紋基本一致,側(cè)面支撐了Sentinel-1數(shù)據(jù)結(jié)果的可靠性.因該軌道震后數(shù)據(jù)獲取于震后3個(gè)月后,該數(shù)據(jù)未用于后續(xù)反演分析.
圖3 反演模擬結(jié)果及殘差對(duì)比(a) S1-T158原始干涉圖; (b) S1-T158模擬結(jié)果; (c) S1-T158 InSAR觀測(cè)與模擬的殘差; (d,e,f) 對(duì)應(yīng)于ALOS2-T158軌道結(jié)果; (g,h,i) 對(duì)應(yīng)于S1-T165干涉相對(duì)結(jié)果. 黑色矩形為斷層在地表的投影,紅線表示斷層面與地表的交線.Fig.3 Comparison of observed and predicted interferograms and their residuals(a) Original InSAR observations of S1-T158; (b) Simulation results of S1-T158; (c) Residuals observed and simulated by S1-T158 InSAR; (d,e,f) are the corresponding results of ALOS2-T158; (g,h,i) are the corresponding results of S1-T165. The black rectangle is the projection of the fault plane on the surface, and the red line represents the intersection line between the fault plane and the surface.
為加快計(jì)算速度,同時(shí)壓制噪聲影響,我們首先對(duì)InSAR升降軌同震圖像進(jìn)行了四叉樹(shù)降采樣處理(Simons et al., 2002),以獲得有限參考點(diǎn)用于反演(表1).我們選擇了兩幅S1 InSAR結(jié)果作為地表約束(表1),采用兩步式的地球物理反演策略(Fukahata and Wright, 2008),即先借助非線性搜索算法確定斷層位置(矩形斷層的上邊界中點(diǎn))以及幾何參數(shù)(走向、傾角和滑動(dòng)角等),再進(jìn)一步開(kāi)展線性反演得到滑動(dòng)空間分布特征.反演工作借助測(cè)地學(xué)反演程序PSOKINV完成,該軟件采用多峰值顆粒群(MPSO)優(yōu)化算法,實(shí)現(xiàn)最優(yōu)擬合條件下的非線性參數(shù)搜索過(guò)程(馮萬(wàn)鵬和李振洪,2010).PSOKINV程序具有收斂效率高和控制參數(shù)少等特點(diǎn),目前已經(jīng)在多個(gè)地震實(shí)例研究中發(fā)揮作用(馮萬(wàn)鵬和李振洪,2010;李永生等,2015;Feng et al., 2018).MPSO算法只需要待解參數(shù)的范圍作為初始參數(shù),本文設(shè)置斷層長(zhǎng)和寬度范圍皆為2至20 km,滑動(dòng)角范圍為-45°至-135°,傾角范圍為-20°至-70°,震級(jí)大小范圍為MW6.25至MW6.35.反演中的正演模型皆以彈性半空間位錯(cuò)理論為基礎(chǔ)(Okada, 1985).
最終反演得到了基于InSAR形變結(jié)果的于田地震的最優(yōu)擬合解(表2),斷層面在地表投影可見(jiàn)圖3.反演結(jié)果表明,于田MW6.3地震發(fā)震斷層震中(矩心)位于82.45°E,35.61°N,斷層面寬約7 km,長(zhǎng)約14 km,近南北展布走向?yàn)?86°,傾角約為59°,滑動(dòng)角為-102° (圖4).考慮當(dāng)?shù)氐貧傂韵禂?shù)為1.0×1022Nm(Masek et al., 1994),可得到震級(jí)為MW6.25.該結(jié)果與USGS、GCMT等震源機(jī)制解以及其他InSAR 反演結(jié)果(He et al., 2020; Yu et al., 2021)近乎一致.
表2 2020年MW6.3于田地震震源參數(shù)Table 2 Source parameters of the 2020 MW6.3 Yutian earthquake
InSAR觀測(cè)的形變誤差(如殘余大氣擾動(dòng))可能會(huì)導(dǎo)致反演結(jié)果中存在不同程度的不確定性.為定量化震源參數(shù)的不確定性水平,我們開(kāi)展了基于反演結(jié)果的一系列正演模擬.考慮到走向和傾角是斷層幾何中非線性程度較高的兩維參數(shù),我們以非線性反演得到的最優(yōu)模型中的走向和傾角為輸入,給定一定范圍建立參數(shù)集合,并對(duì)離散后的每組斷層參數(shù)進(jìn)行了獨(dú)立的InSAR觀測(cè)的模擬,從而得到對(duì)應(yīng)所有參數(shù)集合的殘差分布.以InSAR觀測(cè)遠(yuǎn)場(chǎng)無(wú)形變區(qū)的標(biāo)準(zhǔn)偏差(~0.01 m)為閾值,可大體給出非線性參數(shù)的可信區(qū)間.該方案已經(jīng)成功應(yīng)用于2017年Iran地震的研究實(shí)踐(Feng et al., 2018).最終,我們基于本文InSAR觀測(cè)評(píng)估了斷層參數(shù)的不確定性范圍,其中走向的取值范圍為185°±2°,傾角的取值范圍為59°±3°.
為了理解地震發(fā)生過(guò)程中沿?cái)鄬拥幕瑒?dòng)分布細(xì)節(jié),我們進(jìn)而開(kāi)展了非均勻滑動(dòng)的線性反演.首先,基于上一階段中反演得到的斷層幾何參數(shù)和位置來(lái)固定斷層空間分布,將斷層面的長(zhǎng)度和寬度擴(kuò)展到20 km× 20 km,將斷層面離散成較小的滑動(dòng)單元.考慮到觀測(cè)資料在深度上的分辨能力,我們采用了深度相關(guān)的子斷層離散方案(Fialko, 2004).通過(guò)添加深度相關(guān)的阻尼因子,共劃分了311塊子斷層,空間上深度越淺滑塊越小,反之則滑塊越大(圖4).反演中允許一定滑動(dòng)自由度,將一組沿滑動(dòng)角為-45°至-135°的正交的滑動(dòng)矢量作為待解參數(shù).地表觀測(cè)與滑動(dòng)之間的數(shù)學(xué)表達(dá)為Gs=d,其中G為表征單一滑塊上單位滑動(dòng)矢量的格林矩陣,由Okada模型(Okada,1985)計(jì)算得來(lái),d為InSAR LOS向觀測(cè),s即為待解滑動(dòng)量.同時(shí),拉普拉斯光滑約束作為G矩陣的一部分,也被用于防止出現(xiàn)滑動(dòng)的異常突變.最終利用共軛梯度方法(Ward and Barrientos, 1986)解上述線性方程組,獲取同震滑動(dòng)分布.
于田地震同震滑動(dòng)分布結(jié)果(圖4)表明,地震同震滑動(dòng)的深度范圍為地下4~12 km,未出露到地表.同震的最大滑動(dòng)量和最大滑動(dòng)深度分別約為0.9 m及7 km.滑動(dòng)量大于0.4 m滑塊的平均滑動(dòng)角約為-99°,可見(jiàn)于田地震以正斷分量為主.
圖4 于田地震同震滑動(dòng)分布及精定位后的余震(冉慧敏等, 2020)分布圖Fig.4 Coseismic slip distribution of the 2020 Yutian earthquake and relocated aftershocks (Ran et al., 2020)
基于線性反演結(jié)果,我們開(kāi)展正演模擬,再現(xiàn)了InSAR觀測(cè)同震條紋(圖3).模擬結(jié)果(圖3)表明,由本文確定的斷層滑動(dòng)模型基本可以理論再現(xiàn)S1及ALOS2三軌InSAR觀測(cè)結(jié)果,理論模擬在形變幅度(條紋梯度)和空間分布特征方面均與InSAR同震觀測(cè)趨于一致.兩者殘差顯示,InSAR觀測(cè)結(jié)果與模擬結(jié)果的殘差(圖3c,f和i)整體均處于干涉圖的誤差水平,~2 mm.
為了探究于田地震震后余滑過(guò)程,我們選取了信噪比較高的T165軌道震后干涉圖像,其中第一幅覆蓋震后兩個(gè)月累計(jì)震后形變的干涉像對(duì)(表1)用于震后早期的余滑分析.對(duì)震后6個(gè)月內(nèi)的SAR影像進(jìn)行時(shí)序InSAR處理,以期恢復(fù)于田地震的震后形變趨勢(shì).于田地震震后升降軌的形變幅度差異明顯,表現(xiàn)與同震類似的形變特征(圖5a,b).因升降軌道LOS向位移符號(hào)相同,可以推斷震后形變?nèi)匀灰源瓜蜻\(yùn)動(dòng)為主,而幅度的顯著差異性,揭示了伴隨東西向水平運(yùn)動(dòng)的可能.利用于田地震的同震斷層幾何模型,我們進(jìn)而反演確定了早期余滑的空間分布特征.反演結(jié)果顯示,震后兩個(gè)月時(shí)間內(nèi)于田地震斷層上發(fā)生了較為明顯的震后余滑現(xiàn)象,最大滑動(dòng)量達(dá)0.16 m,滑動(dòng)主要分布在同震周邊(圖6).該現(xiàn)象與其他震后研究類似,表明震后滑動(dòng)主要是同震應(yīng)力驅(qū)動(dòng)的結(jié)果.于田地震同震破裂未延伸至地表,但分布一定淺部無(wú)震滑動(dòng),揭示了斷層淺部區(qū)域具有滑動(dòng)強(qiáng)化屬性,可能存在淺部滑動(dòng)欠滯的現(xiàn)象.
圖6 于田地震同震與震后滑動(dòng)分布紅星代表于田地震震中,色標(biāo)代表同震滑動(dòng)大小,綠線代表余震滑動(dòng)等值線,黑色箭頭代表同震滑動(dòng)方向.Fig.6 Coseismic and postseismic slip distributions of the Yutian earthquakeThe red star is the hypocenter of the Yutian earthquake. The color code represents the magnitude of the coseismic slip, while the green contour lines indicate the early afterslip captured by the postseismic InSAR observations. The black arrows show the directions of coseismic slip vectors on the fault.
以同震反演得到斷層參數(shù)確定的震后滑動(dòng)模型反演中存在顯著的殘差,無(wú)法重現(xiàn)沿剖面(圖7b)的最大震后形變.這里,我們獨(dú)立依托震后形變,類似于同震反演過(guò)程單獨(dú)確定適用震后形變擬合的最優(yōu)斷層幾何.相較于同震幾何參數(shù),震后反演結(jié)果的傾角更大,為75°(表3);滑動(dòng)深度更淺,為3 km.以此斷層模型為輸入,理論模擬與實(shí)際觀測(cè)的LOS向形變與模擬結(jié)果擬合較好(圖7b),斷層上盤持續(xù)下沉.此期間內(nèi),累計(jì)震后形變約為同震的10%(圖7b).根據(jù)同震、震后反演獲取的斷層幾何的特征,可能說(shuō)明了于田地震的實(shí)際發(fā)震斷層由地表向地下延伸時(shí),傾角逐漸變緩,具有鏟形斷層的特征(圖7d).盡管反演中數(shù)據(jù)中的潛在噪聲可能引入一定參數(shù)誤差,但鏟狀特征的存在更可能與實(shí)際斷層結(jié)構(gòu)相符合(Shelton et al., 1984).可見(jiàn),聯(lián)合同震和震后觀測(cè)可能是恢復(fù)鏟狀斷層結(jié)構(gòu)的有效途徑.
大量研究表明,余滑是控制震后形變的重要因素(Wang and Fialko, 2018; Feng et al., 2020; Jing et al., 2020),于田地震震后長(zhǎng)期發(fā)生的無(wú)震蠕滑對(duì)當(dāng)?shù)氐匦蔚乃茉煊兄鴿撛诘挠绊?為了更好地研究震區(qū)震后形變對(duì)裂谷生長(zhǎng)的影響,我們選取本次地震發(fā)震斷層上盤處形變量最大的點(diǎn)O(圖5a,b)
圖5 于田地震震后LOS向形變場(chǎng)(a) S1-T165降軌干涉圖:20200629_20200828; (b) S1-T158升軌干涉圖:20200704_20200821.紅線代表跨越震后余滑作用范圍的剖面BB′,三角形代表用于震后形變衰減曲線模擬的點(diǎn)O.Fig.5 Postseismic InSAR interferograms of the Yutian earthquake(a) S1-T165 descending interferogram of 20200629_20200828; (b) S1-T158 ascending interferogram of 20200704_20200821. Red line across the postseismic deformation area is a profile BB′ and the triangle, O is the selected site, at which the postseismic deformation time series is determined for analyzing postseismic history.
作為特征點(diǎn),對(duì)于田地震同震及震后的地表形變過(guò)程進(jìn)行分析.
震后6個(gè)月InSAR時(shí)序序列揭示了當(dāng)?shù)爻掷m(xù)形變?cè)鲩L(zhǎng)趨勢(shì)(圖7c).基于指數(shù)型余滑模型A(1-e-t/τ),其中A為常數(shù)項(xiàng),t為同震后歷經(jīng)的時(shí)間,τ為松弛時(shí)間.震后余滑被認(rèn)為是主震應(yīng)力驅(qū)動(dòng)的直接結(jié)果,松弛時(shí)間τ大小控制了震后形變的變形趨勢(shì),與斷層摩擦系數(shù)直接相關(guān)(Barbot et al., 2009; Gonzalez-Ortega et al., 2014).我們利用InSAR震后形變序列搜索最優(yōu)擬合參數(shù),嘗試恢復(fù)自同震時(shí)刻開(kāi)始的完整震后余滑衰減趨勢(shì).回歸分析中,采用非線性斷層反演中沿用的MPSO算法(馮萬(wàn)鵬等,2010),得到A和τ分別為-0.037和98.566.模擬結(jié)果表明(圖7c),于田地震的震后余滑可能持續(xù)了將近一年的時(shí)間,其中理論累計(jì)形變曲線中震后第4天至65天造成的總形變量與本文得到的InSAR觀測(cè)相當(dāng),~2 cm.由模型預(yù)測(cè)可知,地震震后余滑過(guò)程在震后的一年內(nèi)可能累積造成了將近3 cm的LOS向形變,達(dá)到了同震LOS向形變的15%.可見(jiàn)于田地震的震后無(wú)震蠕滑不僅在長(zhǎng)時(shí)間內(nèi)持續(xù)控制著裂谷的東西向擴(kuò)張,也對(duì)裂谷生長(zhǎng)及相關(guān)地貌的塑造起到了非常顯著的作用.參考前人對(duì)2008年當(dāng)雄地震震后研究結(jié)果(Bie et al., 2014),可見(jiàn)震后過(guò)程對(duì)裂谷生長(zhǎng)的控制作用具有一定的普遍性.
圖7 (a)沿BB′地形剖面;(b)2020于田地震同震及震后LOS向觀測(cè)和模擬的形變趨勢(shì);(c)由InSAR時(shí)序分析得到的2020于田地震震后累計(jì)最優(yōu)指數(shù)形變曲線;(d)于田地震鏟狀斷層模型Fig.7 (a) Topography along the profile BB′;(b) Observed, modelled co- and post-seismic LOS deformations across the 2020 Yutian earthquake fault; (c) Optimal postseismic accumulated deformation history of the 2020 Yutian earthquake retrieved from the InSAR time series analysis based on an exponential numerical model; (d) The listric fault model of the Yutian earthquake
青藏高原內(nèi)部一系列南北展布的裂谷與青藏高原的東西向伸展及生長(zhǎng)有著密切關(guān)聯(lián).地貌學(xué)證據(jù)表明,青藏高原南北裂谷形態(tài)差異十分巨大,發(fā)育狀況也不盡相同(Harrison et al., 1992).Kapp等(2008)對(duì)青藏高原不同地區(qū)裂谷的構(gòu)造特征進(jìn)行了分析,提出了一種與傾角相關(guān)的裂谷生長(zhǎng)模式.該模式下,裂谷起始以高角度正斷層的構(gòu)造形式啟動(dòng),在擴(kuò)張過(guò)程中,伴隨著拆離及均衡回彈,正斷層傾角逐漸變緩.這一模式提出正斷層的傾角大小可能指示了裂谷的生長(zhǎng)發(fā)育階段,即傾角大的裂谷可能處于裂谷發(fā)育的早期階段,而傾角小的裂谷已發(fā)育至較為成熟的階段.
我們收集了青藏高原不同地區(qū)基于大地測(cè)量學(xué)觀測(cè)的正斷層地震傾角數(shù)據(jù),結(jié)果(表3和圖8)顯示,青藏高原不同地區(qū)正斷事件的傾角有著較為顯著的差異,如藏北羌塘地區(qū)正斷事件的傾角明顯大于藏南地區(qū).從藏南過(guò)渡到藏北,地震斷層傾角表現(xiàn)出較為明顯的遞增趨勢(shì).根據(jù)青藏高原裂谷生長(zhǎng)階段與傾角的關(guān)系(圖8),可以推論,藏北羌塘地區(qū)的裂谷可能處于裂谷生長(zhǎng)的早期階段,而藏南拆離系裂谷的生長(zhǎng)發(fā)育可能已趨于成熟.
表3 青藏高原正斷層地震傾角信息Table 3 Dip angles of normal faulting earthquakes in the Tibetan Plateau collected in this study
有研究表明,藏北地區(qū)的羌塘地塊分布有許多與青藏高原東西向伸展構(gòu)造相關(guān)的埃達(dá)克巖,其低溫?zé)崮甏鷮W(xué)測(cè)年結(jié)果將藏中裂谷可能的啟動(dòng)時(shí)間限定在了47~38 Ma的始新世(Wang et al., 2010;Ou et al., 2017).而藏南裂谷的啟動(dòng)時(shí)間則普遍被相關(guān)測(cè)年結(jié)果限定在10 Ma左右(張佳偉等, 2020),這可能表明藏南裂谷的啟動(dòng)時(shí)間普遍晚于藏北裂谷.結(jié)合年代學(xué)定年結(jié)果及基于傾角大小的藏北藏南發(fā)育階段,可以推論啟動(dòng)時(shí)間較早的藏北羌塘地區(qū)裂谷的活動(dòng)擴(kuò)張速率可能遠(yuǎn)低于啟動(dòng)更晚的藏南裂谷.對(duì)青藏高原裂谷活動(dòng)速率的研究結(jié)果顯示(圖8),藏北的雙湖裂谷的活動(dòng)速率約為0.1~0.3 mm·a-1(李亞林等, 2001),藏中的溫泉盆地附近的斷裂滑動(dòng)速率約為0.2~0.3 mm·a-1(吳中海等, 2005);而藏南的亞?wèn)|裂谷的活動(dòng)速率則約為0.9 mm·a-1,藏南其他裂谷的活動(dòng)速率則更高,約為1~2 mm·a-1(Wang et al., 2020).藏北裂谷活動(dòng)速率隨著緯度的增加而降低,這一變化趨勢(shì)與傾角隨緯度的變化關(guān)系呈現(xiàn)出很好的反相關(guān)現(xiàn)象(圖8).同時(shí),GCMT地震目錄顯示,自1975年至今,青藏高原緯度大于31°的地區(qū)共發(fā)生了102起被記錄的正斷層事件;而高原內(nèi)緯度小于31°的藏南地區(qū)則共發(fā)生了151起被記錄的正斷層事件.這表明藏南地區(qū)的地震活動(dòng)性強(qiáng)于藏北地區(qū),可能與青藏高原南北裂谷活動(dòng)性的差異相關(guān).以上的研究結(jié)論皆與我們推論的藏南裂谷活動(dòng)擴(kuò)張速率可能遠(yuǎn)大于藏北的羌塘地區(qū)相一致.其中,圖8中謝通門地震斷層的傾角大小相較于其他地震顯得較小,考慮到謝通門地震相對(duì)震級(jí)也較小,這可能是因?yàn)樵摰卣鸢l(fā)生在特殊的構(gòu)造環(huán)境上(如次生構(gòu)造)造成的.
圖8 青藏高原自南向北正斷事件的傾角及裂谷活動(dòng)速率變化圖Fig.8 Diagram of dips of normal fault earthquakes vs rift extending rates from south to north Tibetan Plateau
本研究應(yīng)用干涉合成孔徑雷達(dá)技術(shù),厘定了2020年于田MW6.3地震的同震、早期震后位移場(chǎng),并進(jìn)而確定了震源參數(shù),同震以及震后空間滑動(dòng)分布,并利用震后6個(gè)月連續(xù)InSAR觀測(cè)嘗試恢復(fù)了于田地震的震后形變趨勢(shì).同時(shí)結(jié)合已有青藏高原正斷層地震的研究結(jié)果,討論了藏北羌塘地區(qū)裂谷帶的發(fā)育特征.此次工作中得到主要認(rèn)識(shí)如下:
(1)反演結(jié)果表明,于田地震破裂區(qū)域?qū)捈s7 km,長(zhǎng)約14 km,發(fā)震斷層近南北走向,傾角約為59°,滑動(dòng)角約為-102°.同震滑動(dòng)主要發(fā)生的范圍為地下4~12 km,最大滑動(dòng)量為0.9 m,最大滑動(dòng)深度約為7 km,同震破裂未出露至地表.地震矩震級(jí)大小為MW6.3,震中(矩心)位于 82.45°E,35.61°N.同震最大視線向沉降達(dá)20 cm.
(2)InSAR觀測(cè)表明,于田地震震后無(wú)震蠕滑的形變空間特征與同震類似,二者都持續(xù)貢獻(xiàn)裂谷的生長(zhǎng).聯(lián)合同震和震后斷層反演,此次地震的發(fā)震斷層可能具有鏟形特征.震后6個(gè)月連續(xù)的時(shí)序InSAR擬合出的余滑曲線表明,于田地震的震后余滑可能持續(xù)將近一年的時(shí)間,累積可以造成3 cm的LOS向形變.
(3)通過(guò)分析近年來(lái)青藏高原地區(qū)正斷層地震的傾角大小及裂谷擴(kuò)張速率與啟動(dòng)時(shí)間,得出了自裂谷啟動(dòng)以來(lái),藏北羌塘地區(qū)裂谷的擴(kuò)張生長(zhǎng)速率可能遠(yuǎn)小于藏南拆離系裂谷這一推論.
致謝特別感謝新疆地震局冉慧敏高工提供此次地震的重定位余震目錄.本研究采用的Sentinel-1數(shù)據(jù)下載自阿拉斯加衛(wèi)星中心(ASF, https:∥search.asf.alaska.edu/),ALOS2數(shù)據(jù)由日本宇航局(JAXA)通過(guò)RA6課題(3159)提供,大部分圖件使用GMT5.4.5(Wessel et al., 2013)繪制完成,在此一并表示感謝.