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

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

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.     

2010年马乌莱MW8.8地震震后形变三维黏弹性数值模拟

2010年智利马乌莱M_W8.8地震发生在纳斯卡板块与南美板块的板块边界处,引起了显著的同震和震后效应。GPS台网数据显示记录到的同震海向位移最大约5 m,垂向沉降最大约50 cm。在经过对俯冲效应、季节变化等效应的校正后,震后6年的海向最大位移约68 cm,垂向抬升最大约20 cm。马乌莱地震显著的震后形变对该区域的地下三维黏弹性结构有良好的约束。本文建立了智利中南部俯冲带区域的三维有限元模型,黏弹性的地幔楔及海洋地幔均使用伯格斯体材料,并在断层面上设置2 km厚的软弱层以模拟震后余滑。在与GPS台站震后位移数据进行比较后,模拟结果表明,大洋地幔顶部存在约120 km厚,黏度为1×1019 Pa·s的软流圈。模拟震后余滑效应的软弱层黏度为5×1017 Pa·s,其等效震后余滑的最大值在震后前两年接近2 m,且随着时间的增长而快速衰减。      

利用GPS观测资料分析2008年于田Ms7.3地震的同震位移及震后形变

2008年3月21日新疆于田发生Ms7.3级地震.本文通过处理、分析GPS数据,得到破裂断层北侧100 km附近的同震位移及震后形变信息.在观测区域GPS点监测到10 mm左右的同震位移,其中最大为南向14 mm,东向5 mm.同震位移呈现一致性的东南向运动特征,证实于田地震存在显著的左旋走滑分量.震后台站向西南方向运...     

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.     

基于GPS观测研究2010年青海玉树MS7.1地震震后地壳形变特征及其机制

2010年4月14日青海玉树M_S7.1地震发生在青藏高原东南部甘孜-玉树地震带,在震后7~10天内,我们快速建立了由15个GPS测站组成的跨地震破裂带观测剖面,包括1个连续站,3个半连续站和11个流动站,对所有站进行了240多天的观测,获取了该次地震的震后形变时空特征.采用欧拉矢量和位错模型解算了背景速度场,并从GPS观测的形变场中扣除该分量.采用分层黏弹性位错模型计算余震引起的地表形变,结果表明余震对部分测站的位移造成不可忽视的影响.采用对数模型拟合位移时间序列,表明特征衰减时间为6.7±1.2天.利用最速下降法反演震后余滑时空分布,反演结果表明震后断层活动以左旋滑动为主,断层南盘具有少量的抬升.在空间分布上,余滑主要位于同震破裂区的两侧,西北侧的余滑几乎达到地表,而东南区的余滑基本在同震破裂区的下方,余滑最大的区域位于结古镇东南下方10~20 km的深度范围.随着震后离逝时间的增加,2个余滑区在空间上保持不变,余滑区的面积逐渐扩大.余滑的矩释放为（1.5~5.1）×10¹⁸Nm,相当于1个M_W6.1~6.4地震释放的能量.分层岩石圈黏弹性模型计算的地壳孔隙弹性反弹形变与地表观测值相差较大,不能解释观测到的震后变形.采用麦克斯维尔流变体模型计算下地壳和上地幔松弛引起的地表形变,显示出其对地表形变的贡献较小.GPS观测得到的震后形变所具有的快速衰减特征,以及余滑模型能够较好地拟合GPS地表形变,表明2010年玉树M_S7.1地震后早期阶段的地壳形变主要是由余滑机制决定的.     

Postseismic Deformation Following the 2010 $$M = 7.2$$ El Mayor-Cucapah Earthquake: Observations,Kinematic Inversions,and Dynamic Models

Due to its location on a transtensional section of the Pacific-North American plate boundary, the Salton Trough is a region featuring large strike-slip earthquakes within a regime of shallow asthenosphere, high heat flow, and complex faulting, and so postseismic deformation there may feature enhanced viscoelastic relaxation and afterslip that is particularly detectable at the surface. The 2010 \(M = 7.2\) El Mayor-Cucapah earthquake was the largest shock in the Salton Trough since 1892 and occurred close to the US-Mexico border, and so the postseismic deformation recorded by the continuous GPS network of southern California provides an opportunity to study the rheology of this region. Three-year postseismic transients extracted from GPS displacement time-series show four key features: (1) 1–2 cm of cumulative uplift in the Imperial Valley and \(\sim\)1 cm of subsidence in the Peninsular Ranges, (2) relatively large cumulative horizontal displacements \(>\)150 km from the rupture in the Peninsular Ranges, (3) rapidly decaying horizontal displacement rates in the first few months after the earthquake in the Imperial Valley, and (4) sustained horizontal velocities, following the rapid early motions, that were still visibly ongoing 3 years after the earthquake. Kinematic inversions show that the cumulative 3-year postseismic displacement field can be well fit by afterslip on and below the coseismic rupture, though these solutions require afterslip with a total moment equivalent to at least a \(M = 7.2\) earthquake and higher slip magnitudes than those predicted by coseismic stress changes. Forward modeling shows that stress-driven afterslip and viscoelastic relaxation in various configurations within the lithosphere can reproduce the early and later horizontal velocities in the Imperial Valley, while Newtonian viscoelastic relaxation in the asthenosphere can reproduce the uplift in the Imperial Valley and the subsidence and large westward displacements in the Peninsular Ranges. We present two forward models of dynamically coupled deformation mechanisms that fit the postseismic transient well: a model combining afterslip in the lower crust, Newtonian viscoelastic relaxation in a localized zone in the lower crust beneath areas of high heat flow and geothermal activity, and Newtonian viscoelastic relaxation in the asthenosphere; and a second model that replaces the afterslip in the first model with viscoelastic relaxation with a stress-dependent viscosity in the mantle. The rheology of this high-heat-flow, high-strain-rate region may incorporate elements of both these models and may well be more complex than either of them.     

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.     

GPS measurements of postseismic deformation due to October 8, 2005 Kashmir earthquake

Relaxation of the coseismic stresses following an earthquake causes postseismic crustal deformation, which can last for days to years. Continuous monitoring of postseismic deformation facilitates the understanding of the mechanism of deformation and postseismic relaxation and viscous rheology. After the October 8, 2005 Kashmir earthquake, global positioning system data for 8 months, starting from October, 2005 have been analyzed from three continuous sites located at Gulmarg, Amritsar, and Jaipur. The average velocity during the observation period at Gulmarg (8.6 cm/year) is significantly higher than the Indian plate velocity exhibiting postseismic crustal deformation. The velocity at Amritsar (5.9 cm/year) and Jaipur (5.1 cm/year) is comparable to the Indian plate velocity. At Gulmarg, the logarithmic function fits well to the north–south component of postseismic transients (~in the coseismic slip direction). The nature of decay in these transients suggests that the deformation is mainly due to an afterslip, and the second possible contribution may be from the viscous relaxation process. This paper presents the characteristics of postseismic transients and possible contributions from various postseismic mechanisms subsequent to the Kashmir earthquake.     

Spatio-temporal Slip, and Stress Level on the Faults within the Western Foothills of Taiwan: Implications for Fault Frictional Properties

We use preseismic, coseismic, and postseismic GPS data of the 1999 Chi-Chi earthquake to infer spatio-temporal variation of fault slip and frictional behavior on the Chelungpu fault. The geodetic data shows that coseismic slip during the Chi-Chi earthquake occurred within a patch that was locked in the period preceding the earthquake, and that afterslip occurred dominantly downdip from the ruptured area. To first-order, the observed pattern and the temporal evolution of afterslip is consistent with models of the seismic cycle based on rate-and-state friction. Comparison with the distribution of temperature on the fault derived from thermo-kinematic modeling shows that aseismic slip becomes dominant where temperature is estimated to exceed 200° at depth. This inference is consistent with the temperature induced transition from velocity-weakening to velocity-strengthening friction that is observed in laboratory experiments on quartzo-feldspathic rocks. The time evolution of afterslip is consistent with afterslip being governed by velocity-strengthening frictional sliding. The dependency of friction, μ, on the sliding velocity, V, is estimated to be ${{\partial \mu }/{\partial \, {\rm ln}\, V}} = 8 \times 10^{ - 3}$ . We report an azimuthal difference of about 10–20° between preseismic and postseismic GPS velocities, which we interpret to reflect the very low shear stress on the creeping portion of the décollement beneath the Central Range, of the order of 1–3 MPa, implying a very low friction of about 0.01. This study highlights the importance of temperature and pore pressure in determining fault frictional sliding.     

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

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

2001年MW7.8昆仑山地震震后形变观测、机制与青藏高原中北部岩石圈流变:认识与挑战

2001年M_W7.8昆仑山地震是近半个世纪以来青藏高原发生的最大震级地震。同震破裂产生的巨大应力扰动驱动控制着显著震后形变。二十年尺度的大地测量数据记录了地震后长时间、大范围、时空依赖的震后形变演化过程及差异,揭示了昆仑山地震破裂段复杂的断层分段震后运动学特征、分段摩擦性质差异和青藏高原中北部岩石圈流变性质/结构横向各向异性。本文简要回顾昆仑山地震后基于二十年尺度时序InSAR和GPS的震后形变观测方法和时空特征,特别是时空密集的InSAR观测,是该构造区震后GPS观测的重要补充及其不可替代的观测手段;总结大范围震后形变模拟方法及其揭示的震后运动过程、多种动力学机制及其关系。最后总结提出昆仑山地震震后形变20年研究取得的科学认识及尚待深入研究的科学问题,一方面要持续性地对东昆仑断裂带大范围地表形变进行观测研究;另一方面,要不断更新震后形变模型,进一步深化对该断裂带地震周期形变、区域构造对周期形变控制作用、复杂断层运动时空演化机制的认识。      

丽江M_S7.0地震余震深度揭示出的中地壳脆塑性转化特征

1996年丽江MS7.0地震的余震深度分布明显具有时间依赖性,主震发生后短时间内余震震源深度较深,随着时间的延续,余震震源深度变得越来越浅。余震的这种深度分布受地壳脆塑性转化带深度控制,而脆塑性转化带的深度变化与地震前后断层的应变速率有关。由震后GPS地表变形数据得到的地表变形模型表明,震后地表变形主要来自地壳深部,震后滑动与地壳深部弹性松弛有关。根据鲜水河断层地表的滑动数据和按Marone's(1991)给出的方程确定的震后滑动模型,估计的应变速率显示,主震发生后应变速率较高,随时间延续,应变速率逐渐下降。基于地壳P波速度结构和利用热流数据估计的丽江地区地壳温度,采用含水石英的塑性流变参数,估计了中地壳脆塑性转化带深度随震后应变速率的变化。结果表明,主震震源深度与余震深度分布下限与中地壳脆塑性转化带的深度随时间变化趋势一致。由于断层的震后快速滑动致使断层带深部具有很高的应变速率,高应变速率引起断层脆塑性转化带深度下移,主震之后短时间内发生了较深的余震;随着震后时间的延续,断层逐渐进入蠕变阶段,断层滑动速率逐渐减小,地壳应变速率逐渐降低,断层脆塑性转化带也逐渐恢复到间震期的深度,相应余震深度随之变化。因此,余震分布的深度变化是中地壳流变结构和脆塑性转化带深度变化的直接反映。     

基于迭代PCA的GPS时间序列震后形变估计方法和应用

GPS时间序列的震后形变分析对于研究区域震后形变机制和岩石圈流变学特性以及维持国际动态地球参考框架具有重要意义.本文在现有参数估计方法的基础上,提出兼顾震后形变衰减特征空间相关性和整体建模的"迭代PCA参数估计方法",并利用模拟数据证实了新方法可以获取更稳健可靠的震后形变、同震形变和震间速度参数.最后,以37个新西兰GPS连续站坐标时间序列为例,利用迭代PCA方法提取了2016年11月13日Kaikoura地震共性的震后形变时间演化过程和各站点的震后形变,并定量分析了震后形变对地表速度的影响.结果表明各站的震后形变在时间域上以衰减常数τ为4天的对数模型持续松弛;空间域上南岛北东部和北岛最南部震后形变较大,其中,最大震后形变点为cmbl站,截至2017年6月10日NEU方向累计的震后形变分别达到107mm,135mm和187mm,地表速度分别达到133.58mm·a-1,112.05mm·a-1和175.58mm·a-1,仍高于稳定的震间速度.     

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

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

Modeling transient and long-term postseis-mic deformation: A case study of 1995 M_w9.5 Chile earthquake

In this paper, we firstly use finite element method (FEM) with Burgers model to simulate the postseismic viscoe-lastic relaxation taking 1960 Chile earthquake as an example. The postseismic deformation modeled with Burgers model includes co-seismic deformation, transient postseismic deformation and long-term postseismic deformation. So if we apply Burgers model to calculate postseismic deformation of 1960 Chile earthquake, there is no discrep-ancy phenomenon due to different durations of postseismic deforma...     

Modeling transient and long-term postseismic deformation: A case study of 1995 M W9.5 Chile earthquake

In this paper, we firstly use finite element method (FEM) with Burgers model to simulate the postseismic viscoe-lastic relaxation taking 1960 Chile earthquake as an example. The postseismic deformation modeled with Burgers model includes co-seismic deformation, transient postseismic deformation and long-term postseismic deformation. So if we apply Burgers model to calculate postseismic deformation of 1960 Chile earthquake, there is no discrep-ancy phenomenon due to different durations of postseismic deformations that happens in Maxwell model.     

基于2001年MW7.8可可西里地震震后形变模拟研究藏北地区岩石圈流变学结构

青藏高原岩石圈的流变学结构和形变机制是地学界长期争论的重大科学问题.2001年发生在东昆仑断裂带的M_W7.8可可西里地震造成青藏高原北部地区岩石圈构造应力场的很大改变,引起下地壳与上地幔的快速弛豫形变,从而为研究这一问题提供了难得的机会.本研究采用该区域的GPS震后观测,反演这一地区岩石圈的流变学参数并探讨其形变机制.反演所采用的数据来自45个GPS观测点,其中包括一个中国地壳运动观测网络的基准站,数据最长时间跨度达6.4年.大地震震后形变场主要来源于地壳、上地幔的黏弹性松弛与断层面上的震后余滑,因此本研究同时反演介质的黏滞系数和断层的震后余滑.考虑到东昆仑断层南侧的巴颜喀拉-羌塘地区与北侧的柴达木盆地地区具有明显不同的地壳结构,断层南北两侧采用不同的Burgers体流变学结构,其下地壳-上地幔的短期和长期黏滞系数采用网格搜索法获得;断层震后余滑反演则同时施加近似正比于库仑应力的约束.最终结果显示：东昆仑断层北侧柴达木盆地地区下地壳-上地幔短期和长期黏滞系数分别为5×10¹⁸ Pa·s和1.5×10²⁰ Pa·s;东昆仑断层南侧巴颜喀拉-羌塘地区下地壳-上地幔短期和长期黏滞系数分别为1.5×10¹⁸ Pa·s和1.5×10¹⁹ Pa·s.这一结果表明：巴颜喀拉-羌塘地区下地壳-上地幔黏滞系数显著低于柴达木盆地,意味着巴颜喀拉-羌塘地区下地壳可能存在部分熔融,其地壳形变模式更趋近于连续形变,而柴达木盆地形变模式更趋近于块体运动.研究区下地壳长期黏滞系数比下地壳流模型所主张的黏滞系数高2~3个数量级,表明下地壳流在本地区可能不存在.     

Long aseismic slip duration of the 2006 Java tsunami earthquake based on GPS data

The Java earthquake occurred on July 17,2006 with magnitude 7.8 associated to the subduction process of Indo-Australian plate and Sundaland block off southwestern coast of Java. We present postseismic deformation parameters of the 2006 Java earthquake analyzed using campaign GPS observation from 2006 to 2008 and continuous observation from 2007 to 2014. We use an analytical approach of logarithmic and exponential functions to model these GPS data. We find that the decay time in the order of hundreds of days after the mainshock as observed by 8 years' data after the mainshock for magnitude 7 earthquake is longer than a general megathrust earthquake event. Our findings suggest that the 2006 Java earthquake which is considered as ‘‘tsunami earthquake' most probably occurred in the region that has low rigidity and tends to continuously slip for long time periods.     

大地测量约束下的阿尔泰山岩石圈流变结构

本文根据大地测量数据得到过去50年左右阿尔泰山富蕴区域形变场，破裂带中段的相对位移最大，平均速率达6mm/a，总体上表现为沿断裂的走滑运动特征．该形变场可以用1931年富蕴8级地震的震后黏弹性松弛模型进行模拟，反演得到下地壳黏滞系数为1.6×10¹⁹～7.9×10¹⁹Pa·s，上地幔黏滞系数为16×1018～63×1019Pa·s，与华北、Nevada等地区利用震后变形资料推算的黏度基本一致．根据该地区最佳黏弹性分层模型，最近五十年由于岩石圈下部应力松弛引起的地震破裂带两侧最大水平速率约为4mm/a． 我们的研究表明：大陆7～8级大震在几十年后仍可能有可观的地表变形，GPS监测得到的现今变形场可能包含震后变形成分．