2018年11月26日台湾海峡M_S6.2地震发震构造研究

本文利用福建省地震台网、广东省地震台网和台湾"中央"气象局17个台的宽频带记录,使用CAP方法反演了2018年11月26日台湾海峡M_S6.2地震震源机制解,得到节面1走向/倾角/滑动角为89°/82°/-173°,节面2走向/倾角/滑动角为358°/84°/-7°,最佳拟合深度14km,矩震级5.8.使用双差定位获取了94个M_L2.0以上地震的精定位结果,结果显示,主震位于北纬23.36°,东经118.62°,震源深度10.43km.根据小震分布和构造应力场反演得到余震断层面走向和倾角分别为88°和60°.研究认为,台湾海峡6.2级地震发震构造为近EW向的台湾浅滩断裂,受南海板块张裂拉伸发育而成,孕震过程中有东山隆起东缘断裂的参与,推测在菲律宾板块对欧亚板块NW-SE向挤压碰撞背景下,近EW向的台湾浅滩断裂与近NS向的东山隆起东缘断裂交接部位属于强度薄弱区,最终产生高倾角右旋走滑错动而引发地震,余震主要沿台湾浅滩断裂分布.     

GNSS形变监测静态基线网中多基线联合解算两步法

GNSS(global navigation satellite system)静态基线网的解算方法通常采用单基线解算模式或多基线解算模式。这两种方法均有一定的局限性。本文综合两类方法的优点,提出了单基线和多基线联合解算模式的两步法。首先基于单基线模式利用多频率的载波相位观测值解算监测站间的双差整周模糊度,并使用严格的固定准则确定模糊度;然后利用多基线模式进行联合坐标和天顶对流层参数的估计,再结合最小二乘估计解算定位结果。选取国内某水电站静态形变监测数据为研究对象,进行基线解算方法测试和分析。实验结果表明多基线联合解算两步法在精度和可靠性上有较大优势,平面定位精度为亚毫米级,高程定位精度提高1～2倍;并可以明显改善基线网的解算速度,在高精度变形监测应用中具有可行性与通用性。     

基于StarFire星站差分GPS系统的精密静态解研究

针对StarFire星站差分GPS系统静态观测定位结果存在收敛和发散现象,提出了直接平均法、平均中位法和密度中位法以获取其精密静态解。平均中位法考虑了收敛和发散过程,其计算结果要优于直接平均法;但由于没有相应的收敛点和观测点定位模型,平均中位法获取的精密静态解存在一定的不确定性;基于观测值密度分布的密度中位法不需要确定收敛点和发散点,就能得到精密静态解。实例表明,当观测环境相对稳定,观测过程没有明显的收敛和发散现象时,直接平均法、平均中位法和密度中位法的精密静态解基本相同,密度中位法略优;当观测环境差,观测过程出现明显收敛和发散现象时,密度中位法最优,平均中位法次之。     

洞庭湖南荻(Triarrhena lutarioriparia)腐解对水质的影响

洞庭湖区南荻下游产业造纸行业于2019年整体退出,使得大量南荻无法被刈割利用而滞留于湖区腐解.为评估未刈割南荻残体腐解对洞庭湖水环境影响,本研究探究了沉积物湖水模拟体系内不同生物量南荻腐解过程水质变化特征.实验结果表明:在90 d腐解过程中,不同南荻生物量下(0～2 g/L)腐解体系上覆水COD、氨氮(NH3-N)、总氮(TN)、总磷(TP)水质指标浓度均呈现"急速上升快速下降缓慢恢复"趋势,其中,COD和TP浓度均在腐解第5天时急速上升至最高值.腐解体系内南荻生物量对水质影响显著,相同腐解时间下,随着南荻投加量增加,上覆水溶解氧(DO)浓度越低,NH3-N、TN、TP、COD水质指标浓度越高.当南荻投加量为2 g/L时,上覆水COD、NH3-N、TN和TP水质指标浓度可分别高达46.00、4.11、6.84和1.14 mg/L,分别是空白组的2.60、1.58、1.42和8.14倍.当腐解时间为90 d时,南荻投加量为0～1 g/L条件下,腐解系统上覆水水质基本恢复至未腐解时水平;当南荻投加量增加至2 g/L,上覆水NH3-N、TN和TP浓度仍分别为未腐解时的0.59、0.35和2倍,DO仅为初始浓度的16.6%.总体上,南荻腐解可造成短期内营养盐及有机组分浓度显著上升,且其腐解对水质的影响程度随着其生物量的增加而显著增强,以上结果可为造纸业退出后洞庭湖南荻管控提供科学依据.     

联合星载GPS和KBRR星间测速数据反演GRACE时变重力场模型

高精度GRACE卫星时变重力场反演一直是卫星重力测量中的难题.为了恢复高精度的时变地球重力场模型,本文联合GRACE卫星的星载GPS和KBR星间测速观测数据,在对GRACE卫星进行精密定轨的同时,解算出60阶月平均地球重力场模型.通过对GRACE卫星的定轨精度、星载GPS相位和KBR星间测速数据的拟合残差以及时变地球重力场模型解算精度等分析,表明:(1)与美国宇航局喷气推进实验室(JPL)发布的约化动力学精密轨道相比,本文确定GRACE卫星轨道三维位置误差小于5 cm.(2)星载GPS相位数据拟合残差为5~8 mm,KBR星间测速数据拟合残差为0.18~0.30μm·s~(-1).(3)解算的月平均重力场模型与美国德克萨斯大学空间研究中心(CSR)、德国地学研究中心(GFZ)和JPL发布的RL05模型精度接近,时变信号在全球范围内具有很好的空间分布一致性.通过计算亚马逊流域和长江流域的水储量变化,本文与上述三个机构的计算结果无明显差异,且相关系数均达0.9以上.可见,本文建立的卫星轨道与重力场同解算法具有反演高精度GRACE时变重力场能力,为我国卫星重力场反演提供了重要的技术支持.     

盖州青石岭地震序列发震构造初探

利用双差定位方法对盖州青石岭震群2012年2月至2015年8月的地震活动进行了重新定位,并使用CAP方法和P波初动法计算了M_L≥4.0地震的震源机制解,之后结合盖州地区的地震地质资料,分析了青石岭震群的发震构造.结果表明:青石岭震群在平面上呈NW向分布,地震活动主要分布在6 km×3 km的矩形范围内,震源深度为7—10 km;较大地震的震源机制解的走向与精定位后地震的优势分布方向一致;综合分析双差定位结果、震源机制解和发震区的地震地质等资料,初步认为九寨—盖县北段西北侧存在NW向次级铲式正断层,青石岭震群即为该断层在区域应力场作用下不断地左旋走滑-拉张错动造成的.      

日本海沟俯冲带MW9.0地震震源区应力场演化分析

于2011年3月11日发生在日本东北部的M_W9.0级逆冲型板间地震是日本有地震记录以来震级最大的一次地震.本研究基于NIED F-net矩张量解目录中的震源机制解,选取两个长轴相互垂直的矩形区域进行应力场2D反演,获取了日本海沟俯冲带地区应力场的空间及时间分布图像.结果表明：主震前,俯冲带地区应力状态在空间上大体趋于一致,即应力轴（P轴、σ₁轴及S_Hmax轴）系统性地倾向板块汇聚方向,P轴、σ₁轴倾角整体偏缓（<30°）,且远离震源区及日本海沟东侧区域内的应力轴倾角普遍大于主震震源区内应力轴倾角;主震前,受2003年5月26日在宫城县北部发生的M_W7.0地震影响,位于M_W9.0地震震源区西北侧的应力场出现明显扰动,σ₁轴倾向顺时针偏转150°~180°,并于之后大体恢复至震前状态,同期其他地区没有明显变化,这种情况可能和主震断层局部（深部）的前兆性滑动有关;主震后,距离震源区较远处应力场变化不大,主震震源区内应力场发生显著改变,P轴及σ₁轴均以大角度（>60°）倾伏于板块汇聚方向,S_Hmax轴顺时针偏转60°~90°且在日本海沟附近普遍平行于海沟轴.这项研究以时空图像的方式展示了大地震前应力场变化的特点,反映了大地震孕震过程中构造与地震的相互作用,对于理解大地震孕震过程有重要意义.     

采用CAP方法反演2016年广西苍梧M_S 5.4地震震源机制解

北京时间2016年7月31日广西梧州市苍梧县发生M_S 5.4地震,基于海南地震台网数字波形资料,采用CAP方法反演震源机制解。结果显示,此次M_S 5.4地震震源深度较浅,最佳深度为5.1 km,其中节面Ⅰ参数为:走向340°,倾角37°,滑动角-18°;节面Ⅱ参数为:走向85°,倾角79°,滑动角-125°。初步推断苍梧M_S 5.4地震破裂面运动以走滑为主,兼有正断性质,反演参数与中国地震台网中心结果较为一致。     

首都圈地区震源机制解及现今构造应力场时空变化特征研究

基于2009年1月至2017年11月首都圈地区发生的8 061个地震事件的23 293条P波初动极性数据,采用改进的格点尝试法计算了首都圈地区单次地震的震源机制解和小震综合断层面解。在初步分析这些数据的基础上,利用计算得到的单次地震的震源机制解和搜集到的已有历史地震的震源机制解数据,运用线性反演法对首都圈地区构造应力场的时空变化特征进行了研究。结果显示:① 研究区的地震震源机制解类型以走滑型为主,正断型次之,这些地震震源机制解的P轴方位大都为ENE向和近EW向,与该地区的构造应力场方向基本一致,仅有个别地震的P轴方位为NNW向;② 首都圈地区的构造应力场具有较好的一致性和连续性,最大主应力轴方位由西部的ENE向至东部的近EW向呈现顺时针旋转的趋势,应力类型整体上为走滑型,这与以往的研究结果相一致;③ 通过与已有研究结果相比较认为:京西北地区现今构造应力场是相对稳定的,最大主应力轴未呈明显改变;唐山地区和北京地区的构造应力场（最大主应力轴）在1976年唐山地震前后可能发生了变化,唐山地震后一年至今（1977—2017年）是否发生变化,依据现有的计算结果尚不得而知,需要更多的研究来进一步验证.