waveform classification of volcanic low frequency earthquake swarms and its implication at soufriere hills volcano montserrat low frequency volcanic earthquakes are indicators of magma transport and activity within shallow conduit systems at a number of volcanoes these events exhibit a high degree of waveform similarity providing a criterion for classification using cross correlation techniques to quantify the degree of similarity we develop a method to sort events into families containing comparable waveforms events within a family have been triggered within one small source volume from which the seismic wave has then travelled along an identical path to the receiver this method was applied to a series of 16 low frequency earthquake swarms well correlated with cyclic deformation recorded by tiltmeters at soufriere hills volcano montserrat in june 1997 nine waveform groups were identified containing more than 45 events each the families are repeated across swarms with only small changes in waveform indicating that the seismic source location is stable with time the low frequency seismic swarms begin prior to the point at which inflation starts to decelerate suggesting that the seismicity indicates or even initiates a depressurisation process a major dome collapse occurred within the time window considered removing the top 100 m of the dome this event caused activity within some families to pause for several cycles before reappearing this shows that the collapse did not permanently disrupt the source mechanism or the path of the seismic waves

an analytical study for 1 dimensional steady flow in volcanic conduits the effects of magma properties viscosity and volatile content and geological conditions radius and length of conduit on the dynamics of explosive eruptions are analytically investigated on the basis of a simplified model for 1 dimensional steady flow in a volcanic conduit the conduit flow is characterized by bubbly flow before magma fragmentation and by gas pyroclast flow after magma fragmentation the features of the flow are largely determined by the balance of viscous effects of the bubbly flow and expansion of the gas pyroclast flow under a given boundary condition those effects are described by two non dimensional numbers and the parameter represents the ratio of effects of wall friction and gravity in the bubbly flow the length of the bubbly flow region decreases as increases the parameter represents the ratio of pressures at the vent and at fragmentation level the length of the gas pyroclast flow region decreases as increases three regimes of conduit flow are defined according to when is smaller than a certain value determined by initial volatile content and length of conduit the conduit flow is composed mainly of the gas pyroclast flow region regime 1 on the other hand when 10 the conduit flow is composed mainly of the bubbly flow region regime 3 regime 2 is the intermediate regime between regimes 1 and 3 the flow variables such as level of fragmentation and magma discharge rate are sensitive to magmatic properties and geological conditions in regime 2 a graphical method to estimate the effects of these parameters is also proposed

source characterization for an explosion during the 2009 eruption of redoubt volcano from very long period seismic waves abstract the 2009 eruption of redoubt produced several very long period vlp signals associated with explosions we invert for the source location and mechanism of an explosion at redoubt volcano using waveform methods applied to broadband recordings such characterization of the source carries information on the geometry of the conduit and the physics of the explosion process inversions are carried out assuming the volcanic source can be modeled as a point source with mechanisms described by a a set of 3 orthogonal forces b a moment tensor consisting of force couples and c both forces and moment tensor components we find that the source of the vlp seismic waves during the explosion is well described by either a combined moment force source located northeast of the crater and at an elevation of 1 6 km asl or a moment source at an elevation of 800 m to the southwest of the crater the moment tensors for the solutions with moment and force and moment only share similar characteristics the source time functions for both moment tensors begin with inflation pressurization and execute two cycles of deflation reinflation depressurization repressurization although the moment force source provides a better fit to the data we find that owing to the limited coverage of the broadband stations at redoubt the moment only source is the more robust and reliable solution based on the moment only solution we estimate a volume change of 19 000 m3 and a pressure change of 7 mpa in a dominant sill and an out of phase volume change of 5000 m3 and pressure change of 1 8 mpa in a subdominant dike at the source location these results shed new light on the magmatic plumbing system beneath redoubt and complement previous studies on vulcanian explosions at other volcanoes

source mechanism of long period events recorded by a high density seismic network during the 2008 eruption on mount etna 1 one hundred twenty nine long period lp events divided into two families of similar events were recorded by the 50 stations deployed on mount etna in the second half of june 2008 during this period lava was flowing from a lateral fracture after a summit strombolian eruption in order to understand the mechanisms of these events we perform moment tensor inversions inversions are initially kept unconstrained to estimate the most likely mechanism numerical tests show that unconstrained inversion leads to reliable moment tensor solutions because of the close proximity of numerous stations to the source positions however single forces cannot be accurately determined as they are very sensitive to uncertainties in the velocity model constrained inversions for a crack a pipe or an explosion then allow us to accurately determine the structural orientations of the source mechanisms both numerical tests and lp event inversions emphasise the importance of using stations located as close as possible to the source inversions for both families show mechanisms with a strong volumetric component these events are most likely generated by cracks striking sw ne for both families and dipping 70 se family 1 and 50 nw family 2 for family 1 events the crack geometry is nearly orthogonal to the dikelike structure along which events are located while for family 2 the location gave two pipelike bodies that belong to the same plane as the crack mechanism the orientations of the cracks are consistent with local tectonics which shows a sw ne weakness direction the lp events appear to be a response to the lava fountain occurring on 10 may 2008 as opposed to the flank lava flow

shallow s wave attenuation and actively degassing magma beneath taal volcano philippines taal volcano philippines is one of the world s most dangerous volcanoes given its history of explosive eruptions and its close proximity to populated areas a real time broadband seismic network was recently deployed and has detected volcano tectonic events beneath taal our source location analysis of these volcano tectonic events using onset arrival times and high frequency seismic amplitudes points to the existence of a region of strong attenuation near the ground surface beneath the east flank of volcano island in taal lake this region is beneath the active fumarolic area and above sources of pressure contributing inflation and deflation and it coincides with a region of high electrical conductivity the high attenuation region matches that inferred from an active seismic survey conducted at taal in 1993 these features strongly suggest that the high attenuation region represents an actively degassing magma body near the surface that has existed for more than 20 years

satellite observations of a western boundary current in the bay of bengal satellite infrared observations of the bay of bengal during the latter part of february 1985 reveal the existence of two bands of warm water that resemble western boundary currents along the east coasts of india and sri lanka a recently formed elliptical warm core eddy with a major axis of nearly 150 km appears at longitude 90 e and latitude 19 n at the end of the axis of the current color infrared images of sea surface temperature are used to illustrate the sea surface temperature patterns and to estimate their displacements with time

combined electromagnetic geochemical and thermal surveys of taal volcano philippines during the period 2005 2006 since 1572 33 phreatic to phreatomagmatic eruptions have occurred on taal volcano philippines some of them causing several hundred casualties considering the time delay between two consecutive eruptions there is an 88 probability that taal volcano should have already erupted since 1992 several phases of seismic activity have been recorded accompanied by ground deformation opening of fissures and surface activity the volcanic activity of taal appears to be controlled by dike injections and magma supply buffered by a hydrothermal system that releases fluids and heat through boiling and subsequent steaming in early 2005 a multidisciplinary project was launched for studying the hydrothermal activity to map the hydrothermal system combined surveys were carried out to investigate self potential total magnetic field ground temperature and carbon dioxide soil degassing along with satellite thermal imaging of the main crater lake the elevated temperatures and high concentrations of carbon dioxide as well as electromagnetic anomalies indicate large scale hydrothermal degassing this process is enhanced along the tectonic features e g crater rim and faults of the volcano while active fissures opened along the e w northern flank during the 1992 1994 seismic activity heat and fluids from the hydrothermal system are essentially released in the northern part of the crater which is bounded to the south by a suspected nw se fault along which seismicity seems to take place and dikes are thought to be intruded during the january 2005 surveys a new seismic crisis started and the felt earthquakes prompted spontaneous evacuation of hundreds of inhabitants living on the volcano repeated surveys show changes of self potential total magnetic field and ground temperature with time without any noticeable spatial enlargement these observations suggest that the northern flank located between the crater rim and the 1992 1994 fissures is connected with a deep thermal source in main crater and is reactivated during seismic crises this sector could be subjected to flank failure

source mechanism of volcanic tremor low frequency 10 hz volcanic earthquakes originate at a wide range of depths and occur before during and after magmatic eruptions the characteristics of these earthquakes suggest that they are not typical tectonic events physically analogous processes occur in hydraulic fracturing of rock formations low frequency icequakes in temperate glaciers and autoresonance in hydroelectric power stations we propose that unsteady fluid flow in volcanic conduits is the common source mechanism of low frequency volcanic earthquakes tremor the fluid dynamic source mechanism explains low frequency earthquakes of arbitrary duration magnitude and depth of origin as unsteady flow is independent of physical properties of the fluid and conduit fluid transients occur in both low viscosity gases and high viscosity liquids a fluid transient analysis can be formulated as generally as is warranted by knowledge of the composition and physical properties of the fluid material properties geometry and roughness of the conduit and boundary conditions to demonstrate the analytical potential of the fluid dynamic theory we consider a single phase fluid a melt of mount hood andesite at 1250 c in which significant pressure and velocity variations occur only in the longitudinal direction further simplification of the conservation of mass and momentum equations presents an eigenvalue problem that is solved to determine the natural frequencies and associated damping of flow and pressure oscillations for a simple constant pressure reservoir conduit orifice fluid system a change in orifice size can cause the source signal duration due to a single disturbance to change from several seconds to in excess of an hour fluid kinematic viscosity 0 412 m2 s is both included and neglected in the analysis to illustrate its effect upon system damping tremor magnitude is related to the magnitude of the conduit wall motion occurring in response to oscillating fluid pressures the entire volcanic fluid system does not generally experience a transient flow condition at a given time closed end and constant pressure fluid system boundaries reflect pressure waves and isolate components of the fluid system total fluid movement may not be accurately monitored by tremor data alone since steady state and slowly changing flows are aseismic

source mechanisms of explosions at stromboli volcano italy determined from moment tensor inversions of very long period data 1 seismic data recorded in the 2 30 s band at stromboli volcano italy are analyzed to quantify the source mechanisms of strombolian explosions during september 1997 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 elastic medium that takes topography into account two source centroids are identified each representative of the distinct event types associated with explosive eruptions from two different vents the observed waveforms are well reproduced by our inversion and the two source centroids that best fit the data are offset 220 and 260 m beneath and 160 m northwest of the active vents the source mechanisms include both moment tensor and single force components the principal axes of the moment tensor have amplitude ratios 1 1 2 which can be interpreted as representative of a crack if one assumes the rock matrix at the source to have a poisson ratio n 1 3 a value appropriate for hot rock both imaged cracks dip 60 to the northwest and strike northeast southwest along a direction parallel to the elongation of the volcanic edifice and a prominent zone of structural weakness as expressed by lineaments dikes and brittle structures for our data set the volume changes estimated from the moments are 200 m 3 for the largest explosion from each vent together with the volumetric source is a dominantly vertical force with a magnitude of 10 8 n consistent with the inferred movement of the magma column perched above the source centroid in response to the piston like rise of a slug of gas in the conduit index terms 7215 seismology earthquake parameters 7280 seismology volcano seismology 8419 8414 volcanology eruption mechanisms keywords very long period seismicity moment tensor inversions eruption mehanics

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

interpretation of resonance frequencies recorded during hydraulic fracturing treatments hydraulic fracturing treatments are often monitored by strings of geophones deployed in boreholes instead of picking discrete events only we here use time frequency representations of continuous recordings to identify resonances in two case studies this paper outlines an interpretational procedure to identify their cause using a subdivision into source path and receiver side effects for the first case study two main resonances are observed both at depth by the downhole geophones and on the surface by two broadband arrays the two acquisition networks have different receiver and path effects yet recorded the same resonances these resonances are therefore likely generated by source effects the amplitude pattern at the surface arrays indicates that these resonances are probably due to pumping operations in the second case study selective resonances are detected by the downhole geophones resonances coming from receiver effects are either lower or higher frequency and wave propagation modeling shows that path effects are not significant we identify two possible causes within the source area namely eigenvibrations of fractures or non darcian flow within the hydraulic fractures in the first situation 15 10 m long fluid filled cracks could generate the observed resonances an interconnected fracture network would then be required corresponding to mesoscale deformation of the reservoir alternatively systematic patterns in non darcian fluid flow within the hydraulic fracture could also be their leading cause resonances can be used to gain a better understanding of reservoir deformations or dynamic fluid flow perturbations during fluid injection into hydrocarbon and geothermal reservoirs co2 sequestration or volcanic eruptions

source mechanism of vulcanian degassing at popocatepetl volcano mexico determined from waveform inversions of very long period signals 1 the source mechanism of very long period vlp signals accompanying volcanic degassing bursts at popocatepetl is analyzed in the 15 70 s band by minimizing the residual error between data and synthetics calculated for a point source embedded in a homogeneous medium the waveforms of two eruptions 23 april and 23 may 2000 representative of mild vulcanian activity are well reproduced by our inversion which takes into account volcano topography the source centroid is positioned 1500 m below the western perimeter of the summit crater and the modeled source is composed of a shallow dipping crack sill with easterly dip of 10 intersecting a steeply dipping crack northeast striking dike dipping 83 northwest whose surface extension bisects the vent both cracks undergo a similar sequence of inflation deflation and reinflation reflecting a cycle of pressurization depressurization and repressurization within a time interval of 3 5 min the largest moment release occurs in the sill showing a maximum volume change of 500 1000 m3 pressure drop of 3 5 mpa and amplitude of recovered pressure equal to 1 2 times the amplitude of the pressure drop in contrast the maximum volume change in the dike is less 200 300 m3 with a corresponding pressure drop of 1 2 mpa and pressure recovery equal to the pressure drop accompanying these volumetric sources are single force components with magnitudes of 108 n consistent with melt advection in response to pressure transients the source time histories of the volumetric components of the source indicate that significant mass movement starts within the sill and triggers a mass movement response in the dike within a few seconds such source behavior is consistent with the opening of a pathway for escape of pent up gases from slow pressurization of the sill driven by magma crystallization the opening of this pathway and associated rapid evacuation of volcanic gases induces the pressure drop pressure recovery in the magma filling the sill is driven by diffusion of gases from the resulting supersaturated melt into bubbles assuming a penny shaped crack at ambient pressure of 40 mpa the observed pressure and volume variations can be modeled with the following attributes crack radius 100 m crack aperture 5 m bubble number density 1010 1012 m 3 initial bubble radius 10 6 m final bubble radius 10 5 m and net decrease of gas concentration in the melt 0 01 wt

site amplification from s wave coda in the long valley caldera region california following the work of phillips 1985 we have computed site amplification factors for coda waves at many sites in the long valley region in the eastern sierra nevada we computed ratios of coda amplitudes measured at 15 stations in and around long valley caldera relative to a granitic site mmpm for six frequency bands centered at 1 5 3 0 6 0 9 0 12 0 and 15 0 hz all station sites located within the caldera experienced large ground motion amplification at 1 5 and 3 0 hz ranging between five and 17 times that of the reference site however at higher frequencies these same sites exhibited significantly less amplification than the reference granite site this is attributed to the competing effects of an impedance contrast between the basement rock and caldera fill and very high absorption in the caldera fill at high frequencies station mmlm located on top of a volcanic plug displayed the largest amplitudes of all the sites studied for frequencies between 9 0 and 15 0 hz a dike structure attached to the plug couples the basement rock to the surface at high frequencies the resulting large amplitudes at mmlm are not due to amplification resulting from a strong impedance contrast rather the absorption under this site is very low perhaps lower than at the reference site mmpm outside the caldera another hard rock site located at devil s postpile mdpm generally behaved like the reference site for all frequencies the lowest amplifications observed came from a site outside the caldera mdcm located on thin pyroclastic ash deposits overlying granitic basement this can be attributed to a dominance of absorption over the amplification caused by lower impedance of this layer variations among sites on similar surface geology may be due to small local variations in impedance and absorption under and adjacent to the site the range in the spectral decay parameter between caldera and rock sites are comparable to results of anderson and hough 1984 for sites on alluvium and rock in the san fernando region these surprisingly different amplifications support the need for additional site specific studies amplifications determined in this study for the frequency range 1 5 and 3 0 hz correlate remarkably well with eaton s 1990 residuals for duration magnitude fmag and amplitude magnitude xmag for the usgs northern california seismic array further supporting the use of coda waves in determining site specific amplification

a preliminary seismic study of taal volcano luzon island philippines abstract the very active taal volcano lies in the southern part of luzon island only 60 km from manila the capital of the philippines in march 2008 we deployed a temporary seismic network around taal that consisted of 8 three component short period seismometers this network recorded during the period from march to november 2008 about 1050 local events in the early data processing stages unexpected linear drifting of clock time was clearly identified for a number of stations the drifting rates of each problematic station were determined and the errors were corrected before further processing initial location of each event was derived by manually picked p s phases arrival times using hypo71 and a general velocity model based on ak135 since the velocity structure beneath taal is essentially unknown we used travel times of 338 well located events in order to derive a minimum 1d velocity model using velest the resulting locations show that most events occurred at the shallow depth beneath the taal volcano and two major earthquake groups were noticed with one lying underneath the western shore of taal lake and the other one spread around the eastern flank of the taal volcano since there is no reported volcano activities during the operation period of our seismic array we are still not confident to interpret these findings in terms of other natures of volcano at the current stage however our work represents an important pioneer step towards other more advanced seismic studies in taal volcano

shallow conduit system at kilauea volcano hawaii revealed by seismic signals associated with degassing bursts 1 eruptive activity at the summit of kilauea volcano hawaii beginning in march 2008 and continuing to the present time is characterized by episodic explosive bursts of gas and ash from a vent within halemaumau pit crater these bursts are accompanied by seismic signals that are well recorded by a broadband network deployed in the summit caldera we investigate in detail the dimensions and oscillation modes of the source of a representative burst in the 1 10 s band an extended source is realized by a set of point sources distributed on a grid surrounding the source centroid where the centroid position and source geometry are fixed from previous modeling of very long period vlp data in the 10 50 s band the source time histories of all point sources are obtained simultaneously through waveform inversion carried out in the frequency domain short scale noisy fluctuations of the source time histories between adjacent sources are suppressed with a smoothing constraint whose strength is determined through a minimization of the akaike bayesian information criterion abic waveform inversions carried out for homogeneous and heterogeneous velocity structures both image a dominant source component in the form of an east trending dike with dimensions of 2 9 2 9 km the dike extends 2 km west and 0 9 km east of the vlp centroid and spans the depth range 0 2 3 1 km the source model for a homogeneous velocity structure suggests the dike is hinged at the source centroid where it bends from a strike e 27 n with northern dip of 85 west of the centroid to a strike e 7 n with northern dip of 80 east of the centroid the oscillating behavior of the dike is dominated by simple harmonic modes with frequencies 0 2 hz and 0 5 hz representing the fundamental mode 11 and first degenerate mode 12 21 of the dike although not strongly supported by data in the 1 10 s band a north striking dike segment is required for enhanced compatibility with the model elaborated in the 10 50 s band this dike provides connectivity between the east trending dike and the new vent within halemaumau pit crater waveform inversions with a dual dike model suggest dimensions of 0 7 0 7 km to 2 6 2 6 km for this segment further elaboration of the complex dike system under halemaumau does not appear to be feasible with presently available data

spatiotemporal structure of the east india coastal current from satellite altimetry we use a newly processed altimeter data set to present a hitherto unprecedented description of the spatiotemporal structure of the east india coastal current eicc the data set resolves timescales ranging from a few months to a few years and the high along track resolution yields the first description of the cross shore structure of the current the seasonal cycle dominates the variability but the nonannual timescales have similar energy levels all along the eicc path there are short lived intense intraseasonal bursts in contrast to the seasonal cycle the interannual and intraseasonal components are decorrelated along the coast and possible mechanisms for the decorrelation discontinuity in the flow are discussed in the cross shore direction the current is highly correlated at all timescales the eicc is trapped against the shelf with the current offshore flowing in the opposite direction at most locations the eicc appears as an inherently discontinuous flow taking the form of a few recirculating loops along the eicc path with a typical cross shore spatial scale of 150 200 km the loops are highly variable in direction at all timescales from intraseasonal to interannual this discontinuity of the eicc in space and time implies that the basic pathways and advective timescales for the interbasin exchange of water masses between the bay of bengal and the arabian sea are not robust when the full spatiotemporal variability of the eicc is considered

satellite thermal analysis of the 1986 izu oshima lava flows abstract surface temperatures of the eruption products and fumaroles associated with the 1986 eruption of izu oshima volcano japan were analysed using landsat thematic mapper tm imagery acquired between about one and eight months after eruption cessation field observations and video footage of the eruption assisted the analysis throughout the post eruptive period studied the eruptive products showed elevated radiance in the tm thermal infrared channel band 6 10 4 12 5 m equating to pixel integrated temperatures of 6 48 c the tm shortwave infrared channels band 5 1 55 1 75 m and band 7 2 08 2 35 m are sensitive just to higher temperature objects and thermal signatures were found only at the fumarolic areas where temperatures estimated with the dual band method were a maximum of 700 800 c in tm band 6 a systematic surface temperature difference which decreased with time after the eruption was observed between two clastogenic lava flows lb i and lb iii that erupted at shortly spaced intervals on 21 november 1986 16 20 22 00 h local time potential reasons for this temperature difference are explored with the conclusion that the cause is a difference in the explosivity which relates to fountaining height and degrees of magma fragmentation of the two fire fountains rather than differences in magmatic flow temperature surface structure or any non volcanic factor explosivity was greater for the fire fountain feeding lb iii than lb i thereby causing a more extensive heat loss during the airborne passage of the lb iii fragments due to their on average more fragmented nature and longer flight times thus the temperature of the hot core of the lb i flow remained higher than that of lb iii and this was reflected in the surface temperature of the upper clinker of the flow when imaged by landsat tm

origin of the pulse like signature of shallow long period volcano seismicity short duration pulse like long period lp events are a characteristic type of seismicity accompanying eruptive activity at mount etna in italy in 2004 and 2008 and at turrialba volcano in costa rica and ubinas volcano in peru in 2009 we use the discrete wave number method to compute the free surface response in the near field of a rectangular tensile crack embedded in a homogeneous elastic half space and to gain insights into the origin of the lp pulses two source models are considered including 1 a vertical fluid driven crack and 2 a unilateral tensile rupture growing at a fixed sub rayleigh velocity with constant opening on a vertical crack we apply cross correlation to the synthetics and data to demonstrate that a fluid driven crack provides a natural explanation for these data with realistic source sizes and fluid properties our modeling points to shallow sources 1 km while a slow rupture failure provides another potential model for these events the synthetics and resulting fits to the data are not optimal in this model compared to a fluid driven source we infer that pulse like lp signatures are parts of the continuum of responses produced by shallow fluid driven sources in volcanoes

the source of infrasound associated with long period events at mount st helens 1 during the early stages of the 2004 2008 mount st helens eruption the source process that produced a sustained sequence of repetitive long period lp seismic events also produced impulsive broadband infrasonic signals in the atmosphere to assess whether the signals could be generated simply by seismic acoustic coupling from the shallow lp events we perform finite difference simulation of the seismo acoustic wavefield using a single numerical scheme for the elastic ground and atmosphere the effects of topography velocity structure wind and source configuration are considered the simulations show that a shallow source buried in a homogeneous elastic solid produces a complex wave train in the atmosphere consisting of p sv and rayleigh wave energy converted locally along the propagation path and acoustic energy originating from the source epicenter although the horizontal acoustic velocity of the latter is consistent with our data the modeled amplitude ratios of pressure to vertical seismic velocity are too low in comparison with observations and the characteristic differences in seismic and acoustic waveforms and spectra cannot be reproduced from a common point source the observations therefore require a more complex source process in which the infrasonic signals are a record of only the broadband pressure excitation mechanism of the seismic lp events the observations and numerical results can be explained by a model involving the repeated rapid pressure loss from a hydrothermal crack by venting into a shallow layer of loosely consolidated highly permeable material heating by magmatic activity causes pressure to rise periodically reaching the pressure threshold for rupture of the valve sealing the crack sudden opening of the valve generates the broadband infrasonic signal and simultaneously triggers the collapse of the crack initiating resonance of the remaining fluid subtle waveform and amplitude variability of the infrasonic signals as recorded at an array 13 4 km to the nw of the volcano are attributed primarily to atmospheric boundary layer propagation effects superimposed upon amplitude changes at the source

a synthesis and review of historical eruptions at taal volcano southern luzon philippines abstract the philippines is an area of persistent volcanism being located in one of the most tectonically active regions in the world taal volcano in southern luzon is the second most frequently erupting volcano of the 24 active volcanoes in the philippines a comprehensive and critical review of published and unpublished references describing the 33 known historical eruptions of taal may provide answers to knowledge gaps on past eruptive behavior processes and products that could be utilized for hazard and risk assessment of future eruptions data on the prehistoric eruptions and evolution of taal caldera and subsequent deposits are limited only four caldera forming events were identified based on four mapped ignimbrite deposits from oldest to youngest these are the silicic alitagtag ali and caloocan cal pumice flow deposits the dacitic sambong ignimbrite sam and the basaltic andesitic taal scoria flow renamed scoria pyroclastic flow sfl except for sfl with 14c dating yielding 5380 70 to 6830 80 ky there are no age constraints or estimates of extent for the three older deposits a comprehensive review of the historical eruptions of taal volcano is the central element of this paper and includes all eruptions from ad1572 the first known historic event to ad1977 eruption styles and the interplay between processes and products for each eruption are reinterpreted based on the narrative descriptions from all available accounts a change of classification of eruption styles and eruptive products is undertaken for some events at least nine reported eruptions were deemed uncertain including the ad1605 ad1611 event more likely seismic swarms the ad1634 ad1635 and ad1645 may simply be solfataric or hydrothermal activity events and the ad1790 ad1825 ad1842 ad1873 and ad1903 events that were listed in recent published and unpublished documents but do not provide any details to describe and confirm the eruptions except for listing a default vei of 2 pyroclastic density currents brought devastating impacts to the communities around taal during the ad1749 ad1754 ad1911 and ad1965 eruptions and remain the biggest threat in the case of renewed volcanic activity significant implications for aviation are implied by the narrative of tephra fall dispersal towards manila the central gateway of international aviation operation in the philippines during the ad1754 eruptions the dispersal of tephra in the event of an explosive eruption at taal towards metro manila would have catastrophic effects to transport utilities and business activity potentially generating enormous economic losses hazards from earthquake events associated with future volcanic activity may also have localized impacts occurrences of liquefaction phenomena as a consequence of severe ground shaking are interpreted during the ad1749 ad1754 and ad1911 eruptions more work needs to be done to develop a comprehensive understanding of the hazards and risks associated with an explosive eruption at taal volcano especially related to the older quaternary caldera forming eruptions that produced large volume pyroclastic deposits that are extensively distributed and exposed we acknowledge that there may be additional prehistoric eruptions where the eruptive products have not been preserved recognized or reported events that cannot be verified or do not have sufficient details to confirm the eruption have been downgraded to uncertain eruptions that are confirmed with identified dispersal and emplacement of tephra fall and other eruptive deposits as interpreted from narrated records could provide crucial information that may be utilized in hazard assessment

does a monsoon climate exist over south america the climatology and the basic state of the summertime circulation and rainfall over south america are studied using assimilation products from the data assimilation system of goddard earth observing system 1 geos1 and satellite derived rainfall results indicate the existence of a regional summer monsoon circulation regime induced by strong diabatic heating over the subtropical south american highland centered at the altiplano plateau sensitivity of the results to the assimilation scheme is tested by comparing that with the national centers for environmental prediction ncep reanalysis and with satellite rainfall estimates results show general agreement between the model produced rainfall anomaly and the satellite estimates as well as consistency between the basic circulation features in the geos 1 and the ncep reanalyses a case study of 1989 90 south american summer monsoon sasm reveals the following characteristics 1 in late spring the onset of sasm is signaled by an abrupt merging of the upper tropospheric double westerly jets one in the subtropics and the other in the subpolar region into a single jet in the midlatitudes this is followed by the establishment of a vortex to the southeast of altiplano and occurrence of heavy precipitation over subtropical eastern brazil 2 during the mature phase of sasm the heavy rainfall zone moves over the altiplano plateau and the southernmost brazilian highland the fully established sasm features are the following a an enhancement of equatorial north atlantic trade wind which emanates from the sahara high and crosses the equator over the south american continent b a buildup of strong northwesterlies along the eastern side of the tropical andes and c development of the south atlantic convergence zone in the southernmost position with strong convective activity meanwhile the upper tropospheric return flow emerges from an anticyclone formed over the altiplano plateau crosses the equator and sinks over northwestern africa 3 the withdrawal of sasm in late summer is signaled by the resplitting of the midlatitude westerly jet at the same time the low level northwest monsoon flow diminishes reducing the moisture supply and leading to the termination of heavy precipitation over the subtropical highland results also show that the above mentioned characteristics of sasm are clearly linked to the tropospheric temperature changes over the central south american highland sensible versus latent heating over the highland are bound to play an important role in the evolution of sasm to provide further support of presence of a monsoon climate over south america sasm is compared and contrasted to the classic east asian summer monsoon easm many similar features including evolution characteristics between the two systems have been identified contrasting aspects of the sasm from the easm are also discussed it is pointed out that a number of monsoonal characteristics of the climate of south america such as the seasonal reversal of the low level wind become apparent only when the strong annual mean wind is removed based on the characteristic features and their evolution the authors conclude that a monsoon climate does exist over south america

determining the seismic source mechanism and location for an explosive eruption with limited observational data augustine volcano alaska 1 waveform inversions of the very long period components of the seismic wavefield produced by an explosive eruption that occurred on 11 january 2006 at augustine volcano alaska constrain the seismic source location to near sea level beneath the summit of the volcano the calculated moment tensors indicate the presence of a volumetric source mechanism systematic reconstruction of the source mechanism shows the source consists of a sill intersected by either a sub vertical east west trending dike or a sub vertical pipe and a weak single force the trend of the dike may be controlled by the east west trending augustine seldovia arch the data from the network of broadband sensors is limited to fourteen seismic traces and synthetic modeling confirms the ability of the network to recover the source mechanism the synthetic modeling also provides a guide to the expected capability of a broadband network to resolve very long period source mechanisms particularly when confronted with limited observational data

site amplification of coda waves from local earthquakes in central california the goals of this study are two fold 1 to examine the site effect of coda waves believed to be the average site effect of shear waves in order to understand its spatial and frequency dependent behavior and 2 to learn what we can about the processes that generate the coda itself we use digital data from over 90 earthquakes each recorded by a subset of some 150 stations in the coast ranges of central california between san francisco and san luis obispo results from the band 1 5 to 24 hz indicate that site amplification depends strongly on surface geology and frequency at low frequencies 1 5 to 3 0 hz changes in amplification up to a factor of 20 are observed amplification generally varies inversely with deposit age at sediment sites and is lowest at granite and franciscan basement sites at high frequencies 6 to 24 hz the pattern changes many granite sites the majority in the gabilan range exhibit increasing amplification with frequency relative to the average station this behavior differs strikingly from that observed at sediment and franciscan sites which roll off at an intermediate rate and many sediment sites adjacent to the san andreas fault zone which decay away even more rapidly with frequency these results can be explained qualitatively by appealing to variable near site impedance and attenuation however complicating phenomena such as mode excitation in valley sediments must exist at low frequencies peaked behavior at a few hard rock sites at high frequency 12 hz currently lacks explanation at a few sediment sites the low frequency excitation is so strong that the coda envelopes begin to take on different shapes this contradicts the assumption based on previous observations that the coda shape is independent of source receiver location thus some site measurements may be artificially high the working model of coda waves as body wave energy backscattered from randomly distributed inhomogeneities must be modified to include the possibility of efficient near site resonance excited by incident direct and coda wave energy

dynamics of active magmatic and hydrothermal systems at taal volcano philippines from continuous gps measurements 1 a dense network of continuous single and dual frequency gps receivers at taal volcano philippines reveals four major stages of volcanogenic deformation deflation from installation in june 1998 to december 1998 inflation from january to march 1999 deflation from april 1999 to february 2000 and inflation from february to november 2000 the largest displacements occurred during the february november 2000 period of inflation which was characterized by 120 mm of uplift of the center of volcano island relative to the northern caldera rim at average rates up to 216 mm yr deformation sources were modeled using a constrained least squares inversion algorithm the source of 1999 deflation and inflation in 2000 were modeled as contractional and dilatational mogi point sources centered at 4 2 and 5 2 km depth respectively beneath volcano island the locations of the inflationary and deflationary sources are indistinguishable within the 95 confidence estimates modeling using a running 4 month time window from june 1999 to march 2001 reveals little evidence for source migration we suggest that the two periods of inflation observed at taal result from episodic intrusions of magma into a shallow reservoir centered beneath volcano island subsequent deflation may result from exsolution of magmatic fluids and or gases into an overlying unconfined hydrothermal system

improved identification of volcanic features using landsat 7 etm volcanic eruptions can present unpredictable hazards to populations living within regions containing potentially active volcanoes and for people traveling in jet aircraft that intersect with ash laden eruption clouds methods of monitoring volcanic activity include searching for variations in the thermal signal from active fumaroles lava domes lava lakes flows and other features over many active volcanoes in the western hemisphere low spatial resolution 4 km pixel weather satellite data acquired every 15 min are used to identify changes in eruptive activity but are of insufficient spatial resolution to map active volcanic features the enhanced thematic mapper plus etm on landsat 7 can be used to monitor active volcanoes at a higher spatial resolution 15 to 60 m pixels etm also offers improvements over its predecessor the thematic mapper of landsats 4 and 5 by way of a 15 m panchromatic band and higher spatial resolution 60 m pixel thermal infrared ir band with higher spatial resolution panchromatic data we are able to map lava flow fields trace very high temperature lava channels and at lascar volcano identify an arcuate feature associated with a collapsed crater floor a phenomenon that may precede explosive activity with improved spatial resolution in the thermal ir we are able to map the bifurcation and braiding of underground lava tubes at kilauea identifying tube locations and tracking their extension are important because for a given volumetric lava production effusion rate tube fed flows can extend a much greater distance than surface flows at both kilauea and etna we are able to use the thermal data to estimate effusion rates an important parameter for assessing how far a flow is likely to extend and therefore the hazard it poses an improved etm data collection and distribution system includes a well formulated and executed long term acquisition plan and less expensive data that is available much faster than previously possible

eruption dynamics at mount st helens imaged from broadband seismic waveforms interaction of the shallow magmatic and hydrothermal systems 1 the current eruption at mount st helens is characterized by dome building and shallow repetitive long period lp earthquakes waveform cross correlation reveals remarkable similarity for a majority of the earthquakes over periods of several weeks stacked spectra of these events display multiple peaks between 0 5 and 2 hz that are common to most stations lower amplitude very long period vlp events commonly accompany the lp events we model the source mechanisms of lp and vlp events in the 0 5 4 s and 8 40 s bands respectively using data recorded in july 2005 with a 19 station temporary broadband network the source mechanism of the lp events includes 1 a volumetric component modeled as resonance of a gently nnw dipping steam filled crack located directly beneath the actively extruding part of the new dome and within 100 m of the crater floor and 2 a vertical single force attributed to movement of the overlying dome the vlp source which also includes volumetric and single force components is 250 m deeper and nnw of the lp source at the sw edge of the 1980s lava dome the volumetric component points to the compression and expansion of a shallow magma filled sill which is subparallel to the hydrothermal crack imaged at the lp source coupled with a smaller component of expansion and compression of a dike the single force components are due to mass advection in the magma conduit the location geometry and timing of the sources suggest the vlp and lp events are caused by perturbations of a common crack system

very low frequency resistivity self potential and ground temperature surveys on taal volcano philippines implications for future activity abstract taal volcano is one of the most dangerous volcanoes in the philippines thirty three eruptions have occurred through historical time with several exhibiting cataclysmic phases most recent eruptions are confined to volcano island located within the prehistoric taal collapse caldera that is now filled by taal lake the last eruptive activity from 1965 to 1977 took place from mt tabaro about 2 km to the southwest of the main crater center since this time episodes of seismic activity ground deformation gas release surface fissuring fumarole activity and temperature changes are recorded periodically around the main crater but no major eruption has occurred this period of quiescence is the third longest period without eruptive activity since 1572 in march 2010 a campaign based on very low frequency vlf resistivity surveys together with repeated surveys of self potential ground temperature and fissure activity was intensified and the results compared to a large scale electrical resistivity tomography experiment this work fortunately occurred before within and after a new seismovolcanic crisis from late april 2010 to march 2011 the joint analysis of these new data together with results from previous magnetotelluric soundings allows a better description of the electrical resistivity and crustal structure beneath the main crater down to a depth of several kilometers no indication of growth of the two geothermal areas located on both sides of the northern crater rim was apparent from 2005 to march 2010 these areas appear controlled by active fissures opened during the 1992 and 1994 crises that dip downward towards the core of the hydrothermal system located at about 2 5 km depth beneath the crater older mineralized fissures at lower elevations to the north of the geothermal areas also dip downward under the crater repeated self potential and ground temperature surveys completed between 2005 and 2015 show new geothermal and hydrothermal activity in the areas of these older mineralized fissures that occurred during the april 2010 to march 2011 seismovolcanic crisis this dramatically extends the geothermal activity further to the north on the volcano the occurrence of these newly activated fissures after a long period of quiescence and indications of inflation in mid 2010 under the north rim of the main crater suggests that new eruptive activity near the north rim of the crater could occur in the future

evidence of thermal driven processes triggering the 2005 2014 unrest at campi flegrei caldera abstract an accelerating process of ground deformation that began 10 years ago is currently affecting the campi flegrei caldera the deformation pattern is here explained with the overlapping of two processes short time pulses that are caused by injection of magmatic fluids into the hydrothermal system and a long time process of heating of the rock the short pulses are highlighted by comparison of the residuals of ground deformation fitted with an accelerating polynomial function with the fumarolic co 2 ch 4 and he ch 4 ratios which are good geochemical indicators of the arrival of magmatic gases the two independent datasets show the same sequence of five peaks with a delay of 200 days of the geochemical signal with respect to the geodetic signal the heating of the hydrothermal system which parallels the long period accelerating curve is inferred by temperature pressure gas geoindicators referring to a recent interpretation that relates variations in the fumarolic inert gas species to open system magma degassing we infer that the heating is caused by enrichment in water of the magmatic fluids and by an increment in their flux heating of the rock caused by magmatic fluids can be a central factor in triggering unrest at calderas

geochemical evidence of magma intrusion inferred from diffuse co 2 emissions and fumarole plume chemistry the 2010 2011 volcanic unrest at taal volcano philippines we report herein significant increases in the diffuse co2 emission through the lake and soils of the taal volcano main crater philippines during a volcano seismic unrest period in april 2010 june 2011 monitoring of co2 emissions from the taal main crater lake and soils in the fumarole areas started in 2008 with a collaborative project between iter and phivolcs in august 2010 the total diffuse co2 emission from the taal main crater lake increased almost three times from a previously measured baseline emission value 1 000 t day 1 to 2 716 54 t day 1 this increase occurred during anomalous seismic activity which started in april 2010 the highest co2 emission value in the crater lake was measured in march 2011 at 4 670 159 t day 1 when the volcano was still showing signs of unrest significant changes in the co2 so2 co2 h2o and so2 h2s ratios in the fumarole gases suggest the occurrence of a magmatic injection several weeks before the onset of the volcano seismic unrest period these observed changes in the co2 emission rate and the fumarole gas composition as the result of magmatic activity were likely caused by magma moving from a deep reservoir to a shallow reservoir

downflow amplitude decrease of ground penetrating radar reflections in base surge deposits ground penetrating radar gpr data were collected parallel to a flow direction in the base surge deposits of the ubehebe hydrovolcanic field death valley california these data showed a downflow amplitude decrease of the reflections that was displayed by computing fourier transforms and the average square of the sample points of the traces this trend probably reflects some aspects of the lateral facies variation of base surge deposits and could be explained by the concomitant downflow increase in sorting and decrease in grain size with a consequent reduction of the differences between the beds which in turn causes a reduction of the reflected energy these results suggest that although the gpr profiles are traditionally used merely as 2 d images of the subsurface stratigraphy they could be useful in obtaining information about the heterogeneous distribution of subsurface physical properties the downflow decrease of the amplitude of the signal could be used as a flow direction indicator when the position of the vent is unknown although other geologic information and field constraints are probably necessary

waveform inversion of volcano seismic signals assuming possible source geometries 1 we propose an approach to quantify the source of volcano seismic signals assuming possible source geometries such an assumption reduces the number of free parameters in a waveform inversion so we can quantify the source of these signals observed by a small number of seismic stations we test this method by using the synthesized waveforms radiated by a fluid filled crack and with seismic data observed at kusatsu shirane volcano we find that waveform data from 2 to 3 three component seismic stations surrounding the source enable us to estimate the source mechanism and location in the present approach

geochemical and isotopic evidence for seawater contamination of the hydrothermal system of taal volcano luzon the philippines the hydrologic structure of taal volcano has favored development of an extensive hydrothermal system whose prominent feature is the acidic main crater lake ph 3500 mg kg and low so4 cl ratios 0 3 the br cl ratio of most samples suggests incorporation of seawater into the hydrothermal system water and dissolved sulfate isotopic compositions reveal that the main crater lake and spring discharges are derived from a deep parent fluid t approximate to 300 degrees c which is a mixture of seawater volcanic water and lake taal water the volcanic end member is probably produced in the magmatic hydrothermal environment during absorption of high temperature gases into groundwater boiling and mixing of the parent water give rise to the range of chemical and isotopic characteristics observed in the thermal discharges incursion of seawater from the coastal region to the central part of the volcano is supported by the low water levels of the lakes and by the fact that lake taal was directly connected to the china sea until the sixteenth century the depth to the seawater meteoric water interface is calculated to be 80 and 160 m for the main crater lake and lake taal respectively additional data are required to infer the hydrologic structure of taal geochemical surveillance of the main crater lake using the so4 cl na k or mg cl ratio cannot be applied straightforwardly due to the presence of seawater in the hydrothermal system

subevents of long period seismicity implications for hydrothermal dynamics during the 2004 2008 eruption of mount st helens 1 one of the most striking aspects of seismicity during the 2004 2008 eruption of mount st helens msh was the precise regularity in occurrence of repetitive long period lp or drumbeat events over sustained time periods however this precise regularity was not always observed and at times the temporal occurrence of lp events became more random in addition accompanying the dominant lp class of events during the 2004 2008 msh eruption there was a near continuous randomly occurring series of smaller seismic events these subevents are not always simply small amplitude versions of the dominant lp class of events but appear instead to result from a separate random process only loosely coupled to the main lp source mechanism we present an analysis of the interevent time and amplitude distributions of the subevents using waveform cross correlation to separate lp events from the subevents we also discuss seismic tremor that accompanied the 8 march 2005 phreatic explosion event at msh this tremor consists of a rapid succession of lps and subevents triggered during the explosion in addition to broadband noise from the sustained degassing immediately afterward seismicity returned to the pre explosion occurrence pattern this triggering in relation to the rapid ejection of steam from the system and subsequent return to pre explosion seismicity suggests that both seismic event types originated in a region of the subsurface hydrothermal system that was 1 in contact with the reservoir feeding the 8 march 2005 phreatic explosion but 2 not destroyed or drained by the explosion event finally we discuss possible thermodynamic conditions in a pressurized hydrothermal crack that could give rise to seismicity pressure drop estimates for typical lp events are not generally large enough to perturb pure water in a shallow hydrothermal crack into an unstable state however dissolved volatiles such as co2 may lead to a more unstable system increasing the seismogenic potential of a hydrothermal crack subject to rapid heat flux the interaction of hydrothermal and magmatic systems beneath msh in 2004 2008 thus appears able to explain a wide range of observed phenomena including subevents lp events larger md 2 events and phreatic explosions

source process of a long period event at kilauea volcano hawaii summary we analyse a long period lp event observed by a dense seismic network temporarily operated at kilauea volcano hawaii in 1996 we systematically perform spectral analyses waveform inversions and forward modeling of the lp event to quantify its source process spectral analyses identify two dominant spectral frequencies at 0 6 and 1 3 hz with associated q values in the range 10 20 results from waveform inversions assuming six moment tensor and three single force components point to the resonance of a horizontal crack located at a depth of approximately 150 m near the northeastern rim of the halemaumau pit crater waveform simulations based on a fluid filled crack model suggest that the observed frequencies and q values can be explained by a crack filled with a hydrothermal fluid in the form of either bubbly water or steam the shallow hydrothermal crack located directly above the magma conduit may have been heated by volcanic gases leaking from the conduit the enhanced flux of heat raised the overall pressure of the hydrothermal fluid in the crack and induced a rapid discharge of fluid from the crack which triggered the acoustic vibrations of the resonator generating the lp waveform the present study provides further support to the idea that lp events originate in the resonance of a crack

source mechanism of long period events at kusatsu shirane volcano japan inferred from waveform inversion of the effective excitation functions abstract we investigate the source mechanism of long period lp events observed at kusatsu shirane volcano japan based on waveform inversions of their effective excitation functions the effective excitation function which represents the apparent excitation observed at individual receivers is estimated by applying an autoregressive filter to the lp waveform assuming a point source we apply this method to seven lp events the waveforms of which are characterized by simple decaying and nearly monochromatic oscillations with frequency in the range 1 3 hz the results of the waveform inversions show dominant volumetric change components accompanied by single force components common to all the events analyzed and suggesting a repeated activation of a sub horizontal crack located 300 m beneath the summit crater lakes based on these results we propose a model of the source process of lp seismicity in which a gradual buildup of steam pressure in a hydrothermal crack in response to magmatic heat causes repeated discharges of steam from the crack the rapid discharge of fluid causes the collapse of the fluid filled crack and excites acoustic oscillations of the crack which produce the characteristic waveforms observed in the lp events the presence of a single force synchronous with the collapse of the crack is interpreted as the release of gravitational energy that occurs as the slug of steam ejected from the crack ascends toward the surface and is replaced by cooler water flowing downward in a fluid filled conduit linking the crack and the base of the crater lake

magmatic dike resonances inferred from very long period seismic signals very long period vlp signals showing simple decaying harmonic oscillations with periods near 10 seconds and lasting for about 300 seconds were observed in association with an earthquake swarm that occurred beneath hachijo island japan results from the source mechanism analysis and waveform simulation based on a fluid filled crack model consistently point to the resonance of a dike filled with a basaltic magma as the source of the vlp signals thus vlp signals can be used to probe the state of the fluid and dynamic processes within a volcanic system

principles of electrography applied to self potential electrokinetic sources and hydrogeological applications 1 the electrical potential field passively recorded at the ground surface of the earth and termed self potential can be analyzed to determine the shape and the depth of the piezometric surface the coupling between hydraulic flow and electrical current density is electrokinetic in nature the electrokinetic coupling coefficient entering into the integral equation relating the depth of the water table to self potential signals is analyzed for various types of porous materials it is simply related to the electrical conductivity of the pore water in steady state conditions each element of the water table can be seen as an elementary dipole with an inclination locally perpendicular to the water table and strength proportional to the water table elevation then we propose three methods to obtain the shape and range of possible depths of the water table from the study of the self potential distribution recorded at the ground surface the nonuniqueness of the solution is removed if one knows either the electrokinetic coupling coefficient or the water table at one location and under the assumption of the spatial homogeneity of the electrokinetic coupling coefficient two field cases are discussed to show the success of the proposed methods for estimating the shape and depth of the water table at two different scales of investigations they concern the study of self potential signals associated with the shape of the water table in the vicinity of a pumping well and in the flank of the kilauea volcano

theoretical basis of the amplitude source location method for volcano seismic signals the source location method using high frequency seismic amplitudes based on the assumption of isotropic radiation of s waves has been used successfully to locate seismic events and tremor at volcanoes this amplitude source location asl method overcomes the limitations of traditional hypocenter determination methods that use onset arrival times and has great potential to improve volcano seismic monitoring and the investigation of source processes of volcano seismic signals however theoretical justification of the basic assumption used in the asl method has not been provided in previous studies in this study we tested the asl method by using seismic waveforms simulated with a finite difference method in realistic heterogeneous volcanic structures with intrinsic attenuation and topography our results showed that asl determinations were successful when using waveforms characterized by multiple scattering that can be approximated by the diffusion model we found that the energy solution of the diffusion model provided an amplitude distance relation that is similar to that used in the asl method our results suggest that the asl method is applicable to high frequency seismograms in highly heterogeneous media with transport mean free paths of 103 m or smaller and strong intrinsic attenuation these medium parameters are consistent with those estimated at various volcanoes our study validates the basic assumption of the asl method and justifies its use to locate the sources of high frequency seismic signals observed at volcanoes

a large hydrothermal reservoir beneath taal volcano philippines revealed by magnetotelluric resistivity survey 2d resistivity modeling taal volcano located in the southwestern part of luzon island philippines has frequently experienced catastrophic eruptions from both the main crater on volcano island and flank eruptions these eruptions have been magmatic phreatomagmatic and hydrothermal with the latter implying the existence of a large scale hydrothermal system beneath the volcano we conducted an electrical resistivity survey using the magnetotelluric method in order to identify the location and geometry of the hydrothermal reservoir and sealing cap rock two dimensional inversion using the observed data indicates four similar resistivity sections the structure at shallow depths corresponds to volcanic deposits and an aquifer below 1 km the structure features a relatively resistive zone beneath the main crater surrounded by a conductive shell we interpreted these to be a large hydrothermal reservoir with an impermeable cap rock sealing it recent ground deformation detected by gps measurements suggests that the hydrothermal reservoir is active the interpreted cap rock thins just beneath the main crater and could easily be destroyed by an imbalance in the hydrothermal system we conclude that this hydrothermal reservoir plays a significant role in driving catastrophic eruptions that begin with a hydrothermal explosion at the main crater

waveform inversion of oscillatory signatures in long period events beneath volcanoes 1 the source mechanism of long period lp events is examined using synthetic waveforms generated by the acoustic resonance of a fluid filled crack we perform a series of numerical tests in which the oscillatory signatures of synthetic lp waveforms are used to determine the source time functions of the six moment tensor components from waveform inversions assuming a point source the results indicate that the moment tensor representation is valid for the odd modes of crack resonance with wavelengths 2l n 2w n n 3 5 7 where l and w are the crack length and width respectively for the even modes with wavelengths 2l n 2w n n 2 4 6 a generalized source representation using higher order tensors is required although the efficiency of seismic waves radiated by the even modes is expected to be small we apply the moment tensor inversion to the oscillatory signatures of an lp event observed at kusatsu shirane volcano central japan our results point to the resonance of a subhorizontal crack located a few hundred meters beneath the summit crater lakes the present approach may be useful to quantify the source location geometry and force system of lp events and opens the way for moment tensor inversions of tremor

real time monitoring and massive inversion of source parameters of very long period seismic signals an application to stromboli volcano italy 1 we present a comprehensive processing tool for the real time analysis of the source mechanism of very long period vlp seismic data based on waveform inversions performed in the frequency domain for a point source a search for the source providing the best fitting solution is conducted over a three dimensional grid of assumed source locations in which the green s functions associated with each point source are calculated by finite differences using the reciprocal relation between source and receiver tests performed on 62 nodes of a linux cluster indicate that the waveform inversion and search for the best fitting signal over 100 000 point sources require roughly 30 s of processing time for a 2 min long record the procedure is applied to post processing of a data archive and to continuous automatic inversion of real time data at stromboli providing insights into different modes of degassing at this volcano

arabian sea and bay of bengal exchange of salt and tracers in an ocean model exchanges of mass and salt between the mixed layers in the arabian sea and the bay of bengal are examined in a model of the indian ocean using passive tracers as a tool to map the pathways inflow of high salinity water from the arabian sea into the bay of bengal is significant and occurs after the mature phase of the southwest monsoon freshwater transport out of the bay of bengal is southward throughout the year along the eastern boundary of the indian ocean low salinity water transport into the arabian sea occurs in the somali current during the southwest monsoon closing a clockwise path of water mass transport only a small fraction of low salinity water is advected into the eastern arabian sea from the bay of bengal

climate variability in southern south america associated with el nino and la nina events a comprehensive view is given of the precipitation and circulation anomalies associated with the various stages of el nino en and la nina ln events all over southern south america ssa this view comprises the delineation of coherent regions with respect to precipitation anomalies the identification of the seasons of maximum anomalies the indication of their magnitude and the assessment of their consistency during those events in addition the spatial and temporal variability of these anomalies is detailed by calculating the expected precipitation percentiles and the consistency of wet and dry anomalies for each station and each three month running season during en and ln events composites of circulation anomalies and an assessment of their consistency are also presented and their connection with the precipitation anomalies is discussed southern brazil presents the strongest average signal in en events the general behavior toward opposite signals in the precipitation and circulation anomalies over ssa during almost the same periods of the en and ln events indicates a large degree of linearity in the response to these events the timing of the anomalies changes throughout ssa leading to the identification of eight different coherent regions in the en case and six in the ln case this regionalization is mostly caused by different processes leading to precipitation anomalies in ssa during those events all these regions show a significant response in some part of each event the magnitude and consistency of this response show a large spatial variability and some areas present very strong and consistent anomalies sometimes not disclosed when large coherent regions are analyzed in spite of the differences in timing some features of the precipitation anomalies are rather uniform throughout the region during en and ln events in en episode there is a tendency to lower than median precipitation in the year before the event which continues until march of the year of the event in a vast region east of the andes the strongest positive precipitation anomalies occur in spring of this year when the circulation anomalies concur to enhance rainfall over several regions during the summer of the mature stage the positive precipitation anomalies almost disappear and then reappear in some regions in late summer early autumn and in winter of the year following the starting year of the event this description holds partially for the ln event but with opposite signs although there is a larger spatial variability in the ln related anomalies in the following year and some shifts in timing as for precipitation the symmetry of the geopotential height anomaly fields with opposite signs between ln and en cases is also remarkable especially during the year 0

bifurcation of the east india coastal current east of sri lanka the east india coastal current eicc flows equatorward during october december carrying low salinity water from the bay of bengal en route using results from a high resolution ocean general circulation model satellite altimeter data argo float profiles and ocean color images we show that the eicc bifurcates east of sri lanka one part continues along the coast of sri lanka but the major part of the eicc called here as the east sri lanka jet eslj flows eastward into the bay of bengal as a result of this bifurcation there is offshore transport of chlorophyll a rich low salinity water from the coast of sri lanka altimeter data from 1993 2004 show that the bifurcation occurred every year except during the indian ocean dipole iod years of 1994 and 1997 the bifurcation occurs when an anticyclonic eddy that propagates westward ahead of a downwelling rossby wave front impinges on the sri lanka coast this new finding suggests that the main route of the low salinity water from the bay of bengal into the southeastern arabian sea may not be along the coast around sri lanka but through the winter monsoon current

a numerical investigation of choked flow dynamics and its application to the triggering mechanism of long period events at redoubt volcano alaska we use numerical simulations of transonic flow through a crack to study the dynamics of the formation of shock waves downstream from a nozzle like constriction inside the crack the model solves the full set of navier stokes equations in two dimensions via an explicit multifield finite difference representation the crack walls are assumed to be perfectly rigid and elastic coupling to the solid is not considered the simulations demonstrate how the behavior of unsteady shock waves near the walls can produce recurring step like pressure transients in the flow which in turn induce resonance of the fluid filled crack the motion of the shock waves is governed primarily by smooth low amplitude pressure fluctuations at the outlet of the crack the force induced on the walls scales with the amplitude of the shock which is a function of the magnitude of the inlet pressure aperture of the constriction and thickness of the boundary layer the applied force also scales in proportion to the spatial extent of the shock excursion which depends on the fluctuation rate of outlet pressure using the source parameters of long period lp events at redoubt volcano alaska as a guide for our simulations we infer that coupling of the shock to the walls occurs for crack inlet to outlet pressure ratios pi po 2 31 and that the position of the shock front becomes most sensitive to outlet pressure fluctuations for flow regimes with pi po 2 48 for such regimes fluctuations of outlet pressure of up to 0 5 mpa at rates up to 3 mpa s are sufficient to induce pressure transients with magnitudes up to 12 5 mpa over 0 1 2 5 m of the walls within 0 5 s these flow parameters may be adequate for triggering the lp events in the precursory swarm to the december 14 1989 eruption of redoubt according to the flow model the recurrence rate and amplitudes of l p events are inferred to be a manifestation of the response of a shallow hydrothermal reservoir to the sustained injection of superheated steam from a magma column roofing below this reservoir

pressure recovery in magma due to bubble growth 1 pressure recovery in magma subject to sudden depressurization is investigated in terms of the growth of bubbles containing magmatic volatiles assuming conservation of volatile mass and pressure equilibration between the bubbles melt and a surrounding elastic medium the final magma pressure is completely determined by its initial pressure and the magnitude of the pressure drop simulations show that the initial magma pressure is easily recovered or even exceeded when magma containing tiny bubbles is surrounded by a relatively stiff elastic medium under low confining pressure this result suggests that when processes such as magma withdrawal dike intrusion and nearby seismicity decrease the pressure in a magma chamber magma re pressurization can occur without new injection of magma

anomalies in the atmospheric circulation associated with the rainfall excess or deficit in the pampa region in argentina 1 we have investigated the relationship between extreme anomalies in the monthly precipitation of the pampa region pr argentina and large scale circulation anomalies over the south american continent and the neighboring oceans the composite maps of the extreme rainfall anomalies reveal that this atypical behavior extends beyond the pr and that comparable anomalies simultaneously take place in the center and south of chile and in the western south atlantic ocean we have analyzed the anomalies in the water vapor flux and in the circulation of the lower and upper levels of the troposphere that occurred in 85 months with extreme rainfall anomalies between 1955 and 1998 the study combines the use of a complete set of observed precipitation data and several atmospheric variables from the national centers for environmental prediction reanalysis we have found well defined anomalous circulation patterns that embrace the south of the continent and the vicinity of the pacific and atlantic oceans which are very similar but of opposite signs for both rainfall extremes the interaction between these anomalies in the high level flow and the subtropical jet stream generates the necessary conditions for vorticity advection and flow divergence which is consistent with the expected vertical motions the principal component analysis was applied to the pressure anomaly fields to distinguish those aspects of the large scale patterns that are more directly related with the rainfall extremes the amplitudes of some of the dominant pc are significantly correlated with the precipitation anomalies and they highlight the mechanisms leading to rainfall extremes during el nino southern oscillation and neutral periods

characteristics and source modeling of broadband seismic signals associated with the hydrothermal system at satsuma iwojima volcano japan abstract the broadband seismic data observed at satsuma iwojima volcanic island in southwestern japan show a wide variety of interesting features one of them is the periodic and frequency dependent amplitude modulation of volcanic tremors another interesting feature is the occurrence of very long period vlp seismic pulses synchronized with the tremor amplitude modulation both the tremor amplitude modulation and the vlp pulses have fairly regular intervals of 46 50 min suggesting a spatially fixed non destructive repetitive source we performed waveform inversions to reveal the source mechanism of the vlp pulses assuming a point source embedded in a medium that realistically represented the topography of the island numerical tests showed that the green s functions used in the inversion were significantly affected by the topography the observed vlp pulses were well explained by the rapid expansion of an inclined crack 100 m beneath the floor of the summit crater which had a seismic moment of 3 10 10 nm both the characteristic amplitude modulation and the vlp pulses are explained systematically by a model consisting of a water filled crack and a network of gas paths the water filled crack was heated above the boiling temperature of the water in it a sudden break of the unstable superheated state caused sudden vaporization and a pressure increase which was observed as the vlp pulse the network of gas paths consisted of swollen portions and constricted portions which worked as tremor generators when vapor from the water filled crack passed through

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

geodetic observations and modeling of time varying deformation at taal volcano philippines abstract we investigate active crustal deformation patterns associated with magmatic activity at taal volcano an active tholeiitic volcano located in southwestern luzon philippines we present comparisons of elastic and combined elastic viscoelastic modeling results with deformation observed by continuous gps measurements on and near the volcano continuous dual frequency and single frequency l1 gps data between 1998 and 2005 provide evidence for smoothly transitioning periods of inflation and deflation centered under the volcano within deformation phases that typically last 3 9 months with rates exceeding 200 mm yr 1 prominent inflationary phases were observed in february november 2000 and june 2004 march 2005 the 2000 period of inflation was characterized by up to 145 mm yr 1 surface extension and 220 mm yr 1 of uplift of the center of volcano island relative to the northern caldera rim while the 2005 inflation was characterized by as much as 116 mm yr 1 horizontal extension of the volcanic edifice and 180 mm yr 1 uplift inversions of observed surface deformation to determine the source location points to a relatively stable spherical source situated under volcano island s central crater at depths of 5 km based on the preferred mogi model using this source geometry we develop axisymmetric finite element models with annuli of concentric shells around the magma reservoir embedded within a multi layered elastic lithosphere using simply varying pressure histories as input a variety of forward models are fit to the time history of continuously observed deformation for gps stations located on volcano island though the inferred source geometry volume and near source rheological properties are non unique they represent more physically appropriate material properties than those assumed in simple elastic half space models the sources inferred using viscoelastic shell models suggest simpler pressure histories and reduced overall pressure changes relative to equivalent elastic half space models

source process of long period seismic events at taal volcano philippines vapor transportation and condensation in a shallow hydrothermal fissure 1 we analyzed observations of a swarm of more than 40 000 long period lp seismic events at taal volcano philippines in 2010 2011 the event waveforms are strongly correlated to each other consistent with a fixed source location and begin with a dilatational first motion they have a peak frequency around 0 8 hz and a quality factor q of 6 waveform inversion of the events pointed to a tensile crack source dipping 30 60 at a shallow 100 200 m depth a simulation using a fluid filled crack model indicated that the complex frequencies of the waveforms are explained by the fundamental longitudinal mode resonance of a vapor filled crack 188 m long a satellite thermal infrared image acquired during the swarm period suggests that the lp events were not accompanied by surface gas releases we considered a vapor transportation model in which vapor exsolved from magma and rose in a fissure extending to the lp source this model yielded estimates that 105 107 m3 of magma was involved in the lp swarm and that the temperature of vapor in the lp source crack was around 600 k we modeled a triggering mechanism of the crack resonance based on sudden condensation of vapor at the crack tip in a cold aquifer this model explained observed characteristics of the events including the dilatational first motion the total volumetric change and the fixed source location

historical climatic extremes as indicators for typical scenarios of holocene climatic periods in the pampean plain abstract the holocene has been marked by two different climatic periods in the pampas one 8 5 3 5 kyr b p characterized by temperatures and precipitation higher than today and the other 3 5 1 4 kyr b p semiarid generating parabolic dunes the humid period produced intense pedogenesis down to 40 lat s and a mobilization of iron oxides to 30 lat s which means tropical climate with temperatures above 20 c and precipitation higher than 2000 mm yr sand deflation and development of parabolic dunes occur under climates with 300 400 mm yr this contribution is based on geological and physical proxies typical weather scenarios of those climates are sporadically reproduced today during extremely humid or dry periods on that basis meteorological parameters of years beyond the thresholds of 2000 and 400 mm yr were processed dry periods were characterized by a large thermal amplitude frosts and stronger winds reproducing the continental anticyclonic circulation of the late holocene humid extremes were warmer than normal lower thermal amplitude with rains produced by local convection processes typical scenarios with precipitation above 2000 mm yr are 21 c mean temperature ca 1 c higher than the present record 27 c maximal annual temperatures 16 c minimal annual temperatures more than 1 c higher than the present as a consequence thermal amplitude was smaller than today virtually without frosts characteristic parameters of the late holocene were mean annual precipitation ca 350 mm yr 15 c mean temperature 22 c maximal annual temperatures 8 c minimal annual temperatures showing significative shifts in the dry season lower than normal monthly maximal and mean temperatures higher than today the result is a thermal amplitude larger than today probably higher quantity of frosts per year stronger winds and lower relative air humidity in comparison with the humid climate extreme complete the scenario

origins of eddy kinetic energy in the bay of bengal by analyzing satellite observational data and ocean general circulation model experiments this study investigates the key processes that determine the spatial distribution and seasonality of intraseasonal eddy kinetic energy eke within the bay of bengal bob it is revealed that a complicated mechanism involving both local and remote wind forcing and ocean internal instability is responsible for the generation and modulation of eke in this region high level eke mainly resides in four regions east of sri lanka region 1 the western bob region 2 northwest of sumatra region 3 and the coastal rim of the bob region 4 the high eke levels in regions 1 and 2 are predominantly produced by ocean internal instability which contributes 90 and 79 respectively prominent seasonality is also observed in these two regions with higher eke levels in boreal spring and fall due to enhanced instability of the east indian coast current and the southwest monsoon current respectively in contrast ocean internal instability contributes 49 and 52 of the total eke in regions 3 and 4 respectively whereas the atmospheric forcing of intraseasonal oscillations isos also plays an important role isos produce eke mainly through wind stress involving both the remote effect of equatorial winds and the local effect of monsoonal winds equatorial origin wave signals significantly enhance the eke levels in regions 3 and 4 in the form of reflected rossby waves and coastal kelvin waves respectively the local wind forcing effect through ekman pumping also has a significant contribution in regions 3 and 4 24 and 22 respectively

fluid dynamics inside a wet volcano inferred from the complex frequencies of long period lp events an example from papandayan volcano west java indonesia during the 2011 seismic unrest abstract we present results of our study aimed at understanding the dynamics of fluids inside a wet volcano through the analysis of swarms of long period lp events accompanying the 2011 seismic unrest at papandayan volcano west java indonesia prior to this unrest we measured an extremely high percentage 100 of co2 in the ground at the summit crater however with a very low value of so2 flux 6 tons day increase in volcanic activity was also observed from the records of a tiltmeter a long term inflation was followed by an abrupt deflation that took place concurrently with the swarms of lp events thereafter swarms of local tectonic lt and volcano tectonic vt earthquakes started we focus here on analyzing the lp events in the following manner first we estimate the source location of lp events by applying a 3 d non linear hypocenter localization algorithm which includes topography we then study the waveforms and spectral characteristics of lp events recorded at different stations and investigate whether or not these characteristics are due to source effects subsequently we compute the oscillation frequencies f and the decay characteristics q factor in the complex frequency domain of the coda part of the lp events by using the sompi method which is based on a homogeneous auto regressive ar equation the rectangular fluid filled crack model is used to estimate the physical processes related to the observed temporal variations in the complex frequencies we divide the swarms of lp events into two intervals the first interval occurred between june and july 2011 48 lp events while the second interval extended from september to october 2011 36 lp events the frequencies of lp events observed during these intervals range between 1 1 and 6 2 hz while the q factors are widely scattered between 20 and 400 we estimate the compositions of fluids inside the crack during both intervals as either water foam mixtures of water and h2o gas steam or misty gas mixtures of water droplets and h2o gas steam we finally suggest that if an eruption were to have taken place following the 2011 unrest it would have been in phreatic style rather than magmatic style the results of our study therefore contribute to the effort in the prediction of the behavior of future eruptions and to volcanic hazards assessment and therefore to volcanic risk mitigation

electrical resistivity tomography study of taal volcano hydrothermal system philippines taal volcano 311 m in altitude is located in the philippines 14 n 121 e and since 1572 has erupted 33 times causing more than 2 000 casualties during the most violent eruptions in march 2010 the shallow structures in areas where present day surface activity takes place were investigated by dc resistivity surveys electrical resistivity tomography ert lines were performed above the two identified hydrothermal areas located on the northern flank of the volcano and in the main crater respectively due to rough topography deep valleys and dense vegetation most measurements were collected using a remote method based on a laboratory made equipment this allowed retrieval of information down to a depth of 250 m erts results detail the outlines of the two geothermal fields defined by previous self potential co2 soil degassing ground temperature and magnetic mapping harada et al japan acad sci 81 261 266 2005 zlotnicki et al bull volcanol 71 29 49 2009a phys chem earth 34 294 408 2009b hydrothermal fluids originate mainly from inside the northern part of the main crater at a depth greater than the bottom of the crater lake and flow upward to the ground surface furthermore water from the main crater lake infiltrates inside the surrounding geological formations the hydrothermal fluids outlined by gas releases and high temperatures cross the crater rim and interact with the northern geothermal field located outside the main crater

long period seismic source characterization at popocatepetl volcano mexico 1 the seismicity of popocatepetl is dominated by long period and very long period signals associated with hydrothermal processes and magmatic degassing we model the source mechanism of repetitive long period signals in the 0 4 2 s band from a 15 station broadband network by stacking long period events with similar waveforms to improve the signal to noise ratio the data are well fitted by a point source located within the summit crater 250 m below the crater floor and 200 m from the inferred magma conduit the inferred source includes a volumetric component that can be modeled as resonance of a horizontal steam filled crack and a vertical single force component the long period events are thought to be related to the interaction between the magmatic system and a perched hydrothermal system repetitive injection of fluid into the horizontal fracture and subsequent sudden discharge when a critical pressure threshold is met provides a non destructive source process

temporal evolution of a magmatic dike system inferred from the complex frequencies of very long period seismic signals 1 i present a detailed description of temporal variations in the complex frequencies of very long period vlp events observed in association with an earthquake swarm that occurred beneath hachijo island japan the individual vlp signals show simple decaying harmonic oscillations with periods near 10 s and they last for about 300 s the vlp activity began on 20 august 2002 and continued until 2 april 2003 spectral analysis of the vlp waveforms identifies two spectral peaks of which the complex frequencies show systematic temporal variations at the beginning of the vlp activity the q values of both the peaks show increasing trends from roughly 15 to 35 in the time period between 22 and 26 august the frequency of one peak shows a decreasing trend from 0 1 to 0 09 hz while the frequency of the other peak shows an increasing trend from 0 128 to 0 136 hz during this time after these systematic variations the frequencies and q values of the two peaks remain roughly constant waveform simulations based on a fluid filled crack model indicate that the observed systematic temporal variations can be reasonably explained by gradual expansion of a crack containing a basalt gas mixture the roughly constant frequencies and q values of the two peaks after these systematic variations suggest that the magmatic dike persistently existed beneath hachijo island for more than 6 months mass and heat advection is required to sustain the dike system implying magmatic flow and convection in the dike and magma chamber beneath this island

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

global co2 emission from volcanic lakes the global co 2 discharge from subaerial volcanism has been estimated at 300 mt yr 1 however estimates of co 2 emissions from volcanic lakes have not been considered in order to improve this information extensive research on co 2 emissions of volcanic lakes worldwide has been performed the observed normalized average co 2 emission rates increase from alkaline 5 5 t km 2 d 1 to neutral 201 2 t km 2 d 1 to acid 614 2 t km 2 d 1 in volcanic lakes taking into account 1 normalized co 2 emission rates 2 the number of volcanic lakes in the world 769 and 3 the fraction and average areas of the investigated alkaline neutral and acid volcanic lakes the estimated global co 2 emission from volcanic lakes is 117 19 mt yr 1 with 94 17 mt yr 1 as deep seated co 2 this study highlights the importance of a revision of the actual global co 2 discharge from subaerial volcanism

long period seismicity and very long period infrasound driven by shallow magmatic degassing at mount pagan mariana islands long period lp seismicity and very long period infrasound ivlp were recorded during continuous degassing from mount pagan mariana islands in july 2013 to january 2014 the frequency content of the lp and ivlp events and delay times between the two arrivals were remarkably stable and indicate nearly co located sources using phase weighted stacking over similar events to dampen noise we find that the lp source centroid is located 60 m below and 180 m west of the summit vent the moment tensor reveals a volumetric source modeled as resonance of a subhorizontal sill intersecting a dike we model the seismoacoustic wavefields with a coupled earth air 3 d finite difference code the ratios of pressure to velocity measured at the infrasound arrays are an order of magnitude larger than the synthetic ratios so the ivlp is not the result of lp energy transmitting into the atmosphere at its epicenter based on crater shape and dimensions determined by structure from motion we model the ivlp as acoustic resonance of an exponential horn the source of the continuous plume from gas analysis is shallow magmatic degassing which repeatedly pressurized the dike sill portion of the conduit over the 7 months of observation periodic gas release caused the geologically controlled sill to partially collapse and resonate while venting of gas at the surface triggered resonance in the crater lp degassing only accounts for 12 of total degassing indicating that most degassing is relatively aseismic and that multiple active pathways exist beneath the vent

analysis of crustal deformation in luzon philippines using geodetic observations and earthquake focal mechanisms we utilize regional gps velocities from luzon philippines with focal mechanism data from the harvard centroid moment tensor cmt catalog to constrain tectonic deformation in the complex plate boundary zone between the philippine sea plate and eurasia the sundaland block processed satellite imagery and digital elevation models are used with existing gravity anomaly seismicity and geologic maps to define a suite of six elastic blocks geodetic and focal mechanism data are inverted simultaneously to estimate plate rotations and fault locking parameters for each of the tectonic blocks and faults comprising luzon major tectonic structures that were found to absorb the plate convergence include the manila trench 20 100 mm yr 1 and east luzon trough 9 15 mm yr 1 philippine trench 29 34 mm yr 1 which accommodate eastward and westward subduction beneath luzon respectively the left lateral strike slip philippine fault 20 40 mm yr 1 and its northward extensions the northern cordillera fault 17 37 mm yr 1 transtension and the digdig fault 17 27 mm yr 1 transpression the macolod corridor a zone of active volcanism crustal thinning extension and extensive normal and strike slip faulting in southwestern luzon is associated with left lateral transtensional slip of 5 10 mm yr 1 the marikina fault which separates the central luzon block from the southwestern luzon block reveals 10 12 mm yr 1 of left lateral transpression our analysis suggests that much of the philippine fault and associated splays are locked to partly coupled while the manila and philippine trenches appear to be poorly coupled luzon is best characterized as a tectonically active plate boundary zone comprising six mobile elastic tectonic blocks between two active subduction zones the philippine fault and associated intra arc faults accommodate much of the trenchparallel component of relative plate motion 2006 elsevier b v all rights reserved

a 2 ma record of explosive volcanism in southwestern luzon implications for the timing of subducted slab steepening 1 new chemical and isotopic analyses of the tephra layers plus deep sea tephrostratigraphic record from two cores from either side of luzon island philippines have allowed the identification of two periods of explosive volcanic activity originating from the macolod corridor in the southwestern part of the luzon the first period extended from prior to 1355 ka to 1977 ka and the second period extended from 478 ka to the present separated by a period of relative quiescence the time intervals between large explosive eruption events in each period were 31 15 ka and 156 52 ka respectively combined with published chronological and geochemical data from onshore volcanic deposits the tephrostratigraphic record shows that the locus of large explosive eruptions has migrated southwestward from the northeastern section to the middle and southwestern sections of the macolod corridor the period of relative quiescence is characterized by monogenetic volcanism in the central section of the corridor the migration of active volcanism across the southwestern part of luzon during the quaternary is used to infer the evolution of the subducting south china sea crust the period of relative quiescence represents a period of adjustment of the subducted slab by steepening which began around 1355 ka or shortly thereafter and finished at around 478 ka

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

laboratory simulation of volcano seismicity the physical processes generating seismicity within volcanic edifices are highly complex and not fully understood we report results from a laboratory experiment in which basalt from mount etna volcano italy was deformed and fractured the experiment was monitored with an array of transducers around the sample to permit full waveform capture location and analysis of microseismic events rapid post failure decompression of the water filled pore volume and damage zone triggered many low frequency events analogous to volcanic long period seismicity the low frequencies were associated with pore fluid decompression and were located in the damage zone in the fractured sample these events exhibited a weak component of shear double couple slip consistent with fluid driven events occurring beneath active volcanoes

on the functioning of the southern oscillation in the south american sector part i surface climate abstract the functioning of the southern oscillation so over south america and the surrounding tropical atlantic and eastern pacific oceans is studied through analyses of bimonthly patterns of correlation between an so index and the fields of pressure wind temperature and rainfall hydrometeorological anomalies being further ascertained from river discharge and lake level records and measurements of outgoing longwave radiation most pronounced is the anomalously low pressure over the eastern pacific during the negative so phase defined by anomalously low high pressure at tahiti darwin consistently positive so pressure correlations over the eastern tropical pacific extend eastward into the north atlantic region and southern south america during the respective winter semester during the austral winter the tendency for anomalously abundant rainfall in central chile during the negative so phase appears related to an anomalously weak and northward displaced south pacific subtropical high while relativ