上海地震台阵的地震定位方法

上海地震台阵数据处理软件系统的地震定位方法,采用台阵的聚束方法得到地震方位角和视慢度,根据统计得到的视慢度——震中距表推算震中距。并结合了地震台网的定位方法,由单台记录的各类主要震相从JB走时表得到震中距,然后进行地震定位。该定位方法可对近震、远震进行定位处理,并由深震相得到震源深度。     

核试验监测的地震学研究综述

侦察探测核爆事件,评估核爆信息等核试验监测工作,对评估核爆威力、研究核爆炸现象以及各种破坏因素的变化规律有着重要意义.地震学提供了核试验监测中关键的研究手段.核试验地震学监测研究主要包括核试验识别、定位、埋深、震级和当量的估计等内容.根据近半个世纪以来核试验方式的变化和技术的发展,介绍和总结了核试验监测中相关研究进展状况,主要包括事件定位、判别、震级测定、当量和埋深估计等地震学观测、理论和方法.     

近台资料对近震相对定位算法的影响

从地震定位的基本原理出发，分析认为近台资料可用于双差定位算法，但不能用于其他相对定位方法.在此基础上，通过分析和模拟计算表明，应用近台资料时，双差法有可能给出一定精度的地震绝对位置.但仅采用近台资料时，相对位置的误差会比使用远台资料时有所增加.近台和远台资料的联合使用，有利于得到较为精确的定位结果.当震源深度远小于震中距时，如果没有深度震相的参与，只能得到误差较小的震中相对分布，深度的相对位置仍有较大的误差.对2003年新疆伽师地震余震序列中部分余震的重新定位试验，验证了近台资料对双差定位算法的上述影响.     

A Multidisciplinary Study of the DPRK Nuclear Tests

The Democratic People Republic of Korea announced two underground nuclear tests carried out in their territory respectively on October 9th, 2006 and May 25th, 2009. The scarce information on the precise location and the size of those explosions has stimulated various kinds of studies, mostly based on seismological observations, by several national agencies concerned with the Nuclear Test Ban Treaty verification. We analysed the available seismological data collected through a global high-quality network for the two tests. After picking up the arrival times at the various stations, a standard location program has been applied to the observed data. If we use all the available data for each single event, due to the different magnitude and different number of available stations, the locations appear quite different. On the contrary, if we use only the common stations, they happen to be only few km apart from each other and within their respective error ellipses. A more accurate relative location has been carried out by the application of algorithms such as double difference joint hypocenter determination (DDJHD) and waveform alignment. The epicentral distance between the two events obtained by these methods is 2 km, with the 2006 event shifted to the ESE with respect to that of 2009. We then used a dataset of VHR TerraSAR-X satellite images to detect possible surface effects of the underground tests. This is the first ever case where these highly performing SAR data have been used to such aim. We applied InSAR processing technique to fully exploit the capabilities of SAR data to measure very short displacements over large areas. Two interferograms have been computed, one co-event and one post-event, to remove possible residual topographic signals. A clear displacement pattern has been highlighted over a mountainous area within the investigated region, measuring a maximum displacement of about 45 mm overall the relief. Hypothesizing that the 2009 nuclear test had been carried out close to the area where the displacement has been observed through the DInSAR technique, its relation with the epicenter location obtained through seismological processing has been discussed as a possible alternative hypothesis with respect to the preferred solutions reported by the nuclear explosion database (NEDB). The distance of about 10 km between the two places can be considered acceptable in light of the possible systematic location shifts commonly observed in the seismological practice over a global scale. The difference between the m _b magnitudes of the two tests could reflect differences in geological conditions of the two test sites, even if the yield of the two explosions had been the same.     

Sources of Error Model and Progress Metrics for Acoustic/Infrasonic Analysis: Location Estimation

How well can we locate events using infrasound? This question has obvious implications for the use of infrasound within the context of nuclear explosion monitoring, and can be used to inform decision makers on the capability and limitations of infrasound as a sensing modality. This paper attempts to answer this question in the context of regional networks by quantifying current capability and estimating future capability using an example regional network in Utah. This example is contrasted with a sparse network over a large geographical region (representative of the IMS network). As a metric, we utilize the location precision, a measure of the total geographic area in which an event may occur at a 95 % confidence level. Our results highlight the relative importance of backazimuth and arrival time constraints under different scenarios (dense vs. sparse networks), and quantify the precision capability of the Utah network under different scenarios. The final section of this paper outlines the research and development required to achieve the estimated future location precision capability.     

用CAP方法反演2010年6月5日山西阳曲M_S4.6地震震源机制解

利用CAP方法反演了2010年6月5日阳曲MS4.6地震震源机制解,得到震级MW为4.5,节面I走向213°、倾角47°、滑动角-161°,节面II走向109°,倾角76°,滑动角-44°,属于倾滑型;精确定位显示震中处于石岭关隆起区,CAP反演和精定位结果推断本次地震的震源深度为17～20km。震源机制解节面参数与震中附近的山根底断裂和系舟山西麓断裂产状存在差异,这两条断裂不是阳曲地震的发震断裂,由于现场野外地质考察未发现地表断裂,不排除本次地震为隐伏断层活动的结果。     

Use of GSETT-3 gamma data in the Slowness-Azimuth Calibration of IMS primary arrays at regional distances

For the purpose of verifying compliance with the CTBT seismic monitoring is one of the four techniques used by the IDC. In order to improve the accuracy of the automatic and the reviewed bulletin epicenter locations the IDC uses SASC for the IMS seismic stations. SASC determination is a straightforward calculation done by comparing for selected events the azimuth and slowness from the waveform processing using array techniques to the theoretical values based on the event locations and the velocity model.The main problem, however, is to build a set of reference events, whose locations are accurate enough and not based on information from the stations to be calibrated. A reference event list assumed to meet this requirement is the Gamma bulletin, which was collected since 1993 and was compiled in the framework of the GSETT-3. In this work calculation of SASC for regional to teleseismic distances (up to 30 degrees) was performed for 11 IMS primary arrays. The calculation was done using Pg, Pn, P, Sg, Sn, and S phases based on the detection list obtained from the pIDC and the Gamma bulletin for 6 years (1993–1999). The number of Gamma events varies from several hundreds for some arrays (BRAR) to several tens of thousands for others (i.e. ARCES, ILAR). Due to the fact that the Gamma bulletin is purely voluntary, the coverage is non-uniform both in time and in space and the location accuracy is non-uniform. This drawback can be overcome by encouraging signatory states to submit quality Gamma bulletin data to the IDC. The work presented here can be used as a routine procedure for improving IMS array performance, especially at regional distances.     

2015年8月12日天津化学爆炸产生的多模式面波分析及其应用研究

2015年8月12日天津滨海新区发生的强烈化学品爆炸造成了巨大的经济损失和社会影响.天津爆炸产生了清晰的大振幅面波信号,分析结果表明这组信号由基阶和高阶面波组成,可以追踪到约135 km外的远处台站.利用这组面波信号分别开展了以下研究：（1）利用附近三个台站记录的四个单频基阶Rayleigh波信号对爆破事件的绝对位置进行了网格搜索,结果与利用GPS测量的位置相差仅0.498 km;（2）分别利用网格搜索和主事件定位法,对两次子事件的相对位置进行了确定,距离约75 m左右,与前人研究结果吻合;（3）从面波记录中测量到36条基阶Rayleigh波、49条第一高阶Rayleigh波、9条基阶Love波和29条第一高阶Love波的频散曲线,并进一步反演获得研究区域地下4 km内的S波速度结构.反演结果显示地表处S波速度低至0.375 km·s^-1,在小于1 km的浅地表速度梯度较大,符合典型的盆地结构特征.本文的研究结果为类似爆炸等突发事件快速定位提供了新的思路,有助于灾后救援的迅速展开;同时得到天津滨海新区及周边浅层精细的速度结构,对于地震灾害评估有很大帮助.     

Relocation of Wyoming Mine Production Blasts Using Calibration Explosions

—?An important requirement for a comprehensive seismic monitoring system is the capability to accurately locate small seismic events worldwide. Accurate event location can improve the probability of determining whether or not a small event, recorded predominantly by local and regional stations, is a nuclear explosion. For those portions of the earth where crustal velocities are not well established, reference event calibration techniques offer a method of increased locational accuracy and reduced locational bias.¶In this study, data from a set of mining events with good ground-truth data in the Powder River Basin region of eastern Wyoming are used to investigate the potential of event calibration techniques in the area. Results of this study are compared with locations published in the prototype International Data Center's Reviewed Event Bulletin (REB). A Joint Hypocenter Determination (JHD) method was applied to a s et of 23 events. Four of those events with superior ground-truth control (mining company report or Global Positioning System data) were used as JHD reference events. Nineteen (83%) of the solutions converged and the resulting set of station-phase travel-time corrections from the JHD results was then tested. When those travel-time corrections were applied individually to the four events with good ground-truth control, the average locational error reduced the original REB location error from 16.1?km to 5.7?km (65% improvement). The JHD locations indicated reduced locational bias and all of the individual error ellipses enclosed the actual known event locations.¶Given a set of well-recorded calibration events, it appears that the JHD methodology is a viable technique for improving locational accuracy of future small events where the location depends on arrival times from predominantly local and/or regional stations. In this specific case, the International Associ ation of Seismology and the Physics of the Earth's Interior (IASPEI) travel-time tables, coupled with JHD-derived travel-time corrections, may obviate the need for an accurately known regional velocity structure in the Powder River Basin region.     

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.     

公元前179年“齐楚地震”考证与发震构造讨论

对公元前179年(汉文帝前元元年)"齐楚地震"的基本参数一直有不同认识,属于疑难历史地震。在吸收西汉政区研究成果的基础上,深入分析了地震史料隐含的约束条件,对这次地震震中位置进行了考证。提出了以齐楚边界和国都连线为基点,结合史料给出的有关约束条件,综合确定震中位置的技术方案。同时,依据地震地质调查成果,从发震构造角度论证了震中位置的合理性。1)鉴于西汉时期政权更迭频繁,通过追踪西汉初期(尤其是汉文帝前元元年前后)齐、楚政区变化情况,尽力复原当时的行政区划和齐楚两国边界分布。2)认真分析史料记载的"齐楚地震,二十九山同日崩,大水溃出"所指示的地区范围,对比历史上其他震例造成的破坏现象的空间分布以及破坏程度,为震中位置和震级确定提供重要约束条件。经文献考证认为,这次地震可定为公元前179年6月6日平邑南(35.2°N,117.6°E)7级地震。3)结合卫星影像判读和野外地震地质调查资料,分析研究了"齐楚地震"所在地区的地质构造特点和活动断裂发育情况。在排除该地震事件属于郯庐断裂带地震事件后,通过对比分析郯庐断裂带西侧发育的各条NW向断裂的构造地貌以及剖面表现,认为苍尼断裂可能是"齐楚地震"的发震断裂。在震中附近,苍尼断裂控制着白彦断陷盆地的发育,也是断层地貌最为清楚的地方,断错的地层最新,断层运动量也最大。     

基于PHP的地震台站综合管理系统

天津市地震局滨海地震台测震、地球物理等观测手段众多,数据处理量大,处理软件繁多。为此,利用PHP开发网页化的综合管理系统,将工作内容进行分类整理,将台站工作展示、地震观测工作辅助、工作提示等功能整合在系统中,从而减少数据处理错误,提高数据观测质量及工作效率。该平台具有扩展性,可随时添加新的功能模块,适应未来工作变化。     

汶川Ms8.0地震前InSAR垂直形变场变化特征研究

利用D-InSAR技术,选用汶川地震前日本ALOS/PALSAR数据,提取了汶川地震前4个条带的形变场,形变场覆盖了映秀镇、茂县、北川县、平武和青川县.结果显示,汶川地震前形变场沿断层呈“凸”状对称形变分布特征,沿发震断层附近出现隆起形变,而远离发震断层两盘均出现下沉,且随着远离断层距离的增加,沉降幅度逐步增加.在断层两侧附近的平武、北川、安县、江油一带出现了大面积形变隆起,幅度为5～10 cm,且断层西侧抬升范围大于东侧.在映秀镇和汶川地震震中附近,沿断层两侧呈现出小范围斑块状隆起,范围在10～15 cm.该隆起区域与汶川震中位置和破裂最强烈段落基本一致.在远离断层的两盘区域出现沉降,绵阳至成都一带沉降范围在0～ -10 cm之间,松潘至文县一带沉降范围在-5～ -10 cm之间.研究表明对于处于闭锁阶段的逆冲断层,震前的垂直变形可能是其主要变形特征.虽然获得的震前形变变化可能存在5 cm的DEM和对流层大气延迟造成的误差,但这种震前垂直形变场变化趋势仍然存在,可能为地震监测预报提供科学依据.     

The use of teleseismic P-wave complexity for seismic event screening – results determined from GRF and GERESS array data

Seismic data recorded at the broad-band teleseismic GRF array and theshort-period regional GERESS array, which is a designated IMS primarystation, are analyzed to determine the effectiveness of teleseismic P-wave complexity for the purpose of seismic event screening within theframework of Comprehensive Nuclear-Test-Ban Treaty verification. For theGRF array, seismic waveform data from nearly 200 nuclear explosions havebeen recorded since its installation in the late 1970's, which were studiedalong with several thousand earthquakes from the last few years.Additionally, we investigated teleseismic P wave complexity for a similarnumber of earthquakes recorded at GERESS. However, owing to itsoperation starting in 1991, only a limited number of nuclear explosionseismograms are available for study.For nuclear explosions, complexity does not exceed levels of 0.3 except fora number of events from the Nevada Test Site recorded only at the GRFarray and located at a large distance where PcP may interfere with the initialP wavelet. Since all events with complexity at GRF larger than 0.3 areexclusively located on Pahute Mesa within the Nevada Test Site,near-source geology or topography must play a dominant role for theseincreased complexity values, while PcP may not contribute significantly tothe high-frequency energy measured by the complexity parameter.Although many earthquakes show complexities below this level, for morethan 25% of the earthquakes investigated the complexities determined arelarger than 0.7, thus showing distinctly larger values than nuclearexplosions. Therefore, this percentage may be screened as earthquakes fromall seismic events detected. As currently only about half of the eventsdetected by the global IMS network are screened out based on focal depthand the m _b :M _s criterion, teleseismic P-wavecomplexity may contribute significantly to the task of seismic eventscreening.     

Improving Regional Seismic Event Location in China

—?In an effort to improve our ability to locate seismic events in China using only regional data, we have developed empirical propagation path corrections and applied such corrections using traditional location routines. Thus far, we have concentrated on corrections to observed P arrival times for crustal events using travel-time observations available from the USGS Earthquake Data Reports, the International Seismic Centre Bulletin, the preliminary International Data Center Reviewed Event Bulletin, and our own travel-time picks from regional data. Location ground truth for events used in this study ranges from 25?km for well-located teleseimic events, down to 2?km for nuclear explosions located using satellite imagery. We also use eight events for which depth is constrained using several waveform methods. We relocate events using the EvLoc algorithm from a region encompassing much of China (latitude 20°–55°N; longitude 65°–115°E). We observe that travel-time residuals exhibit a distance-dependent bias using IASPEI91 as our base model. To remedy this bias, we have developed a new 1-D model for China, which removes a significant portion of the distance bias. For individual stations having sufficient P-wave residual data, we produce a map of the regional travel-time residuals from all well-located teleseismic events. Residuals are used only if they are smaller than 10?s in absolute value and if the seismic event is located with accuracy better than 25?km. From the residual data, correction surfaces are constructed using modified Bayesian kriging. Modified Bayesian kriging offers us the advantage of providing well-behaved interpolants and their errors, but requires that we have adequate error estimates associated with the travel-time residuals from which they are constructed. For our P-wave residual error estimate, we use the sum of measurement and modeling errors, where measurement error is based on signal-to-noise ratios when available, and on the published catalog estimate otherwise. Our modeling error originates from the variance of travel-time residuals for our 1-D China model. We calculate propagation path correction surfaces for 74 stations in and around China, including six stations from the International Monitoring System. The statistical significance of each correction surface is evaluated using a cross-validation technique. We show relocation results for nuclear tests from the Balapan and Lop Nor test sites, and for earthquakes located using interferometric synthetic aperture radar. These examples show that the use of propagation path correction surfaces in regional relocations eliminates distance bias in the residual curves and significantly improves the accuracy and precision of seismic event locations.     

张北-尚义地震的地震构造环境与宏观破坏特征

阐述了张北－尚义地震的区域地震构造环境,结合本次地震的宏观烈度分布、震源机制解和震区卫星影像的线性构造解释等资料,讨论了地震的孕震构造和可能的发震断层,认为本次地震是张家口－蓬莱断裂带北西端最新活动的结果,北西向的张家口－蓬莱断裂带未来的地震活动趋势应引起重视     

A multi-channel approach for automatic microseismic event association using RANSAC-based Arrival Time Event Clustering (RATEC)

In the presence of background noise, arrival times picked from a surface microseismic data set usually include a number of false picks that can lead to uncertainty in location estimation. To eliminate false picks and improve the accuracy of location estimates, we develop an association algorithm termed RANSAC-based Arrival Time Event Clustering (RATEC) that clusters picked arrival times into event groups based on random sampling and fitting moveout curves that approximate hyperbolas. Arrival times far from the fitted hyperbolas are classified as false picks and removed from the data set prior to location estimation. Simulations of synthetic data for a 1-D linear array show that RATEC is robust under different noise conditions and generally applicable to various types of subsurface structures. By generalizing the underlying moveout model, RATEC is extended to the case of a 2-D surface monitoring array. The effectiveness of event location for the 2-D case is demonstrated using a data set collected by the 5200-element dense Long Beach array. The obtained results suggest that RATEC is effective in removing false picks and hence can be used for phase association before location estimates.     

Seismic phase picking in China Seismic Array using a deep convolutional neuron network*

Seismic phase picking is the preliminary work of earthquake location and body-wave travel time tomography. Manual picking is considered as the most accurate way to access the arrival times but time consuming. Many automatic picking methods were proposed in the past decades, but their precisions are not as high as human experts especially for events with low ratio of signal to noise and later arrivals. As the increasing deployment of large seismic array, the existing methods can not meet the requirements of quick and accurate phase picking. In this study, we applied a phase picking algorithm developed on the base of deep convolutional neuron network (PickNet) to pick seismic phase arrivals in ChinArray-Phase III. The comparison of picking error of PickNet and the traditional method shows that PickNet is capable of picking more precise phases and can be applied in a large dense array. The raw picked travel-time data shows a large variation deviated from the traveltime curves. The absolute location residual is a key criteria for travel-time data selection. Besides, we proposed a flowchart to determine the accurate location of the single-station earthquake via dense seismic array and phase arrival picked by PickNet. This research expands the phase arrival dataset and improves the location accuracy of single-station earthquake.     

Rapid topographic change measured by high-resolution satellite radar at Soufriere Hills Volcano,Montserrat, 2008–2010

High-resolution satellite radar observations of erupting volcanoes can yield valuable information on rapidly changing deposits and geomorphology. Using the TerraSAR-X (TSX) radar with a spatial resolution of about 2 m and a repeat interval of 11 days, we show how a variety of techniques were used to record some of the eruptive history of the Soufriere Hills Volcano, Montserrat between July 2008 and February 2010. After a 15-month pause in lava dome growth, a vulcanian explosion occurred on 28 July 2008 from a vent that was hidden by dense cloud. We were able to show the civil authorities using TSX difference images of surface roughness change that this explosion had not disrupted the dome sufficiently to warrant continuation of a previous, precautionary evacuation. Change difference images also proved to be valuable in mapping new pyroclastic flow deposits: the valley-occupying block-and-ash component tended to increase backscatter and the marginal surge deposits to reduce it, with the pattern reversing after the event due to erosion and deposition. By comparing east- and west-looking images acquired 12 h apart, the deposition of some individual pyroclastic flows can be inferred from change differences. Some of the narrow upper sections of valleys draining the volcano received many tens of metres of rockfall and pyroclastic flow deposits over periods of a few weeks. By measuring the changing radar shadows cast by these valleys in TSX images the changing depth of infill by deposits could be estimated. In addition to using the amplitude data from the radar images we also used their phase information within the InSAR technique to calculate the topography during a period of no surface activity. This enabled areas of transient topography, crucial for directing future flows, to be captured.     

高频GPS和强震仪数据在汶川地震参数快速确定中的初步应用

高频GPS可以实时获取地表位移数据,在地震学中有十分重要的现实应用,比如快速获取震中、震级、地震烈度甚至震源破裂过程.本文以汶川地震为例,首先利用近场7个GPS台站数据反演震中位置,由于高频GPS和测震学确认的震相不一致,两种震中结果相距约15.7 km.然后对高频GPS和强震动数据进行了比较分析,我们的统计结果表明,尽管由于工作原理不同,高频GPS数据中的地震动峰值与强震记录相比存在明显差异,但是高频GPS记录的PGA、PGV和PGD同样可以作为计算地震烈度的指标.进而,使用SMBLOC程序对强震记录进行事后的基线偏移校正,得到与实时高频GPS精度相当的地表位移序列.最后,采用移动平均窗口对这些位移数据作平滑,基于最速下降法和OKADA模型,对汶川地震断层破裂的过程进行了回溯性准实时反演.结果表明,汶川地震主断层由西南向东北方向破裂,以14∶28∶04为基准,在震后20 s提供初始震级M_W7.0,震后70 s震级稳定在M_W7.8,但断层仍在破裂,在震后159 s根据位移波形判断事件基本结束.研究表明,实时地表位移数据可以快速准确获取强震震级和破裂方向,从而使得高频GPS将对现有地震预警系统提供很好的补充.