Blockly earthquake transformer:A deep learning platform for custom phase picking

Deep-learning(DL)algorithms are increasingly used for routine seismic data processing tasks,including seismic event detection and phase arrival picking.Despite many examples of the remarkable performance of existing(i.e.,pre-trained)deep-learning detector/picker models,there are still some cases where the direct applications of such models do not generalize well.In such cases,substantial effort is required to improve the performance by either developing a new model or fine-tuning an existing one.To address this challenge,we present Blockly Earthquake Transformer(BET),a deep-learning platform for efficient customization of deep-learning phase pickers.BET implements Earthquake Transformer as its baseline model,and offers transfer learning and fine-tuning extensions.BET provides an interactive dashboard to customize a model based on a particular dataset.Once the parameters are specified,BET executes the corresponding phase-picking task without direct user interaction with the base code.Within the transfer-learning module,BET extends the application of a deep-learning P and S phase picker to more specific phases(e.g.,Pn,Pg,Sn and Sg phases).In the fine-tuning module,the model performance is enhanced by customizing the model architecture.This no-code platform is designed to quickly deploy reusable workflows,build customized models,visualize training processes,and produce publishable figures in a lightweight,interactive,and open-source Python toolbox.

Determination of focal depth by two waveform-based methods:A case study for the 2008 Panzhihua earthquake

With the 2008 Ms6.1 Panzhihua earthquake as a case study, we demonstrate that the focal depth of the main shock can be well constrained with two approaches： （1） using the depth phase sPL and （2） using full waveform inversion of local and teleseismic data. We also show that focal depths can be well constrained using the depth phase sPL with single broadband seismic station. Our study indicates that the main shock is located at a depth of ii kin, much shallower than those from other studies, confirming that the earthquake occurs in upper crust. Aftershocks are located in the depth range of 11 16 kin, which is consistent with a ruptured near vertical fault whose width is about 10 km, as expected for an Ms6.1 earthquake.

Recent advances in earthquake monitoring II: Emergence of next-generation intelligent systems

Seismic data processing techniques,together with seismic instrumentation,determine our earthquake monitoring capability and the quality of resulting earthquake catalogs.This paper is intended to review the improvement of earthquake monitoring capability from the perspective of data processing.Over the past two decades,seismologists have made considerable advancements in seismic data processing,partly thanks to the significant development of computational power,signal processing,and machine learning techniques.In particular,wide application of template matching and increasing use of deep learning significantly enhance our capability to extract signals of small earthquakes from noisy data.Relative location techniques provide a critical tool to elucidate fault geometries and seismicity migration patterns at unprecedented resolution.These techniques are becoming standard,leading to emerging intelligent software systems for next-generation earthquake monitoring.Prospective improvements in future research must consider the urgent needs in highly generalizable detection algorithms(for both permanent and temporary deployments)and in emergency real-time monitoring of ongoing sequences(e.g.,aftershock and induced seismicity sequences).We believe that the maturing of intelligent and high-resolution processing systems could transform traditional earthquake monitoring workflows and eventually liberate seismologists from laborious catalog construction tasks.

On the generation of deep focus earthquakes in subduction zones

Following a quasi-dynamic scheme proposed by Minear and Toksoz (1970), thermal structures of subduction zonesfor different models by finite element method (FEM) were calculated. Density distribution and p-wave anomaly ofsubduction zones were calculated at the same time. Comparing with seismological evidences and results of laboratories. it is proposed that earthquakes occurred below 400 km depth are probably controlled by anti-crackmechanism.

Coupled large earthquakes in the Baikal rift system: Response to bifurcations in nonlinear resonance hysteresis

The current lithospheric geodynamics and tectonophysics in the Baikal rift are discussed in terms of a nonlinear oscillator with dissipation.The nonlinear oscillator model is applicable to the area because stress change shows up as quasi-periodic inharmonic oscillations at rifting attractor structures （RAS）.The model is consistent with the space-time patterns of regional seismicity in which coupled large earthquakes,proximal in time but distant in space,may be a response to bifurcations in nonlinear resonance hysteresis in a system of three oscillators corresponding to the rifting attractors.The space-time distribution of coupled MLH ＞ 5.5 events has been stable for the period of instrumental seismicity,with the largest events occurring in pairs,one shortly after another,on two ends of the rift system and with couples of smaller events in the central part of the rift.The event couples appear as peaks of earthquake ‘migration＇ rate with an approximately decadal periodicity.Thus the energy accumulated at RAS is released in coupled large events by the mechanism of nonlinear oscillators with dissipation.The new knowledge,with special focus on space-time rifting attractors and bifurcations in a system of nonlinear resonance hysteresis,may be of theoretical and practical value for earthquake prediction issues.Extrapolation of the results into the nearest future indicates the probability of such a bifurcation in the region,i.e.,there is growing risk of a pending M ≈ 7 coupled event to happen within a few years.

Locating the Small 1999 Frenchman Flat, Nevada Earthquake with InSAR Stacking

Due to high interferometric coherence in the Nevada region,Interferometric Synthetic Aperture Radar(InSAR)phase stacking is capable of mapping coseismic signals from the 27 January 1999,M w 4.8 Frenchman Flat earthquake.This is one of the smallest earthquakes yet studied using InSAR with line-of-sight displacements as small as~1.5 cm.Modelling the event as dislocation in an elastic half space suggests that the fault centroid was located at(115.96°W,36.81°N)with a precision of 0.2~0.3 km(1σ)at a depth of 3.4±0.2 km.Despite the dense local seismic network in southern Nevada,differences as large as 2~5 km were observed between our InSAR earthquake location and those estimated from seismic data.The InSAR-derived magnitude appeared to be greater than that from seismic data,which is consistent with other studies,and believed to be due to the relatively long time interval of InSAR data.

The April 24, 2013 Changning M_s4.8 earthquake: a felt earthquake that occurred in Paleozoic sediment

The dense broadband seismic network provides more high-quality waveform that is helpful to improve constraint focal depth of shallow earthquake. Many shal- low earthquakes occurring in sediment were regarded as induced events. In Sichuan basin, gas industry and salt mining are dependent on fluid injection technique that triggers microseismicity. We method with regional records adopted waveform inversion to obtain focal mechanism of an Ms4.8 earthquake at Changning. The result suggested that the Changning earthquake occurred at a ESE thrust fault, and its focal depth was about 3 km. The depth phases including teleseismic pP phase and regional sPL phase shows that the focal depth is about 2 km. The strong, short- period surface wave suggests that this event is a very shallow earthquake. The amplitude ratio between Rayleigh wave and direct S wave was also used to estimate the source depth of the mainshock. The focal depth （2-4 km） is far less than the depth of the sedimentary layer thickness （6-8 km） in epicentral region. It is close to the depth of fluid injection of salt mining, which may imply that this event was triggered by the industrial activity.

The solar and lunar effect of earthquake duration and distribution

Phase folding algorithms are conventionally used in periodicity analyses using X-ray astronomy pulsar. These allow for accurate identification of the cycle and phase characteristics of the physical parameters of the periodic variation. Although periodic variations in earth- quake activity have long been studied, this paper is the first to apply the phase folding algorithm to the analysis of shallow （〈70 km） seismic data for the period 1973-2010. The goal is to study the phase distribution characteristics of earthquake frequencies and we see a connection between earthquake occurrence and solar and lunar cycles. First, the rotation of the Sun may play a significant role in impacting on the occurrence time of earthquakes with magnitudes of less than 6.0. This may be especially pertinent for earth- quakes with magnitudes between 5.0 and 6.0, when the modulation ratio reaches 12 %. The Moon＇s gravity, which is generally thought to have the greatest influence on the global environment, may actually play less of a role on earthquake timing than the rotation of the Sun. Second, when we consider the world to be divided into 72 local regions based on latitude and longitude, we can see that there are more than a dozen regions with significant non- uniform distributions of earthquake occurrence time. In these regions, the ratio of χ2 to the number of degrees of freedom far exceeds five. As a result, we posit that some factors associated with the Sun-Earth-Moon relationship may trigger earthquake activity under certain temporal and spatial conditions.

Numerical simulation of precursors of Gonghe M7.0 earthquake in Qinghai Province

Based on a simplified media model of Gonghe area, the precursor characteristics of Gonghe M 7.0 earthquake in 1990 are simulated in this paper by using the constitutive relationship of binary medium (solid and water). The results show that the simulated state distribution and extension variation of media are identical with the spatial and temporal distribution of reliable anomalies before Gonghe earthquake. The study also suggests that the development of the Gonghe earthquake has experienced a series of processes such as elastic deformation of large scale, early nonelastic dilatation, strain softening and elastic recovery in the neighbour region of Xining, nonelasscale, early nonelastic dilatation, strain softening and elastic recovery in the neighbour region of nonelastic dilatation of high dense block near seismic source and earthquake occurrence. Therefore, it can be concluded that this earthquake development is neither a simple process for fissures developing and linking up with one another, nor a process of nonelastic volume expanion and water flowing into medium around seismic source, conrary, it is a complicated medium saate changing process, and the style of such changing depends on geological structure environment of seismic source and its neighbour region. It is considered that different earthquakes are accompanied by different geological conditions, the spatial and temporal behavior of their precusors are certainly different.

A Study on the Early Warning Time of Strong Earthquakes in Taiwan

The paper collects the records by the Fujian Digital Seismic Network of 40 shallow earthquakes in Taiwan with M_S≥5.0 from 1999 to 2013,analyzes the seismic phase(Pn,Sn phase)characteristics and travel-time rules,determines travel-time models and develops a seismic phase travel-time equation based on the two-step fitting algorithm.With the deduction of processing time and network delay time,this method can provide an accurate estimation of early warning time of Taiwan earthquakes for the Fujian region,and has been officially employed in the earthquake early warning system of Fujian Province.

Rapid Determination of the Centroid Moment Tensor of the 2014 M_S7.3 Yutian,Xinjiang Earthquake

We successfully employ an automatic centroid moment tensor(CMT) inversion system to infer the CMT solutions of the February 12,2014 MS7.3 Yutian,Xinjiang earthquake using near-field seismic waveforms(4° < △ < 12°) observed by the virtual China seismic networks,which have been recently set up.The results indicate that this event occurred on a rupture plane(strike 243°,dip 70°,and rake-18°),showing left-lateral strike-slip faulting with a minor normal-faulting component.The centroid in the horizontal direction is located nearly 13 km east of the epicenter(36.123° N,82.499° E),and the best-fitting centroid depth is around 10 km.The total scalar moment,M0,is retrieved with an average value of 3.05 × 1019N·m,corresponding to moment magnitude MW6.92.Most of the energy is released within about 14 s.Moreover,we discuss about the potential application of this system in earthquake disaster decision.

Physical and chemical changes of water in the deep interior of the Earth―Decrepitation-style mud-volcano-earthquake―A bright lamp to shed light on the mysteries of the deep interior of the Earth

The 2008-05-12 Wenchuan mud-volcano-earthquake was accompanied with eruption of a huge volume of gas and stone,revealing that earthquakes generally result from instant reverse phase explosion of supercritical water(SCW) at the supercritical point.In the deep parts of the crust and mantle there still exists a large amount of supercritical water equivalent in order of magnitude to that of the Earth's hydrosphere.Soft fluids which exist in the MOHO at the top of the upper mantle are the so-called deep supercritical fluids(SCWD).Supercritical water(SCW) has n×103 times strong capability to dissolve gas.Its viscosity is extremely low and its diffusivity is extremely strong.Therefore,it can naturally migrate toward a region with relatively negative pressure.In the steep break zone of the MOHO at the 57-65 km depth beneath the earthquake belt,due to mutation of overburden pressure,SCWD can automatically separate out CaSiO3 and other inorganic salts,evolving into the SCW(H2O-CO2-CH4O system.In going upwards to the 10-20-km depth of the crust SCW will be accumulated as an earthquake-pregnant reservoir in the broken terrain.The phase-transition heat of SCW is estimated at 606.62 kJ/kg and the reverse phasing kinetic energy is 2350.8 kJ/kg.When automatic exhaust at the time of decompression reaches the critical pressure(Pc),the instant explosion reverse phase will be normal-state air water.It will release a huge volume of energy and high-kinetic-energy gas which has been expanded by a factor of 1000,leading to the breaking of the country rocks overlying the earthquake-pregnant reservoir,thus giving rise to a Ms 8.0 earthquake.As a result,there were formed eruptive and air-driven(pneumatic) debris flows whose volumatric flow rate reaches n×1014 m3/s,and their force greatly exceeds the power of INT explosive of the same equivalent value.

2023-08-23辽宁普兰店M 4.6地震震源参数测定及发震构造初判

2023-08-2318:19辽宁省大连市普兰店区发生M 4.6地震,为准确描述本次地震的震源特征,探讨其孕育和发震机理,本文通过测定地震的震源深度,反演震源机制、矩张量及质心深度,给出震源机制中心解;同时分析震源机制与构造应力场的关系,根据小震重定位结果对发震断层面进行拟合,初步确定了本次地震的发震断层。结果表明,普兰店M 4.6地震初始破裂深度为12.0 km,震源机制解参数为节面Ⅰ走向50°,倾伏角75°,滑动角-169°;节面Ⅱ走向317°,倾伏角80°,滑动角-15°,矩震级M W4.8,最优质心深度12 km;地震矩M_(0)为1.796×10^(16)Nm,矩张量解M_(rr)、M_(tt)、M_(pp)、M_(rt)、M_(rp)、M_(tp)分别为-0.004、0.946、-0.942、0.017、-0.305、-0.125;中心解参数为节面Ⅰ走向47.03°,倾伏角79.04°,滑动角-168.15°;节面Ⅱ走向314.75°,倾伏角78.37°,滑动角-11.19°。构造应力场作用在中心解节面Ⅰ上的相对剪应力为0.877,相对正应力为-0.544;投影于节面Ⅱ上的相对剪应力为0.911,相对正应力为0.161。拟合的断层面走向148.91°,倾伏角89.85°,其在构造应力体系下的滑动角为26.47°。综合分析认为,普兰店M 4.6地震发生在NW向普兰店-长海构造带,是沿应力场最优节面以左旋走滑为错动方式的天然地震。

短周期线性密集台阵揭示宣城地区浅部地壳速度结构及断裂发育特征

宣城市位于东南丘陵与长江中下游平原过渡地带,区内经历多期构造活动,地质构造复杂且断裂十分发育,开展宣城地区浅部速度结构和断裂带探测,不仅有助于深入了解区内地质构造和成矿作用,还可以为区域地震危险性和危害性评价提供重要的参考模型.本文在研究区布设了一条由110套三分量地震仪组成的短周期线性密集台阵,利用采集的一个月的三分量背景噪声数据,采用噪声谱比方法(HVSR)探测研究区场地峰值频率及基岩界面结构;通过垂直分量背景噪声数据重建经验格林函数,采用拓距相移法(ERPS)提取相速度频散曲线并反演线性台阵下方浅部地壳精细横波速度结构.基于HVSR计算的场地峰值频率和基岩界面埋深,刻画了测线下方浅部地表松散沉积层的结构,结果显示研究区内皖南山区基岩埋深较浅,盆地区域基岩埋深较深,地表松散沉积层的厚度最大可至地下80 m,此外还评估了区内地震破坏性和建筑的抗震性,认为研究区内平均振幅放大系数相对较高,局部区域可能会产生明显的地震放大效应,并且皖南山区低矮建筑(1~2层)及宣城—南陵盆地区域的高层建筑(7层及以上)在抗震设防上应当特别注意;依据反演的二维横波速度模型,获得了研究区浅部地壳6 km以深的结构,结果显示宣城—南陵盆地在浅部表现为明显的低速特征,其基底深度可达2 km,并且盆地部分区域呈现相对高速的异常特征,指示在区域多期构造活动中产生的岩浆侵入和逆冲推覆构造作用。此外,区内深部总体表现为高低速异常交替分布的“叠瓦状”构造模式,出现的多处低速异常带推断分别是周王断裂、江南断裂和茅山断裂以及清水河—河湾断裂的构造破碎带,指示了研究区内经历的多期次挤压—拉伸的强烈变形改造.综合上述研究成果,本研究为宣城地区地质构造条件的分析、地震危险性及危害性的评价以及区域找矿勘查等相关工作提供了新的依据.

Benchmark on the accuracy and efficiency of several neural network based phase pickers using datasets from China Seismic Network

Seismic phase pickers based on deep neural networks have been extensively used recently,demonstrating their advantages on both performance and efficiency.However,these pickers are trained with and applied to different data.A comprehensive benchmark based on a single dataset is therefore lacking.Here,using the recently released DiTing dataset,we analyzed performances of seven phase pickers with different network structures,the efficiencies are also evaluated using both CPU and GPU devices.Evaluations based on F1-scores reveal that the recurrent neural network(RNN)and EQTransformer exhibit the best performance,likely owing to their large receptive fields.Similar performances are observed among PhaseNet(UNet),UNet++,and the lightweight phase picking network(LPPN).However,the LPPN models are the most efficient.The RNN and EQTransformer have similar speeds,which are slower than those of the LPPN and PhaseNet.UNet++requires the most computational effort among the pickers.As all of the pickers perform well after being trained with a large-scale dataset,users may choose the one suitable for their applications.For beginners,we provide a tutorial on training and validating the pickers using the DiTing dataset.We also provide two sets of models trained using datasets with both 50 Hz and 100 Hz sampling rates for direct application by end-users.All of our models are open-source and publicly accessible.

Stochastic modeling for starting-time of phase evolution of random seismic ground motions

In response to the challenge inherent in classical high-dimensional models of random ground motions, a family of simulation methods for nonstationary seismic ground motions was developed previously through employing a wave-group propagation formulation with phase spectrum model built up on the frequency components’ starting-time of phase evolution. The present paper aims at extending the formulation to the simulation of non-stationary random seismic ground motions. The ground motion records associated with N–S component of Northridge Earthquake at the type-II site are investigated. The frequency components’ starting-time of phase evolution of is identified from the ground motion records, and is proved to admit the Gamma distribution through data fitting. Numerical results indicate that the simulated random ground motion features zeromean, non-stationary, and non-Gaussian behaviors, and the phase spectrum model with only a few starting-times of phase evolution could come up with a sound contribution to the simulation.

Olivine-spinel Phase Boundary Rheology, Strength and Deep Seismicity

In this paper,we discuss (1) whether a significant change in dominant creep mechanism will occur at 400 km discontinuity in some subducting slabs as a result of olivine-spinel phase transition;(2) could the result influence phase boundary strength and deep seismicity? Through this study,we noticed that a transition of creep mechanism from dislocation to diffusion (or superplasticity) could occur at the olivine-spinel phase boundary where temperature effect on phase transition dominates over pressure,which will result in a weakening phase boundary.Triggered by this phase transition weakening,a deep strong earthquake might be generated in the relatively strengthening region above the phase boundary so that the phase boundary is naturally the ceasing boundary of deep seismicity.Contrasting to this,the transition of the creep mechanism from dislocation glide to dislocation climb may be common at the phase boundary where pressure effect on phase transition dominates over temperature.In this case。

Amplitude Ratio of Regional Seismic Phases: An Efficient Indicator of Thermal State of the Lithosphere

By analyzing seismograms of short period records at the Beijing SeismoJogicaJ Observatory, the present paper investigates the amplitude ratio of seismic phases. The results indicate that the amplitude ratio of Sn/Lg is correlated with the lithosphere structure, the thermal state, and strong earthquake occurrence in the region the seismic rays pass through. The significance of such a correlation in the study on the genesis and prediction of strong earthquakes is discussed.

A study on small magnitude seismic phase identification using 1D deep residual neural network

Reliable seismic phase identification is often challenging especially in the circumstances of low-magnitude events or poor signal-to-noise ratio.With improved seismometers and better global coverage,a sharp increase in the volume of recorded seismic data has been achieved.This makes handling seismic data rather daunting by using traditional approaches and therefore fuels the need for more robust and reliable methods.In this study,we develop 1D deep Residual Neural Network(ResNet),for tackling the problem of seismic signal detection and phase identification.This method is trained and tested on the dataset recorded by the Southern California Seismic Network.Results demonstrate that the proposed method can achieve robust performance for the detection of seismic signals and the identification of seismic phases.Compared to previously proposed deep learning methods,the introduced framework achieves around 4%improvement in earthquake detection and a slightly better performance in seismic phase identification on the dataset recorded by Southern California Earthquake Data Center.The model generalizability is also tested further on the STanford EArthquake Dataset.In addition,the experimental result on the same subset of the STanford EArthquake Dataset,when masked by different noise levels,demonstrates the model’s robustness in identifying the seismic phases of small magnitude.

基于多点源的W-phase反演及其在2004年苏门答腊MW9.1大地震中的应用

在大地震发生后,快速准确地获得地震震源信息对应急救援十分重要,但是现有技术方法往往难以在快速准确地获取地震震源机制的同时获得破裂空间的分布特征。文章在传统W-phase反演技术的基础上开发了多点源W-phase反演方法,实现对大震破裂空间尺度上能量释放特征及震源机制的快速测定,并以2004年苏门答腊MW9.1大地震为例,测试程序的有效性。研究中设置了1、2、3、4、5、6个点源来分别测定此次地震的能量释放特征及震源机制。结果显示,震源机制随着空间位置由南向北的变化与俯冲面走向变化一致,与设定点源附近的历史地震震源机制高度吻合。因此,基于多点源的W-phase快速反投影技术将能更好获得大震空间能量释放特征,为震后应急及海啸预警提供科学支持。