GNSS空间大地测量技术在中国大陆活动地块划分中的应用和研究进展

引用本文

郝明, 王庆良. GNSS空间大地测量技术在中国大陆活动地块划分中的应用和研究进展[J]. 地震地质, 2020,42(2): 283-296.
HAO Ming, WANG Qing-liang. PROGRESS IN APPLICATION OF GNSS TO DIVISION OF ACTIVE TECTONIC BLOCKS IN CONTINENTAL CHINA[J]. SEISMOLOGY AND GEOLOGY, 2020,42(2): 283-296. 复制到剪切板

DOI:10.3969/j.issn.0253-4967.2020.02.003
Permissions

《地震地质》编辑部所有

GNSS空间大地测量技术在中国大陆活动地块划分中的应用和研究进展

郝明, 王庆良^*

中国地震局第二监测中心, 西安 710054

*通讯作者: 王庆良, 男, 1963年生, 博士, 研究员, 主要从事地震动力学与地球探测技术、地质资源与灾害地质等方面的研究, 电话: 029-85506503, E-mail:wangql63@163.com。

〔作者简介〕郝明, 男, 1982年生, 2012年于中国地震局地质研究所获固体地球物理专业博士学位, 副研究员, 主要从事地震大地测量数据处理与应用等方面的研究工作, 电话: 029-85506514, E-mail: ha_mg@163.com。

收稿日期: 2020-02-21

基金项目: 国家重点研发计划项目(2017YFC1500102)、第二次青藏高原综合科学考察(2019QZKK0901)和国家自然科学基金(41874017, 41874117)共同资助

摘要

中国科学家提出了“中国大陆的强震受控于活动地块运动与变形”的科学假说, 很好地解释了地震活动在空间上的分区、分带现象。活动地块为形成于晚新生代、晚第四纪的强烈活动的地质构造, 地块边界带的构造活动性最为强烈。以GNSS为代表的空间大地测量技术以其时空分辨率高、覆盖范围广、观测精度高等特点, 被应用于现今的地壳形变监测。基于GNSS建立的中国大陆及周边地区高空间分辨率的地壳运动速度场图像, 清晰地揭示了中国大陆内部不同活动地块的不同运动与变形方式, 以及地块相互作用的构造边界带的强烈变形。文中回顾了用于监测中国大陆地壳运动的GNSS观测数据来源及处理方法, 提出了中国大陆GNSS运动速度场以稳定的华南地块为区域参考框架的优势, 介绍了3种划分活动地块的方法以及检验地块是否具有刚性运动的判别方法。结合GNSS速度场的最新研究成果, 阐述了具有刚性运动的华南、塔里木、鄂尔多斯、阿拉善和东北地块现今的运动特征, 以及青藏高原、天山和华北平原地块的非刚性变形。在利用空间大地测量技术丰富和完善活动地块假说理论方面, 不仅需要地块内部和边界带分布密集的GNSS观测网络, 同时也要充分利用InSAR等其它空间大地测量技术, 并需加强与地震地质、地球物理等多学科的交叉融合。文中的系统总结, 对利用GNSS技术研究活动地块的运动变形特征和活动地块边界带强震孕育的动力学机理具有重要的参考价值。

关键词: 活动地块; GNSS; 地块划分; 刚性运动; 连续变形

中图分类号:P315.2 文献标志码:A 文章编号:0253-4967(2020)02-0283-14

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

HAO Ming, WANG Qing-liang

The Second Monitoring and Application Center, China Earthquake Administration, Xi'an 710054, China

Abstract

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.

Keyword: active tectonic block; GNSS; block division; rigid movement; continuous deformation

文章图片

0 引言

20世纪60年代建立并发展起来的板块构造学说合理解释了全球80%的地震发生在板块边界的现象, 但由于板内地震和板块边界地震在动力环境、构造成因、变形机制和破坏类型等方面存在很多显著差别, 板块构造学说在解释大陆板内地震时遇到了诸多难以解释的问题(张国民等, 1999)。随着板块构造学说被引入中国, 老一辈中国科学家针对中国大陆的板块构造体系和特征提出了断块构造(张文佑, 1984)、活动亚板块(马杏垣, 1987)、构造块体(丁国瑜, 1991)以及活动断块(邓起东等, 1994)等理论或划分方案。

基于国家重点基础研究发展计划(973)项目“ 大陆强震机理与预测” 研究, 在前人得到的大量研究基础和成果下, 中国科学家提出了“ 中国大陆的强震受控于活动地块运动与变形” 的科学假说, 用于解释中国大陆强震的空间分布(张国民等, 1999, 2000; 张培震, 1999; 张培震等, 2003)。张培震等(2003)明确了活动地块的定义, 即活动地块是被形成于晚新生代、晚第四纪至现今强烈活动的构造带所分割和围限的、具有相对统一运动方式的地质单元; 活动地块边界构造活动强烈, 绝大多数强震都发生在边界的活动构造带上。从以上定义可以看出, 该假说在时间上强调是晚第四纪至现今时段的活动状态, 地块边界活动构造带是活动地块间相对运动和相互作用的形变带, 其活动强度和频度都远大于地块内部的构造活动。

近30a来, 以GNSS为代表的空间大地测量技术因其时空分辨率高、覆盖范围广、观测精度高等特点, 被应用于现今地壳形变监测, 为促进中国大陆活动地块的理论研究发挥了重要作用。 Wang 等(2001)首次应用GNSS数据研究了中国大陆现今的地壳运动和变形, 其结果清晰地反映了以活动地块为单元的分块运动特征。随着国家重大科学工程“ 中国地壳运动观测网络” (简称网络工程)和“ 中国大陆构造环境监测网络” (简称陆态网络)项目的开展, 利用丰富的GNSS观测数据精细地刻画出中国大陆现今地壳的水平运动图像(Niu et al., 2005; 李强等, 2012; Liang et al., 2013; Wang et al., 2017; Zheng et al., 2017; Hao et al., 2019; Rui et al., 2019; Wang et al., 2020), 揭示了中国大陆现今的总体构造变形具有刚性和非刚性的运动学特征(张培震等, 2002, 2004, 2005; 李延兴等, 2003; 王敏等, 2003; Zhang et al., 2004; Shen et al., 2005; Gan et al., 2007), 为大陆变形的动力学机理研究提供了重要支撑。以GNSS运动速度场为约束, 采用块体模型可估计出活动地块边界带现今的滑动速率和闭锁强度(Meade, 2007; Thatcher, 2007; Wang et al., 2011; 李煜航等, 2014; Jiang et al., 2015; 郝明等, 2017), 为预测中长期强震的危险性提供了基础数据。

本文首先梳理了用于监测中国大陆地壳运动的GNSS观测网络及数据处理方法, 其次介绍了划分活动地块的方法以及检验地块是刚性变形还是非刚性连续变形的判别方法, 然后结合GNSS速度场的最新研究成果阐述了中国大陆活动地块具有的刚性运动与连续变形特征, 最后对利用空间大地测量技术丰富和完善活动地块假说等工作进行了展望。

1 中国大陆GNSS高精度监测网络体系与现今地壳的运动图像

目前用于监测中国大陆地壳变形的GNSS观测数据主要来自“ 九五” 和“ 十一五” 建设的网络工程(牛之俊等, 2002)和陆态网络(甘卫军等, 2012)项目。网络工程由32个连续观测站、 56个基本站和1 000个区域站组成, 其中基本站于1998— 2009年每年复测1次, 区域站于1999年、 2001年、 2004年和2007年观测。陆态网络项目则是将网络工程项目的GNSS连续站增加至260个, 区域站增加至2 000个, 并于2009— 2017年间每2a开展1次复测。依托中国地震局地震行业科研重点专项“ 中国综合地球物理场” 一、二、三期项目, 由中国地震局第二监测中心牵头并联合多家单位于2012年、 2014年和2016年对位于青藏高原东缘、鄂尔多斯地块周缘以及华北和东北地区的陆态网络区域站进行了加密复测。在Wang 等(2020)提供的中国大陆及周边最新的GNSS运动速度场结果中还使用了一些区域网络GNSS数据, 这些区域网络主要布设在活动构造带, 如华北、汾渭断陷带和川滇地区。图 1 为网络工程和陆态网络项目的流动和连续GNSS测站的分布情况。

	Figure Option View Download New Window
	图 1 网络工程和陆态网络项目GNSS测站分布图(据Wang et al., 2020修改)Fig. 1 Locations of GNSS stations operated by the Crustal Movement Observation Network of China(CMONOC)(after Wang et al., 2020).

GNSS观测数据处理主要采用GAMIT/GLOBK、 BERNESE或GIPSY/OASIS软件完成。前两者采用双差模式获取以24h为1个时段的单日松弛解, 后者则采用非差观测得到单日松弛解, 然后将单日区域网解与全球网解合并为1个单日松弛解, 并通过赫尔默特七参数转换方法将区域网解固定在ITRF全球参考框架下, 以获取各GNSS测站的位移时间序列。中国大陆内部及周边在GNSS观测时段内发生了数次大地震, 而这势必影响了GNSS捕捉到的地壳震间期的变形。 Wang (2020)针对上述问题详细地阐述了消除大地震的同震和震后变形的方法。

为揭示中国大陆的活动构造变形特征, 通常将ITRF框架下的GNSS速度场转换至相对于欧亚板块的速度场。 Wang 等(2020)给出的高时空分辨率的GNSS运动速度场图像整体反映出中国大陆呈顺时针旋转的运动状态, 青藏高原内部的一部分地壳物质沿NE向流动迁移, 一部分沿SE向绕喜马拉雅东构造结发生顺时针旋转, 似乎还有一小部分近E向挤入华南地块, 华北和华南地块的差异性运动并不十分明显。

Hao等(2019)选择中国大陆内部稳定的华南地块作为区域参考基准, 所提供的东亚地区的GNSS运动速度场图像清晰地展示出中国大陆受到印度板块的推挤, 并在太平洋和菲律宾板块的俯冲作用下其内部产生了差异性构造变形特征。因此, 本文将Wang 等(2020)提供的相对于欧亚板块的GNSS速度场旋转为相对于华南地块的GNSS速度场(图 2)。

	Figure Option View Download New Window
	图 2 中国大陆及周边GNSS水平运动速度场(1999— 2017年, 相对于华南地块)Fig. 2 GNSS velocity filed within continental China and vicinities with respect to the South China block(1999— 2017).

2 基于GNSS技术划分活动地块的方法

活动地块假说认为地块内部具有相对统一的运动方式, 但不同地块之间存在相对运动。因此, 利用GNSS观测技术建立的现今高精度、高时空分辨率的构造变形运动图像, 可为研究划分活动地块提供定量约束(张培震等, 2005)。需要注意的是, 在利用GNSS划分地块之前, 需要参考已有的活动地块模型作为初始模型, 然后再结合GNSS速度场进行研究划分。对存在明显形变的梯度带, 还需要结合活动断裂和地震活动性资料, 以判断是否划分为新的地块边界带。

如果活动地块内部不存在相对变形, 可将其视为刚性块体。球面上刚性块体的运动可以用绕欧拉极的旋转(李延兴等, 2003)来描述:

$[\begin{array}{l} V_{e} \\ V_{n} \end{array}] = r [\begin{array}{l} - \sin φ \cos λ & - \sin φ \sin λ & \cos φ \\ \sin λ & - \cos λ & 0 \end{array}] [\begin{array}{l} ω_{x} \\ ω_{y} \\ ω_{z} \end{array}]$ (1)

式中, V_e和V_n为地块上GNSS测站大地坐标(λ , φ )的E向和N向速度, r为地球的平均半径, ω _x、 ω _y、 ω _z为地块的欧拉旋转矢量。对于地块内2个以上的GNSS测站速度, 可根据式(1)利用加权最小二乘方法估计出其欧拉矢量。

如何识别具有同一运动趋势的相对稳定点是划分活动地块的关键技术。下面简要介绍3种方法。

2.1 F检验法

(1)检验单个测站的稳定性。对一个给定的活动地块, 根据式(1)使用地块内所有的GNSS水平运动速度拟合地块的欧拉矢量, 然后计算各GNSS测站的观测值与模型值的残差, 如果残差的绝对值> 2mm/a, 则删除该测站。重复此步骤, 直至所有GNSS测站的拟合残差绝对值均< 2mm/a。

(2)检验每个测站与其它测站的兼容性(Shen et al., 2005)。对执行完上一步后地块内剩余的GNSS测站, 再利用式(1)计算欧拉旋转的拟合残差平方和 $χ_{n}^{2}$ , 然后每次去除1个不同的测站, 计算去除该测站后欧拉旋转的拟合残差平方和 $χ_{n - 1}^{2}$ , 并求

$F = F (χ_{n}^{2}, 2 n - 3; χ_{n - 1}^{2}, 2 n - 5) = \frac{χ_{n}^{2}}{2 n - 3} / \frac{χ_{n - 1}^{2}}{2 n - 5}$ (2)

式中, n为所选用的GNSS测站的个数。若取置信度为95%, 即如该测站与其它测站相互独立的显著性> 95%, 则去除该测站, 否则认为该GNSS测站为属于地块的稳定点。

(3)检验相邻2个地块的独立性(王敏等, 2003; 王敏, 2008)。分别计算2个地块欧拉旋转的拟合残差平方和 $χ_{1}^{2}$ 和 $χ_{2}^{2}$ , 以及假设2个地块合为1个地块的欧拉旋转拟合残差平方和 $χ_{12}^{2}$ , 则

$F = F (χ_{12}^{2}, 2 n_{1} + 2 n_{2} - 3; χ_{1}^{2} + χ_{2}^{2}, 2 n_{1} + 2 n_{2} - 6) = \frac{χ_{12}^{2}}{2 n_{1} + 2 n_{2} - 3} / \frac{χ_{1}^{2} + χ_{2}^{2}}{2 n_{1} + 2 n_{2} - 6}$ (3)

式中, n₁和n₂分别为块体一和块体二内的GNSS测站数。仍取置信度为95%, 即如果2个相邻地块相互独立的显著性 < 95% , 则将2个地块合二为一。通过循环的方式对各相邻地块的独立性进行F检验。

2.2 拟准检定法

拟准检定法(欧吉坤, 1999)是在“ 拟稳平差” 理论的基础上发展起来的, 其在真误差与观测值之间建立解析关系用于探测观测值中的粗差。将该方法应用于地块划分, 即是将地块内的相对稳定点作为拟准观测, 而将不符合运动趋势的点看作粗差并予以剔除(柴艳菊等, 2002)。该方法的关键是确定拟准观测, 可分为初选和复选2个阶段:

(1)初选。将式(1)线性化为观测方程

$V = A \hat{X} - L$ (4)

式中, A为设计矩阵, $\hat{X}$ 为未知参数即欧拉旋转矢量的估值, L为观测值, 即GNSS测站的速度(V_e和V_n), V为观测值的残差。观测值的权矩阵记为P。

未知参数的估值 $\hat{X}$ =(A^TPA)^-1A^TPL, 计算平差因子矩阵J=A(A^TPA)^-¹A^TP以及其正交补投影矩阵R=I-J, I为单位矩阵。

计算所有GNSS测站2个分量的残差V_i([V_e_iV_n_i]^T)以及V_i的中位值σ ₀=med|V_i|, 将其作为单位权中误差的估值 ${\hat{σ}}_{0}$ , 有较强的抗差性。

计算初选指标值u_i=V_i/( ${\hat{σ}}_{0} \sqrt[]{R_{ii}}$ ), R_ii 为投影阵的第i个对角元。如果GNSS测站的残差V_i≤ u_i× f₁, 则将该测站初选为相对稳定点, 否则视为不稳定点。其中, f₁为判别准则因子。

(2)复选。基于选出的相对稳定点重新计算未知参数估值 ${\hat{X}}_{0}$ =( ${A^{T}}_{0}$ P₀A₀)^-1 ${A^{T}}_{0}$ P₀L₀。其中, A₀、 P₀、 L₀为A、 P、 L中与选出的相对稳定点有关的部分。真误差的估值 ${\hat{Δ}}_{0}$ =-(R^TR+P₀G^TGP₀)^-1R^TRL₀, G^T=[0 A₀]^T。计算指标W_i=| ${\hat{Δ}}_{i}$ |/C, C=( ${\hat{Δ}}^{T}_{0}$ P₀ ${\hat{Δ}}_{0}$ /(r-1))^1/2, r为稳定点的个数。

若W_i≤ f₂, 则将该测站归入相对稳定点, 否则作为不稳定点。 f₂为判别准则因子。重复该步骤, 直到没有再找到新的稳定点时结束。

由拟准检定法的计算流程可以看出, f₁和f₂ 2个判别因子的取值大小决定此筛选方法标准的宽或严。当一个块体上的观测点不多时, 较小的判别因子可能会造成没有足够数量的点归入到稳定点的情况; 当观测点较多时, 选取较大的判别因子又可能会使估计出的地块运动和变形参数的误差较大。因此, 需根据数据的质量情况进行多次试算和选取。

2.3 聚类分析法

聚类分析是指将数据对象的集合分组为由类似的对象组成的多个类的分析过程。聚类生成的组称为簇, 而簇是数据对象的集合。簇内部的任意2个对象之间具有较高的相似度, 属于不同簇的2个对象间具有较高的相异度。 Simpson等(2012)将层次凝聚聚类(HAC)方法应用于加州旧金山湾地区的GNSS速度场, 试图找出该地区是否具有类似地块变形的特征。 Savage等(2013a, b, 2015)将k-中心点聚类分析方法应用于圣安地列斯断裂和美国西北部地区的地块划分。

利用聚类分析方法划分地块的流程为: 首先, 需预先给定簇的个数k, 即地块数目。从n个GNSS测站速度值中随机挑选出k个(分别挑选E分量和N分量), 作为初始簇的中心; 根据剩余测站的速度值与各簇中心距离的平方, 将具有最小距离平方的速度值赋给最近的簇; 然后重新计算每个簇的中心点。从剩余测站的速度值中随机挑选出k个, 计算用这个k个值分别代替原中心点形成新聚类的总代价S, 如果S< 0, 则用从剩余速度值中挑选的非中心点代替中心点, 形成新的k个中心点集合。重复上述过程, 直到簇的中心点不发生变化, 认为聚类结束。

使用聚类分析方法划分地块时, 需通过先验信息确定地块的个数。 k-中心点聚类分析法抗噪声(不合群速度值)的能力强, 但运行较为耗时。

3 基于GNSS技术检验活动地块的刚性与非刚性变形

地块在周围板块或地块的作用下不仅会发生刚性运动(平动和旋转), 其内部可能会产生均匀或非均匀的线性应变变形(李延兴等, 2003, 2004)。

如果地块的边界断裂不发生闭锁, 则地块的刚性运动和内部均匀应变变形可用式(5)表示(王敏, 2008; Shen et al., 2015):

$[\begin{array}{l} V_{x} \\ V_{y} \end{array}] = [\begin{array}{l} 1 & 0 & Δ x & Δ y & 0 & Δ y \\ 0 & 1 & 0 & Δ x & Δ y & - Δ x \end{array}] [\begin{array}{l} u_{x} \\ u_{y} \\ τ_{xx} \\ τ_{xy} \\ τ_{yy} \\ ω \end{array}] + Δ$ (5)

式中, V_x、 V_y为GNSS测站的E向和N向速度, Δ x、 Δ y为GNSS测站与参考点之间的相对位置, u_x、 u_y为地块的E向和N向平移速度, τ _xx、 τ _xy、 τ _yy 为应变张量, ω 为旋转率。如果应变张量为零, 即为刚性块体模型, 则式(5)与式(1)等价。

如果地块内部的应变不是均匀的而是随位置线性变化的, 则地块的刚性运动和内部线性应变变形可以用式(6)表示(李延兴等, 2004; Qu et al., 2014):

$[\begin{array}{l} V_{x} \\ V_{y} \end{array}]$ =r $[\begin{array}{l} - \sin φ \cos λ & - \sin φ \sin λ & \cos φ \\ \sin λ & - \cos λ & 0 \end{array}] [\begin{array}{l} ω_{x} \\ ω_{y} \\ ω_{z} \end{array}]$ + $[\begin{array}{l} A_{0} & B_{0} \\ B_{0} & C_{0} \end{array}] [\begin{array}{l} x \\ y \end{array}]$ + $\frac{1}{2} [\begin{array}{l} A_{1} & B_{2} \\ B_{1} & C_{2} \end{array}] [\begin{array}{l} x^{2} \\ y_{2} \end{array}]$ + $[\begin{array}{l} A_{2} + B_{1} \\ B_{2} + C_{1} \end{array}]$ +Δ (6)

式中, A₀、 A₁、 A₂、 B₀、 B₁、 B₂、 C₀、 C₁和C₂分别是待求的应变参数。

此外, 同样可以用F检验方法来判定某个地块为刚性平动、旋转变形或非刚性的连续变形(王敏, 2008)。对于一活动地块, 首先假设其为刚性运动, 用内部GNSS测站的水平速度计算欧拉运动的3个参数(u_x、 u_y和ω ), 并计算模型的拟合残差平方和 $χ_{1}^{2}$ 。然后, 假设该地块内部存在应变变形, 采用式(5)估计欧拉运动的3个参量(u_x、 u_y和ω )以及应变率张量的3个参数(τ _xx、 τ _xy和τ _yy), 并计算拟合后残差平方和 $χ_{2}^{2}$ 。最后通过F检验判断地块内部应变量的显著性:

$F = F (χ_{1}^{2}, 2 n - 3; χ_{2}^{2}, 2 n - 6) = \frac{χ_{1}^{2}}{2 n - 3} / \frac{χ_{2}^{2}}{2 n - 6}$ (7)

当F检验的显著性> 95%才允许地块内部存在均匀应变变形。该方法同样也适用于地块内部的线性应变模型检验。

4 基于GNSS技术获得的中国大陆主要活动地块的运动与变形

张培震等(2003)将中国大陆及其周边地区划分为6个I级地块区和22个Ⅱ 级活动地块(图 3)。随后, 很多学者以此为基础模型, 采用时空分辨率不断提高的GNSS速度场结果, 愈来愈细致地揭示出中国大陆及周边活动地块的运动学特征(Shen et al., 2005; Meade, 2007; Thatcher, 2007; 王敏, 2008; Loveless et al., 2011; Wang et al., 2011; Zhang et al., 2013; 李煜航等, 2014; Wang et al., 2017; Wang et al., 2020)。

	Figure Option View Download New Window
	图 3 中国大陆主要的活动地块与强震分布(张培震等, 2013)Fig. 3 Active tectonic blocks and distribution of large earthquakes in continental China(after ZHANG Pei-zhen et al., 2013).

总体而言, 中国大陆的活动地块有的刚性较好, 内部不发生构造变形和大地震, 有的地块则刚性很差, 内部发生明显的构造变形并伴随强震发生(张培震等, 2013)。 Wang 等(2020)采用刚性和连续应变模型(式(5)), 利用中国大陆的GNSS速度场拟合了东北、华北、华南、鄂尔多斯、阿拉善、塔里木和准噶尔地区的运动变形, 结果表明这几个地区具有刚性地块运动特征, 但地块内部的变形(10^-9/a量级)也不可忽视(图 4), 而青藏高原和天山地区则具有明显的连续变形特征。

Figure Option
View Download New Window

图 4 中国大陆典型活动地块运动模型(据Wang et al., 2020修改)
红色箭头、蓝色扇形和黑色箭头对分别代表地块相对欧亚板块的平动矢量、地块旋转率(逆时针旋转为负)和地块内部的主应变率, 蓝色粗箭头示意印度板块和青藏高原的运动方向Fig. 4 Block-like movement within continental China(modified from Wang et al., 2020).

4.1 刚性活动地块

如式(1)所示, 块体的刚性运动可用球面上刚性块体绕欧拉极的旋转来描述。块体内部各点的运动速度随着旋转半径长度的增加呈有规律的增加; 在以欧拉旋转极为圆心的同一大圆弧上, 各观测点的运动速率相等; 相邻刚性块体之间的相对运动和变形主要集中在块体边界带上, 块体内部不发生大规模变形(张培震等, 2005)。

夹持在构造活动十分强烈的天山和青藏高原之间的塔里木地块的构造活动性相对较弱, 有历史记载以来地块内部没有发生过强震, 也不发育晚第四纪活动断裂(Deng et al., 2003)。 GNSS拟合结果表明塔里木地块的顺时针旋转量最大((-0.630± 0.017)° /Ma), 并以13.8mm/a的速率向N运动, 属典型的刚性块体(张培震等, 2005; Wang et al., 2020)。

青藏高原北侧的阿拉善地块表现为顺时针旋转, 但旋转量仅有(-0.103± 0.023)° /Ma。地块以4mm/a的速率向E运动, 地块内部最大和最小主应变率相近, 约2× 10^-9/a(Wang et al., 2020)。可见, 尽管阿拉善地块内部存在微弱的应变变形, 但刚性块体的运动特征仍十分显著。

鄂尔多斯地块除西南缘受到青藏高原东北部的强烈挤压作用外, 其它边缘被断陷盆地所围限且以剪切变形为特征, 控制了有历史记载以来的7次7级以上强震的发生, 表现出很强的构造活动性。但鄂尔多斯地块内部构造的活动性微弱, 也不发育大规模的活动断层(张培震等, 2005)。 GNSS资料表明鄂尔多斯地块具有逆时针旋转的特点, 旋转量为(0.172± 0.017)° /Ma, 以5.7mm/a的速率向SEE运动, 地块内部的主应变率约1× 10^-9/a(Wang et al., 2020), 说明该地块属于刚性块体。

华南和东南沿海地区在新构造运动上属于比较稳定的地块, 内部不发育明显的活动断裂和褶皱, 地震活动性与华北和西部相比相对较弱(张培震等, 2005)。华南地块内的GNSS速度场表明该地块整体以7.5mm/a的速率向SEE运动, 逆时针旋转量非常小, 地块内部的主应变率< 1× 10^-9/a(Wang et al., 2020)。

位于阿穆尔板块上的东北地块强震活动相对微弱(深震除外), 构造活动相对稳定。地块内部近200个GNSS台站的运动速度表明该地块整体的SEE向运动速率仅约2mm/a, 逆时针旋转量(0.046± 0.006)° /Ma也很小, 地块内部变形微弱(Wang et al., 2020), 符合刚性块体的运动特征。

4.2 非刚性活动地块

青藏高原和天山地区不仅活动断裂十分发育, 且6级以上地震广泛分布, 7级以上的强震沿主要断裂带频繁发生。 Zhang 等(2004)和Gan等(2007)利用青藏高原及其周边的GNSS观测资料, 发现印度和欧亚板块相对运动速率的70%~94%被青藏高原内部的构造变形所调整吸收。高原周边以收缩应变为主, 表现为地壳缩短增厚; 内部以扩张和剪切应变为主。 Zhang 等(2013)利用GNSS资料建立了青藏高原东缘的可变形地块模型, 结果表明青藏高原地块内部的变形明显大于其周边的华南地块和阿拉善地块。地块的内部变形和旋转量非常显著, 青藏高原内部各地块均是可变形的非刚性活动地块(Zhang et al., 2008)。

天山南、北两侧分别被构造稳定的准噶尔和塔里木地块所夹持, 构造变形主要发生在天山山前和山体内部(Zhang, 2004)。强震主要发生在天山两侧边界的前陆盆地的活动褶皱和逆冲断裂带上。横跨天山的SN向GNSS剖面揭示地壳缩短是均匀的, 而不是集中在两侧的前陆冲断带或某一断层附近(张培震等, 2002)。此外, 天山的地壳缩短由西向东逐渐衰减, 这种逐渐减小的连续变形可以用刚性的塔里木块体相对“ 软弱” 的天山发生顺时针旋转所导致的地壳缩短变化来解释(Zhang, 2004)。

华北平原活动地块从长期、平均构造变形的角度来看不具有刚性块体的性质, 其内部发育一系列活动断裂, 并伴有大量历史强震发生(Deng et al., 2003)。现今的GNSS资料表明华北地块整体以5.5mm/a的速率向SEE运动, 逆时针旋转量(0.086± 0.006° )/Ma和约2.7× 10^-9/a的内部主应变率(Wang et al., 2020)均比其南、北两侧的华南和东北地块高。

5 总结与展望

中国大陆晚新生代的构造变形主要表现为活动地块的运动学特征, 不同活动地块的运动方式也不相同, 地块之间相互作用的地块边界带构造变形最为强烈, 也是强震高发区(张培震等, 2003)。近30a发展起来的GNSS观测技术, 能够在短时间内获取大范围、高精度、高空间分辨率的构造变形运动学图像, 为活动地块划分的应用和研究提供了重要基础数据。

(1)利用现代空间大地测量技术GNSS可监测现今地壳变形、划分活动地块、识别活动地块内部的变形特征。在利用GNSS划分活动地块时, 需结合活动断裂探测、地震台阵观测、大地电磁测深等观测成果综合研究。已有的研究结果表明, 中国大陆由2类活动地块所组成(张培震等, 2005, 2013; Zhang et al., 2008): 第一类是刚性变形, 即地块的内部结构完整, 不发生内部构造变形, 地震活动性也比较弱, 华南、塔里木、鄂尔多斯、阿拉善和东北地块属于第一类; 第二类是非刚性变形, 即地块内部发生相对变形, 地块的完整性较差, 内部各测点之间发生相对运动, 地震活动发生在整个块体内部, 变形和运动不能用刚体运动来描述, 青藏高原、天山和华北平原属于该类地块。

(2)如前文所述, 利用GNSS建立活动地块的非刚性连续变形模型时, 假设地块边界断裂带自由滑动, 没有考虑断层闭锁。如果考虑断层闭锁, 块体模型可用块体欧拉旋转、块体内部均匀应变以及块体边界由于断层闭锁产生的地表弹性变形3个分量描述(McCaffrey, 2005)。以GNSS运动速率作为约束, 反演活动地块边界断裂的滑动亏损速率和断层面上闭锁强度的时空演化过程, 无疑可为地震危险性预测提供重要支撑。

(3)为刻画中国大陆构造变形的精细特征, 需要在地块内部和边界带上合理建立足够空间密度的GNSS观测网络。如美国板块边界观测计划(PBO)在美国西部建立了由1 200个GNSS连续站组成的密集观测网络, 日本国土地理院在日本境内建立的GEONET网络由个1 400个GNSS连续站组成。中国地震科学试验场、国家重点研发计划等重大项目也已在青藏高原东南缘建立了较为密集的GNSS连续观测网络, 但在青藏高原内部, 如柴达木盆地、羌塘地块和巴颜喀拉地块东部等仍存在监测盲区, 今后需加强形变监测网络的建设。

(4)利用空间大地测量技术丰富和完善活动地块假说, 需要充分利用其它空间大地测量技术(如InSAR观测)和常规大地测量(如精密水准观测)的成果, 实现对重点及局部加密区域的有效补充监测, 使多种监测技术达到优势互补及有效融合。

致谢所有参加网络工程和陆态网络项目同仁的辛勤付出, 使得针对中国大陆地壳变形的精细化研究得以逐步发展和完善。评审专家提供了建设性意见, 使文章内容得到了大幅提高。在此一并表示感谢!本文是在大量前人研究基础上的概括和总结, 限于篇幅, 还有很多科研价值较高的相关学术文章未能参考引用, 敬请谅解。

参考文献

文献列表

[1]	柴艳菊, 欧吉坤, 独知行. 2002. 拟准检定法用于划分不同运动趋势块体[J]. 地震学报, 24(6): 579—586. CHAI Yan-ju, OU Ji-kun, DU Zhi-xing. 2002. A new approach for distinguishing different deformation trend blocks with displacement observations[J]. Acta Seismolgoica Sinica, 24(6): 579—586(in Chinese). [本文引用:1]
[2]	邓起东, 徐锡伟, 于贵华. 1994. 中国活动断裂分区特征 [A]∥ 中国地震学会地震地质专业委员会. 中国活动断层研究. 北京: 地震出版社: 1—14. DENG Qi-dong, XU Xi-wei, YU Gui-hua. 1994. Regional characteristics of active faults in China [A]∥ Research of Active Faults in China. Seismological Press, Beijing: 1—14(in Chinese). [本文引用:1]
[3]	丁国瑜(编). 1991. 中国岩石圈动力学概论 [M]. 北京: 地震出版社. DING Guo-yu(ed). 1991. Lithospheric Dynamics of China [M]. Seismological Press, Beijing(in Chinese). [本文引用:1]
[4]	甘卫军, 李强, 张锐, 等. 2012. 中国大陆构造环境监测网络的建设与应用[J]. 工程研究: 跨学科视野中的工程, 4(4): 324—331. GAN Wei-jun, LI Qiang, ZHANG Rui, et al. 2012. Construction and application of tectonic and environmental observation network of mainland China[J]. Journal of Engineering Studies, 4(4): 324—331(in Chinese). [本文引用:1]
[5]	郝明, 李煜航, 秦姗兰. 2017. 基于GPS数据的海原-六盘山断裂带滑动速率亏损时空分布[J]. 地震地质, 39(3): 471—484. HAO Ming, LI Yu-hang, QIN Shan-lan. 2017. Spatial and temporal distribution of slip rate deficit across Haiyuan-Liupan Shan fault zone constrained by GPS data[J]. Seismology and Geology, 39(3): 471—484(in Chinese). [本文引用:1]
[6]	李强, 游新兆, 杨少敏, 等. 2012. 中国大陆构造变形高精度大密度GPS监测: 现今速度场[J]. 中国科学(D辑), 42(5): 629—632. LI Qiang, YOU Xin-zhao, YANG Shao-min, et al. 2012. High precision and density of GPS velocity field in continental China[J]. Science in China(Ser D), 42(5): 629—632(in Chinese). [本文引用:1]
[7]	李延兴, 李智, 张静华, 等. 2004. 中国大陆及周边地区的水平应变场[J]. 地球物理学报, 47(2): 222—231. LI Yan-xing, LI Zhi, ZHANG Jing-hua, et al. 2004. Horizontal strain field in the Chinese mainland and its surrounding areas[J]. Chinese Journal of Geophysics, 47(2): 222—231(in Chinese). [本文引用:2]
[8]	李延兴, 杨国华, 李智, 等. 2003. 中国大陆活动地块的运动与应变状态[J]. 中国科学(D辑), 33(S1): 65—81. LI Yan-xing, YANG Guo-hua, LI Zhi, et al. 2003. Active tectonic block movement and strain state in continental China[J]. Science in China(Ser D), 33(S1): 65—81(in Chinese). [本文引用:3]
[9]	李煜航, 郝明, 季灵运, 等. 2014. 青藏高原东缘中南部主要活动断裂滑动速率及其地震矩亏损[J]. 地球物理学报, 57(4): 1062—1078. LI Yu-hang, HAO Ming, JI Ling-yun, et al. 2014. Fault slip rate and seismic moment deficit on major active faults in mid and south part of the eastern margin of Tibet Plateau[J]. Chinese Journal of Geophysics, 57(4): 1062—1078(in Chinese). [本文引用:2]
[10]	马杏垣. 1987. 中国岩石圈动力学概要[J]. 地质学报, 61(2): 113—125. MA Xing-yuan. 1987. Summary of the lithospheric dynamics in China[J]. Acta Geologica Sinica, 61(2): 113—125(in Chinese). [本文引用:1]
[11]	牛之俊, 马宗晋, 陈鑫连, 等. 2002. 中国地壳运动观测网络[J]. 大地测量与地球动力学, 22(3): 88—93. NIU Zhi-jun, MA Zong-jin, CHEN Xin-lian, et al. 2002. Crustal movement observation network of China[J]. Journal of Geodesy and Geodynamics, 22(3): 88—93(in Chinese). [本文引用:1]
[12]	欧吉坤. 1999. 粗差的拟准检定法(QUAD法)[J]. 测绘学报, 28(1): 15—20. OU Ji-kun. 1999. Quasi-accurate detection of gross errors(QUAD)[J]. Acta Geodaetica et Cartographica Sinica, 28(1): 15—20(in Chinese). [本文引用:1]
[13]	王敏. 2008. GPS观测结果的精化分析与中国大陆现今地壳形变场研究[D]. 北京: 中国地震局地质研究所. WANG Min. 2008. Analysis of GPS data with high precision and study on present-day crustal deformation in China [D]. Institute of Geology, China Earthquake Administration, Beijing(in Chinese). [本文引用:4]
[14]	王敏, 沈正康, 牛之俊, 等. 2003. 现今中国大陆地壳运动与活动块体模型[J]. 中国科学(D辑), 33(S1): 21—32. WANG Min, SHEN Zheng-kang, NIU Zhi-jun, et al. 2003. Present-day crustal movement and active block model of mainland China[J]. Science in China(Ser D), 33(S1): 21—32(in Chinese). [本文引用:2]
[15]	张国民, 张培震. 1999. 近年来大陆强震机理与预测研究的主要进展[J]. 中国基础科学, 2-4: 47—58. ZHANG Guo-min, ZHANG Pei-zhen. 1999. Main recent development in the research on the mechanism and forecast for continental strong earthquakes[J]. China Basic Science, 2-4: 47—58(in Chinese). [本文引用:2]
[16]	张国民, 张培震. 2000. “大陆强震机理与预测”中期学术进展[J]. 中国基础科学, 10: 4—10. ZHANG Guo-min, ZHANG Pei-zhen. 2000. Academic progress on “The Mechanism and Forecast for Continental Strong Earthquake”in the first two years[J]. China Basic Science, 10: 4—10(in Chinese). [本文引用:1]
[17]	张培震. 1999. 中国大陆岩石圈最新构造变动与地震灾害[J]. 第四纪研究, 39(5): 404—413. ZHANG Pei-zhen. 1999. Late Quaternary tectonic deformation and earthquake hazard in continental China[J]. Quaternary Sciences, 39(5): 404—413(in Chinese). [本文引用:1]
[18]	张培震, 邓起东, 张国民, 等. 2003. 中国大陆的强震活动与活动地块[J]. 中国科学(D辑), 33(S1): 12—20. ZHANG Pei-zhen, DENG Qi-dong, ZHANG Guo-min, et al. 2003. Active tectonic blocks and strong earthquakes in the continent of China[J]. Science in China(Ser D), 33(S1): 12—20(in Chinese). [本文引用:4]
[19]	张培震, 邓起东, 张竹琪, 等. 2013. 中国大陆的活动断裂、地震灾害及其动力过程[J]. 中国科学(D辑), 43(10): 1607—1620. ZHANG Pei-zhen, DENG Qi-dong, ZHANG Zhu-qi, et al. 2013. Active faults, seismic hazard and geodynamics of mainland China[J]. Science in China(Ser D), 43(10): 1607—1620(in Chinese). [本文引用:2]
[20]	张培震, 甘卫军, 沈正康, 等. 2005. 中国大陆现今构造作用的地块运动和连续变形耦合模型[J]. 地质学报, 79(6): 748—756. ZHANG Pei-zhen, GAN Wei-jun, SHEN Zheng-kang, et al. , 2005. A coupling model of rigid-block movement and continuous deformation: Patterns of the present-day deformation of China's continent and its vicinity[J]. Acta Geologica Sinica, 79(6): 748—756(in Chinese). [本文引用:7]
[21]	张培震, 沈正康, 王敏, 等. 2004. 青藏高原及周边现今构造变形的运动学[J]. 地震地质, 26(3): 367—377. ZHANG Pei-zhen, SHEN Zheng-kang, WANG Min, et al. 2004. Kinematics of present-day tectonic deformation of the Tibetan plateau and its vicinities[J]. Seismology and Geology, 26(3): 367—377(in Chinese). [本文引用:1]
[22]	张培震, 王琪, 马宗晋. 2002. 中国大陆现今构造运动的GPS速度场与活动地块[J]. 地学前缘, 9(2): 430—441. ZHANG Pei-zhen, WANG Qi, MA Zong-jin. 2002. GPS velocity field and active crustal blocks of contemporary tectonic deformation in continental China[J]. Earth Science Frontiers, 9(2): 430—441(in Chinese). [本文引用:2]
[23]	张文佑. 1984. 断块构造导论 [M]. 北京: 石油工业出版社. ZHANG Wen-you. 1984. An Introduction to Fault-block Tectonics [M]. Petroleum Industry Press, Beijing(in Chinese). [本文引用:1]
[24]	Deng Q D, Zhang P Z, Ran Y K, et al. 2003. Basic characteristics of active tectonics of China[J]. Science in China(Ser D), 46(4): 356—372. [本文引用:2]
[25]	Gan W, Zhang P, Shen Z, et al. 2007. Present-day crustal motion within the Tibetan plateau inferred from GPS measurements[J]. Journal of Geophysical Research: Solid Earth, 112(B8): B08416. [本文引用:2]
[26]	Hao M, Li Y H, Zhuang W Q. 2019. Crustal movement and strain distribution in East Asia revealed by GPS observations[J]. Scientific Reports, 9: 16797. [本文引用:2]
[27]	Jiang G, Xu X, Chen G, et al. 2015. Geodetic imaging of potential seismogenic asperities on the Xianshuihe-Anninghe-Zemuhe fault system, southwest China, with a new 3-D viscoelastic interseismic coupling model[J]. Journal of Geophysical Research: Solid Earth, 120(3): 1855—1873. [本文引用:1]
[28]	Liang S M, Gan W J, Shen C Z, et al. 2013. Three-dimensional velocity field of present-day crustal motion of the Tibetan plateau derived from GPS measurements[J]. Journal of Geophysical Research: Solid Earth, 118(10): 1—11. [本文引用:1]
[29]	Loveless J P, Meade B J. 2011. Partitioning of localized and diffuse deformation in the Tibetan plateau from joint inversions of geologic and geodetic observations[J]. Earth Planetary Science Letters, 303(1-2): 11—24. [本文引用:1]
[30]	McCaffrey R. 2005. Block kinematics of the Pacific-North America plate boundary in the southwestern United States from inversion of GPS, seismological, and geologic data[J]. Journal of Geophysical Research: Solid Earth, 110(B7): B07401. [本文引用:1]
[31]	Meade B J. 2007. Present-day kinematics at the India-Asia collision zone[J]. Geology, 35(1): 81—84. [本文引用:2]
[32]	Niu Z J, Wang M, Sun H R, et al. 2005. Contemporary velocity field of crustal movement of Chinese mainland from Global Positioning System measurements[J]. Chinese Science Bulletin, 50(9): 939—941. [本文引用:1]
[33]	Qu W, Lu Z, Zhang Q, et al. 2014. Kinematic model of crustal deformation of Fenwei Basin, China based on GPS observations[J]. Journal of Geodynamics, 75: 1—8. [本文引用:1]
[34]	Rui X, Stamps D S. 2019. A geodetic strain rate and tectonic velocity model for China[J]. Geochemistry, Geophysics, Geosystems, 20(3): 1280—1297. [本文引用:1]
[35]	Savage J C, Simpson R W. 2013a. Clustering of GPS velocities in the Mojave Block, southeastern California[J]. Journal of Geophysical Research: Solid Earth, 118(4): 1747—1759. [本文引用:1]
[36]	Savage J C, Simpson R W. 2013 b. Clustering of velocities in a GPS network spanning the Sierra Nevada Block, the northern Walker Lane belt, and the central Nevada seismic belt, California-Nevada[J]. Journal of Geophysical Research: Solid Earth, 118(9): 4937—4947. [本文引用:1]
[37]	Savage J C, Wells R E. 2015. Identifying block structure in the Pacific Northwest, USA[J]. Journal of Geophysical Research: Solid Earth, 120(11): 7905—7916. [本文引用:1]
[38]	Shen Z K, Lü J, Wang M, et al. 2005. Contemporary crustal deformation around the southeast borderland of the Tibetan plateau[J]. Journal of Geophysical Research: Solid Earth, 110(B11): B11409. [本文引用:3]
[39]	Shen Z K, Wang M, Zeng Y H, et al. 2015. Optimal interpolation of spatially discretized geodetic data[J]. Bulletin of Seismological Society of America, 105(4): 2117—2127. [本文引用:1]
[40]	Simpson R W, Thatcher W, Savage J C. 2012. Using cluster analysis to organize and explore regional GPS velocities[J]. Geophysical Research Letter, 39(18): L18307. [本文引用:1]
[41]	Thatcher W. 2007. Microplate model for the present-day deformation of Tibet[J]. Journal of Geophysical Research: Solid Earth, 112(B1): B01401. [本文引用:2]
[42]	Wang H, Liu M, Cao J, et al. 2011. Slip rates and seismic moment deficits on major active faults in mainland China[J]. Journal of Geophysical Research: Solid Earth, 116(B2): B02405. [本文引用:2]
[43]	Wang M, Shen Z K. 2020. Present-day crustal deformation of continental China derived from GPS and its tectonic implications[J]. Journal of Geophysical Research: Solid Earth, 125(2): e2019JB018774. [本文引用:13]
[44]	Wang Q, Zhang P Z, Freymueller J T, et al. 2001. Present-day crustal deformation in China constrained by Global Positioning System(GPS)measurements[J]. Science, 294(5542): 574—577. [本文引用:1]
[45]	Wang W, Qiao X J, Yang S M, et al. 2017. Present-day velocity field and block kinematics of Tibetan plateau from GPS measurements[J]. Geophysical Journal International, 208(2): 1088—1102. [本文引用:2]
[46]	Zhang P Z. 2004. Late Cenozoic tectonic deformation in the Tianshan Mountain and its foreland basins[J]. Chinese Science Bulletin, 49(4): 311—313. [本文引用:3]
[47]	Zhang P Z, Gan W J. 2008. Combined model of rigid-block motion with continuous deformation: Patterns of present-day deformation in continental China[J]. Special Paper of the Geological Society of America, 444: 59—71. [本文引用:2]
[48]	Zhang P Z, Shen Z K, Wang M, et al. 2004. Continuous deformation of the Tibetan plateau from Global Positioning System data[J]. Geology, 32(9): 809—812. [本文引用:1]
[49]	Zhang Z Q, McCaffrey R, Zhang P Z. 2013. Relative motion across the eastern Tibetan plateau: Contributions from faulting, internal strain and rotation rates[J]. Tectonophysics, 584: 240—256. [本文引用:2]
[50]	Zheng G, Wang H, Wright T J, et al. 2017. Crustal deformation in the India-Eurasia collision zone from 25 years of GPS measurements[J]. Journal of Geophysical Research: Solid Earth, 122(11): 9290—9312. [本文引用:1]

2002

0.0

... 将该方法应用于地块划分, 即是将地块内的相对稳定点作为拟准观测, 而将不符合运动趋势的点看作粗差并予以剔除(柴艳菊等, 2002) ...

1994

0.0

... 随着板块构造学说被引入中国, 老一辈中国科学家针对中国大陆的板块构造体系和特征提出了断块构造(张文佑, 1984)、活动亚板块(马杏垣, 1987)、构造块体(丁国瑜, 1991)以及活动断块(邓起东等, 1994)等理论或划分方案 ...

1991

0.0

2012

0.0

... 建设的网络工程(牛之俊等, 2002)和陆态网络(甘卫军等, 2012)项目 ...

2017

0.0

... 郝明等, 2017), 为预测中长期强震的危险性提供了基础数据 ...

2012

0.0

... 李强等, 2012 ...

2004

0.0

李延兴, 李智, 张静华, 等. 2004. 中国大陆及周边地区的水平应变场[J]. 地球物理学报, 47(2): 222—231.

Yan-xing

, LI

Zhi

, ZHANG

Jing-hua

, et al. 2004. Horizontal strain field in the Chinese mainland and its surrounding areas[J]. Chinese Journal of Geophysics, 47(2): 222—231(in Chinese).

AbstractIn this paper, two kinds of equations of block elastic motion are established. One is the equation of rotation in the entire block and homogeneous strain, and the other is the equation of rotation in the entire block and linear strain. Using the velocities of GPS sites located in the Western Region, QinghaiXizang and North China blocks, these two kinds of elastic motion equations and the rigid motion equation are tested statistically for their unbiasedness and validity to simulate the site velocities of the blocks. The results indicate that the equation with rotation in the entire block and linear strain is the optimal model to describe the block motion. Dividing the Chinese mainland into 10 blocks, the rotation and strain parameters of each block are estimated

1 First Monitoring Center, China Seismological Bureau, Tianjin 300180, China 2 Technological College of Equipments and Commands, Beijing 101416, China 3 Shanghai Astronomical Observatory, Chinese Academy of Sciences, Shanghai 200030, China 4 Center for Analysis and Prediction, China Seismological Bureau, Beijing 100036, China

推导并建立了块体的两种弹性运动方程：块体的整体旋转与均匀应变方程和块体的整体旋转与线性应变方程. 应用统计学原理，使用西域、青藏和华北块体上的GPS站速度数据，对这两种弹性运动方程与刚体运动方程模拟块体站速度的无偏性和有效性进行了统计检验. 检验结果表明，块体的整体旋转与线性应变方程是描述块体运动的最优模型. 将中国大陆划分为10个块体，应用块体的整体旋转与线性应变方程和块体上的GPS站速度估计了各个块体上的旋转与应变参数，按照1°×1°的间距计算了中国大陆及周边地区上1005个点的应变参数，分析了中国大陆及周边地区应变场的基本特征. 用本文方法得到的主压应变方向与地质学方法和测震学方法得到的主压应力轴方向具有很好的一致性(华南块体除外).

... 3 基于GNSS技术检验活动地块的刚性与非刚性变形地块在周围板块或地块的作用下不仅会发生刚性运动(平动和旋转), 其内部可能会产生均匀或非均匀的线性应变变形(李延兴等, 2003, 2004) ...

... 如果地块内部的应变不是均匀的而是随位置线性变化的, 则地块的刚性运动和内部线性应变变形可以用式(6)表示(李延兴等, 2004 ...

2003

0.0

... 李延兴等, 2003 ...

... 球面上刚性块体的运动可以用绕欧拉极的旋转(李延兴等, 2003)来描述: ...

2014

0.0

李煜航, 郝明, 季灵运, 等. 2014. 青藏高原东缘中南部主要活动断裂滑动速率及其地震矩亏损[J]. 地球物理学报, 57(4): 1062—1078.

Yu-hang

, HAO

Ming

, JI

Ling-yun

, et al. 2014. Fault slip rate and seismic moment deficit on major active faults in mid and south part of the eastern margin of Tibet Plateau[J]. Chinese Journal of Geophysics, 57(4): 1062—1078(in Chinese).

As one of the most important pathways for southeastward extrusion of Tibet plateau, the mid and south part of the eastern margin of Tibet plateau feature widespread active faults, intensive and frequent seismic activity. In this paper, we use the spherical linear block theory constrained by geodetic observations from Global Positioning System (GPS) stations to estimate fault slip rates in the mid and south part of the eastern margin of Tibet plateau. The 3D geometric block model is based on the previous researches of active block, considering the latest research works in this field. Using these slip rates and optimal locking depths, the rates of moment accumulation on the major faults are estimated. And then we compared them with the seismic moment released on each fault using the historical earthquake catalog. The result shows that there are six faults or fault segments with large moment deficits in the research area; it may be helpful for the middle-long term strong earthquake hazard study.

[Li Yu-Hang] China Earthquake Adm, Inst Geol, Beijing 100029, Peoples R China.;[Li Yu-Hang; Hao Ming; Ji Ling-Yun; Qin Shan-Lan] China Earthquake Adm, Crust Monitoring & Applicat Ctr 2, Xian 710054, Peoples R China.

作为青藏高原南东向“挤出逃逸”的重要通道，青藏高原东缘中南部具有大型走滑断裂广泛发育和地震活动强烈而频繁的特征.本文使用线性球面块体模型理论，在前人活动地块研究的基础上吸收新近研究成果，建立研究区三维块体几何模型，使用1999—2007年的GPS数据反演得到青藏高原东缘中南部主要活动断裂滑动速率.使用反演得到的滑动速率和最优断层闭锁深度估算了川滇菱形块体主要边界和其内部断裂的地震矩积累，并利用历史强震目录估算了地震矩释放，在比较两者差异的基础上得到了研究区地震矩亏损（未释放的地震矩）较大的断层和断层段，该结果可以作为研究区强震中长期危险性研究的参考.

... 李煜航等, 2014 ...

1987

0.0

2002

0.0

... 建设的网络工程(牛之俊等, 2002)和陆态网络(甘卫军等, 2012)项目 ...

1999

0.0

... 2 拟准检定法拟准检定法(欧吉坤, 1999)是在#cod#x0201c ...

2008

0.0

... 王敏, 2008) ...

... 如果地块的边界断裂不发生闭锁, 则地块的刚性运动和内部均匀应变变形可用式(5)表示(王敏, 2008 ...

... 此外, 同样可以用F检验方法来判定某个地块为刚性平动、旋转变形或非刚性的连续变形(王敏, 2008) ...

... 王敏, 2008 ...

2003

0.0

... 王敏等, 2003 ...

... (3)检验相邻2个地块的独立性(王敏等, 2003 ...

1999

0.0

... 0 引言20世纪60年代建立并发展起来的板块构造学说合理解释了全球80%的地震发生在板块边界的现象, 但由于板内地震和板块边界地震在动力环境、构造成因、变形机制和破坏类型等方面存在很多显著差别, 板块构造学说在解释大陆板内地震时遇到了诸多难以解释的问题(张国民等, 1999) ...

... 的科学假说, 用于解释中国大陆强震的空间分布(张国民等, 1999, 2000 ...

2000

0.0

... 的科学假说, 用于解释中国大陆强震的空间分布(张国民等, 1999, 2000 ...

1999

0.0

... 张培震, 1999 ...

2003

0.0

... 张培震等, 2003) ...

... 张培震等(2003)明确了活动地块的定义, 即活动地块是被形成于晚新生代、晚第四纪至现今强烈活动的构造带所分割和围限的、具有相对统一运动方式的地质单元 ...

... 4 基于GNSS技术获得的中国大陆主要活动地块的运动与变形张培震等(2003)将中国大陆及其周边地区划分为6个I级地块区和22个#cod#x02161 ...

... 5 总结与展望中国大陆晚新生代的构造变形主要表现为活动地块的运动学特征, 不同活动地块的运动方式也不相同, 地块之间相互作用的地块边界带构造变形最为强烈, 也是强震高发区(张培震等, 2003) ...

2013

0.0

... 总体而言, 中国大陆的活动地块有的刚性较好, 内部不发生构造变形和大地震, 有的地块则刚性很差, 内部发生明显的构造变形并伴随强震发生(张培震等, 2013) ...

... 已有的研究结果表明, 中国大陆由2类活动地块所组成(张培震等, 2005, 2013 ...

2005

0.0

... , 2020), 揭示了中国大陆现今的总体构造变形具有刚性和非刚性的运动学特征(张培震等, 2002, 2004, 2005 ...

... 因此, 利用GNSS观测技术建立的现今高精度、高时空分辨率的构造变形运动图像, 可为研究划分活动地块提供定量约束(张培震等, 2005) ...

... 相邻刚性块体之间的相对运动和变形主要集中在块体边界带上, 块体内部不发生大规模变形(张培震等, 2005) ...

... 8mm/a的速率向N运动, 属典型的刚性块体(张培震等, 2005 ...

... 但鄂尔多斯地块内部构造的活动性微弱, 也不发育大规模的活动断层(张培震等, 2005) ...

... 华南和东南沿海地区在新构造运动上属于比较稳定的地块, 内部不发育明显的活动断裂和褶皱, 地震活动性与华北和西部相比相对较弱(张培震等, 2005) ...

... 已有的研究结果表明, 中国大陆由2类活动地块所组成(张培震等, 2005, 2013 ...

2004

0.0

张培震, 沈正康, 王敏, 等. 2004. 青藏高原及周边现今构造变形的运动学[J]. 地震地质, 26(3): 367—377.

ZHANG

Pei-zhen

, SHEN

Zheng-kang

, WANG

Min

, et al. 2004. Kinematics of present-day tectonic deformation of the Tibetan plateau and its vicinities[J]. Seismology and Geology, 26(3): 367—377(in Chinese).

As the most prominent example of large scale continental deformation,Tibetan Plateau offers an ideal natural laboratory for quantifying such deformation and understanding the relevant dynamic processes. Global Positioning System (GPS) provides a powerful means to directly measure the kinematics of present day deformation. Our synthesis of GPS velocities from 553 stations in Tibetan Plateau and its margins quantitatively show that most of the relative India/Eurasia motion has been accommodated primarily by crustal shortening along the margins,strike slip and normal faulting in the plateau interior,and clockwise rotation around the eastern end of Himalayas. Taken 36～40mm/yr as total relative motion between India and Eurasia,the eastern Tibetan Plateau and its margins accommodates 85%～94% of the total motion,whereas western Tibet absorbs 70%～91% of total convergence and the rests are taken up by shortening across the Tianshan in the north. The NNE SSW shortening of the plateau interior is accommodated by conjugate strike slip faulting and orthogonal normal faulting,which do not require crustal thickening or thrust faulting. The eastward extrusion of Tibetan Plateau out of India's northward pass is carried out by roughly eastward flow of crustal material rather than by rigid block rotation. The flow of Tibetan crustal material rotates around the eastern Himalayan syntaxis,causing southeastward to southward and even southwestward velocities observed in southern and western Yunnan Province of China. To the east,the eastward flow of crustal material causes shortening across the eastern margin of the plateau and clockwise rotations of the region where resistance to such flow is weak. To the west,the westward motion of the western margin of the plateau is observed with only 4mm/yr slip rate. Components of velocity of Tibetan Plateau in both parallel and perpendicular directions to the relative motion between India and Eurasia can not be attributed to slips along a few faults either strike slip or thrust because of their distributed nature. Thus,the present day tectonics in the Tibetan Plateau is best described as deformation of a continuous medium,at least when averaged over distances of ～100km.

1. State Key Laboratoray of Earthquake Dynamics, Institute of Geology, China Earthuake Administration, Beijing 100029, China; 2. Institute of Earthquake Prediction, China Earthquake Administration, Beijing 100036, China

青藏高原现今构造变形的定量化研究是理解其动力过程的基础,近年来高速发展的GPS(全球定位系统)技术为测量大尺度现今构造变形提供了最有效的手段。我们利用青藏高原及周边的553个GPS观测数据给出了其现今构造变形的速度场,表明印度和欧亚板块之间的相对运动主要被青藏高原周边的地壳缩短和内部的走滑剪切所调整吸收。其中,喜马拉雅山系吸收了青藏高原总缩短量的44%～53%,北部的阿尔金山、祁连山和柴达木盆地吸收了15%～17%,高原内部吸收了32%～41%。青藏高原的“向东挤出”实际上是地壳物质的向东流动而不是刚性地块的挤出。这一地壳物质流动带在高原西部以地表张性正断层和共轭剪切走滑断层为特征,到高原中东部转换为巨型的弧形走滑断裂带,再到高原东北缘转换为地壳缩短和绕东喜马拉雅构造结的顺时针旋转。青藏高原的大尺度现今构造变形以连续变形为特征

... , 2020), 揭示了中国大陆现今的总体构造变形具有刚性和非刚性的运动学特征(张培震等, 2002, 2004, 2005 ...

2002

0.0

... , 2020), 揭示了中国大陆现今的总体构造变形具有刚性和非刚性的运动学特征(张培震等, 2002, 2004, 2005 ...

... 横跨天山的SN向GNSS剖面揭示地壳缩短是均匀的, 而不是集中在两侧的前陆冲断带或某一断层附近(张培震等, 2002) ...

1984

0.0

2003

0.0

... 夹持在构造活动十分强烈的天山和青藏高原之间的塔里木地块的构造活动性相对较弱, 有历史记载以来地块内部没有发生过强震, 也不发育晚第四纪活动断裂(Deng et al ...

... 华北平原活动地块从长期、平均构造变形的角度来看不具有刚性块体的性质, 其内部发育一系列活动断裂, 并伴有大量历史强震发生(Deng et al ...

2007

0.0

... Gan et al ...

... Zhang 等(2004)和Gan等(2007)利用青藏高原及其周边的GNSS观测资料, 发现印度和欧亚板块相对运动速率的70%~94%被青藏高原内部的构造变形所调整吸收 ...

2019

0.0

... Hao et al ...

... Hao等(2019)选择中国大陆内部稳定的华南地块作为区域参考基准, 所提供的东亚地区的GNSS运动速度场图像清晰地展示出中国大陆受到印度板块的推挤, 并在太平洋和菲律宾板块的俯冲作用下其内部产生了差异性构造变形特征 ...

2015

0.0

... Jiang et al ...

2013

0.0

... Liang et al ...

2011

0.0

... Loveless et al ...

2005

0.0

... 如果考虑断层闭锁, 块体模型可用块体欧拉旋转、块体内部均匀应变以及块体边界由于断层闭锁产生的地表弹性变形3个分量描述(McCaffrey, 2005) ...

2007

0.0

... 以GNSS运动速度场为约束, 采用块体模型可估计出活动地块边界带现今的滑动速率和闭锁强度(Meade, 2007 ...

... Meade, 2007 ...

2005

0.0

... (简称陆态网络)项目的开展, 利用丰富的GNSS观测数据精细地刻画出中国大陆现今地壳的水平运动图像(Niu et al ...

2014

0.0

... Qu et al ...

2019

0.0

... Rui et al ...

2013

0.0

... Savage等(2013a, b, 2015)将k-中心点聚类分析方法应用于圣安地列斯断裂和美国西北部地区的地块划分 ...

2013

0.0

... Savage等(2013a, b, 2015)将k-中心点聚类分析方法应用于圣安地列斯断裂和美国西北部地区的地块划分 ...

2015

0.0

... Savage等(2013a, b, 2015)将k-中心点聚类分析方法应用于圣安地列斯断裂和美国西北部地区的地块划分 ...

2005

0.0

... Shen et al ...

... (2)检验每个测站与其它测站的兼容性(Shen et al ...

... 随后, 很多学者以此为基础模型, 采用时空分辨率不断提高的GNSS速度场结果, 愈来愈细致地揭示出中国大陆及周边活动地块的运动学特征(Shen et al ...

2015

0.0

... Shen et al ...

2012

0.0

... Simpson等(2012)将层次凝聚聚类(HAC)方法应用于加州旧金山湾地区的GNSS速度场, 试图找出该地区是否具有类似地块变形的特征 ...

2007

0.0

... Thatcher, 2007 ...

2011

0.0

... Wang et al ...

2020

0.0

... Wang et al ...

... 在Wang 等(2020)提供的中国大陆及周边最新的GNSS运动速度场结果中还使用了一些区域网络GNSS数据, 这些区域网络主要布设在活动构造带, 如华北、汾渭断陷带和川滇地区 ...

... Wang (2020)针对上述问题详细地阐述了消除大地震的同震和震后变形的方法 ...

... Wang 等(2020)给出的高时空分辨率的GNSS运动速度场图像整体反映出中国大陆呈顺时针旋转的运动状态, 青藏高原内部的一部分地壳物质沿NE向流动迁移, 一部分沿SE向绕喜马拉雅东构造结发生顺时针旋转, 似乎还有一小部分近E向挤入华南地块, 华北和华南地块的差异性运动并不十分明显 ...

... 因此, 本文将Wang 等(2020)提供的相对于欧亚板块的GNSS速度场旋转为相对于华南地块的GNSS速度场(图 2) ...

... Wang et al ...

... Wang 等(2020)采用刚性和连续应变模型(式(5)), 利用中国大陆的GNSS速度场拟合了东北、华北、华南、鄂尔多斯、阿拉善、塔里木和准噶尔地区的运动变形, 结果表明这几个地区具有刚性地块运动特征, 但地块内部的变形(10^-9/a量级)也不可忽视(图 4), 而青藏高原和天山地区则具有明显的连续变形特征 ...

... Wang et al ...

... 10^-9/a(Wang et al ...

... /Ma也很小, 地块内部变形微弱(Wang et al ...

... 10^-9/a的内部主应变率(Wang et al ...

2001

0.0

... Wang 等(2001)首次应用GNSS数据研究了中国大陆现今的地壳运动和变形, 其结果清晰地反映了以活动地块为单元的分块运动特征 ...

2017

0.0

... Wang et al ...

2004

0.0

... Zhang et al ...

... 天山南、北两侧分别被构造稳定的准噶尔和塔里木地块所夹持, 构造变形主要发生在天山山前和山体内部(Zhang, 2004) ...

... 的天山发生顺时针旋转所导致的地壳缩短变化来解释(Zhang, 2004) ...

2008

0.0

... 地块的内部变形和旋转量非常显著, 青藏高原内部各地块均是可变形的非刚性活动地块(Zhang et al ...

... Zhang et al ...

2004

0.0

... Zhang 等(2004)和Gan等(2007)利用青藏高原及其周边的GNSS观测资料, 发现印度和欧亚板块相对运动速率的70%~94%被青藏高原内部的构造变形所调整吸收 ...

2013

0.0

... Zhang et al ...

... Zhang 等(2013)利用GNSS资料建立了青藏高原东缘的可变形地块模型, 结果表明青藏高原地块内部的变形明显大于其周边的华南地块和阿拉善地块 ...

2017

0.0

... Zheng et al ...