汶川地震前龙门山断裂带闭锁程度和滑动亏损分布研究

利用1999—2007期GPS水平速度场数据,采用Defnode负位错反演程序估算了龙门山断裂在汶川地震前的闭锁程度和滑动亏损分布,结合龙门山断裂带附近地表水平应变率场结果,综合分析了震前地壳变形特征.反演结果表明,震前龙门山断裂中北段处于完全闭锁状态,闭锁深度达到21 km(闭锁比例0.99)左右,垂直断层方向的挤压滑动亏损速率约为2.2 mm/a,平行断层方向的右旋滑动亏损速率约为4.6 mm/a.龙门山断裂南段只有地表以下12 km闭锁程度较高(闭锁比例0.99),垂直断层方向滑动亏损速率约为1.4 mm/a,平行断层方向滑动亏损速率约为4.6 mm/a;在12~16 km处闭锁比例约为0.83,垂直断层方向滑动亏损速率约为1.2 mm/a,平行断层方向滑动亏损速率约为3.8 mm/a;在16~21 km处闭锁比例约为0.75,垂直断层方向滑动亏损速率约为1.1 mm/a,平行断层方向滑动亏损速率约为3.5 mm/a.在21~24 km处整条断裂均逐步转变为蠕滑.上述反演结果与区域应变计算获得的龙门山断裂带中北段整体应变积累速率较低、南段应变积累速率较高相一致,均表明中北段闭锁程度高、南段闭锁程度稍低,该特征可以较好地解释汶川地震时从震中向北东向单向破裂现象.     

汶川地震前龙门山断裂带闭锁程度和滑动亏损分布研究

利用1999—2007期GPS水平速度场数据,采用Defnode负位错反演程序估算了龙门山断裂在汶川地震前的闭锁程度和滑动亏损分布,结合龙门山断裂带附近地表水平应变率场结果,综合分析了震前地壳变形特征.反演结果表明,震前龙门山断裂中北段处于完全闭锁状态,闭锁深度达到21 km(闭锁比例0.99)左右,垂直断层方向的挤压滑动亏损速率约为2.2 mm/a,平行断层方向的右旋滑动亏损速率约为4.6 mm/a.龙门山断裂南段只有地表以下12 km闭锁程度较高(闭锁比例0.99),垂直断层方向滑动亏损速率约为1.4 mm/a,平行断层方向滑动亏损速率约为4.6 mm/a;在12~16 km处闭锁比例约为0.83,垂直断层方向滑动亏损速率约为1.2 mm/a,平行断层方向滑动亏损速率约为3.8 mm/a;在16~21 km处闭锁比例约为0.75,垂直断层方向滑动亏损速率约为1.1 mm/a,平行断层方向滑动亏损速率约为3.5 mm/a.在21~24 km处整条断裂均逐步转变为蠕滑.上述反演结果与区域应变计算获得的龙门山断裂带中北段整体应变积累速率较低、南段应变积累速率较高相一致,均表明中北段闭锁程度高、南段闭锁程度稍低,该特征可以较好地解释汶川地震时从震中向北东向单向破裂现象.     

鲜水河断裂带炉霍段的震后滑动与形变

1973年2月在鲜水河断裂带炉霍段发生了M7.6地震破裂.自那以来,先后在炉霍县虾拉沱布设了若干横跨该地震断层(1973年破裂带)的地壳形变观测系统,包括断层近场的短基线、短水准、蠕变仪、人工构筑物等,以及断层近-远场的GPS观测站.利用这些观测系统的长期观测资料,本文分析了鲜水河断裂带炉霍段的震后滑动/变形及其时、空变化特征,并建立起解释这些特征的动力学模式.研究表明:(1)1973年地震后的头5年,地震断层在虾拉沱场地表现为开放性质,近场的断层震后滑动以无震左旋蠕滑为主,速率达到10.27 mm/a,且伴有微量的拉张性蠕动作用;1979年以来,左旋蠕滑速率由5.3 mm/a逐渐减小到2.27 mm/a,减小的过程呈对数函数型,反映此阶段断层面已逐渐重新耦合、正朝闭锁的方向发展,并伴有部分应变积累.(2)1999年以来,地震断层两侧远场的相对左旋位移/变形速率为10 mm/a,远大于同时期断层近场(跨距40~144 m)的左旋蠕滑速率0.66~2.52 mm/a;远-近场位移/形变速率的显著变化发生在地震断层两侧各宽约30 km的范围,显示出这是与大地震应力应变积累—释放相关的断裂带宽度.(3)结合动力学背景与深部构造信息,本文对这里断层的震后位移/变形及其时、空变化的机理进行初步解释,要点是:震后约5年之后,由于逐渐增大的断层滑动/摩擦阻抗,上地壳脆性层中的断层面由震后初期的开放性质逐渐转向重新耦合、并朝闭锁的方向发展,但其两侧地块深部持续的延性相对运动拖拽着浅部脆性层发生相应的弹性位移/变形.(4)可估计再经历15~25年,研究断裂段将完全闭锁,即进入积累下一次大地震应力应变的震间闭锁阶段.     

Deformable backstop as seaward end of coseismic slip in the Nankai Trough seismogenic zone

Wide-angle seismic surveys performed in the last decade have clarified the 3-D crustal structure along the Nankai Trough. The geometry and velocity structure of the southwestern Japan subduction zone are now well constrained. Comparing these observations with the rupture distribution of historic great thrust earthquakes, it appears that the coseismic rupture occurred along plate boundaries deeper than the wedge/backstop boundary (the boundary between the Neogene-Quaternary accretionary wedge and the crust forming the backstop). From the view of spatial relationship, both rupture distributions of the last two large events and the crust forming the backstop are considerably retreated from the trough axis in the west and east off the Kii Peninsula. In both areas, seamount or ridge subduction is apparent in seismic results, geomorphological data and geomagnetic data. The landward indentation of the deformable backstop, which corresponds to the crustal block of old accreted sediments, may be formed by seamount subduction according to published results of sandbox modeling. In particular, the subducted seamount may be a structural factor affecting the recession of the crustal block forming the backstop.     

川滇菱形块体东边界断层闭锁程度与滑动亏损动态特征研究

利用1999—2007期和2009—2013期中国大陆GPS速度场数据,采用DEFNODE负位错反演程序估算了川滇菱形块体东边界——鲜水河—安宁河—则木河—小江断裂带在汶川地震前后的断层闭锁程度和滑动亏损空间分布动态变化特征,讨论了汶川地震对该断裂系统的影响范围和程度,并结合b值空间分布和地震破裂时-空结果分析了断裂系统的强震危险段.结果表明,汶川地震前鲜水河断裂最南端为完全闭锁(闭锁深度25 km),中南段地表以下10~15 km深度为强闭锁状态,中北段基本处于蠕滑状态;安宁河断裂最南端闭锁很弱,其余位置闭锁深度为10~15 km;则木河断裂除最南端闭锁较弱以外,其余位置基本为完全闭锁;小江断裂在巧家以南、东川以南、宜良附近、华宁以北等四处位置闭锁较弱,其余位置为强闭锁.10年尺度的GPS速度场反演所得断层闭锁程度所指示的强震危险段,主要为鲜水河断裂道孚—八美段、安宁河断裂中段、则木河断裂中北段、小江断裂北段东川附近、小江断裂南段华宁—建水段,该结果与地质尺度的断层地震空区和30年尺度的b值空间分布所指示的危险段落具有一致性.汶川地震后断裂带远、近场速度分布和块体运动状态发生变化,这种区域地壳运动调整使得负位错模型反演得到的断裂带闭锁情况发生一定变化.汶川地震前后川滇菱形块体东边界平行断层滑动亏损速率均为左旋走滑亏损,且在安宁河断裂北端、则木河断裂中北段滑动亏损速率最大;除鲜水河断裂中南段与最南端和小江断裂东川附近以外,其余断裂震后滑动亏损速率均有所增加.垂直断层滑动亏损速率既有拉张亏损也有挤压亏损,且鲜水河断裂最南端由震前挤压转变为震后拉张,其余断裂除了安宁河断裂和小江断裂中段与最北端存在挤压滑动亏损速率外均为拉张速率.     

Development of dynamic slip model for a shear crack

<正>The objective of this paper is to develop a dynamic slip model for a shear crack under constant stress drop.This crack problem is formulated by a traction boundary integral equation(BIE) in the frequency domain and then solved by the hyper-singular boundary element method as well as the regularization technique proposed in this paper.Based on the spectral integral form of the kernel function,the unbounded term can be isolated and extracted from the hyper-singular kernel function by using the method of subtracted and added back in wave number domain.Finally,based on the inverse transformation from the frequency domain to the time domain,the time histories of crack opening displacement under constant stress drop can be determined.Three rupture models(simultaneous rupture model,symmetric bilaterally-propagating model and unilaterally propagating model) with specified time histories of stress drop are considered in this paper.Even though these three models will cause the same final slip shapes because of the same constant stress drop,the associated slip time functions differ significantly from each other during the rupture process.     

由重复微震估算的龙门山断裂带深部滑动速率

使用断层附近两个数字地震台网记录的重复地震数据研究了发生2008年汶川M7.9灾难性地震的龙门山断裂带的深部滑动速率。通过波形互相关分析,识别出波形高度相似的231组相似地震对和224组多重相似对。大部分相似地震活动表现为非周期性,复发间隔从数分钟到数百天不等。针对每一组多重相似对,我们首先依据在同一地震台站记录的波形来塑造参考地震的波形,用以量取每个地震相对于参考地震的S与P波的相对到时差,并确定每个地震与多重相似对质心的相对距离,从而挑选出每组地震中的重复事件。我们识别出了位于M7.9汶川地震震源区的12组重复微震,大部分重复微震分布在2008汶川地震同震大破裂闭锁区的边缘,揭示了微震活动与未来破坏性大震在空间位置上的密切关系。由重复微震直接获取的滑动速率随深度增加而增大,在4～18km深度范围滑动速率为3.5～9.6mm/a,约为GPS和地质等浅表观测的滑动速率值的2倍。研究结果对于理解断层活动性和地震危险性分析具有重要意义。     

Study of fault slip modes

We present the data of the laboratory experiments on studying the regularities of gradual transition from the stick-slip behavior to aseismic creeping on the interblock boundary. The experiments show that small variations in the material composition in the principal slip zones of the faults may cause a significant change in the fraction of seismic energy radiated during the dynamic unloading of the adjacent segment of the rock mass. The experiments simulate interblock sliding regimes with the values of the scaled kinetic energy differing by a few orders of magnitude and relatively small distinctions in the strength of the contacts and in the amplitude of the released shear stresses. The results of the experiments show that the slip mode and the fraction of the deformation energy that goes into the seismic radiation are determined by the ratio of two parameters—the stiffness of the fault and the stiffness of the enclosing rock mass. An important implication of the study for solving the engineering tasks is that for bringing a stressed segment of a fault or a crack into a slip mode with low-intensity radiation of seismic energy, the anthropogenic impact should be aimed at diminishing the stiffness of the fault zone rather than at releasing the excessive stresses.     

Spatiotemporal patterns of fault slip rates across the Central Sierra Nevada frontal fault zone

Patterns in fault slip rates through time and space are examined across the transition from the Sierra Nevada to the Eastern California Shear Zone–Walker Lane belt. At each of four sites along the eastern Sierra Nevada frontal fault zone between 38 and 39° N latitude, geomorphic markers, such as glacial moraines and outwash terraces, are displaced by a suite of range-front normal faults. Using geomorphic mapping, surveying, and ¹⁰Be surface exposure dating, mean fault slip rates are defined, and by utilizing markers of different ages (generally, ~ 20 ka and ~ 150 ka), rates through time and interactions among multiple faults are examined over 10⁴–10⁵ year timescales.At each site for which data are available for the last ~ 150 ky, mean slip rates across the Sierra Nevada frontal fault zone have probably not varied by more than a factor of two over time spans equal to half of the total time interval (~ 20 ky and ~ 150 ky timescales): 0.3 ± 0.1 mm year^? 1 (mode and 95% CI) at both Buckeye Creek in the Bridgeport basin and Sonora Junction; and 0.4 + 0.3/?0.1 mm year^? 1 along the West Fork of the Carson River at Woodfords. Data permit rates that are relatively constant over the time scales examined. In contrast, slip rates are highly variable in space over the last ~ 20 ky. Slip rates decrease by a factor of 3–5 northward over a distance of ~ 20 km between the northern Mono Basin (1.3 + 0.6/?0.3 mm year^? 1 at Lundy Canyon site) to the Bridgeport Basin (0.3 ± 0.1 mm year^? 1). The 3-fold decrease in the slip rate on the Sierra Nevada frontal fault zone northward from Mono Basin is indicative of a change in the character of faulting north of the Mina Deflection as extension is transferred eastward onto normal faults between the Sierra Nevada and Walker Lane belt.A compilation of regional deformation rates reveals that the spatial pattern of extension rates changes along strike of the Eastern California Shear Zone-Walker Lane belt. South of the Mina Deflection, extension is accommodated within a diffuse zone of normal and oblique faults, with extension rates increasing northward on the Fish Lake Valley fault. Where faults of the Eastern California Shear Zone terminate northward into the Mina Deflection, extension rates increase northward along the Sierra Nevada frontal fault zone to ~ 0.7 mm year^? 1 in northern Mono Basin. This spatial pattern suggests that extension is transferred from more easterly fault systems, e.g., Fish Lake Valley fault, and localized on the Sierra Nevada frontal fault zone as the Eastern California Shear Zone–Walker Lane belt faulting is transferred through the Mina Deflection.     

Large slip rate detected at the seismogenic zone of the 2008 M_W7.9 Wenchuan earthquake

Repeating microearthquakes were identified along the edge of the rupture area of the 2008 MW7.9 Wenchuan earthquake.Slip rates at depths derived from seismic moments and recurrence intervals are found to be systematically larger than those observed at surface.This large deep slip rate may explain the odds about the occurrence of this unanticipated event.Our observations here suggested that seismic hazard could be underestimated if surface measurements alone are employed.     

Microseismicity but no tremor accompanying slow slip in the Hikurangi subduction zone,New Zealand

Geodetically-detected episodes of slow slip appear in several subduction zones to be accompanied by bursts of low-frequency coherent noise known as seismic tremor, but whether a single physical process governs this association or even whether slow slip is invariably accompanied by tremor remains unresolved. Detailed analysis of broadband seismic data spanning a slow slip episode in the Hikurangi subduction zone, New Zealand, reveals that slow slip was accompanied by distinct reverse-faulting microearthquakes, rather than tremor. The timing, location, and faulting style of these earthquakes are consistent with stress triggering down-dip of the slow slip patch, either on the subduction interface or just below it. These results indicate that tremor is not ubiquitous during subduction zone slow slip, and that slow slip in subduction zone environments is capable of triggering high-frequency earthquakes near the base of the locked subduction thrust. In this and other locations (Hawaii, Boso Peninsula) where slow slip is accompanied by triggered microseismicity, the estimated upper extent of the slow slip is shallower (less than ~ 20 km) than in those locations from which tremor has been reported. This suggests that ambient temperature- or pressure-dependent factors govern the character of the seismic response to slow slip on subduction thrusts and other large faults, with rheological or lithological conditions at shallow depths triggering high-frequency microearthquakes and those at greater depths triggering seismic tremor.     

Determination of the direction and magnitude of recent tectonic stress in the Xianshuihe fault zone using fault slip data

ＤｅｔｅｒｍｉｎａｔｉｏｎｏｆｔｈｅｄｉｒｅｃｔｉｏｎａｎｄｍａｇｎｉｔｕｄｅｏｆｒｅｃｅｎｔｔｅｃｔｏｎｉｃｓｔｒｅｓｓｉｎｔｈｅＸｉａｎｓｈｕｉｈｅｆａｕｌｔｚｏｎｅｕｓｉｎｇｆａｕｌｔｓｌｉｐｄａｔａ谢富仁，李宏Ｆｕ－ＲｅｎＸＩＥａｎｄＨｏｎｇ...     

The slip rate and strong earthquake recurrence interval on the Qianning-Kangding segment of the Xianshuihe fault zone

IntroductionLocated in the western part of Sichuan Province, China, the Xianshuihe fault zone is a notable strong earthquake zone in the eastern Tibetan Plateau. At its northwestern end, the Xianshuihe fault zone overlaps the Ganzi-Yushu fault in a left-stepping pattern. The fault zone extends southeastwards through Luhuo, Daofu, Kangding, and Moxi and disappears near Shimian, with a total length of 400 km. The fault trends N40(-50°W in the north, and N20(-30°W to the south of Kangdi…     

利用重复微震估算红河断裂带北段深部滑动速率

利用云南数字地震台网1999-2015年的波形资料,通过识别在断层同一位置发生的重复微震来估算红河断裂带北段不同孕震深度的断层滑动速率.针对研究区台站分布稀疏的客观情况,采用基于子采样条件下S-P相对到时差约束地震相对位置,从而确保震源位置一致性的方法,在红河断裂带北段及其附近区域识别出23组重复微震,其复发间隔变化明显.重复微震的时空变化初步分析表明,研究区频繁的背景地震活动和相邻较近的重复微震之间的相互作用都会对重复微震的复发间隔产生影响.利用重复微震的震级和复发间隔,估算得到红河断裂带北段6.0~13.4 km孕震深处的滑动速率为2.3~10.0 mm·a^-1,显示不同孕震深度的滑动速率存在差异.

     

海原断裂带老虎山段晚第四纪滑动速率精确厘定与讨论

如何准确测定断裂滑动速率是近年来活动构造研究的前沿与热点.随着高精度地形数据获取手段与第四纪测年方法的不断进步,位错量和地貌面年龄的精度均得到大大提高.在进行滑动速率计算时还要考虑地质过程是否合理,蒙特卡洛方法为获取更加符合地质过程的滑动速率提供了重要工具.本文以滑动速率研究程度较低的海原断裂带老虎山段为例,基于LiDAR高精度地形数据,测得T1—T4阶地面年龄分别为1~3 ka,9~11 ka,15~17 ka,40~45 ka,陡坎前缘的位错分别为7~14 m,28~36 m,59~66 m,180~190 m.综合多地点的左旋走滑位错量及不同时代的地貌面年龄数据,并考虑滑动历史,利用蒙特卡洛模拟方法,将位错-时间两个参数的不确定性定量化,限定老虎山断裂45 ka以来平均滑动速率为4.3±0.16 mm·a^-1,17ka以来的平均滑动速率为4.0±0.15 mm·a^-1,与前人研究得到的狭义海原断裂滑动速率4.5±1.0 mm·a^-1基本一致.综合整个海原断裂带滑动速率,本文结果更支持低滑动速率变化趋势,即海原断裂带整体滑动速率趋于稳定,向东至六盘山断裂,滑动速率开始降低,推测海原断裂带的左旋走滑在尾端主要为马东山—六盘山隆起所吸收.结合老虎山断裂历史地震资料和深部锁闭浅部蠕滑的动力学特征,推测老虎山断裂具备与相邻断裂一起触发强震的能力.

     

甘孜—玉树断裂当江段晚第四纪滑动速率

甘孜—玉树断裂带是青藏高原中东部的一条大型左旋走滑断裂带,同时也是羌塘地体和巴颜喀拉地体的重要地质边界.当江断裂位于甘孜—玉树断裂带的西北段,沿线发育当江荣、当江和哲达等一系列串珠状第四纪断层谷地.通过遥感影像解译和数字高程地形模型(DEM)数据分析,结合野外构造地貌调查,以及断错地貌面的光释光年代测定,发现断裂沿线冲沟、河流阶地和洪积扇等断错地貌发育,反映了该断裂晚第四纪左旋走滑活动性强烈.该断裂最新活动时代为全新世晚期,距今约3.04 ka.当江断裂晚更新世以来的左旋滑动速率为7±3mm·a~(-1).研究结果为该区的地震危险性分析和高原东北部的运动学特征探讨提供了基础资料.     

综合利用多种地壳形变观测资料计算鲜水河断裂带现今滑动速率

本文首先沿走向将鲜水河断裂带划分为炉霍、道孚、乾宁、康定和磨西五个断裂段,利用沿断裂带布设的跨断层短基线、短水准场地测量资料计算了近场的断层活动参数,利用覆盖断裂带相对较大区域的重力、GPS观测资料计算了重力场动态变化、GPS速度场.基于重力场动态变化和GPS速度场采用蚁群算法和粒子群算法（具有全局优化的优势）分别反演计算了五个断裂段断层活动参数,将结果中的走滑分量作为五个断裂段的现今走滑速率.通过对以上三类现今走滑速率及五个断裂段的地质平均滑动速率进行融合与对比分析,将重力资料反演计算结果作为断裂带整体走滑速率,与跨断层短基线、短水准测量计算的断层滑动速率结果进行对比分析,初步判定了各跨断层短基线、短水准场地所跨断裂的性质,最终给出了五个断裂段的现今整体左旋走滑速率和部分分支断裂左旋走滑速率,结果为：①炉霍段为9.13 mm·a^-1,虾拉沱区域西支断裂为2.46 mm·a^-1,东支断裂为5.84 mm·a^-1.②道孚段为8.57 mm·a^-1,东南段沟普区域西支断裂为1.78 mm·a^-1,东支断裂为6.79 mm·a^-1.③乾宁段为7.67 mm·a^-1.④康定段为6.14 mm·a^-1.⑤磨西段为4.41 mm·a^-1.本文还定性讨论了断裂带两侧重力、GPS测点覆盖范围内活动地块的三维弹塑性变形和古地震、历史地震造成的永久位错.

     

台湾—吕宋会聚带的地壳运动特征及其动力学机制

南海东部的台湾-吕宋会聚带是南海四个边界中构造背景最为复杂、构造活动最为活跃.本文收集该区的GPS速度场资料,通过对速度场进行样条插值获得了该区连续的速度场、主应变率场、最大剪应变率场等结果.研究发现,该区的地壳运动受西北侧华南陆缘基底隆起和西南侧巴拉望岛阻挡、台湾北部24°N俯冲极转换、琉球海沟弧后扩张,以及菲律宾大断裂高速走滑等关键构造因素的影响,可分为(1)受华南地块运动影响区;(2)台湾岛上俯冲极转换部位至宜兰盆地之间受弧后扩张影响区;(3)俯冲极转换部位以南受菲律宾海板块俯冲和北港隆起阻挡影响区;(4)吕宋岛菲律宾大断裂以东受菲律宾海板块运动影响区;(5)大断裂以西受菲律宾海板块运动和巴拉望阻挡影响区;(6)巴拉望受巽他块体运动影响区6个地壳运动与变形区.     

Prevalent thickening and local thinning of the mantle transition zone beneath the Baikal rift zone and its dynamic implications

The Baikal rift is the most seismically active continental rift in the world and is significant for studying the dynamics of continental rifts, although its precise dynamic mechanisms remain controversial. We calculated receiver functions (1748) from Global Seismographic Network seismic stations TLY and ULN and stacked receiver functions in different bins. Here we present discontinuities at depths of 410km and 660km and thickness of the mantle transition zone (MTZ) beneath the study area. The MTZ structure shows an obvious thickening (292km) in the Baikal rift zone except for an area of limited thinning (230km), whereas it is basically normal (250km) beneath the Mongolian area, to the southeast of the Baikal rift. Combining these results with previous findings, we propose that the large-scale thickening beneath the Baikal rift zone is likely to be caused by the Mesozoic collision between the Siberian Platform and the Mongolia-North China Block or magmatic intrusion into the lower crust, which would result in crust and lithosphere thickening. Thus, the lower crust becomes eclogitized and consequently detached into the deep mantle because of negative buoyancy. The detachment not only induces asthenosphere upwelling but also accelerates mantle convection of water detached from the subducted slab, which would increase mantle melting, while both processes promote the development of the rift. Our preliminary results indicate that the detachment and the consequent hot upwelling have an important influence on the development of the Baikal rift, and a small-scale mantle upwelling indicated by the located thinning may have destroyed the lithosphere and promoted this development.