Long-period Events Research Articles

Time reverse location of seismic long-period events recorded on Mt Etna

We present the first application of a time reverse location method in a volcanic setting, for a family of long-period (LP) events recorded on Mt Etna. Results are compared with locations determined using a full moment tensor grid search inversion and cross-correlation method. From 2008 June 18 to July 3, 50 broad-band seismic stations were deployed on Mt Etna, Italy, in close proximity to the summit. Two families of LP events were detected with dominant spectral peaks around 0.9 Hz. The large number of stations close to the summit allowed us to locate all events in both families using a time reversal location method. The method involves taking the seismic signal, reversing it in time, and using it as a seismic source in a numerical seismic wave simulator where the reversed signals propagate through the numerical model, interfere constructively and destructively, and focus on the original source location. The source location is the computational cell with the largest displacement magnitude at the time of maximum energy current density inside the grid. Before we located the two LP families we first applied the method to two synthetic data sets and found a good fit between the time reverse location and true synthetic location for a known velocity model. The time reverse location results of the two families show a shallow seismic region close to the summit in agreement with the locations using a moment tensor full waveform inversion method and a cross-correlation location method.

Generation of Vulcanian activity and long-period seismicity at Volcán de Colima, Mexico

During the current episode, which commenced in 1998, activity at Volcán de Colima has been characterised by daily Vulcanian events, several effusive phases and a number of larger dome destroying explosions. The upper edifice comprises of an intricate array of fractures and within this system, variations in the magma ascent rate, rheology and volatile-contents complete a complexity which controls the style of activity. Subtle variations in one or more of these factors can trigger a transition. A model is presented of the Vulcanian explosion mechanism, which is reflected in the associated seismicity: first the breaching of an impermeable cap and the initial gas loss after the rupture (low-frequency signal), followed by fragmentation (high-frequency signal). In 2005, a series of larger Vulcanian explosions associated with ascending magma, represent the period of activity with the highest production rate in recent years. Pyroclastic flows were produced by column collapse, with the absence of vesicularity amongst the products pointing to a deep source of gas driving the eruption. The appearance of swarms of long-period (LP) events associated with the explosions provided a great opportunity for analysis and an insight into the processes within the fracture system, which control the eruptive style. Cross-correlation of the LP waveforms produced a series of ten families which reappeared in subsequent swarms suggesting a consistent source. Brittle fracture associated with the enhanced stain-rate found along the conduit margins is suggested as the source of this seismicity. This is supported by a linear relationship between log event rate and relative amplitude of the events within each swarm. An analysis of the temporal distribution of LP events revealed a variation in conditions between swarms. For some swarms, the conditions within the fracture system meant that a conflict between the processes associated with magma ascent was revealed, whilst for others, it was a situation of failure and backup, which can be interpreted as the failing and activation of different fractures within the upper edifice.

Long-period seismicity during magma movement at Volcán de Colima

During the period from February to September 2005, Volcan de Colima produced 30 Vulcanian explosions of sufficient magnitude to produce pyroclastic flows of variable size, with a total volume of at least 2.5 × 106 m3. Swarms of long-period events were associated with each event, their duration ranging from about 6 h to 3 days and each swarm containing up to 886 events. The characteristics of the swarms have been studied to understand the source mechanism and their relationship with the Vulcanian explosions. In total, 12,548 long-period events were analysed using various comparative and statistical methods. Patterns were not apparent in the data with no correlation between different properties of the swarms (duration, magnitude or frequency of occurrence of LP events) and the magnitude of the associated Vulcanian explosion, whether recorded by seismicity, volume of pyroclastics or altitude of the eruption column. This, along with other characteristics of the swarms, such as the continuation of the swarm after the explosion, with an increase in long-period event amplitude in some cases, suggests that the mechanism is not merely associated with the pressurization under an impermeable cap and resulting pressure differentials between adjacent volumes within the system. It is more likely that the production of long-period events is dominated by brittle fracturing on the margins of an ascending magma body. A model is proposed whereby the unloading above the ascending magma column produced by a Vulcanian explosion resulted in an increase in ascent rate, reflected in the increasing amplitude of long-period events. The results reflect the complexity of non-linear processes involved during magma ascent, degassing, crystallization and rupture of the impermeable plug during the Vulcanian process. At Volcan de Colima, as at many volcanoes, long-period events represent a useful precursor for eruptive activity. For monitoring, this paper highlights some useful analyses that can be carried out, which could illustrate certain characteristics of an eruptive episode. A preliminary model is presented of the conduit processes at work during the cyclic extrusive and explosive activity during 2005.

High precision locations of long-period events at La Fossa Crater (Vulcano Island, Italy)

Since the last eruption in 1888-90, the volcanic activity on Vulcano Island (Aeolian Archipelago, Italy) has been limited to fumarolic degassing. Fumaroles are mainly concentred at the active cone of La Fossa in the northern sector of the island and are periodically characterized by increases in temperature as well as in the amount of both CO2 and He. Seismic background activity at Vulcano is dominated by micro-seismicity originating at shallow depth (<1-1.5 km) under La Fossa cone. This seismicity is related to geothermal system processes and comprises long period (LP) events. LPs are generally considered as the resonance of a fluid-filled volume in response to a trigger. We analyzed LP events recorded during an anomalous degassing period (August-October 2006) applying a high precision technique to define the shape of the trigger source. Absolute and high precision locations suggest that LP events recorded at Vulcano during 2006 were produced by a shallow focal zone ca. 200 m long, 40 m wide and N30-40E oriented. Their occurrence is linked to magmatic fluid inputs that by modifying the hydrothermal system cause excitation of a fluid-filled cavity.

A study of tilt change recorded from July to October 2006 at the Phlegraean Fields (Naples, Italy)

The tiltmetric dataset of Phlegraean Fields area showed a discrete correlation with the volcanic dynamics, suggesting that tiltmetric monitoring is important for the surveillance of active volcanic areas. Tilt data recorded in 2006 at 2 stations belonging to the monitoring network of the Osservatorio Vesuviano (INGV, National Institute for Geophysics and Volcanology, Italy) in the Phlegraean Fields are discussed in this paper. The acquired signals have shown a strong tiltmetric inversion that took place from the end of July 2006. After correcting tilt variations to eliminate the influence of temperature (influencing 90% of the signal at OLB station, hereafter OLB) a significant value of the tilt still remains. This change is related to a local inflation episode lasting 3 months, during an unrest phase that started 2 years before. It is interesting to note that tilt amplitude is much greater at OLB than the slope of the displacement field predicted by the theoretical inflation models, but data show that this field is not homogeneous and in some areas very tilted. Moreover, in the last days before the end of tilt inversion, a low energy seismic swarm happened at about 1 km of distance from the tiltmetric station by hundreds of VT (Volcano-Tectonics) and LP (Long-Period) events.

Long period and very long period events at Mt. Etna volcano: Characteristics, variability and causality, and implications for their sources

Almost 50,000 long period (LP) events were recorded at Mt. Etna from November 2003 to May 2006. We analysed these events, as well as very long period (VLP) events which were associated with some of them. During some intervals the spectral and wavefield features of LP events remained steady, with significant changes occurring between the intervals. Based on the times of the changes, we distinguish five different sub-periods. In particular, during sub-period III (June–November 2005) the wavefield at the stations nearest to the summit area was composed of P waves. Locations for 150 LP events occurring in sub-period III, determined using radial semblance, changed, at the same time as the events' spectral features changed. It was during this sub-period that many of the LP events were associated with VLP events. Based on similarity of the waveforms, we distinguished two families of VLP events, with gradually evolving waveforms. The two families are located in slightly different places, but near the sources of the LP events. The change between the families occurred at the same time as the spectra of the LP events changed. Finally the source of the VLP events was investigated by performing complete waveform moment tensor inversion of stacks of the two families. Synthetic Green's functions for the full moment tensor were calculated for a homogeneous half-space. For both families, the source region with the highest variance reduction lies approximately beneath the active craters, at 500 m below the altitudes of the stations. The solutions for both families are very similar with sources that are between 60 and 70% isotropic. Attempts to determine deviatoric moment tensors produced consistently poorer fits. The remaining energy is poorly constrained and is likely to be noise. In conclusion, these results highlight changes in the LP and VLP events at Mt. Etna over time, and the causal relationship between them.

Source mechanism of a very-long-period event at Mt Ontake, central Japan: Response of a hydrothermal system to magma intrusion beneath the summit

A very-long-period (VLP) volcanic event, which occurred on 25 January 2007 beneath Mt Ontake, central Japan, was recorded by nearby broadband seismometers as well as the Japanese nationwide Hi-net and F-net seismic networks. The VLP waveforms show a pulse-like shape with a duration of approximately 1 min. The VLP event was accompanied by a burst sequence of long-period (LP) events. We carried out waveform inversion to determine the source mechanism and location of the VLP event. Our analysis indicates volumetric changes in an inclined crack (seismic moment of 10 14 Nm) at the VLP source, which is located at 0.6 km above sea level beneath the summit of Mt Ontake and above dikes determined from GPS data. The source-time functions of the VLP event showed a pulse-like triangular shape, indicating inflation followed by deflation of the crack at the source. We interpreted the result of the waveform inversion as follows. Gases released from the intruded magma heated an aquifer system beneath the summit of Mt Ontake. A break of unstable superheated water in a hydrothermal crack in the system caused vaporization of water and associated crack inflation, which was followed by discharge of steam from the crack resulting in its collapse. These processes caused the VLP event, which may have been a precursor of a small phreatic eruption that occurred in late March 2007.

Using WiMAX to Improve Volcano Monitoring during an Emergency: A Feasibility Study of the WEIRD System Applied to the Campi Flegrei Volcanic Area Scenario (Southern Italy)

IEEE Standard for Local and Metropolitan Area Networks (hereafter IEEE 802.16; online at http://standards.ieee.org/getieee802/802.16.html) is one of the most promising mobile and fixed broadband wireless access technologies for next-generation all-IP networks in the 3.5 GHz band (European spectrum profile). Commonly known as Universal WiMAX (worldwide inter-operability for microwave access), this access technology reaches a high bit rate and covers large areas with a single base station, making it possible to offer connectivity to end users in a cost-effective way. A further useful property of the WiMAX technology is that the transmission can be used both in line-of-sight (LOS) and non-line-of-sight (NLOS) environments, allowing highly feasible communications (WiMAX Forum 2004). Thanks to these features, IEEE 802.16 opens the way to the use of wireless technologies in the environmental monitoring of areas such as seismic and volcanic zones. Volcanic eruptions, in particular, are among the highest natural risks in the world. About 600 volcanoes have erupted in historical times, and many are located near densely populated areas. An adequate monitoring system is essential to understanding the physical processes preceding an eruption and for mitigating the risk of such an event. Active volcano surveillance is primarily based on the analysis of geophysical and geochemical parameters gathered by monitoring networks. Seismological methods are among the most useful methods for monitoring volcanoes. In fact, several types of seismic signals, e.g. , volcano-tectonic earthquakes (VT), long-period events (LP), and volcanic tremors, can occur before and during an eruption (McNutt 2005). The analysis and interpretation of these seismic signals is a very important task for volcanic eruption forecasting. WiMAX technology can be applied to provide broadband wireless access in volcano monitoring scenarios to solve some of the problems that prevent us from improving surveillance of volcanic activities during emergencies. In particular, this paper focuses on a novel …

The source of infrasound associated with long‐period events at Mount St. Helens

During the early stages of the 2004–2008 Mount St. Helens eruption, the source process that produced a sustained sequence of repetitive long‐period (LP) seismic events also produced impulsive broadband infrasonic signals in the atmosphere. To assess whether the signals could be generated simply by seismic‐acoustic coupling from the shallow LP events, we perform finite difference simulation of the seismo‐acoustic wavefield using a single numerical scheme for the elastic ground and atmosphere. The effects of topography, velocity structure, wind, and source configuration are considered. The simulations show that a shallow source buried in a homogeneous elastic solid produces a complex wave train in the atmosphere consisting of P/SV and Rayleigh wave energy converted locally along the propagation path, and acoustic energy originating from the source epicenter. Although the horizontal acoustic velocity of the latter is consistent with our data, the modeled amplitude ratios of pressure to vertical seismic velocity are too low in comparison with observations, and the characteristic differences in seismic and acoustic waveforms and spectra cannot be reproduced from a common point source. The observations therefore require a more complex source process in which the infrasonic signals are a record of only the broadband pressure excitation mechanism of the seismic LP events. The observations and numerical results can be explained by a model involving the repeated rapid pressure loss from a hydrothermal crack by venting into a shallow layer of loosely consolidated, highly permeable material. Heating by magmatic activity causes pressure to rise, periodically reaching the pressure threshold for rupture of the “valve” sealing the crack. Sudden opening of the valve generates the broadband infrasonic signal and simultaneously triggers the collapse of the crack, initiating resonance of the remaining fluid. Subtle waveform and amplitude variability of the infrasonic signals as recorded at an array 13.4 km to the NW of the volcano are attributed primarily to atmospheric boundary layer propagation effects, superimposed upon amplitude changes at the source.

GRiD MT (grid-based real-time determination of moment tensors) monitoring the long-period seismic wavefield

We have developed and implemented a new grid-based earthquake analysis system that continuously monitors long-period seismic wavefield of 20–50 s recorded by broadband seismometers. The new analysis system automatically and simultaneously determines the origin time, location and seismic moment tensor of seismic events within 3 min of their occurrence. This system has been in operation since 2003 at the Earthquake Research Institute ( http://wwweic.eri.u-tokyo.ac.jp/GRiD MT/), and the locations and origin times are usually obtained within 3 s and 20 km away from the earthquake catalog values determined by the Japan Meteorological Agency (JMA). In addition, moment tensor solutions are comparable to the network solutions manually obtained. This new system should enable us to monitor long-period seismic wavefield continuously further to help identifying long-period (or low-frequency) events which are undetectable by the conventional monitoring of short-period seismic wavefields.

An efficient method to compute the dynamic response of a fluid-filled crack

SUMMARY The vibration of a fluid-filled crack is considered to be one of the most plausible source mechanisms for the long-period events and volcanic tremors occurring around volcanoes. As a tool for the quantitative interpretation of source process of such volcanic seismic signals, we propose a method to numerically simulate the dynamic response of a fluid-filled crack. In this method, we formulate the motions of the fluid inside and the elastic solid outside of the crack, using boundary integrals in the frequency domain and solve the dynamic interactions between the fluid and the elastic solid using the point collocation method. The present method is more efficient compared with the time-domain finite difference method, which has been used in simulations of a fluid-filled crack and enables us to study the dynamics of a fluid-filled crack over a wide range of physical parameters. The method also allows us direct calculation of the attenuation quality factor of the crack resonance, which is an indispensable parameter for estimating the properties of the fluid inside the crack. The method is also designed to be flexible to many applications, which may be encountered in volcano seismology, and thus, extensions of the method to more complicated problems are promising.

Automatic estimation of optimal autoregressive filters for the analysis of volcanic seismic activity

Abstract. Long-period (LP) events observed on volcanoes provide important information for volcano monitoring and for studying the physical processes in magmatic and hydrothermal systems. Of all the methods used to analyse this kind of seismicity, autoregressive (AR) modelling is particularly valuable, as it produces precise estimations of the frequencies and quality factors of the spectral peaks that are generated by resonance effects at seismic sources and, via deconvolution of the observed record, it allows the excitation function of the resonator to be determined. However, with AR modelling methods it is difficult to determine the order of the AR filter that will yield the best model of the signal. This note presents an algorithm to overcome this problem, together with some examples of applications. The approach described uses the kurtosis (fourth order cumulant) of the deconvolved signal to provide an objective criterion for selecting the filter order. This approach allows the partial automation of the AR analysis and thus provides interesting possibilities for improving volcano monitoring methods.

Analysis of infrasonic and seismic events related to the 1998 Vulcanian eruption at Sakurajima

We present results from a detailed analysis of seismic and infrasonic data recorded over a four day period prior to the Vulcanian eruptive event at Sakurajima volcano on May 19, 1998. Nearly one hundred seismic and infrasonic events were recorded on at least one of the nine seismic–infrasonic stations located within 3 km of the crater. Four unique seismic event types are recognized based on the spectral features of seismograms, including weak seismic tremor characterized by a 5–6 Hz peak mode that later shifted to 4–5 Hz. Long-period events are characterized by a short-duration, wide spectral band signal with an emergent, high-frequency onset followed by a wave coda lasting 15–20 s and a fundamental mode of 4.2–4.4 Hz. Values of Q for long-period events range between 10 and 22 suggesting that a gas-rich fluid was involved. Explosive events are the third seismic type, characterized by a narrow spectral band signal with an impulsive high-frequency onset followed by a 20–30 second wave coda and a peak mode of 4.0–4.4 Hz. Volcano-tectonic earthquakes are the fourth seismic type. Prior to May 19, 1998, only the tremor and explosion seismic events are found to have an infrasonic component. Like seismic tremor, infrasonic tremor is typically observed as a weak background signal. Explosive infrasonic events were recorded 10–15 s after the explosive seismic events and with audible explosions prior to May 19. On May 19, high-frequency impulsive infrasonic events occurred sporadically and as swarms within hours of the eruption. These infrasonic events are observed to be coincident with swarms of long-period seismic events. Video coverage during the seismic–infrasonic experiment recorded intermittent releases of gases and ash during times when seismic and acoustic events were recorded. The sequence of seismic and infrasonic events is interpreted as representing a gas-rich fluid moving through a series of cracks and conduits beneath the active summit crater.

Effects of topography and crustal heterogeneities on the source estimation of LP event at Kilauea volcano

SUMMARY The main goal of this study is to improve the modelling of the source mechanism associated with the generation of long period (LP) signals in volcanic areas. Our intent is to evaluate the effects that detailed structural features of the volcanic models play in the generation of LP signal and the consequent retrieval of LP source characteristics. In particular, effects associated with the presence of topography and crustal heterogeneities are here studied in detail. We focus our study on a LP event observed at Kilauea volcano, Hawaii, in 2001 May. A detailed analysis of this event and its source modelling is accompanied by a set of synthetic tests, which aim to evaluate the effects of topography and the presence of low velocity shallow layers in the source region. The forward problem of Green’s function generation is solved numerically following a pseudo-spectral approach, assuming different 3-D models. The inversion is done in the frequency domain and the resulting source mechanism is represented by the sum of two time-dependent terms: a full moment tensor and a single force. Synthetic tests show how characteristic velocity structures, associated with shallow sources, may be partially responsible for the generation of the observed long-lasting ringing waveforms. When applying the inversion technique to Kilauea LP data set, inversions carried out for different crustal models led to very similar source geometries, indicating a subhorizontal cracks. On the other hand, the source time function and its duration are significantly different for different models. These results support the indication of a strong influence of crustal layering on the generation of the LP signal, while the assumption of homogeneous velocity model may bring to misleading results.

Eruption dynamics at Mount St. Helens imaged from broadband seismic waveforms: Interaction of the shallow magmatic and hydrothermal systems

The current eruption at Mount St. Helens is characterized by dome building and shallow, repetitive, long‐period (LP) earthquakes. Waveform cross‐correlation reveals remarkable similarity for a majority of the earthquakes over periods of several weeks. Stacked spectra of these events display multiple peaks between 0.5 and 2 Hz that are common to most stations. Lower‐amplitude very‐long‐period (VLP) events commonly accompany the LP events. We model the source mechanisms of LP and VLP events in the 0.5–4 s and 8–40 s bands, respectively, using data recorded in July 2005 with a 19‐station temporary broadband network. The source mechanism of the LP events includes: 1) a volumetric component modeled as resonance of a gently NNW‐dipping, steam‐filled crack located directly beneath the actively extruding part of the new dome and within 100 m of the crater floor and 2) a vertical single force attributed to movement of the overlying dome. The VLP source, which also includes volumetric and single‐force components, is 250 m deeper and NNW of the LP source, at the SW edge of the 1980s lava dome. The volumetric component points to the compression and expansion of a shallow, magma‐filled sill, which is subparallel to the hydrothermal crack imaged at the LP source, coupled with a smaller component of expansion and compression of a dike. The single‐force components are due to mass advection in the magma conduit. The location, geometry and timing of the sources suggest the VLP and LP events are caused by perturbations of a common crack system.

Temporal evolution of long-period seismicity at Etna Volcano, Italy, and its relationships with the 2004–2005 eruption

Between December 2004 and August 2005, more than 50,000 long-period events (LP) accompanied by very-long period pulses (VLP) were recorded at Mt. Etna, encompassing the effusive eruption which started in September 2004. The observed activity can be explained by the injection of a gas slug formed within the magmatic column into an overlying cavity filled by either magmatic or hydrothermal fluids, thus triggering cavity resonance. Although a large number of LP events exhibit similar waveforms before the eruption, they change significantly during and after the eruption. We study the temporal evolution of the LP-VLP activity in terms of the source movement, change of the waveforms, temporal evolution of the dominant resonance frequencies and the source Q factor and changes in the polarization of the signal. The LP source locations before and after the eruption, respectively, do not move significantly, while a slight movement of the VLP source is found. The intensity of the LP events increases after the eruption as well as their dominant frequency and Q factor, while the polarization of the signals changes from predominantly transversal to pure radial motion. Although in previous studies a link between the observed LP activity and the eruption was not found, these observations suggest that such a link was established at the latter end of the eruptive sequence, most likely as a consequence of a reestablishment of the pressure balance in the plumbing system, after it was undermined due to the discharge of large amounts of resident magma during the eruption. Based on the polarization properties of the signal and geological setting of the area, a fluid-filled crack is proposed as the most likely source geometry. The spectral analysis based on the autoregressive-models (SOMPI) is applied to the signals in order to analyse the resonance frequencies and the source Q-factors. The results suggest water and basalt at low gas volume fraction as the most likely fluids involved in the source process. Using theoretical relations for the “slow waves” radiated from the fluid-filled crack, we also estimate the crack size for both fluids, respectively.

Moment‐tensor inversion of LP events recorded on Etna in 2004 using constraints obtained from wave simulation tests

The persistent occurrence of long period (LP) events at Mt Etna became apparent with the installation of the first fixed broad‐band seismic network in late 2003. Repeating similar LP events from Nov. ‘03 to Sept. ‘04 indicate a non‐destructive source process. We perform moment tensor (MT) inversions on a stacked high S/N ratio representative LP signal, conducting a grid search for the source geometry and L2‐inversion for the source time function. Results indicate a NNW‐SSE oriented resonating sub‐vertical crack as the most probable source. This result is consistent with deformation and GPS observations. Crucial to this result are constraints imposed by detailed 3D full waveform numerical simulations in a heterogeneous tomographic model with topography, and in particular a detailed assessment of the influence of very near surface velocity structure on LP signals. Pulsating gas injection is hypothesised as the most likely LP trigger.

Broadband seismic measurements of degassing activity associated with lava effusion at Popocatépetl Volcano, Mexico

From November 1999 through July 2000, a broadband seismic experiment was carried out at Popocatépetl Volcano to record seismic activity over a wide period range (0.04–100 s). We present an overview of the seismicity recorded during this experiment and discuss results of analyses of long-period (LP) and very-long-period (VLP) seismic signals recorded at stations nearest to the crater over a four-month interval December 1999–March 2000. Three families of LP signals (Types-I, II, and III) are identified based on distinctive waveform features observed periods shorter than 1 s, periods longer than 15 s, and within the period range 0.5–2.5 s. Type-I LP events have impulsive first arrivals and exhibit a characteristic harmonic wave train with dominant periods in the 1.4–1.9 s range during the first 10 s of signal. These events are also associated with a remarkable VLP wavelet with period near 30 s. Type-II LP events represent pairs of events occurring in rapid succession and whose signatures are superimposed. These are typically marked by slowly emergent first arrivals and by a characteristic VLP wave train with dominant period near 30 s, made of two successive wavelets whose shapes are quasi-identical to those of the VLP wavelets associated with Type-I events. Type-III LP events represent the most energetic signals observed during our experiment. These have an emergent first arrival and display a harmonic signature with dominant period near 1.1 s. They are dominated by periods in the 0.25–0.35 s band and contain no significant energy at periods longer than 15 s. Hypocentral locations of the three types of LP events obtained from phase picks point to shallow seismic sources clustered at depths shallower than 2 km below the crater floor. Observed variations in volcanic eruptive activity correlate with defined LP families. Most of the observed seismicity consists of Type-I events that occur in association with 1–3-min-long degassing bursts (“exhalations”). Eruptive activity increased in intensity in February, coinciding with an increasing occurrence of Type-II LP events. Type-III events were first observed at the end of February and during March, in coincidence with the formation of a new lava dome. Vulcanian eruptions occurred in April and May. These events typically exhibit broadband signatures extending over the full period range of the sensors and lasting 30–80 min.

Source process of very-long-period events accompanying long-period signals at Cotopaxi Volcano, Ecuador

Renewed seismic activity of Cotopaxi, Ecuador, began in January 2001 with the increased number of long-period (LP) events, followed by a swarm of volcano-tectonic (VT) earthquakes in November 2001. In late June 2002, the activity of very-long-period (VLP) (2 s) events accompanying LP (0.5–1 s) signals began beneath the volcano. The VLP waveform was characterized by an impulsive signature, which was accompanied by the LP signal showing non-harmonic oscillations. We observed temporal changes of both the VLP and LP signals from the beginning until September 2003: The VLP signal gradually disappeared and the LP signal characterized by decaying harmonic oscillations became dominant. Assuming possible source geometries, we applied a waveform inversion method to the observed waveforms of the largest VLP event. Our inversion and particle motion analyses point to volumetric changes of a sub-vertical crack as the VLP source, which is located at a depth of 2–3 km beneath the northeastern flank. The spectral analysis of the decaying harmonic oscillations of LP events shows frequencies between 2.0 and 3.5 Hz, with quality factors significantly above 100. The increased VT activity and deformation data suggest an intrusion of magma beneath the volcano. A release of gases with small magma particles may have repetitively occurred due to the pressurization, which was caused by sustained bubble growth at the magma ceiling. The released particle-laden gases opened a crack above the magma system and triggered the resonance of the crack. We interpret the VLP and LP events as the gas-release process and the resonance of the crack, respectively.

Swarms of repeating long-period earthquakes at Shishaldin Volcano, Alaska, 2001–2004

During 2001–2004, a series of four periods of elevated long-period seismic activity, each lasting about 1–2 months, occurred at Shishaldin Volcano, Aleutian Islands, Alaska. The time periods are termed swarms of repeating events, reflecting an abundance of earthquakes with highly similar waveforms that indicate stable, non-destructive sources. These swarms are characterized by increased earthquake amplitudes, although the seismicity rate of one event every 0.5–5 min has remained more or less constant since Shishaldin last erupted in 1999. A method based on waveform cross-correlation is used to identify highly repetitive events, suggestive of spatially distinct source locations. The waveform analysis shows that several different families of similar events co-exist during a given swarm day, but generally only one large family dominates. A network of hydrothermal fractures may explain the events that do not belong to a dominant repeating event group, i.e. multiple sources at different locations exist next to a dominant source. The dominant waveforms exhibit systematic changes throughout each swarm, but some of these waveforms do reappear over the course of 4 years indicating repeatedly activated source locations. The choked flow model provides a plausible trigger mechanism for the repeating events observed at Shishaldin, explaining the gradual changes in waveforms over time by changes in pressure gradient across a constriction within the uppermost part of the conduit. The sustained generation of Shishaldin's long-period events may be attributed to complex dynamics of a multi-fractured hydrothermal system: the pressure gradient within the main conduit may be regulated by temporarily sealing and reopening of parallel flow pathways, by the amount of debris within the main conduit and/or by changing gas influx into the hydrothermal system. The observations suggest that Shishaldin's swarms of repeating events represent time periods during which a dominant source is activated.