combining local and remote infrasound recordings from the 2009 redoubt volcano eruption abstract the explosive phase of the 2009 redoubt volcano eruption produced predominantly short duration high amplitude infrasound signals recorded up to 4500 km away all 19 numbered explosive events were recorded at a local microphone dfr 12 km as well as at an infrasound array in fairbanks alaska i53us 547 km most with high signal to noise ratios the local microphone provides an estimate of the source parameters and comparison between the two datasets allows the unique opportunity to evaluate acoustic source term estimation at a remote array high waveform similarity between dfr and i53us occurs during much of the explosive phase due to strong stratospheric ducting permitting accurate source constraints inferred from i53us data cross correlation analysis after applying a hilbert transform to the i53us data shows how the acoustic energy has passed through a single caustic as predicted by ray theory similar to previous studies significant low frequency infrasound from redoubt recorded at i53us is coincident with high altitude ash emissions the largest events also produced considerable energy at greater than 50 s periods likely related to the initial oscillations of the volcanic plume or jet many of the explosive events have emergent onsets somewhat unusual for explosive short duration eruptions comparison of the satellite derived so2 emissions with the relative amount of acoustic energy at i53us shows a very high statistically significant correlation this study reiterates the utility of using remote infrasound arrays for detection of hazardous emissions and characterization of large volcanic eruptions and demonstrates how under typical meteorological conditions remote infrasound arrays can provide an accurate representation of the acoustic source

broadband recording of strombolian explosions and associated very long period seismic signals on mount erebus volcano ross island antarctica in december 1996 and january 1997 broadband seismometers were deployed on the summit plateau of mount erebus at radial distances of 0 7 1 4 and 1 9 km from the central crater and lava lake strombolian explosions at erebus previously have been observed to produce seismic and acoustic energy between 1 and 6 hz new observations document significant energy with spectral peaks as grave as 20 s nearly identical very long period vlp signals begin 1 5 s prior to explosions have dilatational onsets and persist for up to 150 s similar vlp waveforms were recorded at all three stations indicating that the seismograms are essentially source dominated particle motions suggest an initial depth for the vlp source of up to several hundred meters migrating deeper in the course of 15 s such explosion associated vlp signals may indicate a nondestructive lossy resonance or nonlinear fluid flow excitation within the shallow magmatic system

very long period oscillations of mount erebus volcano as the signal decays vlp scalar moments up to 5 10 11 n m exceed sp moments by an order of magnitude or more suggesting distinct though genetically related sp and vlp source mechanisms we conclude that vlp signals arise from excitation of a quasi linear resonator that is intimately associated with the conduit system and is excited by gravity and inertial forces associated with gas slug ascent eruption and magma recharge vlp signal stability across hundreds of eruptions spanning 5 years the persistence of the lava lake and the rapid posteruptive lava lake recovery indicate a stable near summit magma reservoir and vlp source process index terms 4544 oceanography physical internal and inertial waves 7280 seismology volcano seismology 8419 8414 volcanology eruption mechanisms 8419 volcanology eruption monitoring 7280 keywords strombolian very long period volcano seismology

source mechanism of very long period signals accompanying dome growth activity at merapi volcano indonesia 1 very long period vlp pulses with period of 6 7s displaying similar waveforms were identified in 1998 from broadband seismographs around the summit crater these pulses accompanied most of multiphase mp earthquakes a type of long period event locally defined at merapi volcano source mechanisms for several vlp pulses were examined by applying moment tensor inversion to the waveform data solutions were consistent with a crack striking 70 and dipping 50 sw 100m under the active dome suggest pressurized gas transport involving accumulation and sudden release of 10 60 m3 of gas in the crack over a 6s interval

new insights into the 1999 eruption of shishaldin volcano alaska based on acoustic data data collected by a pressure sensor provide new insights into the 1999 eruption of shishaldin volcano unimak island alaska on 19 april 1999 after 3 months of unrest and an extended period of low level strombolian activity shishaldin experienced a subplinian eruption ash plume to 16 km followed by several episodes of strong strombolian explosions acoustic data from the pressure sensor allow us to investigate the details of an eruption which was instrumentally well recorded but with few visual observations in the 12 h prior to the subplinian phase the pressure sensor detected a series of small repeated pulses with a constant spectral peak at 2 3 hz the amplitude and occurrence rate of the pulses both grew such that the signal became a nearly continuous hum just before the subplinian eruption this humming signal may represent gas release from rising magma the main subplinian phase was heralded by 1 the abrupt end of the humming signal 2 several pulses of low frequency sound interpreted as ash bursts and 3 a dramatic increase in seismic tremor amplitude the change in acoustic signature at this time allows us to precisely time the start of the subplinian eruption previously approximated as the time of strongest tremor increase the 50 min subplinian phase actually contained several bursts of signal each of which may represent a discrete volume of magma passing through the system following the subplinian event the pressure sensor recorded four discrete episodes of strombolian gas explosions on 19 20 april and another on 22 23 april four of the five episodes were accompanied by strong seismic tremor the fifth has not been previously recog nized and was not associated with anomalous tremor amplitudes in time series these events are similar to explosions recorded at other volcanoes but in general they are much larger with maximum amplitudes of 65 pa at 6 5 km from the vent and they have low 0 7 1 5 hz peak frequencies these large explosions occurred at rates of 3 20 per minute for 1 5 h in each episode the explosions were accompanied by a small 5 km above sea level ash plume and only minor amounts of ejecta were produced thus the explosion activity was domi nated by gas release

theory of the airborne sound field generated in a resonant magma conduit abstract explosive sources triggered inside a magma conduit may excite the conduit into acoustic resonance the acoustic field in the conduit can propagate into the atmosphere through an open vent and ensonify the overlying atmosphere the character of the airborne sound field is determined by a combination of propagation and source effects the resonance of the magmatic conduit and the diffraction of the sound field at the volcanic vent are acoustic propagation effects whereas the explosion pressure signature and the firing rate of the explosions define the source characteristics for wavelengths larger than the conduit radius only the longitudinal resonances of the magmatic conduit are relevant and the open vent radiates like a piston surrounded by an infinite baffle in this case the fluid particle velocities are directed along the axis of the conduit and the sound field may propagate into the surrounding bedrock through the conduit wall displacement induced by the fluid overpressure this coupling may produce seismic signals with banded spectra such as volcanic tremor and long period events the airborne pressure field retains the modal structure of the sound field in the magmatic conduit which contains information on the conduit geometry and geo acoustic properties of the magma the seismic wavefield is driven by the acoustic field in the magma and also contains this information although it may be filtered by propagation effects in the bedrock the theoretical sound field in the magma conduit is used to interpret seismic tremor signals recorded at mt spurr volcano alaska rapid variations in the acoustic impedance of the magma conduit terminations can create systematic changes in the tremor spectra which can be used to monitor changes in the magmatic system the results of the modelling illustrate the difference between source effects and conduit resonance as well as the value of seismoacoustic measurements in volcanic environments

relative partitioning of acoustic and seismic energy during strombolian eruptions the relative elastic energy partitioning during strombolian eruptions can be estimated from acoustic pressure and seismic velocity records we outline methods for determining acoustic and seismic energies for sensors deployed within several kilometers of erupting vents we use these energies to introduce the concept of the volcano acoustic seismic ratio g or vasr which is the ratio of elastic energy propagated through the atmosphere and into the earth eruption vasr is a physical diagnostic of explosive degassing that is appropriate for comparing eruption mechanisms at individual and between various volcanoes here we assess acoustic and seismic energies and corresponding vasr for discrete strombolian explosive events at karymsky and erebus volcanoes we attribute the relatively high and stable vasr at erebus g 8 standard deviation 41 to repeatable source conditions occurring at the surface of a persistent lava lake with accompanying strong coupling to the atmosphere lower and more variable vasr at karymsky g 0 18 in 1998 to g 1 51 in 1991 with standard deviations of 93 and 313 respectively is attributed to changing conditions within a narrow partially choked conduit variable seismoacoustic energy partitioning for karymsky as manifested by the large vasr standard deviation suggests that conduit conditions affecting vasr which include magma properties conduit obstruction or fragmentation depths can evolve both during the course of an explosion and between successive events d 2005 elsevier b v all rights reserved

dynamics of coupled diffusive and decompressive bubble growth in magmatic systems bubble growth in an ascending parcel of magma is controlled both by diffusion of oversaturated volatiles and decompression as the magma rises we have developed a numerical model which explores the processes involved in water exsolution from basaltic and rhyolitic melts rising at a constant rate from magma chamber depths of 4 and 1 km while the model does not attempt to simulate natural eruptions it sheds light on the processes which control eruptive behavior under various conditions ascent rates are defined such that a constant rate of decompression dp dt is maintained a variety of initial ascent rates are considered in the model from 1 m s to 100 m s for basalts and from a few centimeters per second to 10 m s for rhyolite at the base of the conduit the model results indicate that for any reasonable ascent rate basaltic melt degasses at a rate sufficient to keep the dissolved volatile concentration at equilibrium with the decreasing ambient pressure rhyolitic melt reaches the surface at equilibrium if its ascent rate is less than 1 m s but it can erupt with high oversaturation at greater ascent rates the latter may lead to explosive eruptions if the ascent rate of rhyolite is 10 m s or more then melt barely degasses at all in the conduit and erupts with the highest oversaturation possible for the case of slow magma rise bubble growth is limited by decompression for the case of rapid magma rise bubble growth is limited by diffusion the results of our simple model do not accurately simulate natural volcanic eruptions but suggest that subsequent more complex models may be able to simulate eruptions using the insights regarding diffusive and decompressive bubble growth processes explored in this study numerical modeling of volcanic degassing may eventually lead to better prediction of eruption timing energetics and hazards of active volcanoes

on the relation between infrasound seismicity and small pyroclastic explosions at karymsky volcano 1 ground and atmospheric propagating elastic waves are integrated with digitized video to reconstruct important parameters during the onset of eruptions at karymsky volcano in 1999 muzzle velocities which range from 10 to 55 m s for 34 analyzed explosive events are found to be very well correlated with eruption velocities synthesized from infrasonic trace data acoustic modeling proposed here invokes a linear acoustic source consistent with a 3 m piston that accelerates upwards timing of acoustic arrivals relative to video records of material emission is used to constrain fragmentation sources depths which are assumed to be coincident with infrasound generation these sources are generally shallow with an observed vertical range of about 15 m systematic changes in source depth are observed during two nights of data collection and are attributed to variations in the depth of the magma free surface relative timing of acoustic and seismic arrivals is used to identify the presence of seismicity prior to fragmentation this seismicity is emergent and low intensity but significant because it precedes the presumed fragmentation events by variable amounts up to 4 s the primary short period seismic phases appear coincident with the initiation of fragmentation and infrasound generation because of the complexity of the seismic green s function it is not possible to correlate amplitude from a specific seismic pulse with observed eruption intensity

study of seismic site effects using h v spectral ratios at arenal volcano costa rica by using data obtained with a linear array at arenal volcano we show that the h v spectral ratio method can be profitably applied to detect site effects on volcanoes similar results are obtained when calculating spectral ratios with different types of seismo volcanic signals tremor ambient noise explosion quakes lp events we compare the h v ratios with theoretical s wave transfer functions calculated using velocity models obtained from seismic refraction studies there is a good agreement when the h v ratios display sharp peaks indicating a close relationship between the ratios and the transfer function of the shallow structure furthermore the main peaks of the spectral ratios are consistent with local amplification of seismic waves observed at the corresponding frequencies

model of a strong volcanic blast and a method of estimating the mass ejected summary a hydrodynamic model of a strong volcanic blast associated with large calderaforming ash gas flow eruptions is presented we show that the mass erupted from a vent influences considerably the damping of the shock wave if the ejected mass flux is constant the shock amplitude as a function of distance will decrease slower in accordance with the inverse 4 3 power law than that of an atomic explosion the inverse cube power law from the modelling we extract a relation connecting the evective energy of a blast with the mass ejected during the blast phase the mass released during the largest explosion of the 1883 eruption of krakatau and the lateral blast energy at mount st helens on 1980 may 18 are estimated on the basis of the equation obtained

a seismically constrained mass discharge rate for the initiation of the may 18 1980 mount st helens eruption we calculate the vertical mass discharge rate from mount st helens for the first few minutes of the may 18 1980 cataclysmic eruption using a new method based on seismic constraints the observed seismic waves indicate that the seismic source is a series of single forces we model these forces as thrusts due to a combination of the momentum flux of the erupted products and the pressure of the eruptive jet the momentum discharge rate is converted to a mass discharge rate based on estimates of the velocity and jet pressure as constrained by a simple fluid dynamical model only two parameters are necessary for the calculation mach number and sound velocity the calculated mass ejected in the first 100 s is 1 6 10 11 4 6 10 11 kg since the total blast deposit is 3 2 10 11 4 1 10 11 kg one possible interpretation is that the directed blast had a significant 30 vertical component

seismicity associated with eruptive activity at langila volcano papua new guinea abstract during periods of eruptive activity at langila volcano various combinations of low frequency earthquakes and high frequency air waves have been recorded ranging from earthquakes alone to air waves alone it is shown that these different types of events are part of a continuous range of events which are all interpreted to be explosions at or near the summit continuous tremor is also often recorded for durations of several minutes to hours and is classified into two types harmonic and non harmonic spectral analysis shows that the non harmonic tremor resembles a string of low frequency events in contrast the spectra of harmonic tremor show simple modal characteristics which indicate some type of resonating source this source is thought to include a gas filled cavity above the magma body within the volcano

highlights from a seismic broadband array on stromboli volcano an array of nine three component broadband seismometers was deployed in two different configurations on stromboli volcano the analysis of the seismic wavefield related to volcanic explosions revealed some observations which offer a completely new insight into the internal dynamics of a volcano these new observations are restricted to the low frequency range below 1 hz and underline therefore the superiority of broadband recordings over conventional short period observations surprisingly simple wavelets indicate an initially contracting source mechanism gas jets that could not be seen in a short period seismic record at all generate a clear dilatational wavelet in a broadband recording suggesting the same contracting source mechanism the analysis of particle motion and seismic array techniques permits a location of the seismic source we find low frequency signals of 3s and 6s period that are not related to eruptions and do not share a common source with the eruption related events a video recording of visible volcanic activity at the crater region allows one to correlate precisely eruptive features with seismic signals

mineral equilibrium modelling and calculated chemical potential relations of reaction textures in the ultrahigh temperature in ouzzal terrane in hihaou area western hoggar algeria

deep volcanic tremor and magma ascent mechanism under kilauea hawaii deep harmonic tremor originating at depths around 40 km under kilauea was studied using records accumulated since 1962 at the hawaii volcano observatory of the u s geological survey the deep source of the tremor was determined by onset times and confirmed by the relative amplitude across the island wide network of seismometers the period of tremor was conclusively shown to be determined by the source effect and not by the path or station site effect because the period would change considerably in time but maintained uniformity across the seismic net during the tremor episode the tremor appeared to be primarily composed of p waves we interpret the observed period and amplitude in terms of the stationary crack model of aki et al 1977 and find that the seismic moment rates for deep tremors are considerably larger than those for shallow tremors suggesting more vigorous transport for the former we propose a kinematic source model which may be more appropriate for deep tremor according to this model a measurable quantity called reduced displacement is directly proportional to the rate of magma flow a systematic search for deep tremor episodes was made for the period from 1962 through 1979 and the amplitude period and duration of the tremor were tabulated we then constructed a cumulative reduced displacement plot over the 18 year period the result shows a generally steady process which does not seem to be significantly affected by major eruptions and large earthquakes near the surface the total magma flow estimated from the reduced displacement is however one order of magnitude smaller than that estimated by swanson 1972 it may be that most channels transport magma aseismically and only those with strong barriers generate tremor

long range acoustic observations of the eyjafjallajokull eruption iceland april may 2010 1 the april may 2010 summit eruption of eyjafjallajokull iceland was recorded by 14 atmospheric infrasound sensor arrays at ranges between 1 700 and 3 700 km indicating that infrasound from modest size eruptions can propagate for thousands of kilometers in atmospheric waveguides although variations in both atmospheric propagation conditions and background noise levels at the sensors generate fluctuations in signal to noise ratios and signal detectability array processing techniques successfully discriminate between volcanic infrasound and ambient coherent and incoherent noise the current global infrasound network is significantly more dense and sensitive than any previously operated network and signals from large volcanic explosions are routinely recorded because volcanic infrasound is generated during the explosive release of fluid into the atmosphere it is a strong indicator that an eruption has occurred therefore long range infrasonic monitoring may aid volcanic explosion detection by complementing other monitoring technologies especially in remote regions with sparse ground based instrument networks

recent eruptions at bezymianny volcano a seismological comparison abstract for the past few decades bezymianny volcano has erupted once to twice per year here i examine eight eruptive events between 2006 and 2010 this is the first time period for which proximal or broadband seismic data have been recorded at bezymianny several recurring patterns are demonstrated in advance of eruptions eruptions are generally preceded by 12 36 h of tremor energy elevated by 2 to 3 orders of magnitude locatable earthquake activity is quite erratic in the days before eruptions for eruptions of juvenile magma however the cumulative moment magnitude increases with the repose time since the previous eruption though tenuous this relationship is statistically significant and could improve forecasts of bezymianny eruptions the most energetic eruptions demonstrate increasing multiplet activity in the run up followed by a rapid cessation at the time of eruption when present this behavior marks increasing pressure in the conduit system as degassing eclipses the capacity for venting very long period seismicity 20 s periods accompanies some eruptions these tend to be the same short lived high energy eruptions that exhibit multiplet precursors four eruptions are examined in detail to illustrate the variety in eruption mechanisms lava dome collapses sustained eruptions singular paroxysmal explosions and post explosion lava flows occur in different combinations demonstrating that more than one eruption trigger is regulating bezymianny compared to bezymianny s fifty year modern history recent eruptions have been shorter lived and separated by longer repose times some evidence suggests that these eruptions may be increasingly explosive a speculation that carries significant hazard implications if true however this threat is tempered by solid evidence that the most explosive eruptions are preceded by the clearest precursors suggesting an ability to improve the already excellent eruption forecasts available for bezymianny

canonical model of volcano acoustics a full wave theoretical model is developed for the acoustic field generated by an explosive point source embedded in a magma column that is open to the atmosphere the green s functions for the field in the magma and the atmosphere are derived on the basis of several simplifying assumptions concerning the geometry of the conduit the boundary conditions and the geoacoustic properties of the magma a kinematic model for the acoustic signature of the explosive source is proposed which when combined with the green s functions provides full analytical expressions for the acoustic field in the magma and in the atmosphere as complex functions of frequency by fourier inversion the airborne pressure spectrum is transformed into a pressure time series the predicted sound pulse in the atmosphere and its energy spectrum are highly dispersive showing complicated structure that arises from the coherent addition of many normal modes of oscillation i e depth and radial resonances in the magma column at high frequencies for which the aperture of the vent is many wavelengths across each mode is launched into the atmosphere as a parallel sided beam of sound with a characteristic angle of elevation which through snell s law is determined by the speed of sound in the magma relative to that in air at somewhat lower frequencies the modal beams of sound undergo angular spreading due to diffraction at the edge of the vent in the lowest frequency regime where the wavelength is comparable with the aperture the airborne field shows little angular structure a comparison between airborne acoustic data that we recorded in july 1994 at the western vent of stromboli volcano and the predictions of the theory using parameters that are characteristic of stromboli show compelling agreement the theoretical and observed power spectra both display the following features 1 a concentration of energy below 20 hz associated with the first four longitudinal resonances 2 radial resonances between 35 and 65 hz and 3 a broad minimum around 30 hz arising because the source lies near nulls in longitudinal modes that would otherwise be excited the conclusion is that the airborne sound signature from an explosive volcanic event may be inverted to provide estimates of the depth and radius of the magma conduit the depth spectral shape and peak shock wave pressure of the source and the viscosity of the magma

very long period signals associated with vulcanian explosions at popocatepetl volcano mexico very long period vlp seismic signals associated with large summit eruptions of popocatepetl volcano in the last four years are investigated using data from a broadband seismometer 0 04 120 s deployed on the north flank of the volcano at a radial distance of 5 km from the active crater the vlp signals associated with individual eruptions share similar waveforms discrete vlp signals accompanying long period lp events also share similar signatures and have dominant periods that are nearly identical to those observed in the vlp waveforms of explosions the vlp particle motions for eruption onsets consistently point to a source located a few km beneath the crater the vlp ground displacement response to each explosion is marked by a compression followed by a dilatation and terminating with another compression suggesting a sequence of pressurization depressurization repressurization of the conduit the repetitive nature of the waveforms points to a non destructive source process which has remained active in the magmatic system of popocatepetl at least since april 1997

degassing explosions at karymsky volcano kamchatka 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

explosion mechanisms at arenal volcano costa rica an interpretation from integration of seismic and doppler radar data 1 we execute an integrated analysis of broadband seismic and doppler radar data to gain insights into the subsurface mechanisms that drive repetitive mildly explosive activity of arenal volcano costa rica we find large variability of both seismic and radar waveforms and nonsystematic relationships between the two seismic recordings display long lasting tremor sequences and numerous explosion quakes radar measurements show that tephra emissions are poorly correlated in both time and energy to the seismic activity tephra emissions were found in association with explosion quakes but also during episodes of tremor and seismic quiescence moreover the exit velocity mass loading and kinetic energy of the emissions show no clear relationship with the coeval seismic amplitude and frequency content we propose a conceptual source model whereby degasing is controlled by opening and closing of fractures that crosscut a rigid cap atop the conduit when the fracture s strength is overcome by the gas pressure building up below it suddenly opens and high velocity gas escapes producing high frequency elastic waves typical of explosion quakes gas release also occurs in relation to periodic opening and closure of the fractures to produce repetitive pressure pulses this being the source of tremor in both cases varying quantities of fragmented material may be carried by the gas which can be detected by the radar if their concentration is high enough moreover the highly variable constantly changing state of lava cap e g thickness fracture network and gas permeability results in nonrepeatable source conditions and explains the complex relationship between tephra emissions and associated seismic signals

the effects of a decompression on seismic parameter profiles in a gas charged magma seismic velocities in a gas charged magma vary with depth and time relationships between pressure density exsolved gas content and seismic velocity are derived and used in conjunction with expressions describing diffusive bubble growth to find a series of velocity profiles which depend on time an equilibrium solution is obtained by considering a column of magma in which the gas distribution corresponds to the magmastatic pressure profile with depth decompression events of various sizes are simulated and the resulting disequilibrium between the gas pressure and magmastatic pressure leads to bubble growth and therefore to a change of seismic velocity and density with time bubble growth stops when the system reaches a new equilibrium the corresponding volume increase is accommodated by accelerating the magma column upwards and an extrusion of lava a timescale for the system to return to equilibrium can be obtained the effect of changes in magma viscosity and bubble number density is examined

seismic monitoring and modeling of an active volcano for prediction we describe seismic monitoring and modeling of an active volcano for the purpose of predicting the magma transport process in real time we selected the piton de la fournaise which has been monitored by an observatory since the early 1980s and is active with 28 eruptions from 1980 until august 1992 it erupted on march 9 1998 after an unusually long quiet period the eruption lasted until september 21 1998 before the march eruption we completed the initial construction of a model of a magma system based primarily on long period lp events and coda localization we found that lp events are associated with the lateral movement of magma to the rift zone and not with the vertical movement to the summit we also found that the source of 2 hz lp events is located below that of 1 hz lp events the coda localization is a newly discovered phenomenon that helped to locate a magma body this model was used for interpreting the incoming information from the precursory seismic crisis and the eruption tremor during the 6 months of the eruption the simultaneous monitoring and modeling of an active volcano led us to a quantitative simulation of the eruption process computing excess pressure and flow rate in a system of reservoirs connected by channels

observations of volcanic tremor at mount st helens volcano digital recordings of ground motion during tremor episodes accompanying eruptions at mount st helens volcano in the state of washington on august 7 and october 16 18 1980 are studied the spectra of the vertical component waveforms contain at least two dominant peaks at 1 0 and 1 3 hz for all events recorded during both eruptions that were studied spectra of horizontal ground motion show peaks at 0 9 and 1 1 hz the relative amplitude of the two peaks changes between tremor episodes and during single tremor episodes and shows no consistent relation to amplitude of ground motion spectra of long period earthquakes are very similar to those of tremor events suggesting that tremor is composed of many long period earthquakes that occur over a period of time the unique waveform of tremor events observed at mount st helens must be due to a source effect since the relative amplitude of the two dominant peaks changes during tremor episodes the path effect on tremor waveforms is small since there are no peaks in the spectra of waveforms recorded during tectonic earthquakes occurring in the vicinity of mount st helens the consistency of the location of the spectral peaks for the wide range of tremor amplitudes means that there must be a physical length at the source that is constant independent of the amplitude of motion at the source variations in amplitude of motion generated during tremor events must be due to variations in the force driving the tremor amplitudes of ground motion varies between 0 11 and 4 7 m seismic moment rates during the two eruptions are calculated using the model of aki et al 1977 and found to vary between 6 1018 and 1 1020 dynes cm s 1 which are larger than values found by aki et al 1977 who studied amplitudes of shallow tremor events recorded during the october 1963 eruption of kilauea volcano in hawaii study of tremor amplitudes recorded at corvallis oregon leads to the conclusion that tremor accompanying the cataclysmic may 18 1980 eruption was at least one order of magnitude larger in amplitude than tremor during august and october

excitation mechanism of atmospheric pressure waves from the 1980 mount st helens eruption summary atmospheric pressure waves from the mount st helens eruption 1980 may 18 have been clearly recorded by a sensitive microbarograph at berkeley california the record shows three types of waves with different group velocities the pressure waves can be interpreted in terms of direct waves a1 antipodean travelling waves a2 and circumnavigating waves a3 all of which are composed of several acoustic gravity modes propagated in the lower atmosphere synthetic barograms appropriate to the berkeley station have been calculated on the basis of the dynamic response of the lower atmospheric structure together with various assumptions of source properties comparisons between synthetic and observed barograms provide estimates for ranges of the time history of upward particle velocity at the source source dimensions and the velocity of the source spreading over the blast zone as well as for the average dissipation effects over the circumferential path the results suggest that two major compression pulses on the a1 record correlate with the arrival of pressure waves from the first lateral blast and second vertical blast although the inferred interblast time interval is not consistent with that estimated from seismic observations

statistical evidence of transitioning open vent activity towards a paroxysmal period at volcan santiaguito guatemala during 2014 2018 abstract long term eruptive activity at the santiaguito lava dome complex guatemala is characterised by the regular occurrence of small to moderate size explosions from the active caliente dome between november 2014 and december 2018 we deployed a seismo acoustic network at the volcano which recorded several changes in the style of eruption including a period of elevated explosive activity in 2016 here we use a new catalogue of explosions to characterise changes in the eruptive regime during the study period we identify four different phases of activity based on changes in the frequency and magnitude of explosions at the two ends of the spectrum of repose times we find pairs of explosions with near identical seismic and acoustic waveforms recorded within 1 10 min of one another and larger explosions with recurrence times on the order of days to weeks the magnitude frequency relationship for explosions at santiaguito is well described by a power law we show that changes in b value between eruptive regimes reflect temporal and spatial changes in rupture mechanisms likely controlled by variable magma properties we also demonstrate that the distribution of inter explosion repose times between and within phases is well represented by a poissonian process the poissonian distribution describing repose times changes between and within phases as the source dynamics evolve we find that changes in source properties restrict the extrapolation of explosive behaviour to within a given eruptive phase limiting the potential for long term assessments of anticipated eruptive behaviour at santiaguito

high frequency analysis of seismic background noise as a function of wind speed and shallow depth we used a deep 1500 m cased borehole near the town of datil in west central new mexico to study high frequency 1 hz seismic noise characteristics the remote site had very low levels of cultural noise but strong winds winter and spring made the site an excellent candidate to study the effects of wind noise on seismograms along with a three component set of surface sensors teledyne geotech gs 13 a vertical borehole seismometer gs 28 was deployed at a variety of depths 5 43 and 85 m to investigate signal and noise variations wind speed was measured with an anemometer event triggered and time triggered data streams were recorded on a reftek 72 02 data acquisition system located at the site our data show little cultural noise and a strong correlation between wind speed and seismic background noise the minimum wind speed at which the seismic background noise appears to be influenced varies with depth 3 m sec at the surface 3 5 m sec at 43 m in depth and 4 m sec at 85 m in depth for wind speed below 3 to 4 m sec we observe omni directional background noise that is coherent at frequencies below 15 hz this coherence is destroyed when wind speeds exceed 3 to 4 m sec we use a test event md 1 6 and superimposed noise to investigate signal to noise ratio snr improvement with sensor depth for the low q valley fill of the datil borehole dbh site we have found that snr can be improved by as much as 20 to 40 db between 23 and 55 hz and 10 to 20 db between 10 and 20 hz by deploying at a 43 m depth rather than at the surface at the surface there is little signal above noise in the 23 to 55 hz frequency band for wind speeds greater than 8 m sec thus high frequency signal information that is lost at the surface can be recorded by deploying at the relatively shallow depth of 40 m because we observe only minor further reductions in seismic background noise sbn at deeper depths 40 m is likely to be a reasonable deployment depth for other high frequency monitoring sites in similar environmental and geologic conditions

origin of the sound generated by strombolian explosions during strombolian eruptions large bubbles break at the surface of the lava column and produce sound acoustic pressure recorded during several explosions on the eastern vent of stromboli volcano shows a pattern consistent between explosions the well marked oscillation contains only very low frequencies around 7 hz and is followed by a signal containing now both higher frequencies and lower frequency less than 4 hz the generation of sound is explained by the vibration 7 hz and bursting higher frequencies of large bubbles followed by the drainage of magma down in the conduit less than 4 hz the interpretation suggests that 1 the bubble vibrates before bursting and 2 the vibration prior to bursting yields more energy than the bursting itself the radius of bubbles responsible for explosions is determined to be close to 1 meter in agreement with observations of a few conduit sizes

image processing of explosive activity at stromboli abstract a new way to study the dynamics of explosive behaviour at different vents of stromboli volcano is proposed the method is based on computer processing of images taken by an infrared film camera at a rate of 4 pictures per second each picture has been digitized by use of an image scanner with a 75 d p i dots per inch resolution and 256 levels of gray then processed and the total mass of the ejected material has been calculated for each photograph with this procedure the flux of ejected matter throughout every single explosive event has been reconstructed the analysis of trajectories made on adequate samples of ejecta has allowed us to evaluate the ejection velocity the mode of the mass flux distribution and the mean value of the velocities calculated for individual particles quantify the kinetic energy of an explosion vents nos 1 and 3 which remained active after the last effusive eruption of 1985 1986 show several explosive patterns the energy of which is measured the observed flux distributions have been checked both numerically and experimentally by modelling impulsive events the computer simulations of the mass distributions have allowed us to speculate on the possible explosive dynamisms two different models have been used a single impulse model confirms the behaviour actually shown by one of the two vents while the behaviour of the other is better depicted by a cumulative trend of successive impulses in rhythmic sequence moreover through the experimental approach we have had the opportunity to test the overall resolution of the method the observed matter flux distribution has then been associated to the seismic signals in order to correlate the dynamics involved in the explosions with the seismic source mechanisms

seismic and acoustic observations at mount erebus volcano ross island antarctica 1994 1998 volcanic activity at mount erebus is dominated by eruptive activity within a phonolitic summit lava lake common eruption styles range from passive degassing to strombolian explosions which typically occur several times daily and occasionally in swarms of up to 900 per day shallow explosions although generally the result of steady exsolution of volatiles from depth can be triggered by surficial input of h 2o through mass wasting of rock snow and ice from the crater walls broadband observations of strombolian explosions document very long period vlp signals with strong spectral peaks near 20 12 and 7 s which are polarized in the vertical radial plane these signals precede lava lake surface explosions by 1 5 s are highly repeatable and persist for up to 200 s first motions indicate a deflationary source with any precursory inflation being below the 30 s passband of our instruments particle motions suggest a vlp source residing up to 800 m below the lava lake surface however this depth could be exaggerated by near field radial tilt seismic and acoustic signals associated with lava lake explosions commonly show evidence for multiple bubble bursts in corresponding complexity features resulting from varying time delays and relative sizes of superimposed bursts a systematic decrease in seismic acoustic ratio for smaller surface explosions suggests that either the seismic energy from the smallest shallowest bubble bursts experiences much greater seismic attenuation than energy arising from larger events which may involve a deeper less attenuative portion of the magma column and or that the shallowest layer is seismically isolated from deeper parts of a stratified magma column which are not excited by the smallest explosions due to sharp impedance contrasts across distinct layers tremor at erebus is uncommon with only a few isolated instances identified in five years of monitoring some tremor events are nearly monochromatic and some exhibit numerous gliding harmonic spectral lines q 2000 elsevier science b v all rights reserved

strombolian explosions 2 eruption dynamics determined from acoustic measurements strombolian activity consists of a series of explosions due to the breaking of a large overpressurized bubble at the surface of the magma column acoustic pressure due to sound waves has been measured and analyzed at stromboli for more than 50 explosions three parts can be distinguished in the acoustic pressure waveform which are related to the behavior of the bubble before during and after its bursting before the sharp rise in acoustic pressure the signal is dominated by waves with a frequency of 2 hz which develop on the nose of the bubble they produce sound in air by imposing a rapid motion to the interface and one could detect a bubble travelling in the uppermost 30 m of the magma column when the bubble reaches the air magma interface its strong vibration driven by a large overpressure inside the gas generates the main event with a frequency around 9 hz after the bubble has burst kinematic waves of frequency around 4 5 hz are the main source of sound they develop at the surface of the magma left on the conduit side the three types of motion although determined independently give consistent results furthermore combining the results obtained for the two types of kinematic waves the magma viscosity is estimated to be of 300 65 pa s which is in good agreement with petrological constraints and corroborates the validity of our analysis this suggests that acoustic measurements constitute a powerful tool in the understanding of eruption dynamics

interpretation and utility of infrasonic records from erupting volcanoes abstract 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

air pressure waves from mount st helens eruptions weather station barograph records as well as infrasonic recordings of the pressure wave from the mount st helens eruption of may 18 1980 have been used to estimate an equivalent explosion airblast yield for this event pressure amplitude versus distance patterns in various directions compared with patterns from other large explosions such as atmospheric nuclear tests the krakatoa eruption and the tunguska comet impact indicate that the wave came from an explosion equivalent of a few megatons of tnt the extent of tree blowdown is considerably greater than could be expected from such an explosion and the observed forest damage is attributed to outflow of volcanic material the pressure time signature obtained at toledo washington showed a long 13 min duration negative phase as well as a second hour long compression phase both probably caused by ejacta dynamics rather than standard explosion wave phenomenology the peculiar audibility pattern with the blast being heard only at ranges beyond about 100 km is explicable by finite amplitude propagation effects near the source compression was slow taking more than a second but probably less than 5 s so that it went unnoticed by human ears and susceptible buildings were not damaged there was no damage asmore toledo 54 km where the recorded amplitude would have broken windows with a fast compression an explanation is that wave emissions at high elevation angles traveled to the upper stratosphere where low ambient air pressures caused this energetic pressure oscillation to form a shock wave with rapid nearly instantaneous compression atmospheric refraction then returned part of this wave to ground level at long ranges where the fast compressions were clearly audible copyright american geophysical union 1987 less

acoustic measurements of the 1999 basaltic eruption of shishaldin volcano alaska 1 origin of strombolian activity the 1999 basaltic eruption of shishaldin volcano alaska usa displayed both classical strombolian activity and an explosive subplinian plume strombolian activity at shishaldin occurred in two major phases following the subplinian activity in this paper we use acoustic measurements to interpret the strombolian activity acoustic measurements of the two strombolian phases show a series of explosions that are modeled by the vibration of a large overpressurised cylindrical bubble at the top of the magma column results show that the bubble does not burst at its maximum radius as expected if the liquid film is stretched beyond its elasticity but bursting occurs after one cycle of vibration as a consequence of an instability of the air magma interface close to the bubble minimum radius during each strombolian period estimates of bubble length and overpressure are calculated using an alternate method based on acoustic power we estimate gas velocity to be 30 60 m s in very good agreement with synthetic waveforms although there is some variation within these parameters bubble length and overpressure for the first strombolian phase are found to be c82f11 m and 0 083 mpa for the second strombolian phase bubble length and overpressure are estimated at 24f12 m and 0 15 mpa for the first 17 h after which bubble overpressure shows a constant increase reaching a peak of 1 4 mpa just prior to the end of the second strombolian phase this peak suggests that at the time the magma in the conduit may contain a relatively large concentration of small bubbles maximum total gas volume and gas fluxes at the surface are estimated to be 3 3 10 7 and 2 9 10 3 m 3 s for the first phase and 1 0 10 8 and 2 2 10 3 m 3 s for the second phase this gives a mass flux of 1 2 10 3 and 8 7 10 2 kg s respectively for the first and the second strombolian phases

permanent tremor of masaya volcano nicaragua wave field analysis and source location the masaya volcano nicaragua is a basaltic caldera in a subduction zone the permanent source of the volcanic tremor was located inside santiago crater at the lava lake s position and 400 m below the ne rim and therefore corresponds to superficial magma activity we used two tripartite arrays 90 m side one semicircular array r 120 m in 1992 and two semicircular arrays r 60 m and a 2500 m long linear array radiating out from the source and on the flank of the crater in 1993 we used both a cross spectrum method and a correlation method to determine the wave delay time between the reference station and the other stations of an array and to quantify the wave field using the delays therefore by intersecting the back azimuth wave directions from the arrays we could pinpoint the source additionally the correlation coefficients obtained as functions of frequency for the three components of motion confirm the inferred position of the source of tremor the tremor s wave field is composed of comparable quantities of dispersed rayleigh and love surface waves whose phase velocities lie in the ranges 730 1240 m s at 2 hz and 330 550 m s at 6 hz the dispersive phase velocities were inverted to obtain crustal structures with a minimal number of layers the resulting velocity models are similar for the northern and southern parts of the volcano after geometrical spreading corrections q2hz 14 and q3hz 31 were determined along the northern linear array the typical low velocities and low q corresponding to the cone structure and are similar to those of other basaltic volcanoes like puu oo hawaii and klyuchevskoy kamchatka

volcanic eruptions observed with infrasound 1 infrasonic airwaves produced by active volcanoes provide valuable insight into the eruption dynamics because the infrasonic pressure field may be directly associated with the flux rate of gas released at a volcanic vent infrasound also enhances the efficacy of volcanic hazard monitoring and continuous studies of conduit processes here we present new results from erebus fuego and villarrica volcanoes highlighting uses of infrasound for constraining quantitative eruption parameters such as eruption duration source mechanism and explosive gas flux

impacts of cascade reservoirs on the longitudinal variability of fine sediment characteristics a case study of the lancang and nu rivers abstract the lancang river and nu river are parallel rivers located in southwestern china while the lancang river has several dams the nu river has none to explore the impacts of cascade reservoirs on the grain size and mineral composition of fine sediments 9 fine sediment samples from the nu river and 29 samples from the lancang river were collected and tested the results indicate that in the free flowing reaches ffrs of the lancang river the sediment grain sizes fined systematically and the mineral compositions varied slightly and these results were similar to those of the nu river without dams along the whole mainstream in the regulated reaches rrs of the lancang river the sediment grain sizes and mineral compositions indicated discontinuous fluctuations due to the dams trapping in the reservoir head areas the sediments were composed of fine silt or medium silt with 50 clay minerals but the sediments immediately downstream of the dams were fine sand or medium sand with 71 quartz and feldspar the relative abundance of each mineral was altered differently by passing through the rrs quartz chlorite and illite increased by 13 33 and 20 respectively while feldspar and calcite decreased by 26 and 80 respectively and dolomite and pyrite were completely trapped in the cascade reservoirs all these common minerals may be replenished except pyrite which originates from mining areas in the lancang subbasin

acoustic measurements of the 1999 basaltic eruption of shishaldin volcano alaska 2 precursor to the subplinian phase the 1999 eruption of shishaldin volcano alaska usa displayed both strombolian and subplinian basaltic activity the subplinian phase was preceded by a signal of low amplitude and constant frequency c2 hz lasting 13 h this bhumming signalq is interpreted as the coalescence of the very shallow part of a foam building up in the conduit which produces large gas bubbles before bursting the acoustic waveform of the hum event is modelled by a helmholtz resonator gas is trapped into a rigid cavity and can only escape through a tiny upper hole producing sound waves at shishaldin the radius of the hole c 5m is close to that of the conduit c6 m the cavity has a length of c60 m and gas presents only a small overpressure between c1 2 10 3 and 4 5 10 3 mpa such an overpressure is obtained by the partial coalescence of a foam formed by bubbles with a diameter from c2 3 mm at the beginning of the episode towards c0 64 mm very close to the end of the phase the intermittency between hum events is explained by the ripening of the foam induced by the h2o diffusion through the liquid films the two extreme values from 600 to 10 s correspond to a bubble diameter from 2 2 to 0 3 mm at the beginning and end of the pre subplinian phase respectively the extremely good agreement between two independent estimates of bubble diameters in the shallow foam reinforces the validity of such an interpretation the total gas volume lost at the surface during the humming events is at most 5 9 10 6 m 3 at the very end of the presubplinian phase there is a single large bubble with an overpressure of c0 42 mpa the large overpressure suggests that it comes from significant depth unlike other bubbles in the pre subplinian phase this deep bubble may be responsible for the entire foam collapse resulting in the subplinian phase d 2004 elsevier b v all rights reserved

gas bubble dynamics model for shallow volcanic tremor at stromboli volcanic tremor at stromboli aeolian islands italy is correlated to small infrasonic transients ripepe et al 1996 which repeat almost rhythmically in time in a range between 0 8 and 1 2 s we demonstrate that infrasonic transients are associated to small gas bubble 0 5 m burstings which produces no transients in the seismic signal tremor ground displacement attenuates with the inverse of the distance from the craters indicating that the source is shallow short term energy release shows that infrasonic and seismic signals are linked to the same dynamical process while at the long term scale it is evident that the two signals are controlled by two distinctive mechanisms we suggest that the possible physical model acts in two steps first gas coalescence and then gas bursting in our model the seismic signal is related to the coalescence of a gas bubble from a layer of small bubbles while the infrasonic signal is linked to the bursting of the bubble when it reaches the magma surface gas bubbles could form by free coalescence in magma or could be forced to coalesce by a structural barrier we calculate that forced coalescence induces in magma a pressure change 104 pa 2 orders of magnitude higher than free coalescence and it explains best the tremor ground displacement 10 5 m recorded at stromboli moreover forced coalescence evidences the role of a structural barrier such as a dike in volcanic tremor source dynamics in this gas dynamic process the delay time of 1 2 s between infrasonic pulses could reflect the gas nucleation interval of basaltic magma thomas et al 1993 manga 1996 we propose that the source function for the shallow volcanic tremor at stromboli could be the viscoelastic reaction of the magma to the pressure decrease induced by gas bubble growth rate under constant depressurization the spectrum of our source function is controlled by the time duration of the pressure pulse which represents the viscoelastic relaxation time of the magma and gas bubble growth rate the predicted asymptotic decay of the frequency contents fits the spectral behavior of the volcanic tremor ground displacement recorded at stromboli we show that the same spectral behavior can be found in ground displacement spectra of volcanic tremor recorded on different volcanoes

generation and propagation of infrasonic airwaves from volcanic explosions abstract analysis of infrasonic pressure waves generated by active volcanoes is essential to the understanding of volcanic explosion dynamics unlike seismic waves propagating in the earth infrasonic airwaves offer a relatively unfiltered representation of source motions at the vent during an eruption time varying acoustic propagation filters caused by changeable atmospheric conditions are minimal for microphones deployed at intermediate distances johnson et al j volcanol geotherm res in press

three distinct regimes of volcanic tremor associated with the eruption of shishaldin volcano alaska 1999 tremor signals associated with the eruption of shishaldin volcano on 19 and 23 april 1999 were the strongest recorded anywhere in the aleutian arc by the alaska volcano observatory avo in its 10 year history reduced displacements dr reached 23 cm2 on 19 april and 43 cm2 on 23 april during the activity dr and spectral data with a frequency resolution of 0 1 hz were computed and put on the world wide web every 10 min these data are analyzed here the general temporal patterns of seismicity of these eruption events were similar but the eruptions and their effects quite different the 19 april event is known to have culminated in a sub plinian phase which ejected ash to an altitude of 16 km despite higher amplitudes and the largest hotspot from satellite data the 23 april event produced little ash reaching only 6 km altitude for several hours prior to the sub plinian phase on 19 april tremor with a peak frequency of 1 3 hz intensified during the sub plinian phase the peak frequency increased to 4 8 hz however in 15 h after the eruption three episodes of stronger tremor occurred with a lower 1 0 hz peak alternating with weaker tremor with a 1 3 hz peak these transitions correspond to dr 8 cm2 although these strong tremor episodes produced higher dr levels than the sub plinian phase data from a pressure sensor show that only strong strombolian explosions occurred the suite of observations suggests three distinct tremor regimes that may correspond to slug flow bubbly flow and sustained strong eruptions or a cyclic change in source parameters e g geometry sound speed or ascent rate this behavior occurred at shishaldin only during the april 1999 sequence and we are not aware of similar behavior at other volcanoes