用全球定位系统（GPS）研究中国大陆现今地壳运动模式

摘要

中国大陆受周缘板块的挤压、各块体间的相互作用影响，致使其内部聚积了特别强大的构造作用力，历史上成为板内构造最活跃、板内地震最强烈的地区，构造变形覆盖了印欧碰撞带以北两到三千公里的范围。正确理解大陆构造变形运动过程及其动力机制，是探求地震成因和进行中、长期或短临地震预报的首要前提。因此，必须首先获取分辨精细、精度一致的岩石圈构造变动速度场，而岩石圈构造变动速度场的获取必须以质量均匀、分布广泛、密度适当的地壳变形和应变的实测资料作基础。迄今为止，对此运动、动力过程的认识基本来自第四纪以来主要活动断层大量观测调查及百余年来地震矩张量反演。中国大陆的晚新生代构造变形分布十分广泛，获取第四纪活动断裂和活动褶皱的定量数据需要漫长的时间和大量的投入，短时间内很难完成。地质学、地震学对地壳整体运动研究仍是一种间接推算，受资料来源及分析方法上制约，其准确性与可靠性在短期内难有根本提高，对现今地壳运动形态的精细观测更难有作为。与地质、地震资料相比，GPS观测结果具有量化清晰、精度高，时间尺度一致，动力学意义明确等优点，逐步成为岩石圈大尺度构造运动与动力学研究的重要基础资料。 自上世纪八十年代开始，中国大陆现今地壳运动的GPS监测资料日益丰富，一张全面反映中国大陆地壳运动的速度图像已经形成，使利用GPS测定的站点速率，研究构造变形的模式速度场成为可能。本文继续了前期工作，主要用“中国地壳运动观测网络”实测速率值（1200个测站），通过数值模拟技术构建了南至青藏高原南侧的印欧板块碰撞边界和南中国海，北至贝加尔湖，西至咸海以东的中亚地区，东以日本海以西为界的中国大陆构造变形的整体速度场。为此本文针对构造变形是以连续还是以分块的方式进行了对比研究。 首先在地块细分的多样性和块体边界的不确定性条件下，中国大陆的构造变形运动虽然块体或断裂带之间的非连续性是绝对的，可当假定中国大陆运动在重力驱动下类似流体缓慢流动，认为可将变形描述为连续粘性介质运动。基于此，本文利用球面上双三次样条函数数值模拟方法，将离散的值在空间域上连续化，来求解中国大陆的构造变形场，应变场。 然后本文根据块体、非连续变形假定反演中国大陆GPS速度场，结合中国大陆活动构造背景，建立地块运动模型。以整个变形区域内的主要走滑断层、正断层和逆断层为边界细致划分块体来同时模拟弹性岩石圈块体的旋转和边界断层闭锁产生的应变。给出了33个块体的细致模型和29个块体的简化模型；得到了各主要断层的活动速率、地块运动速率、欧拉矢量以及应变率。 两种模拟GPS观测结果显示： (1) 印度板块与欧亚板块的碰撞、挤压是构成中国大陆内部岩石层水平形变的主要驱动力；对中国大陆及周边地区构造活动的连续变形分析还表明GPS推算的数年尺度应变场的平均最大、最小水平主压应变率方向与长期地质资料、百年尺度地震资料展示的应力状态基本吻合； (2) 绕东喜马拉雅构造结，模型速度场和主应力方向有近似180o的顺时针旋转，说明中国东南部的应力场主要受到了周边特定板块边界的作用； (3) 太平洋板块和菲律宾板块作用力虽然比印度板块的作用力小，但仍然是影响中国大陆东部的重要因素； (4) 以玛尼-玉树-甘孜-理塘断裂带为中心构成了一条现今地壳物质东流带，并且青藏高原物质向东的挤出或流动没有超出其东边界，而是在东边界一带转化成绕东喜马拉雅构造结的顺时针旋转； (5) 青藏高原内部以连续地壳变形或均匀地壳缩短为现今构造变形的主要特征，印度和欧亚板块之间相对运动量的90％被青藏高原的现今地壳缩短速率所吸收和调节,且西藏中部东西存在东西向的拉张,其拉张速率为16mm/yr左右； (6) 阿尔金断裂带的左旋走滑速率是5.8 ± 1.5 mm/yr ，祁连山和阿拉善之间的左旋走滑速率是6.5±2.0mm/yr，横穿喜马拉雅和青藏的南北向压缩速率为19.0+/-2.0 mm/yr，穿过西天山的南北向压缩速率为14.0 ± 1.1 mm/yr，穿过中天山的南北向压缩速率为7.9 ±1.2mm/yr； (7) 塔里木块体围绕欧拉极（37.39°N, 96.4°E,-0.533°/Ma）相对西伯利亚作顺时针旋转，华南块体则围绕欧拉极（63.38°N,159.87°W）作轻微的顺逆时针方向旋转，鄂尔多斯块体围绕欧拉极（49.06°N,118.51°E，0.213°/Ma）相对西伯利亚作逆时针方向旋转。 从统计结果和GPS残差图来看两种模型在中国大陆范围内都得到了较好的模拟效果。但连续模型似乎得到更好拟合效果，但这并不表明分块模型不能适用于大陆构造研究。实际上，本文得到的块体滑动速率、活动断裂滑动速率说明无论青藏高原，还是整个中国大陆，现今活动构造变形仍可以视为有限制的较低速率岩石圈和块体运动，不过进一步比较有待今后更多的GPS数据和地质、地震资料。