昆仑8.1级地震对青藏块体东北缘地壳水平运动的影响初探

为了探讨昆仑8.1级地震引起的同震形变是否影响了青藏块体东北缘的地壳水平运动,以震后GPS观测资料为数据源,依Okada位错模型原理,计算了地震引起的同震水平位移场,分析了地震对区域地壳水平运动的影响。结果表明:垂直于断层沿震中方向的GPS点受地震影响大,远离断层的测点的同震位移受断层位错的影响小。     

超导重力数据检测到的2011年日本东北大地震(Mw 9.0)震前重力异常及同震重力变化

高精度的超导重力数据已广泛应用于地球动力学的研究,对大地震前重力异常和同震重力变化的探测有助于震源机制和地震预警的研究,同时高精度的同震观测数据可用于断层滑动模型的反演。本文利用日本、中国及欧洲7个超导台站2011年3月的秒采样数据研究日本东北大地震（M_w 9.0）产生的重力变化,经潮汐、气压、漂移等预处理改正后,得到改正后的重力变化。对比同时段日本岛附近发生的2861次M_b≥4级地震,分析滤波后的重力变化数据,排除非主震的影响,发现0.12 Hz≤f≤0.18 Hz频段显示了明显的震前重力异常扰动信息,所有台站在震前89 h均开始出现重力异常扰动,Medicina站的震前最大异常扰动振幅达到28×10^-8 m/s²。另外,利用球形位错理论及地震的CMT解计算得到的地表同震重力变化理论值与超导重力观测值非常接近,同震重力变化绝对值与震中距成反比,近场站点的观测值比远场站点更接近于理论计算值,前者可作为断层反演的约束条件。     

GPS测定的2011年日本9.0级地震的中国大陆地区同震位移场

利用陆态网络东北地区的6个GPS基准站地震期间的观测数据反演了日本本州东海岸9.0级强震区域地表的瞬时形变过程,应用160多个陆态网络基准站以及IGS跟踪站的GPS观测数据解算获得中国大陆地区的远场同震位移。将测站按照中国大陆主要构造带(区)进行分配统计可知,东北地区同震位移最大,平均接近20mm;首都圈与郯-庐带及周围地区平均同震位移也在10mm左右;距离震中比较远的鄂尔多斯及周围、阿尔金-祁连山-阿拉善、新疆等地区也受到了该地震的影响,在东方向发生了5mm左右的位移,并且以上地区的优势方向都为东方向;而本次地震对中国的华南(包括南海)以及青藏等地区基本没有影响。这些解算结果对进一步研究地震动力学特征及判定未来地震趋势提供了详实的实测依据。     

尼泊尔Mw 7.9级地震同震垂直位移与断层运动模型

2015年4月25日尼泊尔地区发生了Mw 7.9级地震,发震断层位于印度板块与欧亚板块碰撞边界带,此次地震是一次典型的板块逆冲型事件。利用中国境内加密的GPS同震观测资料,融合ALOS-2卫星L波段的InSAR（interferometric synthetic aperture radar）同震形变数据,基于最小二乘方法获得了此次地震的同震垂直位移场。同震垂直位移结果表明,此次地震造成尼泊尔加德满都地区抬升约0.95 m,珠穆朗玛峰地区受地震的影响有所下降,其主峰的沉降量为2~3 cm,中国境内的希夏邦马主峰沉降约为20 cm。地区利用改进的二维弹性半空间位错模型反演了发震断层运动参数,本文模型显示此次地震的断层面破裂宽度约为60 km,平均滑动量达到4 m,相当于Mw 7.89级。     

利用InSAR研究西藏拉孜地震同震形变

采用欧空局ERS—1卫星雷达数据获取了1993年Mw 6.2级西藏拉孜地震同震形变场,并采用弹性位错模型反演了该地震的断层参数。结果显示该地震产生的最大地表位移约为12 cm,断层倾角为42°,不支持该断层为低角度正断层的论断。     

日本2011 Tohoku-Oki Mw 9.0级地震的同震形变及其滑动分布反演:GPS和InSAR约束

利用覆盖整个日本东北部的3个条带的Envisat/ASAR降轨数据和6个条带的ALOS/PALSAR升轨数据,通过二通差分干涉处理,首先得到了2011年日本Tohoku-Oki Mw 9.0级地震的初步同震形变干涉结果,然后,利用GPS同震观测值对InSAR结果进行校正和基准统一。在此基础上,采用弹性半空间矩形位错模型,联合GPS和InSAR观测结果对同震滑动分布进行反演,获取了发震断层的断层滑动分布。结果显示,此次地震的滑动分布主要发生在40～50km深度范围内,最大滑移量为50.3m,释放的能量为3.20×1022 N·m(相当于Mw 8.94级)。     

GNSS观测的2021年青海玛多地震（Mw 7.4）同震形变及其滑动分布

北京时间2020年5月22日2时,中国青海省果洛州玛多县发生Mw 7.4地震。收集震中附近9个连续运行卫星定位基准站(continuously operating reference stations,CORS)观测数据,基于模糊度解算的精密单点定位(precise point positioning with ambiguity resolution,PPP-AR)技术处理了2 h时段的高频(1 Hz采样)全球导航卫星系统(global navigation satellite system,GNSS)数据,快速确定本次地震的近场同震位移,其中水平方向最大约为0.6 m。结合远场12个CORS震前、震后各3 d低频(30 s采样)数据的非差PPP解算同震形变场,基于弹性位错模型反演了玛多地震断层几何参数和滑动分布。结果显示:玛多地震属于典型走滑事件,发震断层走向278.49°,倾角为64.38°,滑动角为-10.90°,破裂长度约为138.72 km,宽度为4.82 km;滑动量超过3 m的滑移主要集中在东部小于18 km深度的区域,最大破裂可达4.2 m。反演的地震矩为1.85×1020Nm,相当于矩震级7.45,比美国地质调查局利用地震波的反演结果略大。     

利用GPS资料反演汶川Mw7.9级地震滑动分布

2008年5月12日，青藏高原东缘龙门山断裂带发生Mw 7.9级汶川地震，该地震使得北川-映秀断裂、灌县-江油断裂和汉旺-白鹿断裂带发生了同震破裂。本文主要利用“中国地壳运动观测网络”项目组提供的GPS同震位移资料，采用基于敏感度迭代拟合法(SBIF)及分段发震断层模型，分别应用均匀半空间和分层地壳结构两种方法反演分析了该次地震的精细同震滑动分布。结果表明：（1）SBIF反演方法应用于相对较小的GPS数据集仍然可以获得理想的同震滑动分布;（2）采用考虑分层地壳结构模型的反演结果整体上要优于均匀半空间模型的结果；（3）同震滑动主要发生在10-19千米深度以上，且存在五个滑动峰值区域及北川断裂南西端存在明显的深部滑动区，与震后野外考察结果及余震分布相吻合，同时能够较好地解释GPS观测数据；（4）青川断裂存在明显的右旋走滑分量，其平均滑动量为1.99米，北川断裂主要发生逆冲滑动，平均滑动量为3.35米，灌县断裂的平均滑动量为0.65米；（5）反演得到的地震矩为8.74×1020Nm(Mw 7.90)，与地震学结果吻合。     

基于GPS观测数据的汶川地震断层形变反演分析

通过高精度数据处理方法获取了2008年5月12日汶川地震(Mw=8.O)可靠的同震形变位移,较大形变主要在汶川、青川和北川县,其中最大同震形变不在震源处,而在北川县,水平方向达到2.3 m,垂向为0.7 m.并进一步采用分布滑动模型,利用弹性半空间均匀位错理论反演得到同震断层滑动分布.研究结果表明,不约束断层滑动特性的反演结果显示汶川地震断层不仅仅是右旋、逆断层,还伴随着少许左旋走滑及正断层变化.其地震矩张量为2.38X1021N·M,基于反演的断层滑动分布正演的区域同震形变场主要变形区域分别在都江堰-汶川县,北川县,青川县等一带,地表最大变形达到东方向4.5m,北方向0.8m,具有极强的右旋走势,均与野外地质调查结果相符合.     

GPS解算的日本Mw9.0级地震的远场同震地表位移

利用日本Mw9.0地震前后各3d的77个GPS连续站观测资料和JPL精密星历,计算了该地震对周缘各国及相邻板块运动的同震影响。结果显示,日本地震对日本、中国东北部、朝鲜半岛均有较大影响,其同震水平位移影响范围达2 000km,地震造成中国东北东移1～3cm,华北东移3～8mm,朝鲜半岛东移2cm。中国境内的GPS测站的垂直同震位移均呈上升趋势,但是上升幅度非常小。结合日本岛测定的GPS同震位移和约20m高的海啸可以推断,本次日本地震的同震垂直位移具有波状起伏特征,符合逆冲挤压型地震的粘弹性理论模型。鉴于近年地震活动趋势,应该密切关注横贯中国东部的深大断裂带（郯庐断裂带）的断层活动,借助中国陆态网络密集覆盖提供的有利条件,跟踪研究郯庐断裂带附近300km范围的区域地壳现今变形、应力应变场的时变过程等,分析郯庐断裂带的断层活动,捕捉异常信息。     

基于GPS观测的Baja California地震地壳变形分析

2010-04-04墨西哥Baja California(32.259°N,115.287°W)Ms 7.2级地震发生在GPS连续观测站及钻孔应变仪相当密集的区域,利用GPS监测数据计算出了若干站点地震前后的位移量。地壳变形分析结果表明,离震中近的观测站变形量普遍较大,最大可达0.2m,大部分观测站明显向东南方向运动。通过选择较稳定的观测站作为参考站,计算了4观测站的动态变形序列,同时得到地震波的传播速度为3.2km/s,与应变观测计算结果比较吻合。     

GPS时间序列中同震和震后形变的自动识别和估计

地震发生时,附近GPS观测站会捕捉到同震和震后形变,同震和震后形变对地学研究和维持动态参考框架具有重要作用,识别和估计同震和震后形变也是GPS时间序列处理的重要内容。提出综合检校法识别和估计GPS时间序列中的同震和震后形变,首先利用该方法对GPS站建模,综合考虑多源的同震阶跃及兼容性和均方根RMS减小率来判断测站是否受到地震影响,然后用试错法搜索最优衰减常数τ后,利用最小二乘估计包含同震和震后形变的所有参数,最后以真实震例验证了该方法的有效性。     

台湾集集大地震断层非均匀滑动分布的反演

利用台湾1999集集大地震前后GPS测量获得的地表位移数据,反演断层面上的非均匀滑动分布.反演中采用二次样条函数作为滑动分布的基函数,通过引入阿卡克贝叶斯准则,获得了反演问题的惟一解.反演得到的地表位移与实测位移的差在几个厘米,较好地解释了观测到的GPS位移.另一方面,用反演结果计算得到的地表形变场同INSAR得到的地表实测位移场较相近.     

Deriving 3D coseismic deformation field by combining GPS and InSAR data based on the elastic dislocation model

The density of GPS measurements is usually one of the key issues in resolving 3-D coseismic deformation field from integrating GPS and Interferometric Synthetic Aperture Radar (InSAR) measurements with pure mathematic interpolation methods An approach that combines the elastic dislocation model with the Best Quadratic Unbiased Estimator (BQUE) or a robust estimation method named IGG (Institute of Geodesy and Geophysics) is proposed to reconstruct 3-D coseismic deformation field, in which only a small amount of GPS data is needed to produce a reasonable initial 3-D coseismic deformation. Then the BQUE and IGG are used to weight the InSAR and GPS measurements to avoid computational issues caused by the negative variance problem and to decrease the impact from gross errors. The Wenchuan earthquake is used to test the proposed method. We find that the developed method makes it possible to use only a few GPSs and InSAR data to recover the 3-D coseismic deformation field, which offers extensive future usage for measuring earthquake deformation, particularly in some tectonically active regions with sparse GPS measurements.     

联合地震位错模型和InSAR数据构建2017年九寨沟Mw 6.5地震同震三维形变场

利用哨兵（Sentinel）-1A卫星升、降轨影像,在地震位错模型约束下获取了2017年九寨沟Mw 6.5地震的高质量三维形变场。首先,利用合成孔径雷达干涉测量技术（interferometric synthetic aperture radar,InSAR）提取九寨沟地震升、降轨同震形变场;然后,通过“两步法”反演获取该地震发震断层的几何参数和分布式滑动模型,以此为约束,采用方差分量估计算法联合解算九寨沟地震三维形变场。结果表明,九寨沟地震同震三维形变场以水平位移为主,垂向形变较弱;南北向形变呈拉张趋势,断层上盘向南、下盘向北滑动,最大位移分别为－19.81 cm和14.38 cm;东西向形变不对称性明显,断层上盘西北部向东水平运动,最大位移为18.37 cm,下盘东南部向西运动,最大位移不足8 cm。将南北、东西向形变与6个全球导航卫星系统（global navigation satellite system,GNSS）台站观测数据进行比较,两者一致性较好且均方根误差较小,分别为1.44 cm和1.77 cm,表明联合升、降轨InSAR观测和地震位错模型约束构建同震三维形变场方法具有较高可行性,显著降低了大地测量数据不足、InSAR观测对南北向形变不敏感等问题的影响。     

唐山7.8级地震三维地表形变和断层运动特征

通过系统收集整理1976年唐山7.8级地震前后区域水准和边角网观测资料,在研究发展动态平差程序的基础上,解算得到了唐山地震震前和同震三维形变定量结果与误差分布,并反演得到了断层面同震滑动分布。结果显示,震前阶段丰润-唐山-唐海-昌黎-迁安一带整体隆升,隆升高值区位于唐山附近,震中区的水平形变量不大,唐山菱形地块相对于周围地块逆时针扭转;同震阶段唐山断裂以右旋错动为主,并呈现西升东降的垂向变形特征;反演结果表明多段断层参与了唐山7.8级地震的应变释放,主震发震断层的南段破裂主要发生在深度约6~18 km、长轴约50 km范围,北段破裂主要发生在深度约7~17 km、长轴约30 km范围,滦县地震断层也发生了少量滑动;主震断层最大水平位错为7.82 m、最大垂直位错为2.04 m,等效震级为Mw 7.58。挖掘得到的唐山震前和同震三维形变场定量结果为深入研究该地震的孕育发生机制及地震影响提供了珍贵资料。     

2008年汶川地震前的形变异常及机理解释

2008年5月12日发生了汶川8.0级大地震,震前的定点形变观测值曲线异常很少。震后通过对定点形变观测台站的震前数据进行了深入分析和处理,发现在汶川8.0级地震前存在形变异常,且距震中较近的台站,出现的异常信息的频率较高;距震中稍远的台站,出现的异常信息的频率较低。对这些现象进行了初步的解释。