STUDY ON CHARACTERISTICS OF GRAVITY VARIATION BEFORE AND AFTER HUTUBI MS6.2 EARTHQUAKE

The Hutubi M_S6.2 earthquake of December 8, 2016 is a pure thrust event in the northern Tianshan thrust fold belt. The earthquake occurred between the Qigu Fault and the Junggar southern margin fault, which are both thrust faults. Based on mobile gravity measurements from 2013 to 2018 in the northern Tianshan, the gravity net adjustment was accomplished using Urumqi absolute gravity observation point as the datum, and the absolute gravity value of surface observation points were obtained. In order to eliminate the seasonal effect on gravity variation, the paper uses the observation data in May per annual as studying objects and obtains the temporal-spatial dynamic evolution images of gravity field differences in the northern Tianshan at different time scales as well as the time series of gravity variation of some points in the adjacent area of the epicenter. The characteristics of regional gravity variation before and after the Hutubi M_S6.2 earthquake on December 8, 2016 and their relations are analyzed systematically. The results show that: 1)The gravity variations in the study area are dramatic in generally, and the contours of gravity variation are consistent with the main faults basically. There was a four-quadrant distribution near the epicenter before the earthquake, and the Hutubi M_S6.2 earthquake occurred near the center of the four-quadrant distribution and at the turn of the gravity variation contour. The three years' cumulative gravity variation before the earthquake and the gravity variation after the earthquake are inversed, and the variation amplitudes are equivalent, suggesting that the M_S6.2 earthquake has released the stress and the energy accumulated before the earthquake. 2)This paper focuses on the analysis of gravity variation at the observation points on both sides of the Junggar southern margin fault near the epicenter. Regional gravity variation and gravity time series show that gravity variations at the same side of the Junggar southern margin fault are basically consistent, however, gravity variations at the different sides of the Junggar southern margin fault are different from each other obviously, indicating the difference of material migration laws in different structural regions. In addition, the strain energy accumulated in the epicenter is basically released after the earthquake, and the area nearby the epicenter tends to be stable. 3)The Hutubi M_S6.2 earthquake occurred near the center of the four-quadrant and at the turn of the high-gradient zone of gravity variation, reflecting the location of strong earthquake is related to the distribution of four-quadrant of regional gravity variation, the high-gradient zone of regional gravity variation and its turn. It has a unique advantage in determining the location of strong earthquake using gravity variation results. The regional spatial-temporal gravity variation before the earthquake is manifested as a systematic evolution process of “steady state→regional gravity anomaly→four-quadrant distribution→earthquake occurring in the reverse process”. Studying the temporal-spatial evolution characteristics of gravity field before and after Hutubi M_S6.2 earthquake has important practical significance for understanding the occurrence law of large earthquakes and capturing the precursory information of earthquakes.     

GRAVITY VARIATION BEFORE AND AFTER THE MS5.4 EARTHQUAKE IN CANGWU IN 2016

Based on the mobile gravity observation data in 2014-2016 in Guangxi and its adjacent areas, this paper systematically analyzed the changes of regional gravity field and its relation to the M_S5.4 Cangwu, Guangxi earthquake on July 31, 2016, and combined with GPS observation data and seismic geological survey results, discussed the temporal and spatial distribution characteristics of the changes of regional gravity field and its mechanism. The results show that:(1) Before and after the M_S5.4 Cangwu earthquake, the gravity anomaly changes near the earthquake area were closely related to the major faults in space, which reflects the crustal deformation and tectonic activities that caused the surface gravity change along the seismogenic fault in the period of 2014-2016; (2) The gravity changes near the epicenter before and after the M_S5.4 Cangwu earthquake showed an evolution process in which the positive gravity anomaly zone changed to the negative gravity anomaly zone, a gravity gradient belt appeared along NNE direction and the earthquake occurred in its reverse change process; (3) The epicenter of the M_S5.4 Cangwu earthquake located both near the gravity gradient belt and in the zero transition zone of the surface strain gradient and the edge of the high maximum shear strain rate area, the observational fact further proved that the dynamic image of gravitational field and deformation field have important instruction significance to the location prediction of strong earthquakes; (4) in recent years, the gravity dynamic change in northwestern Guangxi presented a four-quadrant distribution pattern, and there is the risk of generating earthquake of magnitude about 5 in the center of the quadrants.     

地震前后重力场变化的典型震例分析

系统总结了中国地震局地球物理勘探中心利用流动重力资料进行地震预报的实践经验,研究了8次地震的孕震环境、重力场变化、成因分析和预报实况等。结果表明,一次较强的地震前,处于不同地段的重力测点和测段都要出现上升-下降-上升或下降-上升-下降的较明显的中、短期异常变化;重力场的变化图象较清晰地反映了地震前后重力场的演化过程。     

门源6.4级地震前后加卸载响应比的异常变化

研究了门源６．４级地震前后加卸载响应比的异常变化．门源地震前响应比呈明显高值异常．门源地区响应比的动态变化反映了该区域的地震孕育过程．     

RESEARCH OF SEISMOGENIC STRUCTURE OF THE MENYUAN MS6.4 EARTHQUAKE ON JANUARY 21, 2016 IN LENGLONGLING AREA OF NE TIBETAN PLATEAU

On January 21 2016, an earthquake of M_S6.4 hit the Lenglongling fault zone(LLLFZ)in the NE Tibetan plateau, which has a contrary focal mechanism solution to the Ms 6.4 earthquake occurring in 1986. Fault behaviors of both earthquakes in 1986 and 2016 are also quite different from the left-lateral strike-slip pattern of the Lenglongling fault zone. In order to find out the seismogenic structure of both earthquakes and figure out relationships among the two earthquakes and the LLLFZ, InSAR co-seismic deformation map is constructed by Sentinel -1A data. Moreover, the geological map, remote sensing images, relocation of aftershocks and GPS data are also combined in the research. The InSAR results indicate that the co-seismic deformation fields are distributed on both sides of the branch fault(F₂)on the northwest of the Lenglongling main fault(F₁), where the Earth's surface uplifts like a tent during the 2016 earthquake. The 2016 and 1986 earthquakes occurred on the eastern and western bending segments of the F₂ respectively, where the two parts of the F₂ bend gradually and finally join with the F₁. The intersections between the F₁ and F₂ compose the right-order and left-order alignments in the planar geometry, which lead to the restraining bend and releasing bend because of the left-lateral strike-slip movement, respectively. Therefore, the thrust and normal faults are formed in the two bending positions. In consequence, the focal mechanism solutions of the 2016 and 1986 earthquakes mainly present the compression and tensional behaviors, respectively, both of which also behave as slight strike-slip motion. All results indicate that seismic activity and tectonic deformation of the LLLFZ play important parts in the Qilian-Haiyuan tectonic zone, as well as in the NE Tibetan plateau. The complicated tectonic deformation of NE Tibetan plateau results from the collisions from three different directions between the north Eurasian plate, the east Pacific plate and the southwest Indian plate. The intensive tectonic movement leads to a series of left-lateral strike-slip faults in this region and the tectonic deformation direction rotates clockwise gradually to the east along the Qilian-Haiyuan tectonic zone. The Menyuan earthquake makes it very important to reevaluate the earthquake risk of this region.     

DISCUSSION ON THE SEISMOGENIC STRUCTURE OF THE 2016 MENYUAN M6.4 EARTHQUAKE IN MENYUAN,QINGHAI

On January 21, 2016, a M6.4 earthquake occurred in Menyuan county, Qinghai Province. Its epicenter is located in the Qilian-Hexi Zoulang tectonic zone, which records several moderate-large historical earthquakes. Previous studies on this event are based on geology, remote sensing data and focal mechanism solutions, lacking analysis on its seismogenic structure. In order to study seismogenic fault plane and seismoteconic style of the earthquake, this work uses data of seismic intensity, aftershocks, and geology to address this issue. Furthermore, we calculate Coulomb stress changes imposed by the 1927 Gulang M8 and 1986 Menyuan M6.4 earthquake on the fault plane of the 2016 Menyuan M6.4 earthquake. The results indicate the early two events have posed distinct impacts on two nodal planes:loading or triggering on nodal plane Ⅰ, and unloading or delay on Ⅱ. In some cases such triggering stress is approaching or up to the threshold value of 0.01 MPa. Combining isoseismals, aftershock distribution, geological structure and different Coulomb stress changes aforementioned, the nodal plane Ⅱ of the source model is considered the seismogenic feature. In conjunction with geophysical data, we establish the seismogenic model of the Menyuan earthquake, which is a positive flower structure in a profile, gentle in the upper and steep in the lower, characterized by thrusting in a strike slipping fault system. This is a possible model for thrusting earthquakes generated by strike-slip faults in a compressional tectonic regime.     

THE SEISMOGENIC ENVIRONMENT ANALYSIS OF LUDIAN MS6.5 EARTHQUAKE USING GRAVITY DATA

The main rupture of Ludian M_S6.5 earthquake is directed to the northwest, which occurred in the east of Xianshuihe-Xiaojiang fault zone. The epicenter is in the transitional zone of the Sichuan-Yunnan block and the South China block, where there are many slip and nappe structures. Some controversy still remains on the earthquake tectonic environment. So, Bouguer gravity anomalies calculated by EGM2008 were broken down into 1-5 ranks using the way of Discrete Wavelet Transform(DWT), then we get the lateral heterogeneity in different depths of the crust. The distribution characteristics of Bouguer gravity anomaly are analyzed using measured gravity profile data. We also get its normalized full gradient(NFG)picture, and study the differences between different depths in crust. The results show that: (1)the characteristic of Buoguer gravity anomaly in southwest to northeast is high-low-high between the Lianfeng Fault(LFF)and Zhaotong-Ludian Fault(ZLF). The mainshock and aftershocks are distributed in the middle of the low-value zone, which means that the east moving materials of Qinghai-Tibet plateau broke through the southern section of Lianfeng Fault(LFF), moving along the Baogunao-Xiaohe zone(low-value belt)to the southeast, stopped by the Zhaotong-Ludian Fault(ZLF), and then earthquake occurred.(2)The third-order discrete wavelet transform(DWT)details show that: there is a good consistency between the negative gravity anomaly in upper crust and the distribution of major faults, which reflects that the rupture caused by the movements of the faults in crust has reduced gravity anomaly. There is a NW-trending negative anomaly belt near the epicenter, which may has some relationship to the southward development of the Daliangshan Fault(DLSF). So we speculate that the southward movement of Daliangshan Fault is the main direct force source of Ludian earthquake.(3)In the picture of the fourth-order DWT details, there is an obvious positive gravity anomaly under the epicenter of Ludian earthquake, which confirms the presence of a high-density body in the middle crust. While the fifth-order DWT details show that: A positive anomaly belt is below the epicenter too, which may be caused by mantle material intruding to the lower crust. Tensile force in crust caused by mantle uplift and extrusion-torsion force caused by Indian plate push are the main force source in the tensile and strike slip movement of the Ludian earthquake.(4)The normalized total gradient of Bouguer gravity anomalies of Huili-Ludian-Zhaotong profile shows that: there is obvious ‘deformation’ in the Xiaojiang fault zone which dips to the east and controls the local crust movement. There is a local ‘constant body’ at the bottom of the epicenter. The stable constant body in density has limiting effects to the earthquake rupture, which is the reason that the earthquake rupture' scale in strike and in depth are limited.(5)The ability of earthquake preparation in Zhaotong-Ludian Fault is lower than the Xianshuihe-Xiaojiang fault zone, and the maximum earthquake capacity in this area should be around magnitude 7.     

包头西6.4级地震前地球物理场特征的研究

通过对包头６．４级地震前较大范围前兆资料的研究,认为一个强烈地震在发生之前,围绕震中区不同区域内,将出现具有不同特点的前兆场异常群。其近场（Δ≤１００ｋｍ）的主要特征是具有较长时间尺度的趋势性异常;中场（１００ｋｍ〈Δ≤２００ｋｍ）的主要特征是趋势和突发性异常并重;而远场（Δ〉２００ｋｍ）则一般无趋势性异常,只有急剧变化的短临异常。     

临猗ML5.0级地震前后的重力场变化及预报过程

讨论了气压、地形变、降雨和地下水位变化对重力观测资料的影响,论证了观测资料的可靠程度,在此基础之上结合地质构造,地壳深浅构造等资料,讨论了１９９８年７月１１日临猗地震前后重力场的变化;依据临猗地震前重力场的变化,结合该地区的地震活动性,对这次地震作出了预报。     

2016年门源6.4级地震前连续重力频率特征分析

利用兰州地震台及高台地震台连续重力资料,分析研究了2016年1月21日门源地震前连续重力变化特征。对震前5 d及震后2 d秒采样数据使用快速傅里叶变换（FFT）变换方法进行功率谱分析,发现其频率幅度持续增大,临震前一天减小到正常值后发生地震,同时对震前6 d秒采样数据进行滤波处理,发现在临震前48 h观测到重力扰动信号,周期为4~8 s。门源地震前观测到的频率特征及重力扰动信号分析,可能反映了门源地震前孕震的一个整体过程,这为揭示门源地震产生机理提供了参考。     

1976年唐山地震前的重力位二次微商变化

本文介绍了唐山大地震前在北京通县西集地震台用扭秤观测的重力位二次微商变化。这个资料与唐山地震前的重力异常符合得比较好。扭秤测量精度较高，仪器性能稳定，两秤臂的一致性好，它对唐山地震前观测的重力异常是一个很好的旁证和补充。说明唐山地震前的重力变化，尤其是１９７５年的加速上和临震前的得力下降是可信的。西集台第一次观测到与地震孕育过程有关的重力位二次微商异常，这对研究唐山地震的前兆，特别是重力场的变化原     

2013年芦山MS7.0地震和2008年汶川MS8.0地震井水位同震变化的比较分析

以芦山、汶川地震震中为中心,选择了8个井点,对2个地震引起的井水位同震变化进行了比较,重点分析了3口同震变化反向的井点,即重庆的荣昌华江井、北碚柳荫井和四川的泸沽湖井。为使对比研究结果可靠,进一步收集了2011年3月11日日本9.0级地震、2012年4月11日苏门答腊北部海域8.6级地震引起的井水位同震变化资料。对汶川、芦山地震引起的同震体应变的计算结果显示：汶川地震时,荣昌华江井、北碚柳荫井和泸沽湖井水位同震变化与井点位置处的同震体应变一致,即位于同震体应变压缩区的井水位上升,位于体应变膨胀区的井水位下降;芦山地震时,3口井的水位同震变化与同震体应变不一致,表现出与日本、苏门答腊地震时相似的同震变化特征。     

晋中南地区的重力场变化

本文对晋中南地区的重力复测资料、重力场变化及地震活动特征进行了分析,讨论了重力观测对山西断陷带地震监测的意义。     

RELOCATION OF THE HUTUBI MS6.2 EARTHQUAKE SEQUENCE ON 8 DECEMBER 2016 AND ANALYSIS OF THE SEISMOGENIC STRUCTURE

Based on the digital waveforms of Xinjiang Seismic Network, the Hutubi M_S6.2 earthquake sequence (M_L ≥ 1.0) was relocated precisely by HypoDD.The best double-couple focal mechanisms of the main shock and aftershocks of M_L ≥ 4.0 were determined by the CAP method. We analyzed the characteristics of spatial distribution, focal mechanisms and the seismogenic structure of earthquake sequence. The results show that the main shock is located at 43.775 9°N, 86.363 4°E; the depth of the initial rupture and centriod is about 15.388km and 17km. The earthquake sequence extends unilaterally along NWW direction with an extension length of about 15km and a depth ranging 5~15km. The characteristics of the depth profiles show that the seismogenic fault plane dips northward and the faulting is dominated by thrusting. The nodal planes parameters of the best double-couple focal mechanisms are:strike 292°, dip 62° and rake 80° for nodal plane I, and strike 132°, dip 30° and rake 108° for nodal plane Ⅱ, indicating that the main shock is of thrust faulting. The dip of nodal planeⅠis consistent with the dip of the depth profile, which is inferred to be the fault plane of seismogenic fault of this earthquake. According to the comprehensive analysis of the relocation results, the focal mechanism and geological structure in the source region, it is preliminarily inferred that the seismogenic structure of the Hutubi M_S6.2 earthquake may be a backthrust on the deeper concealed thrust slope at the south of Qigu anticline. The earthquake is a "folding" earthquake taking place under the stress field of Tianshan expanding towards the Junggar Basin.     

POISSON'S RATIO VARIATIONS OF CRUSTAL MEDIA BEFORE AND AFTER XINYUAN-HEJING MS6.6 EARTHQUAKE IN 2012

We analyzed the variation characteristics of Poisson's ratio in crustal media from January 2009 to December 2012 at 11 fixed seismic stations(for station SCH, it is from January 2006 to December 2012)within an epicenter distance of 200km of the Xinyuan-Hejing M_S6.6 earthquake in Xinjiang on June 30, 2012 using the methods of P wave receiver functions, H-κ stacking of receiver functions, and time sliding window, and obtained the following conclusions: (1)The crustal media's Poisson ratio of five stations in an epicenter distance less than 130km showed a significant and long-lasting decline about 2~3 years before Xinyuan-Hejing M_S6.6 earthquake. Taking the crustal Poisson ratio mean value as reference, the decrease ranges between 0.003 and 0.014, the decrease in 4 stations are more than twice the mean error. The variations of the Poisson's ratio in crust are characterized by "V" shape or "double V" shape. Earthquakes occur at the end of the formation of "V" shape. After the occurrence of earthquakes, the Poisson's ratio continues to rise. The earliest initial fall appeared in July 2009 at WUS station which has the minimum epicentral distance(77km). The Poisson ratio of the crustal media of 6 stations with epicentral distance more than 150km fluctuated up and down around the mean value, and there is no significant decline or persistent low value. (2)We analyzed the arrival-time variations of the quasi-repetitive receiver functions Ps converted wave(t_Ps)of the 3 stations WUS, SCH and XNY and found that the travel times of Ps converted waves became smaller in the crust before the earthquake and increased after the earthquake. (3)Through the comprehensive analysis on the descending process, decline ranges, variations process, duration of Poisson' ratio, the Ps converted waves arrival time variations, the original time of earthquake, and the number of stations, it is inferred that the cause for Poisson's ratio anomalous variations is the change of physical properties of crustal media in the process of earthquake preparation and occurrence. Since the variation characteristics of crustal media may be related to the earthquake magnitude, the size of seismogenic area, the medium properties under stations, and the focal distance, whether the medium variation characteristics exist before and after Xinyuan-Hejing M_S6.6 earthquake will need more earthquake cases analyses. (4)The H-κ stacking of receiver functions is used to calculate the velocity ratio. Because P-wave velocity is given, this method can only be applied when the Ps converted wave velocity of Moho surface of receiver functions changes before an earthquake. With the application of receiver functions to the analysis of more earthquake cases, we can gain more insights into the variation of crustal medium parameters during the seismogenic process. This observation indicates that the receiver function method may become a new approach to detect the Poisson's ratio change of the crustal media before strong earthquake under the condition of high seismic network density.     

STUDY ON THE SEISMOGENIC FAULT CHARACTERSTICS OF 2016 MW5.9 MENYUAN EARTHQUAKE BASED ON Sentinel -1A DATA

In this paper, we processed and analyzed the Sentinel-1A data by "two-pass" method and acquired the surface deformation fields of Menyuan earthquake. The results show the deformation occurred mainly in the south wall of fault, where uplift deformation is dominant. The uplift deformation is significantly larger than the subsidence and the maximum uplift of ascending and descending in the LOS is 6cm, 8cm respectively. Meanwhile, based on the Okada model, we use the ascending and descending passes data as constraints to invert jointly the fault distribution and source parameters through constructing fault model of different dip directions. The optimum fault parameters are:The dip is 43°, the strike is 128°with the mean rake of 85°. The maximum slip is about 0.27m. The inverted seismic moment M₀ is 1.13×10¹⁸N·m, and the moment magnitude M_W is 5.9. The SW-dipping Minyue-Damaying Fault is possibly the seismogenic fault, based on the comprehensive analysis of the focal mechanisms, aftershocks relocation results and the regional tectonic background. The focus property is dominated by thrust movement with a small amount of dextral strike-slip component. The earthquake is the result of local stress adjustment nearby the Lenglongling Fault under the background of northeastward push and growth of Tibet Plateau.     

1990年共和7.0级地震水化学前兆场特征

The features of Hydrochemical precursory field before the Gonghe M_S7.0 earthquake on April 20,1990 are studied.The results show that proportion of hydrochemistry observation sites in which trend anomalies appear to all sites in focal region of the earthquake is higher than that far away from the focal region and the proportion reduces gradually as the epicentral distance increases.In the focal area the trend anomalies in which radon content in groundwater goes up gradually emerge earlier than far away from foci.The formation of hydrochemical precursory field has something to do with regional stress field,but has closer relation with focal stress field.Focal stress field may play a main role in the precursory process.     

2009年阿合奇M_S5.5、阿图什M_S5.0地震前后重力场变化分析

计算分析了近2期流动重力资料,对2009年4月19日阿合奇发生的MS5.5地震和2009年4月22日阿图什发生的MS5.0地震进行了初步分析。结果表明,这两次地震的孕育与发生,都伴随着重力异常的出现,且地震发生在重力场变化的零线附近和重力场反向回跳的转折线附近。     

2014年云南鲁甸6.5级地震前重力变化特征与3维反演

利用西昌测网2012-2014年6期流动重力观测资料,分析了区域重力场差分动态变化和累积动态变化特征,以及近震源区测点的时序变化曲线,认为鲁甸6.5级地震前2年并没有出现明显的趋势性累积异常,而可能的前兆信号与断层滑动面两侧的测点2014-03-2014-06时段的重力点值差异性变化相关.针对此异常采用等效源模型反演至震源深度范围内的密度值变化在10-5g/cm3量级,推测该质量源在短期内的变化可能与地壳内部介质孔隙流体的充填或运移相关,在短期构造运动等动力条件下,孕震区可能由于流体物质运移、充填并诱发断层滑动,导致地震的发生.     

1986年门源地震(MS 6.4)过程地形变演化特征及块体模型解析

分析了1986年门源6．4级地震震源周边地区不同尺度跨断层形变在地震过程中的演化特征,依据非连续变形数值分析方法(DDA)模拟了震源破裂过程引起的地表垂直位移速率变化．分析表明,地表跨越断层的剖面形变变化是垂直和水平栽荷共同作用的结果,而垂直力源的影响可能是引起本次地震形变变化的主要因素．这与地质构造分析和震源机制结果是吻合的．