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

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

2013年芦山Ms 7.0级地震断层参数模型反演

2013年4月20日四川芦山地区发生了Ms 7.0级地震。利用GPS三维同震形变数据获取地表形变场，基于位错模型反演芦山地震的断层几何参数及滑动分布。首先采用多峰值颗粒群算法（multiple peak particle swarm optimization，MPSO）得到断层几何参数，其中断层走向206.47°，倾角44.11°，长度21.94 km，离地表最浅处为7.66 km，最深处为17.84 km。为了反演断层面的精细滑动分布，分析地震所在的龙门山断裂西南段破裂面具有的铲状型特征，将芦山地震破裂面确立为铲状模型，即将断层的倾角预设为上陡下缓，倾角变化范围为21°~50°。结果显示，断层破裂面在不同深度区域出现了两个滑动峰值，其中最大滑动量为0.68 m，深度位于13 km。地震释放的能量为1.47×1019 N·m，对应的矩震级为Mw 6.74，与地震学的研究结果一致。     

利用GNSS基线分析芦山Ms7.0级地震前后应变演变特征

利用中国地壳运动观测网络6个GNSS基准站2011~2014年的观测数据,解算了芦山Ms7.0级地震近场区域的GNSS基线长度变化时间序列。从投影后的平面基线出发,给出了由基线长度变化求解面应变的方法。利用基线线应变、基线所得的最大剪应变、最大面应变和第一、第二剪应变等5个物理量分析了地震前后地壳运动的动态演变过程。线应变序列结果显示,整体地壳形变表现为NW-SE向挤压,S-N向拉张状态,震前线应变异常并不明显,震后经过调整又恢复为原线性变化趋势,总的线应变仍在持续积累。面应变序列结果显示,震前出现了不同程度的非正常偏离,尤其出现了明显的闭锁状态和反向加速过程,整个变化曲线类似抛物线的弧形,与岩石力学形变理论相吻合。由第一剪应变的显著异常变化,推断发震断层南缘在这次地震活动中产生了强大的左旋剪切构造力。     

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。通过对震后光学影像解译,发现此次地震诱发的滑坡多集中分布在发震断层西侧,该现象与余震主要集中在断层西侧的结果相一致,可认为是地震上盘效应的体现。     

基于多源SAR数据的2022年门源Ms 6.9地震同震破裂模型反演研究

2022年1月8日青海省海北州门源县发生Ms 6.9地震,震中位于青藏高原东北缘祁连-海原断裂中段,属历史地震空区,基于多源合成孔径雷达（synthetic aperture radar,SAR）遥感数据研究该地震的破裂模式对理解青藏高原东北缘构造变形机制、应变释放过程以及地震危险性评估具有重要意义。首先利用Sentinel-1数据和合成孔径雷达差分干涉测量（differential interferometry synthetic aperture radar,D-InSAR）技术获取了门源地震的同震形变场,视线（line of sight,LOS）向形变场显示此次地震造成了约20 km长的地表破裂,最大形变约0.75 m;然后基于Sentinel-2卫星数据,利用光学影像配准和相关技术获取了本次地震的东西向同震形变场,最大同震位移达2.5 m;最后基于均匀弹性半无限位错模型,以LOS向形变场为约束反演了断层的滑动分布模型。结果显示,门源地震是一次典型的左旋走滑型地震,地震破裂主要集中在0~10 km深度范围,最大滑动量3.25 m,滑动角10.44°,对应深度4.89 km;反演给出的矩震量为1.07×1019 N·m,对应矩震级Mw 6.6。结合野外考察和地质资料,初步判定发震断裂为冷龙岭断裂,并引起托莱山断裂发生同震滑动。同震库仑应力结果显示,冷龙岭断裂东段和托莱山断裂西段应力状态为加载,未来具有发生强震的风险。     

九寨沟Ms7.0地震的同震形变场误差研究及断层滑动分布反演

北京时间2017年8月8日四川省九寨沟发生Ms7.0级地震。为了提高断层滑动分布反演的可靠性,本文利用升降轨InSAR数据联合反演断层滑动分布,再根据模拟形变值二次反演,对比反演结果来探讨升降轨形变场误差和确定断层滑动分布的结果。结果表明,九寨沟升轨同震形变场的质量较好,降轨形变场受余震、震后粘弹性松弛效应影响明显,利用模拟形变值二次反演确定的断层滑动分布可靠性更高。     

跨断层形变的慢时变结构分析与异常识别初探

将自适应控制理论引入跨断层流动形变分析。利用甘肃及与宁夏、青海交界的监测区内石灰窑口、红柳峡、窝子滩、扁都口、南营5个场地典型观测曲线,构建慢时变前兆系统结构,研究其与6级左右及6级以上地震孕育-发生过程的可能关系。可显示震前中短期阶段最大特征根模1.0以上异常,异常识别较原始观测曲线分析减少了人为性与经验性,也反映出断层形变前兆系统结构的动态变化和自适应性。     

利用遥感技术对鲁甸6.5级地震发震断层的判定

2014-08-03鲁甸发生6.5级地震,由于其上有北西方向的包谷垴-小河断裂和北东方向的昭通-鲁甸断裂,且震中无明显地表破裂带,故而对发震断层存在一定争议。由往年地震数据分析可知,滑坡分布与发震带在方向上有着一致性,所以选取同一季节时间相近的地震前后的遥感影像,排除季节对植被的影响,分析地震前后的影像中植被覆盖度,间接确定滑坡分布方向。研究表明,植被覆盖度变化主要集中在北西方向,可以判定鲁甸地震的发震断层是北西方向的包谷垴-小河断裂。     

芦山地震前后川滇地区地壳水平形变特征研究

为获取芦山地震前后川滇地区地壳的形变特征,利用陆态网2009—2013年和2014—2016年两期全球卫星导航系统(GNSS)水平速度场资料分别计算并对比分析了地震发生前后主应变率场、最大剪应变率场、面膨胀率场及基线长度的变化情况。结果显示芦山地震之前龙门山断裂带主要以压缩应变为主,面压缩应变和最大剪应变均处于高值区,震后能量部分释放,压缩形变程度减弱,但其西南方向的安宁河断裂带和鲜水河断裂带南段出现明显的压缩应变高值区,同时南汀河断裂带附近也呈现较明显的压缩应变。岷江断裂带附近区域的压缩应变虽然减弱,但仍然没有改变它的应变状态,这可能促使了2017年九寨沟Ms7.0地震的发生。     

GNSS网形变化的区域地壳形变信息提取方法

针对GNSS三维无约束平差不依赖基准且拥有高精度几何网形的特点,以GNSS监测网作为一个整体的时空观测单元,由GNSS网测站间相互关联的所有基线长度和基线夹角变化综合衡量网形的变化,用网形的变化集中反映区域地壳形变信息。利用地壳形变高空间相关性的特点,采用主成分时空响应分析的方法实现了区域地壳形变信息的快速提取。通过对2013年芦山Ms7.0地震近场区域陆态网GPS基准站2010—2014年观测数据的解算,提取并分析了对应地震前后近场区域地壳形变时空分布特征及形变的动态演变过程。     

青藏高原中北部地貌变形的遥感表征及其生态地质学意义

遥感地质特征显示青藏高原北部,在澜沧江—温泉与东昆仑南缘—阿尼玛卿两条岩石圈断裂之间区域存在着4个正交叠加在早期北西—南东向构造地貌上的北东—南西向第四纪坳陷。其中西羌塘北部坳陷和长江源坳陷的地堑式裂陷地貌特征明显,而黄河源坳陷和若尔盖坳陷发育程度较低,但它们均是在青藏高原大陆构造变形的统一动力学背景下,地貌发生东西向伸展变形的产物。进一步研究表明,它们的形成不仅与青藏高原中南部所发育的近南北向裂谷具有同样的动力学意义,而且已成为影响该地区生态系统的主要因素。     

Research on the deformations in Qinghai-Tibet plateau and its margins by inverting seismic moment tensors and GPS velocities

We have determined approximate average rates of deformation in the Qinghai-Tibet plateau and its margins from the GPS data for last 10 years and the moment tensors from earthquakes between 1900 and 1999. We also determined the strain rate (seismic strain rate) associated with the seismic deformation using 254Mω≥5.0 earthquakes, and estimated the shortening and extension rates for every block in the area as well. We also estimated the strain rate (geodetic strain rate) by 80 GPS sites’ velocity vectors and analyzed characteristic of kinematics by two kinds of strain rates and discussed earthquake potential in the area. As a result, the deformation rates from seismic moment tensors and from GPS velocities are basically agreed with each other. It is feasible to analyze seismic risk by comparing geodetic strain rate with seismic strain rate based on the opinion that strain energy will be released through earthquake. It is concluded that there is no strong earthquake potential (>M7) in the Qinghai-Tibet plateau and its margins, but there is earthquake potential (>M5) in middle Tibet in a few years.     

2013年Ms 7.0级中国芦山地震断层曲面模型的构建及其滑动分布的大地测量反演

根据中国芦山地震的余震序列分布和地质构造情况,联合GPS同震位移和水准测量数据,采用格网搜索法和方差分量估计（variance covariance component,VCE）方法建立了芦山地震的断层曲面模型。以建立的断层曲面模型为基础进行滑动分布反演,反演结果表明此次地震以明显的逆冲滑动为主,存在左旋趋势,断层的滑动主要发生在4.56~18.25 km深度范围内,最大滑动量为0.70 m,在深度12.64 km附近,释放的能量为8.77×1018 N/m（相当于Mw 6.57级）。对曲面模型进行棋盘检验,结果表明其能够较好地分辨出断层深度浅于13 km的凹凸体,其面积最小为8 km×8 km。     

2014年鲁甸6.5级地震震前TEC异常分析

2014年8月3日云南省鲁甸发生了6.5级地震.针对震前电离层异常问题,利用IGS提供的全球电离层数据,采用滑动四分位法对鲁甸地震震中及其周围格网点的TEC进行异常探测,发现震前3d、6d和13 d电离层有明显的异常,其中震前3d、6d的异常最为显著.在进一步分析了太阳地磁活动及这几天的全球TEC异常分布后,认为这3d的电离层异常可能与震前的孕育地震有关.另外,这3d的TEC异常分布区域都满足在经度方向跨度和纬度方向跨度长度比约为3:1,且异常的峰值点并不位于震中的正上方,而是位于震中偏向赤道方向,同时赤道共轭区域也存在异常,但范围和幅度都较小.     

Research on the deformations in Qinghai-Tibet plateau and its margins by inverting seismic moment tensors and GPS velocities

We have determined approximate average rates of deformation in the Qinghai-Tibet plateau and its margins from the GPS data for last 10 years and the moment tensors from earthquakes between 1900 and 1999.We also determined the strain rate (seismic strain rate) associated with the seismic deformation using 254 Mw≥5.0 earthquakes,and estimated the shortening and extension rates for every block in the area as well.We also estimated the strain rate (geodetic strain rate)by 80 GPS sites' velocity vectors and analyzed characteristic of kinematics by two kinds of strain rates and discussed earthquake potential in the area.As a result,the deformation rates from seismic moment tensors and from GPS velocities are basically agreed with each other.It is feasible to analyze seismic risk by comparing geodetic strain rate with seismic strain rate based on the opinion that strain energy will be released through earthquake.It is concluded that there is no strong earthquake potential (＞M7) in the Qinghai-Tibet plateau and its margins,but there is earthquake potential (＞M5) in middle Tibet in a few years.