中国大陆地壳水平运动速度场与应变场

收集了中国大陆及周边地区GPS网的有关数据，提出了GPS网速度场的不同融合方法；经过融合建立了中国大陆及周边地区统一的地壳运动速度场，该速度场使用的有效GPS站共423个，其覆盖面积为1200万km^2；初步总结出中国大陆及周边地区地壳水平运动空间分布的基本特征；建立了板内块体的刚性弹塑性运动应变模型，对其进行了块体应变参数唯一性与速度残差中误差最小检验；根据中国大陆及周边地区的速度场，估计了8个块体的应变参数，分析了这些块体的应变状态，估计出的各个块体的应变状态与地质学、地球物理学方法估计的结果具有很好的一致性。用喜马拉雅块体主压应变方向估计的印度板块向欧亚板块碰撞力的主方向为北东7．1度。     

板内块体的刚性弹塑性运动模型与中国大陆主要块体的应变状态

建立了板内块体的刚性弹塑性运动应变模型，并对其进行了块体应变参数唯一性与速度残差中误差最小检验．根据中国大陆及周围地区的速度场，估计了8个块体的应变参数，分析了这些块体的应变状态．本文估计的各个块体的应变状态与地质学、地球物理学方法估计的结果具有很好的一致性．由喜马拉雅块体主压应变方向估计的印度板块向欧亚板块碰撞力的主方向为北东7.1．      

由GPS资料讨论中国大陆块体运动与应变模型的建立及检验

目前很多文献基于弹性力学假设,以小形变、均匀应变场为前提,利用GPS资料建立了中国大陆块体的整体旋转与均匀应变模型,并对模型进行了相关检验.但上述对模型的检验侧重从模型的拟合效果上进行分析和比较,忽略了对模型的适用性和假设条件合理性的检验.实际上,不同块体由于自身构造和外部条件的不同,其运动状况与应变状态并不完全一致,模型的适用性需进一步的验证.本文提出了两步检验方法对所建立的模型在中国大陆6个主要块体上进行了检验：第一步对应变参数进行显著性检验,第二步对应变均匀性进行检验,结果表明：根据所用数据精度,华南块体的刚性程度在95% 的置信水平上可以认为显著的,建立模型时不必考虑内部形变;其他块体内部形变显著,整体旋转与均匀应变模型并不能反映真实的地壳运动状态,必须考虑其形变的不均匀性.     

基于GPS速度场采用RELSM模型分析青藏块体东北缘的形变—应变特征

利用青藏块体东北缘地区1999～2001年GPS观测获得的地壳水平运动速率场,通过对该地区进行块体划分,将该地区划分为9个块体,应用块体的整体旋转线性应变模型(RELSM)估计了各个块体的旋转与应变参数,以及计算了该地区内143个GPS站点的应变参数,以此分析了该地区的应变场的基本特征,结果表明:①阿拉善块体s较稳定,其旋转角为0.630×10-8,运动速率为0.688 mm/a,②相比其他块体,共和块体旋转角最大达到了6.589×10-8 ,运动速率达到了7.296 mm/a,③应变高值区主要集中在祁连山断裂,海原断裂等,在这些地区最大剪应变率达到了7.5×10-8、面膨胀率达到了-2.5×10-8、主压应变达到了-6×10-8.     

巴颜喀拉块体北东地区现今水平运动与变形

本文利用GPS数据研究了巴颜喀拉块体北东地区现今水平运动与变形特征。 在球坐标系中解算了各应变分量, 分析了应变率场的空间分布特征, 并与地球物理学和地震地质学研究结果进行了综合对比分析。 最新的GPS速度场结果表明, 巴颜喀拉块体北东地区与高原整体运动性质一样具有顺时针向南东方向旋转的特征, 自西向东和北东方向测站水平运动速度呈现明显的衰减特征。 应变场结果显示, 研究区以北东向的主压应变为主, 伴随着近北西向的张性应变。 应变较强的区域主要分布在活动块体的边界断裂东昆仑断裂带的东段塔藏段和龙门山断裂带上。 东昆仑断裂带东段塔藏段的主压应变明显, 结合地震地质和活动构造资料, 认为东昆仑断裂带东段塔藏段的运动性质自西向东发生了改变, 水平滑动速率逐渐减小, 垂向运动逐渐增强。 研究区GPS速度场和应变场的这一变形特征表明, 青藏高原内部的块体运动特征较为明显, 变形主要集中在作为活动块体边界的活动断裂带上, 边界断裂带的运动特征在调节活动块体间的相互运动中起着重要作用。     

昆仑山口西8.1级地震前后山西及邻区速度场和应变场的变化

在山西及邻区,1999~2004年几乎所有的GPS站速度都小于1999~2001年,鄂尔多斯西部大多数GPS站1999~2004年的运动方向相对于1999~2001年发生了逆时针旋转。鄂尔多斯和山西中南部的NWW向运动与昆仑山口西8.1级大震断裂带以北地区的运动方向是一致的,而鄂尔多斯以南的向东运动与断裂带以南地壳的运动方向也是一致的。鄂尔多斯的逆时针旋转形成了109.0°E以西地区E-W方向的压应变。而它的逆时针旋转和太行山东部的向东伸展运动形成了109.0°E以东地区E-W方向的张应变。鄂尔多斯东部和山西地区在8.1级大震后的增量应变场减弱了基本应变场E-W方向压应变的积累,而鄂尔多斯西部的增量应变场增加了E-W方向压应变的积累。     

玉树MS7.1级地震前的形变应变场变化特征

基于1999～2007年和2007～2009年的GPS站点观测数据计算的速度场,分析了青藏亚板块的运动特征.利用块体(旋转与线性应变)弹性运动公式计算可可西里-巴颜喀拉块体和羌唐块体的应变率、可可西里-玉树断裂带的运动速度和应变率,研究了玉树7.1级地震前的形变应变特征.结果表明:可可西里-巴颜喀拉块体处于青藏亚板块运...     

台湾地区地壳形变的弹性块体位错模型

在经典的非震形变位错模型中,地壳形变被认为是活动块体刚性运动和上部断层锁定影响的叠加,本文对此模型进行了改进: (1) 用活动块体整体运动和内部线性应变、旋转的贡献代替活动块体刚性运动的贡献;(2) 用分层介质地壳模型代替半无限介质模型计算断层锁定的影响. 利用改进后的非震形变位错模型,拟合了台湾地区1990～1995年间GPS观测资料. 结果显示,在东部海岸山脉区,约有30 mm·a^-1的汇聚率被奇美断层消耗掉,运动速度从奇美断层向北迅速衰减. 在西部平原地区,南部断层是岛内锁定最为强烈的断层,该地区相应的也是史上灾害性地震多发的地区. 根据反演结果计算出的应变率与旋转率分布与前人结果在大部分地区一致,主应变率场显示台湾大部分地区存在近NW-SE方向的主压应变,主压应变方向呈扇形分布. 旋转率场显示台湾东部和南部地区存在着逆时针旋转率,而西部和北部地区则为顺时针旋转率.     

三峡重庆库区现今地壳形变特征分析

基于三峡重庆库区2008-2012年GPS观测资料,利用块体的整体旋转线性应变模型,首次获取了该地区的地壳水平运动及应变场特征.结果表明:三峡重庆库区主应变以压应变为主,方向近似为北西向;最大剪应变高值区集中分布在该地区的华蓥山断裂与金佛山断裂,量值均达到1.4×10-8;区域主压应力方向与大范围构造应力场基本一致,但方位角在华蓥山断裂与金佛山断裂两侧具有明显偏转的特征,反映出上述地区是三峡重庆区内的地壳构造活动区域.     

中国大陆活动地块的运动与应变状态

从地壳运动与应变的角度给出了活动地块的定义,根据中国大陆及周边地区最近几年GPS观测得到的由1598个GPS站速度组成的统一速度场,估计了各个活动地块的运动与应变参数,分析了各个活动地块的运动与应变状态。中国大陆各地块存在一致的向东运动分量,但其南北分量是不一致的。西部地块存在一致的向北运动分量,东部地块存在一致的向南运动分量。在90°E以东,从喜马拉雅地块向NE方向,各地块的运动方向按顺时针方向旋转,各地块的运动速率是不相同的。从总体上看是西部大、东部小,南部大、北部小,西部大约是东部的3～4倍。各地块主压应变方向的空间分布是不相同的。在90°E以西各地块的主应变方向基本上为SN向,在青藏高原的东北部各地块的主压应变方向基本为NE向,在青藏高原东南部各地块的主压应变方向绕喜马拉雅构造东端顺时针方向旋转。各地块的主应变与剪应变率也是不同的,其中喜马拉雅、天山地块的主压和最大剪应变率最高,其次是拉萨、羌塘、滇西南、祁连与川滇地块。东部各地块的应变率较小。根据应变状态推测,喜马拉雅地块南北向的缩短速率为(15.2±1.5)mm/a,仍然是现今构造活动最强烈的地区,其次是天山地块,天山地块南北向的缩短速率为(10.1±0.9)mm/a。这两个地块目前仍处于隆升状态,从面应变看,面膨胀在中国大陆占优势,东部基本都是膨胀区,在西部面压缩与面膨胀从南向北相间分布。中国大陆的大多数东西向或近东西向断裂两侧的相对运动都是左旋或类似左旋走滑型的,大多数南北向断裂两侧的相对运动都是右旋或类似右旋走滑型的。GPS测定的阿尔金断裂中部的左旋走滑速为(4.8±1.3)mm/a,鲜水河断裂的左旋走滑速为(9.8±2.2)mm/a。地块边界断裂带的运动为地块运动创造了条件,地块及其边界的运动是协调一致的统一的,各个地块的活动程度是不相同的,统计检验结果表明,大多数地块之间的相对运动是显著的与非常显著的,这证明活动地块是客观存在的,喜马拉雅、拉萨、天山、羌塘和滇西南是活动最强烈的地块,中蒙、中朝西、阿拉善和华南是较稳定的地块,印度、太平洋、菲律宾板块与欧亚板块的互相作用力是中国大陆地块运动的主要驱动力。青藏高原地壳物质在印度板块NNE向的强烈推挤下,向NNE和NE方向运动,由于受到北部、东北部和东部地块的阻挡,经高原的东南部向印度洋方向运移,     

利用GPS资料研究郯庐带现今运动及变形状态

本文利用华北地区1999—2001年, 2001—2004年, 2004—2007年及2007—2009年四期GPS速度场资料, 基于块体的整体旋转与均匀应变模型, 分析了郯庐断裂带中南段的运动及变形特征。 结果显示: 潍坊—郯城段主要为右旋走滑的变形特征, 而郯城—庐江段则为左旋滑动的变形特征, 两段垂直断层方向上的变形表现为“张压交替”的特征。 基于刚体运动模型, 计算了扣除环渤海湾区域整体刚性运动的华北地区GPS速度场, 并分析了环渤海湾区域块体的变形状态, 结果显示环渤海区虽然各期的变形特征不同, 特别是郯庐带附近, 各期的运动特征差异较大, 但基本可反映燕山—渤海地震带是运动特征差异的分界线, 且每期郯庐带各站点的一致性运动明显。     

THE RELATIONSHIP BETWEEN REGIONAL SEISMIC ENERGY RELEASE AND RELATIVE MOTION BETWEEN BLOCKS ON BOTH SIDES OF YISHU FAULT ZONE

Following the 11 March 2011 Japan M_W9.0 earthquake, frequent moderate and small events occurred on the Yishu fault zone and its either side. Using continuous GPS data and a sliding block model, this work studies the relationship between the energy release of these shocks and the block relative motion of either side of the Yishu fault zone. The results show that(1)the equivalent magnitude M from released energy and the two blocks' relative motion are well correlated when earthquakes are selected in a retrieval circle(whose center is the midpoint of the Yishu fault zone)with a radius of 250~500km and using a sliding time window of 3~10 months. The best correlation coefficient between M and the two blocks' relative motion is 0.74 and the T test shows a significant linear correlation between them.(2)Spatial distribution of the correlation coefficients shows that the relative motion of the blocks on both sides affects the energy release in the area from the north part of Yishu fault zone to the Jiaodong Peninsula area and southwest Shandong-Henan border area obviously.(3)Since June 2014, the relative motion of the two blocks on both sides of the Yishu fault zone presents a wave of change, which may be an expression of the accumulation of seismic strain energy in the Yishu fault zone and its two sides. The linear relationship between the equivalent magnitude M from released energy and two blocks' relative motion V can be fitted by linear equation M=0.51*V+3.9, showing that strain energy accumulation could be released by the moderate and small earthquakes in a timely manner, which may favorable to delay the seismic risk in the study area. It also shows, on the other hand, that earthquake energy was not released so completely in the study area since the end of 2015 to 2016, which is likely associated with the Changdao earthquake swarm in 2017.     

Rotation of vertically oriented objects during earthquakes

Vertically oriented objects, such as tombstones, monuments, columns, and stone lanterns, are often observed to shift and rotate during earthquake ground motion. Such observations are usually limited to the mesoseismal zone. Whether near-field rotational ground motion components are necessary in addition to pure translational movements to explain the observed rotations is an open question. We summarize rotation data from seven earthquakes between 1925 and 2009 and perform analog and numeric rotation testing with vertically oriented objects. The free-rocking motion of a marble block on a sliding table is disturbed by a pulse in the direction orthogonal to the rocking motion. When the impulse is sufficiently strong and occurs at the ‘right’ moment, it induces significant rotation of the block. Numeric experiments of a free-rocking block show that the initiation of vertical block rotation by a cycloidal acceleration pulse applied orthogonal to the rocking axis depends on the amplitude of the pulse and its phase relation to the rocking cycle. Rotation occurs when the pulse acceleration exceeds the threshold necessary to provoke rocking of a resting block, and the rocking block approaches its equilibrium position. Experiments with blocks subjected to full 3D strong motion signals measured during the 2009 L’Aquila earthquake confirm the observations from the tests with analytic ground motions. Significant differences in the rotational behavior of a monolithic block and two stacked blocks exist.     

青藏块体东北缘近期水平运动与变形

利用青藏块体东北缘地区13、1年GPS观测资料，给出了本区地壳水平运动速度场及视应变场分布图，提出了由位移观测值直接求解块体旋转和变形参数的方法，初步研究了本区构造块体运动与变形特征.结果表明：①本区存在整体性向东-东南方的运动(速率约mm/a)；②南部的甘肃-青海块体的运动较快，而北部的阿拉善块体的运动较慢，二者运动速率相差近6mm/a，祁连-海原断裂带左旋走滑运动显著.③自西向东存在北北东-北东东向压性运动；④阿拉善块体、甘肃-青海块体内部存在北西西向张性变形，阿拉善块体的整体张性变形更显著，鄂尔多斯块体西侧的块体交接地带为压性运动.     

On the seismic performances of rigid block‐like structures coupled with an oscillating mass working as a TMD

In this paper, the effects of a mass damper on the rocking motion of a non‐symmetric rigid block‐like structure, subject to different seismic excitation, are investigated. The damper is modelled as a single degree of freedom oscillating mass, running at the top of the block and connected to it by a linear visco‐elastic device. The equations of rocking motion, the uplift and the impact conditions are derived. A nondimensionalisation of the governing equations is performed with the aim to obtain an extensive parametric analysis. The results are achieved by numerical integration of these equations. The slenderness and the base of the rigid block, and the eccentricity of the centre of mass are taken as variable parameters in the analyses. The main objective of the study is to check the performance of the damper versus the spectral characteristics of the seismic input. Three earthquake registrations with different frequency contents are used in the analyses. The results show that the presence of the mass damper leads to different levels of improvement of the response of the system, depending on the spectral characteristics of the seismic input. Curves providing the overturning slenderness of blocks of specific sizes versus the characteristics of the TMD are obtained. Copyright © 2017 John Wiley & Sons, Ltd.     

青藏高原东北隅似三联点构造特征

在大陆内部,对于一定尺度的构造块体,似三联点构造是一种普通的构造形式。由于大陆内部块体旋转的普遍性,似三联点构造可形成顺旋型和逆旋型两种构造型式。     

PROGRESS IN APPLICATION OF GNSS TO DIVISION OF ACTIVE TECTONIC BLOCKS IN CONTINENTAL CHINA

Chinese scientists proposed that large earthquakes that occurred in mainland China are controlled by the movement and deformation of active tectonic blocks. This scientific hypothesis explains zoned phenomenon of seismicity in space. The active tectonic blocks are intense active terranes formed in late Cenozoic and late Quaternary, and the tectonic activity of block boundaries is the intensest. Global Navigation Satellite System(GNSS)has advantages of high spatio-temporal resolution, broad coverage, and high accuracy, and is utilized to monitor contemporary crustal deformation. High accuracy and resolution of GNSS velocity field within mainland China and vicinities provided by previous studies clearly demonstrate that different active tectonic blocks behave as different patterns of movement and deformation, and block interaction boundaries have intense tectonic deformation. The paper firstly introduces the GPS networks operated by the Crustal Movement Observation Network of China(CMONOC)since 1999, and GNSS data processing methods, including GAMIT, BERNESE and GIPSY/OASIS, and discusses the advantages of using South China block as a regional reference frame for GNSS velocity field, then proposes three strategies of block division, F-test, quasi-accurate detection(QUAD), and clustering analysis. Furthermore, we introduce rigid and non-rigid block motions. Rigid block motion can be denoted by translation and rotation, while non-rigid block motion can be described by rigid motion and internal strain deformation. Internal strain deformation can be divided into uniform and linear strains. We also review the usage of F-test to distinguish whether the block acts as rigid deformation or not. In addition, combining with recent GNSS velocity results, we elaborate the characteristics of present movement of rigid block, such as the South China, Tarim, Ordos, Alashan, and Northeast China, and that of non-rigid block, such as the Tibetan plateau, Tian Shan, and North China plain. Especially, the Tibetan plateau and Tian Shan seem to deform continuously with significant internal deformation. In order to enrich and perfect the active tectonic block hypothesis, we should carefully design dense GNSS networks in inner blocks and block boundaries, optimize utilizing other space geodesy technologies such as InSAR, and strengthen combining study of geodesy, seismogeology and geophysics. Through systematic summary, this paper is very useful to employing GNSS to investigate characteristics of block movement and dynamics of large earthquakes happening in block interaction boundaries.     

地壳介质的各向异性弹性本构关系讨论

论述了地壳介质中存在的几类各向异性弹性本构关系，推导了应变能变换公式、EDA薄层各向异性弹性常数及非对称面上EDA构成的单斜对称弹性常数的计算公式，讨论了几种各向异性对称系统的特征方程，并给出了弹性波的角散表达式，最后选取合理的参数对波的角散现象进行了数值模拟。