Search

Modeling harmonic tremor recorded at volcanoes is an essential practice in improving eruption forecasting methods and warning systems. We model the conduit dynamics of quasi-periodic tremor (chugging) recorded at Karymsky, Tungurahua, and Fuego volcanoes to estimate its source characteristics. Chugging mechanisms are estimated using two theoretical models originally derived in Garcés [1997, doi: https://doi.org/10.1029/1999JB900096] and Girona et al. [2019, doi: https://doi.org/10.1029/2019JB017482], respectively. Comparisons of the conduit and fluid output parameters suggests that chugging is primarily limited to near-surface oscillations and outgassing due to an accumulation of gas between eruptive episodes. The modeled results indicate clustered release of volatiles exsolved from a deeper magma conduit region, triggered by an initial explosion. This interpretation is consistent with both infrasonic and seismic observations at each volcano.

The monitoring of seismic activity at active glacier-hosting volcanoes is challenging as volcanic and glacial earthquakes (i.e. icequakes) can have overlapping characteristics (i.e. frequencies, waveform shape and magnitude). Here we present results from the first study to target glacial activity at active ice-covered volcanoes in the Southern Chile. The primary focus so far has been on Llaima volcano, one of the largest and most active volcanoes in the region while hosting >14 km2 of glacial ice on the flanks. We use a combination of automatic multi-station event detection and waveform cross-correlation to find candidate repeating icequakes in seismic data from the permanent volcano monitoring network recorded in early 2019. We identified dozens of low magnitude families of repeating seismic events across two months, the largest of which included over 200 events. These findings are comparable to results from analysis of seismic data recorded at Llaima volcano during the same time period in 2015. The persistent, repetitive nature of these events combined with their waveform characteristics and source locations suggest they originated from multiple sub-glacial stick-slip sources around the upper flanks of the volcano. We also deployed a network of seismo-acoustic sensors at Villarrica volcano in early 2020 to record glacial activity in concurrence with the lava lake and strombolian activity at the summit. We conclude that icequakes at Llaima volcano may be more common than previously thought and has implications for how seismic data at ice-covered volcanoes may be used for assessing future volcanic and glacial hazard potential.

Innovative technological solutions are required to access and observe subglacial hydrological systems beneath glaciers and ice sheets. Wireless sensing systems can be used to collect and return data without the risk of losing data from cable breakage, which is a major obstacle when studying fast flowing glaciers and other high-strain environments. However, the performance of wireless sensors in deep and fast-moving ice has yet to be evaluated formally. We report experimental results from Cryoegg: a spherical probe that can be deployed along an ice borehole and either remain fixed in place or potentially travel through the subglacial hydrological system. The probe makes measurements in-situ and sends them back to the surface via a wireless link. Cryoegg uses very high frequency (VHF) radio to transmit data through up to 1.3 km of cold ice to a surface receiving array. It measures temperature, pressure and electrical conductivity, returning all data in real time. This transmission uses Wireless M-Bus on 169 MHz; we present a simple “radio link budget” model for its performance in cold ice and confirm its validity experimentally. Power is supplied by an internal battery with sufficient capacity for two measurements per day for up to a year, although higher reporting rates are available at the expense of battery life. Field trials were conducted in 2019 at two locations in Greenland (the EastGRIP borehole and the RESPONDER project site on Sermeq Kujalleq/Store Glacier) and on the Rhone Glacier in Switzerland. Our results from the field demonstrate Cryoegg’s utility in studying englacial channels and moulins, including estimating moulin discharge through salt dilution gauging with the instrument deployed deep within the moulin. Future iterations of the radio system will allow Cryoegg to transmit through up to 2.5 km of ice.

In terms of terrestrial animal communication, elephants are especially noteworthy because they have been documented to produce some of the loudest sounds at frequencies between 10-35 Hz. Their vocalisations, or ‘rumbles’, with fundamental frequencies in the infrasonic range (≤ 20 Hz) can have amplitudes as high as 117 dB. Here we present our efforts to evaluate the performance of the Raspberry Shake and Boom, a low-cost seismic and acoustic sensor, for identifying and monitoring the movements of African Bush Elephants (Loxodonta africana). The test area was the Adventures with Elephants elephant reserve in South Africa which includes a herd of 7 elephants (3 females, 2 males, and 2 juveniles). Within the reserve we deployed a local network of 5 Raspberry Shake and Boom units in October 2019 in order to record seismo-acoustic waves generated by the herd. The network also included other seismic and acoustic sensors of different sensitivities which were used to assess the performance Raspberry Shake and Boom units. We show that the Raspberry Shake and Boom units performed well during the deployment, with clear recordings of elephant movement and rumbles. The acoustic data also suggests that we may be able to discriminate between individual elephants due to the distinct frequencies of their rumbles. This presentation will provide general information on the potential use of low-cost sensor units for the purpose of unobtrusively monitoring vulnerable wildlife such as elephants.

Volcanoes and glaciers are both productive sources of seismic activity which may be easily confused for each other, leading to potential missed warnings or false alarms. This presents a challenge for organizations monitoring active volcanoes with glaciers on or near the edifice. Cryogenic earthquakes (i.e. icequakes) have been studied at only a few volcanoes around the world and there is a ready need to develop robust methods for efficiently differentiating them from volcanic events. Here we present results from an ongoing study of icequakes at active ice-covered volcanoes in the Southern Chilean Volcanic Zone. The primary focus of the project so far has been on seismo-acoustic data collected at Llaima volcano, one of the largest and most active volcanoes in the region. The data, recorded in 2015 and 2019, was analysed using a combination of automatic multi-station event detection and waveform cross-correlation to find candidate repeating icequakes. We identified 11 persistent families of repeating events in 2015, and over 50 families in 2019; the largest family containing over 1000 events from January to April 2019. The persistent, repetitive nature of these events combined with their waveform characteristics and source locations suggest they originated from multiple sub-glacial sources on the upper flanks of the volcano. Low levels of volcanic activity at Llaima volcano since 2009 have allowed this clear discrimination of icequake events. We are also targeting Villarrica volcano in early 2020 with a network of seismo-acoustic sensors and to record icequake activity in concurrence with the ongoing eruptive activity at the summit. Altogether, the results from this project so far suggest icequakes may be more common than previously thought and has strong implications for how seismic data at ice-covered volcanoes may be interpreted.

Ambient noise tomography is used to image Greenland's lithosphere, which passed over the Iceland plume between ∼70 and ∼40 Ma. Cross-correlations from 21 stations from GLISN seismic network were used to invert for 2-D Rayleigh wave phase velocity maps for 14 periods between 8 and 40 s. We find that Rayleigh wave phase velocities substantially vary across Greenland, with slow velocities coinciding with NW-SE trending Iceland plume track. In east Greenland the detected velocity reduction at longer periods (33–40s) reflects substantially thinned lithosphere, thermally ablated by the plume. From the east, the reduced velocities shift NW at shorter periods (12–20s), indicating shallowing of the plume-related slow anomaly. In north-central Greenland, the reduced velocities appear in the proximity of the plume ∼60 Ma, reflecting lithospheric weakening in the presence of residual heat that still persists within the lithosphere. Our results provide important new constraints on variations in the seismic velocity structure of Greenland's crust and uppermost mantle, revealing prolonged effects of the mantle plume on the overpassing craton.

Llaima Volcano in Chile is one of the largest and most active volcanoes in the southern Andes, with over 50 eruptions since the 1600s. After years of persistent degassing, Llaima most recently erupted in a series of violent Strombolian eruptions in 2007–2009. This period had few precursory signals, which highlights the need to obtain accurate magma storage information. While petrologic advancements have been made in understanding magma degassing and crystallization trends, a comprehensive seismic study has yet to be completed. Here, we present results of a receiver function survey utilizing a dense seismic array surrounding Llaima volcano. Application of H-κ stacking and common conversion point stacking techniques reveals a new Moho estimate and two structural anomalies beneath Llaima Volcano. We interpret a low velocity zone between 8 and 13 km depth as a newly imaged magma body.

Glacially- and magmatically-derived seismic events have been noted to heavily overlap in characteristics, thus there exists the potential for false-alarms or missed warnings at ice-covered volcanoes. Here we present the first study to specifically target icequakes at an ice-covered volcano in Southern Chile. Two months of broadband seismic data collected at Llaima volcano in 2015 were analyzed in order to quantify, characterize, and locate glacially-derived seismic events at one of the most active ice-covered volcanoes in the region. We find over 1,000 repeating seismic events across 11 families, the largest of which contains 397 events. Approximate locations and characteristics of the largest families lead us to conclude that these events were derived from persistent stick-slip motion along the ice-rock interface at the base of a glacier near the volcano summit. These results have implications for future seismic monitoring at Llaima volcano and other ice-covered active volcanoes in the region.

Current geophysical techniques for visualizing seismic activity employ image reconstruction methods that rely on a centralized approach for processing the raw data captured by seismic sensors. The data is either gathered manually, or relayed by expensive broadband stations, and then processed at a base station. This approach is time-consuming (weeks to months) and hazardous as the task involves manual data gathering in extreme conditions. Also, raw seismic samples are typically in the range of 16-24 bit, sampled at 50-200?Hz and transferring this high fidelity sample from large number of sensors to a centralized station results in a bottleneck due to bandwidth limitations. To avoid these issues, a new distributed method is required which processes raw seismic samples inside each node and obtains a high-resolution seismic tomography in real time. In this paper, we present a component-averaged distributed multi-resolution evolving tomography algorithm for processing data and inverting volcano tomography in the network while avoiding centralized computation and costly data collection. The algorithm is first evaluated for the correctness using a synthetic model in a CORE emulator. Later, our proposed algorithm runs using the real data obtained from Mt. St. Helens, WA, USA. The results validate that our distributed algorithm is able to obtain a satisfactory image similar to centralized computation under constraints of network resources, while distributing the computational burden to sensor nodes. We present a distributed method to obtain seismic image on wireless sensor networks.We perform in-network computing and avoid data gathering.High resolution image can be obtained using proposed multi-resolution scheme.We evaluate the method using synthetic seismic data in a CORE emulator.Our method is able to recover velocity model of Mt. St. Helens.

We present an acoustic waveform inversion technique for infrasound data to estimate volume fluxes from volcanic eruptions. Previous inversion techniques have been limited by the use of a 1-D Green's function in a free space or half space, which depends only on the source-receiver distance and neglects volcanic topography. Our method exploits full 3-D Green's functions computed by a numerical method that takes into account realistic topographic scattering. We apply this method to vulcanian eruptions at Sakurajima Volcano, Japan. Our inversion results produce excellent waveform fits to field observations and demonstrate that full 3-D Green's functions are necessary for accurate volume flux inversion. Conventional inversions without consideration of topographic propagation effects may lead to large errors in the source parameter estimate. The presented inversion technique will substantially improve the accuracy of eruption source parameter estimation (cf. mass eruption rate) during volcanic eruptions and provide critical constraints for volcanic eruption dynamics and ash dispersal forecasting for aviation safety. Application of this approach to chemical and nuclear explosions will also provide valuable source information (e.g., the amount of energy released) previously unavailable.

Several volcanic centers, active and extinct, are located within Iceland. One of them is the Hengill volcanic center, which lies on the plate boundary between the North American and the European crustal plates in southwestern Iceland. The rifting of the two plates has opened a north‐northeast‐trending system of normal faults with frequent magma intrusions. The Hengill central volcano and its transecting fissure swarm, extending 70–80 km long from the coast south of Hengill to north of Lake Thingvallavatn with an associated graben structure, form the Hengill volcanic system. The Hengill central volcano is currently active and is the main volcanic production focus of the area associated with a high‐temperature geothermal field. In Nesjavellir, in the northern part of the Hengill area, a 400 MW geothermal power plant has been in operation since 1987. Another active but less pronounced volcanic system, the Hromundartindur volcanic system, lies at the eastern edge of the Hengill system, outside the Hengill fissure swarm. The area near Mount Hromundartindur can be classified as the central volcano of this system; it is a separate focus of volcanic production with high geothermal activity. Fracture detection by shear‐wave splitting (SWS) is an exploration method of proven reliability and unique imaging power. The method is based on the fact that a shear wave propagating through rocks with stress‐aligned microcracks will split into two waves: a fast one polarized parallel to the predominant crack direction, and a slow one polarized perpendicular to it (Crampin, 1981, 1984; Babuska and Cara, 1991). Recent major studies of SWS include Peng and Ben‐Zion (2004) and Johnson et al. (2011), among many others. Two important parameters are associated with an SWS event: the polarization direction of the fast shear …

Volcano tomography provides valuable information concerning the internal structure of a volcano. During times of increased volcano activities, real-time images of the interior would allow seismologists to better understand volcanic dynamics. However, real-time tomography requires spatial coverage of the volcano utilizing large number of seismic sensors. The current broadband seismic sensors used in volcanology are expensive and difficult to install in large quantities. In this paper, we present the design and deployment of the first wireless sensor network that is capable of in-network processing of seismic signals and performing real-time volcano tomography. Specifically, our system design addresses several major challenges in real-time volcano tomography, including high data intensity and compute overhead, dynamic nature of seismic activities, and extended operational lifetime. We systematically analyze the system delay and lifetime, and present a dynamic scheme that assigns compute tasks to different system tiers to meet energy and delay constraints. Lastly, we present learned lessons from deploying the system on two volcanoes in Ecuador and Chile in 2012 and 2015. Our study shows the feasibility of achieving in-network volcanic event detection and real-time tomography using a sensor network that is two orders of magnitude less expensive than traditional seismic equipment.

Because sources of volcano infrasound are shallow or aerial, infrasound observations can provide valuable information on eruption dynamics (Fee and Matoza, 2013). In addition to extracting eruption source parameters, higher precision source localization becomes feasible for volcano infrasound, because characteristic wavelengths of volcano infrasound are short (34–340 m for 10–1 Hz sound) compared with associated seismic wavelengths. Such high‐resolution source location makes it possible to differentiate active vents in a multivent system (Ripepe et al. , 2007; Jones et al. , 2008; Johnson et al. , 2013) and to track surface flow phenomena such as pyroclastic flow or rock avalanche (Yamasato, 1997; Ripepe and Marchetti, 2002; Moran et al. , 2008; Ripepe et al. , 2009). High precision accuracy in source location can be achieved by deploying sensors in a local range (e.g.

Volcan Santiaguito (Guatemala) “breathes” with extraordinary regularity as the edifice's conduit system accumulates free gas, which periodically vents to the atmosphere. Periodic pressurization controls explosion timing, which nearly always occurs at peak inflation, as detected with tiltmeters. Tilt cycles in January 2012 reveal regular 26 ± 6 min inflation/deflation cycles corresponding to at least ~101 kg/s of gas fluxing the system. Very long period (VLP) earthquakes presage explosions and occur during cycles when inflation rates are most rapid. VLPs locate ~300 m below the vent and indicate mobilization of volatiles, which ascend at ~50 m/s. Rapid gas ascent feeds pyroclast-laden eruptions lasting several minutes and rising to ~1 km. VLPs are not observed during less rapid inflation episodes; instead, gas vents passively through the conduit producing no infrasound and no explosion. These observations intimate that steady gas exsolution and accumulation in shallow reservoirs may drive inflation cycles at open-vent silicic volcanoes.

Infrasound and high-speed imaging during a series of field-scale buried explosions suggest new details about the generation and radiation patterns of acoustic waves from volcanic eruptions. We recorded infrasound and high-speed video from a series of subsurface explosions with differing burial depths and charge sizes. Joint observations and modeling allow the extraction of acoustic energy related to the magnitude of initial ground deformation, the contribution of gas breakout, and the timing of the fallback of displaced material. The existence and relative acoustic amplitudes of these three phases depended on the size and depth of the explosion. The results motivate a conceptual model that relates successive contributions from ground acceleration, gas breakout, and spall fallback to the acoustic amplitude and waveform characteristics of buried explosions. We place the literature on infrasound signals at Santiaguito Volcano, Guatemala, and Sakurajima and Suwonosejima Volcanoes, Japan, in the context of this model.

Source mechanisms of explosive volcanic eruptions are critical for understanding magmatic plumbing systems and magma transport. Tungurahua is a large andesitic stratovolcano where seismoacoustic data has been recorded over several years. In May 2010, an energetic eruption cycle began with a midsize Vulcanian explosion followed by swarms of explosive eruptions. The five-station seismoacoustic network recorded significant seismic and infrasonic signals from the explosions. Source mechanisms of 50 explosion earthquakes associated with Vulcanian explosions during this eruptive period are investigated here. The source centroid locations and geometries of explosive signals in the 10–2 s band were determined by full-waveform moment tensor inversion. The observed waveforms are well explained by a combination of volumetric moment tensor components and a single, vertical, downward force component. The source centroids are positioned about 1.5 km below and about 320 m north of the active crater. Eigenvalue and eigenvector analysis indicates that the source geometries can be described by a subhorizontal, thin ellipsoid representing a sill-like magma accumulation. Resultant source time histories show a repetitive sequence of inflation and deflation from event to event, indicating identical source processes frequently occurred over the period. The inflation/deflation in the deep source region may be the result of crack opening. Volatile bubble growth at depth opens a pathway for gases to escape and triggers shallow explosions at the summit crater. The associated downward single force is interpreted as an exchange of linear momentum between the source and the surrounding region during the escaping gas flow.

[1] Switching between seismic-only harmonic tremor (SHT) and seismo-acoustic harmonic tremor (SAHT) has been reported at few volcanoes worldwide, but its occurrence may indicate important changes in shallow conduit conditions. Switching was simulated in a laboratory experiment in which harmonic signals were produced with a flow-driven valve and compressed air. The harmonic signals were passed through a tank of shear-thinning viscoelastic fluid, and the resulting signals were recorded. At high fluid stiffness, a stable, open conduit was produced, and the harmonic signals generated within the experimental apparatus were efficiently transmitted into the atmosphere. At lower fluid stiffness, bubbling dominated the activity, stable pathways were not generated in the fluid, and HT was not recorded in the atmosphere. These results are compared to observations of switching at Fuego volcano, Guatemala. We conclude that at intermediate magma viscosities, the development of stable degassing pathways open to the atmosphere will allow HT generated in the conduit to be transmitted into the atmosphere. Further, subtle changes in magma properties and supply rate may control whether SHT or SAHT is recorded, providing information about the state of the shallow conduit and vent at active volcanoes.

SUMMARY Volcanic explosions are accompanied by strong acoustic pressure disturbances in the atmosphere. With a proper source model, these acoustic signals provide invaluable information about volcanic explosion dynamics. Far-field solutions to volcanic infrasound radiation have been derived above a rigid half-space boundary, and a simple inversion method was developed based on the half-space model. Acoustic monopole and dipole sources were estimated simultaneously from infrasound waveforms. Stability of the inversion procedure was assessed in terms of variances of source parameters, and the procedure was reliable with at least three stations around the infrasound source. Application of this method to infrasound observations recorded at Tungurahua volcano in Ecuador successfully produced a reasonable range of source parameters with acceptable variances. Observed strong directivity of infrasound radiation from explosions at Tungurahua are successfully explained by the directivity of a dipole source model. The resultant dipole axis, in turn, shows good agreement with the opening direction of the vent at Tungurahua, which is considered to be the origin of the dipole source. The method is general and can be utilized to study any monopole, dipole or combined sources generated by explosions.

Data from a December 2008–January 2009 seismometer deployment show substantial variation in seismic signal characteristics among different sites and eruptive events at Santiaguito volcano. A station on one of the inactive domes typically records higher amplitudes and a more peaked spectrum than does another station closer to the vent on lower ground. Amplitude ratios between the stations and spectral peakedness of the station on the dome vary substantially and depend on vertical signal polarity, which is controlled by seismic source mechanism. Finite‐difference models considering two likely source mechanisms and several hypothetical subsurfaces show that the topography of the free surface and subsurface units can increase or decrease amplitudes by focusing or obstructing body waves. A homogeneous velocity structure is shown not to be a close approximation to the actual shallow subsurface, while a subsurface including a dipping reflective layer produces results that fit observations better. These models also indicate that path effects act more strongly on waves originating from an implosion than from a downward force, possibly due to differences in body and surface wave proportions and differing susceptibilities of body and surface waves to these path effects. Thus, recorded data and models both indicate that the level of wave‐field distortion associated with propagation effects is controlled by the source mechanism.

We implement an infrasound semblance technique to identify acoustic sources originating from volcanic vents and apply the technique to the generally low-amplitude infrasound (

Broadband seismic observations made near the active vent at Santiaguito provide constraints on changing pressures inside the edifice. During January 2009 we recorded episodic ultra-long period signals (ULPs, > 30 s) that were coincident with eruptive activity centered on the active El Caliente vent. These signals lasted 2–5 min and corresponded to depressurization of a volume beneath the dome. Coincident shorter-period seismic signals (> 0.5 Hz) appeared to be associated with shallow level pyroclastic emissions and dome surface motions. These shorter period signals varied significantly from event to event, unlike the ULP waveforms, which were self-similar. As opposed to the short period signals, which were recorded across the network, the ULPs were only recorded on the horizontal channels of 3 stations that were within 1.1 km of the vent. This and other factors suggest that they were a result of long period rotational motion rather than horizontal translation. ULP source locations derived from horizontal particle motions point to the northwestern part of El Caliente 150 m from the vent axis, the general locus of activity from 1967 until flank collapses in 1989/1990. We modeled the ULP signal as a Mogi source and solved for three dimensional location and equivalent source strength. Through analysis of 78 events, we identify a mean Mogi source that is 200 m west of and 250 m beneath the centre of the vent with the peak volume changes accompanying eruptions ranging from 125 to 1825 m3. These inferred volume losses coincide with the onset of gas venting and/or explosive eruptions at the surface of Caliente. Based upon these observations we develop a model to describe the pre-eruptive, eruptive and inter-eruptive physical processes occurring at the surface and inside Santiaguito.

Mount Fuji is the focus of intense study because of its potential hazard signaled by seismic, geologic and historical activity. Based on extensive seismic data recorded in the vicinity of Mt Fuji, coda quality factor () using a single scattering model hypothesis, and intrinsic and scattering quality factor ( and ) using the Multiple Lapse Time Window (MLTW) method were measured. To focus the study on the magmatic structure below Mt Fuji, the data were separated into two groups: a near-Fuji region of rays traversing an area with radius 5 km around the summit ( < 5 km), and a far-Fuji region of rays beyond a radius of 20 km around the summit ( > 20 km). This classification shows the largest discrepancy of at a range of sampling volumes corresponding to overlapped sampling depth of about 80 km. Further, the spatial division shows significant difference of and at hypocentral distance of 80 km. The large difference of in bandwidths 2–4, and 4–8 Hz indicates lithospheric heterogeneity beneath Mt Fuji with a characteristic heterogeneity scale length of about 1 km. The results have a small error range due to the large data sample, showing that all values in the near-Fuji area are greater than those of the far-Fuji area, and for both the near and far-Fuji areas is higher than at high frequencies. The and values for far-Fuji are in the range of values for typical non-volcanic areas. The values of the near-Fuji area are lower than those of other volcanic areas considered, where as values of are not. The low for the volcanic region of near-Fuji suggests that the magmatic activity, indicated by percent partial melt, in the vicinity of Mt Fuji is not as active as hot spot volcanoes, such as Kilauea, Hawaii.

Dome growth and explosive degassing are fundamental processes in the cycle of continental arc volcanism. Johnson et al. present novel optical geodetic observations to constrain the sources responsible for long-period volcanic earthquakes, which are coincident with frequent explosive eruptions at Santiaguito volcano, Guatemala. They find that acceleration in deformation of the volcanic dome, extracted from high-resolution optical image processing, is coincident with recorded long-period seismic sources. Based on these observations, abrupt mass shift of solidified domes, conduit magma or magma pads seem to be part of the mechanism responsible for generating long-period earthquakes at silicic volcanic systems. This study presents optical geodetic observations to constrain the sources responsible for long-period volcanic earthquakes, which are coincident with frequent explosive eruptions at Santiaguito Volcano, Guatemala. It is found that acceleration in deformation of the volcanic dome, extracted from the high-resolution optical image processing, is coincident with recorded long-period seismic sources. On the basis of these observations, abrupt mass shift of solidified domes, conduit magma or magma pads seem to be part of the mechanism responsible for generating long-period earthquakes at silicic volcanic systems. Dome growth and explosive degassing are fundamental processes in the cycle of continental arc volcanism1,2. Because both processes generate seismic energy, geophysical field studies of volcanic processes are often grounded in the interpretation of volcanic earthquakes3. Although previous seismic studies have provided important constraints on volcano dynamics4,5,6, such inversion results do not uniquely constrain magma source dimension and material properties. Here we report combined optical geodetic and seismic observations that robustly constrain the sources of long-period volcanic earthquakes coincident with frequent explosive eruptions at the volcano Santiaguito, in Guatemala. The acceleration of dome deformation, extracted from high-resolution optical image processing, is shown to be associated with recorded long-period seismic sources and the frequency content of seismic signals measured across a broadband network. These earthquake sources are observed as abrupt subvertical surface displacements of the dome, in which 20–50-cm uplift originates at the central vent and propagates at ∼50 m s-1 towards the 200-m-diameter periphery. Episodic shifts of the 20–80-m thick dome induce peak forces greater than 109 N and reflect surface manifestations of the volcanic long-period earthquakes, a broad class of volcano seismic activity that is poorly understood and observed at many volcanic centres worldwide7. On the basis of these observations, the abrupt mass shift of solidified domes, conduit magma or magma pads may play a part in generating long-period earthquakes at silicic volcanic systems.

We report on the evaluation of in situ measurements of eruption velocities and relative mass flux, collected by Doppler radar, as well as acoustic and infrared data, that were recorded at Stromboli volcano. Doppler radar observations were made alongside thermal, acoustic and seismic measurements, to (1) further investigate the complex waveforms of the SW crater by combining infrasonic, infrared and Doppler radar measurements, (2) establish a relationship between infrared, acoustic and Doppler radar measurements and (3) verify that all instruments observe the same behavior in terms of relative mass eruption rate. We also explore the relationship between kinetic and acoustic energy released during an eruption. Comparing the different methods to each other we are able to show that the mass erupted can be estimated either from the total reflected Doppler radar energy or from infrared observations. However, neither thermal nor reflected energy can provide a value in terms of absolute mass in kg. The erupted masses of different eruptions can only be evaluated relative to each other. Using the combined three data sets we especially focus on the eruption dynamics of the SW crater of Stromboli, namely its fluctuations in eruption strength. These pulses in one eruptive event, so called ’pulsations’ dominate more than 40% of the eruptions of SW crater. Previous models that explain pulsations to be generated by multiple consecutively exploding bubbles are supported by our combined analysis of infrasound, thermal data, particle velocities, and the reflected energy.

A detailed analysis of discrete degassing pulses, chugs, at Sangay volcano, was performed on seismic and infrasonic records to determine the physics of the conduit. Infrasonic chugging signals appear as repetitive pulses with small variations in amplitude and time lag. An automated time-domain analysis was developed to measure with high precision time intervals and amplitudes at different wave arrivals, reducing the possibility error associated with hand picking. Using this automated method, a strong positive correlation of acoustic amplitude with repose time between individual pulses on chugging signals of Sangay was found on numerous oscillating sequences. Frequency gliding of apparent harmonic frequencies generally trends from high to low frequency at Sangay, in contrast to trends at Karymsky Volcano, Russia. A new description of chugging events using wavelet transform methods, appropriate for non-stationary signals, shows subtle changes in the waveforms relate to physical processes in the volcano. A system of non-linear feedback, based on choked flow at the vent, is postulated as the most likely source of this volcanic tremor.

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.

We studied regional ![Graphic][1] at 1 Hz (![Graphic][2] ) around the Korean Peninsula based on broadband vertical component seismic records of six Incorporated Research Institutions for Seismology (IRIS) Global Seismographic Network stations and 19 Korean stations of the Korea Institute of Geoscience and Mineral Resources (KIGAM). Using 177 seismic events with magnitudes between 5.3 and 5.7 and depths less than 50 km, the reversed two station method was applied, and 94 high quality interstation paths were selected from 869 possible pairs. These results show high ![Graphic][3] paths around the Sea of Japan (East Sea) reflecting the typical oceanic structure and low ![Graphic][4] paths around northeastern China related to inactive seismicity. Assigning these path values into 193 cells around South Korea with a size of 1° by 1°, we observed that the regional ![Graphic][5] decreased gradually from east to west between 2 and 1×10-3. [1]: /embed/inline-graphic-2.gif [2]: /embed/inline-graphic-3.gif [3]: /embed/inline-graphic-4.gif [4]: /embed/inline-graphic-5.gif [5]: /embed/inline-graphic-6.gif

Seismic inversion for three-dimensional variations of velocity and attenuation are often used to delineate magma bodies in the crust and upper mantle. Problems related to spatial resolution and data noise can obscure details relevant to investigating magma chambers, and the introduction of smoothing constraints, or damping, causes blurring. Tomographic inversions for P- and S-wave velocity/attenuation are summarized including large calderas, rift zones and smaller scale subduction zone volcanoes. While results vary considerably from place to place, most anomalies are found to be in the range of ± 10% perturbation, a range often controlled by the method of smoothing or regularization imposed during analysis. At many volcanoes high velocity anomalies are observed in the shallow regions below active areas where conduits, dykes or sills are expected to be present. At other locations low velocity perturbations are seen and interpreted as magma accumulation. Resolution limitations and regularization play a significant role in determining the level of perturbation observed in tomographic studies, although there may be regions where diffuse accumulations of magma do not exhibit strong anomalies and their identification will be elusive.

The most recent eruptive cycle of Tungurahua volcano began in May 2004, and reached its highest level of activity in July 2004. This activity cycle is the last one of a series of four cycles that followed the reawakening and major eruption of Tungurahua in 1999. Between June 30 and August 12, 2004, three temporary seismic and infrasonic stations were installed on the flanks of the volcano and recorded over 2,000 degassing events. The events are classified by waveform character and include: explosion events (the vast majority, spanning three orders of pressure amplitudes at 3.5 km from the vent, 0.1–180 Pa), jetting events, and sequences of repetitive infrasonic pulses, called chugging events. Travel-time analysis of seismic first arrivals and infrasonic waves indicates that explosions start with a seismic event at a shallow depth (

Pacific plate subduction at the Aleutian–Kamchatka juncture, or corner, could be accommodated by either a large bend or a tear in the oceanic lithosphere. In this paper, we describe a number of observations which suggest that the Pacific plate terminates abruptly at the Bering transform zone (TZ). Seismicity shoals along the subduction zone from Southern Kamchatka (600 km) to relatively shallow depths near the Kamchatka–Bering Fault intersection (100–200 km). This seismicity shoaling is accompanied by an increase in the heat flow values measured on the Pacific plate. Moreover, unusual volcanic products related to adakites are erupted on Kamchatka peninsula at the juncture. Simple thermal modeling shows that a slab torn and thinner along the northern edge of the Pacific plate would be compatible with the observations. Delayed thickening of the lithosphere due to the Meiji–Hawaiian hotspot may be responsible for the required thinning. D 2004 Published by Elsevier B.V.

Eruptive activity at Stromboli Volcano was significantly elevated over background levels in May 2001. During 63 h of observation, eight vents produced, on average, 17 explosions/h with an average repose interval of 3 min. During this period, the Stromboli vents exhibited consistent seismic and acoustic signatures, based on cross-correlation cluster analysis. Dendrogram clustering, based on waveform cross-correlation, was used to illustrate the complexity of the near surface plumbing system of Stromboli's multi-vent edifice. Cross-correlations of displacement seismograms produced by explosions at specific craters, such as the Northeast crater (NEC), form dense waveform clusters with correlation coefficients between 0.96 and 0.99, while displacement waveforms from other craters, such as the Southwest crater (SWC), exhibit loose clusters with correlations between 0.88 and 0.96. The inconsistency of SWC events, as compared to the NEC, suggests that the vent system there is more heterogeneous. Cluster linkage distance between the NEC cluster and the Central crater (CC) cluster is shorter than the linkage distance between the NEC and SWC clusters, indicating that NEC and CC are more closely related. Infrasonic observations were used to locate vent explosions confirming that the clusters of events are associated with specific vents or craters. Qualitative analyses of acoustic waveforms from approximately 500 explosions reveal that impulsive acoustic signals were associated with short, mechanically simple ground displacement responses. These events may correspond to the bursting of an individual gas slug. Similar degassing mechanisms from vents in the NEC and the CC show common characteristics in their displacement waveforms.

SUMMARY The explosions of Karymsky volcano often produce signals containing a sequence of repeating pulses recorded on acoustic and seismic sensors, known as chugging. The amplitudes of these pulses correlate with the time interval between pulses. For a given measured acoustic pressure, seismic amplitudes take on arbitrary values up to a specific, empirically determined threshold. Conversely, events with a small seismic amplitude yielded acoustic waves with large variations and large-amplitude seismic events corresponded to large acoustic waves. These observations are not consistent with a source modelled by a resonating conduit. Rather, a model consisting of a sequence of discrete pulses explains the data and provides a framework for understanding the dynamics of degassing at the vent. The physical model for chugging involves a time-varying narrowing vent where gasses are released in a series of oscillations which appear to be harmonic but instead are modelled as short-term transients, or discrete pulses, suggestive of choked flow.

In the most basic seismo–acoustic studies at volcanoes, infrasound monitoring enables differentiation between sub-surface seismicity and the seismicity associated with gas release. Under optimal conditions, complicated degassing signals can be understood, relative explosion size can be assessed, and variable seismo–acoustic energy partitioning can be interpreted. The extent to which these points may be investigated depends upon the quality of the infrasonic records (a function of background wind noise, microphone sensitivity, and microphone array geometry) and the type of activity generated by the volcano (frequency of explosions, bandwidth of the signals, and coupling efficiency of the explosion to elastic energy). To illustrate the features, benefits, and limitations of infrasonic recordings at volcanoes, we showcase acoustic and seismic records from five volcanoes characterized by explosive degassing. These five volcanoes (Erebus in Antarctica, Karymsky in Russia, and Sangay, Tungurahua, and Pichincha in Ecuador) were the focus of seismo–acoustic experiments between 1997 and 2000. Each case study provides background information about the volcanic activity, an overview of visual observations during the period of monitoring, and examples of seismo–acoustic data. We discuss the benefits and utility of the infrasound study at each respective volcano. Finally, we compare the infrasound records and eruptive activity from these volcanoes with other volcanoes that have been the focus of previous seismo–acoustic experiments.

Teleseismic receiver functions (RFs) from a yearlong broadband seismological experiment in Kamchatka reveal regional variations in the Moho, anisotropy in the supra-slab mantle wedge, and, along the eastern coast, Ps converted phases from the steeply dipping slab. We analyze both radial- and transverse-component RFs in bin-averaged epicentral and backazimuthal sweeps, in order to detect Ps moveout and polarity variations diagnostic of interface depth, interface dip, and anisotropic fabric within the shallow mantle and crust. At some stations, the radial RF is overprinted by near-surface resonances, but anisotropic structure can be inferred from the transverse RF. Using forward modeling to match the observed RFs, we find Moho depth to range between 30 and 40 km across the peninsula, with a gradational crust–mantle transition beneath some stations along the eastern coast. Anisotropy beneath the Moho is required to fit the transverse RFs at most stations. Anisotropy in the lower crust is required at a minority of stations. Modeling the amplitude and backazimuthal variation of the Ps waveform suggests that an inclined axis of symmetry and 5–10% anisotropy are typical for the crust and the shallow mantle. The apparent symmetry axes of the anisotropic layers are typically trench-normal, but trench-parallel symmetry axes are found for stations APA and ESS, both at the fringes of the central Kamchatka depression. Transverse RFs from east-coast stations KRO, TUM, ZUP and PETare fit well by two anisotropic mantle layers with trench-normal symmetry axes and opposing tilts. Strong anisotropy in the supraslab mantle wedge suggests that the mantle ‘‘lithosphere’’ beneath the Kamchatka volcanic arc is actively deforming, strained either by wedge corner flow at depth or by trenchward suction of crust as the Pacific slab retreats. D 2002 Elsevier Science B.V. All rights reserved.

A balloon borne infrasound sensor circumnavigating Antarctica over the Southern Ocean shows significant broad band signals recorded on a simple microphone. A 10 hour portion of the flight traversed South American Andes, in southern Chile just north of the active Villarrica Volcano. We present signals recorded during this expedition and suggest that the high frequency signals represent transmissions from the earth, or possible bolide activity. While the source of the high frequency signals is yet unknown we provide a statistical analysis of the content during the flight. Verification of signals using chemical explosions in New Mexico demonstrate that infrasound signals are well recorded on stratospheric stations floating hundreds of km from the source. Comparison of stratospheric floating balloon platforms with ground based (IMS) infrasound stations shows considerable improvement in the noise floor, suggesting that balloon borne signals offer a high fidelity alternative to traditional acoustic installations.

Most island-arc magmatism appears to result from the lowering of the melting point of peridotite within the wedge of mantle above subducting slabs owing to the introduction of fluids from the dehydration of subducting oceanic crust. Volcanic rocks interpreted to contain a component of melt (not just a fluid) from the subducting slab itself are uncommon, but possible examples have been recognized in the Aleutian islands, Baja California, Patagonia and elsewhere. The geochemically distinctive rocks from these areas, termed 'adakites, are often associated with subducting plates that are young and warm, and therefore thought to be more prone to melting. But the subducting lithosphere in some adakite locations (such as the Aleutian islands) appears to be too old and hence too cold to melt. This implies either that our interpretation of adakite geochemistry is incorrect, or that our understanding of the tectonic context of adakites is incomplete. Here we present geochemical data from the Kamchatka peninsula and the Aleutian islands that reaffirms the slab-melt interpretation of adakites, but in the tectonic context of the exposure to mantle flow around the edge of a torn subducting plate. We conclude that adakites are likely to form whenever the edge of a subducting plate is warmed or ablated by mantle flow. The use of adakites as tracers for such plate geometry may improve our understanding of magma genesis and thermal structure in a variety of subduction-zone environments.

Frequent degassing explosions, occurring at intervals of minutes to tens of minutes, are common at many active basaltic and andesitic volcanoes worldwide. In August 1997, April 1998, and September 1998 we recorded seismic and acoustic signals generated at two andesitic volcanoes with ‘Strombolian-type’ activity. Despite variations in explosion frequency (5‐15 h 21 at Karymsky as opposed to 1‐3 h 21 at Sangay), the signatures of the explosions are remarkably similar at these two, diverse field sites. In all explosions, gas emission begins rapidly and is correlated with an impulsive acoustic pressure pulse. Seismic waveforms are emergent and begin 1‐2 s before the explosion. We classify explosion events at the two volcanoes as either short-duration (less than 1 min) simple impulses or long-duration (up to 5 min) tremor events. Many tremor events have harmonic frequency spectra and correspond to regular 1 s acoustic pulses, often audible, that sound like chugging from a locomotive. Chugging events are intermittent, suggesting that the geometry or geochemistry of the process is variable over short time scales. We attribute the 1 Hz periodic chugs to a resonant phenomenon in the upper section of the conduit. q 2000 Elsevier Science B.V. All rights reserved.

Two recent earthquake sequences near the Coso geothermal field show clear evidence of faulting along conjugate planes. We present results from analyzing an earthquake sequence occurring in 1998 and compare it with a similar sequence that occurred in 1996. The two sequences followed mainshocks that occurred on 27 November 1996 and 6 March 1998. Both mainshocks ruptured approximately colocated regions of the same fault system. Following a comparison with the background seismicity of the Coso region, we have detected evidence of stress loading within the geothermal field that appears to be in response to the 1998 earthquakes. The ML = 5.2 mainshock in the 1998 sequence occurred at 5:47 a.m. UTC and was located approximately 45 km north of the town of Ridgecrest in the Coso range. The mainshock of the 1996 sequence had an ML magnitude of 5.3. There have been no observable surface ruptures associated with either of these sequences. Though the mainshocks for both sequences were located about 900 m apart and have nearly the same local magnitudes, the sequences differ in both their temporal and spatial characteristics. An analysis of the fault-plane solutions of the mainshocks and the aftershock locations suggests that the two sequences ruptured fault planes that are perpendicular to one another. We observe a much faster temporal decay of the 1998 sequence compared to the one in 1996; moreover, while the 1996 sequence was not followed by any sizeable (i.e., ML > 4.0) aftershocks, the 1998 sequence had four such events. From an estimate of the tectonic stressing rate on the fault that produced the 1998 sequence, we infer a repeat cycle of 135 years for an earthquake of comparable magnitude at Coso.

We present results from a local earthquake tomographic imaging experiment in the greater Mount Rainier area. We inverted P wave arrival times from local earthquakes recorded at permanent and temporary Pacific Northwest Seismograph Network seismographs between 1980 and 1996. We used a method similar to that described by Lees and Crosson [1989], modified to incorporate the parameter separation method for decoupling the hypocenter and velocity problems. In the upper 7 km of the resulting model there is good correlation between velocity anomalies and surface geology. Many focal mechanisms within the St. Helens seismic zone have nodal planes parallel to the epicentral trend as well as to a north-south trending low-velocity trough, leading us to speculate that the trough represents a zone of structural weakness in which a moderate (M 6.5–7.0) earthquake could occur. In contrast, the western Rainier seismic zone does not correlate in any simple way with anomaly patterns or focal mechanism fault planes, leading us to infer that it is less likely to experience a moderate earthquake. A ∼10 km-wide low-velocity anomaly occurs 5 to 18 km beneath the summit of Mount Rainier, which we interpret to be a signal of a region composed of hot, fractured rock with possible small amounts of melt or fluid. No systematic velocity pattern is observed in association with the southern Washington Cascades conductor. A midcrustal anomaly parallels the Olympic-Wallowa lineament as well as several other geophysical trends, indicating that it may play an important role in regional tectonics.

During the summer of 1997, Karymsky Vol- cano produced summit explosions about six times each hour. Typical explosive episodes lasted between 30 seconds and three minutes, produced gas and ash columns several hun- dred meters high, and ejected some incandescent material. To better understand the physical source mechanisms re- sponsible, we recorded hundreds of explosions with a three component broad-band seismometer and microphone lo- cated 1650 meters from the active vent. Nearly every explo- sion is recorded as an emergent yet identical seismic wavelet which is followed 4.15 s later by an impulsive acoustic ar- rival. We interpret the signals as a near-surface gas volume burst which fractures the vent 'plug,' lowers the lithostatic pressure within the magma column, and often induces fur- ther degassing. When degassing continues, it is generally manifested as either a series of regular one second 'chug- ging 1 explosions, steady higher frequency 'jetting', or a hy- brid combination. We believe that the seismic signature for 'chugs,' short duration harmonic tremor with integer over- tones, is the result of repeated gas volume bursts at the vent. In contrast, seismograms for jetting are non-harmonic and contain higher frequencies. We believe that the competing degassing behaviors are influenced by the gas flux as well as the plug/conduit characteristics. We propose that a plug exists due to a viscosity gradient caused by volatile depletion in the upper conduit.

The Kamchatka peninsula in Russian East Asia lies at the junction of a transcurrent plate boundary, aligned with the western Aleutian Islands, and a steeply-dipping subduction zone with near-normal convergence. Seismicity patterns and P-wave tomography argue that subducting Pacific lithosphere terminates at the Aleutian junction, and that the downdip extension (>150km depth) of the slab edge is missing. Seismic observables of elastic anisotropy (SKS splitting and Love-Rayleigh scattering) are consistent with asthenospheric strain that rotates from trench-parallel beneath the descending slab to trench-normal beyond its edge. Present-day arc volcanism is concentrated near the slab edge, in the Klyuchevskoy and Sheveluch eruptive centers. Loss of the downdip slab edge, whether from thermo-convective or ductile instability., and subsequent "slab-window" mantle return flow is indicated by widespread Quaternary volcanism in the Sredinny Range inland of Klyuchevskoy and Sheveluch, as well as the inferred Quaternary uplift of the central Kamchatka depression. The slab beneath Klyuchevskoy has shallower dip (35°) than the subduction zone farther south (55°) suggesting a transient lofting of the slab edge, either from asthenospheric flow or the loss of downdip load. Such lofting may induce pressure-release melting to supply the Klyuchevskoy and Sheveluch eruptive centers. Petrologic indicators of high magma-peridotite equilibrium temperatures, long residence times for the hydrous arc-volcanic component, and weak expression of subducted sediment flux support the lofting hypothesis, and discourage an alternate interpretation in terms of accelerated slab rollback and/or a heightened influx of subducted volatiles. Over the late Cenozoic, the Komandorsky Basin subducted beneath northern Kamchatka and produced arc volcanics in the Sredinny Range. Several lines of evidence suggest the northeast migration of a plate triple junction (North America/Pacific/Komandorsky) along the eastern Kamchatka coast in Oligocene-Miocene times. Three "Cape terranes" (Shipunsky, Kronotsky, Kamchatka) along the coastline are exotic, with geologic similarities to present-day western Aleutian islands, and may have accreted in a "caulking-gun" process as the triple junction migrated NE. The late Cenozoic transfer of arc volcanism from the Sredinny Range to the eastern volcanic front of Kamchatka may have been facilitated by the progressive replacement of a shallow-dipping Komandorsky slab with a steeply-dipping Pacific slab.

Unzen volcano, located in southwestern Japan, erupted on November 17, 1990 after 198 years of dormancy and has been ejecting lava since May 20, 1991. In this paper, we present three-dimensional P- and S-wave velocity variations below Unzen volcano using 22,473 P- and 14,349 S-wave arrival times from 3457 local earthquakes recorded on a local network of twelve seismic stations. The model was parameterized by 24,000 approximately 2.0 km cubic blocks, targeting a volume of 120 × 100 × 20 km. A prominent low-velocity anomaly (greater than 4% slowness perturbation) starting from 2.5 to 5.0 km depth beneath the volcano is observed on both the P- and S-wave inversions. Below Unzen volcano slightly lower amplitude low-velocity anomalies are further observed to a depth of 12.5 km for P- and 7.5 km for S-waves. Shallow low-velocity anomalies are observed below the Chijiwa bay area (northwest of Unzen), and deeper anomalies are seen below Shimabara bay (southwest of Shimabara peninsula). This NW-SE-trending feature represents an elongated system of dykes which supply the shallow magma reservoirs.

[1] Ground tilt is measured from broadband seismic records prior to frequent explosions at Fuego volcano, Guatemala. We are able to resolve tilt beginning 20–30 minutes prior to explosions, followed by a rapid reversal in deformation coincident with explosion onsets. The tilt amplitude and polarity recorded on the horizontal channels vary from station to station such that the steep and unusual topography of the upper cone of Fuego appears to affect the ultra-long-period signals. We account for the effect of topography and attempt to constrain the tilt source depth and geometry through finite-difference modeling. The results indicate a shallow spherical pressure source, and that topography must be considered when attempting to model tilt sources at volcanoes with steep topography. The tilt signals are interpreted as pressurization of the shallow conduit beneath a crystallized plug followed by elastic deflation concurrent with explosive pressure release.

[1] Numerical modeling of waveform diffractions along the rim of a volcano vent shows high correlation to observed explosion signals at Karymsky Volcano, Kamchatka, Russia. The finite difference modeling assumed a gaussian source time function and an axisymmetric geometry. A clear demonstration of the significant distortion of infrasonic wavefronts was caused by diffraction at the vent rim edge. Data collected at Karymsky in 1997 and 1998 were compared to synthetic waveforms and variations of vent geometry were determined via grid search. Karymsky exhibited a wide range of variation in infrasonic waveforms, well explained by the diffraction, and modeled as changing vent geometry. Rim diffraction of volcanic infrasound is shown to be significant and must be accounted for when interpreting source physics from acoustic observations.

[1] Vertical inflation of the Caliente lava dome at Santiaguito (Guatemala) occurs coincidentally with the onset of explosive eruptions and produces infrasound that is generally peaked between 0.5 and 2 Hz with amplitude of up to 5 Pa (at ∼1 km from vent). Inflation of up to ∼1 m progresses rapidly (within ∼2 s) and can encompass the entire surface of the ∼200-m-diameter dome. We use particle image velocimetry to quantify the time history of dome uplift and demonstrate that deflection of a volcano's solid surface can generate infrasound waves with amplitudes of a few Pa. The volumetric displacement history is used to model linear volumetric acoustic sources, both compact and of finite extent. Synthetic waveforms match recorded infrasound in terms of timing and frequency content. Amplitude fit is very good for a station located 3.3 km from the vent and less good for stations shielded by near vent topography.

Fifteen earthquakes (Mw 4.1–6.4) occurring at ten major segments of the Sumatran Fault Zone (SFZ) were analyzed to identify their respective fault planes. The events were relocated in order to assess hypocenter uncertainty. Earthquake source parameters were determined from three-component local waveforms recorded by IRIS-DMC and GEOFON broadband lA networks. Epicentral distances of all stations were less than 10°. Moment tensor solutions of the events were calculated, along with simultaneous determination of centroid position. Joint analysis of hypocenter position, centroid position, and nodal planes produced clear outlines of the Sumatran fault planes. The preferable seismotectonic interpretation is that the events activated the SFZ at a depth of approximately 14–210 km, corresponding to the interplate Sumatran fault boundary. The identification of this seismic fault zone is significant to the investigation of seismic hazards in the region.

High-resolution, three-dimensional images of P-wave velocity anomalies below Mt. St. Helens, Washington, were derived using tomographic inversion. The model is a 27.5 × 21 × 20 km target volume parameterized by blocks 0.5 km per side. The area included 39 stations and 5454 local events leading to 35,475 rays used in the inversion. To diminish the effects of noisy data, the Laplacian was constrained to be small within horizontal layers, providing smoothing of the model. Non-linear effects were compensated for by iterating three-dimensional ray tracing (using pseudo-bending) between inversions and relocating earthquakes relative to the updated three-dimensional model. The structural differences between the linear and non-linear inversions appear to be insignificant, although the amplitudes of the anomalies are larger in the non-linear models. Results indicate a low-velocity anomaly (> 7%), approximately 1 km in lateral extent, from 1.5 to 3 km depths. Between 3 and 6 km depth the anomaly appears to spread out. Below 6 km depth the low-velocity feature changes to a higher-velocity perturbation with lower-velocity perturbations flanking around the perimeter of the volcano. The higher-velocity material, which correlates with the higher seismicity at that depth, is interpreted as being a plug capping the low-velocity magma chamber which begins below 9 km depth.