the uwekahuna ash member of the puna basalt product of violent phreatomagmatic eruptions at kilauea volcano hawaii between 2800 and 210014c years ago abstract kilauea volcano s reputation for relatively gentle effusive eruptions belies a violent geologic past including several large phreatic and phreatomagmatic eruptions that are recorded by holocene pyroclastic deposits which mantle kilauea s summit area and the southeast flank of adjacent mauna loa volcano the most widespread of these deposits whose original distribution can be reconstructed is the uwekahuna ash member of the puna basalt a basaltic surge and fall deposit emplaced during two or more eruptive episodes separated by a few decades to several centuries the first episode occurred between 2770 70 and 2265 50 14 c yr ago it included two major pyroclastic surges each preceded by unusually vigorous lava fountaining from a vent near the volcano s summit before the second eruptive episode 2110 120 14 c yr ago plants had re colonized the rainforest environment northeast of the summit and at least two lava flows from mauna loa had buried parts of the first episode deposits the second episode also began with vigorous lava fountaining followed by widespread lithic ashfall a third major surge and finally a fourth fountaining event before the final pumice deposit could be significantly reworked it was partly buried by picritic basalt flows that are unusual in kilauea s summit area in proximal areas the uwekahuna ash member is more than 1 m thick locally 5 m and includes lithic blocks up to 0 8 m in diameter coarse primarily lithic debris was deposited around the vent by laterally expanding surges fallout deposits accumulated preferentially to the northeast under the influence of high altitude counter tradewinds the area devastated by surges and originally buried by at least 15 cm of the uwekahuna was about 420 km 2 the bulk volume of the deposits was approximately 0 3 km 3 including less than 0 1 km 3 of juvenile material juvenile constituents are olivine tholeiitic basalts similar in major element composition to typical kilauea summit lava flows but variations in both major elements and trace elements suggest that the eruptions tapped more than a single uniform source region we infer that the eruptions which produced the uwekahuna were driven by water interacting with a fluctuating magma column magma withdrawal episodes may have been accompanied by large volume submarine effusive eruptions and by summit collapse the volume extent and character of the uwekahuna deposits underscore the hazards posed by relatively infrequent but potentially devastating explosive eruptions at kilauea as well as at other basaltic volcanoes

time dependent spatial amplitude patterns of harmonic tremor at arenal volcano costa rica seismic wave interferences seismograms recorded at the receivers of a small aperture seismic array usually display very similar waveforms and amplitudes as a consequence of their close proximity during the analysis of the volcanic tremor wave field at arenal volcano costa rica we detected significant differences in the amplitudes of harmonic tremor recorded at the stations of a small aperture 210 m seismic array the amplitude distributions are geometrically complex and characterized by strong gradients they occur just for harmonic tremors any other type of seismic event produces nearly uniform amplitudes across the array this suggests some relation with harmonic frequency content moreover the spatial amplitude patterns change with time some of these observations could be explained by an extreme combination of source path and site effects but they also could be produced by interference of different components of the seismic wave field we use numerical calculations to investigate the amplitude pattern generated by two interfering plane waves and are able to reproduce the main features of the observed amplitude patterns we propose mechanisms that might generate seismic wave fields with multiple components and conclude that interference can explain the complexity and variability of the harmonic tremor wave field at arenal volcano online material wave field animation

velocity and spectral characteristics of the volcanic tremor at etna deduced by a small seismometer array abstract a small l shaped short period seismic array was set up at etna volcano to measure the wave velocity and the incoming azimuth of the volcanic tremor generated by the volcanic activity at etna in non eruptive conditions cross correlation techniques were used to measure the phase shifts between the array stations an approximately n s direction of the propagation direction was estimated coincident with the azimuth of the crater area of etna an apparent velocity of less than 500 m s was estimated suggesting a surface propagation of the tremor waves spectral analysis confirming previous results shows velocity spectra peaked at approximately 2 hz a localized source coincident with the crater area agrees with our data and could generate the volcanic tremor at etna in non eruptive stages

shallow velocity structure of stromboli volcano italy derived from small aperture array measurements of strombolian tremor abstract the properties of the tremor wave field at stromboli are analyzed using data from small aperture arrays of short period seismometers deployed on the north flank of the volcano the seismometers are configured in two semi circular arrays with radii of 60 and 150 m and a linear array with length of 600 m the data are analyzed using a spatiotemporal correlation technique specifically designed for the study of the stationary stochastic wave field of rayleigh and love waves generated by volcanic activity and by scattering sources distributed within the island the correlation coefficients derived as a function of frequency for the three components of motion clearly define the dispersion characteristics for both rayleigh and love waves love and rayleigh waves contribute 70 and 30 respectively of the surface wave power the phase velocities of rayleigh waves range from 1000 m sec at 2 hz to 350 m sec at 9 hz and those for love waves range from 800 to 400 m sec over the same frequency band these velocities are similar to those measured near puu oo on the east rift of kilauea volcano hawaii although the dispersion characteristics of rayleigh waves at stromboli show a stronger dependence on frequency such low velocities are consistent with values expected for densely cracked solidified basalt the dispersion curves are inverted for a velocity model beneath the arrays assuming those dispersions represent the fundamental modes of rayleigh and love waves

seismic source mechanism of degassing bursts at kilauea volcano hawaii results from waveform inversion in the 10 50 s band the current march 2008 to february 2009 summit eruptive activity at kilauea volcano is characterized by explosive degassing bursts accompanied by very long period vlp seismic signals we model the source mechanisms of vlp signals in the 10 50 s band using data recorded for 15 bursts with a 10 station broadband network deployed in the summit caldera to determine the source centroid location and source mechanism we minimize the residual error between data and synthetics calculated by the finite difference method for a point source embedded in a homogeneous medium that takes topography into account the vlp signals associated with the bursts originate in a source region 1 km below the eastern perimeter of halemaumau pit crater the observed waveforms are well explained by the combination of a volumetric component and a vertical single force component for the volumetric component several source geometries are obtained which equally explain the observed waveforms these geometries include 1 a pipe dipping 64 to the northeast 2 two intersecting cracks including an east striking crack dike dipping 80 to the north intersecting a north striking crack another dike dipping 65 to the east 3 a pipe dipping 58 to the northeast intersecting a crack dipping 48 to the west southwest and 4 a pipe dipping 57 to the northeast intersecting a pipe dipping 58 to the west southwest using the dual crack model as reference the largest volume change obtained among the 15 bursts is 24 400 m3 and the maximum amplitude peak to peak of the force is 20 gn each burst is marked by a similar sequence of deflation and inflation trailed by decaying oscillations of the volumetric source the vertical force is initially upward synchronous with source deflation then downward synchronous with source reinflation followed by oscillations with polarity opposite to the volumetric oscillations this combination of force and volume change is attributed to pressure and momentum changes induced during a fluid dynamic source mechanism involving the ascent expansion and burst of a large slug of gas within the upper 150 m of the magma conduit as the slug expands upon approach to the surface and more liquid becomes wall supported by viscous shear forces the pressure below the slug decreases inducing conduit deflation and an upward force on the earth the final rapid slug expansion and burst stimulate vlp and lp oscillations of the conduit system which slowly decay due to viscous dissipation and elastic radiation consideration of the fluid dynamic arguments leads us to prefer the dual crack vlp source model as it is the only candidate model capable of producing plausible values of length scales and pressure changes the magnitudes of the vertical forces observed in the 15 bursts appear consistent with slug masses of 104 to 106 kg

spatial extent of a hydrothermal system at kilauea volcano hawaii determined from array analyses of shallow long period seismicity 2 results we present a probabilistic method to locate the source of seismic events using seismic antennas the method is based on a comparison of the event azimuths and slownesses derived from frequency slowness analyses of array data with a slowness vector model several slowness vector models are considered including both homogeneous and horizontally layered half spaces and also a more complex medium representing the actual topography and three dimensional velocity structure of the region under study in this latter model the slowness vector is obtained from frequency slowness analyses of synthetic signals these signals are generated using the finite difference method and include the effects of topography and velocity structure to reproduce as closely as possible the behavior of the observed wave fields a comparison of these results with those obtained with a homogeneous half space demonstrates the importance of structural and topographic effects which if ignored lead to a bias in the source location we use synthetic seismograms to test the accuracy and stability of the method and to investigate the effect of our choice of probability distributions we conclude that this location method can provide the source position of shallow events within a complex volcanic structure such as kilauea volcano with an error of 200 m

iron isotope geochemistry of biogenic magnetite bearing sediments from the bay of vidy lake geneva article dissimilatory microbial iron oxide reduction dir has been hypothesized to be an important respiratory pathway on early earth potentially generating significant quantities of fe ii that have been preserved in proterozoic and archean sedimentary rocks in particular dir has been implicated in the formation of magnetite inprecambrian marine sediments todate however onlyone modern sedimentary environmentexists wherein situ magnetite formation has been linked to dir the bay of vidy in lake geneva switzerland previous work at this locality has characterized a magnetic susceptibility anomaly that reflects the presence of fine grained magnetite produced viamicrobial reduction of amorphous fe iii oxidesthat enterthebayof vidy from a nearby sewage treatment plant in this study we report on the fe isotope composition of aqueous and solid phase fe in the bay of vidy sediments extensive fe iii reduction has occurred resulting in the conversion of nearly all reactive non silicate fe iii to a variety of fe ii bearing phases with mixed fe valence magnetite being a minor but easily detectable component 0 5 8 wt very little fe isotope variation was observed in any solid phase fe components including magnetite although significant fractionation was observed between aqueous and solid phase fe ii because fe mass balance was dominated by the solid phase little net change in 56 fe

mechanism of explosive eruptions of kilauea volcano hawaii a small explosive eruption of kilauea volcano hawaii occurred in may 1924 the eruption was preceded by rapid draining of a lava lake and transfer of a large volume of magma from the summit reservoir to the east rift zone this lowered the magma column which reduced hydrostatic pressure beneath halemaumau and allowed groundwater to flow rapidly into areas of hot rock producing a phreatic eruption a comparison with other events at kilauea shows that the transfer of a large volume of magma out of the summit reservoir is not sufficient to produce a phreatic eruption for example the volume transferred at the beginning of explosive activity in may 1924 was less than the volumes transferred in march 1955 and january february 1960 when no explosive activity occurred likewise draining of a lava lake and deepening of the floor of halemaumau which occurred in may 1922 and august 1923 were not sufficient to produce explosive activity a phreatic eruption of kilauea requires both the transfer of a large volume of magma from the summit reservoir and the rapid removal of magma from near the surface where the surrounding rocks have been heated to a sufficient temperature to produce steam explosions when suddenly contacted by groundwater

location of long period events below kilauea volcano using seismic amplitudes and accurate relative relocation 1 we present methods for improving the location of long period lp events deep and shallow recorded below kilauea volcano by the permanent seismic network lp events might be of particular interest to understanding eruptive processes as their source mechanism is assumed to directly involve fluid transport however it is usually difficult or impossible to locate their source using traditional arrival time methods because of emergent wave arrivals at kilauea similar lp waveform signatures suggest the existence of lp multiplets the waveform similarity suggests spatially close sources while catalog solutions using arrival time estimates are widely scattered beneath kilauea s summit caldera in order to improve estimates of absolute lp location we use the distribution of seismic amplitudes corrected for station site effects the decay of the amplitude as a function of hypocentral distance is used for inferring lp location in a second stage we use the similarity of the events to calculate their relative positions the analysis of the entire lp seismicity recorded between january 1997 and december 1999 suggests that a very large part of the lp event population both deep and shallow is generated by a small number of compact sources deep events are systematically composed of a weak high frequency onset followed by a low frequency wave train aligning the low frequency wave trains does not lead to aligning the onsets indicating the two parts of the signal are dissociated this observation favors an interpretation in terms of triggering and resonance of a magmatic conduit instead of defining fault planes the precise relocation of similar lp events based on the alignment of the high energy low frequency wave trains defines limited size volumes

temporal and spatial characteristics of volcanic tremor wave fields abstract detailed volcanic tremor studies were carried out with portable seismic stations installed both in the near and distant zones at klyuchevskoy volcano in 1984 1986 and 1987 it has been determined that the coherence of wave field decreases with distance to 0 5 km and does not change at greater distances based on coherence and phase spectra the frequency dependence of the phase velocities has been determined the dispersion of phase velocity points to the surface character of the wave field the behaviour of volcanic tremor characteristics through time is stationary which implies stability of the parameters of the volcanic tremor source a detailed spectral study has demonstrated that when the intensity of volcanic tremor varies by a few orders of magnitude parts of the low frequency and high frequency spectra change which is likely due to the variations in source dimensions the detailed spectral characteristics of tremor may be used for both determining the source parameters and studying the structure of the medium

three dimensional p and s wave velocity structure of redoubt volcano alaska the three dimensional p and s wave structure of redoubt volcano alaska and the underlying crust to depths of 7 8 km is determined from 6219 p wave and 4008 s wave first arrival times recorded by a 30 station seismograph network deployed on and around the volcano first arrival times are calculated using a finite difference technique which allows for flexible parameterization of the slowness model and easy inclusion of topography and source receiver geometry the three dimensional p wave velocity structure and hypocenters are determined simultaneously while the three dimensional s wave velocity model is determined using the relocated seismicity and an initial s wave velocity model derived from the p wave velocity model assuming an average vp vs ratio of 1 78 convergence is steady with approximately 73 and 52 reduction in p and s wave arrival time rms respectively after 10 iterations the most prominent feature observed in the three dimensional velocity models derived for both p and s waves is a relative low velocity near vertical pipelike structure approximately 1 km in diameter that extends from 1 to 6 km beneath sea level this feature aligns axially with the bulk of seismicity and is interpreted as a highly fractured and altered zone encompassing a magma conduit the velocity structure beneath the north flank of the volcano between depths of 1 and 6 km is characterized by large lateral velocity variations high velocities within this region are interpreted as remnant dikes and sills and low velocities as regions along which magma migrates no large low velocity body suggestive of a magma chamber is resolved in the upper 7 8 km of the crust

three dimensional p wave velocity structure of mt etna italy the three dimensional p wave velocity structure of mt etna is determined to depths of 15 km by tomographic inversion of first arrival times from local earthquakes recorded by a network of 29 permanent and temporary seismographs results show a near vertical low velocity zone that extends from beneath the central craters to a depth of 10 km this low velocity region is coincident with a band of steeply dipping seismicity suggesting a magmatic conduit that feeds the summit eruptions the most prominent structure is an approximately 8 km diameter high velocity body located between 2 and 12 km depth below the southeast flank of the volcano this high velocity body is interpreted as a remnant mafic intrusion that is an important structural feature influencing both volcanism and east flank slope stability and faulting

modeling the three dimensional structure of macroscopic magma transport systems application to kilauea volcano hawaii we report the results of modeling the three dimensional internal structure of kilauea s magmatic passageways the approach uses a clear plexiglass model containing equally spaced levels upon which well located seismic hypocenters are plotted application of constraining geologic and geophysical criteria to this distributed volume of earthquakes permits the interpretation of seismic structures produced by fracturing in response to locally high fluid pressures four magma transport and storage structures have been identified within and beneath kilauea 1 primary conduit the conduit transporting magma into kilauea s summit storage reservoir rises from the model base 14 6 km to the 6 5 km depth level it is a zone of intense fracturing and inferred intrusion whose horizontal sections are elliptical in planform over its height the average major axis of the component horizontal sections is 3 3 km with an average minor axis of 1 7 km this yields an aspect ratio of 0 52 at the 14 6 km level the strike of the major axis is n67 e during passage from the upper mantle through the oceanic crust this axis rotates in a right handed sense until the strike is n41 w at the 6 5 km level 2 magma chamber complex floor the interval from 6 5 to 5 7 km immediately over the primary conduit is aseismic this suggests differentially high fluid to rock ratios and relatively weak pathways for further vertical transport into higher levels of the storage complex as well as lateral leakage eastward into the mauna ulu staging area for later vertical ascent beneath the upper east rift zone seismicity within the immediately subjacent rocks that form the top of the primary conduit at 6 5 km suggests that this inferred magma rich horizon forms the effective floor of the summit storage complex 3 magma chamber crown intense seismicity over the 1 1 1 9 km depth interval defines an elliptical region in plan view the top of this region has a broad apex with an average major axis of 1 5 km and an average minor axis of 1 1 km producing an aspect ratio of 0 73 this region coincides with the epicentral position of summit vertical displacement maxima during inflationary deformation and is interpreted as a region of intense diking crowning the summit storage complex 4 upper east rift zone pipe from a depth of 5 7 km a cylindrical region of seismicity rises beneath the upper east rift zone to the 1 9 km level where it merges with the subhorizontal upper east rift zone duct this pipe like zone has a mean diameter of 1 4 km and a vertical axis that pierces the surface of the volcano at the intersection of the koae fracture system and the upper east rift zone we suggest this pipe forms a principal connection between the base of the summit storage reservoir and near surface storage compartments and ducts in the upper east rift zone in this role the pipe is suggested to have contributed heavily to mauna ulu s magma supply during 1969 1974 the model of magma transport and storage thus constructed provides a context in which other types of data e g ground deformation and electrical conductivity may be interpreted as well as a means of refining petrological inferences based on lava chemistry

a double seismic antenna experiment at teide volcano existence of local seismicity and lack of evidences of volcanic tremor data analyzed in the present work correspond to a 40 days field experiment carried out in teide volcano canary islands spain with two short period small aperture dense seismic antennas in 1994 the objective of this experiment was to detect analyze and locate the local seismicity we analyzed also the background seismic noise to investigate the possible presence of volcanic tremor from a set of 76 events we selected 21 of them in base of their good signal to noise ratio and their possibility to locate their seismic source by using the seismic antennas a visual classification based on the s p time and seismogram shape has permitted to establish three groups of events local seismicity s p time between 3 and 5 s very local earthquakes s p time smaller than 3 s and artificial explosions these earthquakes have been located by applying the zero lag crosscorrelation technique and the inverse ray tracing procedure those earthquakes that were recorded simultaneously by both seismic antennas were also located by intersecting both back azimuths the analysis of the seismicity has revealed that the amount of seismicity in teide volcano is moderate this seismicity could be distributed in three main areas inside the caldera edifice below the teide pico viejo complex in the eastern border of the caldera edifice and offshore of the island at present this activity is the only indicator of the volcano dynamics the analysis of the back ground seismic noise has revealed that at frequencies lower than 2 hz the oceanic load signal is predominant over other signals even over local earthquakes with a magnitude of 2 0 due to this although if in the teide area were present a weak volcanic tremor or other volcanic signals with predominant peaks below 2 hz to observe them would be a very difficult task q 2000 elsevier science b v all rights reserved

shallow structure of mt vesuvius volcano italy from seismic array analysis data from a portable dense seismic array deployed on mt vesuvius in may 1994 during a 2d seismic tomography experiment are analyzed in the present paper the array consisted of two groups of short period geophones 4 5 hz natural frequency formed by 16 and 25 vertical components plus two horizontal components distributed along an arc like shape along the summit crater stacks of later arrivals interpreted as reflected phases provide a significant constraint to a boundary layer located in the depth range 1 5 2 2 km beneath the summit crater with average velocity v 1 8 2 2 km s interpreted as the top of the limestone basement the correlation methods applied on microtremor records allowed to infer the shallow velocity structure up to 400 m beneath the crater rim

observations of loma prieta aftershocks from a dense array in sunnyvale california a small aperture 300 m four station array was deployed in sunnyvale california for 5 days to record aftershocks of the loma prieta earthquake of october 1989 the purpose of the array was to study the seismic response of the alluvium filled santa clara valley and the role of surface waves in the seismic shaking of sedimentary basins strong motion records of the loma prieta mainshock indicate that surface waves produced the peak velocities and displacements at some sites in the santa clara valley we use the recordings from the dense array to determine the apparent velocity and azimuth of propagation for various arrivals in the seismograms of four aftershocks with magnitudes between 3 6 and 4 4 apparent velocities are generally observed to decrease with increasing time after the s wave in the seismograms phases arriving less than about 8 sec after the s wave have apparent velocities comparable to the s wave and appear to be body waves multiply reflected under the receiver site or reflected by crustal interfaces for times 10 to 30 sec after the direct s wave we observe long period 1 to 6 sec arrivals with apparent velocities decreasing from 2 5 to 0 8 km sec we interpret these arrivals to be surface waves and conclude that these surface waves produce the long duration of shaking observed on the aftershock records much of the energy in the 40 sec after the s wave is coming approximately from the direction of the source although some arrivals have backazimuths as much as 60 different from the backazimuths to the epicenters two of the aftershocks show arrivals coming from 30 to 40 more easterly than the epicenters this energy may have been scattered from outcrops along the southeastern edge of the basin in contrast the deepest aftershock studied d 17 km displays later arrivals with backazimuths 30 to 40 more westerly than the epicenter a distinct arrival for one of the aftershocks propagates from the southwest possibly scattered from the western edge of the basin synthetic seismograms derived from a plane layered crustal model do not produce the long period love waves observed in the waveforms of the ml 4 4 aftershock these love waves may be generated by the conversion of incident s waves or rayleigh waves near the edge of the basin

array analysis using circular wave front geometry an application to locate the nearby seismo volcanic source summary the zero lag cross correlation technique used for array analysis in the hypothesis of plane waves has been modified to allow the wave front to be circular synthetic tests have been performed to check the capability of the method which returns the input test data when the source array distances are not greater than two or three times the array aperture for this distance range the method furnishes an estimate of the apparent velocity and the epicentral coordinates of the source for more distant sources the method becomes equivalent to that based on the planar wave approximation which gives an estimate of the backazimuth to the source and the apparent velocity the method has been applied to seismic data recorded at the active volcano located at deception island antarctica 35 volcanic long period events occurring in a small swarm were selected results show that the epicentres are close to the array between 0 4 and 2 km and aligned in a sw direction in agreement with one of the main directions of the fracture system of deception volcano

characteristics of seismic waves composing hawaiian volcanic tremor and gas piston events observed by a near source array a correlation method specifically designed for describing the characteristics of a complex wave field is applied to volcanic tremor and gas piston events recorded by a semicircular array of geos instruments set at the foot of the puu oo crater on the east rift of kilauea volcano hawaii the spatial patterns of correlation coefficients obtained as functions of frequency for the three components of motion over the entire array are similar for gas piston events and tremor and clearly depict dispersive waves propagating across the array from the direction of puu oo the wave fields are composed of comparable amounts of rayleigh and love waves propagating with similar and extremely slow phase velocities ranging from 700 m s at 2 hz to 300 m s at 8 hz the highly cracked solidified lava flow on which the array was deployed and subjacent structure of alternating lava and ash layers formed during repeated eruptions of puu oo since 1983 appear to be responsible for the low velocities observed the results from puu oo stand in sharp contrast to those obtained in an experiment conducted in 1976 on the partially solidified lava lake of kilauea iki rayleigh waves were not observed in kilauea iki but well developed trains of love waves were seen to propagate there with velocities twice as high as those observed near puu oo these differences in the propagation characteristics of surface waves at the two sites may be attributed to the presence of a soft horizontal layer of molten rock in kilauea iki which may have lowered the phase velocity of rayleigh waves more drastically than that of love waves resulting in severe scattering of the rayleigh wave mode on the other hand the thin superficial pahoehoe flow under our array at puu oo may have favored the development of vertical columnar joints more extensively at this location than at kilauea iki which may have reduced the shear moduli controlling the love wave mode the average phase velocities in the frequency band from 2 to 5 hz found at puu oo are roughly similar to those determined for tremor at klyuchevskoy volcano kamchatka but the frequency dependence appears much stronger at puu oo

three dimensional velocity structure of the kilauea caldera hawaii high resolution velocity models 0 5 km resolution of the kilauea caldera region are obtained by the tomographic inversion of both p and s wave arrival times data are from the permanent hawaiian volcano observatory hvo seismic network a broadband seismic network and a temporary array of stations centered on the southern boundary of the caldera a low velocity p wave anomaly is imaged centered on the southeastern edge of the caldera with a velocity contrast of about 10 and a volume of 27 km the vp vs model mimics the spatial extent of the p wave anomaly but is partitioned into two discrete anomalous volumes centered on the southern boundary of the caldera and on the upper east rift of the volcano the corresponding poisson s ratio in these zones is high 0 25 0 32 which is consistent with a densely cracked hot volume which may contain partial melt the large scale features of the models are consistent with results obtained from an earlier larger scale 2 km resolution tomographic image of kilauea volcano based on hvo network data

seismovolcanic signals at deception island volcano antarctica wave field analysis and source modeling the seismovolcanic signals associated with the volcanic activity of deception island antarctica recorded during three antarctic summers 1994 1995 1995 1996 and 1996 1997 are analyzed using a dense small aperture 500 m seismic array the visual and spectral classification of the seismic events shows the existence of long period and hybrid isolated seismic events and of low frequency quasi monochromatic and spasmodic continuous tremors all spectra have the highest amplitudes in the frequency band between 1 and 4 hz while hybrids and spasmodic tremors have also significant amplitudes in the high frequency band 4 10 hz the array analysis indicates that almost all the well correlated low frequency signals share similar array parameters slowness and back azimuth and have the same source area close to the array site the polarization analysis shows that phases at high frequency are mostly composed of p waves and those phases dominated by low frequencies can be interpreted as surface waves no clear shear waves are evidenced from the energy evaluation we have found that the reduced displacement values for surface and body waves are confined in a narrow interval volcano tectonic seismicity is located close to the array at a depth shallower than 1 km the wave field properties of the seismovolcanic signals allow us to assume a unique source model a shallow resonating fluid filled crack system at a depth of some hundreds of meters all of the seismic activity is interpreted as the response of a reasonably stable stationary geothermal process the differences observed in the back azimuth between low and high frequencies are a near field effect a few episodes of the degassification process in an open conduit were observed and modeled with a simple organ pipe

deterministic frequency wavenumber methods and direct measurements of rupture propagation during earthquakes using a dense array data analysis data recorded by the smart 1 array in taiwan are used to make the first direct measurements of two dimensional earthquake rupture propagation using frequency wavenumber techniques and ray theory we obtain estimates of the spatial extent duration and directions of rupture propagation during two earthquakes located offshore taiwan our estimates include uncertainties due to errors in fault plane location and orientation velocity structure and slowness ray parameter measured at the array we find that the january 29 1981 ml 6 3 taiwan earthquake ruptured unilaterally updip and toward the west on a 60 dipping 109 striking reverse fault using p waves we find its fault length and duration to be 25 18 km and 7 4 3 4 s respectively s waves indicate a fault length and duration of 27 15 km and 9 4 3 6 s respectively we suggest that the november 14 1986 ml 7 0 hualien taiwan earthquake was triggered by a shallow h 14 km foreshock ml 5 4 which occurred on a subparallel or splay fault approximately 5 s before the mainshock based on p waves rupture propagation during the foreshock is unilateral towards the northeast and slightly downdip on a 58 dipping 38 striking reverse fault its spatial extent and duration are approximately 23 km and 5 4 s the mainshock initiated at greater depth h 34 km than the foreshock with its primary direction of rupture unilateral toward the northeast and slightly updip on a 57 dipping 43 striking reverse fault using p waves we found a spatial extent and duration of 68 15 km and 12 4 2 2 s s waves indicate a similar extent duration and direction of propagation for the mainshock however we found it difficult to distinguish the foreshock and mainshock 5 waves we show by example that these results represent an improvement over what can be obtained using long period t 5 s teleseismic methods comparison of our results with those obtained using teleseismic data suggests either a short duration slip time function or that most slip occurs near the rupture front our results indicate that much of the complexity in the seismograms recorded during these events by the smart 1 array is due to spatial and temporal variations of the source this suggests that at least in some cases predictions of both the phase and amplitude of strong ground motion as a function of time would require a priori knowledge of source properties finally we recommend siting small 10 station dense arrays on land above offshore subduction zones which are considered seismic gaps because they can provide more detailed measurements of rupture propagation than is presently possible

fluctuations in surface height of active lava lakes during 1972 1974 mauna ulu eruption kilauea volcano hawaii an active lava lake is one that forms directly on top of a magma column as opposed to an inactive one formed by passive ponding of lava within a preexisting crater two active lava lakes operating within summit craters of newly constructed volcanic shields at mauna ulu and alae were the longest lived of any observed in a hawaiian rift zone to date small scale fluctuations a few to 20 m in the height of lake surfaces are typically associated with degassing gas piston activity larger fluctuations as great as 200 m reflect major changes in eruptive activity measurements of surface height fluctuations of mauna ulu lake show gross positive correlation with variation in summit tilt which in turn reflects the inflation deflation state of kilauea s summit magma reservoir a more tenuous conclusion is that a sympathetic variation possibly exists between surface heights of the mauna ulu and alae lakes collectively these observations 1 require the existence of a fully engorged connected open and hydraulically efficient magma conduit system linking the summit reservoir with the lava lakes 2 provide additional evidence for the postulated regime of quasi steady state magma transfer that characterized much of the eruption 3 imply that fluctuations in lake surface reflect variation in magma supply rate to the reservoir and or outflow rate from the mauna ulu alae vents and 4 suggest undetected episodes of magma intrusion or storage in kilauea s east rift zone during 1973

imaging the crustal magma sources beneath mauna loa and kilauea volcanoes hawaii three dimensional seismic p wave traveltime tomography is used to image the magma sources beneath mauna loa and kilauea volcanoes hawaii high velocity bodies 6 4 km s in the upper 9 km of the crust beneath the summits and rift zones of the volcanoes correlate with zones of high magnetic intensities and are interpreted as solidified gabbro ultramafic cumulates from which the surface volcanism is derived the proximity of these high velocity features to the rift zones is consistent with a ridge spreading model of the volcanic flank southeast of the hilina fault zone along the south flank of kilauea low velocity material

shallow velocity models at the kilauea volcano hawaii determined from array analyses of tremor wavefields the properties of the surface wavefield at kilauea volcano are analysed using data from small aperture arrays of short period seismometers deployed in and around the kilauea caldera tremor recordings were obtained during two japan us cooperative experiments conducted in 1996 and 1997 the seismometers were deployed in three semi circular arrays with apertures of 300 300 and 400 m and a linear array with length of 1680 m data are analysed using a spatio temporal correlation technique well suited for the study of the stationary stochastic wavefields of rayleigh and love waves associated with volcanic activity and scattering sources distributed in and around the summit caldera spatial autocorrelation coefficients are obtained as a function of frequency and are inverted for the dispersion characteristics of rayleigh and love waves using a grid search that seeks phase velocities for which the l 2 norm between data and forward modelling operators is minimized within the caldera the phase velocities of rayleigh waves range from 1400 to 1800 m s 1 at 1 hz down to 300 400 m s 1 at 10 hz and the phase velocities of love waves range from 2600 to 400 m s 1 within the same frequency band outside the caldera rayleigh wave velocities range from 1800 to 1600 m s 1 at 1 hz down to 260 360 m s 1 at 10 hz and love wave velocities range from 600 to 150 m s 1 within the same frequency band the dispersion curves are inverted for velocity structure beneath each array assuming these dispersions represent the fundamental modes of rayleigh and love waves the velocity structures observed at different array sites are consistent with results from a recent 3 d traveltime tomography of the caldera region and point to a marked velocity discontinuity associated with the southern caldera boundary

multiscale matched field processing for noise source localization in exploration geophysics abstractin the study we investigate the use of ambient noise data to locate microseismic sources at the exploration scale we develop a multiscale matched field processing mfp approach to localize seismic sources at frequencies below 10 hz an application to an actual data set acquired over a hydrocarbon field is presented to determine the reliability of the mfp procedure the data used were continuously recorded over five days at a total of 397 stations on a 1 km per side square seismic network the mfp results show 1 a dominant and stable surface source associated with human activities road and exploration platforms around the reservoir and 2 weaker sources at depth below the seismic network that are related to the injection extraction process

wavefield properties of a shallow long period event and tremor at kilauea volcano hawaii the wavefields of tremor and a long period lp event associated with the ongoing eruptive activity at kilauea volcano hawaii are investigated using a combination of dense small aperture 300 m and sparse large aperture 5 km arrays deployed in the vicinity of the summit caldera measurements of azimuth and slowness for tremor recorded on the small aperture array indicate a bimodal nature of the observed wavefield at frequencies below 2 hz the wavefield is dominated by body waves impinging the array with steep incidence these arrivals are attributed to the oceanic microseismic noise in the 2 6 hz band the wavefield is dominated by waves propagating from sources located at shallow depths 1 km beneath the eastern edge of the halemaumau pit crater the hypocenter of the lp event determined from frequency slowness analyses combined with phase picks appears to be located close to the source of tremor but at a shallower depth 0 1 km the wavefields of tremor and lp event are characterized by a complex composition of body and surface waves whose propagation and polarization properties are strongly affected by topographic and structural features in the summit caldera region analyses of the directional properties of the wavefield in the 2 6 hz band point to the directions of main scattering sources which are consistent with pronounced velocity contrasts imaged in a high resolution three dimensional velocity model of the caldera region the frequency and q of the dominant peak observed in the spectra of the lp event may be explained as the dominant oscillation mode of a crack with scale length 20 100 m and aperture of a few centimeters filled with bubbly water the mechanism driving the shallow tremor appears to he consistent with a sustained excitation originating in the oscillations of a bubble cloud resulting from vesiculation and degassing in the magma

seismic source dynamics of gas piston activity at kilauea volcano hawai i since 2008 eruptive activity at the summit of kilauea volcano hawai i has been confined to the new overlook pit crater within the halema uma u crater among the broad range of magmatic processes observed in the new pit are recurring episodes of gas pistoning the gas piston activity is accompanied by seismic signals that are recorded by a broadband network deployed in the summit caldera we use raw data recorded with this network to model the source mechanism of representative gas piston events in a sequence that occurred on 20 25 august 2011 during a gentle inflation of the kilauea summit to determine the source centroid location and source mechanism we minimize the residual error between data and synthetics calculated by the finite difference method for a point source embedded in a homogeneous medium that takes topography into account we apply a new waveform inversion method that accounts for the contributions from both translation and tilt in horizontal seismograms through the use of green s functions representing the seismometer response to translation and tilt ground motions this method enables a robust description of the source mechanism over the period range 1 10 000 s most of the seismic wavefield produced by gas pistoning originates in a source region 1 km below the eastern perimeter of the halema uma u pit crater the observed waveforms are well explained by a simple volumetric source with geometry composed of two intersecting cracks featuring an east striking crack dike dipping 80 to the north intersecting a north striking crack another dike dipping 65 to the east each gas piston event is marked by a similar rapid inflation lasting a few minutes trailed by a slower deflation ramp extending up to 15 min attributed to the efficient coupling at the source centroid location of the pressure and momentum changes accompanying the growth and collapse of a layer of foam at the top of the lava column assuming a simple lumped parameter representation of the shallow magmatic system the observed pressure and volume variations can be modeled with the following attributes foam thickness 10 50 m foam cell diameter 0 04 0 10 m and gas injection velocity 0 01 0 06 m s 1 gas piston activity occurs in a narrow pipe with diameter of 6 m connecting the halema uma u pit crater to the subjacent dike system the height of the magma column is estimated at 104 m at the start of the sequence based on the period of very long period vlp oscillations accompanying the onset of the gas piston signal based on the change in the period of vlp oscillations and tilt evidence the height of the magma column is inferred to have risen by up to 23 m by the end of the 5 day long sequence a penny shaped crack model of the dike geometry yields effective diameters of 1 2 2 9 km for the east dike and 0 7 km for the north dike the shallower north dike segment is embedded in a relatively weak medium compatible with expected mechanical properties in the hydrothermal environment of this dike

array tracking of the volcanic tremor source at deception island antarctica we have found experimental evidence which shows that the volcanic tremor recorded at deception island south shetland islands antarctica is a superposition in time of overlapping hybrid events we studied data from a small aperture seismic array data analysis for tremor and hybrids included 1 spectral analysis 2 apparent slowness and back azimuth determination by using the zero lag cross correlation method and 3 polarization analysis both types of events share these common features a two dominant spectral bands at frequencies 1 3 hz the most energetic and 4 8 hz b several coherent phases with the same back azimuth to the source and apparent slowness along the whole signal c in the high frequency band the apparent slowness is very low around 0 17 s km indicating the propagation of body waves d in the low frequency band the apparent slowness is high around 1 6 s km consistent with the presence of surface waves and e clear p wave onset followed by a complex pattern of rayleigh waves therefore both types of events are strongly related because they share the same source region the same wave propagation properties and the same wave composition moreover several arrivals that resemble a single hybrid event have been found along the tremor signals due to these reasons we hypothesize that volcanic tremor of deception island is a superposition of hybrid type events the source of both types could be the interaction between thaw water and hot materials in a shallow aquifer

top down inflation and deflation at the summit of kīlauea volcano hawai i observed with insar 1 we use interferometric synthetic aperture radar insar to study deformation of the summit caldera at k lauea volcano during 2000 2008 which spanned both an east rift zone eruptive event in 2007 and the start of the ongoing summit eruption in 2008 the data set consists of small baseline subset sbas time series generated from 270 acquisitions on three separate beam modes from the radarsat 1 satellite we identify 12 time periods with distinct patterns of displacement that we attribute until late 2003 to secular tectonic driven deformation and from 2004 to 2008 to four different sources in the summit area we model the shallow magmatic system as a spherical reservoir at 1 9 0 1 km depth below the surface to the northeast of halemaumau source 1 and three vertically stacked sills at greater depths in the southern caldera area source 2 at the southern edge of the caldera at 2 9 0 2 km depth source 3 to the south southeast of the caldera at 3 4 0 5 km depth and source 4 south of the caldera at 3 6 0 4 km depth the sequence for filling of and withdrawal from these reservoirs reveals a top down process with sequences of both inflation and deflation initiating in the shallowest source inflation of source 3 is coincident with elevated seismic activity in the upper east rift zone in february 2006 and may 2007 source 4 is elongated toward the southwest rift zone and also shows elevated seismicity that extends toward the southwest rift zone

precursory swarms of long period events at redoubt volcano 1989 1990 alaska their origin and use as a forecasting tool abstract during the eruption of redoubt volcano from december 1989 through april 1990 the alaska volcano observatory issued advance warnings of several tephra eruptions based on changes in seismic activity related to the occurrence of precursory swarms of long period lp seismic events dominant period of about 0 5 s the initial eruption on december 14 occurred after 23 years of quiescence and was heralded by a 23 hour swarm of lp events that ended abruptly with the eruption after a series of vent clearing explosions over the next few days dome growth began on december 21 another swarm with lp events similar to those of the first began on the 26th and ended in a major tephra eruption on january 2 eruptions continued over the next two weeks and then ceased until february 15 when a large eruption initiated a long phase of repetitive dome building and dome destroying episodes that continued into april warnings were issued before the major events on december 14 and january 2 but as the eruptive sequence continued after january 2 the energy of the swarms decreased and forecasting became more difficult a significant but less intense swarm preceded the february 15 eruption which was not forecast this eruption destroyed the only seismograph on the volcanic edifice and stymied forecasting until march 4 when the first of three new stations was installed within 3 km of the active vent from march 4 to the end of the sequence on april 21 there were eight eruptions six of which were preceded by detectable swarms of lp events although weak these swarms provided the basis for warnings issued before the eruptions on march 23 and april 6 the initial swarm on december 13 had the following features 1 short duration 23 hours 2 a rapidly accelerating rate of seismic energy release over the first 18 hours of the swarm followed by a decline of activity during the 5 hours preceding the eruption 3 a magnitude range from 0 4 to 1 6 4 nearly identical lp signatures with a dominant period near 0 5 s 5 dilatational first motions everywhere and 6 a stationary source location at a depth of 1 4 km beneath the crater this occurrence of long period events suggests a model involving the interaction of magma with groundwater in which magmatic gases steam and water drive a fixed conduit at a stationary point throughout the swarm the initiation of that sequence of events is analogous to the failure of a pressure relief valve connecting a lower supercharged magma dominated reservoir to a shallow hydrothermal system a three dimensional model of a vibrating fluid filled crack recently developed by chouet is found to be compatible with the seismic data and yields the following parameters for the lp source crack length 280 380 m crack width 140 190 m crack thickness 0 05 0 20 m crack stiffness 100 200 sound speed of fluid 0 8 1 3 km s compressional wave speed of rock 5 1 km s density ratio of fluid to rock 0 4 and ratio of bulk modulus of fluid to rigidity of rock 0 03 0 07 the fluid filled crack is excited intermittently by an impulsive pressure drop that varies in magnitude within the range of 0 4 to 40 bar such disturbance appears to be consistent with a triggering mechanism associated with choked flow conditions in the crack

seismological analysis of conduit dynamics in fragmentation experiments we simulate gas burst and volcanic explosions under controlled laboratory conditions inducing fragmentation of volcanic rocks by rapid depressurization a series of experiments were performed in a shock tube apparatus at room temperature and a pressure range of 4 to 20 mpa using argon ar gas and particles or pumice samples of different porosities the instrumentation of this system with high precision piezoelectric sensors enabled us to capture elastic waves and to recognize their characteristic signatures by relating these signals to physical processes in the wave field we have been able to characterize the conduit mechanism and the source dynamics we compare and discuss conspicuous features of the waveforms and frequency spectra of these experimental signals with those of volcanic origin despite the fact that these signals are different in amplitude resulting from different scale conditions our observations indicate that the physical processes that occur during simulated explosions and those that occur during volcanic eruptions yield comparable signatures in their respective records the effects of the source receiver configuration and resonance also have significant implications all this suggests that the physical processes e g pressurization and depressurization of a system involve a system response that causes similar distinctive effects independent of the system size reflecting its intrinsic dynamics these similarities imply that powerful constraints on the source mechanisms of volcanic seismicity can emerge from seismic investigations of experimental simulations of volumetric sources such constraints may yield significant advances in the understanding of volcanic conduit dynamics and in the interpretation of seismic unrest at volcanic centers

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

influence of near surface volcanic structure on long period seismic signals and on moment tensor inversions simulated examples from mount etna 1 long period lp seismicity on volcanoes is thought to be associated with moving fluids or resonating fluid filled conduits hence lp moment tensor mt source inversions might have a direct bearing on our understanding of the plumbing system using 3 d full wavefield simulations and 2 d sensitivity kernels in a digital elevation model of mount etna we investigate the influence of near surface volcanic structure on lp signals and on moment tensor inversions contrary to common wisdom in crustal seismology we find that despite their relatively long wavelengths lps are severely distorted by near surface structures including layering and topographic features in particular near surface low velocity structures which are commonly observed on volcanoes play a critical role in controlling the nature of lp signals if not accounted for these path effects leak into the source solution leading to the emergence of erroneous source geometries spurious forces and incorrect source time functions this is particularly problematic if one adopts an unconstrained solution space for the source with many free parameters hence there is a fine balance in the trade off between the velocity model and the source in the absence of high resolution near surface velocity control we demonstrate the importance of employing a priori source information from other fields e g structural geology for shallow lps constraining the number of free parameters in the inversion a probabilistic approach should then be taken as the model with the best fit is not necessarily the true solution

earthquake classification location and error analysis in a volcanic environment implications for the magmatic system of the 1989 1990 eruptions at redoubt volcano alaska abstract determination of the precise locations of seismic events associated with the 1989 1990 eruptions of redoubt volcano posed a number of problems including poorly known crustal velocities a sparse station distribution and an abundance of events with emergent phase onsets in addition the high relief of the volcano could not be incorporated into the hypoellipse earthquake location algorithm this algorithm was modified to allow hypocenters to be located above the elevation of the seismic stations the velocity model was calibrated on the basis of a posteruptive seismic survey in which four chemical explosions were recorded by eight stations of the permanent network supplemented with 20 temporary seismographs deployed on and around the volcanic edifice the model consists of a stack of homogeneous horizontal layers setting the top of the model at the summit allows events to be located anywhere within the volcanic edifice detailed analysis of hypocentral errors shows that the long period lp events constituting the vigorous 23 hour swarm that preceded the initial eruption on december 14 could have originated from a point 1 4 km below the crater floor a similar analysis of lp events in the swarm preceding the major eruption on january 2 shows they also could have originated from a point the location of which is shifted 0 8 km northwest and 0 7 km deeper than the source of the initial swarm we suggest this shift in lp activity reflects a northward jump in the pathway for magmatic gases caused by the sealing of the initial pathway by magma extrusion during the last half of december volcano tectonic vt earthquakes did not occur until after the initial 23 hour long swarm they began slowly just below the lp source and their rate of occurrence increased after the eruption of 01 52 ast on december 15 when they shifted to depths of 6 to 10 km after january 2 the vt activity migrated gradually northward this migration suggests northward propagating withdrawal of magma from a plexus of dikes and or sills located in the 6 to 10 km depth range precise relocations of selected events prior to january 2 clearly resolve a narrow steeply dipping pencil shaped concentration of activity in the depth range of 1 7 km which illuminates the conduit along which magma was transported to the surface a third event type named hybrid which blends the characteristics of both vt and lp events originates just below the lp source and may reflect brittle failure along a zone intersecting a fluid filled crack the distribution of hybrid events is elongated 0 2 0 4 km in an east west direction this distribution may offer constraints on the orientation and size of the fluid filled crack inferred to be the source of the lp events

bayesian statistics applied to the location of the source of explosions at stromboli volcano italy abstract we present a method for determining the location and spatial extent of the source of explosions at stromboli volcano italy based on a bayesian inversion of the slowness vector derived from frequency slowness analyses of array data the method searches for source locations that minimize the error between the expected and observed slowness vectors for a given set of model parameters the conditional probability density function of slowness vectors is approximated by a gaussian distribution of expected errors the method is tested with synthetics using a five layer velocity model derived for the north flank of stromboli and a smoothed velocity model derived from a power law approximation of the layered structure application to data from stromboli allows for a detailed examination of uncertainties in source location due to experimental errors and incomplete knowledge of the earth model although the solutions are not constrained in the radial direction excellent resolution is achieved in both transverse and depth directions under the assumption that the horizontal extent of the source does not exceed the crater dimension the 90 confidence region in the estimate of the explosive source location corresponds to a small volume extending from a depth of about 100 m to a maximum depth of about 300 m beneath the active vents with a maximum likelihood source region located in the 120 to 180 m depth interval

array measurements and modeling of sources of shallow volcanic tremor at kilauea volcano hawaii we show that resonances in near surface layers path effects or a combination of near surface resonances and source resonances source effects could be responsible for the sharp spectral peaks found in the records of gas piston activity and volcanic tremor recorded near puu oo crater on kilauea volcano hawaii two types of sources are found to be compatible with observations and can not be distinguished using our data in one model impulsive explosion point sources are embedded in the puu oo structure in the other model impulsive point sources of pressure excite the resonances of a fluid filled crack embedded in that structure both models require shallow sources z 100 m and display a strong dependence of the radiated spectrum on source depth over a depth range comparable to the overall thickness of the surficial layers based on these results it may be possible to track changes in source depth through temporal changes in observed spectra a potentially useful tool for monitoring volcanic activity at this site we estimate the depth and spatial extent of volcanic tremor and gas piston activity using two dense arrays with respective apertures of 800 and 120 m located near puu oo measurements of slowness ray parameter and azimuth as a function of time clearly indicate that the sources of volcanic tremor and gas piston activity are located beneath or in close proximity to the puu oo crater at depths shallower than approximately 1 km based on slowness and particle motion analyses we find that the records of volcanic tremor and gas piston events at puu oo consist of a complex combination of body and surface waves

observations of high frequency scattered waves using dense arrays at teide volcano abstract a study of the seismic coda of seven small earthquakes recorded on the teide volcano canary islands spain was carried out using a temporary small aperture 12 station seismic array the purpose was to measure backazimuth apparent velocity and the type of waves that compose the coda in the frequency range 4 to 6 hz we used the zero lag cross correlation zlc method to obtain the components of the wave vector and three component analysis techniques based on the covariance matrix of the signal in the time domain to investigate the polarization properties of the signals the results show that a great part of the coda signals in the analyzed frequency range are almost uncorrelated while a low number of isolated correlated arrivals show an apparent slowness between 2 10 3 and 2 5 10 4 sec m and an almost random distribution of backazimuths the correlated arrivals have been interpreted as generated by strong scatterers probably related to the presence in the study area of surface topography irregularities the volcanic cone and the caldera rim the wave type varies from pure s type waves to mixed surface waves with some rayleigh components