利用GPS观测结果反演中国大陆及邻区现代构造变形运动.pdfVIP

第45卷 增刊 地 球 物 理 学 报 Vo1．45，Supp1． 2002焦 CHINESE JOURNAL OF GEOPHYSICS 2oo2 [文章编号】 0001—5733(2002)增一0123—13 [中图分类号】 P546 利用 GPS观测结果反演中国大陆及邻区 现代构造变形运动 任金卫 (中国地震局地质研究所，北京 100029) [摘 要】 通过模拟水平应变率场获得了中国大陆地区相对于稳定欧亚大陆的现今地 壳运动速度场．用此方法得出的应变率解相当于一个牛顿体薄壳在水平作用力下的反 应．在模型中印度板块和阿拉伯板块相对于欧亚板块的速率采用了由NUVEL一1A模型的 欧拉矢量所定义的速度值作为约束，对每一个单元格中的模型应变率张量通过地震矩张 量和第四纪活动断层滑动速率加以约束．模拟结果表明，在各向同性的不均匀体中由NU． VEL一1A模型所定义的印度板块相对于欧亚板块的运动绝大部分都被近南北向的地壳缩 短变形所吸收，走滑运动只占了非常小的比例．在各向异性模型中，走滑运动分量增大， 但最大不超过总量的25％，说明在调节整个青藏高原构造变形的过程中，逆断和地壳增 厚起了主要的作用． [关键词】 地壳变形，GPS，中国大陆． 1 引 言 由印度大陆向北楔入引起的南北向缩短主要是通过地壳增厚和走滑断层作用这两种方 式来吸收的，在这一点上研究青藏高原的科学家们基本上取得了一致的看法，但在这两种作 用谁占主导地位的问题上仍然存在着较大的分歧．与此同时，在具体工作中上述争论的焦点 集中在青藏高原主体侧向挤出量的估计上和这种侧向挤出在调节构造变形中的作用．本文 试图以中国大陆GPS观测结果为基础，通过数值方法求解中国大陆及邻区地壳水平运动速 度梯度张量场．研究的范围南至青藏高原南侧的板块碰撞边界和南中国海，北至贝加尔湖以 北，西至咸海以东的中亚地区，东侧以日本海以西为界．同时包括所涉及的亚洲、印度和阿拉 伯板块的边界带，目的是要确定由印度板块向北推挤所产生的内陆构造变形是如何分布的． 本文将利用不规则网格来勾画板块边界带、主要构造带和主要地质构造单元，使其尽可能符 合实际的几何形状，以便能合理地确定板块边界和内部变形的分布，从而更好地了解不同构 造单元在调节板块运动变形方面所起的作用，在此基础上回答与青藏高原物质侧向挤出量 估计有关的一系列问题，并探讨其可能的动力学驱动机制． [收稿日期】 2001—11—21收到，2002—05—18收到修定稿． [基金项目】 国家重点基础研究发展规划项目(G1998040703)，国家攀登计划项目(970231003)和地震联合基金课题资 助(198015)． 【作者简介】任金卫，男，1957年生，1982年南京大学地质学系毕业，研究员，主要从事活动构造研究和利用 GPS观测 资料研究现代地壳运动．E—tlulilI jwren@163bj．eOIll 地 球 物 理 学 报 45卷 2 方法和资料 Haines和Holt…提出了利用地震资料反演地壳运动速度场的方法．本项研究采用改 进的Haines 数学反演方法，它利用不规则网格勾画构造带和板块边界带，使其几何形状 得到尽可能符合实际的描述．此外，改进后的方法还容许在黏性体中加入刚性块体的运 动，这不仅使研究区内相对刚性的块体的运动得到恰当的描述，也大大减少了计算量．在 整个中国大陆及邻近地区确定其地壳运动速度场时，需要在求解整个模型区应变率分布 最佳拟合的过程中，按照NUVEL一1A模型 J，加入印度板块和阿拉伯板块与欧亚板块之 间的相对运动速度矢量作为边界速度条件．因为所获得的连续的速度场是自恰的，在所拟 合的应变率场中，通过所有可能的积分路径，自稳定欧亚地块到相应的印度和阿拉伯板块 边界，所获得的速度矢量应该与该板块由NUVE