华北地区GPS地壳应变能密度变化率场及其构造运动分析

利用GAMIT软件将华北地区GPS地壳形变监测网1995、1996及1999年3期观测数据同全球IGS跟踪站资料进行了统一处理,并采用基线向量单天解及其全协方差矩阵和武汉大学PowerAdj Ver3.0软件在ITRF97框架下进行了GPS网整体平差,得出了相对于欧亚板块的水平运动速度场;进一步利用GPS速度场计算华北地区GPS地壳应变率场及应变能密度变化率场. GPS结果表明:华北地区现以(3-12mm/a)±3.4mm/a速度相对于稳定的欧亚板块向东或东南方向运动;唐山-河间-磁县是一个压缩活动边界,大同-太原则是一拉张活动边界,呈略带右旋的拉张运动;北京-天津-渤海湾地区,以及邢台地区是高剪应变率地区;张家口-北京-天津-渤海湾地区,以及济南地区是高应变能密度变化率地区,地壳积累能量高,且渤海湾地区处于面膨胀变化率正、负值过渡带,那里可能存在与地震危险性有关的应变能积累背景.     

MATLAB在测量中的应用

MATLAB是一种功能十分强大的科技应用软件和编程语言，它集数值计算、可视化和编程于一体，易学易用。本文结合实例介绍了MATLAB在测量中的应用，给出了程序的源代码及详细的注释，使读者能深刻的体会到MATLAB编程的语言风格、使用的简洁方便，以及其强大的功能。     

用大地测量资料检验青藏高原隆升的两种构造演化模式

对当前地质学和地球物理学方面提出的两种青藏高原构造演化模式进行了分析，着重讨论了用大地测量资料检验构造演化模式的原理和方法，以实测资料初步验证了回流模式的合理性。     

InSAR数据处理中的误差分析

针对星载雷达 ,从InSAR的基本原理出发 ,综合推导并讨论了相位观测误差、卫星轨道误差、地形误差等各类误差源对不同工作模式 (包括重复飞行InSAR生成DEM ,Two passDInSAR监测地表形变 ,以及Three passDInSAR监测地表形变 )的影响     

利用PALSAR数据研究长白山火山活动性

利用2006-10~2011-03期间共24景ALOS卫星的PALSAR数据对长白山天池火山区的近年形变进行了研究,基于PSInSAR方法进行了数据处理,获取了长白山近年火山的时序形变场和形变速率结果。结果表明,长白山火山近年火山活动较弱,仅在火山口附近存在10mm/a的形变速率。同时利用Mogi模型对火山区的岩浆变化体积速率进行了反演,显示该区域的岩浆以~1×106 m3/a的速率补充。     

升降轨InSAR数据约束下的2007年阿里地震反演分析

2007年5月5日,西藏阿里地区发生了Mw6.1级地震。本文采用Envisat卫星的升、降轨SAR数据获取地震的精确同震地表位移,然后基于弹性半空间位错模型,分别采用MPSO非线性和最小二乘线性反演算法确定了断层几何参数和滑动分布模型。结果表明,分布式滑动分布模型能较好地解释观测到的InSAR地表形变场。地震发震断层为走向158°、倾角43°的西南倾断层,主要的滑动量集中在7~12km的深度,最大滑动量约0.3m,位于9km的深度。反演给出的地震矩为1.24×1018 Nm,与地震学结果相一致。     

联合Envisat和ALOS卫星影像确定 L'Aquila地震震源机制

2009年4月6日意大利L'Aquila地区发生了M_w6.3级地震,该地震造成了300余人的人员死亡. 本文联合不同波长、不同入射倾角的升降轨Envisat和ALOS卫星的差分干涉数据对该地震进行震源机制解的反演研究. 研究首先对卫星雷达影像进行二通差分干涉处理,获取了覆盖L'Aquila地震震区的完整InSAR同震形变场,然后结合四叉树和均匀采样方法对原始观测数据进行降采样. 在此基础上,联合GPS形变观测数据,利用弹性半空间矩形和三角位错模型,以及断层自动剖分技术对断层面进行最优离散剖分,反演获取了发震断层的精确几何参数和最优断层滑动分布,结果显示分布式三角位错滑动模型能够很好地解释观测到的地表形变场. 反演结果表明发震断层是一个以正倾滑为主兼有少量右旋走滑的盲断层;基于观测数据最优确定的断层单元的最短边长为0.4 km,最长边长为6.3 km;此次地震的滑动分布主要发生在5～14 km深度的范围内,最大滑移量为1.07 m,释放的能量为3.43×10¹⁸ N·m(M_w6.32),与地震学的研究结果非常一致.     

岗纳楼冰川表面流速时空变化特征提取及分析

冰川表面流速是反映冰川动态变化的重要指标,能够为冰川物质平衡提供重要信息.利用2016年的13景Sentinel-1A影像和合成孔径雷达(synthetic aperture radar,SAR)偏移量追踪法测定岗纳楼冰川表面流速场,并根据地表温度是否大于0℃将其分为冰封期(10月—次年3月)和消融期(4月—9月).其...     

华北地区地壳形变的GPS及地震矩张量反演分析

给出了在椭球面上利用ＧＰＳ速度场计算地壳应变率张量的公式,利用华北地区ＧＰＳ监测网１９９２年、１９９５年和１９９６年观测数据所得的速度场在椭球面上进行了应变率计算,得到了其最大剪切变态为率、主应变率的图像;并利用１９６６年以来２５次大于５．０级地震的地震矩张量,由Ｋｏｓｔｒｏｖ公式计算了地震应变率,对两种数据所得的应变率地比较分析;结合地质构造背景探讨了华北地区现今形变特征及地震危险性地区的分布。     

基于方差分量估计的2015年尼泊尔MW7.8地震同震滑动分布

基于三角位错的断层自动剖分方法,联合ALOS-2 InSAR和GPS同震位移数据反演了2015年尼泊尔M_W7.8地震的发震断层参数及滑动分布.反演中基于赫尔默特方差分量估计方法确定InSAR、GPS水平向与垂直向的3种观测位移数据的相对权比.反演结果显示,方差分量估计方法定权得到的最大滑动量5.5 m,大于验前方差定权给出的5.1 m的结果,同时前者显示同震破裂沿东南方向传播至接近M_W7.3最大余震区域,然后分别继续向北和东南传播约30 km,并形成一个滑动空区（slip gap）,而后者的结果没有表现出这种显著的滑动空区特征.方差分量估计定权得到的滑动模型能更好地解释观测到的InSAR同震形变场,其拟合残差标准差较验前方差定权的结果减小了2.5 cm,GPS的东和北分量的残差标准差分别减小了0.4 cm和0.9 cm.与两种不同的基于矩形位错的方法相比,利用方差分量估计的三角位错断层自动剖分方法的滑动误差均值和标准差相对最小,其中均值为0.037 m,标准差为0.028 m,表明该方法得到的滑动模型具有相对更好的无偏性和有效性.