Acceleration and evolution of faults: An example from the Hunter Mountain–Panamint Valley fault zone,Eastern California

We present new space geodetic data indicating that the present slip rate on the Hunter Mountain–Panamint Valley fault zone in Eastern California (5.0 ± 0.5 mm/yr) is significantly faster than geologic estimates based on fault total offset and inception time. We interpret this discrepancy as evidence for an accelerating fault and propose a new model for fault initiation and evolution. In this model, fault slip rate initially increases with time; hence geologic estimates averaged over the early stages of the fault's activity will tend to underestimate the present-day rate. The model is based on geologic data (total offset and fault initiation time) and geodetic data (present day slip rate). The model assumes a monotonic increase in slip rate with time as the fault matures and straightens. The rate increase follows a simple Rayleigh cumulative distribution. Integrating the rate-time path from fault inception to present-day gives the total fault offset.     

Effects of fault geometry and slip style on near-fault static displacements caused by the 1999 Chi-Chi,Taiwan earthquake

We investigated the fault geometry effects and the corresponding coseismic slip distribution using various proposed earthquake fault models for the Chi-Chi earthquake of 21 September 1999. The types of fault geometries are threefold: a simple planar fault plane, a two segmented planar fault plane and a three dimensional (3D) curved fault surface rupture propagation model. We derived the estimated spatial slip distribution from an inversion analysis of GPS coseismic displacement data and show that the 3D fault model is the preferred solution. The simple and segmented fault models lead to significant artificial slip distributions associated with the pre-defined fault geometry and the spatial distribution of GPS stations. The spatial distribution of coseismic slip deduced from the 3D fault model has three observable features: (1) the overall slip is concentrated at depth of less than 12 km, which may well correspond to a shallow-dipping detachment; (2) the maximum slip of about 10 m is located 45 km to the north of the epicenter; and (3) the slip vector is dominated by the dip-slip component. In addition, the results from the inversion of GPS data are consistent with those from the inversion analysis of teleseismic broadband data. A resolution analysis, further, demonstrates that the results are highly correlated with field GPS data studies when we used synthetic test data. The inversion of spatially distributed GPS data is highly sensitive to fault geometry. We conclude that the use of the 3D fault model is not only necessary but also certainly competent enough to well explain the inferred slip style and the observed static coseismic displacements.     

Co-seismic fault geometry and slip distribution of the 26 December 2004, giant Sumatra–Andaman earthquake constrained by GPS,coral reef,and remote sensing data

We analyze co-seismic displacement field of the 26 December 2004, giant Sumatra–Andaman earthquake derived from Global Position System observations,geological vertical measurement of coral head, and pivot line observed through remote sensing. Using the co-seismic displacement field and AK135 spherical layered Earth model, we invert co-seismic slip distribution along the seismic fault. We also search the best fault geometry model to fit the observed data. Assuming that the dip angle linearly increases in downward direction, the postfit residual variation of the inversed geometry model with dip angles linearly changing along fault strike are plotted. The geometry model with local minimum misfits is the one with dip angle linearly increasing along strike from 4.3oin top southernmost patch to 4.5oin top northernmost path and dip angle linearly increased. By using the fault shape and geodetic co-seismic data, we estimate the slip distribution on the curved fault. Our result shows that the earthquake ruptured *200-km width down to a depth of about 60 km.0.5–12.5 m of thrust slip is resolved with the largest slip centered around the central section of the rupture zone78N–108N in latitude. The estimated seismic moment is8.2 9 1022 N m, which is larger than estimation from the centroid moment magnitude(4.0 9 1022 N m), and smaller than estimation from normal-mode oscillation data modeling(1.0 9 1023 N m).     

基于InSAR同震形变观测反演2010年新西兰南岛Mw7.1 Darfield地震同震破裂分布

2010年9月4日新西兰南岛Canterbury平原发生了M_w7.1地震,震源深度约为10 km.本次地震发生在一条震前不为人所知的断层上.我们利用覆盖整个震区的合成孔径雷达(SAR)观测资料,通过干涉处理分析获得雷达视线向(LOS)同震形变场;以此资料为约束反演了断层的几何参数以及同震破裂分布.结果显示,该地震造成四条相对独立断层的破裂.大部分的地震矩释放发生在Greendale断层(编号1-4),其错动以右旋走滑为主,最大破裂约为8.5 m.其它三条断层中,经过震源的逆冲断层最大破裂为5.1 m (编号6),位于Greendale断层以西的逆冲断层最大破裂为3.5 m (编号5),位于Greendale断层北面的走滑断层最大破裂为1.9 m(编号7).反演的Greendale断层地表滑动与地质调查得到的地表破裂在形态和数值上均吻合较好.本次地震释放的地震矩为5.0×10¹⁹N·m,矩震级为7.1.板块边界带形变场分析表明,Darfield地震的发生受边界带应变分配在该地区残留构造应力场控制,其复杂性体现了区域构造应力场的特点.地震对其周围地区的应力场影响较大,库仑应力增加区与余震分布有一定对应关系,并在2011年Christchurch 6.3级地震发震断层区域造成约0.1bar的库仑应力增加,对此地震有一定的触发作用.     

The Grouping of the Strike-Slip Offset Data and Its Relation With the Seismicity Along Middle Segment of the Xiaojiang Active Fault

The data of the strike-slip offset along the Xiaojiang active fault can be obviously grouped.The groups of small orders of magnitude data within 100 m show clear linear characteristics of increments between 8 m and 12 m,which indicates that the segments of the Xiaojiang active fault is of characteristic seismicity and the distribution of the values of each group indicates that there are smaller earthquakes and creep between two large earthquakes along each segment of the Xiaojiang active fault.The interval between two characteristic large earthquakes can be calculated with the increments for two groups of slip data and the slip rate of the fault.Furthermore,the frequency of smaller earthquakes can also be estimated by comparing the distributions of the displacements of the large earthquakes with the distributions of the values of each group of data.The groups of large slip displacements show that there is close relationship between the records of the displacements of the fault and the changes of the cli     

Effect of non-uniform station coverage on the inversion for earthquake rupture history for a Haskell-type source model

In order to determine how reliably one can invert accelerograms to determine the rupture process details, when the station configuration is less than optimal, we use the vertical component of synthetic accelerograms for a Haskell-type earthquake rupture model, at stations in the vicinity of a dip-slip fault and solve the inverse problem. Of the various station configurations used, one is a uniform distribution and the others are very non-uniform. Faults of two different aspect ratios are considered. We mainly use much larger spatial and temporal cell sizes in the inversion than we use to construct the artificial data. The fault mechanism and the fault area are taken as known in the inversions. To solve the inverse problem, we use the method of linear programming and stabilize the solution by the use of physical constraints. The constraints of positivity of the slip rates on the fault is used in all cases in this study. In some cases, additional physical constraints such as preassigning the final moment, the rupture speed, and so on, are also used. We find that using a cell size almost double the wavelength of interest, we are able to reproduce the solution of the problem, even when we add a small amount of random noise to the artificial data, provided the source medium structure is known. We show that the best station configuration is when the stations are on the hanging wall, due to the fact that they provide the best illumination of the fault surface. This provides an incentive to install permanent ocean bottom strong ground motion stations in subduction zones. We also analyzed the effect of the rupture propagation direction on the results of the inversion showing that even four stations are sufficient to retrieve the rupture process if they are in the forward direction of the rupture propagation; the results for this case are better than when the four stations are placed in the backward direction, even when their positions are such that they illuminate the fault in exactly the same way as the four stations in the forward direction. Thus azimuthal distribution and the resulting illumination of the fault as well as the relation of the position of the stations to the direction of rupture propagation are more important than simply the number of stations. Finally, we find that proper knowledge of source medium structure is essential to recover the source process details reliably and that poor knowledge of crustal structure cannot be compensated by adding stations or by additional constraints.     

Combination of Precise Leveling and InSAR Data to Constrain Source Parameters of the Mw = 6.5, 26 December 2003 Bam Earthquake

We used new precise leveling data acquired 40 days after the Bam earthquake in combination with radar interferometry observations from both ascending and descending orbits to investigate static deformation associated with the 2003 Bam earthquake. We invert this geodetic data set to gain insight into the fault geometry and slip distribution of the rupture. The best-fitting dislocation model is a steeply east-dipping right-lateral strike-slip fault that has a size of 11 by 8 km and strikes N2°W. We find that such smooth geometry fits available geodetic data better than previously proposed models for this earthquake. Our distributed slip model indicates a maximum strike slip of 3 m occurring about 3 to 5 km deep. The slip magnitude and depth of faulting taper to the north, where the fault approaches the Bam city. Inclusion of crustal layering increases the amount of maximum slip inferred at depth by about 4%.     

青藏高原东北缘主要断裂闭锁特征和滑动亏损研究

利用1999—2015年GPS水平速度场, 基于块体-位错模型, 反演了青藏高原东北缘4条主要断裂(海原断裂, 六盘山断裂, 陇县—宝鸡断裂, 西秦岭北缘断裂)的闭锁程度和滑动亏损速率的空间分布, 并分析了各断裂的地震危险性。 结果显示, 六盘山断裂南段、 陇县-宝鸡断裂北段、 西秦岭北缘断裂东段闭锁程度最强, 闭锁深度达到24 km左右; 西秦岭北缘断裂东段滑动亏损速率最大, 平均值达到3 mm/a; 六盘山断裂南段、 陇县—宝鸡断裂北段滑动亏损速率平均值达到1.9 mm/a, 稍弱于西秦岭北缘断裂东段; 海原断裂闭锁程度和滑动亏损速率相对较小, 闭锁程度和滑动亏损都仅分布在浅部。 我们认为现阶段海原断裂的地震危险性相对较小, 六盘山断裂南段、 陇县—宝鸡北段、 西秦岭北缘断裂东段地震危险性高于这些断裂的其他段落。 这些结果对于青藏高原东北缘地震危险性判定和地震灾害评估具有参考意义。     

断层交汇区附近的变形特点与声发射特点的实验研究

把昌平花岗岩切割成20多种含两组断层交汇情况的模型.断层物质用2毫米厚的石膏模拟.模型受单轴应力的作用.实验中得出四种破裂组合情况:1.滑动破裂型 发生在两组断层与主压应力轴1的交角一组大于45,一组小于45的情形下.其中,当断层物质石膏湿度为10——15%时,先在大角度面上产生剪切破裂,后沿小角度面产生最终磨擦滑动.在应力位移曲线上表现为两次应力降.当断层物质石膏湿度小于5%时,则只沿小角度面上产生磨擦滑动.应力位移曲线为单应力降。2.压楔型 发生在两组断层与1夹角均大于45的情况下.其应力位移曲线为阶梯式上升型.断层位移方向反复.3.张楔型 发生在两组断层与1夹角均不大于30时.应力位移曲线为单应力降.4.双滑或双破型 发生在两组断层与1夹角均在30和45之间.应力位移曲线为连续的两次应力降.根据应力位移曲线和应力声发射累积数曲线的相似性,主破裂前的声发射累积数与主破裂时峰值应力的相似性以及断层组合与声发射组合的关系中,看出声发射过程对构造条件是敏感的.文中还对不同断层组合下可能出现的地震前兆进行了讨论.      

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

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

2000年日本鸟取县7.3级地震断层滑动分布的研究

利用Hori(2001)的反演方法,由同震GPS观测数据和水准观测数据,计算了2000年日本鸟取县西部地区7.3级地震断层滑动分布,并与余震分布进行了比较。结果显示,其最大滑动量达到5m,断层滑动分布主要集中在浅部,与余震分布并不一致。     

近断层应力变化快速计算方法及对强余震空间分布的指示意义:以M_w7.9汶川地震为例

构造地震一般由断层摩擦失稳所致.断层内部及周边所累积的剪切形变则通过同震滑动位移得到局部释放.因此,震后断层内部及近断层周边的静态剪切应力变化量的空间分布可通过断层面上的滑动位移分布计算得到.本文采用傅氏变换方法(FTM)计算单一有限断层同震滑移场所形成的静态剪切应力变化分布,近断层两侧的应力变化计算可由波数域内应力近似算法获得.结果表明,FTM快速有效、易于实现,有效地避免了常规应力计算中奇异值的出现.以2008年Mw7.9中国汶川大地震为例,采用前人所得有限断层滑动位移模型,得到了断层面和近断层周边准3D剪切应力分布解,并同主震后中强余震的空间分布特征作了比较.结果表明,大部分的中强余震震源位置处于剪切应力变化值为正的区域,由同震位移所产生的局部应力降峰值和均值大小同应力变化的正值大致相当,从而表明了快速且有效地计算断层内部及近断层附近的应力变化分布可以为主震后强余震发生的潜在区域提供指示意义.需要强调的是,应力变化空间分布特征的获取强烈地依赖于有限断层滑移模型解.有关滑动位移反演解的多解性对应力变化计算结果的影响,本文作了必要的讨论.     

Monte Carlo Inversion of DInSAR Data for Dislocation Modeling: Application to the 1997 Umbria-Marche Seismic Sequence (Central Italy)

— The study of surface deformation due to seismic activity is often made using dislocations with uniform slip and simple geometries. A better modeling of coseismic and postseismic surface displacements can be obtained by using dislocations with variable slip and nonregular shapes. This is consistent with the asperity model of fault surfaces, assuming a friction distribution on faults made of locked zones with much higher friction than surrounding zones. In this paper we consider the 1997–1998 Colfiorito seismic sequence. The coseismic surface displacements in the Colfiorito zone are used in order to infer the slip distribution on the fault surface at different stages of the sequence. The displacement field has been modeled varying the slip distribution on the fault, and comparing the deformation observed by SAR and GPS techniques with model results. The slip distribution is calculated by Monte Carlo simulations on a normal fault with the dip angle equal to 40°. A good approximation is obtained by using square asperity units of 1.5×1.5 km². In the first stage, we employed a simplified model with uniform slip, in which each asperity unit is allowed to slip a constant amount or not to slip at all, and in the second stage, we evaluate the slip distribution in the dislocation area determined by the Monte Carlo inversion: in this case we allow unit cells to undergo different values of slip in order to refine the initial dislocation model. The results show that the 1997 seismic events of the sequence can be modeled by irregular dislocations, obtaining a good fit to the DInSAR and GPS observations. The model also confirms the results of previous studies by a different methodology, defining the distribution of asperities on the fault plane using the fault geometry, the geodetic data and the seismic moment of the 1997–1998 Colfiorito seismic sequence. Furthermore, the analysis of 1997 aftershocks in the seismogenic region shows a strong correlation between most events and the asperity distribution, which can be considered as an independent test of the validity of the model.     

用形变资料反演1976年唐山地震序列的破裂分布

1976年唐山发生了7.8级地震,相继又发生了两次大余震——滦县7.1级地震和宁河6.9级地震.地震发生在观测条件比较好的地区,水准测量和三角测量测得了地震的同震位移场.本研究采用原始水准测量数据,而不是采用根据水准数据处理的地面沉降图像,和三角测量数据反演了该地震序列的破裂分布.模型构建中考虑了滦县地震和宁河地震的断层形态和大小.结果表明,唐山地震主震断层有明显的右旋走滑性质,最大走向滑动错距>6 m,位于断层南段,北段的走滑分量明显小于南段.主震总地震矩达2.58×10²⁰N·m,与地震波反演得到的地震矩的量级相当;滦县地震断层总体表现为左旋正断层,释放地震矩达4.95×10¹⁹N·m;宁河地震断层总体表现为右旋正断层,释放地震矩达3.94×10¹⁹N·m,比地震波反演的地震矩大一个量级.据此可以推测唐山地震的无震滑移主要发生在宁河地震断层的西部上,滑动性质以正断层为主.该结果对于唐山地震序列后的动力学演变过程及余震发生机理有一定参考.     

2016年台湾美浓MW6.4地震震源参数的InSAR和GPS反演

2016年2月6日台湾西南部高雄市美浓区发生了M_W6.4地震.本文结合ALOS2卫星升降轨、Sentinel-1A升轨SAR数据,采用两轨差分干涉技术获取了该区域的同震形变场,形变结果表明震中西北部以抬升为主,最大视线向形变量约为11.2 cm.基于均匀位错模型和多峰值粒子群（MPSO）算法,利用InSAR和GPS形变数据联合反演了美浓地震的断层几何参数,结果表明震源中心位于22.920°N,120.420°E,深度约12 km,发震断层长度约15 km,走向角307°,倾角16.5°,平均滑动角为51.5°,此次地震是以逆冲倾滑兼左旋走滑的破裂模式.利用格网迭代搜索法得到最优倾角为15.7°,GPS和InSAR最优权比为18：1,最优平滑因子为0.06.基于非均匀位错模型,利用非负最小二乘方法进行线性反演,结果显示最大倾滑和走滑量分别为51.7 cm和55.3 cm,对应矩震级为M_W6.38,略小于GCMT （M_W6.4）的结果.通过与已有文献的比较和对该区域断层构造的分析,发现美浓地震的发震断层为单一断层的解释更为合理,我们推测发震断层是位于左镇、后甲里等断层之间的一条东南-西北走向往东北倾斜的盲断层,并初步推测2010年M_W6.3甲仙地震也同该断层有关.     

联合采用定性和定量断层资料的应力张量反演方法及在乌鲁木齐地区的应用

地壳应力场是地球动力学研究和地震孕育环境研究的重要基础资料.求解应力场的一种重要方法是根据断层擦痕资料反演,然而擦痕往往难以观测,但断层滑动性质也与应力场有紧密的关系,这部分资料十分丰富.本研究按照应力张量在断层面上的剪应力方向与擦痕滑动方向及断层滑动性质一致为准则提出了联合采用擦痕数据和定性断层滑动数据求解应力场的网格搜索法,并给出了反演参数的一定置信度下的置信区间.采用假定的走滑型、挤压型和拉张型应力张量产生由7个擦痕数据和80个定性断层滑动数据组成的三组数据集.采用这些人工合成数据求解应力张量的实验结果表明,该方法可以更为准确地反演应力张量参数,定性断层滑动资料的参与使得应力张量的主应力方向更加接近假定的主应力方向,并且,应力张量参数的不确定性大大减小.将该方法应用于乌鲁木齐附近的定量擦痕和定性断层滑动性质数据,得到该地区较为精确的应力张量.结果表明,乌鲁木齐主压应力方向为近南北向、主张应力方向近垂直向,应力形因子为0.03.该地区占主导地位的近东西向的逆冲断层是近南北向挤压和近垂直向拉张形成的.较低的应力形因子表明乌鲁木齐地区在南北向挤压并辅以博罗科努山和博格达山的东、西向挤压的应力状态下处于隆升状态.     

Modeling coseismic displacements resulting from SAR interferometry and GPS measurements during the 1997 Umbria-Marche seismic sequence

In this study we analyse coseismic GPS displacements and DInSAR data to constrain a dislocation model for the three largest earthquakes of the 1997 Umbria-Marche seismic sequence. The first two events, which occurred on September 26 at 00:33 GMT (Mw 5.7) and 09:40 GMT (Mw 6.0) respectively, are investigated using both GPS displacements and DInSAR interferograms. We discuss and compare the results of previous studies which separately modeled a smaller subset of geodetic data. We provide a dislocation model for these two earthquakes which fits well both GPS and DInSAR data and agrees with the results of seismological and geological investigations. The first event consists of a unilateral rupture towards the southeast with a uniform dislocation. The strike, rake and dip angles are those resulting from the CMT solution. The second event consists of an unilateral rupture towards the northwest and a variable slip distribution on the fault plane. The strike and the rake are consistent with the CMT solution, but the dip angle has been slightly modified to improve the simultaneous fit of GPS and DInSAR data. While the second rupture (09:40 GMT) arrived very close to the surface, the fit to geodetic data shows that the first rupture (00:33 GMT) is deeper (2 km), despite the more evident surface geological effects. The analysis of new SAR interferograms allows the identification of a 5–6 cm additional displacement caused by the October 3 (Mw 5.2) and 6 (Mw 5.4) seismic events.We use data from a new DInSAR interferogram to model the displacement field of the Sellano earthquake of October 14, 1997. For this event significant GPS measurements were not available. We tested two different fault plane geometries: a blind, planar fault (top depth = 2.4 km), and a curved (listric) fault reaching the surface. The two models provide a generally similar fit to the data, and show that most of the slip was released at depths greater than 2.4 km along a gently dipping (40^°–45^°) fault surface. They also show that a unilateral rupture does not allow fitting the interferometric fringes since there is evident surface deformation to the northwest of the hypocenter. Moreover, we suggest that the concentration of high residuals in the southern part of our uniform slip model may in fact indicate a certain slip variability in this area.We conclude that, despite the moderate magnitudes and the lack of significant surface faulting, the space geodetic data allowed to constrain dislocation models giving new insights in the rupture process of the three largest events of the sequence.     

Ms8.0汶川地震断层的应力状态以及对余震危险性的影响

利用中国地震局在汶川地震前后对四川盆地以及龙门山断裂进行的水压致裂绝对应力测量数据、近断层强震记录、以及由美国USGS公布的包括地震矩和地震波能量等在内的远场震源参数解,从简单断层模型出发,应用地震能量分配原理和库仑摩擦准则,初步估算和判断了2008年Ms8.0汶川地震断层破裂过程和震源参数,以及滑移弱化模型下应力变化...     

利用InSAR技术研究2003年西藏班戈MS6.1地震发震构造

2003年7月7日发生在西藏与青海交界处的班戈M_S6.1地震, 由于缺乏余震分布等可靠资料, 其发震构造及其活动性质等问题一直认识不清。 几家机构利用远场波动资料给出的震源机制解差异很大。 本文利用ENVISAT卫星ASAR数据和D-InSAR技术, 计算获得了该地震LOS方向的同震形变场图像, 并以此为约束反演获得了该地震的断层几何参数和同震滑动分布。 结果表明, 班戈M_S6.1地震的发震构造为控制唐古拉山西边界的波涛湖—土门断裂, 其发震断层为走向161°的高角度右旋斜滑正断层, 破裂长度约10 km, 滑动量主要集中在3~7 km深度范围, 最大滑动量0.26 m, 矩震级M_W5.6。     

THE VITRINITE REFLECTANCE OF FAULT ROCKS FROM THE WENCHUAN EARTHQUAKE FAULT ZONE: CONSTRAINTS ON FRICTIONAL PROPERTIES OF THE FAULT DURING THE EARTHQUAKE

The temperature rise caused by frictional heating during seismic slip is able to indicate dynamic frictional properties of the seismic fault,which provides an approach to understand the dynamic process and energy budget of an earthquake.The residual indicators of frictional heating within the fault zone also can be taken as an evidence for seismic events.The vitrinite reflectance is a commonly-used geothermometer in the coal,oil and gas industries.It also has some potential applications in the studies of fault rock and fault mechanics.We studied vitrinite reflectance (VR) of fault rocks collected from surface outcrops of the Wenchuan earthquake fault zone in this paper.The measured data reveal that the VR of fault rocks are affected by fault motion,and there is a trend that the VR increases towards the fault core,which indicates the effects of frictional heating.The VR of fault rocks from the Bajiaomiao outcrop is much higher than those from the Shenxigou outcrop,which probably suggests the difference in fault activity at the two outcrops.Our study also suggests that systematic measurement of VR across the fault zone is helpful in identifying slip zones and determining their widths.From the VR measurement on an oriented specimen containing the slip surface of the Wenchuan earthquake from the Shenxigou outcrop,we observed anomalous high VR values in two black-colored slip zones of about 2mm in width near the slip surface.The numerical calculation shows that the maximum temperature rise on the fault plane near Shenxigou was probably less than 162℃ during the Wenchuan earthquake,which means the dynamic fault strength was quite low.These estimations are roughly in accord with the results from the high-velocity frictional experiments.