Stick-slip and strain waves in the physics of earthquake rupture: experiments and models

The study of stick-slip is directly connected with the formulation of the consistent concept of seismicity migration and explanation of the nature of strain waves observed at the contacts of blocks and plates. We propose a phenomenological model to describe the initiation of stick-slip at the rough contact of blocks of rocks. The model contains the leading factors of stick-slip process (friction, roughness of contact surfaces, and asperity). The model reproduces the universal profile of slip velocity and displacement, velocity of dynamic rupture observed in experiments. The stick-slip motion at a rough surface of a fault is shown to be a nonlinear process and can be described by the generalized sine-Gordon equation. The results of calculations and a comparison with experimental data testify to the existence of waves of a new type — solitary waves of sliding. The analogies are found between the strain waves generated due to stick-slip at the contact of the blocks of rocks, strain waves in the crustal faults, and those within the lithosphere.     

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

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

汶川8.0级地震前GPS与跨断层资料反映的运动与变形演化特征

为了研究与总结2008年5月12日汶川8.0级地震前GPS与跨断层资料反映的龙门山断裂带及其周边地区的运动、构造变形、应变积累演化过程,以及汶川地震临震阶段可能的物理机制,本文综合1999~2007期GPS速度场、1999~2008年大尺度GPS基线时间序列、1985~2008年跨断层短水准等资料进行了相关分析与讨论。结果表明:(1)GPS速度剖面结果显示,宽达500km的川西高原在震前有明显的连续变形,而四川盆地一侧和跨龙门山断裂带基本没有变形趋势,表明震前川西高原在持续不断地为已经处于闭锁状态的龙门山断裂带提供能量积累。(2)GPS应变率结果显示,震前龙门山断裂带中北段的NW侧EW向挤压变形明显,变形幅度从远离断裂带较大到靠近断裂带逐渐减小,而断裂带变形微弱;龙门山断裂带西南段周边形成了显著的EW向挤压应变集中区,应变积累速率明显大于中北段。(3)断层闭锁程度反演结果显示,除了汶川地震的震源位置闭锁相对较弱,且西南段有大概20km宽度断层在12~22.5km深度为蠕滑状态以外,震前整条龙门山断裂基本处于强闭锁状态。(4)大尺度GPS基线结果显示,跨南北地震带区域的NE向基线从2005年开始普遍出现压缩转折,反映NE向地壳缩短的相对运动增强。(5)跨断层短水准场地结果显示,震前年均垂直变化速率和形变累积率很低,表明断层近场垂向活动很弱、闭锁较强。通过以上分析认为,在相对小尺度的地壳变形中,震前龙门山断裂带深浅部均处于强闭锁状态,断裂带水平与垂直变形都很微弱,这可能经历了一个缓慢的过程,而且越是临近地震的发生,微弱变形的范围可能越大;在相对大尺度的地壳变形中,震前龙门山断裂带西侧的巴颜喀拉块体东部地区经历了地壳缓慢且持续的缩短挤压变形,为龙门山断裂带应变积累持续提供了动力支持。     

饱和土体-地下综合管廊结构地震响应分析

将土体视为固-液两相介质,基于饱和土体有效应力原理,建立饱和土体-地下综合管廊结构体系相互作用动力模型:在地应力平衡的静力状态下采用Duncan-Chang非线性弹性本构模型,在地震波作用的动力状态下采用Davidenkov非线性黏弹性本构模型;考虑饱和土体黏弹性动力人工边界条件,将地震动作用转化为作用在人工边界节点上的动力荷载。模型考察不同地震波时程、地震波加速度峰值、入射角度、孔隙率以及地应力场的影响,得出如下结论:(1)地震波的卓越周期与场地卓越周期相近时引起结构上的变形最大;随着地震波加速度峰值的增大结构变形增大;随着地震波入射角度的增加结构变形增大,地震波斜入射情况下产生的行波效应使得结构变形最大。(2)土体材料的孔隙水压力是影响地震中结构变形的主要因素之一。(3)将土体材料考虑为单相介质时结构上的变形要比考虑为固-液两相介质时大得多,直接将饱和土体场地中得到的地震波等效荷载施加到单相土介质-结构动力相互作用模型上,能够得到与完全基于有效应力法一致的结果。     

有缝重力坝动力非线性数值计算模型及其应用

对于有缝重力坝的分析，基于连续介质力学的界面单元-有限元方法难于比较精确地模拟缝间的接触应力，从而无法合理地估算坝体应力与变形。作者对于多体系统和分区连续介质所发展的非连续变形计算力学模型能够根据接触界面的本构特性及其力学和运动学约束条件精确地再现受力过程中界面相互作用力的传递与非连续变形状态，本文将其应用于有缝重力坝的动力分析。实例数值分析表明该模型的计算结果从定性上讲是合理的，并且为判断坝体缝隙的工作状态与界面应力提供了有力的依据。     

Deformation Evolution Characteristics Revealed by GPS and Cross-fault Leveling Data before the MS8.0 Wenchuan Earthquake

Based on GPS velocity during 1999-2007, GPS baseline time series on large scale during 1999-2008 and cross-fault leveling data during 1985-2008, the paper makes some analysis and discussion to study and summarize the movement, tectonic deformation and strain accumulation evolution characteristics of the Longmenshan fault and the surrounding area before the M_S8.0 Wenchuan earthquake, as well as the possible physical mechanism late in the seismic cycle of the Wenchuan earthquake. Multiple results indicate that:GPS velocity profiles show that obvious continuous deformation across the eastern Qinghai-Tibetan Plateau before the earthquake was distributed across a zone at least 500km wide, while there was little deformation in Sichuan Basin and Longmenshan fault zone, which means that the eastern Qinghai-Tibetan Plateau provides energy accumulation for locked Longmenshan fault zone continuously. GPS strain rates show that the east-west compression deformation was larger in the northwest of the mid-northern segment of the Longmenshan fault zone, and deformation amplitude decreased gradually from far field to near fault zone, and there was little deformation in fault zone. The east-west compression deformation was significant surrounding the southwestern segment of the Longmenshan fault zone, and strain accumulation rate was larger than that of mid-northern segment. Fault locking indicates nearly whole Longmenshan fault was locked before the earthquake except the source of the earthquake which was weakly locked, and a 20km width patch in southwestern segment between 12km to 22.5km depth was in creeping state. GPS baseline time series in northeast direction on large scale became compressive generally from 2005 in the North-South Seismic Belt, which reflects that relative compression deformation enhances. The cross-fault leveling data show that annual vertical change rate and deformation trend accumulation rate in the Longmenshan fault zone were little, which indicates that vertical activity near the fault was very weak and the fault was tightly locked. According to analyses of GPS and cross-fault leveling data before the Wenchuan earthquake, we consider that the Longmenshan fault is tightly locked from the surface to the deep, and the horizontal and vertical deformation are weak surrounding the fault in relatively small-scale crustal deformation. The process of weak deformation may be slow, and weak deformation area may be larger when large earthquake is coming. Continuous and slow compression deformation across eastern Qinghai-Tibetan Plateau before the earthquake provides dynamic support for strain accumulation in the Longmenshan fault zone in relative large-scale crustal deformation.     

GPS应变率场计算方法研究进展

基于大地测量资料获取地壳运动与应变积累定量结果一直是国内外重视的地震中长期预测的技术途径。 针对地震变形过程的准确描述问题, 国内外学者基于GPS资料, 发展了多种应变率场解算方法。 本文首先简要介绍了GPS应变率计算的基本原理, 然后系统梳理了国内外多种计算方法的优势和不足, 结果表明: 应变率计算的数学方法只考虑几何关系, 其中整体方法主要适合数据密度和分布较好条件下获取区域地壳变形分布与趋势, 局部方法主要适用于数据较为稀疏情况下描述构造块体的变形特征; 应变率计算的物理方法既考虑几何关系又考虑物理关系, 其中, 位错方法根据主要适合于研究区域存在主控断层的情况(研究区域的变形主要由少数断层控制); 数值模拟方法(如有限元法)主要适用于区域地质、 地球物理的资料比较完备的情况。     

Focal process of the great Chilean earthquake May 22, 1960

Long-period strain seismogram recorded at Pasadena is used to determine the focal process of the 1960 Chilean earthquake. Synthetic seismograms computed for various fault models are matched with the observed strain seismogram to determine the fault parameters. A low-angle (～ 10°) thrust model with rupture length of 800 km and rupture velocity of 3.5 km/sec is consistent with the observed Rayleigh/Love wave ratio and the radiation asymmetry. A seismic moment of 2.7 · 10³⁰ dyn · cm is obtained for the main shock. This value, together with the estimated fault area of 1.6 · 10⁵ km², gives an average dislocation of 24 m. The strain seismogram clearly shows unusually long-period (300–600 sec) wave arriving at the P time of a large foreshock which occurred about 15 minutes before the main shock, suggesting a large slow deformation in the epicentral area prior to the major failure. A simple dislocation model shows that a dislocation of 30 m, having a time constant of 300–600 sec, over a fault plane of 800 × 200 km² is required to explain this precursory displacement. The entire focal process may be envisaged in terms of a large-scale deformation which started rather gradually and eventually triggered the foreshocks and the “main” shock. This mechanism may explain the large premonitory deformations documented, but not recorded instrumentally, for several Japanese earthquakes. The moments of the main shock and the precursor add to 6 · 10³⁰ dyn · cm which is large enough to affect the earth's polar motion.     

Seismomagnetic response of a fault zone

Based on the results of instrumental observations of geomagnetic variations caused by the propagation of seismic waves through a fault zone, the dependences between the amplitudes of the induced seismomagnetic effect and seismic signal as a function of distance r to the midline of the fault are obtained. For the first time, it is shown that the amplitude of the seismomagnetic effect is maximal in the fault damage zone. The phenomenological model describing the generation of magnetic signals by seismic waves propagating through the crushed rock in the tectonic fault zone is suggested. It is assumed that geomagnetic variations are generated by the changes in the electrical conductivity of the fragmented rocks as a result of the deformation of the rock pieces contacts. The amplitudes of the geomagnetic variations calculated from the model agree with the instrumental observations.     

基于不同结构损伤模型的RC框架地震易损性研究

针对目前RC框架结构地震易损性分析中整体损伤模型研究的薄弱性以及广泛采用的层间位移角方法不能准确反应结构在地震作用下损伤机理的现状,本文基于现有损伤模型的对比分析,提出了一种较准确反映地震破坏机理同时便于应用的最大变形和滞回耗能非线性组合的双参数损伤模型。以8层RC框架结构为例,进行50条地震波作用下的结构增量动力分析,分别绘制了变形和能量2种单参数模型以及牛荻涛模型和本文模型两种双参数模型的结构损伤曲线与易损性曲线,并进行了模型的对比分析和检验评估。分析结果表明:仅以层间位移角作为结构整体损伤指标会高估结构的抗倒塌性能,仅以能量作为结构整体损伤指标会低估结构损伤的超越概率。本文模型能较好地平衡最大变形和累积损伤对结构损伤的影响程度。     

我国西部近期应变强度比动态演化特征

探索由断层形变资料提取应变强度比并获得动态图像的方法，结合2000年景泰（Ms5.9）等震例，给出由我国西部甘、宁、青及川滇地区最新跨断层形变资料所得的应变强度比时空分布动态图像。结果表明：其具有压性为主的中期背景及中短期前兆（一年尺度的绝对值2.0以上高值区）等异常的动态演化特征可作为上述地区中短期强震危险区的判定依据。     

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

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年,研究断裂段将完全"闭锁",即进入积累下一次大地震应力应变的震间闭锁阶段.     

A fault model of the 1946 Nankaido earthquake derived from tsunami data

A fault model of the 1946 Nankaido earthquake (M = 8.2) is determined by the use of tsunami records of Uwajima, Shimotsu and Hososhima which were located within or near the area of major coseismic crustal deformation. Synthetic tsunamis computed for various fault models are matched with the observed tsunamis to determine the fault parameters. A low-angle thrust model slightly revised from a previous model by Ando is consistent with the observed tsunamis. The duration of faulting is constrained as less than 10 min based upon the tsunami. The fault is divided into an eastern and a western segment corresponding to areas associated with and without aftershocks, respectively. The fault area and dislocation for the western segment are 150 × 70 km² and 6 m, and those for the eastern segment are 150 × 70 km² and 3 m, respectively. The total seismic moment is 4.7 × 10²⁸ dyn·cm, significantly smaller than that obtained from a geodetic model by Fitch and Scholz, but still larger than that of the seismic model by Kanamori. The discrepancy in seismic moment between the seismic and the present models (RAN2) could be interpreted in terms of a slow dislocation on the fault, but this interpretation does not match the seismic intensity distribution and damage pattern, and the slow-slip model for the Nankaido earthquake is rejected. The discrepancy between the two seismic moments is considered insignificant within error involved in data and modeling assumptions. If the revised geodetic model (RAN2) is modified, the seismic moment required to explain the observed tsunamis would be reduced further by ～30%. If we consider the uncertainties involved in the fault model of Kanamori and the fault-finiteness effect affecting the amplitude of seismic waves, the seismic moment required to interpret the seismic-wave data could be increased, possibly being more than twice that of Kanamori. Thus, the two seismic moments from the different data sets could be close to each other within allowable tolerance. This implies that the rise time of the Nankaido earthquake was short enough to generate short-period seismic waves from both the western and the eastern fault segments.     

GNSS在中国大陆的地震预测应用研究进展与展望

GNSS技术的快速发展为地震预测研究提供了前所未有的大尺度、高精度的观测结果,为强震变形模型的发展提供了可靠的观测约束。本文针对GNSS技术在中国大陆地震预测中的应用,系统梳理了断裂带滑动特征描述、变形场动态演化解析、应变集中过程识别、潜在震源危险程度判断等方面的研究进展。通过典型震例总结了GNSS资料在长、中、短临不同的地震预测阶段的应用。针对地震中长期预测,基于构造动力过程给出了强震危险性时空逼近的科学思路,即“板块边界动力作用—大-中尺度动态形变场—应力应变增强/集中区—孕震危险段中短期危险性的时空逼近”的过程。在此基础上,针对GNSS监测能力提升、地震孕育过程相关的多尺度地壳形变动态信息获取、GNSS多参量动力学模型构建及产出等问题进行了讨论和发展展望。总体而言,GNSS技术的应用显著增强了我国地震预测的地壳形变观测基础支撑,丰富了对大陆地震孕育发生物理过程的科学认识,推动了大陆地震预测科学思路和预测方法的发展,并促进了地震预测由经验预测向物理预测的拓展。     

The genesis of seismic activity on faults in Central Asia in real time and its variations

We used the Digital Faults geoinformation system that we developed to propose an algorithm for quantitative estimation of seismic activity on faults. The resulting technique was used to study the spatiotemporal patterns in the present-day activity of faults in Central Asia. Fault activity was found to vary at frequencies of a few years and cannot be explained by changes in the regional stress fields. We studied the tendency of seismic events to be localized to areas of dynamic influence due to faults. The active faults were grouped by the criteria of seismicity organization in the influence areas of these faults. It was shown that fault activity and its comparatively high frequency on real time scales are caused by strain waves, which may be generated by interplate and interblock movements in the brittle lithosphere. Judging by the speed of strain waves, the active faults are classified into groups that differ in their geological and geophysical parameters. They can be used to estimate the directions of strain wave fronts and to identify areas of dominant fault activation over intervals of real (geologically speaking) time. We give a map showing active faults in Central Asia, plots of a quantitative index of their seismic activity, and the directivity vectors of strain waves that excite fault activity. The methods we developed for classifying active faults by the quantitative index of seismic activity and for determining the vectors of strain waves that excite fault activity are all tools that significantly expand our possibilities when developing tectonophysical models of the seismic process in earthquake-generating zones of the lithosphere and open new methods for attacking problems in intermediate-term earthquake prediction.     

GNSS在断层形变研究中的应用

大地震多发生在主要活动断裂带及其附近, 因此断层形变的监测与研究非常重要。 多年来采用跨断层短水准、 短基线、 跨断裂带综合观测剖面等手段监测断层活动, 进而研究断层活动与变形规律。 GNSS多用于大面积地壳水平运动场和主要地块活动与变形的研究, 随着GNSS数据的积累和应用研究的深入, 已有学者用GNSS从不同角度开展断层形变的监测与研究。 本文介绍和论述了GNSS在断层形变研究中已有的部分应用和新的探索, 主要介绍以下几个方面: ① GNSS用于断层走滑活动的监测与研究; ② GNSS用于断层垂直活动的监测与研究; ③ 通过研究断层两边地块变形规律并建立应变模型, 可分别依据各自模型推测断层的活动与形变, 将推测结果与跨断层观测到的实际情况进行研究与比较。 研究分析认为, 由块体推测的变形与实测变形差异越大, 说明块体相互作用越强烈。 这可能对强地震地点预测、 地震危险性判定及地震趋势预测均有重要意义。     

亚失稳准动态及同震过程变形场时空演化特征 ——实验与分析

断层亚失稳模型指出,在临震亚失稳阶段中各种物理量存在规律性的时空演化特征,控制这些物理参数变化的根本原因是震源的力学过程.为深入观测和分析该过程,文中介绍了一套自主研发的64通道、16位分辨率、4MHz采样频率、可并行连续采集的超动态变形场观测系统(Ultra-HiDAM),首次实现了在4MHz频率下对应变信号和声发射...     

Characteristics of Recent Horizontal Crustal Movement and Tectonic Deformation in the Northwest China Region

INTRODUCTIONThe vector diagramof the crustal horizontal movement of the north margin of Qinghai-Xizang(Tibet) block obtained from GPS measurements of 2001 ~2003 shows that the western part of theregion presents a southwestward movement oppositetothetrend beforethe Nov.14 ,2001 west KunlunMountains PassMS8 .1 earthquake ,which mayreflectthe post-earthquake relaxationand adjustment ofthe region (Zhang Xi ,2004) .It was the variation of regional dynamics that led to instability of thepo…