Driving Galactic Turbulence by Supernova Explosions

We investigate the general properties of supernova driven interstellar turbulence using local three-dimensional MHD simulations under Galactic conditions. Our model includes the effects of large-scale shear due to Galactic differential rotation, density stratification, compressibility, magnetic fields, heating via supernova explosions and parameterized radiative cooling of the interstellar medium. In addition to investigating isolated supernova explosions we allow for multiple supernovae distributed randomly in the Galactic disc and exponentially in the vertical direction. Single supernova explosions drive a strong shock, the lifetime of which is approximately 2 Myr in our model. This stage is found to be characterized by a kinetic energy spectrum in the diffuse gas with spectral index consistent with k = –2. Large-scale shear and Coriolis force act on the supernova remnant producing some vorticity inside it, but this process was found to be very weak. In the case of multiple supernova explosions, older remnants have an important role causing density fluctuations in the interstellar medium. In this clumpy medium, the propagation velocity of the shock fronts changes due to the changing density, and vorticity is generated. In the absence of these supernova interactions the kinetic energy spectrum shows a relatively wide shock spectrum with spectral index k = –2, but when the supernova interactions become dominant the classical k = –5/3 spectrum is observed.     

Underground Nuclear Explosions Recorded At Rarotonga: Estimation of MB From T-Phase Amplitude

Summary. More than seventy underground nuclear explosions detonated in the Tuamotu Archipelago of French Polynesia have been detected on the short-period vertical Benioff seismograph at Rarotonga in the Cook Islands. Early data suggested that the duration of the T -phase provides a better estimate of the magnitude of the events than does the amplitude, but with more data now available it has been shown that amplitude is in fact the more reliable parameter. Data prior to 1980 are more scattered than those since. This appears to be related to the substantially different station network in those early years, compared with the later period, and suggests some unreliability in early data. Station corrections have been extracted and improved estimates of the magnitudes of the events obtained.     

Acoustic-Gravity Waves from Bolide Sources

We have developed a new approach to modeling the acoustic-gravity wave (AGW) radiation from bolide sources. This first effort involves entry modeling of bolide sources that have available satellite data through procedures developed in ReVelle (Earth Moon Planets 95, 441–476, 2004a; in: A. Milani, G. Valsecchi, D. Vokrouhlicky (eds) NEO Fireball Diversity: Energetics-based Entry Modeling and Analysis Techniques, Near-earth Objects: Our Celestial Neighbors (IAU S236), 2007b). Results from the entry modeling are directly coupled to AGW production through line source blast wave theory for the initial wave amplitude and period at (at 10 blast wave radii and perpendicular to the trajectory). The second effort involves the prediction of the formation and or dominance of the propagation of the atmospheric Lamb, edge-wave composite mode in a viscous fluid (Pierce, J. Acoust. Soc. Amer. 35, 1798–1807, 1963) as a function of the source energy, horizontal range and source altitude using the Lamb wave frequency that was deduced directly during the entry modeling and that is used as a surrogate for the source energy. We have also determined that Lamb wave production by bolides at close range decreases dramatically as either the source energy decreases or the source altitude increases. Finally using procedures in Gill (Atmospheric-Ocean Dynamics, 1982) and in Tolstoy (Wave Propagation, 1973), we have analyzed two simple dispersion relationships and have calculated the expected dispersion for the Lamb edge-wave mode and for the excited, propagating internal acoustic waves. Finally, we have used the above formalism to fully evaluate these techniques for four large bolides, namely: the Tunguska bolide of June 30, 1908; the Revelstoke bolide of March 31, 1965; the Crete bolide of June 6, 2002 and the Antarctic bolide of September 3, 2004. Due to page limitations, we will only present results in detail for the Revelstoke bolide.     

天津滨海爆炸事件检测与相对定位研究

2015年8月12日发生在天津滨海新区的化学品爆炸事故造成了严重的人员伤亡和财产损失。 基于区域地震台网记录的数字化连续地震波形资料, 对爆炸发生过程的精细分析是深入调查爆炸事故发生原因、 快速估计事故所造成灾害与损失的重要手段之一。 本文首先利用震相到时和质点运动轨迹信息对记录波形中能量较强的震相性质进行分析, 再通过模板匹配算法检测爆炸过程中是否存在两次大爆炸以外的小爆炸事件, 并对检测方法的效能进行检验, 最后依据波形互相关得到的震相到时差估计两次较大爆炸的相对位置。 计算结果显示区域台网记录波形中能量最大的体波震相可能为沿沉积层顶部传播的Pg震相, 垂直向能量最强的面波震相为Rg震相, 而能量较弱的初至震相可能为结晶基底首波; 爆炸过程中除了两次能量较大的爆炸以外还能检测到两次较小的爆炸事件, 其当量相当于约M_L0.5的地震; 基于两次大爆炸Pg震相相对到时信息得到第二次大爆炸发生在第一次大爆炸的北西侧, 两次爆炸位置距离约50~55 m。 本文研究结果为爆炸事故的调查分析提供了地震学依据。     

平顶山地区爆破识别判据分析

对2008年1月1日-2012年12月31日河南数字地震台网记录的平顶山地区的地震与爆破事件,从初动方向、振幅比、周期比和记录波形对比等进行分析,结果表明：振幅比和记录波形对比可以作为平顶山地区爆破识别判据,初动方向可以作为该地区爆破识别的辅助判据.     

地下储液罐抗爆炸地冲击作用的流固耦合有限元分析

储液罐的抗震性能虽然得到了较为深入的研究,但对其抗爆炸产生的地冲击作用的研究还十分罕见。本文通过建立流固耦合的数值分析模型,研究爆炸地冲击作用下地下立式储液罐的动力特性,包括储液晃动波高、罐壁的应力和应变、底板的提离和浮放储液罐的“象足”变形产生原因及破坏机理等。研究结果表明,储液罐在爆炸地冲击荷载和地震荷载作用下的动力响应有明显区别,地震荷载作用下的储液罐抗震验算方法不适用于爆炸地冲击荷载。     

The Influence of Path Corrections and a Three-dimensional Global <Emphasis Type="Italic">P</Emphasis>-wave Velocity Model on Seismic Event Location in Kazakhstan

We relocate 81 large nuclear explosions that were detonated at the Balapan and Degelen Mountain subregions of the Semipalatinsk test site in Kazakhstan during the years 1978 to 1989. The absolute locations of these explosions are available, as well as very accurate estimates of their origin times. This ground truth information allows us to perform a detailed analysis of location capability. We use a sparse network of stations with highly accurate first arrival picks measured using a waveform cross-correlation method. These high quality data facilitate very accurate location estimates with only a few phases per event. We contrast two different approaches: 1) a calibration-based approach, where we achieve improved locations by using path corrections, and 2) a model-based approach, where we achieve improved locations by relocating in a recently published global 3-D P-velocity model. Both methods result in large improvements in accuracy of the obtained absolute locations, compared to locations obtained in a 1-D reference earth model (ak135). The calibration-based approach gives superior results for this test site, in particular when arrival times from regional stations are included. Estimated locations remain well within a 1000 km² region surrounding the ground truth locations when the path corrections for the Balapan and Degelen Mountain subregions are interchanged, but even for the short separation between these two regions, we find variations in the path corrections that cause systematic mislocations. The model-based approach also results in substantially reduced mislocation distances and has the distinct advantage that it is, in principle, transportable to other source regions around the world.