不同卫星导航系统组合监测的玛多地震同震形变与精度分析

为全面分析不同卫星导航系统(GNSS)组合监测地震形变的精度,本文首先基于青海玛多地震(Mw 7.4)周边近场及远场共8个区域CORS站的高频(采样率为1 Hz)GNSS观测数据,采用精密单点定位技术进行高频数据解算,初步分析了多系统组合相对单GPS系统的同震形变精度提升,然后对具有相同卫星系统数/不同卫星系统组合的同震形变误差进行对比分析,最后基于GREC2(GPS/GLONASS/Galileo/BD-2)和GREC3(GPS/GLONASS/Galileo/BD-3)两种组合情况分析了北斗二号与北斗三号的同震形变精度差异。研究结果表明,多系统组合观测同震形变的精度优于单系统,其中GEC2+3(GPS/Galileo/BD-2/BD-3)组合观测同震形变精度最优,相对于单系统GPS而言,在水平方向上的精度提升可达30%～60%,在垂向上的精度平均提升可达30%;在三系统卫星导航系统组合中,含北斗的多系统组合在水平方向上的精度提升最高可达38.75%,垂向上的精度提升最高可达41.46%;北斗三号观测的同震形变精度比北斗二号精度更高、效果更好。     

哨兵-2号光学影像地表形变监测:以2016年Mw7.8新西兰凯库拉地震为例

光学影像已经广泛地应用于地表形变监测研究。哨兵-2号光学影像作为一种新的对地观测数据,具有重访周期短、空间分辨率高、影像覆盖范围大以及数据免费等优点,因此该数据在地表形变监测上有广泛的应用潜力。本文以COSI-Corr软件包为数据处理平台,基于亚像素的频率域相关性匹配技术,处理多时相的哨兵-2号数据获取地表形变。本文选取2016年M_w7.8新西兰凯库拉（Kaikoura）地震覆盖区域为例,对哨兵-2号影像地表形变中存在的各种系统误差源进行了系统分析和改正,并提出了改进的均值相减法去除形变场中的卫星姿态角误差。另外,还对哨兵-2号4个10 m空间分辨率的波段（Band 2/3/4/8）中可用于地表形变监测的最佳波段进行了分析,统计结果显示Band8的地表形变监测效果最好。最后,利用哨兵-2号光学影像获取了2016年11月14日M_w7.8新西兰凯库拉地震的同震形变场;分析了沿地震主要断层的滑移分布,结果表明最大水平滑移量达10 m;并与同期Landsat8全色影像的同震形变监测结果进行对比分析,结果表明哨兵-2号结果精度更高。本文的研究成果可以为哨兵-2号光学影像的应用提供参考。     

基于GNSS的青海玛多M7.4级地震的地表形变分析

2021年5月22日,青海省玛多县发生M7.4级地震,地震对区域基准站造成了影响. 本文收集了距震中400 km以内的青海省全球卫星导航系统(GNSS)连续运行参考站(CORS)的高频观测数据,以国际GNSS服务 (IGS)站和陆态网为参考基准,分别采用静态模式和动态模式分析了震后地表三维形变. 结果表明:站点位移量级随震中距的加大而呈现衰减趋势,其中距震中约35 km的江多站(JDUO)产生东偏南28.0 cm的永久性位移. 整体来讲,地震对100 km内的震源区的影响表现为走滑型运动,与震源机制相吻合,区域运动特征整体表现为拉张运动,北西-南东(NW-SE)向拉张和北东-南西(NE-SW)向挤压.      

青海玉树地震差分干涉雷达同震形变测量

2010-04-14青海玉树发生7.1级地震后,作者利用震前和震后获取的日本ALOS卫星PALSAR遥感数据,开展了差分干涉雷达(D-InSAR)地震同震形变测量与分析。结果表明:玉树地震引起较大范围地表变形,地震变形沿玉树—甘孜断裂带向南东东方向扩展,在N33.7°,E96.81°附近达到最大形变量,D-InSAR监测到雷达视向上的最大形变量为35cm。地表形变特征对于评价玉树地震破坏程度、推断断层性质、研究地震形变和地震孕育特征具有重要的参考价值。     

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,比美国地质调查局利用地震波的反演结果略大。     

2021年青海玛多Mw 7.4地震前后b值的变化特征研究

古登堡-里克特频度-震级关系式中的b值常被用作估算区域应力大小的一个重要指标,对于评估地震发生概率起着至关重要的作用。北京时间2021-05-22T02:04:11青海省玛多县附近发生Mw 7.4地震,此次地震发生在巴颜喀拉块体内部,东昆仑断裂以南的次级断裂昆仑山口-江错断裂上。为了揭示这次玛多地震之前震源区的孕震环境和震后余震的演化,利用中国地震台网中心提供的2009-01-01-2021-05-20地震目录,对玛多地震发震前后的b值分布和变化进行了研究。研究发现,在2021年5月22日之前,玛多地震主震区b值最低,表明震源区存在明显的应力积累;从近十年震源区b值的时间序列来看,b值也呈现下降趋势,该发现可以作为今后判定块体内部大地震主震临近的重要依据。利用震后的余震序列计算了断层周边的b值变化,发现早期的b值变化对后期余震的发生位置预测也具有一定的指导作用。     

DInSAR技术监测玉树地震同震形变场的研究

近二十年来,合成孔径雷达干涉测量技术（DInSAR）作为一种监测地表形变的有效技术得到广泛应用,其视线向精度可达厘米级甚至毫米级。该技术尤其是对监测快速、激烈的地表形变更为有效,如地震、火山等。利用DInSAR技术,选用日本ALOS卫星PALSAR数据,获取了2010—04—14玉树Mw6．9地震的同震形变场。结果表明：地表至少发生过3次地表破裂,沿断层走向的形变分布范围远远大于垂直断层分布范围,形变分布特征以左旋走滑为主。     

2022年泸定Mw 6.6地震InSAR同震形变与滑动分布

2022-09-05,青藏高原东缘的鲜水河断裂上发生了泸定Mw 6.6地震,该地震是鲜水河断裂上40年来发生的最大地震,研究该地震的运动学和同震破裂模式对理解青藏高原东缘构造形变机制和评估鲜水河断裂以及安宁河断裂的地震危险性具有重要意义。利用Sentinel-1和ALOS-2卫星雷达影像,采用合成孔径雷达干涉技术获取了泸定地震的同震形变场,进而基于弹性半空间的位错模型,确定了本次地震发震断层的几何参数和滑动分布。结果表明,泸定地震是一次典型的左旋走滑事件,发震断层西倾,倾角约为72°,走向沿NNW-SSE方向,约为167°;断层破裂主要集中在0~10 km深度,最大滑动发生在约5.8 km深度,约为2.23 m;同震释放的地震矩约为8.74×1018 N·m,相当于矩震级Mw 6.59。通过对震后光学影像解译,发现此次地震诱发的滑坡多集中分布在发震断层西侧,该现象与余震主要集中在断层西侧的结果相一致,可认为是地震上盘效应的体现。     

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

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

芦山Ms 7.0级地震震前及同震地表形变

对四川地区GNSS观测数据进行处理,获得如下结果:①汶川Ms8.0级地震之前龙门山断裂带及其以北近150km的范围内几乎不存在可分辨的右旋走滑活动和压性形变,四川盆地以西-鲜水河安宁河断裂带以东约100km的范围内几乎也不存在可辨的形变;②汶川Ms8.0级地震之后芦山Ms7.0级地震之前,龙门山断裂带中、北段存在约5mm/a右旋走滑活动,芦山震源及附近地区的三维形变一直处在闭锁状态,鲜水河安宁河断裂带及其以东的左旋形变有所增强;③芦山地震可变的同震形变场基本分布于以震源为中心的数十km范围内,大致以震中为界东侧呈右旋形变,西侧呈左旋形变;④距震中较近约12km的LS05震时\     

2008年大柴旦M_w6.3级地震的InSAR同震形变观测及断层参数反演

2008年11月10日,青海省大柴旦地区发生了Mw6.3级地震。本文利用EnviSat卫星升降轨SAR数据和差分干涉测量技术提取了同震形变场,基于均匀位错模型反演确定了地震断层参数,然后利用格网迭代搜索法确定了较优断层倾角,同时基于非均匀位错模型获得了精细滑动分布。结果表明,地震使得上盘区域沿降轨、升轨视线向分别产生最大约8.5cm、10cm的抬升;较优断层倾角为47.9°;地震滑动未延伸至地表,主要发生在地下8.2~23.7km深度范围内,最大和平均滑动量为0.5m和0.19m,平均滑动角为104.9°。反演的地震矩为3.74×1018 N·m,矩震级为Mw6.35。     

日本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级)。     

利用MSMPCA去噪的CEEMD方法监测高频GNSS同震形变

考虑到全球导航卫星系统(global navigation satellite system,GNSS)地震信号的非线性和非平稳性,利用一种多尺度多方向主成分分析(multiscale multiway principal component analysis,MSMPCA)去噪的完备总体经验模态分解(complete...     

2015年尼泊尔地震同震滑动及震后余滑的三角位错模型反演

2015年尼泊尔Mw 7.9级地震发生在印度板块向欧亚板块低角度俯冲的喜马拉雅断裂带上。对该地震的滑动模型进行精化,对于理解喜马拉雅断裂带的孕震模式具有重要意义。采用三角形位错元构建主喜马拉雅断裂“双断坡”几何模型,联合全球定位系统（Global Positioning System,GPS）和合成孔径雷达干涉测量（interferometric synthetic aperture radar,InSAR）资料反演2015年尼泊尔地震同震滑移及震后余滑。结果表明,尼泊尔地震最大同震滑移达到7.8 m,深度为15 km,位于中地壳断坡和浅层断坪的接触部位。不考虑中地壳断坡结构会使反演的最大滑移量偏低。震后余滑主要分布在同震破裂区北侧,释放的地震矩为1.02×1020 N·m,相当于一次Mw 7.3级地震,约占主震释放地震矩的12%。同震库伦应力变化和震间断层闭锁分布均表明,尼泊尔地震破裂区南部宽约60 km的区域仍具有较高的地震危险性。     

基于GIS的InSAR结果分析方法及在汶川Mw7.9级地震同震解释中的应用

利用GIS对InSAR结果进行分析可以更全面地解释地表形变机制。介绍了利用GIS对SAR影像进行裁剪、切割和拼接处理、剖面和等值分析的方法。在此基础上,以汶川Mw7.9级地震为例介绍了基于GIS的InSAR结果分析方法在地震形变机制解释方面的应用:①根据InSAR距离向偏移量和同震形变场,利用GIS技术提取了汶川Mw7.9级地震发震断层的地表迹线;②将基于位错模型正演得到的视线向形变场与InSAR同震形变场进行比较,得到了汶川地区InSAR同震形变场的上盘和整体改正数,进一步借助GIS进行修正得到了更为真实的InSAR同震形变场;③利用2D/3DInSAR同震形变场、干涉纹图和剖面分析结果对汶川Mw7.9级地震同震形变特征进行了具体分析和解释。     

尼泊尔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级。     

震后松弛过程的粘弹性模型在1997年Mw7.6玛尼地震中的应用研究

采用麦克斯韦体模型、标准线性体模型和伯格斯体模型,具体研究了1997年玛尼地震在衰减时间1~100 a间的震后位移场和库仑应力场。结果表明:①在1~5 a的短期内,采用标准线性体和伯格斯体模型的震后位移场具有时空上的一致性;而在10 a以上的中长期内,麦克斯韦单元开始起主要作用,表现为采用麦克斯韦体和伯格斯体模型的震后位移场具有时空上的一致性。②虽然应力阴影区和应力触发区随时间而同步扩展,但短期内,采用标准线性体和伯格斯体模型的粘弹性库仑应力场的空间结构具有一致性;而中长期内,采用麦克斯韦体和伯格斯体模型的粘弹性库仑应力场的空间结构具有一致性。为更好地解释震后形变及粘弹性库仑应力触发响应,流变模型的选择需要根据研究的时间段和观测数据源具体确定。     

东日本Mw9.0大地震前、同震及震后地壳水平运动

利用“中国地壳运动观测网络(二期)”245个GPS连续观测站及其他站的大量观测结果,采用离震中很远的7个GPS站作为位错参考框架,得到2011年3月11日东日本Mw9.0大地震的同震水平位移和应变。此种参考框架的位移场解,符合弹性位错模型离震中很远且水平位移接近于0的要求,不仅有利于位错模型反演,也有利于合理分析此次大地震同震水平位移的影响范围。地震前后区域参考框架下远场GPS连续观测位移时间序列和同震位移弹性位错模型反演结果表明,东日本大地震同震水平位移既是地震断层破裂的结果,也是远场震前数月或数年水平位移的弹性回跳,据此可以进一步确认这些站震前水平位移前兆异常。由于这些震前异常位移中亦包含了同震位移量(方向相反),因此,弹性位错模型为定量研究大地震前远场地壳运动异常提供了一种理论模型。位错参考框架中的同震位移及其位错模型和区域参考框架的位移时间系列的综合研究有助于对此次大地震地壳运动异常的认识。这种持续时间为数月至数年的弹性形变异常,可以为大地震的中短期预测提供依据。