天然源面波勘探台阵对比试验

为了对比天然源面波勘探不同台阵布局的探测效果, 筛选出探测成果可靠、 效率高和便于野外施工的天然源面波勘探台阵阵形, 在天水市黄土覆盖区的同一场地分别用4种常见的阵形进行数据采集试验, 并对各种阵形数据使用空间自相关法或扩展空间自相关法提取相应的频散曲线, 通过反演得到了试验点地下的浅层速度结构模型. 分析对比试验结果表明: 4种台阵提取的频散曲线数值很相近; 频散谱能量集中度较高的是嵌套式等边三角形和圆形台阵, L形和直线形台阵相对分散; L形台阵低频段(4—8 Hz)比直线形台阵差, 其高频段(8—40 Hz)比直线形台阵好. 针对直线形台阵在高频段信噪比较低的情况, 在确保探测成果可靠性的前提下, 为了提高探测效率, 提出了在同一直线形台阵开展天然源与人工源面波联合勘探的数据采集方法. 实验结果证实, 这种联合勘探方法不仅可弥补直线形台阵高频段的不足, 确保探测精度和结果的可靠性, 而且还能实现“高低”频兼顾, 即“深浅”兼顾.      

Identification of blasting sources in the Dobrogea seismogenic region,Romania using seismo-acoustic signals

In order to discriminate between quarry blasts and earthquakes observed in the Dobrogea seismogenic region, a seismo-acoustic analysis was performed on 520 events listed in the updated Romanian seismic catalogue from January 2011 to December 2012. During this time interval, 104 seismo-acoustic events observed from a distance between 110 and 230 km and backazimuth interval of 110–160° from the IPLOR infrasound array were identified as explosions by associating with infrasonic signals. WinPMCC software for interactive analysis was applied to detect and characterize infrasonic signals in terms of backazimuth, speed and frequency content. The measured and expected values of both backazimuths and arrival times for the study events were compared in order to identify the sources of infrasound. Two predominant directions for seismo-acoustic sources’ aligning were observed, corresponding to the northern and central parts of Dobrogea, and these directions are further considered as references in the process of discriminating explosions from earthquakes. A predominance of high-frequency detections (above 1 Hz) is also observed in the infrasound data. The strong influence of seasonally dependent stratospheric winds on the IPLOR detection capability limits the efficiency of the discrimination procedure, as proposed by this study.     

Seismo-acoustic signals associated with degassing explosions recorded at Shishaldin Volcano,Alaska, 2003–2004

In summer 2003, a Chaparral Model 2 microphone was deployed at Shishaldin Volcano, Aleutian Islands, Alaska. The pressure sensor was co-located with a short-period seismometer on the volcano’s north flank at a distance of 6.62 km from the active summit vent. The seismo-acoustic data exhibit a correlation between impulsive acoustic signals (1–2 Pa) and long-period (LP, 1–2 Hz) earthquakes. Since it last erupted in 1999, Shishaldin has been characterized by sustained seismicity consisting of many hundreds to two thousand LP events per day. The activity is accompanied by up to ∼200 m high discrete gas puffs exiting the small summit vent, but no significant eruptive activity has been confirmed. The acoustic waveforms possess similarity throughout the data set (July 2003–November 2004) indicating a repetitive source mechanism. The simplicity of the acoustic waveforms, the impulsive onsets with relatively short (∼10–20 s) gradually decaying codas and the waveform similarities suggest that the acoustic pulses are generated at the fluid–air interface within an open-vent system. SO₂ measurements have revealed a low SO₂ flux, suggesting a hydrothermal system with magmatic gases leaking through. This hypothesis is supported by the steady-state nature of Shishaldin’s volcanic system since 1999. Time delays between the seismic LP and infrasound onsets were acquired from a representative day of seismo-acoustic data. A simple model was used to estimate source depths. The short seismo-acoustic delay times have revealed that the seismic and acoustic sources are co-located at a depth of 240±200 m below the crater rim. This shallow depth is confirmed by resonance of the upper portion of the open conduit, which produces standing waves with f=0.3 Hz in the acoustic waveform codas. The infrasound data has allowed us to relate Shishaldin’s LP earthquakes to degassing explosions, created by gas volume ruptures from a fluid–air interface.     

The Finnish seismic array

This paper deals with the Finnish seismic array, which in its present state consists of two tripartite arrays, HESA and JYSA. The first of these arrays, situated near Helsinki, has been in continuous operation since 1964 and uses analog recording by frequency-modulated telemetry. The JYSA array is sited in central Finland near the town of Jyväskylä. This array consists of three substations equipped with vertical-component SP seismometers. The signals are transmitted by digital FM telemetry in the UHF band and recorded at Helsinki on magnetic tape. Events are detected visually from monitoring drums. The digital data of the events detected are handled by a microprocessor with a graphic display and fed to a Burroughs 6700 analysing computer.Comparison with the NORSAR array indicates that at frequencies below 2 Hz the absolute noise level at JYSA is about 6 dB lower, but above 2 Hz the noise level for single instruments at the two arrays is about the same. As regards detection capability, the incremental threshold is about the same for JYSA as for Hagfors and 0.5 magnitude units higher than for NORSAR, but 0.5 magnitude units lower for JYSA than for the WWSSN station NUR in southern Finland.     

Locating the Tohoku-Oki 2011 tsunami source using acoustic–gravity waves

The giant Tohoku-Oki earthquake of 11 March 2011 in offshore Japan did not only generate tsunami waves in the ocean but also infrasound (or acoustic–gravity) waves in the atmosphere. We identified ultra-long-period signals (>500 s) in the recordings of infrasound stations in northeast Asia, the northwest Pacific, and Alaska. Their source was found close to the earthquake epicenter. Therefore, we conclude that in general, infrasound observations after a large offshore earthquake are evidence that the surface and the floor of the sea have been significantly vertically displaced by the earthquake and that a tsunami must be expected. Since infrasound is traveling faster than the tsunami, such information may be used for tsunami early warnings.     

Characterizing complex eruptive activity at Santiaguito,Guatemala using infrasound semblance in networked arrays

We implement an infrasound semblance technique to identify acoustic sources originating from volcanic vents and apply the technique to the generally low-amplitude infrasound (< 3 Pa at 1 km) signals produced by Santiaguito dome in Guatemala. Semblance detection is demonstrated with data collected from two-element miniature arrays with ~ 30 m spacing between elements. The semblance technique is effective at identifying a range of eruptive phenomena, including pyroclastic-laden eruptions, vigorous degassing events, and rockfalls, even during periods of high wind contamination Many of the detected events are low in amplitude (tens of mPa) such that they are observed only by select arrays positioned with proximity and line-of-sight to the source. Larger events, such as the pyroclastic-laden eruptions, which occurred bi-hourly in 2009, were detected by all five arrays and produced an infrasonic signal that was correlated across the network. Network correlated events can be roughly located and map to the summit of the Caliente Vent where pyroclastic-laden eruptions originate. In general, the degree of Santiaguito infrasound event correlation is poor across the network, suggesting that complex source geometry contributes to asymmetric sound radiation.     

Detection of Impact Points of Fragments of Spent Launch Vehicle Stages Using Infrasound Direction-Finding Methods

The paper describes the principles and techniques used to detect signals propagating in the atmosphere in the infrasonic frequency range. Such signals can be generated by different sources: ground and atmospheric explosions, as well as objects moving in the atmosphere at supersonic speed (aircraft, rockets, bolides, fragments of spent stages of launch vehicles). Portable infrasound monitoring stations are described, each of which includes three spaced infrasonic microphones. Each such station makes it possible to determine three basic parameters of the detected infrasound signal, which are subsequently used to solve the direction- finding problem: the time of arrival of an infrasonic wave, the azimuth to the source in the horizontal plane, and the wave approach angle from the source of infrasonic waves to the Earth’s surface in the vertical plane. An acoustic detector used to extract useful signals against a noise background is described. The detector is based on an algorithm similar to the STA/LTA detection algorithm known in seismology. Examples of the operation of an acoustic detector with data obtained during real measurements are given. Passive infrasound direction-finding technology is described. It is based on mathematical modeling of the of infrasonic wave propagation in the atmosphere, which are generated by objects moving along possible trajectories; comparison of theoretical signals with real ones recorded by monitoring stations; and determination of the realized trajectories. The paper gives examples of experimental verification of the effectiveness of passive infrasound direction-finding technology for determining the impact points of the first and second stages of launch vehicles. It is shown that infrasound direction-finding systems makes it possible to reduce the estimated search area for launch vehicle fragments that fall to the Earth, significantly decrease the time and costs for their search and utilization, and mitigate the negative environmental impact of the rocket and space industry.     

Some results of microseismic monitoring in the Elbrus Region

The results of microseismic monitoring in the adit of the neutrino laboratory (Institute of Nuclear Research, Russian Academy of Sciences) in the Elbrus Region are presented. Three-component computer recording was carried out in the frequency band 10–120 Hz with 5-second averaging of the amplitudes of vibrations. We consider the results of processing the data yielded by the monitoring in the adit during nine months since March 2008. The seismic field at the observation point is composed of three main contributions: the local earthquakes associated with two main fluid-magmatic sources, which are located nearby the adit; a slowly varying wave field, which is probably related to the same sources; and the industrial noise, namely, the vibrations caused by the traffic of carriages and a weak stationary background noise.     

Source location variability and volcanic vent mapping with a small-aperture infrasound array at Stromboli Volcano,Italy

Stromboli Volcano in Italy is a persistently active, complex volcanic system. In May 2002 activity was confined to 3 major summit craters within which several active vents hosted multiple explosions each hour. During a 5-day field campaign an array of 3 low-frequency microphones was installed to investigate the coherent infrasound produced by degassing from these vents. Consistent phase lags across the 3 stations indicate distinct sources that are subsequently investigated to determine the associated vent location, apparent depth, and origin time. The cross-correlation routine allows for variations in comparison window length, waveform filtering bandwidth, and correlation and consistency thresholds, allowing for improved detection of certain types of degassing sources. Identification of activity at the various vents could be subsequently corroborated with 3 channels of synchronously acquired thermal data and video. During the May 2002 experiment persistent, energetic infrasound was observed from a passive degassing source within the Central Crater (CC) and transient infrasound, produced by discrete Strombolian explosions, was identified at 4 additional vents. The continuous infrasound produced by the CC exhibits variable frequency-dependent correlation lag times that are interpreted as a diffraction effect due to the acoustic radiators recessed location within a steep-walled crater. Such dispersion has important implications for accurate eruption source modeling because it indicates that infrasonic waveforms may be significantly filtered during propagation. Transient explosion signals from the Northeast Crater (NEC) and Southwest Crater (SWC) vents also exhibit dynamic correlation lag times, but this scatter may be more reasonably attributed to variable epicentral locations. Explosions from the NEC west vent, for instance, appear to emanate from a diffuse zone with a lateral extent in excess of 10 m.Editorial responsibility: R. Cioni     

Analysis of Signals from an Unique Ground-Truth Infrasound Source Observed at IMS Station IS26 in Southern Germany

Quantitative modeling of infrasound signals and development and verification of the corresponding atmospheric propagation models requires the use of well-calibrated sources. Numerous sources have been detected by the currently installed network of about 40 of the final 60 IMS infrasound stations. Besides non-nuclear explosions such as mining and quarry blasts and atmospheric phenomena like auroras, these sources include meteorites, volcanic eruptions and supersonic aircraft including re-entering spacecraft and rocket launches. All these sources of infrasound have one feature in common, in that their source parameters are not precisely known and the quantitative interpretation of the corresponding signals is therefore somewhat ambiguous. A source considered well-calibrated has been identified producing repeated infrasound signals at the IMS infrasound station IS26 in the Bavarian forest. The source results from propulsion tests of the ARIANE-5 rocket’s main engine at a testing facility near Heilbronn, southern Germany. The test facility is at a range of 320 km and a backazimuth of ~280° from IS26. Ground-truth information was obtained for nearly 100 tests conducted in a 5-year period. Review of the available data for IS26 revealed that at least 28 of these tests show signals above the background noise level. These signals are verified based on the consistency of various signal parameters, e.g., arrival times, durations, and estimates of propagation characteristics (backazimuth, apparent velocity). Signal levels observed are a factor of 2–8 above the noise and reach values of up to 250 mPa for peak amplitudes, and a factor of 2–3 less for RMS measurements. Furthermore, only tests conducted during the months from October to April produce observable signals, indicating a significant change in infrasound propagation conditions between summer and winter months.     

Ambient noise levels and detection threshold in Norway

Ambient seismic noise is caused by a number of sources in specific frequency bands. The quantification of ambient noise makes it possible to evaluate station and network performance. We evaluate noise levels in Norway from the 2013 data set of the Norwegian National Seismic Network as well as two temporary deployments. Apart from the station performance, we studied the geographical and temporal variations, and developed a local noise model for Norway. The microseism peaks related to the ocean are significant in Norway. We, therefore, investigated the relationship between oceanic weather conditions and noise levels. We find a correlation of low-frequency noise (0.125–0.25 Hz) with wave heights up to 900 km offshore. High (2–10 Hz) and intermediate (0.5–5 Hz) frequency noise correlates only up to 450 km offshore with wave heights. From a geographic perspective, stations in southern Norway show lower noise levels for low frequencies due to a larger distance to the dominant noise sources in the North Atlantic. Finally, we studied the influence of high-frequency noise levels on earthquake detectability and found that a noise level increase of 10 dB decreases the detectability by 0.5 magnitude units. This method provides a practical way to consider noise variations in detection maps.     

首都圈数字地震台网对微弱爆破信号的检测能力

利用首都圈数字地震台网接收人工地震信号,进行地下结构研究具有重要意义.但人工震源释放的能量小,激发的地震波以短周期为主,因此本文较全面地研究了地震台网对短周期微弱信号（1～20 Hz）的检测能力：(1) 分析了台网的背景噪声,结果表明基岩台址的地震台噪声比沉积盖层台址的地震台噪声低约13 dB,这相当于近1个震级的检测阈值;夜间的噪声比白天低约5 dB;噪声有逐年增高的趋势,2006年比2001年噪声提高约4 dB.(2 )分析了在台网内进行的药量为25 kg的陆地井下爆破实验,一次爆破相当于0.69级（M_L）的天然地震,有18个地震台可辨认爆炸产生的Pg、Pm或Pc波;离爆破点218 km的基岩台,仍可以接收到振幅只有1.6 nm 的Pm波,这个结果可为地震勘探实际工作提供参考.(3) 研究了台网外核爆试验的信号特征,2006年发生在朝鲜的地下核试验是一次检验台网检测微弱信号能力的好机会.波形记录经1～5Hz滤波后,台网中噪声小的18个基岩台可以清晰辨认核爆破产生的P波或Lg波,P波平均振幅为16 nm,计算的平均震级为m_b4.3,和NEIC给出的震级相同;分析还表明背景噪声是影响台站信号检测能力的主要因素之一.     

Infrasound Strain Perturbations Caused by Typhoons

Izvestiya, Physics of the Solid Earth - The paper addresses the excitation of high-power infrasound disturbances in the 7–9 Hz frequency band associated with the passage of tropical cyclones...     

Fine structure of Pi2-type geomagnetic pulsations

An analysis of the BEAR experiment and Polar satellite data showed that several ionospheric sources could act almost simultaneously in the area of Pi2-type geomagnetic pulsations. The aim of this study is to locate these sources and determine their coherence in the narrow-frequency band of (6–10) × 10⁻³ Hz. By using gradient analysis methods, we revealed the local sources of oscillations coherent (phased), that allows one to treat their contributions as manifestations of the fine structure of the field of Pi2 pulsations.     

Seismic activity related to the 2000 eruption of the Hekla volcano, Iceland

The 2000 Hekla eruption took place from February 26 to March 8. Its seismic expressions were a swarm of numerous small earthquakes related to its onset, and low-frequency volcanic tremor that continued throughout the eruption. A swarm of small earthquakes was observed some 80 min before the onset of the eruption, and the size of the events increased with time. Low-frequency volcanic tremor, with a characteristic frequency band of 0.5–1.5 Hz and dominant spectral peak(s) at 0.7–0.9 Hz, became visible at 18:19 GMT on February 26, marking the onset of the eruption. The tremor amplitude rose quickly and was very high in the beginning of the eruption. However, it soon began to decrease after about an hour. In general, the seismic activity related to the 2000 Hekla eruption was very similar to what was observed in the previous eruption in 1991. Based on knowledge gained from seismicity and strain observations from 1991, this was the first time that a Hekla eruption was predicted.Editorial responsibility: J Stix     

Swelling of a lava plug associated with a Vulcanian eruption at Sakurajima Volcano,Japan, as revealed by infrasound record: case study of the eruption on January 2, 2007

In order to clarify the time relation of the expansion of a gas pocket and failure of its overlying plug of lava during Vulcanian eruptions, infrasound records and video images of the Vulcanian eruption that occurred at Sakurajima volcano on January 2, 2007 were analyzed with respect to their origin times. Weak (≤3 Pa) and slowly increasing air pressure preceded the impulsive compression phase by 0.25–0.32 s, and a longer-period rarefaction phase of infrasound waves was recognized at all microphone stations. The velocity of the compression phase was assumed to be supersonic (ca. 400 m/s) up to 850 m above the crater bottom from other recent explosions. On the other hand, the propagation velocity of the preceding weak signal was regarded to be similar to the air sound velocity because the lack of impulsiveness is unlikely to be related to the main compression phase. Therefore, the estimated origin time of the main compression phase was delayed by 0.5–0.7 s from the preceding phase. The origin time of the preceding phase coincided with the onset of the isotropic expansion process of the pressurized gas pocket, which was obtained by the waveform inversion of the explosion earthquake. In contrast, the origin time of the main impulsive phase coincided with the time when the expansion rate reached its peak. This observation suggests that the volumetric increase of the gas pocket caused swelling of the surface of the crater bottom and its subsequent failure. When the expansion velocity exceeded a threshold level, the main impulsive compression phase radiated with a high velocity by the sudden releases of the pressurized gases. The volumetric change at the source was estimated to be 280–560 m³ from the preceding phase of the infrasound. This volume change indicates that the vertical displacement of the swelling ground was on the order of 1.0 m, assuming the radius of the lava plug was ca. 10 m.     

Comment on “Small-scale heterogeneities below the Lanzhou CTBTO seismic array,from seismic wave field fluctuations” by Shen and Ritter

A so-called teleseismic fluctuation wavefield method (TFWM) was adopted by Shen and Ritter to study small-scale heterogeneities beneath the Lanzhou array, a small-aperture seismic array locating in the northeastern margin of the Tibetan plateau. Through this method, they speculated a correlation length of 2.4–8.2 km and a root-mean-squared (RMS) velocity perturbation of 1–3% for the crust under the array. However, for such a tectonically active and heavily deformed region where great earthquakes up to M7 occurred, such a low RMS velocity perturbation proposed seems to be unconvincing and conflicts with other studies. To confirm the suspicion, the acoustic finite difference method is used to model scattering wavefields in the random medium created through parameters suggested by Shen and Ritter. TFWM is applied to the synthetic data from both the small- and large-aperture arrays. For synthetic data from arrays with similar aperture as the Lanzhou array, it fails to reproduce the parabolic fit for real data, implying that the statistical parameters proposed by Shen and Ritter are incorrect. Another test presuming a correlation length of 2.4 km and a RMS velocity perturbation of 5% can reproduce the fit for real data. The numerical experiments prove that it is improper to apply TFWM to small-aperture arrays like the Lanzhou array, or else the heterogeneity will be heavily underestimated. The aperture of an array must be much larger than the correlation length and wavelength to get unbiased measurements.     

1D velocity structure of the Po River plain (Northern Italy) assessed by combining strong motion and ambient noise data

Strong ground motions recorded on the sedimentary deposits of the Po River alluvial plain during the Emilia (Northern Italy) Mw 5.7 earthquake of May 29, 2012 are used to assess the vertical profile of shear-wave velocity above the limestone basement. Data were collected by a linear array installed for site effect studies after the Mw 5.9 mainshock of May 20, 2012. The array stations, equipped with both strong and weak motion sensors, are aligned in the South–North direction, at distances ranging from 1 to 26 km from the epicenter. The vertical components of ground motion show very distinctive, large-amplitude, low-frequency dispersive wave trains. Wavelet analysis yields group-velocity dispersion curve in the 0.2–0.7 Hz frequency band. The availability of a long ambient noise record allows estimates of the site resonance frequency along with its stability among stations. The joint inversion of dispersion of surface waves and ellipticity curves derived from ambient noise H/V allows extending investigations down to the sediment-limestone interface, at a depth of about 5,000 m. Our results add new information about the velocity structure at a scale that is intermediate between the local scale already investigated by other authors with small-aperture arrays using ambient noise and the regional scale inferred from modeling of seismogram waveforms recorded at hundreds of kilometers from the source.     

Application of seismo-acoustic signals to the study of local site effects

The Nevada Seismic Array (NVAR) is a small-aperture seismic array designed for monitoring an eventual nuclear test ban treaty. In spite of the 4 km aperture, large amplitude variations are recorded due to the complicated local geology. This study takes advantage of the collocated infrasound and seismic sensors to discuss the use of air-to-ground coupled waves to characterize the shallow geological structure existing beneath the array.     

2003年日本北海道8.0级地震次声波特征研究

利用在北京昌平安装的次声三点阵, 记录了2003年9月26日日本北海道地震的前兆次声波和震后次声波。 实际记录的P-t曲线及波速波向图, 经快速傅立叶变换得到其三维动态频谱。 结果表明: ① 强震前2～7天能测到振幅很强、 方向可测的地震前兆次声波。 其频谱特点是: 振幅由弱(10～15 Pa)逐渐加强到(50～80 Pa)或更大; 先为短周期波(10～30 min), 紧跟着长周期波(30～50 min), 然后又出现更长的周期波, 最后长短周期的波一起出现连成一片(2～62 min); ② 地震前兆次声波的产生和记录较容易, 而当地次声波和震中次声波记录和识别比较困难; ③ 若两地的三点阵波向相交, 可预测震中位置。 故地震前兆次声波的测量研究, 有可能发展成为临震预报中一种有效的新方法