GEODETIC AND TELESEISMIC CONSTRAINTS ON SLIP DISTRIBUTION OF 2015 MW6.4 PISHAN EARTHQUAKE

On July 3^rd, 2015, a M_W6.4 earthquake occurred on Pishan County, Xinjiang, located in the front of western Kunlun thrust belt, which is the largest earthquake(M_W6.0~7.0)in the past 40 years in this region. In this study, we collected both the near-filed geodetic coseismic deformation observations including 4 GPS sites and one high-resolution ALOS-2 InSAR imagery, and far-field teleseismic P waveforms from 25 stations provided by IRIS/USGS, to invert the fault parameters(strike and dip)and coseismic rupture model of 2015 M_W6.4 Pishan earthquake. Using the finite fault theory, a non-linear simulated annealing algorithm was employed to resolve our joint inversion problem. The strike (120°~130°) and dip angle(35°~40°)of optimal models are different from that of some previous studies, and the dip change is strongly constrained by combined data than that of strike. In fixing the geometric parameters of optimal fault model, we also considered data weight(5)(geodetic data/teleseismic P waveforms)and constrained weight from moment and smooth factor(2.5). Clearly, our results indicate that the slip distribution mainly concentrates in the depth range from 9 to 16km and a length range of 20km along the strike direction, which is similar to the spatial distribution of the relocated aftershocks. The maximum slip is~95cm. The seismic moment release is 5.45×10¹⁸N·m, corresponding to M_W6.42. Compared with the single data set, geodetic data or teleseismic waveform, our joint inversion model could simultaneously constrain the seismic moment and slip distribution well, thus avoiding effectively a lower-resolution rupture distribution determined by teleseismic-only inversion and a bias released moment estimated by the geodetic-only inversion. Importantly, we should consider both the near-field geodetic data and far-field teleseismic data in retrieving the rupture model for accurately describing the seismogenic structure of active fault in western Kunlun region.     

2011年日本东北大地震（M_W=9.0）震间与震前变形场特征及其对强震预测的启示

强震震前(preseismic)动力学过程的研究对于地震预测具有十分重要的意义,但由于观测资料的限制,目前对强震前孕震区力学状态及其演化过程的认识还非常有限.2011年日本东北9.0特大地震(Tohoku-Oki)发生在GPS观测台站最为密集的地区,为研究特大地震震间(interseismic)与震前的变形状态提供了难得的机会.文中将利用日本东北大地震之前连续的GPS观测资料,分别计算震间与震前的速度场与变形场.通过对比分析发现,日本东北地区(Tohoku)震前的应变状态与震间的有很大的不同,震间的变形主要受到太平洋板块向日本海沟北西西向的俯冲挤压作用所控制,其主压应变以近东西向压缩为主,日本东北地区的运动方向与太平洋板块的运动方向大体一致.但是,临近地震前(震前)日本东北地区的运动方向发生了很大变化,震前30天的连续GPS观测结果显示,速度场的优势方向经常变换,间歇性地出现与太平洋板块运动方向相反的情况.这意味着震前孕震区的力学状态发生了很大的改变.这种变化可能与震前破裂成核或慢滑移及慢地震等过程有关,这些过程将加速或促进大地震的发生,从而为大地震的发生准备了力学条件.值得特别强调的是,这些现象都是可以通过直接观测能够发现的大地震之前的异常现象.由此可见,加密GPS站点进行连续观测,寻找震前变形异常区以及探索异常的物理机制对于地震预测预报有重要的科学意义.     

基于GPS应变、地震应变率与震源应力场对帕米尔高原现今构造变形特征的分析

研究帕米尔高原的构造变形特征对于理解印度板块向北推挤过程中的应变分配方式以及应力转换模式具有重要的意义.本文利用区域GPS应变场、地震应变场与震源应力场分析帕米尔高原的构造形变特征.主要结论为:(1)该区域变形主要以NNW-SSE或近N-S向的挤压为主,在高原内部伴有明显的近ENE-WSW或E-W向拉张,应力方向在帕米尔高原与塔吉克盆地区域呈现逆时针旋转的趋势,而在塔里木盆地则显示几乎与帕米尔高原的一致的应力状态,这可能与两侧盆地块体的强度差异有关.(2)安德森断层参数A∅显示帕米尔高原北缘与西侧区域为逆断层应力状态,在高原内部为正断层应力状态,这与GPS应变的结果显示的应变主要集中在主帕米尔断裂与阿莱谷地附近而在高原内部应变较低是一致的,另外应力在喀喇昆仑断裂北段的方向基本平行于断层走向,以及断层北端较低的滑动速率,这说明了地壳挤压缩短可能是帕米尔高原主要的的构造变形特征,并不支持由于边界走滑断裂导致的应变分异或者块体挤出的模式.(3)综合考虑地震应变方向与SHmax从帕米尔北部NNW-SSE方向到天山北部的近N-S方向的转换,GPS应变方向在帕米尔高原两侧盆地都存在不同程度的旋转,应力场安德森参数A∅显示的应力状态以及SKS的结果显示的近ENE-WSW方向,我们认为印度板块向北推挤与天山造山带碰撞导致帕米尔高原不对称的径向逆冲是帕米尔高原现今构造变形的主要成因与构造模式.     

基于稀疏约束反演的三维混采数据分离

混合震源采集技术相对于传统的地震数据采集,在极大提高采集效率的同时引入了混叠噪声,很大程度上影响了成像结果的精度.二维混采数据中,我们通常利用混叠噪声在非共炮域呈非相干分布这一特点来压制混叠噪声,从而实现混合震源数据分离.相对于二维混采数据,三维混采数据具有数据量巨大,构建混合震源算子困难,混合度的增加引入了高强度混叠噪声的特点.针对上述问题,本文采用稀疏约束反演方法在Radon域实现混采数据分离,混叠噪声强度比较大的情况下,稀疏约束反演方法能够得到更高精度的分离结果;利用震源激发的GPS时间通过长记录的方式在共接收点道集对上一次迭代分离结果做混合、伪分离,实现了单个共接收点道集自身混合、伪分离,避免了对整个数据做运算,同时不需要构建混合震源算子.通过模拟数据和实际数据计算来验证上述方法的适用性.     

COSEISMIC DISPLACEMENT FIELD OF THE WENCHUAN EARTHQUAKE DERIVED FROM STRONG MOTION RECORDS AND APPLICATION IN SLIP INVERSION

The development of high-rate GNSS seismology and seismic observation methods has provided technical support for acquiring the near-field real-time displacement time series during earthquake. But in practice, the limited number of GNSS continuous stations hardly meets the requirement of near-field quasi-real-time coseismic displacement observation, while the macroseismographs could be an important complement. Compared with high-rate GNSS, macroseismograph has better sensitivity, higher resolution(100~200Hz)and larger dynamic range, and the most importantly, lower cost. However, baseline drift exists in strong-motion data, which limits its widespread use. This paper aims to prove the feasibility and reliability of strong motion data in acquiring seismic displacement sequences, as a supplement to high-rate GNSS. In this study, we have analyzed the strong-motion data of Wenchuan M_S8.0 earthquake in Longmenshan fault zone, based on the automatic scheme for empirical baseline correction proposed by Wang et al., which fits the uncorrected displacement by polynomial to obtain the fitting parameters, and then the baseline correction is completed in the velocity sequence. Through correction processing and quadratic integration, the static coseismic displacement field and displacement time series are obtained. Comparison of the displacement time series from the strong motions with the result of high-rate GPS shows a good coincidence. We have worked out the coseismic displacement field in the large area of Wenchuan earthquake using GPS data and strong motion data. The coseismic displacement fields calculated from GPS and strong motions are consistent with each other in terms of magnitude, direction and distribution patterns. High-precision coseismic deformation can provide better data constraint for fault slip inversion. To verify the influence of strong-motion data on slip distribution in Wenchuan earthquake, we used strong motion, GPS and InSAR data to estimate the stress drop, moment magnitude and coseismic slip model, and our results agreed with those of the previous studies. In addition, the inversion results of different data are different and complementary to some extent. The use of strong-motion data supplements the slip of the fault in the 180km segment and the 270~300km segment, thus making the inversion results of fault slip more comprehensive. From this result, we can draw the following conclusions:1)Based on the robust baseline correction method, the use of strong motion data, as an important complement to high-rate GNSS, can obtain reliable surface displacement after the earthquake. 2)The strong motion data provide an effective method to study the coseismic displacement sequence, the surface rupture process and quick seismogenic parameters acquisition. 3)The combination of multiple data can significantly improve the data coverage and give play to the advantages of different data. Therefore, it is suggested to combine multiple data(GPS, strong motion, InSAR, etc.)for joint inversion to improve the stability of fault slip model.     

CHARACTERISTICS OF SATELLITE GRAVITY FIELD AND SEISMOGENIC MECHANISM IN THE LONGMENSHAN FAULT ZONE

Longmenshan fault zone is a famous orogenic belt and seismic zone in the southeastern Tibetan plateau of China. The Wenchuan M_S8.0 earthquake on May 12, 2008 and the Ya'an M_S7.0 earthquake on April 20, 2013 occurred in the central-southern part of Longmenshan fault zone. Because of its complex geological structures, frequent earthquakes and special geographical locations, it has attracted the attention of many scholars around the world. Satellite gravity field has advantages in studying gravity field and gravity anomaly changes before and after earthquake. It covers wide range, can be updated regularly, without difficulty in terms of geographical restrictions, and is not affected by environmental factors such as weather, terrain and traffic. Therefore, the use of high-precision Earth satellite gravity field data inversion and interpretation of seismic phenomena has become a hot topic in earth science research. In order to understand satellite gravity field characteristics of the Longmenshan earthquake zone in the southeastern Tibetan plateau and its seismogenic mechanism of earthquake disasters, the satellite gravity data was used to present the terrain information of the study area. Then, by solving the regional gravity anomaly of the Moho surface, the crustal thickness of the study area was inverted, and the GPS velocity field data was used to detect the crustal deformation rate and direction of the study area. Combining the tectonic setting of the Longmenshan fault zone and the existing deep seismic sounding results of the previous researchers, the dynamic characteristics of the gravity time-varying field after the earthquake in the Longmenshan earthquake zone was analyzed and the mechanism of the earthquake was explored. The results show that the eastward flow of deep materials in the eastern Tibetan plateau is strongly blocked at the Longmenshan fault zone. The continuous collision and extrusion process result in a "deep drop zone" in the Moho surface, and the long-term stress effect is conducive to the formation of thrust-nappe and strike-slip structures. The Longmenshan earthquake zone was in the large-scale gradient zone of gravity change before the earthquake, the deep plastic fluid material transport velocity differed greatly, the fluid pressure was enhanced, and the rock mechanical strength in the seismic source region was weakened, which contributed to the intrusion of crustal fluid and the upwelling of the asthenosphere. As a result, the continuous accumulation of material and energy eventually led to continuous stress imbalance in the deep part and shear rupture of the deep weak structure, causing the occurrence of the thrust-nappe and strike-slip earthquake.     

GPS掩星技术低轨卫星计划的现状及进展

主要介绍了几个国际上已经实施的GPS掩星计划（GPS／MET，Φrsted,Sunsat,SAC-C,CHAMP)和在研的（COSMIC，ACE）低轨卫星计划的现状及取得的进展。通过总结可知，下一代GPS掩星接收器应具有体积小、质量轻（几百克）、低功耗的特点；能提供实时的卫星物理状态参数；飞行器上的自主计算和控制；星上数据通信和命令解释；所有跟踪数据应符合厘米级精密定轨的要求；能提供掩星和海洋反射实验的测量。     

基于GIS的中山站至Dome A冰盖考察数据处理和信息提取

中山站至DomeA冰川学考察断面是国际ITASE计划的核心断面之一。首次利用GIS开展了该断面的数据处理与信息提取 ,采样点的布设与管理。介绍了ADD数字地图拼接与裁剪 ,野外数据和BEDMAP数据地理坐标投影变换的方法。利用ArcView内插等高线方法对GPS导航仪高程数据进行了校正。沿考察路线提取了BEDMAP的冰下基岩高程和冰厚度数据 ,以及冰盖表面的坡度与坡向数据。指出应进一步开展该断面以及兰伯特冰川盆地GIS应用研究。     

基于高斯距离加权格网法研究华北地区地壳形变特征

分析了高斯距离加权格网法计算应变率的原理和特征,研制了相应的计算程序,并用于华北地区现今地壳应变特征研究。基于华北地区1999～2009年的GPS观测结果,研究了确定高斯平滑半径的原则和方法,并给出华北地区地壳形变多个应变率物理量的分布。结果表明,基于150 km高斯平滑半径计算水平速度梯度能够较好地分辨华北地区的应变率空间变化特征,首都圈地区地壳以张渤带上的左旋剪切为主,其中唐山-秦皇岛地区同时存在南北向的拉张变形;山西带整体表现为沿断裂的右旋剪切变形,同时兼具“南挤北张”特性;郯庐断裂带整体形变特征不显著,在南北两端局部区域存在较小的剪切形变;鄂尔多斯块体内部比较稳定,北缘以拉张和左旋剪切变形为主,西缘表现为左旋剪切以及EW向的挤压特征,南缘的秦岭块体南北两侧分别具有左旋剪切和右旋剪切的特性,东西两侧分别表现为正断拉张与逆断挤压特征。     

气象中使用统计检验的几个问题

针对目前在气象和气候要素变化分析中,对相关系数和合成分析使用统计检验存在的若干问题进行评论。认为对相关系数等统计检验,不应该称为“信度检验”,其检验水平也不应该称为“95％置信度”,而应该称为“显性检验”,“显性水平为5％”。