Study of dislocation model and evolution characteristics of vertical deformation field of gonghe M S=6.9 earthquake

The preliminary research results of vertical deformation dislocation model of GongheM _S =6.9 earthquake show that, the causative structure is a hidden fault with strike N60°W, dipping S47°W, which lies near the current subsidence center of Gonghe basin. The rupture length and width are 30km and 14km, the upper and lower bound depth of the fault in width direction are 3km and 17km respectively. The maximum coseismic and preseismic vertical deformation of GongheM _S =6.9 earthquake are 247mm and about 100mm. The reasons why there existed rapid postseismic uplift are also given a tentative discussion. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 289–295, 1993.     

1990年青海共和7.0级地震震后垂直形变研究

根据1990年青海共和地震震后地表垂直形变,通过模型拟合得到了支配共和地区震后形变场时空演化的形变源及其力学机制.分析穿过断层的震前1期和震后6期水准数据,结果表明震后垂直形变具有以下特征:①震后震区上盘继续发生继承性的大幅度上升,其中震后头一年上升速率最大;②震后上升区范围显著,范围随时间变化不大,但较同震形变上升区范围增大;③震后相邻测站高差观测值的时间序列明显具有对数衰减特征或指数衰减特征,衰减特征时间分别为0.165年和1.344年.本文还发展了一个利用水准数据与连续介质位错模型研究震后形变机制的新方法.该方法用相邻水准点之间的原始高差观测值而非它们相对参照点的积分值来约束连续介质位错模型,可以有效减少误差累积带来的偏差并充分利用观测数据.利用这一方法的初步分析结果表明,断层震后滑移和介质黏弹性松弛共同导致了共和地震震后形变.前者表现为发生在断层面及其延伸部分的滑移,特别是位于主破裂上方沉积层内的滑移;后者则表现为下地壳与上地幔内的黏弹性松弛,黏滞系数为1020Pa.s量级.     

1990年共和7.0级地震的发震构造讨论

根据地质、地震及形变资料研究了１９９０年４月青海省共和ＭＳ７．０地震的发震构造．认为ＮＷＮＷＷ向的哇玉香卡拉干逆冲盲断裂是该次地震的发震构造．其地表形变以褶皱隆起为主．因此,共和地震是一次典型的“盲断裂褶皱地震”．     

Application of rock creep experiment in calculating the viscoelastic parameters of earth medium

The viscosity of earth medium can explain many geophysical phenomena. The uplifting velocity of Ti- betan Plateau may be caused by extrusion of the low viscosity channel in middle and lower crust[1]. The post-deformation after earthquake also relates to the rheology of earth medium[2]. The aftershock sequence is considered a process of strain release after stress dropped[3]. There are two kinds of models in the study of rheological character of crustal medium, elastoplas- only when the forced…     

Mechanisms of Transient Postseismic Deformation Following the 2001 Mw 7.8 Kunlun (China) Earthquake

Using global positioning system (GPS) technology, significant postseismic surface displacements were observed within the first 4 months after the 2001 Mw 7.8 Kunlun earthquake which occurred in China. In this study, we investigated the mechanisms that may have possibly contributed to the postseismic deformations that have been observed. Based on the modeling results, we find that an afterslip model can interpret postseismic displacements in the near field even when the fault plane is extended to the bottom of the crust (~70 km). Models based on the viscoelastic relaxation theory showed a large discrepancy in the spatial pattern of the deformation compared with what has been observed. Thus, we infer that both mechanisms cannot interpret the observed postseismic deformation independently. A combination of afterslip and viscoelastic relaxation can further improve the data fit, especially at sites far from the fault. With maximum afterslip of ~0.4 m occurring at a depth of 10 km in the central section, the combined model shows that the estimated afterslip occurred mostly on and below the coseismic rupture plane, as well as on its eastern extension. The estimated moment released by the afterslip in the first 4 months is almost 40% of that released by the coseismic slip. The best-fitting viscoelastic relaxation model shows a “weak” upper mantle with a viscosity of ~1.0 × 10¹⁸ Pa s. The combined model also suggests the existence of a lower crust with viscosity larger than 1.0 × 10¹⁸ Pa s, although it cannot be constrained accurately.     

Study of dislocation model and evolution characteristics of vertical deformation field of gongheM S=6.9 earthquake

The preliminary research results of vertical deformation dislocation model of GongheM _S =6.9 earthquake show that, the causative structure is a hidden fault with strike N60°W, dipping S47°W, which lies near the current subsidence center of Gonghe basin. The rupture length and width are 30km and 14km, the upper and lower bound depth of the fault in width direction are 3km and 17km respectively. The maximum coseismic and preseismic vertical deformation of GongheM _S =6.9 earthquake are 247mm and about 100mm. The reasons why there existed rapid postseismic uplift are also given a tentative discussion.     

昆仑山Ms 8.1级地震震后变形场数值模拟与成因机理探讨

2001年11月14日,在青海和新疆交界处发生了昆仑山M_s8.1级强烈地震,GPS后观测显示,此次地震震后形变不仅在断裂南北两侧存在很大的差异,而且在短时间调整后断裂南北两侧表现为同向运动.本文以观测的地震形变为约束,通过有限元数值模拟分析昆仑山地震震后形变的物理机制.建立有关的有限元虚功方程,通过有限元数值方法模拟震后形变,从理论上分析介质的非均匀性、黏滞性松弛、流体调整对震后形变的影响.采用网格搜索确定昆仑断裂南北两侧下地壳的黏滞系数分别为5.0×10¹⁷Pa·s, 9.0×10¹⁸Pa·s左右,正是这十余倍的差异引起了断裂两侧震后形变的非对称性和同向运动,这一差异既是长期地质作用的结果,又是现代地球动力学环境的决定因素之一.通过数值模拟定性讨论了断裂北侧地表形变在震后短期内的调整,对于靠近断裂附近的测点可能是黏弹性松弛和孔隙流体调整共同作用的结果,所以在分析短期震后形变时综合考虑黏弹松弛和孔隙流体调整是很有必要的.     

1999年台湾集集地震震后地表变形的力学机制

1999年台湾集集地震震后450天的GPS观测资料显示了几十到几百毫米的地表位移.下地壳的震后黏性松弛和断层无震蠕变产生的震后滑动是用来解释地表震后变形的两个主要机制.本文利用接触问题的黏弹性有限元(LDDA)方法,以GPS观测数据作为约束,分别考察了黏性松弛和震后滑动机制对地表震后变形的影响.计算结果表明,黏性松弛机制产生的地表位移与观测数据吻合较好,通过试错法由震后GPS观测约束得到的下地壳黏度为10¹⁷Pa·s,而上地幔黏度对计算结果影响不大.考察震后滑动机制对地表变形的影响时,在LDDA方法中结合了速率状态摩擦定律,结果显示震后滑动机制不能很好地解释震后450天的观测数据,它产生的地表变形只在震后50天内与观测大致吻合,之后位移值基本不随时间变化.这些结果有助于增进对集集地震震后变形机制的认识.     

Viscoelastic relaxation of the upper mantle and afterslip following the 2014 MW8.1 Iquique earthquake

An improved understanding of postseismic crustal deformation following large subduction earthquakes may help to better understand the rheological properties of upper mantle and the slip behavior of subduction interface. Here we construct a three-dimensional viscoelastic finite element model to study the postseismic deformation of the 2014 M_W8.1 Iquique, Chile earthquake. Elastic units in the model include the subducting slab, continental and oceanic lithospheres. Rheological units include the mantle wedge, the oceanic asthenosphere and upper mantle. We use a 2 ​km thick weak shear zone attached to the subduction fault to simulate the time-dependent stress-driven afterslip. The viscoelastic relaxation in the rheological units is represented by the Burgers rheology. We carry out grid-searches on the shear zone viscosity, thickness and viscosity of the asthenosphere, and they are determined to be 10¹⁷ ​Pa ​s, 110 ​km and 2 ​× ​10¹⁸ ​Pa ​s, respectively. The stress-driven afterlsip within the first two years is up to ~47 ​cm and becomes negligible after two years (no more than 5 ​cm/yr). Our results suggest that a thin, low-viscosity oceanic asthenosphere together with a weak shear zone attached to the fault are required to better reproduce the observed postseismic deformation.     

Dynamic mechanisms of the post-seismic deformation following large events: Case study of the 1999 Chi-Chi earthquake in Taiwan of China

The mechanism of postseismic deformation related to strong earthquakes is important in geodynamics, and presumably afterslip or viscoelastic relaxation is responsible for the postsesimic deformation. The 1999 Chi-Chi, Taiwan of China, earthquake occurred in the region where GPS observation station is most densely deployed in the world. The unprecedented GPS data provides a unique opportunity to study the physical processes of postseismic deformation. Here we assume that the interactions of viscoelastic relaxation, afterslip, fault zone collapse, poroelastic rebound, flow of underground fluids, and all these combined contribute to the surface displacements following the main shock. In order to know the essence of the postseismic deformation after the strong event, fault zone collapse, poroelastic rebound, flow of underground fluids, and so on, are represented equivalently by the variations of the focal medium properties. Therefore, the viscoelastic relaxation, afterslip, and the variations of the equivalent focal medium properties are inverted by applying the GPS temporal series measurement data with viscoelastic finite element method. Both the afterslip rate distribution along the fault and the afterslip evolution with time are obtained by means of inversion. Also, the preliminary result suggests that viscosities of the lower crust and the upper mantle in Taiwan region is 2.7×10¹⁸ and 4.2×10²⁰ Pa·s, respectively. Moreover, the inversion results indicate that the afterslip contributing to postseismic deformation of 44.6% in 450 days after the Chi-Chi earthquake, with 34.7% caused by the viscous relaxation and 20.7% by other factors such as fault zone collapse, poroelastic rebound, and the flow of liquids.     

Surface rupture zone of the 1303 Hongtong M=8 earthquake, Shanxi Province

Based on the latest displacement of Huoshan piedmont fault, Mianshan west-side fault and Taigu fault obtained from the beginning of 1990‘s up to the present, the characteristics of distribution and displacement of surface rupture zone of the 1303 Hongtong M = 8 earthquake, Shanxi Province are synthesized and discussed in the paper. If Taigu fault, Mianshan west-side fault and Huoshan piedmont fault were contemporarily active during the 1303 Hongtong M = 8 earthquake, the surface rupture zone would be 160 km long and could be divided into 3 segments, that is, the 50-km-long Huoshan piedmont fault segment, 35-km-long Mianshan west-side fault segment and 70-km-long Taigu fault segment, respectively. Among them, there exist 4 km and 8 km step regions. The surface rupture zone exhibits right-lateral features. The displacements of northern and central segments are respectively 6～7 m and the southern segment has the maximum displacement of 10 m. The single basin-boundary fault of Shanxi fault-depression system usually corresponds to M ≈ 7 earthquake, while this great earthquake (M = 8) broke through the obstacle between two basins. It shows that the surface rupture scale of great earthquake is changeable.     

1976年唐山地震震时和震后变形的模拟

本文采用三维粘弹性有限元方法拟合唐山地区１９７６-１９８５年观测到的地震震时和震后的水平与垂直地形变，反演华北板块下方深部物质的流变学性质．模型采用多层弹性覆盖层与线性粘弹性层的有限块体，发展断层面上存在着位错运动，并用正交设计法拟合观测数据．模拟计算表明，华北板块下方软流层粘度为７．１×１０^１８Ｐａ·ｓ；上地幔粘度为２．１×１０^１９Ｐａ·ｓ．     

Combined GPS and InSAR Models of Postseismic Deformation from the Northridge Earthquake

--Models of combined Global Positioning System (GPS) and Interferometric Synthetic Aperture Radar (InSAR) data collected in the region of the Northridge earthquake indicate that significant afterslip on the main fault occurred following the earthquake. Additional shallow deformation occurred to the west of the main rupture plane. Both data sets are consistent with logarithmic time-dependent behavior following the earthquake indicative of afterslip rather than postseismic relaxation. Aftershocks account for only about 10% of the postseismic motion. The two data sets are complimentary in determining the postseismic processes. Fault afterslip and shallow deformation dominate the deformation field in the two years following the earthquake. Lower crustal deformation may play an important role later in the earthquake cycle.     

Inversion of long-period body-wave data for the source process of the Gonghe, Qinghai, China earthquake

Long period body waves data recorded by the China Digital Seismograph Network (CDSN) are inverted for the seismic moment tensors of the April 26, 1990, Gonghe, QinghaiM _S=6.9 earthquake and itsM _S=5.0 after-shock occurred on May 7, 1990. In the inversion, the generalized reflection-transmission coefficient matrix method is used to generate Green’s function. From the inversion it is obtained that the rupture process of theM _S=5.0 aftershock is relatively simple, and that of the main shock is rather complex. There are at least two events during main shock rupture process with an interval about 35 seconds. The focal mechanisms of two events are roughly the same as that of the aftershock, all of them were mainly reverse dip-slipping faulting with minor left-lateral strike-slip motion. These results indicate that the Gonghe earthquake was the result of the farther extension of one NWW-SEE striking buried fault on the southern margin of Gonghe basin from shallower depth to deeper depth and from NW to SE under the action of a nearly horizontal NE direction compressive stress. Contribution No. 95A0111, Institute of Geophysics, SSB, China.     

GPS测定的四川九寨沟MW6.6地震同震形变

2017年8月8日四川阿坝州九寨沟发生M_W6.6地震,震源机制解显示该地震为左旋走滑型地震。对震中周围的GPS连续站观测资料进行处理,获得高频GPS动态形变和静态同震水平位移。震中100km范围内四川松潘和甘肃武都站观测到1 Hz动态形变。距离震中约69km的松潘站观测的同震水平位移为7.4mm。根据少量的GPS静态同震位移反演的同震破裂模型显示本次地震的最大滑动量为376mm,地震矩为7.25×1018 N·m,等效矩震级为M_W6.6。正演计算的同震三维形变场显示本次地震的最大水平位移可达4~5cm,垂直位移呈四象限分布,最大可达1.5cm,区域内10个流动GPS站可观测到同震形变。     

The features of multiple fractures associated with the great Haiyuan 8.5 magnitude earthquake of 1920

It is deduced on the basis of field investigation that the total length of the stratigraphic fault associated with the great Haiyuan 8.5 magnitude earthquake of 1920 was 225 km. This fault was formed by 6 secondary faults with different geometric parameters, which align regularly inen echelon arrangement. Each secondary fault can be divided into three segments with different characteristics of deformation where the middle segment was mainly of the horizontal strike—slip fault while another two segments the vertical deformation as shown by the features of reverse or normal faults. It is also shown by the data of vertical and horizontal displacements that the horizontal displacement approached a maximum at the middle segment for each secondary fault and gradually decreased toward and finally disappeared at both ends of each segment while in contrast the vertical displacement was minimum at the middle and became large at both ends of the segment. The feature of the multiple peaks appeared in the deformation as shown by the earthquake displacements along the whole fault. This feature indicates that the 6 secondary faults associated with the great Haiyuan earthquake were the horizontal interrupted planes (i.e., dislocation surface) which were independent on each other, and hence each dislocation surface may represent an independent secondary fracture event of the earthquake. We thus think that the 6 relatively independent secondary events which occurred successfully might result in the great 8.5 magnitude Haiyuan earthquake. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 21–31, 1991. This project is sponsored by the Chinese Joint Seismological Science Foundation.     

2010年青海玉树地震同震-震后形变场特征及演化过程

基于星载合成孔径雷达差分干涉测量技术(DInSAR)和4期ENVISAT/ASAR雷达数据,获得了不同时间基线的三个同震干涉形变场和两个震后干涉形变场,并对这五个在时段上互有重叠的形变场进行了综合分析.结果表明,玉树地震同震形变场为围绕发震断层NW展布的椭圆形干涉条纹,覆盖范围约89 km×59 km.断层运动性质为左旋走滑.两盘最大视线向相对形变量至少达45 cm,最大形变出现在结古镇附近.时间基线不同的同震形变场总体上基本一致,但两盘最大相对形变量和局部形变存在差异.震后时间较长的干涉对反映的最大形变量反而减小;在震后时间较短的干涉对上于结古镇西南侧观测到的局部形变,在震后时间较长的干涉对上却没有出现.分析认为在形变量最大的结古镇附近可能出现了震后快速弹性回弹,导致随震后时间延续,形变量反而减小的现象.玉树地震震后形变主要出现在断层附近、震后不久的时段内,形变量在8 cm以下,具有与同震方向一致和相反的两种震后形变方式.在结古镇西南观测到一个与同震形变相反的局部沉降,应为震后弹性回弹.在微观震中处的断层附近观测到与同震方向一致的震后形变,可能是震后余滑.通过对地震前后不同时间基线的多个干涉对的联合对比分析,可以在一定程度上区分同震形变与震后形变,更好地研究地震引起的变形过程,特别是地震断层附近短期震后形变场的演化过程,为进一步研究断层带的岩性特征、物理力学及运动特性提供约束.     

Source process of the 1990 Gonghe,China, earthquake and tectonic stress field in the northeastern Qinghai-Xizang (Tibetan) plateau

TheM _s =6.9 Gonghe, China, earthquake of April 26, 1990 is the largest earthquake to have been documented historically as well as recorded instrumentally in the northeastern Qinghai-Xizang (Tibetan) plateau. The source process of this earthquake and the tectonic stress field in the northeastern Qinghai-Xizang plateau are investigated using geodetic and seismic data. The leveling data are used to invert the focal mechanism, the shape of the slipped region and the slip distribution on the fault plane. It is obtained through inversion of the leveling data that this earthquake was caused by a mainly reverse dip-slipping buried fault with strike 102°, dip 46° to SSW, rake 86° and a seismic moment of 9,4×10¹⁸ Nm. The stress drop, strain and energy released for this earthquake are estimated to be 4.9 MPa, 7.4×10^–5 and 7.0×10¹⁴ J, respectively. The slip distributes in a region slightly deep from NWW to SEE, with two nuclei, i.e., knots with highly concentrated slip, located in a shallower depth in the NWW and a deeper depth in the SEE, respectively.Broadband body waves data recorded by the China Digital Seismograph Network (CDSN) for the Gonghe earthquake are used to retrieve the source process of the earthquakes. It is found through moment-tensor inversion that theM _s =6.9 main shock is a complex rupture process dominated by shear faulting with scalar seismic moment of the best double-couple of 9.4×10¹⁸ Nm, which is identical to the seismic moment determined from leveling data. The moment rate tensor functions reveal that this earthquake consists of three consecutive events. The first event, with a scalar seismic moment of 4.7×10¹⁸ Nm, occurred between 0–12 s, and has a focal mechanism similar to that inverted from leveling data. The second event, with a smaller seismic moment of 2.1×10¹⁸ Nm, occurred between 12–31 s, and has a variable focal mechanism. The third event, with a sealar seismic moment of 2.5×10¹⁸ Nm, occurred between 31–41 s, and has a focal mechanism similar to that inverted from leveling data. The strike of the 1990 Gonghe earthquake, and the significantly reverse dip-slip with minor left-lateral strike-slip motion suggest that the pressure axis of the tectonic stress field in the northeastern Qinghai-Xizang plateau is close to horizontal and oriented NNE to SSW, consistent with the relative collision motion between the Indian and Eurasian plates. The predominant thrust mechanism and the complexity in the tempo-spatial rupture process of the Gonghe earthquake, as revealed by the geodetic and seismic data, is generally consistent with the overall distribution of isoseismals, aftershock seismicity and the geometry of intersecting faults structure in the Gonghe basin of the northeastern Qinghai-Xizang plateau.Contribution No. 96 B0006 Institute of Geophysics, State Seismological Bureau, Beijing, China.     

岩石圈结构对大地震震后形变的影响——以1976年唐山大地震和2001年昆仑山大地震为例

大地震发生之后通常会诱发一系列的余震序列,对比1976年M_S7.8唐山大地震和2001年M_S8.1昆仑山大地震周边区域的地震事件可以看出,唐山大地震余震活动时间要明显长于昆仑山大地震余震活动时间.余震序列往往与震后形变密切相关,而影响震后形变的因素不仅与地震发震断层和震级有关,同时与岩石圈的结构有关.考虑到唐山大地震的发震区华北地块和昆仑山大地震的发震区青藏高原有着较大的岩石圈结构差异,本文采用PSGRN/PSCMP软件计算了岩石圈分层模型的大地震同震和震后形变,分析了地壳弹性模量、弹性厚度以及黏滞性系数对同震和震后形变的影响,进而讨论了影响唐山地震和昆仑山地震余震序列差异的原因.计算结果显示,震后形变会在黏弹性效应的作用下逐渐调整,震后形变的持续时间与地壳弹性模量、地壳弹性厚度和下地壳黏滞性系数有关.上地壳和下地壳弹性模量越大,震后形变达到稳定值的时间越短,弹性模量对震后形变稳定值影响很小.地壳弹性厚度越大,震后形变达到稳定值的时间越短,当断层面底端深度小于地壳弹性厚度时,地壳弹性厚度的增加会引起震后形变稳定值的减小;下地壳厚度对震后形变达到稳定值的时间和稳定值基本无影响.下地壳黏滞性系数越大,震后形变达到稳定值的时间越长,反之亦然.结合唐山地震区的华北地块和昆仑山地震的青藏高原深部结构发现,两者之间的上地壳弹性模型差别不大,唐山地震区地壳弹性厚度略大于昆仑山地震区,但昆仑山地震区下地壳黏滞性系数明显低于唐山地震区.这些因素均决定了昆仑山地震的震后形变持续时间短（余震时间序列短）而唐山地震的震后形变持续时间长（余震时间序列长）.由此可见,岩石圈结构差异可能是导致唐山地震和昆仑山地震余震序列差异的主要因素之一.     

Fault deformation anomaly and intermediate and short-term prediction of the Jingtai M s=5.9 earthquake

The time-space distribution characteristics of fault deformation anomaly in the near-source region and its outlying zone in the seismogenic process of the Jingtai M _s=5.9 earthquake occurred on June 6, 2000 in Gansu Province is studied preliminarily. The distribution scope of fault deformation anomaly before the earthquake is wide, the anomaly shape is complicated and the pattern anomalous zone of fault deformation (strain) information index is obvious. The shape and amplitude of fault deformation anomaly in different regions differ significantly, which is closely related with the tectonic location of anomaly. The fault deformation anomaly of α, β, and γ phases along the western segment of Haiyuan fault zone shows the process from the quasi-linearity to non-linearity of fault movement in the near-source region, matches the high-value anomalous area of fault deformation (strain) information index, and reflects the high strain accumulation in the seismogenic region. However, the anomaly of abrupt jump and cusp with a large amplitude occurred in the areas far from the earthquake, such as Liupanshan fault zone which is the tectonic convergent section does not reflect the strain accumulation of its location, maybe it is a sign that the regional tectonic stress field is strengthened in the seismogenic process. Based on the above-mentioned facts and combined with the preliminary summary of experiences and lessons in the intermediate and short-term prediction of the Jingtai M _s=5.9 earthquake, we study and explore the application of fault deformation anomaly to earthquake judgment. Foundation item: National Key Basic Research Development Program (G1998040703 and G1998040705), and State Scientific and Technological Project of the “Ninth Five-Year Plan” (96-913-09-01-02-03 and 96-913-09-02-02-03), China.