Statistical analysis of intermittent volcanic tremor associated with the September 1989 summit explosive eruptions at Mount Etna, Sicily

The pattern of volcanic tremor accompanying the 1989 September eruption at the south-east summit crater of Mount Etna is studied. In specific, sixteen episodes of lava fountaining, which occurred in the first phase of the eruption, are analysed. Their periodic behaviour, also evidenced by autocorrelation, allows us to define the related tremor amplitude increases as intermittent volcanic tremor episodes. Focusing on the regular intermittent behaviour found for both lava fountains and intermittent volcanic tremors, we tried an a posteriori forecast using simple statistical methods based on linear regression and the Student’ t-test. We performed the retrospective statistical forecast, and found that several eruptions would have been successfully forecast. In order to focus on the source mechanism of tremor linked to lava fountains, we investigated the relationship between volcanic and seismic parameters. A mechanism based on a shallow magma batch ‘regularly’ refilled from depth is suggested.     

Volcanic tremor at Mt. Etna: Inferences on magma dynamics during effusive and explosive activity

During 1999, the volcanic activity at Mt. Etna was both explosive and effusive at the summit craters: Strombolian activity, lava fountains and lava flows affected different areas of the volcano, involving three of the four summit craters. Results from analysis of the 1999 volcanic tremor features are shown at two different time scales. First, the long-term time variation of the features of the volcanic tremor (including spectral and polarization parameters), during the entire year, was compared with the evolution of the eruptive activity. This approach demonstrated the good agreement between tremor data and observed eruptive activity; the activation of different tremor sources was suggested. Then, a more refined analysis of the volcanic tremor, recorded during 14 lava fountain eruptions, was performed. In particular, a shift of the dominant frequencies towards lower values was noted which corresponds with increasing explosive activity. Similar behaviour in the frequency content has already been observed in other explosive eruptions at Mt. Etna as well as on other volcanoes. This behaviour has been explained in terms of either an increase in the tremor source dimension or a decrease in the sound speed in the magma within the conduit. These results confirm that the volcanic tremor is a powerful tool for better understanding the physical processes controlling explosive eruptions at Mt. Etna volcano.     

Seismic recordings of subterranean volcanic activity at Ruapehu during 1964

Volcanic vibrations from Ruapehu volcano which is situated in the centre of the North Island of New Zealand have been recorded on high magnification slow motion tape seismographs and drum seismographs at the Chateau Volcanological Observatory since December, 1960. In 1964, volcanic microtremor with a dominant frequency of 2 c/s commenced late in March, and reached a peak seismic power level of 100 KW for short periods in May. During the maximum phase, the tremor completely ceased for up to 20 minutes, and recommenced with an explosion or burst of strong tremor. The sequence and timing were very similar to the sequence of ash discharge-stoppage-explosion-ash discharge observed during the 1945 eruption of Ruapehu. In 1964 it is thought that the eruption of ash clouds was prevented by the Crater Lake, so that visible activity was limited to a few fumaroles, steam rising from the lake, and turbulence in the water. The lake temperature increased from about 25°C (a normal temperature) on 20 March to 50°C on 26 May. The corresponding rates of heat loss from the lake are 200 and 700 MW respectively, and an additional 150 MW or more was required to heat up the lake, giving a total average heat output of 850 MW for about 4 weeks. The corresponding average seismic power was about 3 KW, which is 0.0005 per cent of the thermal power. If this relationship is constant, the peak thermal power over a period of 6 minutes was of the order of 20,000 MW. The explosions initiating the tremor were mostly identical except in amplitude, and their magnitudes ranged up to 2.3, corresponding to 10¹⁴ erg. This was slightly greater than the energy conserved during the preceding stoppage of tremor. The dominant tremor frequency was normally 2.2 c/s and sometimes 1.2 c/s. Individual bands of frequency within the spectrum varied in power, but coherent migrations of frequency bands occurred on a few occasions. Near the end of the maximum phase, the entire spectrum migrated downwards by half an octave and the tremor power decreased to zero. Tremor power and frequency increased again, and a violent period of explosions and tremor of changing frequency occurred. Quite frequently the explosions initiating the tremor were multiple, and a few explosions occurred which were followed by tremor lasting only a matter of seconds. Explosions apparently unrelated to tremor were very rare and minor. The explosions were located by temporary installation of portable slow motion tape recorders around the volcano. The epicentre is very close to Crater Lake, but the depth cannot be determined from the data. Graphs of tremor power against time covering the active period from April to September, and graphs of cumulative energy against time, and frequency spectrum against time, for explosions and tremor are presented, and slow motion tape recordings will be played many times faster than the recording speed so that the tremor and explosion vibrations can be heard as audible noises.     

New insights into the 1999 eruption of Shishaldin volcano,Alaska, based on acoustic data

Data collected by a pressure sensor provide new insights into the 1999 eruption of Shishaldin volcano, Unimak Island, Alaska. On 19 April 1999, after 3 months of unrest and an extended period of low-level Strombolian activity, Shishaldin experienced a Subplinian eruption (ash plume to >16 km), followed by several episodes of strong Strombolian explosions. Acoustic data from the pressure sensor allow us to investigate the details of an eruption which was instrumentally well recorded, but with few visual observations. In the 12 h prior to the Subplinian phase, the pressure sensor detected a series of small, repeated pulses with a constant spectral peak at 2–3 Hz. The amplitude and occurrence rate of the pulses both grew such that the signal became a nearly continuous hum just before the Subplinian eruption. This humming signal may represent gas release from rising magma. The main Subplinian phase was heralded by (1) the abrupt end of the humming signal, (2) several pulses of low-frequency sound interpreted as ash bursts, and (3) a dramatic increase in seismic tremor amplitude. The change in acoustic signature at this time allows us to precisely time the start of the Subplinian eruption, previously approximated as the time of strongest tremor increase. The 50-min Subplinian phase actually contained several bursts of signal, each of which may represent a discrete volume of magma passing through the system. Following the Subplinian event, the pressure sensor recorded four discrete episodes of Strombolian gas explosions on 19–20 April and another on 22–23 April. Four of the five episodes were accompanied by strong seismic tremor; the fifth has not been previously recognized and was not associated with anomalous tremor amplitudes. In time series these events are similar to explosions recorded at other volcanoes but in general they are much larger, with maximum amplitudes of >65 Pa at 6.5 km from the vent, and they have low (0.7–1.5 Hz) peak frequencies. These large explosions occurred at rates of 3–20 per minute for 1–5 h in each episode. The explosions were accompanied by a small (<5 km above sea level) ash plume and only minor amounts of ejecta were produced. Thus, the explosion activity was dominated by gas release.     

Volcanic tremor related to the 1991 eruption of the Hekla volcano,Iceland

Volcanic tremor at the Hekla volcano is directly related to eruptive activity. It starts simultaneously with the eruptions and dies down at the end of them. No tremor at Hekla has been observed during non-eruptive times. The 1991 Hekla eruption began on 17 January, after a short warning time. Local seismograph stations recorded small premonitory earthquakes from 16:30 GMT on. At 17:02 GMT, low-frequency volcanic tremor became visible on the seismograph records, marking the onset of the eruption. The initial plinian phase of the eruption was short-lived. During the first day several fissures were active but, by the second day, the activity was already limited to a segment of one principal fissure. The eruption lasted almost 53 days. At the end of it, during the early hours of 11 March, volcanic tremor disappeared under the detection threshold and was followed by a swarm of small earthquakes. At the start of the eruption, the tremor amplitude rose rapidly and reached a maximum in only 10 min. The tremor was most vigorous during the first hour and started to decline sharply during the next hour, and later on more gently. During the eruption as a whole, the tremor had a continuous declining trend, with occasional increases lasting up to about 2 days. Spectral analysis of the tremor during the first 7 h of the eruption shows that it settled quickly, within a couple of minutes, to its characteristic frequency band, 0.5–1.5 Hz. The spectrum had typically one dominant peak at 0.7–0.9 Hz, and a few subdominant peaks. Hekla tremor likely has a shallow source. Particle motion plots suggest that it contains a significant component of surface waves. The tremor started first when the connection of the magma conduit with the atmosphere was reached, suggesting that degassing may contribute to its generation.     

Observations of volcanic tremor during the January–February 2005 eruption of Mt. Veniaminof,Alaska

Mt. Veniaminof, Alaska Peninsula, is a stratovolcano with a summit ice-filled caldera containing a small intracaldera cone and active vent. From January 2 to February 21, 2005, Mt. Veniaminof erupted. The eruption was characterized by numerous small ash emissions (VEI 0 to 1) and accompanied by low-frequency earthquake activity and volcanic tremor. We have performed spectral analyses of the seismic signals in order to characterize them and to constrain their source. Continuous tremor has durations of minutes to hours with dominant energy in the band 0.5–4.0 Hz, and spectra characterized by narrow peaks either irregularly (non-harmonic tremor) or regularly spaced (harmonic tremor). The spectra of non-harmonic tremor resemble those of low-frequency events recorded simultaneously with surface ash explosions, suggesting that the source mechanisms might be similar or related. We propose that non-harmonic tremor at Mt. Veniaminof results from the coalescence of gas bubbles while low-frequency events are related to the disruption of large gas pockets within the conduit. Harmonic tremor, characterized by regular and quasi-sinusoidal waveforms, has duration of hours. Spectra containing up to five harmonics suggest the presence of a resonating source volume that vibrates in a longitudinal acoustic mode. An interesting feature of harmonic tremor is that frequency is observed to change over time; spectral lines move towards higher or lower values while the harmonic nature of the spectra is maintained. Factors controlling the variable characteristics of harmonic tremor include changes in acoustic velocity at the source and variations of the effective size of the resonator.     

Volcanic Tremor at Mt. Etna,Italy, Preceding and Accompanying the Eruption of July – August, 2001

The July 17 – August 9, 2001 flank eruption of Mt. Etna was preceded and accompanied by remarkable changes in volcanic tremor. Based on the records of stations belonging to the permanent seismic network deployed on the volcano, we analyze amplitude and frequency content of the seismic signal. We find considerable changes in the volcanic tremor which mark the transition to different styles of eruptive activity, e.g., lava fountains, phreatomagmatic activity, Strombolian explosions. In particular, the frequency content of the signal decreases from 5 Hz to 3 Hz at our reference station ETF during episodes of lava fountains, and further decreases at about 2 Hz throughout phases of intense lava emission. The frequency content and the ratios of the signal amplitude allow us to distinguish three seismic sources, i.e., the peripheral dike which fed the eruption, the reservoir which fed the lava fountains, and the central conduit. Based on the analysis of the amplitude decay of the signal, we highlight the migration of the dike from a depth of ca. 5 km to about 1 km between July 10 and 12. After the onset of the effusive phase, the distribution of the amplitude decay at our stations can be interpreted as the overall result of sources located within the first half kilometer from the surface. Although on a qualitative basis, our findings shed some light on the complex feeding system of Mt. Etna, and integrate other volcanological and geophysical studies which tackle the problem of magma replenishment for the July–August, 2001 flank eruption. We conclude that volcanic tremor is fundamental in monitoring Mt. Etna, not only as a marker of the different sources which act within the volcano edifice, but also of the diverse styles of eruptive activity. An erratum to this article is available at .     

Three distinct regimes of volcanic tremor associated with the eruption of Shishaldin Volcano, Alaska 1999

Tremor signals associated with the eruption of Shishaldin Volcano on 19 and 23 April 1999 were the strongest recorded anywhere in the Aleutian Arc by the Alaska Volcano Observatory (AVO) in its 10-year history. Reduced displacements (D_R) reached 23 cm² on 19 April and 43 cm² on 23 April. During the activity, D_R and spectral data with a frequency resolution of 0.1 Hz were computed and put on the World Wide Web every 10 min. These data are analyzed here. The general temporal patterns of seismicity of these eruption events were similar, but the eruptions and their effects quite different. The 19 April event is known to have culminated in a sub-Plinian phase, which ejected ash to an altitude of 16 km. Despite higher amplitudes and the largest hotspot from satellite data, the 23 April event produced little ash reaching only 6 km altitude. For several hours prior to the sub-Plinian phase on 19 April, tremor with a peak frequency of 1.3 Hz intensified. During the sub-Plinian phase the peak frequency increased to 4-8 Hz. However, in 15 h after the eruption, three episodes of stronger tremor occurred with a lower 1.0-Hz peak, alternating with weaker tremor with a 1.3-Hz peak. These transitions correspond to D_R=~8 cm². Although these strong tremor episodes produced higher D_R levels than the sub-Plinian phase, data from a pressure sensor show that only strong Strombolian explosions occurred. The suite of observations suggests three distinct tremor regimes that may correspond to slug flow, bubbly flow, and sustained strong eruptions, or a cyclic change in source parameters (e.g., geometry, sound speed, or ascent rate). This behavior occurred at Shishaldin only during the April 1999 sequence, and we are not aware of similar behavior at other volcanoes.     

以颗粒物理原理认识汶川MW 7.9地震发生前临震预滑和震颤现象

发生在孕震区周边地块上的临震预滑和震颤现象,对破坏性地震预测有一定前兆意义,是值得地震学界关注的问题。选取2008年5月12日汶川M_W 7.9地震发生前,临夏和湟源地震台分量应变仪记录与临夏、恩施和西安地震台数字地震仪记录以及临夏和周至地震台深井水位仪记录,分析发现,在临震前数天至数小时,上述各地震台不同学科观测仪器均记录到一些"跃变"和"震颤"震相。文中试图以颗粒物理原理,来认识不同距离、不同台站、不同学科的观测仪器在临震前相近时间段内记录的低频和高频震相,可能是不同地块在临震前发生预滑错动后激发的预滑震相Xp和地下气体在裂隙内流动激发的震颤震相Tp。观测结果表明：2008年5月8日03时至主震发生,各地震台所处地块在相近时段内逐次发生次数不等的预滑错动,其中1-2次较大错动可在噪声背景中被识别;各地震台预滑错动方向指向或背向主震震中。据此认为：汶川M_W 7.9地震前,上述各地震台所处地块在不同大小、不同方向的力链驱动下,发生指向或背向主震震中的临震预滑现象。     

以颗粒物理原理认识汶川M_W 7.9地震发生前临震预滑和震颤现象

发生在孕震区周边地块上的临震预滑和震颤现象,对破坏性地震预测有一定前兆意义,是值得地震学界关注的问题。选取2008年5月12日汶川M_W 7.9地震发生前,临夏和湟源地震台分量应变仪记录与临夏、恩施和西安地震台数字地震仪记录以及临夏和周至地震台深井水位仪记录,分析发现,在临震前数天至数小时,上述各地震台不同学科观测仪器均记录到一些"跃变"和"震颤"震相。文中试图以颗粒物理原理,来认识不同距离、不同台站、不同学科的观测仪器在临震前相近时间段内记录的低频和高频震相,可能是不同地块在临震前发生预滑错动后激发的预滑震相Xp和地下气体在裂隙内流动激发的震颤震相Tp。观测结果表明:2008年5月8日03时至主震发生,各地震台所处地块在相近时段内逐次发生次数不等的预滑错动,其中1—2次较大错动可在噪声背景中被识别;各地震台预滑错动方向指向或背向主震震中。据此认为:汶川M_W7.9地震前,上述各地震台所处地块在不同大小、不同方向的力链驱动下,发生指向或背向主震震中的临震预滑现象。     

Reventador Volcano 2005: Eruptive activity inferred from seismo-acoustic observation

Reventador Volcano entered an eruptive phase in 2005 which included a wide variety of seismic and infrasonic activity. These are described and illustrated: volcano-tectonic, harmonic tremor, drumbeats, chugging and spasmodic tremor, long period and very long period events. The recording of this simultaneous activity on an array of three broadband, seismo-acoustic instruments provides detailed information of the state of the conduit and vent during this phase of volcanic eruption. Quasi-periodic tremor at Reventador is similar to that observed at other volcanoes and may be used as an indicator of vent aperture. Variations in the vibration modes of the volcano, frequency fluctuations and rapid temporal fluctuations suggest the influx of new material, choking of the vent and possible modification of the conduit geometry during explosions and effusion over a period of six weeks.     

River flow regimes in a changing climate

Abstract

A river flow regime describes an average seasonal behaviour of flow and reflects the climatic and physiographic conditions in a basin. Differences in the regularity (stability) of the seasonal patterns reflect different dimensionality of the flow regimes, which can change subject to changes in climate conditions. The empirical orthogonal functions (EOF) approach can be used to describe the intrinsic dimension of river flow regimes and is also an adopted method for reducing the phase space in connection to climate change studies, especially in studies of nonlinear dynamic systems with preferred states. A large data set of monthly river flow for the Nordic countries has been investigated in the phase space reduced to the first few amplitude functions to trace a possible signature of climate change on the seasonal flow patterns. The probability density functions (PDF) of the weight coefficients and their possible change over time were used as an indicator of climate change. Two preferred states were identified connected to stable snowmelt-fed and rainfed flow regimes. The results indicate changes in the PDF patterns with time towards higher frequencies of rainfed regime types. The dynamics of seasonal patterns studied in terms of PDF renders it an adequate and convenient characterization, helping to avoid bias connected to flow regime classifications as well as uncertainties inferred by a modelling approach.     

Volcanic tremor at Pavlof Volcano,Alaska, October 1973–April 1986

In thirteen years (1973–1986) of seismic monitoring of Pavlof Volcano, 488 episodes of volcanic tremor have been recorded, only 26 of which have been previously described in the literature. This paper tabulates and describes all the tremor episodes and reports on the results of all analyses to date. Pavlof tremor durations range from 2 minutes to greater than 1 week; episodes accompanying magmatic eruptions have durations greater than 1 hour, and sustained amplitudes of greater than 6 mmP-P (=54 nanometers at 1.5 Hz) on station PVV, 8.5 km from the vent. Digital data provide much better amplitude resolution than helicorders do. Helicorders, however, provide continuous coverage, whereas digital data are intermittent. Correlations of tremor with visual eruption observations shows that tremor amplitudes are roughly correlated with heights of lava fountains, but the correlation of tremor amplitudes with plume heights is more problematic. Fast Fourier Transform (FFT) spectra show that Pavlof tremor is quite statinary for the entire time period, 1973–1983. All principal spectral peaks lie between 0.8 and 3.0 Hz, and may be caused by resonance of magma and gas, and resonance of the volcanic pile. Preliminary analysis of 2-and 3-component data shows thatP, S, PL, and Rayleigh waves may be present in Pavlof volcanic tremor. Other waveforms can be misidentified as tremor, most commonly those caused by storms orS-waves of regional earthquakes. A strategy is proposed to distinguish tremor from noise using automatic seismic data acquisition and analysis systems. Pavlof's volcanic tremor is briefly compared with a preliminary sample of over 1100 cases of tremor from 84 volcanoes worldwide. Finally, several recommendations for monitoring and reporting volcanic tremor are discussed.     

Geological and seismological evidence of increased explosivity during the 1986 eruptions of Pavlof volcano,Alaska

We present results of study of the best-documented eruptions of Pavlof volcano in historic time. The 1986 eruptions were mostly Strombolian in character; a strong initial phase may have been Vulcanian. The 1986 activity erupted at least 8×10⁶ m³ of feldspar-phyric basaltic andesite lava (SiO₂=53–54%), and a comparable volume of wind-borne tephra. During the course of the eruption, 5300 explosion earthquakes occurred, the largest of which was equivalent to an M _L =2.5 earthquake. Volcanic tremor was recorded for 2600 hours, and the strongest tremor was recorded out to a distance of 160 km and had an amplitude of at least 54 cm² reduced displacement. The 1986 eruptions modified the structure of the vent area for the first time in over two decades. A possible pyroclastic flow was observed on 19 June 1986, the first time such a phenomenon has been observed at the volcano. Overall, the 1986 eruptions were the strongest and longest duration eruptions in historic time, and changed a temporal pattern of activity that had persisted from 1973–1984.     

Dynamics of hourly sea level at Hillarys Boat Harbour,Western Australia: a chaos theory perspective

Water level forecasting using recorded time series can provide a local modelling capability to facilitate local proactive management practices. To this end, hourly sea water level time series are investigated. The records collected at the Hillarys Boat Harbour, Western Australia, are investigated over the period of 2000 and 2002. Two modelling techniques are employed: low-dimensional dynamic model, known as the deterministic chaos theory, and genetic programming, GP. The phase space, which describes the evolution of the behaviour of a nonlinear system in time, was reconstructed using the delay-embedding theorem suggested by Takens. The presence of chaotic signals in the data was identified by the phase space reconstruction and correlation dimension methods, and also the predictability into the future was calculated by the largest Lyapunov exponent to be 437 h or 18 days into the future. The intercomparison of results of the local prediction and GP models shows that for this site-specific dataset, the local prediction model has a slight edge over GP. However, rather than recommending one technique over another, the paper promotes a pluralistic modelling culture, whereby different techniques should be tested to gain a specific insight from each of the models. This would enable a consensus to be drawn from a set of results rather than ignoring the individual insights provided by each model.     

昆仑山强震前的震颤波并非源自慢地震

2001年11月14日的昆仑山M_s8.1级地震前几天,中国地震台网多个台站都观测到了持续数天的低频震颤波信号.由于这些震颤波发生在强震前,所以备受关注.多年来研究人员对该震颤波的产生原因进行过多方探讨,但没有定论.该震颤波信号是否源自强震区的慢地震?是否是地震前兆?或为其它因素?为了回答这些问题,我们从多方面分析和研究了昆仑山强震前中国大陆宽频地震仪所观测的震颤波信号的特征、持续时间、震颤波强度变化与大规模大气运动的关系、信号强度随观测空间的衰减变化特征.结果表明:中国大陆宽频地震仪在昆仑山强震前观测到的震颤波由两个信号组成,其中11月10日开始出现,主要频率范围0.15~0.22 Hz (周期约4~7 s)、持续时间在10-13日的震颤波,主要由同时间段内发生在西太平洋的强台风玲玲(Ling Ling)引发;而11月11日开始出现,主要频率范围0.1~0.13 Hz (周期7~10 s)、持续时间在11-12日的震颤波,不是来自昆仑山强震区的慢地震,而是由来自欧洲北部及欧亚大陆的强温带气旋引发.     

Seismological studies of the Raoul Island eruption, 1964

Raoul Island, the largest of the Kermadec Group, is about 8 km across and is situated about 1100 km north-east of Auckland, New Zealand. An eruption of steam and mud occurred on the island on 1964 November 21 (local date); this was preceded by an earthquake swarm which started on 10 November. Records made with the Willmore seismograph at the Meteorological Station on the north coast of the island show that within four hours of the start of the swarm, the frequency of shocks had risen to a peak of more than 80 an hour. Most shocks had sharp beginnings with S-P intervals of 1 to 2 seconds. Tremor was noticeable on 11 November and increased until by 12 November it was continuous, masking all but the largest of the discrete shocks. The level of tremor and the number of recorded shocks then decreased, until by the time the largest earthquake occurred on November 14, only 30 to 40 shocks were being recorded per hour. The largest earthquake was assigned a magnitude of 5.7 from recordings made in New Zealand, and was felt on Raoul Island at intensity 7 on the Modified Mercalli scale. Further shocks and tremor were recorded, and on 15 November the tremor was particularly active for several hours. On 21 November, an eruption occurred, throwing steam, mud, and rocks to an estimated height of 800 m from a crater on the edge of Green Lake, about 2 km from the Meteorological Station. During the eruption the seismograph recorded a peak vibration 30 to 40 times the amplitude of the normal background level. By two hours after the eruption, the level of vibration had stabilised at double that before the eruption. The island was evacuated from 23 November until 6 December, during which time the seismograph was inoperative. From 6 December to 11 December three portable seismographs were recording in addition to the permanent station. By this time the frequency of recorded shocks had dropped to about two per hour. The earthquakes were located in the vicinity of Denham Bay, some kilometres to the west of the main crater. Volcanic activity has been observed previously in Denham Bay, and it is thought that the Bay may be a former crater.     

Stochastic simulation of volcanic tremor from Ruapehu

The most common volcanic tremor produced by Ruapehu is a continuous signal with a dominant frequency of about 2 Hz. This signal has a sharply peaked spectrum, and an autocorrelation function with a high degree of coherence, even for lags of over 20 seconds. These characteristics strongly suggest that the cause of this tremor is a single resonator, probably a fluid-filled cavity resonating in an “organ-pipe” mode.The stochastic simulation of such a resonator uses the equation of motion of a Simple Harmonic Oscillator, which applies to an “organ-pipe” fundamental resonance, with either the characteristics of the oscillator, or the forcing function, containing a random element. A “white noise” forcing function, which would be appropriate for excitation of the cavity by a high pressure gas input, gave good agreement with the observed spectra and autocorrelation functions. Another possible model used an oscillator with a damping factor which varied randomly, and was sometimes negative, so oscillations built up, rather than decayed. This also gave a reasonable simulation of Ruapehu tremor.The third excitation model used a Poisson process, in which during each time interval there was a certain probability of applying a fixed impulse to the resonator. It was found that the impulses had to be frequent, i.e. several times a second, to match the characteristics of Ruapehu tremor.It has been suggested that tremor is composed of a succession of low-frequency (“B-type”) earthquakes. The results of this simulation show that at Ruapehu tremor could be produced by a resonator with positive feedback just sustaining oscillation, or by a resonator excited by external impulses. The most promising model for low-frequency earthquakes describes them as the result of a major external disturbance of the resonator.     

临夏地震台观测到的临震预滑和震颤震相

对临夏地震台的YRY-4分量应变仪、水位仪和地震仪记录数据分析后发现:2008年5月12日汶川MW 7.9地震前,在3种不同学科的观测记录上,在相近时间段内均记录到了预滑震相Xp和震颤震相Tp.总结Xp震相和Tp震相记录特征的基础上,试图用实验室做的小尺度黏滑实验结果来佐证和解释所记录到的临震预滑和震颤震相的物理机制....     

Source location using time‐reverse imaging

We present the chain of time‐reverse modeling, image space wavefield decomposition and several imaging conditions as a migration‐like algorithm called time‐reverse imaging. The algorithm locates subsurface sources in passive seismic data and diffractors in active data. We use elastic propagators to capitalize on the full waveforms available in multicomponent data, although an acoustic example is presented as well. For the elastic case, we perform wavefield decomposition in the image domain with spatial derivatives to calculate P and S potentials. To locate sources, the time axis is collapsed by extracting the zero‐lag of auto and cross‐correlations to return images in physical space. The impulse response of the algorithm is very dependent on acquisition geometry and needs to be evaluated with point sources before processing field data. Band‐limited data processed with these techniques image the radiation pattern of the source rather than just the location. We present several imaging conditions but we imagine others could be designed to investigate specific hypotheses concerning the nature of the source mechanism. We illustrate the flexible technique with synthetic 2D passive data examples and surface acquisition geometry specifically designed to investigate tremor type signals that are not easily identified or interpreted in the time domain.