两次大柴旦Mw6.3地震间地表形变的InSAR观测及与同震破裂的联合分析

基于2008-11-10、2009-08-28两次大柴旦Mw6.3地震期间的6景Envisat SAR影像,利用小基线集技术提取相应的InSAR地表形变场。结果显示,震中周缘区域未观测到显著的震后形变空间模式,6个特征点区域的形变序列也未显示由快到慢的时变特征。基于震后形变理论对获取的地表形变进行余滑解释,分析两次地震同震破裂的主要参数,建立震时及期间地下断层滑移过程。两次地震分别在不同深度发生破裂且平均滑移量差异较大,前次地震未产生显著的浅层震后余滑。地震区域岩石圈8.2~23.7km深度内在同震阶段的滑移亏损可能以震间无震滑动的形式予以补偿。     

台湾美浓MW6.4地震InSAR形变场初步分析

利用欧洲航天局最新发射的宽刈幅、高分辨率Sentinel-1A卫星,第一时间获取2016-02-06台湾美浓MW6.4地震干涉像对,使用D-InSAR技术获取美浓地区的同震形变场。利用震中附近39个同震GPS观测进行对比验证后显示,InSAR获取的同震形变场精度优于1 cm（3σ）。形变发生在宏观震中30 km×30 km范围内,主要表现为沿雷达视线向抬升,最大抬升约12 cm。从形变场空间分布特征可以看出,空间连续性较好,说明宏观震中附近地表未发生明显破裂。宏观震中并不与震中位置重合,而是位于震中西部约15 km处。     

新疆精河2017年6.6级地震对2018年5.4级地震的应力影响分析

基于2017-08-09新疆精河6.6级地震同震破裂位错模型和区域地壳速度结构,利用粘弹性模型计算其地表位移、同震及震后应力场扰动变化,分析其对2018-10-16精河5.4级地震的应力触发关系。结果表明,精河6.6级地震造成的同震位移主要为地表垂向位移,最大位移量为84 mm,地表水平位移最大位移量为25 mm;地震引起发震断裂及其与博阿断裂交会区、库松木楔克断裂东段及分支断裂、喀什河断裂东段库仑应力增加,这些断裂靠近精河6.6级地震震中部分库仑应力增加超过10 kPa,但5~10 a尺度内的应力扰动基本为弹性响应,并未随时间显著变化。精河6.6级地震造成的精河5.4级地震震中位置同震库仑应力变化显著低于10 kPa的触发阈值,两者不存在同震或震后应力触发关系。     

2015年尼泊尔Mw7.9地震震前断层闭锁与同震位移特征分析

对尼泊尔震前的GPS水平速度场进行融合处理,获得跨喜马拉雅中东段沿N15°E方向的GPS水平速度场剖面。采用弹性半空间矩形位错模拟反演,结果表明,尼泊尔地震前喜马拉雅主前缘断裂带存在浅部闭锁,深部无震蠕滑,闭锁深度为22.5km,蠕滑段滑动速率为19.0mm/a,断层倾角为10°。GPS观测的同震形变场揭示了尼泊尔地震引起的地壳形变特征,最大同震位移位于尼泊尔境内的KKN4,该站向南移动1.89m;地震还在我国藏南地区造成最大0.54m的永久形变。距离震中400km以外GPS观测到的同震形变微弱,在误差范围之内。依据GPS同震水平位移在N15°E方向的位移剖面,采用矩形位错模型简单模拟了尼泊尔地震的同震破裂。结果显示,GPS同震位移剖面可以用喜马拉雅冲断带前缘主断裂(MFT)以北的基底低角度逆断层引起的弹性位错模拟,浅部的主前缘逆冲断裂、主边界逆冲断裂和主中央逆冲断裂等分支并没有破裂。     

中国玉树 Mw6.9地震InSAR地表形变特征分析

利用ALOS PALSAR和ENVISAT ASAR雷达遥感干涉像对,获取了2010年4月14日玉树 M w6.9地震的InSAR同震形变场.对形变结果的分析表明:至少有3次地表破裂;沿断层走向的形变分布范围远大于垂直断层分布范围,且对断层两侧的影响范围在10 km左右;最大的视线向形变量达54 cm,最大水平位移达180 cm;形变分布特征与左旋走滑断层特征吻合.对比ALOS PALSAR和ENVISAT ASAR的形变测量结果,发现两种模式获取的结果非常一致.另外基于InSAR获取的宏观震中与野外考察结果也非常一致,表明InSAR是宏观震中快速定位的简单快捷方法.     

2020年土耳其MW5.7地震InSAR同震形变特征与破裂滑动分布

为准确认识2020-06-14土耳其MW5.7地震的发震位置、构造特点以及地震危险性,利用D-InSAR技术对Sentinel-1A数据进行处理,基于GACOS进行大气校正,获得视线向（line of sight,LOS）同震形变场。降轨LOS向同震形变场显示,断层北侧抬升,最大抬升形变量约8.87 cm;南侧沉降,最大沉降量约-7.75 cm。以LOS向同震形变为约束,先采用贝叶斯自举优化法反演发震断层几何参数,然后使用有限断层方法反演地震破裂滑动分布。结果显示,断层走向约257.48°±0.65°、倾角约79.69°±0.98°、滑动角约154.2°±3.8°,震中位置为（40.754°E,39.389°N）,破裂区域长度约8 km、宽度约6 km,破裂的最浅埋深约0.8 km、最大埋深约8.9 km,最大滑动量约0.57 m,对应深度约4.278 km。地震释放地震矩约4.54×1017 Nm,对应矩震级MW5.7,与土耳其灾害和紧急情况管理局公布的结果一致。此次土耳其地震受近东西向的潜伏断层控制,以右旋走滑为主兼具少量的逆冲特性。地震造成部分地区的库伦破裂应力增量超过0.1 bar,这些区域未来的地震危险性值得关注。     

利用GPS观测资料分析2015年尼泊尔MS8.1地震震前及震后形变

通过对尼泊尔MS8.1地震前后附近区域GPS台站记录到的观测数据进行处理,获得了震区以及中国青藏高原地区地震前后GPS站点速度场以及震后形变场。震前速度场显示,喜马拉雅构造带整体呈现出约16 mm/a的压缩特征。同时,震前喜马拉雅构造带根据形变特征可分为东、中、西3段,其地震发生在中段,主要以北向挤压为主,而东西两段分别具有逆时针旋转和顺时针旋转的特征。震后GPS站点形变场显示,此次地震对中国新疆、青海、西藏等地区的影响较大,其最大震后位移达20 mm左右。震后速度场显示,本次地震对尼泊尔地区以及中国藏南地区的构造形变影响较大,主要表现为喜马拉雅构造带的年推挤速度减小,藏南地区的南北向运动速率减小,而东西向速度有增大的现象。这一现象可能对藏南地区的走滑断层有较大影响。     

荣县地震InSAR同震与时间序列形变初步分析

借助Sentinel-1A卫星升降轨数据，利用合成孔径雷达差分干涉测量(DInSAR)和Okada均匀滑动模型反演获得四川荣县地区三维同震形变场及荣县地震震源参数。结果表明，此次地震造成的断裂现象大体沿NS向，两侧沉降区域中间伴随抬升现象，最大垂直形变量约18 mm；荣县地震发震断层符合逆断层特性，断层长约9.5 km，宽约2.1 km。采用时间序列分析技术获取该地区2018-11～2019-03各时间段的累计形变，结果发现，该区域在震前形变波动较小，地震是造成累计形变较大的主要原因，震后地表变化趋于稳定。     

俯冲带地震循环的数值模拟——以日本Tohoku MW9.0地震为例

基于滑移弱化摩擦准则,以日本2011年Tohoku MW9.0地震为例,建立一个以物理规律控制地震循环过程的二维有限元数值模型。结果显示,参考模型在1 000 a间的6次大地震表现出特征地震的规律,地震重复周期约161 a,单位破裂长度地震矩为1.13×1020 Nm/km,在两次大地震中间会发生一次5.62×1018 Nm/km的小地震。参考模型的数值结果与地表同震GPS位移、震间GPS速度分布具有较好的一致性。模型弹性参数的不确定性对同震以及震间形变的影响有限,模型粘性参数主要影响震间形变场。数值计算也显示,假设震间形变仅由断层运动规律所决定,那么在一个地震周期内,模型空间重力异常基本上随时间均匀变化,在大陆一侧距海沟100 km处可达-370 μGal;速度场的变化主要发生在震后约5 a的时间内,此后基本保持稳定增加。     

高频GNSS资料的强震震相识别与震中确定

利用2011日本东北Mw9.0地震和2008汶川Mw7.9地震若干台站的高频（1 Hz）GNSS观测资料,通过精密单点定位方法解算出地震发生时的位移形变波形。对位移波形一阶差分得到震时高频GNSS速度时间序列,采用S变换谱的方法对该速度时序进行地震波到时的拾取,即震相识别。将地震波拾取结果用于震中位置和发震时刻反演,并与美国地质调查局USGS公布的资料对比。结果发现,两者日本地震震中位置相差约16 km,发震时刻相差约0.7 s;汶川地震震中相差约4 km,发震时刻相差约0.4 s。     

震前突变——水准资料反映的2013年岷县-漳县MS6.6地震孕震特征

利用兰（州）天（水）武（都）水准监测网1973～2014年水准资料，获取2013年岷县-漳县M_S6.6地震前后震中及其所在区域的垂直形变场。结果表明：1）从1999～2006年区域垂直形变场来看，本次地震位于NNW-NWW向隆起带至近EW向相对下沉带的转换位置，也是反映岷山与秦岭之间构造活动的差异带附近。2）从震中以西附近沿NNE向的陇西-岷县水准剖面演化过程来看，垂直形变具有明显的阶段性变化特征，最终表现为震前突变。由分时段结果来看，1973～1993年，在漳县至岷县之间出现最大上升幅度近40 mm的压性隆起区；1993～2008年，隆起区消失且转平，期间剖面整体呈波动状态，幅度在±10 mm之间，西秦岭断裂、临潭-宕昌断裂两侧无明显的垂直差异运动；2008～2011年，相对陇西盆地，漳县至岷县之间跨临潭-宕昌断裂两侧出现近40 km长的下沉梯度带，其中距震中16 km以内的武定80、武定78、武定77等3个水准点同步出现向下突变，接近或大于1973年以来正常变化幅度的2倍；2011～2014年，漳县、岷县间相对陇西盆地恢复以逆冲隆升为主，但震前临潭-宕昌断裂两侧出现的形变梯度带已不存在，而上述3个水准点也恢复至震前波动区间。3）在总结岷县-漳县M_S6.6地震前形变前兆的基础上，针对目前关于震前是否出现形变加速或突变特征的争议认为，只要形变资料的时空尺度足够，至少对某些构造区域的中强地震还是可以捕捉到震前的突变现象。     

GPS�ⶨ���Ĵ�«ɽMw6.6����ͬ���α�

??????????????Χ??GPS??????????????з?????????t???????λ?????????α?????50 km??Χ????GPS????α???????????????????70 mm????????????30 s???????λ??????????????????α????????λ???????????У???С?????о???????????????29 km????????λ???21 mm????????150 km???λ???1??2 mm????GPS???????α????????????????????????????????????????????λ?????и???20 km????????廬?????0.6 m??GPS????????α??о????????????????????????????????ε?????????????????     

����Ms6.6�����InSAR�۲⼰�ϲ�λ����

????3??Envisat?????????,?????2008??10??6?????????Ms6.6??????????????α?,????????????4????????????δ?????????????α?,?????????????????20 km??20 km??Χ??????0.3 m????????Ρ??α??????????????????????????????????????t???????????λ???????????????????????????λ?????????????????????????????????????????????λ???3 m, ????λ??90.374??E, 29.745??N????Mw6.35??????????????189?????60????????????????′???????ε?μ?????????????12 km,??11 km,???9.5 km???????????????????????????????????????????????й??     

安徽省CORS坐标时间序列噪声和速度场分析

利用GAMIT/GLOBK软件解算安徽省卫星定位综合服务系统（AHCORS）2013-01～2018-06数据,得到参考站坐标时间序列,并对其进行噪声分析。结果表明,安徽省境内参考站的最佳噪声模型为白噪声+闪烁噪声（WN+FN）模型。在顾及有色噪声影响时,AHCORS在ITRF2008框架下的平均运动速度为31.72 mm/a,方向为E22.76°S;以欧亚板块为参考的平均运动速度为6.28 mm/a,方向为E2.65°N。垂直方向上整体平均运动速率为-0.71 mm/a,以淮河为界呈现出南方隆升、北方沉降的趋势。其中,大别山地区的隆起速率大于江淮丘陵区和黄山地区,江淮丘陵区抬升速率相对平缓,淮河以南一带平均抬升速率为2.94 mm/a。淮北地区沉降明显,平均沉降速率为10.92 mm/a,最大沉降速率为32.82 mm/a,发生在宿州砀山（SZDS）站,可能是近年来开采矿产资源所致。     

基于GPS的2021年玛多MW7.4地震同震静、动态形变场

利用青海省CORS网GPS连续观测数据获取2021年玛多M_W7.4地震的同震形变场以及震时地表运动状态。结果表明,断层南北两侧的近场测站分别展现出南东向和北西向运动,符合左旋走滑机制。同震形变集中在震中距300 km范围内,震中距150 km以内的站点均能反映出cm级的位移,最大水平向位移为28.3 cm(JDUO站)。高频GPS动态形变提取出的永久位移与静态解算结果相当,其动态波形最大峰值为49.9 cm(KANQ站)。依据震中距和波形初动时间估计得到地震波速为2.8~4.9 km/s,断层东端的站点估算速度大于其他站点,可能与此次破裂的方向性或者破裂传播速度有关。依据震级经验公式,动态波形估算震级在M6.8~7.6区间,拟合平均震级为M7.35±0.15。若在实时条件下,震后70 s内可得到稳定的震级。     

����PS-InSAR�������ݱ��������ر��������

???27??ENVISAT ??????????????????????????PS??InSAR??????????????д??? ????????????????2006-2010?????α????о??????????1??????????????????????????????????????????????????????7??????????????г??????????????????????????????????????????????????????????LOS??????-27.19±5.66 mm/a??-24.65±5.14 mm /a??2?????????????????????????230 km²??3???????????????????????10.92±3.24 mm/a????????20????80????????????????????ε???????????????????4?????????????????????????????????????????????????     

山东省黄河北煤炭矿区地表沉陷预测与综合治理研究

黄河北煤炭矿区为鲁西煤炭基地最后一块整装后备资源,阳谷茌平煤田和黄河北煤田其开发建设具有重要的意义。该文采用实际观测资料利用回归分析法获得的岩移参数,对山东省黄河北矿区内上组煤开采后引起的地表沉陷进行了预测,预测结果表明:地表沉陷总面积达到3 724.5 km~2,最大下沉深度为5 455 mm;其中,阳谷茌平煤田最大下沉深度为5 455 mm,地表沉陷面积为1 188.0 km~2;黄河北煤田最大下沉深度为2 976 mm,地表沉陷面积为2 536.5 km~2。地表沉陷预测结果与多年水位观测数据进行对照,煤炭开采会造成1 091.5 km~2季节性积水区和62.3km~2的常年积水区。地表沉陷对地面、地物产生不同的影响和危害,其中对农业的影响最大。地表沉陷形成的塌陷裂缝、不同坡度耕地、大型塌陷坑等可分别采取黄土填堵、大型设备平整、退耕还林、煤矸石充填、修复为鱼塘等治理措施。该文对下一步黄河北矿区开发过程中,相关政府部门采取的规划控制措施提出了建议。