南北地震带近期地壳水平运动变化分析摘要：对2009、2011及2013年3期陆态网络流动GPS观测资料进行了处理，获得了2009～2011及2011～2013年两个时间段南北地震带地壳水平运动速度图像。参考1999～2007年网络工程流动GPS速度场结果，分析了南北地震带近期地壳水平运动场与应变场变化特征，认为近期地壳运动变化受汶川地震震后影响较大。2013年芦山7.O级地震和岷县漳县6.6级地震是2008年汶川8.0级地震震后应力场调整触发的结果。关键词：南北地震带；GPS；应变场；汶川8.O级地震；芦山7.O级地震；岷县漳县6.6级地震：P315.7文献标识码：A：1000-0666（2015）02-0229-080引言中国大陆中部发育着著名的南北向构造带，由于其控制着一系列历史强震的发生也称为南北地震带，该带北起鄂尔多斯地块西缘，跨越秦岭、穿过龙门山，再沿鲜水河一安宁河一小江断裂带向南延伸至缅甸境内，成为分隔中国大陆东部相对稳定的鄂尔多斯高原、四川盆地和华南地块与西部强烈隆升的青藏高原之间的边界活动构造带（张培震，2008）。深部地球物理探测和人地测量及地形地貌表明，该带不仪是地壳厚度梯度带和重力梯度带（马杏垣，1978），也是中国大陆现今地壳运动和历史构造活动最为强烈的区域之一。自2008年南北地震带中段的龙门山断裂发生汶川8.0级地震后，该区域地壳形变发生了显著的变化。许多学者从不同的时间和空间尺度对该区域近期地壳运动特征进行了分析与研究（蒋锋云等，2013；唐红涛等，2014；方颖等，2014）。在分析该区域构造形变背景特征的基础上，本文主要分析了利用流动GPS观测获得的汶川地震震后两个时段地壳水平速度场及计算的应变场的变化特征。1资料及处理资料主要于2009、2011及2013年“中国大陆构造环境网络”项目（简称“陆态网络”）在南北地震带观测的GPS数据。采用GAMIT-GLOBK软件对数据进行处理：首先利用GAMIT软件获得GPS测站坐标、卫星轨道及台站对流层天顶延迟等参数及其方差协方差矩阵的单日松弛解（H文件），解算时采用SOPAC数据中心提供的精密星历文件，并将中国人陆及周边分布均匀、稳定性较好的19个国际IGS永久跟踪站一并参与解算，对这些IGS站坐标根据实际情况给出合理的坐标约束，对区域站则进行较为宽松的约束。其次，利用GLOBK软件（卡尔曼滤波）进行多时段综合解，获得网平差结果。在此过程中，我们把解算获得的区域站单日松弛解和全球IGS站单日松弛解捆绑在一起，并选取全球均匀分布的90多个稳定站为基准点，将区域站速度场同定到ITRF2008框架下。最后通过基准变化，获得相对稳定欧亚板块的2009～2011年及2011～2013年两个时间段速度场。1999～2007年速度场主要采用郑文俊等（2009）给出的结果。2水平运动速度场及应变率场2.1水平运动特征分析为了便于后续的分析，将区域的主要断裂及发生的主要地震在图1中标出。断裂分布位置以邓起东等（2002）给出的断裂为主要依据。图2为南北地震带3个不同时间段GPS水平运动速度场，总体来看其运动趋势较为一致。从构造上来看，由于印度板块向欧亚板块俯冲，使得青藏高原向北、向东扩展，在这一过程中，由于受到周边稳定阿拉善地块、鄂尔多斯地块及华南地块的阻挡，其运动主要被高原内部及边缘的地壳缩短增厚（Englandetal.，1986）和一系列火型走滑断裂切割形成的块体的挤出、旋转（Tap-ponnieretal，1982）所调整和吸收。鉴于“网络T程”建立后，1999～2007年南北地震带没有发生过Ms≥7.0、尤其是8级特大地震，因此可将汶川地震前的1999～2007年较长时间段GPS水平运动速度场近似看作研究区水平运动的本底趋势形态来分析其运动特征，而用震后的2009～2011年及2011～2013年的运动状态来分析与大地震发生有关的的水平运动变化的时间阶段性特征。从图2a可以看出，现今水平运动的空间差异性主要体现在运动强度和运动方式上：（1）不同构造部位运动强度的差异性主要表现在青藏高原与华南地块、鄂尔多斯地块和阿拉善地块交界地带地壳水平运动速率较小，而在高原内部运动速率较人，反映了青藏高原的强烈活动与周边相对稳定地块活动的差异。如：青藏高原东北缘的柴达木一祁连山地区GPS速度矢量由高原内部的15～20mm/a减小到高原边界地带的5...