基于形变观测分析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级左右地震复发周期比较一致.在实际孕震能力判定的工作中,由于不同性质的断层段在同震过程中会表现更多的组合形式,断层发震能力判定结果存在更多的不确定性,但利用区域形变观测等资料给出震间断层运动特征的研究工作对于断层强震发震能力的判定具有非常重要的实际意义.     

Retardations in fault creep rates before local moderate earthquakes along the San Andreas fault system,central California

Records of shallow aseismic slip (fault creep) obtained along parts of the San Andreas and Calaveras faults in central California demonstrate that significant changes in creep rates often have been associated with local moderate earthquakes. An immediate postearthquake increase followed by gradual, long-term decay back to a previous background rate is generally the most obvious earthquake effect on fault creep. This phenomenon, identified as aseismic afterslip, usually is characterized by above-average creep rates for several months to a few years. In several cases, minor step-like movements, called coseismic slip events, have occurred at or near the times of mainshocks. One extreme case of coseismic slip, recorded at Cienega Winery on the San Andreas fault 17.5 km southeast of San Juan Bautista, consisted of 11 mm of sudden displacement coincident with earthquakes ofM _L =5.3 andM _L =5.2 that occurred 2.5 minutes apart on 9 April 1961. At least one of these shocks originated on the main fault beneath the winery. Creep activity subsequently stopped at the winery for 19 months, then gradually returned to a nearly steady rate slightly below the previous long-term average.The phenomena mentioned above can be explained in terms of simple models consisting of relatively weak material along shallow reaches of the fault responding to changes in load imposed by sudden slip within the underlying seismogenic zone. In addition to coseismic slip and afterslip phenomena, however, pre-earthquakeretardations in creep rates also have been observed. Onsets of significant, persistent decreases in creep rates have occurred at several sites 12 months or more before the times of moderate earthquakes. A 44-month retardation before the 1979M _L =5.9 Coyote Lake earthquake on the Calaveras fault was recorded at the Shore Road creepmeter site 10 km northwest of Hollister. Creep retardation on the San Andreas fault near San Juan Bautista has been evident in records from one creepmeter site for the past 5 years. Retardations with durations of 21 and 19 months also occurred at Shore Road before the 1974 and 1984 earthquakes ofM _L =5.2 andM _L =6.2, respectively.Although creep retardation remains poorly understood, several possible explanations have been discussed previously. (1) Certain onsets of apparent creep retardation may be explained as abrupt terminations of afterslip generated from previous moderate-mainshock sequences. (2) Retardations may be related to significant decreases in the rate of seismic and/or aseismic slip occurring within or beneath the underlying seismogenic zone. Such decreases may be caused by changes in local conditions related to growth of asperities, strain hardening, or dilatancy, or perhaps by passage of stress-waves or other fluctuations in driving stresses. (3) Finally, creep rates may be lowered (or increased) by stresses imposed on the fault by seismic or aseismic slip on neighboring faults. In addition to causing creep-rate increases or retardations, such fault interactions occasionally may trigger earthquakes.Regardless of the actual mechanisms involved and the current lack of understanding of creep retardation, it appears that shallow fault creep is sensitive to local and regional effects that promote or accompany intermediate-term preparation stages leading to moderate earthquakes. A strategy for more complete monitoring of fault creep, wherever it is known to occur, therefore should be assigned a higher priority in our continuing efforts to test various hypotheses concerning the mechanical relations between seismic and aseismic slip.     

Influence of 2008 Wenchuan earthquake on earthquake occurrence trend of active faults in Sichuan-Yunnan region

The Wenchuan earthquake coseismic deformation field is inferred from the coseismic dislocation data based on a 3-D geometric model of the active faults in Sichuan-Yunnan region. Then the potential dislocation displacement is inverted from the deformation field in the 3-D geometric model. While the faults' slip velocities are inverted from GPS and leveling data, which can be used as the long-term slip vector. After the potential dislocation displacements are projected to long-term slip direction, we have got the influence of Wenchuan earthquake on active faults in Sichuan-Yunnan region. The results show that the northwestern segment of Longmenshan fault, the southern segments of Xianshuihe fault, Anninghe fault, Zemuhe fault, northern and southern segments of Daliangshan fault, Mabian fault got earthquake risks advanced of 305, 19, 12, 9.1 and 18, 51 years respectively in the eastern part of Sichuan and Yunnan. The Lijiang-Xiaojinhe fault, Nujiang fault, Longling-Lancang fault, Nantinghe fault and Zhongdian fault also got earthquake risks advanced in the western part of Sichuan-Yunnan region. Whereas the northwestern segment of Xianshuihe fault and Xiaojiang fault got earthquake risks reduced after the Wenchuan earthquake.     

历史大地震断层滑动模型的建立及其对同震数值计算的影响——以1920年宁夏海原MS8.5大地震为例

大地震的发生会引起区域位移场和应力场发生变化,进而改变区域内及临近断层的应力状态和地震活动性.目前,研究学者可据已有的断层滑动模型来计算分析大地震同震应力变化,同时采用库仑应力触发理论来进一步分析震后余震分布和断层危险性.然而,历史上曾经发生过不少大地震,例如,1920年的海原M_S8.5大地震,是全球范围内少见的特大地震之一.局限于无确切的地震台站地震波等资料,前人在研究历史地震的影响时往往给出一些简单的断层滑动模型,将断层面上错动量视为均匀分布.为更准确地了解历史地震对后续地震的影响,基于前人研究和一般地震滑动形态分布规律及地震反射剖面等资料,以海原M_S8.5大地震为例,探讨了如何建立海原大地震断层滑动模型,并分别搭建了简单断层滑动模型和复杂断层滑动模型的全球同震横向不均匀并行椭球型地球模型.通过对海原M_S8.5地震的同震位移场和应力场的计算,发现采用复杂断层滑动模型比简单断层滑动模型地表位错分布更切合实际.同时,进一步计算和分析了此次大地震对青藏高原东北缘近100年历史地震和周围断层的应力触发作用,得出断层滑动模型对同震计算结果的影响集中在发震断层附近而对远场影响较小.     

2015尼泊尔M_S8.1地震中等余震震源机制研究

本研究利用西藏台网记录的波形数据,采用gCAP方法反演了2015年4月25日尼泊尔MS8.1大震5次中等余震(5.0≤MS≤6.5)及西藏定日MS5.9地震震源机制解.结果显示,6次地震包含2个正断、2个走滑及2个逆冲型地震.其中2个正断型地震位于主震的东北方向,即发震断层的上盘,表明该区域受到主震同震位移的影响,表现出应力拉张的变化特征;2个走滑型地震在主震破裂的东南方向上,说明随着破裂往东南方向延伸,余震的走滑分量增强;另外2个逆冲型地震位于5月12日MS7.5强余震区域,与MS7.5地震的滑移状态一致,可能与主震同震位移引起该区域处于应力挤压状态密切相关.这些结果表明,尼泊尔MS8.1主震发生后,由于同震位移的影响,不同区域处于不同的应力状态,从而使中等余震表现出不同的震源类型.     

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

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

Ductile creep and compaction: A mechanism for transiently increasing fluid pressure in mostly sealed fault zones

A simple cyclic process is proposed to explain why major strike-slip fault zones, including the San Andreas, are weak. Field and laboratory studies suggest that the fluid within fault zones is often mostly sealed from that in the surrounding country rock. Ductile creep driven by the difference between fluid pressure and lithostatic pressure within a fault zone leads to compaction that increases fluid pressure. The increased fluid pressure allows frictional failure in earthquakes at shear tractions far below those required when fluid pressure is hydrostatic. The frictional slip associated with earthquakes creates porosity in the fault zone. The cycle adjusts so that no net porosity is created (if the fault zone remains constant width). The fluid pressure within the fault zone reaches long-term dynamic equilibrium with the (hydrostatic) pressure in the country rock. One-dimensional models of this process lead to repeatable and predictable earthquake cycles. However, even modest complexity, such as two parallel fault splays with different pressure histories, will lead to complicated earthquake cycles. Two-dimensional calculations allowed computation of stress and fluid pressure as a function of depth but had complicated behavior with the unacceptable feature that numerical nodes failed one at a time rather than in large earthquakes. A possible way to remove this unphysical feature from the models would be to include a failure law in which the coefficient of friction increases at first with frictional slip, stabilizing the fault, and then decreases with further slip, destabilizing it.     

断层横向构造在逆冲型地震破裂中的作用——以汶川地震小鱼洞断层为例

大地震破裂大多由横向构造(如阶区、弯曲和分叉)所分割的多个段落组成.2008年5·12汶川地震破裂沿北东走向上穿过了多个横向构造部位,特别在震中北东45 km的位置,小鱼洞断层、北川断层和彭灌断层三者之间呈现复杂的断裂切割相交关系.复杂断层几何结构对破裂的扩展是有抑制还是促进的作用？在相交的断裂段之间是否存在最优的破裂顺序？本文以库仑应力分析为手段,探讨在汶川同震破裂初始30 s内,破裂在多分支断裂中选择扩展路径时的可能应力相互作用.库仑应力分析显示:如果北川断层先发生破裂,其滑动对小鱼洞断层和彭灌断层均产生强烈负应力的抑制作用,而彭灌断层的滑动却反而对小鱼洞断层和北川断层浅部有强烈正应力的促进作用.因此,从准静态应力分析角度,彭灌断层先于北川断层发生破裂的可能性较大,这一破裂顺序与小鱼洞断层参与同震破裂过程的事实相符.此外,小鱼洞断层在链接北川和彭灌断层的同震位移中可能起到桥梁作用,但非静态应力的影响.横向构造在逆冲型地震破裂扩展过程中起到的牵引作用使得逆冲型地震破裂能够比走滑型地震跨越更宽的阶区.横向构造是逆冲断裂带内广泛发育的构成单元,因此在地震危险性分析的最大潜在震级测算中应该考虑其作用.     

Instability model for recurring large and great earthquakes in southern California

The locked section of the San Andreas fault in southern California has experienced a number of large and great earthquakes in the past, and thus is expected to have more in the future. To estimate the location, time, and slip of the next few earthquakes, an earthquake instability model is formulated. The model is similar to one recently developed for moderate earthquakes on the San Andreas fault near Parkfield, California. In both models, unstable faulting (the earthquake analog) is caused by failure of all or part of a patch of brittle, strain-softening fault zone. In the present model the patch extends downward from the ground surface to about 12 km depth, and extends 500 km along strike from Parkfield to the Salton Sea. The variation of patch strength along strike is adjusted by trial until the computed sequence of instabilities matches the sequence of large and great earthquakes sincea.d. 1080 reported by Sieh and others. The last earthquake was theM=8.3 Ft. Tejon event in 1857. The resulting strength variation has five contiguous sections of alternately low and high strength. From north to south, the approximate locations of the sections are: (1) Parkfield to Bitterwater Valley, (2) Bitterwater Valley to Lake Hughes, (3) Lake Hughes to San Bernardino, (4) San Bernardino to Palm Springs, and (5) Palm Springs to the Salton Sea. Sections 1, 3, and 5 have strengths between 53 and 88 bars; sections 2 and 4 have strengths between 164 and 193 bars. Patch section ends and unstable rupture ends usually coincide, although one or more adjacent patch sections may fail unstably at once. The model predicts that the next sections of the fault to slip unstably will be 1, 3, and 5; the order and dates depend on the assumed length of an earthquake rupture in about 1700.     

A prototype earthquake warning system for strike-slip earthquakes

A prototype expert system has been developed to provide rapid warning of earthquakes while they are occurring. Warning times of up to 100 seconds will be possible. In the complete system, several accelerometers are distributed at intervals within a few kilometers of a known fault; data are telemetered to a central computer which implements the expert system. The expert system incorporates specific information about the type of fault to be monitored, and includes simple rules for estimating the fault slip, rupture length, and seismic moment, all in real time. If the seismic moment exceeds a preset value, an alarm may be issued. The prototype is designed for deployment on near-surface strike-slip faults such as the San Andreas and has been successfully tested with data from the 1979 Imperial Valley and 1984 Morgan Hill earthquakes. Crucial concepts have also been tested using synthetic data calculated for a model of the 1857 Fort Tejon earthquake. Parkfield, California, could be used as a test site.     

Visco-elastic stress triggering model of Tangshan earthquake sequence

We calculated the Coulomb failure stress change generated by the 1976 Tangshan earthquake that is projected onto the fault planes and slip directions of large subsequent aftershocks.Results of previous studies on the seismic fail-ure distribution,crustal velocity and viscosity structures of the Tangshan earthquake are used as model constraints.Effects of the local pore fluid pressure and impact of soft medium near the fault are also considered.Our result shows that the subsequent Luanxian and Ninghe earthquakes occurred in the regions with a positive Coulomb fail-ure stress produced by the Tangshan earthquake.To study the triggering effect of the Tangshan,Luanxian,and Ninghe earthquakes on the follow-up small earthquakes,we first evaluate the possible focal mechanisms of small earthquakes according to the regional stress field and co-seismic slip distributions derived from previous studies,assuming the amplitude of regional tectonic stress as 10 MPa.By projecting the stress changes generated by the above three earthquakes onto the possible fault planes and slip directions of small earthquakes,we find that the "butterfly" distribution pattern of increased Coulomb failure stress is consistent with the spatial distribution of follow-up earthquakes,and 95% of the aftershocks occurred in regions where Coulomb failure stresses increase,indicating that the former large earthquakes modulated occurrences of follow-up earthquakes in the Tangshan earthquake sequence.This result has some significance in rapid assessment of aftershock hazard after a large earthquake.If detailed failure distribution,seismogenic fault in the focal area and their slip features can be rapidly determined after a large earthquake,our algorithm can be used to predict the locations of large aftershocks.     

1997年以来巴颜喀拉块体周缘强震之间的黏弹性触发研究

本文利用考虑黏弹性地壳结构和精确震源参数的震后形变模型,计算了玛尼、昆仑山口西、汶川地震所形成的同震和震后形变场的变化过程与特征,结果显示玛尼地震的震后形变场有利于昆仑山口西地震的能量积累,玛尼地震与昆仑山口西地震的震后形变场对汶川地震的作用不明显,而这三次地震总的震后形变场在玉树断裂带附近形成了一个明显带有左旋性质的...     

A viscoelastic relaxation model for post-seismic deformation from the San Francisco earthquake of 1906

We present a realistic model of the San Andreas fault zone. We propose that aseismic ground displacement is a sum of visco-elastic relaxation following large earthquakes, transient fault slip, steady fault slip and a large-scale relative plate motion. We used the model to explain the aseismic ground displacements observed after the San Francisco earthquake of 1906.The data do not resolve the question of which is the dominant mechanism, but viscoelastic relaxation can contribute a significant fraction of the displacement if the elastic plate thickness is 50 km or less. If the relative plate motion is taken to be 5.5 cm/yr, as found from plate rotation pole studies, then the zone of significant shearing in the mantle is probably at least 100 km thick beneath California.     

Stress accumulation and release on the San Andreas fault

Summary The San Andreas fault can be divided into locked and free sections. On the locked sections accumulated slip is released in great earthquakes. On the free sections slip is occurring continuously either aseismically or during smaller earthquakes. Stress drops during earthquakes can be estimated from the ratio of short to long period amplitudes and from surface strain. Surface heat flow may provide an upper bound on the absolute stress. The failure or yield stress must reach a maximum at some depth on the fault. This maximum may occur in the near-surface brittle zone or deeper in the plastic zone of the fault. The historic distribution of seismic activity provides information on the stress level. The accumulation of strain and stress on the fault can be predicted using elastic theory. It is necessary, however, to include the viscous coupling of the lithosphere to the asthenosphere in order to fully model the problem.     

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.     

龙门山断裂带深部构造变形的黏弹性模拟及其与强震活动的关联性探讨

根据龙门山断裂带地区的主要构造特征,建立该地区的有限元模型,同时考虑地下深处的黏弹性蠕动和不同部位间的接触关系,模拟计算了研究区在强震轮回活动中的时间演化历程.模拟结果表明:龙门山断裂带深处的滑动速率比浅表的滑动速率大,龙门山断裂带周围是相对容易发生应变积累的地区,其5~19km深度也是高应力聚集成核区,随时间推移的应力集中程度加剧而引发强震.本模拟分析证实了重复地震观测所揭示的龙门山断裂带存在深浅活动速率差异的现象,这在一定程度上可以解释出乎预料的汶川MW7.9地震的孕育机理.综合分析研究提出:应充分利用重复地震这一天然的"地下蠕变计(subsurface creepmeter)"来探测深部构造变形的活动,为强震危险性分析提供必要的"原位(in situ)观测"约束信息.     

2014年鲁甸6.5级地震成因、破裂特征和余震分布特征的库仑应力作用

根据2014年鲁甸M_S6.5地震的区域构造和余震共轭分布特征,本文首先计算了1733年小江断裂带北段M7.75地震,1850年则木河断裂带M7.5地震和1974年马边M_S7.1地震对鲁甸M_S6.5地震震源机制解两个节面的黏弹性库仑应力作用,结果显示NNW向发生主破裂的包谷垴-小河断裂受到这3次地震,尤其是1850年M7.5地震明显的库仑应力作用,我们认为则木河断裂的高速左旋走滑运动以及7级以上强震的重复发生对于包谷垴-小河断裂的强震孕育和断裂演化方面具有促进作用;然后分析了鲁甸6.5级地震的共轭破裂与余震分布特征,并分别计算了两个共轭破裂面单独破裂对另一破裂面的库仑应力作用,结果显示NEE向破裂促进NNW向破裂的发生,而NNW向破裂后则阻碍了NEE向破裂的进一步发展,最终发展成以NNW向破裂为主的共轭破裂事件;最后计算了共轭破裂产生的库仑应力变化对余震的影响,认为位于NEE向破裂面西侧的余震集中分布主要是由于应力触发作用而形成.     

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观测的1999-2007年汶川震前3期地表变形数据和2008年汶川同震地表变形数据,结合地震位错理论,通过高斯变换和坐标旋转建立断层模型,运用遗传算法,反演了龙门山断裂带断层震前3期和同震滑动参数。结果表明龙门山断层震前3期平均走滑位移为-5.39mm,倾向位移为2.66mm,与同震断层滑移相比较,发现震前断层的滑移趋势与同震断层滑移一致,均为逆冲兼右旋的挤压运动。比较震前3期逆冲方向的滑移量,发现逆冲滑移有加速的现象。并根据震前和同震的断层滑动量估算了汶川地震复发周期。     

2001年昆仑山口西MS8.1地震对周围断层的应力影响数值分析

大地震的发生往往会引起周围区域形变场和应力场变化,且对临近断层上的应力状态也有影响.2001年11月4日,昆仑山口西发生了半个世纪以来中国最大的M_S8.1级地震.本文基于已有的滑动模型,建立了三维含地形高程的横向不均匀性椭球型地球有限元模型,采用等效体力方法,分析了此次M_S8.1地震产生的全球同震位移和应力场变化.与解析方法相比,该模型考虑了地形、Moho面起伏和地球介质横向不均匀性;与一般的有限元数值模拟相比,该模型考虑了地球曲率和椭率,合理地规避了有限块体模型假定边界位移为零所引入的误差.计算得出同震位移与GPS观测数据可以很好地吻合.据库仑破裂应力准则和震源参数,计算得出昆仑山口西M_S8.1地震的发生造成了汶川、芦山、改则和当雄地震的发震断层上库仑应力增加,对这些地震的发生起促进作用;而造成玉树和德令哈地震发震断层上的库仑应力变化为负值,在一定程度上抑制了这些断层的地震活动性.此外,计算结果显示地球地形高程、介质非均匀性和椭率对昆仑山口西M_S8.1地震同震变化计算有一定的影响,其中地形和椭率造成的同震位移场相对误差约10%.