基于形变观测分析2011年日本9.0级地震与断层运动间关系

2011年3月11日日本发生9.0级地震,本文以此次地震的震间、同震和震后形变观测为约束,依据不同时段断层运动空间分布特征分析日本海沟地区强震与断层运动间关系.震间日本海沟地区,断层运动闭锁线深度约为60km,闭锁线以上从深到浅依次为断层运动强闭锁段、无震滑移段和弱闭锁段.由同震位错反演结果,2011年日本9.0级地震同震存在深浅两个滑移极值区,同震较浅的滑移极值区(同震位错量10~50m,深度小于30km)震间为断层弱闭锁段;同震较深的滑移极值区(同震位错量10~20m,深度在40km左右)震间为断层强闭锁段;而在两者之间的过渡带同震位错相对较小,震间断层运动表现为无震滑移.震后初期断层运动主要分布在在闭锁线以上的同震较深滑移极值区,而同震较浅的滑移极值区能量释放比较彻底,断层震后余滑量相对较小.依据本文同震和震间断层运动反演结果,震间强闭锁段积累10m同震位错需要100多年时间,与该区域历史上7级地震活动复发周期相当;震间弱闭锁段积累30~50m同震位错约需要300~600年时间,与相关研究给出的日本海沟9级左右地震复发周期比较一致.在实际孕震能力判定的工作中,由于不同性质的断层段在同震过程中会表现更多的组合形式,断层发震能力判定结果存在更多的不确定性,但利用区域形变观测等资料给出震间断层运动特征的研究工作对于断层强震发震能力的判定具有非常重要的实际意义.     

Coseismic Fault Slip of the September 16, 2015 <Emphasis Type="Italic">Mw</Emphasis> 8.3 Illapel,Chile Earthquake Estimated from InSAR Data

The complete surface deformation of 2015 Mw 8.3 Illapel, Chile earthquake is obtained using SAR interferograms obtained for descending and ascending Sentinel-1 orbits. We find that the Illapel event is predominantly thrust, as expected for an earthquake on the interface between the Nazca and South America plates, with a slight right-lateral strike slip component. The maximum thrust-slip and right-lateral strike slip reach 8.3 and 1.5 m, respectively, both located at a depth of 8 km, northwest to the epicenter. The total estimated seismic moment is 3.28 × 10²¹ N.m, corresponding to a moment magnitude Mw 8.27. In our model, the rupture breaks all the way up to the sea-floor at the trench, which is consistent with the destructive tsunami following the earthquake. We also find the slip distribution correlates closely with previous estimates of interseismic locking distribution. We argue that positive coulomb stress changes caused by the Illapel earthquake may favor earthquakes on the extensional faults in this area. Finally, based on our inferred coseismic slip model and coulomb stress calculation, we envision that the subduction interface that last slipped in the 1922 Mw 8.4 Vallenar earthquake might be near the upper end of its seismic quiescence, and the earthquake potential in this region is urgent.     

联合GNSS和InSAR观测位移反演2008年汶川大地震断层位错模型参数

利用现代空间大地测量技术,尤其是卫星合成孔径雷达干涉测量,能够获取高精度、高空间分辨率的同震和孕震形变,为地震断层形变和破裂机制研究提供了前所未有的机遇。本文介绍了利用大地测量观测数据反演地震断层位错模型参数的贝叶斯反演方法。联合运用2008汶川大地震前后GNSS和InSAR技术观测获得的同震位移,反演了地震断层的几何参数和滑动位错分布。研究结果表明,汶川地震的断层滑动主要集中在倾角较陡的浅部,同时包含逆冲和右旋走滑,其中最大逆冲6.1m,最大右旋6.5m。根据断层滑动分布正演计算得到的上盘同震位移明显小于下盘,预示该断层两侧孕震形变可能存在较大的不对称性。     

Analysis of Coseismic Fault Slip Models of the 2012 Indian Ocean Earthquake: Importance of GPS Data for Crustal Deformation Studies

Based on continuous GPS data, we analyze coseismic deformation due to the 2012 Indian Ocean earthquake. We use the available coseismic slip models of the 2012 earthquake, derived from geodetic and/or seismic waveform inversion, to calculate the coseismic displacements in the Andaman-Nicobar, Sumatra and Java. In our analysis, we employ a spherical, layered model of the Earth and we find that Java Island experienced coseismic displacements up to 8 mm, as also observed by our GPS network. Compared to coseismic offsets measured from GPS data, a coseismic slip model derived from multiple observations produced better results than a model based on a single type of observation.     

近场位移数据约束的2013年芦山地震破裂模型及其构造意义

以往的研究显示了2013年芦山MS7.0级地震发震断层的隐伏逆冲断层基本特征,但是破裂深部细节差异较大.本文以近场密集的同震形变数据约束芦山地震破裂面几何形状及滑动分布,结果显示芦山地震破裂面具有铲状结构,上部16km为43°~50°高角度断层,深部16~25km为小于27°的低角度断层,破裂深度与重定位的余震分布深度一致.破裂分布模型清楚显示上下两个断层上各有一个滑动幅度大于0.5m的峰值破裂区,最大滑动量1.5m位于13km深处.重定位的余震分布基本都落在最大滑动量等值线外部库仑应力增加的区域.芦山地震破裂面几何形状和滑动分布特征与2008年汶川MS8.0级地震映秀—北川破裂相似,支持龙门山冲断带发育大规模的近水平滑脱层,是青藏高原东缘地壳缩短增厚、龙门山挤压隆升的重要证据.     

龙门山断层地震周期及其动力学过程模拟研究

在断层面上引入速率-状态相依摩擦本构关系、考虑铲形逆冲断层几何结构特征、断层下盘和上盘中下地壳及上地幔为黏弹性介质、上盘上地壳为弹塑性介质,本文用二维有限元动力学模型模拟了龙门山断层上大震准周期复发行为、分析了断层上地震孕育位置、地震周期不同阶段的应力/应变场演化特征.不同于近垂直走滑断层上的地震周期行为,大陆铲形逆冲断层上的构造应力的积累和释放过程更复杂、有其独特性.我们得到如下认识：（1）铲形逆冲断层上的地震复发是准周期行为.（2）龙门山断层最大库仑应力位于断层17~20 km深处,应力长期积累和同震释放都在此深度最大,说明地震会在此处孕育、发动.（3）在断层破裂的深部和浅部,同震滑动大小和构造应力释放大小并非同步,而是差异悬殊.（4）地震仅部分释放区域积累的应变能,断层上盘上地壳顶部和底部的褶皱、破裂等永久变形形式也是释放应变能的重要形式.（5）应变能密度增量的演化图像分为：震间、同震、震后期,清晰反应了龙门山断层附近的地震动力学过程.（6）地震发生除释放能量外,同时也对近断层的中地壳和断层底部有很大的应变能加载;这些加载,在震后期可能通过震后滑移、余震或中下地壳乃至上地幔的驰豫形变用几十年时间释放.以上对大陆内铲形逆冲断层上变形特征的了解,有助于我们在其地震周期行为中评估地震危险性.     

Interplate Earthquake Fault Slip During Periodic Earthquake Cycles in a Viscoelastic Medium at a Subduction Zone

--A 2-D finite-element-method (FEM) numerical experiment of earthquake cycles at a subduction zone is performed to investigate the effect of viscoelasticity of the earth on great interplate earthquake fault slip. We construct a 2-D viscoelastic FEM model of northeast Japan, which consists of an elastic upper crust and a viscoelastic mantle wedge under gravitation overlying the subducting elastic Pacific plate. Instead of the dislocation model prescribing an amount of slip on a plate interface, we define an earthquake cycle, in which the plate interface down to a depth is locked during an interseismic period and unlocked during coseismic and postseismic periods by changing the friction on the boundary with the master-slave method. This earthquake cycle with steady plate subduction is periodically repeated to calculate the resultant earthquake fault slip.¶As simulated in a previous study (Wang, 1995), the amount of fault slip at the first earthquake cycle is smaller than the total relative plate motion. This small amount of fault slip in the viscoelastic medium was considered to be one factor explaining the small seismic coupling observed at several subduction zones. Our simulation, however, shows that the fault slip grows with an increasing number of repeated earthquake cycles and reaches an amount comparable to the total relative plate motion after more than ten earthquake cycles. This new finding indicates that the viscoelasticity of the earth is not the main factor in explaining the observed small seismic coupling. In comparison with a simple one-degree-of-freedom experiment, we demonstrate that the increase of the fault slip occurs in the transient state from the relaxed initial state to the stressed equilibrium state due to the intermittent plate loading in a viscoelastic medium.     

InSAR约束下的2008年汶川地震同震和震后形变分析

2008年5月12日,青藏高原东缘的龙门山断裂带上发生了M_w7.9级汶川地震.本文通过分析覆盖汶川地震震中区域的ALOS/PALSAR像对的方位向偏移量来选择无明显电离层扰动影响的像对进行干涉处理,获取了高精度、连续的InSAR地表形变场.在此基础上,结合高精度GPS同震形变数据,采用同震、黏弹性松弛震后形变联合反演模型同时确定了汶川地震的同震滑动分布和龙门山地区的流变结构参数.研究结果表明,汶川地震是一个断层破裂非常复杂的地震事件,其中,北川段、岳家山段、虹口段和汉旺段的滑动以逆冲为主,而青川段以右旋走滑为主.滑动主要发生在10 km深度以上的区域,最大滑动量位于虹口段的东北端,达10.7 m.地震释放的总能量为9.28×10²⁰ N·m（M_w7.91）,与地震学的结果一致.联合反演模型确定的龙门山地区中下地壳的黏性系数为2×10¹⁸ Pa·s,为青藏高原东部地区的黏性系数提供了一个可靠的下限值.如果有更长时间的震后形变观测时间序列,将为该区域提供更为可靠的流变结构.     

GNSS观测的九寨沟7.0级地震同震位移初步结果

2017年8月8日四川省九寨沟县发生了7.0级地震,中国大陆构造环境网络与北斗地基增强系统的GNSS连续观测共同监测到了此次地震的同震位移（坐标：东向为正,北向为正）,结果显示：3个站点记录到了明显的同震位移,距离震中43 km的九寨沟台站（SCJZ）在东西向的位移为-9.8±1.5 mm,在南北向的位移为3.3±0.7 mm;距离震中65 km的松潘站（SCSP）在东西向的位移为-1.8±0.7 mm,在南北向的位移为-7.7±0.6 mm;距离震中77 km的舟曲站（GSZQ）在东西向的位移为0.4±1.2 mm,在南北向的位移为3.6±0.8 mm.通过同震位移分布特征,可以推测此次地震为一次左旋走滑型事件,引起水平向同震位移大致不超过150 km范围,地震对东南侧的龙门山断裂带影响非常小,对北侧的塔藏断裂和西侧的岷江断裂处引起的同震位移为厘米级.同震位移的反演结果显示：断层面上滑动量主要集中在7 km深度,最大量值约为0.4 m,平均滑动角为-15°,利用滑动分布计算的相应矩震级为M_W6.4,与地震波反演结果相当.结合同震滑动分布、同震主应变分布、余震分布和震源机制解等特征,推测此次地震破裂极值区累积的能量得到较充分释放,进一步分析得出此次地震在塔藏断裂、岷江断裂和虎牙断裂处产生了一定的应力变化,值得持续关注.     

基于GNSS与InSAR约束的九寨沟MS7.0地震滑动模型及其库仑应力研究

2017年8月8日的九寨沟M_S7.0地震发生在岷江断裂、塔藏断裂及虎牙断裂交汇地区,地处青藏高原东北部的川甘交界地区,位于巴颜喀拉地块的东缘,地质构造复杂,对于九寨沟地震震中位置和发震断层的确定,存在不同意见.本文利用GNSS及升降轨InSAR观测,在获取九寨沟地震同震形变场的基础上,基于均匀弹性半无限位错模型,联合反演了发震断层的滑动分布模型,并计算了同震库仑应力变化.InSAR同震形变场显示,视线向最大沉降量和抬升量分别为0.21 m和0.16 m,形变场长轴为NW向,形变主要集中在断层西侧.距震中40 km和65 km的九寨和松潘两县,水平向的GNSS同震位移分别达14.31 mm和8.22 mm.联合GNSS和InSAR同震形变场反演得到的滑动分布主要集中在沿走向5~33 km,倾向2~20 km的范围内,平均滑动量为0.18 m,最大滑动量为0.91 m.发震断层长40 km,宽30 km,走向155°,倾角81°,滑动角-9.56°.同震位移场及滑移分布模型表明此次地震为一次左旋走滑为主的地震事件,地震破裂并未完全到达地表,与虎牙断裂北段的几何产状和运动学性质更为接近,结合精定位余震的分布,我们确定虎牙断裂北段为此次地震的发震断层,震中位于北纬33.25°,东经103.82°,震源深度10.86 km,矩震量为7.754×10¹⁸ Nm,相应的矩震级为M_W6.5,与美国地调局和哈佛大学给出的震源机制解基本一致.同震库仑应力导致了虎牙断裂北段延长线的东北和西南两端应力增强,其中塔藏断裂的罗叉段和马磨段未来强震的危险性值得关注.     

GPS观测到的智利及其邻区多次大地震前后地壳运动

2001年以来,智利及其邻区发生了3次7.1—7.7级和4次8级以上大地震,其震中附近至少有一个甚至多个GPS连续观测站观测到地震前后的地壳运动。从http://geodesy.unr.edu网站可获得GeoffreyBlewitt教授用GIPSY软件处理得到的南美大量GPS连续观测站南美板块（SA）区域参考框架位移时间序列,获得的这些大地震同震、震前位移积累和震后位移,特别是同震水平位移和震前水平位移积累,为探索地震预测,增添了更多有意义的震例。研究表明,这些大地震的同震水平位移也是震前水平位移积累的回跳或弹性回跳,同样也证明了震前存在前兆地壳形变;这些大地震前震中及其附近也无明显的垂直位移积累,由此证明了板块运动或地壳水平运动就是地震成因。尽管东日本和智利近海大地震的构造环境不同,日本2011年9级和智利2010年8.8级巨大地震前的地壳运动都清楚显示太平洋海底扩张。这些地震的同震水平位移回跳或弹性回跳的规律一致,地震成因都是水平挤压。智利多次大地震GPS观测到的最特殊现象是,在2015年8.3级地震震中以北,2010年8.8级地震的同震水平位移量值偏小,且方向异常一致向北,可认为是8.3级地震的前兆形变现象。临近智利的南美地区应是全球最利于地震预测探索的地区之一。      

九寨沟M_s7.0地震的InSAR观测及发震构造分析

2017年8月8日四川省九寨沟县发生M_s7.0地震.本文基于Sentinel-1 SAR影像,利用InSAR技术获取了此次地震的同震形变场,反演获得同震滑动分布,计算了同震位错对余震分布和周边断层的静态库仑应力变化,并对发震构造进行了分析讨论.结果表明:①InSAR同震形变场显示,九寨沟地震造成地表形变最大量级约为20 cm(雷达视线方向),同震形变存在非对称性分布特征.②同震位错以左旋走滑为主,主要发生在4~16 km深度,最大滑动量约为77 cm,位于9 km深处.反演得到的矩震级为Mw6.46.同震错动未破裂到地表.③大部分余震发生在库仑应力增加区.此次地震增加了震中周边地区一些断裂的库仑应力,如东昆仑断裂带东段、龙日坝断裂、虎牙断裂等.④东昆仑断裂东段的未来地震危险性值得关注.⑤九寨沟地震的发震断层为树正断裂,可能是虎牙断裂的北西延伸隐伏部分,此次地震是巴颜喀拉块体南东向运动受到华南块体的强烈阻挡过程中发生的一次典型构造事件.     

Interplate Coupling in the Kanto District, Central Japan, and the Boso Peninsula Silent Earthquake in May 1996

— I studied crustal deformation in the Kanto district, central Japan, based on continuous GPS data. Horizontal as well as vertical displacement rate demonstrate significant interaction between the landward Kanto block and the Philippine Sea plate. Although the subduction effect of the Pacific plate is not apparent, it is reasonable to consider the entire Kanto district is displaced westward due to the interaction with the Pacific plate. The GPS velocity data were inverted to estimate the slip deficit distribution on the Sagami Trough subduction zone. The result delineates a strongly coupled region on the plate interface, part of which corresponds to the 1923 Kanto earthquake. The strongly coupled region is located shallower than 20 km. In addition, the plate interaction is laterally heterogeneous even in the same depth range, implying thermal structure is not the only factor controlling interplate coupling. The GPS data also detected a silent earthquake event on the interface of the Philippine Sea slab east of the Boso Peninsula in the middle of May, 1996. The silent rupture propagated over a 50 km * 50 km wide area during about a week. The maximum slip was approximately 50 mm and the released seismic moment was 4.7*10¹⁸Nm (M _w 6.4). There was a small seismicity triggered by this silent event. The silent slip was located in the peripheral of the strongly coupled area, suggesting that frictional properties and/or stress conditions are inhomogeneous on the plate boundary interface.     

南美西海岸(智利)地震科学研究回顾及其对地震监测预报的启示

智利俯冲带是全球地震活动最为活跃的地区之一,2010年以来已经发生了3次8级强震。本文在调研了该地区地震地质、历史地震、运动学特征等方面研究成果的基础上,梳理了其中与地震预测预报有关的观测现象,总结了相关科学问题的研究进展和认识过程。结果显示:(1)智利俯冲带历史地震记录丰富,基于现代大地测量的断层震间闭锁研究深入,近年来多次8级地震均发生在历史地震破裂空段和滑动亏损较大的强闭锁区;(2) 2010年Maule M_W8.8地震和2014年Iquique M_W8.1地震震前均有不同程度的前兆变化,Maule M_W8.8地震前主要是"一类空区"和"二类空区"的嵌套;而Iquique M_W8.1地震前还伴有b值持续降低,前震活动频繁,GPS时间序列瞬态变化等一系列短临异常变化;(3)同震位错和震间闭锁的空间分布显示智利俯冲带断层摩擦性质复杂;(4)孕震物理力学机制复杂,不同学者不同方法之间还存在一些认识争议。     

5·12汶川地震地表破裂基本参数的再论证及其构造内涵分析

地震地表破裂基本参数是反演地震破裂过程的基本约束条件和预测其他活动断层地震危险性不可缺少的物理量.以野外地表破裂带重要观测点全站仪或差分GPS仪实测数据为基础,结合高分辨率遥感资料解译、先存断层陡坎构造地貌标志的识别、以及地形测绘资料的考证等,重新论证了5·12汶川地震地表破裂带展布样式、长度、最大同震位移值等基本参数.结果表明,地震地表破裂带长度可达240 km,最大垂直位移为6.5±0.5 m,最大右旋走滑位移4.9 m,基于倾角向下变缓逆断层模型推测汶川地震在龙门山推覆构造带中段产生了最大~7 m的地壳缩短量,说明青藏高原东缘横向逆断层为将高原内部东向水平运动转换为高原隆升的转换构造,这一研究结果有助于深化认识青藏高原东缘隆升机理.     

Frequency-dependent rupture process,stress change,and seismogenic mechanism of the 25 April 2015 Nepal Gorkha <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> 7.8 earthquake

On 25 April 2015, an M _w 7.8 earthquake occurred on the Main Himalaya Thrust fault with a dip angle of ~ 7° about 77 km northwest of Kathmandu, Nepal. This Nepal Gorkha event is the largest one on the Himalayan thrust belt since 1950. Here we use the compressive sensing method in the frequency domain to track the seismic radiation and rupture process of this event using teleseismic P waves recorded by array stations in North America. We also compute the distribution of static shear stress changes on the fault plane from a coseismic slip model. Our results indicate a dominant east-southeastward unilateral rupture process from the epicenter with an average rupture speed of ~3 km s^?1. Coseismic radiation of this earthquake shows clear frequency-dependent features. The lower frequency (0.05–0.3 Hz) radiation mainly originates from large coseismic slip regions with negative coseismic shear stress changes. In comparison, higher frequency (0.3–0.6 Hz) radiation appears to be from the down-dip part around the margin of large slip areas, which has been loaded and presents positive coseismic shear stress changes. We propose an asperity model to interpret this Nepal earthquake sequence and compare the frequency-dependent coseismic radiation with that in subduction zones. Such frequency-dependent radiation indicates the depth-varying frictional properties on the plate interface of the Nepal section in the main Himalaya thrust system, similar to previous findings in oceanic subduction zones. Our findings provide further evidence of the spatial correlation between changes of static stress status on the fault plane and the observed frequency-dependent coseismic radiation during large earthquakes. Our results show that the frequency-dependent coseismic radiation is not only found for megathrust earthquakes in the oceanic subduction environment, but also holds true for thrust events in the continental collision zone.     

2014年美国加州纳帕M_W6.1地震断层参数的Sentinel-1A InSAR反演

2014年8月24日,在美国加州旧金山海湾北部的纳帕地区发生了MW6.1地震.发震断层是西纳帕断裂系统中的一部分,但是该断层之前并未被足够重视.本文利用欧洲空间局最近发射成功并刚刚投入使用的Sentinel-1A卫星获取的第一对同震干涉像对(20140807-20140831),得到了该地震的地表同震形变场,结合震后24h内区域GPS同震形变资料作为约束条件,反演了纳帕地震的断层几何参数以及滑动分布.Sentinel-1A干涉结果表明,此次地震造成了明显的地面形变,视线向最大抬升和最大沉降量均达到了10cm.联合反演结果表明,该发震断层的走向为344°,倾角为80°.主要破裂以右旋走滑为主,平均倾滑角为-146.5°,最大倾滑量达到了1.1m,位于地表下约4km,存在明显的滑动亏损现象.此次地震,累计释放地震矩达1.5×1018 N·m,约合矩震级MW6.1.该结果略小于InSAR单独约束结果,可能与Sentinel-1A像对中包含的快速震后形变分量有关.     

利用GPS观测资料分析2008年于田Ms7.3地震的同震位移及震后形变

2008年3月21日新疆于田发生Ms7.3级地震.本文通过处理、分析GPS数据,得到破裂断层北侧100 km附近的同震位移及震后形变信息.在观测区域GPS点监测到10 mm左右的同震位移,其中最大为南向14 mm,东向5 mm.同震位移呈现一致性的东南向运动特征,证实于田地震存在显著的左旋走滑分量.震后台站向西南方向运...     

InSAR数据约束的2021年5月21日云南漾濞MS6.4地震发震构造研究

运用Sentinel-1A卫星数据和D-InSAR技术,获取2021-05-21云南漾濞M_S6.4地震的同震形变场。结果显示,漾濞地震同震形变场长轴近NW展布升降轨形变场符号相反,视线向最大沉降量和抬升量为0.1 m。InSAR同震形变场反演的滑动分布主要集中在沿走向2～12 km,倾向1～9 km的范围内,最大滑动量0.35 m,发震断层长9.8 km、宽4 km,滑动量主要集中在地下3～6 km范围内,滑动角-146.7°。同震位移场及滑动分布模型反映本次地震为发震断层的右旋走滑事件,地震破裂未达到地表。断层模型反演结果显示,矩震级为M_W6.1,发震断层以北西走向右旋走滑运动为主,初步认为本次M_W6.1地震发震断裂可能是一条NW向的维西—乔后断裂西侧的隐伏次生断裂。     

2008年汶川8.0级地震前横跨龙门山断裂带的震间形变

利用区域GPS和水准测量资料,结合地震构造背景的分析,本文研究2008年汶川8.0级地震前横跨龙门山断裂带地区的震间地壳形变,探讨引起这种形变的活动构造与动力学模式,并由此认识汶川地震的孕育与成因机制.主要结果表明：1997～2007年期间,自龙门山断裂带中段朝北西约230 km的地带内存在垂直于断裂的水平缩短变形、以及平行于断裂的水平右旋剪切变形,缩短率为1.3×10^-8/a （即：0.013 mm/km/a）,角变形速率为2.6×10^-8/a;同一地带在1975～1997年期间还表现出垂直上隆变形,上隆速率在龙门山前山断裂与中央断裂之间仅0.6 mm/a,而至龙门山后山断裂及其以西达2～3 mm/a.这些反映了在汶川地震之前至少10～30余年,龙门山断裂带中段的前山与中央断裂业已闭锁、并伴有应变积累.造成这种形变的主要原因是：以壳内的低速层为“解耦”带,巴颜喀拉地块上地壳朝南东的水平运动在四川盆地西缘受到华南地块的阻挡、转换成龙门山断裂带中段的逆冲运动;由于该断裂段的震间闭锁,致使西侧的巴颜喀拉地块的上地壳发生横向缩短以及平行断裂的右旋剪切变形.然而,龙门山断裂带北段在1997～2007年期间除了有大约0.9 mm/a的右旋剪切变形外,横向的缩短变形极微弱,这可能与该断裂段西侧的岷江、虎牙、龙日坝等断裂带吸收了巴颜喀拉地块朝东水平运动的大部分有关.另外,汶川地震前,横跨龙门山断裂带中段与北段的地壳形变特征的差异,与汶川地震时能量释放的空间分布吻合.