利用GPS 位移反演日本Mw 9.0 仙台地震及Mw 7.9 强余震静态位错模型

基于GPS 观测的M_w 9.0 日本仙台地震同震形变及震后8 h 的地表形变, 考虑分层介质模型的影响, 反演了主震及震后8 h 的断层的位错分布. 结果显示, 主震位错对应的矩震级为M_w 8.98, 最大同震位错量可达23.3 m, 在震中区显示出纯逆冲特性, 震中两侧断层错动显示出一定的走滑分量. 另外, 约90%的地震能量发生在40 km 深度以上. 而由震后8 h 的地表形变反演的结果显示, 这一期间释放能量约相当于一个M_w 8.13 的地震, 位错主要分布在主震破裂的西南区域, 最大位错量约1.5 m, 破裂峰值区与M_w 7.9 强余震有明显的对应关系, 暗示震后8 h的断层错动主要是由M_w 7.9 强余震引起的. 另外, 主震破裂的深部延伸面在震后0.2~0.4 m 的滑动可能主要是无震余滑.     

四川九寨沟M_s7.0地震主震及其余震序列精定位

采用基于三维速度模型的定位方法,测定了九寨沟M_s7.0级地震的主震位置,震中为103.806°E,33.201°N,震源深度为20.4 km.使用最近一个强震台(51JZB)的S-P到时差,估算了主震初始破裂点的深度不浅于14.3 km.利用流动观测获得定位精度较高的余震,对早期余震位置进行校正,并采用双差方法对震后一个月的余震进行了重定位,获得了3030个地震的位置.余震呈北西西向的条带状分布,长约42 km,余震北临塔藏断裂,南接虎牙断裂北端.主震位于余震带的中央,其两侧各有长约20 km的余震带,西北侧有一长约5 km的余震稀疏段.西北段余震深度较浅,余震带宽度约6 km;东南段余震深度较深,余震带较窄,约4 km.余震震源深度的优势分布范围在4~20 km之间.发震断层倾角较陡,平均值约为84°.断层倾向和倾角沿走向方向有明显变化,断层在浅部向西南倾斜,深部略向北东倾斜.主震初始破裂点深度大于矩心深度和余震平均深度,地震破裂由深向浅传播.余震沿断层走向有明显的时空扩展特征,震源区可能存在余滑.     

GPS观测的2015年尼泊尔M_S 8.1级地震震前应变积累及同震变形特征

2015年4月25日,在印度板块与欧亚板块交界区的喜马拉雅地震带上发生了尼泊尔MS 8.1级大地震.震前GPS速度场和应变率场显示,喜马拉雅地震带整体表现为15.94±1.82 mm/a的压缩特征,同时还具有分段活动特征.此次地震发生在速度场顺时针旋转和逆时针旋转的分界带,该处最大主压应变率的量值在喜马拉雅地震带并非最大.GPS观测的同震位移场揭示了尼泊尔MS 8.1级地震引起的地壳变形特征,分别有9和6个测站观测到明显的水平向和垂向同震位移,其水平分量的运动方向整体表现为南向运动,位于震中东南侧的3个测站垂直分量表现为上升,其余测站为下降.中国境内距离震中最近的5个测站的垂向同震位移显示,此次地震造成珠穆朗玛峰的沉降量约为4mm.依据GPS观测到的同震位移场,利用非负最小二乘方法反演震源断层面上的滑动分布.反演结果表明最大滑动量为6.84 m,滑动量较大的区域分布在加德满都附近及其以北区域的下方,这可能是造成加德满都地区具有较大破坏的原因之一,该滑动分布模型能够很好地解释GPS观测到的同震位移.利用此滑动分布模型计算的地震矩为8.21×1020 N m,对应的矩震级为MW7.9.     

2015年智利8.3级地震构造背景及对俯冲带地震的启示

2015年发生在智利Coquimbo的M_w8.3级地震是一个典型的俯冲带地震,其大小、破裂滑移分布与利用震间期GPS观测得到的智利南部俯冲带闭锁模型一致.利用反投影方法得到的破裂过程初步解显示,这次地震的破裂由震中(约25 km深)向浅部传播,并且这次地震破裂的能量辐射与频率明显相关,与2010年发生在Coquimbo地震破裂南部约50 km的M_w8.8级地震的观测结果一致,表明这部分俯冲带的断层性质随深度存在着系统变化.这次2015年8.3级地震的震前及同震观测资料将对地震破裂动力学的数值模拟提供一个难得的模型约束信息.     

西藏米林M6.9级地震及其余震序列地震定位

利用布设在南迦巴瓦构造结及其周边地区流动地震台的观测数据,测定了2017年11月18日西藏米林M6.9级地震的主震位置,主震震中位于29.87°N,95.05°E,震源深度12 km.采用Geiger法和双差定位法对余震序列进行了地震定位.定位结果表明余震沿着主震的北西和南东向两侧扩展分布,其展布长度约40 km,宽度约10 km.米林地震发生在西兴拉断裂带最南侧的一条次级断裂上,震源深度剖面和主震震源机制研究结果揭示发震断层以高角度、北东倾向为主要特征,米林地震触发了靠近主震发震断层北东侧一条近平行断层的地震活动.米林地震表明,南迦巴瓦构造结顶部目前仍处于较高的构造挤压状态,南迦巴瓦变质体持续向北东方向推进,导致构造结内部块体缩短变形,其东西两侧的墨脱断裂和米林断裂存在未来发生大震的可能性.     

GPS测定的尼泊尔M_w7.9和M_w7.3级地震同震形变场

2015年尼泊尔M w7.9和M w7.3级地震致灾范围包括尼泊尔、印度北部、巴基斯坦、孟加拉和中国藏南地区,地震应变应力调整将对震区和邻区的地震活动产生不同程度的影响.尼泊尔境内的GPS连续观测数据、"中国大陆构造环境监测网络"GPS基准站和地震应急流动站观测数据计算结果揭示了尼泊尔地震的静态和动态同震形变场.对M w7.9级地震,尼泊尔境内近场静态同震GPS水平位移最大为1.89 m,距震中100~400 km的中国藏南地区观测到几毫米到几十厘米的静态水平位移,毗邻尼泊尔的聂拉木县最大形变54.0 cm.M w7.3级地震静态同震形变最大为2 cm,局限在震中附近200 km范围内.距震中约2000 km内的GPS基准站均记录到M w7.9大震明显的动态形变信号,高频GPS动态形变幅度与地震破裂方向有关,位于破裂方向上的测站动态形变幅度明显大于其他方向的测站.用弹性半空间位错模型正演模拟了震区和青藏高原南部格网点上的同震形变,并分析了地震应变影响,认为尼泊尔地震对中国藏南地区产生一定程度的拉张型应力变化,需要持续关注.     

汶川8.0级地震地表破裂白鹿镇段的变形特征*

汶川8.0级地震在彭县-灌县断裂上形成的地表破裂带主要位于彭州磁峰与安县睢水之间,全长约90 km,白鹿镇位于该断裂地表破裂带的南段。应用全站仪和GPS对地表破裂带的垂直位错量、水平位错量、构造缩短率进行精确测量,结果显示,沿彭县-灌县断裂,白鹿镇地表破裂的垂直位错量和水平位错量都是最大的,其中最大垂直位错量位于白鹿老街,为(2.7±0.2) m,最大水平位错量位于白鹿河堤岸,为(0.7±0.2) m。究其原因,可能是沿彭县-灌县断裂,白鹿地区煤层含量最为丰富,煤层厚度最大,因而该地区软弱岩层的屈服强度相对最低,从而导致该地区断层两侧岩体的滑动量最大。另外,通过白鹿中学附近垂直于断层走向的探槽,甄别出两条呈叠瓦状排列的分支断层,该探槽所揭示的地表破裂面倾角度数及倾角随深度的变化趋势与通过震源机制解推断的结论一致。     

2022年泸定MS6.8地震同震形变特征及周边强震危险性

2022年四川泸定6.8级地震发生于鲜水河断裂磨西段,为揭示此次地震的同震形变特征,收集了震中200 km范围内的连续全球定位系统(global positioning system, GPS)测站和震中50 km范围内的强震动数据,进行了高精度处理以提取测站的同震水平位移.此外,还收集了覆盖震中区域的Sentinel和ALOS-2降轨数据,通过干涉差分技术处理并获得了卫星视线向同震形变场.结果显示:水平同震形变呈四象限分布,表明泸定地震的震源机制为左旋走滑;距震中16 km的强震动台站记录到了明显的位移波形,其峰值位移达14 cm、永久变形约为12 cm;合成孔径雷达干涉测量(interferometric synthetic aperture radar, InSAR)显示断层近场30 km×30 km区域具有明显的同震变形,断层西侧卫星视线向最大位移达15 cm,挖角乡附近的局部形变超过15 cm;野外调查推测在“二台子-爱国村”间的断层段发育地表破裂,但从InSAR形变场上难以判定,同震是否破裂到地表仍需详细的野外考察予以确定.反演了川滇块体北东边界主干断裂在泸定地震前的闭锁分...     

玉树M_S7.1级地震地表破裂与历史大地震

野外调查表明,青海玉树MS7.1级地震较清晰的同震地表破裂带由3条主破裂左阶组成,走向310°～320°,总长约31km,左旋走滑性质.另在隆宝镇东侧一带见有长约2km的雁列式张裂缝带,如以该点为破裂带的北端点,则破裂带总长约51km.地表破裂带由一系列挤压鼓包与张裂缝相间排列或雁列式裂缝组成,实测最大水平位错约1.8m.地表破裂带沿甘孜-玉树断裂展布,显示该断裂是此次地震的发震构造.甘孜-玉树断裂历史上记载过多次7级左右地震,古地震遗迹明显,具有短周期的大地震重复特征.玉树地震的孕育机制与汶川地震一样,都是青藏高原东扩、地块边界应力积累和释放的结果,不同的是玉树地震为巴颜喀拉地块与川滇块体向东不均匀挤出产生的左旋走滑型地震.     

GPS,OBP和GRACE估计地球动力学扁率变化及其地球物理激发

行星地球是一个旋转的扁椭球体,其动力学扁率(即J2)变化主要由地球系统物质流动和各圈层相互作用引起.目前国际上地球动力学扁率测定主要用卫星激光测距(SLR)资料得到,然而SLR地面观测台站少,南北半球分布不均,且不是连续观测,以及受动力学模型和常数等影响.尽管新一代重力卫星(GRACE)观测大大提高了地球重力场球谐系数低阶项一两个数量级,但对二阶项C20不敏感.本文利用全球连续GPS观测得到的地表负荷位移估计地球动力学扁率J2,并联合GPS与海底气压(OBP)以及OBP和GRACE资料分别估计地球动力学扁率J2,比较和分析地球动力学扁率多尺度变化特征及其机理.结果发现GPS单独估计结果振幅偏小,GPS+OBP,GPS+OBP+GRACE和GRACE估计J2周年变化与SLR结果非常接近.而J2半周年变化,GRACE估计相对较差,主要由于GRACE资料处理没有很好扣除约161天的S2潮汐影响.另外,GPS+OBP和GPS+OBP+GRACE估计J2季节内和年季变化与SLR结果较一致,而GRACE和GPS单独估计结果与SLR结果偏差较大.并进一步利用地球物理模型资料研究对J2变化的贡献,结果表明,地球动力学扁率J2季节性、季节内和年季变化主要由大气、海洋和陆地水地表流体质量重新分布和迁移激发.     

基于SAR差分干涉测量的张北-尚义地震震源参数反演

基于欧洲遥感卫星（ERS）合成孔径雷达差分干涉测量技术（D-InSAR）获取1998年1月10日张北-尚义地震形变图,通过弹性半空间位移模型模拟形变场,由最佳参数解得到张北-尚义地震的震源参数.模拟结果表明,张北-尚义地震的发震构造为走向95°,倾向西南,倾角30°的具右族性质逆断层,断层面长12 km,宽14 km,震源（断层面中心）深度7.5km,震中位置为北纬40°58＇,东经114°21＇.断层滑动角105.95°,滑动方向为北西I3.26°,位移量为0.728m,地震矩M0为2.69×1018N·m.根据震源参数和区域构造背景分析,断裂活动可能与张家口-蓬莱断裂北西端的活动有关.     

2022年青海门源Mw 6.6地震的发震断层及孕震构造模式

2022年1月8日青海门源盆地北缘发生Mw 6.6地震,震源机制反演表明此次地震属于左旋走滑事件.震后10 d内,近600个余震被检测到,最大余震为M 5.1级.此次地震发生在祁连-海原左旋走滑断裂系统的冷龙岭段,该断裂段全长127 km,由古地震研究确定的特征地震大小在Mw 7.3～7.5.为了更为全面理解此次地震的震源机制以及当地孕震模式,我们分析了地震波形,获取了主震和17个Ms≥3.0余震的震源机制与矩心深度.利用升、降轨道SAR数据获取的像元偏移数据和同震干涉相位(interferometric synthetic aperture radar, InSAR)确定了两条地表破裂带的位置,并利用InSAR数据反演了主震的滑动模型.研究发现,此次地震破裂带对应于冷龙岭断裂西段和托莱山断裂的阶区,发震断层存在3个形变中心,最大滑动量约为4 m,出现在冷龙岭断裂上,形变中心深度为4 km.滑动模型显示释放了累计能量～1.58×10¹⁹Nm,约合矩震级Mw 6.68,与本文利用地震学方法得到的Mw 6.58接近.结合区域活动构造特征、1986和2016年两次门源地...     

全球定位系统(GPS)测定的印尼苏门达腊巨震的远场同震地表位移

以中国的连续GPS观测资料为基础, 结合其他全球GPS观测资料, 计算了2004年12月和2005年3月发生在苏门达腊岛9.3和8.7级巨大地震的同震位移场. 2004年地震的同震永久位移的影响范围达6000~7000 km, 并呈现波状起伏的特征. 2005年地震的同震位移主要局限在震中区附近, 可能只是2004年地震诱发的强余震.     

致密砂岩层系多尺度波耗散规律及双重分形结构模型

了解地层岩石的多尺度非均质性特征对不同频率弹性波响应的影响,对基于地震资料解释与反演地下介质属性至关重要.本文联合采用超声波实验(550 kHz)、声波测井(约10 kHz)和地震勘探(约30 Hz)手段,对鄂尔多斯盆地延长组致密岩石开展观测,估算了不同频率下致密砂岩中的纵波传播速度和衰减.超声波观测结果显示,随着围压增大,孔隙度对衰减的影响有变小的趋势,而黏土含量的影响随围压增大而增大.不同频率范围内,波速和衰减随孔隙度或黏土含量的变化趋势基本一致.测量结果中值及分布特征显示,从超声波到地震频段,纵波速度趋于下降,而衰减却有增大趋势.从超声波到声波频段,速度减小更为显著,中值减小了362 m/s;而测井速度与地震观测结果较为接近.衰减在声波和地震频率之间的差异更为明显.为解释多尺度观测数据,提出了描述孔隙介质中裂隙和黏土尺寸具有自相似分布特征的双重分形结构模型.模型预测与实验结果对比揭示了相关储层黏土包体和裂隙的分形特征.本研究可为进一步解释复杂岩石宽频带波频散和非弹性现象提供理论模型.     

利用初至压缩波对2004年苏门答腊-安达曼地震破裂直接成像

2004年12月26日在印度尼西亚苏门答腊岛西侧海域发生的Mw9.3级地震, 是现代数字地震和GPS观测台网布设以来记录到的最大地震. 高质量的地震和形变测量资料, 为准确估计引发巨大海啸的地震破裂特征提供了难得的可相互验证的资料. 利用中国国家数字地震台网(CNDSN)记录到的初至压缩波反演这次地震破裂过程, 并通过对方法的改进, 克服了CNDSN稀疏的台站分布带来的影响. 所获得的地震破裂过程与全球和区域地震台网及大地测量得到的结果较为一致, 即这次地震起始于苏门达腊西北海岸并以(2.7±0.2) km/s的破裂速度向北北西方向单侧破裂, 破裂持续时间达420 s或更长, 破裂沿着安达曼海槽延伸1200~1300 km, 并且出现两次强能量脉冲.     

全球定位系统测定的尼泊尔M_w7.8级地震同震位移

据尼泊尔境内的连续GPS观测资料显示,2015年4月25日的尼泊尔M_w7.8级地震造成尼泊尔向南移动,最大位移量达到了1.88 m.利用国家重大科技基础设施项目"中国大陆构造环境监测网络"、中国气象局中日合作JICA项目及加州理工学院在尼泊尔境内所建的连续GPS观测资料,分析了此次地震对中国大陆构造形变场的震时动态响应及同震影响.结果显示,地震造成中国西藏地区毫米至厘米级的同震水平位移,最大值在震中北东方向,震中距220 km的珠峰站,约为30 mm,波及范围远至1500 km之外的高原东北部的祁连山地区.1 Hz的高频GPS数据分析表明,震中距1200 km内的站点记录了弹性震波和永久变形的叠加和破裂方向效应等丰富信息.通过分析中国境内地震活动频度并对比2011年日本宫城地震发现,本次尼泊尔地震并未造成中国境内微震活动显著增加,同震过程对中国西藏部分地区、川滇地区及南北地震带断裂带的应力扰动和构造加载有限,与其震级相对较小有关.     

上海GPS网层析水汽三维分布改善数值预报湿度场

水汽三维分布信息是数值天气预报模式的初始条件之一, 大气中的水汽具有较强的时空变化特性, 传统观测手段受时空分辨率的限制, 对模式初始水汽场的描述远远不够. 利用地基GPS信号斜路径观测值探测水汽垂直剖面信息是目前国际上的研究前沿, 利用上海GPS综合应用网GPS信号斜路径水汽观测值, 在国内首次层析了上海地区水汽三维分布. 层析中采用高斯加权函数进行水平约束, 利用数值预报结果做先验值, 考虑了边界条件. 针对无法准确给出各类观测信息先验权的问题, 在层析水汽三维分布问题中首次采用了方差分量估计, 进行验后定权; 为了降低观测噪声对参数估计的影响, 采用了抗差估计; 针对水汽密度随高度减小的特点, 垂直方向上采用随高度变化的非等权约束方法. 层析结果与MM5数值预报结果比较, GPS信号斜路径水汽观测值对数值模式的预报湿度场有明显改善, 利用GPS信号斜路径水汽观测值准实时地探测水汽的三维分布, 应用于短期数值预报中, 可改善短期预报的初始湿度场, 提高数值预报的准确度.     

2015年尼泊尔M_s8.1与M_s7.5级地震活动特征

2015年尼泊尔8.1级地震发生在喜马拉雅地震带上,为低角度逆冲型单侧破裂.余震区呈WNW-ESE展布,长轴约170 km,短轴约60 km.余震空间分布不均匀,主震和强余震分布在余震区两端,中部余震稀疏,这与8.1级地震矩释放主体区一致.7.5级地震发生在8.1级地震余震区的东部边缘,8.1级地震对其具有显著的触发作用.8.1级和7.5级地震发生在尼泊尔1505年和1934年两次大震之间的8级地震破裂空段上,1870年以来至本次地震前该破裂空段内没有发生过6级以上地震,存在6级地震背景空区.这次地震前13年,形成长约590 km的5级地震空区,震前19个月空区被打破.8.1级地震序列发生在5级地震空区的中部,其东、西两侧仍有较大范围没有发生地震,库伦应力计算表明8.1级地震对其东西两侧断层具有明显的触发作用.考虑到历史地震的离逝时间与复发周期,认为1934年地震破裂区再次发生大震的危险性较小,而1505年地震破裂区发生大震的危险性增大.2005年巴基斯坦M_w7.6和2015年尼泊尔8.1级地震的发生,表明喜马拉雅地震带已经进入了一个7级以上地震相对活跃的时段.     

汶川大地震余震序列震源机制的空间分段特征

利用汶川地震序列周围450 km范围内6 个台站记录的宽频带波形数据, 通过时间域矩张量反演方法, 得到该地震序列88 个地震的矩张量解. 结果显示震源机制解类型具有明显的空间分段特征. 沿着主破裂带方向自西南至东北共分为6 段, 震源机制类型由开始的逆冲为主, 经过各种震源机制类型交替出现的过渡带, 转变为以右旋走滑为主. 在小鱼洞镇附近出现与主破裂正交的左旋走滑型地震. 利用FMSI程序分别反演各段的应力场, 发现沿着主破裂带方向自西南至东北, 最大主应力方向由近EW向逐渐转变为NW-SE向, 最后又转变为近EW向.     

全球定位系统测定的2011 年日本宫城MW 9.0 级地震远场同震位移

据覆盖日本全境的GEONET 网络GPS 观测资料显示, 2011 年3 月11 日的日本宫城M_W9.0 级地震造成日本半岛向东移动, 最大达到了5.3 m. 利用国家重大科技基础设施项目“中国大陆构造环境监测网络”的GPS 观测资料, 分析此次地震对中国大陆构造形变场的同震影响,结果显示, 地震造成我国东北和华北地区产生毫米至厘米级的同震水平位移, 最大值为35 mm.通过应变分析发现, 地震导致东北和华北地区一系列北北东走向的断裂产生了不同程度的张性应变. 虽然在东北地区张性应变相对比较明显, 最大处约为40 nano-strain, 但对断裂带的静态库仑应力加载有限, 不会对区域地震活动产生大的影响.