fractal properties of tremor and gas piston events observed at kilauea volcano hawaii we study the fractal properties of shallow volcanic tremor and gas piston events associated with magma degassing at kilauea volcano hawaii using data from two dense short baseline arrays of seismographs deployed near the active crater of puu oo on the east rift of the volcano we found an upper bound on the fractal dimension of a strange attractor common to phase portraits of both tremor and gas piston activities in the range of 3 1 4 5 with a mean of 3 75 over a scale length of the wave field of the order of 100 m the dimension of the attractor fluctuates over the course of tremor and gas piston episodes within a range comparable to the spatial variations of fractal dimensions found in the wave field the existence of a categorically stable attractor characterizing both types of activities strongly suggests that the excitation mechanism of tremor is similar to that of gas piston events which in turn are correlated with visual observations at the volcanic vent the low value derived for the dimension of the attractor in phase space points to significant self organization in the process of generation of tremor and offers general constraints on the dimensionality of attractors derived from models of acoustic emission associated with magma flow vesiculation and degassing

spectral and dynamical hints on the time scale of preparation of the 5 april 2003 explosion at stromboli volcano stromboli volcano is well known for its continuous strombolian activity moreover the volcano occasionally shows effusive phases the latest in 1985 1986 on 28 december 2002 stromboli entered a new effusive phase accompanied by different paroxysmal events that led to considerable hazards for inhabitants and tourists on the island of stromboli on 30 december 2002 a major sector collapse affected the sciara del fuoco slope and initiated a tsunami on 5 april 2003 a powerful explosion which can be compared in size with the most recent explosion in 1930 covered a large part of the normally tourist accessible summit area with bombs as this explosion was not forecasted although the island was by then effectively monitored by a dense deployment of instruments in this paper we tackle the problem of highlighting the time scale of preparation of this event and conduct a search for possible precursors for this purpose we analyze the seismic data preceding the paroxysm with spectral and dynamical methods

strombolian explosions 1 a large bubble breaking at the surface of a lava column as a source of sound strombolian activity consists of a series of explosions caused by the breaking of large overpressurized bubbles at the surface of the magma column acoustic pressure has been measured for 36 explosions at stromboli we propose that sound is generated by the vibration of the bubble before it bursts oscillations are driven by an initial overpressure inside the bubble assumed to be initially at rest just below the magma air interface inertia effects cause the bubble to overshoot its equilibrium radius then the bubble becomes underpressurized and contracts because of gas compressibility these oscillations are only slightly damped by viscous effects in the magma layer above the bubble the bubble cannot complete more than one cycle of vibration because of instabilities developing on the magma layer that lead to its breaking near the minimum radius assuming a simple geometry we model this vibration and constrain the radius and length of the bubble and the initial overpressure by fitting a synthetic waveform to the measured acoustic pressure the fit between synthetic and observed waveforms is very good both for frequency 60 rad s 1 and amplitude the initial bubble radius is 1 m and the length varies between several and a few tens of meters from the initial overpressure approximately 105 pa we calculate the maximum radial velocity of ejecta 30 m s 1 the generally good agreement between data and predictions of our model leads us to suggest that acoustic measurements are a powerful tool for the understanding of eruption dynamics

infrasonic waves and volcanic tremor at stromboli the origin of the volcanic tremor is still under debate many theories have been proposed in the last years but none has yet been completely accepted in 1993 highly sensitive pressure sensors 2 175 pa volt used to monitor the explosive activity at stromboli have revealed unexpected correlation between small spike shaped pressure signals 1 2 pa and volcanic tremor these pressure pulses repeat regularly in time with a recurrent period of ca 1 s video camera images allowed us to correlate the pressure pulses with small gas bursts occurring at one of the active vents the striking correlation 0 971 between infrasonic and seismic energy fluctuations is particularly meaningful in the frequency domain infrasonic and seismic signal share the same spectral content 3 hz for every station within a range of 700 m around the craters correlations in time and frequency domain remained unaltered during the 1994 field experiments moreover during 1994 the increased degassing activity has been followed by an increase in pressure release 7 8 pa and by a shift towards higher frequencies 8 hz both in the infrasonic and seismic records infrasonic waves and volcanic tremor show similar energy fluctuations and frequency contents appearing therefore to be produced by the same dynamical process on this basis we claim that volcanic tremor at stromboli originates by continuous outbursting of small gas bubbles in the upper part of the magmatic column

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

a m gurnell and m j clark eds 1987 glacio fluvial sediment transfer an alpine perspective chichester etc john wiley and sons ltd

analysis of the temporal occurrence of seismicity at deception island antarctica a nonlinear approach deception island is characterized by small magnitude local events with constant energy flux and very low stress drop to obtain information about its origin an interevent time series of 546 events corresponding to an observational period of two month has been analyzed from a statistical point of view data satisfies a weibull distribution and presents clustering a rescaled range analysis reveals that data are not independent i e have memory and the correlation dimension saturates at 2 2 as a consequence the system can be modeled as a nonlinear iterative equation with three degrees of freedom that presents chaotic behavior taking into account that the average interevent time is of the order of 130 minutes too short to be only due to tectonic activity the above results indicate that some other mechanism may coexist with the regional tectonic one according to several geological and geophysical observations we suggest that most of the local events may be originated by pressure waves generated by a sudden change of phase of sea and fresh water infiltrated into the main fractures and faults and also from shallow and confined water saturated layers

automated mapping of glacial overdeepenings beneath contemporary ice sheets approaches and potential applications awareness is growing on the significance of overdeepenings in ice sheet systems however a complete understanding of overdeepening formation is lacking meaning observations of overdeepening location and morphometry are urgently required to motivate process understanding subject to the development of appropriate mapping approaches high resolution subglacial topography data sets covering the whole of antarctica and greenland offer significant potential to acquire such observations and to relate overdeepening characteristics to ice sheet parameters we explore a possible method for mapping overdeepenings beneath the antarctic and greenland ice sheets and illustrate a potential application of this approach by testing a possible relationship between overdeepening elongation ratio and ice sheet flow velocity we find that hydrological and terrain filtering approaches are unsuited to mapping overdeepenings and develop a novel rule based gis methodology that delineates overdeepening perimeters by analysis of closed contour properties we then develop gis procedures that provide information on overdeepening morphology and topographic context limitations in the accuracy and resolution of bed topography data sets mean that application to glaciological problems requires consideration of quality control criteria to a remove potentially spurious depressions and b reduce uncertainties that arise from the inclusion of depressions of nonglacial origin or those in regions where empirical data are sparse to address the problem of overdeepening elongation potential quality control criteria are introduced and discussion of this example serves to highlight the limitations that mapping approaches and applications of such approaches must confront we predict that improvements in bed data quality will reduce the need for quality control procedures and facilitate increasingly robust insights from empirical data

evidence for gas influence on volcanic seismic signals recorded at stromboli abstract the meaning of the large variety of seismic waveforms showing different frequencies associated with volcanic activity is still uncertain no definitive model for source dynamics has been proposed yet at present seismic models explaining the origin of the spectral content of the seismic signal are mainly linked to geometrical features of the volcanic conduit and to resonant effects induced by pressure fluctuations in a fluid filled conduit such models assume the physical system to be in a steady state over a long interval of time at stromboli the seismicity produced by each of the three active vents is not stable and can suddenly change in time therefore the application of stable resonator models appears to be unwarranted on the basis of infrared image analyses atmospheric pressure and seismic wavefield signals we show that different frequency contents and different explosive styles occur at stromboli at the same time at different vents moreover we give evidence to support the idea that seismicity and explosivity are both controlled by a variable gas flow regime released during explosions at the top of the magma column we have recognized two main pressure regimes controlling the explositivity at stromboli one is impulsive short and possibly produced by a gas accumulation beneath a cap rock the other has an intermittent regime feature and lasts longer we demonstrate how the first one generates a monochromatic low frequency wavefield while the second shows a high frequency spectrum where the frequency content depends on pressure fluctuations and on pressure gradient we suggest different pressure growth and gas flow in the magma as the common source for both explosive style and seismic wavefield

on the spectral peaks of volcanic tremor at stromboli abstract seismic activity at stromboli volcano is characterized by a variety of signals emanating from three vents for a long time the northwest vent has been in constant activity periodically large explosions occur and material is ejected beyond the crater walls these large explosions are accompanied by sonic and infrasonic pressure waves in the atmosphere and explosion quakes apart from large explosions there is constant activity in the form of continuous gas bursts which are related to low infrasonic pulses in the atmosphere and volcanic tremor we assume that volcanic tremor and low pressure infrasonics are generated by gas bubbles inside the volcanic conduit and accordingly we compute synthetic tremor by modeling the source function as a pressure variation in a spherical cavity that propagates through a finely layered medium by means of haskell s formalism to simulate a tremor we superpose in time domain a large number of such pulses of varying amplitudes and time delays according to the observed infrasonic series in addition to the spectral similarity the observed and synthetic tremor display the same autocorrelation and hurst exponents implying similar long term correlation we present strong evidence in favour of an interpretation of the spectral peaks of the volcanic tremor at stromboli in terms of resonances of the layered structure hence as a path effect rather than a source effect

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

the origin of harmonic tremor at old faithful geyser volcanic tremor the seismic signal which is often associated with eruptions is among the least understood phenomena in seismology large variations in the physical properties of volcanic rocks and fluids along with the complex geological structure make it difficult to obtain reliable source models old faithful geyser yellowstone national park was chosen as a laboratory for studying fluid flow induced seismicity the geyser s cyclic behavior on an hourly time scale and its accessibility make it suitable for a detailed study of seismic behavior between eruptions simultaneously with underwater pressure measurements inside the geyser we observe that sharp pressure pulses inside the water column attributed to steam bubble collapse are followed by distinct seismic events with a sharp onset and a harmonic horizontal motion whose frequency varies spatially but not temporally a superposition of these events creates the appearance of continuous harmonic tremor since the pressure in the water column exhibits no sign of resonance the harmonic motion must be caused by elastic waves reverberating in the solid medium a near surface soft layer is demonstrated to be a possible elastic model the harmonic tremor observed at old faithful seems analogous to some examples of shallow volcanic tremor

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

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

interactions between bubbles in magmas and lavas effects of bubble deformation the interactions between two deformable bubbles are studied in order to determine the effects of deformation on bubble dynamics and to determine the limits in which the effects of deformation can be ignored deformation due to hydrodynamic interactions with other bubbles leads to alignment of horizontally offset bubbles and thus an enhanced rate of coalescence bubble alignment may produce spatial inhomogeneities in bubble concentrations in magmas which implies temporal variations in the concentration and size of bubbles in erupting magmas and spatial variations in solidified lavas an approximate quantitative model for the rate of coalescence of deformable bubbles is developed based on a series of experiments the effects of deformation can be ignored if the largest bubbles have radii less than about 5 mm in silicate magmas

rapid oceanic and atmospheric changes during the younger dryas cold period the younger dryas event which began approximately 12 900 years ago was a period of rapid cooling in the northern hemisphere driven by large scale reorganizations of patterns of atmospheric and oceanic circulation 1 3 environmental changes during this period have been documented by both proxy based reconstructions 3 and model simulations 4 but there is currently no consensus on the exact mechanisms of onset stabilization or termination of the younger dryas 5 8 here we present high resolution records from two sediment cores obtained from lake krakenes in western norway and the nordic seas multiple proxies from lake krakenes are indicative of rapid alternations between glacial growth and melting during the later younger dryas meanwhile reconstructed sea surface temperature and salinity from the nordic seas show an alternation between sea ice cover and the influx of warm salty north atlantic waters we suggest that the influx of warm water enabled the westerly wind systems to drift northward closer to their present day positions the winds thus brought relatively warm maritime air to northern europe resulting in rising temperatures and the melting of glaciers subsequent input of this fresh meltwater into the ocean spurred the formation of sea ice which forced the westerly winds back to the south cooling northern europe we conclude that rapid alternations between these two states immediately preceded the termination of the younger dryas and the permanent transition to an interglacial state the impact of freshwater pulses released from northern hemisphere ice sheets on the north atlantic meridional overturning circulation has been suggested as a likely mechanism for the younger dryas cooling 7 11 in addition a responsive sea ice cover has been proposed as a mechanism that could rapidly alter the heat exchange between the ocean and atmosphere during the last glacial period 12 13 the roles of the ocean and the atmosphere in triggering the climate shifts around the younger dryas are much discussed 9 these hypotheses are difficult to resolve partly because the temporal resolution of most palaeoclimate archives is too coarse to enable the anatomy of the younger dryas stadial to be adequately described 14 to address these questions we here present a lake sediment record at sub annual resolution together with new data on the sea surface conditions of the norwegian atlantic current nwac throughouttheyoungerdryasstadial fig 1

modeling of volcanic tremor as explosive point sources in a single layered elastic half space i investigate path effects on the waveform of a volcanic tremor using a model that consist of an impulsive dilatational point source embedded in a single layered elastic half space such a source may arise during a variety of volcanic conditions such as a sequence of rapid explosions at a volcanic vent as well as gas piston activity observed in basaltic lava lakes using acoustic signals recorded during the 1975 tolbachik eruption the 1983 klyuchevskoy eruption and the 1989 tokachi dake eruption i propose a general impulsive time function for the tremor source consisting either of a single impulse or a sequence of identical impulses computed seismograms illustrate the influence of a low velocity surface layer on the shape and spectra of the waveforms for the same layered model a single impulse produces a dispersive wave while a sequence of identical impulses produces a wave of multiple harmonics for a single impulse the peak in the seismograms shifts to higher frequencies when the surface layer is very thin for multiple impulses spectral peaks shift to higher frequencies as the time between successive impulses increases indicating an increasing importance of higher order harmonics

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

interaction of seismic and air waves recorded at stromboli volcano explosion quake seismograms recorded at stromboli show that seismic phases with a high amplitude and high frequency content propagate with a velocity of approximately 330 m s the sound speed the analysis of seismograms recorded at a distance of 500 m from one of the three active vents shows for the first onset a low frequency and particle motion characteristic of a p wave which loses its longitudinal polarization with the onset of the air wave recording the explosion quakes simultaneously with a microphone we could ascertain that the high frequency onset coincided with the air wave s in order to better understand the seismic wavefield generated by the atmospheric pressure we performed a controlled source experiment at stromboli using a seismic gun seismograms with the same two phases and particle motions comparable with the volcano seismic data were obtained a second experiment demonstrated that the air wave propagates at least in the uppermost 1 m of the ground we suggest that the seismic source of the corresponding seismograms is an explosion at the top of the magma column and conclude that the p and air waves are both generated in the same point and at the same time

formalising an inversion methodology for reconstructing ice sheet retreat patterns from meltwater channels application to the british ice sheet glacial meltwater channels yield important information regarding patterns of palaeo ice sheet retreat however it is important to distinguish between channels formed at different positions in relation to the ice mass e g lateral subglacial in britain whilst meltwater channels have been widely reported in the literature such categorisation is often lacking here we develop an inversion methodology for extracting ice margin retreat patterns from meltwater channels diagnostic geomorphological criteria are presented derived from published descriptions of meltwater channel form by which a channel of unknown origin may be assigned to its appropriate category these criteria are applied to channels recorded in the britice gis database of glacial landforms and we find most channels to be of lateral or subglacial origin the categorisation of channels enables palaeo margin positions and the direction of retreat stage ice flow to be approximated and a preliminary reconstruction is made of fragments of the ice sheet retreat pattern further mapping and integration of melt traces with other geomorphic data will permit a more complete palaeo glaciological reconstruction of the last british ice sheet copyright 2007 john wiley sons ltd

continuous measurement of sediment transport in the erlenbach stream using piezoelectric bedload impact sensors we report on bedload transport observations using piezoelectric bedload impact sensors pbis an indirect method of estimating the volume of bedload transport of coarse sediment the pbis device registers vibrations produced by bedload particle diameter 20 mm and records the signal as a sum of the number of impulses per time sediment transport at the erlenbach stream has been continuously monitored with a pbis array starting in 1986 the sensor array spans the width of an entire cross section and is mounted flush with the surface of a check dam immediately upstream of a sediment retention basin we compare pbis data with long term sedimentation records obtained from repeated surveys of material stored in the sediment retention basin with artificial sediment input under controlled conditions in the field and also with laboratory experiments the rate of bedload transport is proportional to the number of impacts on the sensor per unit time the reliability of the calibration relationship increases with the length of the observation period e g for higher numbers of impacts and larger bedload volumes sediment volumes for individual flood events estimated with the pbis method are in agreement with volumes estimated using an independent empirical method based on the effective runoff volume of water the peak water discharge and the critical discharge for the onset of sediment transport copyright 2007 john wiley sons ltd

sedimentation in arctic proglacial lakes mittivakkat glacier south east greenland several sediment cores were collected from two proglacial lakes in the vicinity of mittivakkat glacier southeast greenland in order to determine sedimentation rates estimate sediment yields and identify the dominant sources of the lacustrine sediment the presence of varves in the ice dammed icefall lake enabled sedimentation rates to be estimated using a combination of x ray photography and down core variations in 137 cs activity sedimentation rates for individual cores ranged between 0 52 and 1 06 g cm 2 year 1 and the average sedimentation rate was estimated to be 0 79 g cm 2 year 1 despite considerable down core variability in annual sedimentation rates there is no significant trend over the period 1970 to 1994 after correcting for autochthonous organic matter content and trap eaciency the mean fine grained minerogenic sediment yield from the 3 8k m 2 basin contributing to the lake was estimated to be 327 t km 2 year 1 cores were also collected from the topset beds of two small deltas in icefall lake the deposition of coarse grained sediment on the delta surface was estimated to total in excess of 15 cm over the last c 40 years in the larger lake kuutuaq which is located about 5 km from the glacier front and for which the glacier represents a smaller proportion of the contributing catchment sedimentation rates determined for six cores collected from the centre of the lake based on their 137 cs depth profiles were estimated to range between 0 05 and 0 11 g cm 2 year 1 and the average was 0 08 g cm 2 year 1 the longer term c 100 150 years average sedimentation rate for one of the cores estimated from its unsupported 210 pb profile was 0 10 0 13 g cm 2 year 1 suggesting that sedimentation rates in this lake have been essentially constant over the last c 100 150 years the average fine grained sediment yield from the 32 4k m 2 catchment contributing to the lake was estimated to be 13 t km 2 year 1 the 137 cs depth profiles for cores collected from the topset beds of the delta of lake kuutuaq indicate that in excess of 27 cm of coarse grained sediment had accumulated on the delta surface over the last approximately 40 years caesium 137 concentrations associated with the most recently deposited uppermost fine grained sediment in both icefall lake and lake kuutuaq were similar to those measured in fine grained sediment collected from steep slopes in the immediate proglacial zone suggesting that this material rather than contemporary glacial debris is the most likely source of the sediment deposited in the lakes this finding is confirmed by the 137 cs concentrations associated with suspended sediment collected from the mittivakkat stream which are very similar to those for proglacial material copyright 2000 john wiley sons ltd

kinematic parameters of the tremor wave field at mt etna sicily the tremor wave field at mt etna was recorded with an array of three 3 component broad band seismometers the analysis of its polarisation and dispersion properties leads to the conclusions that 1 the wave field consists mainly of surface waves and that 2 two of the observed spectral peaks are caused by a wave guide effect along the propagation path

analysis and interpretation of patterns within and between hydroclimatological time series in an alpine glacier basin this paper explores patterns within and between climatological and hydrological time series from an alpine glacier basin time series recorded in the basin of the haut glacier d arolla over the 1989 ablation season are subdivided into five subperiods box jenkins arima autoregressive integrated moving average and tfn transfer function noise models are estimated for each of the five subperiods and differences between the models are interpreted in the context of changing glacier hydrology particularly the changing nature and extent of the glacier drainage network

geochemical weathering at the bed of haut glacier d arolla switzerland a new model waters were sampled from 17 boreholes at haut glacier d arolla during the 1993 and 1994 ablation seasons three types of concentrated subglacial water were identified based on the relative proportions of ca 2c h co 3 and so4 2 to si type a waters are the most solute rich and have the lowest relative proportion of si they are believed to form in hydrologically inefficient areas of a distributed drainage system most solute is obtained from coupled sulphide oxidation and carbonate dissolution so cd it is possible that there is a subglacial source of o2 perhaps from gas bubbles released during regelation because the high so4 2 levels found up to 1200 eq l are greater than could be achieved if sulphides are oxidized by oxygen in saturated water at 0 c c 414 eq l a more likely alternative is that sulphide is oxidized by fe 3c in anoxic environments if this is the case exchange reactions involving fe iii and fe ii from silicates are possible these have the potential to generate relatively high concentrations of hco3 with respect to so4 2 formation of secondary weathering products such as clays may explain the low si concentrations of type a waters type b waters were the most frequently sampled subglacial water they are believed to be representative of waters flowing in more efficient parts of a distributed drainage system residence time and reaction kinetics help determine the solute composition of these waters the initial water rock reactions are carbonate and silicate hydrolysis and there is exchange of divalent cations from solution for monovalent cations held on surface exchange sites hydrolysis is followed by so cd the so4 2 concentrations usually are 414 eq l although some range up to 580 eq l which suggests that elements of the distributed drainage system may become anoxic type c waters were the most dilute yet they were very turbid their chemical composition is characterized by low so4 2 h co3 ratios and high ph type c waters were usually artefacts of the borehole chemical weathering environment true type c waters are believed to flow through sulphide poor basal debris particularly in the channel marginal zone the composition of bulk runoff was most similar to diluted type b waters at high discharge and was similar to a mixture of type b and c waters at lower discharge these observations suggest that some supraglacial meltwaters input to the bed are stored temporarily in the channel marginal zone during rising discharge and are released during declining flow little of the subglacial chemical weathering we infer is associated with the sequestration of atmospheric co2 t he progression of reactions is from carbonate and silicate hydrolysis through sulphide oxidation by first oxygen and then fe iii which drives further carbonate and silicate weathering a crude estimate of the ratio of carbonate to silicate weathering following hydrolysis is 4 1 we speculate that microbial oxidation of organic carbon also may occur both sulphide oxidation and microbial oxidation of organic carbon are likely to drive the bed towards suboxic conditions hence we believe that subglacial chemical weathering does not sequester significant quantities of atmospheric co2 and that one of the key controls on the rate and magnitude of solute acquisition is microbial activity which catalyses the reduction of fe iii and the oxidation of fes2 copyright 2002 john wiley sons ltd

an isotropic source of volcanic tremor observation with a dense seismic network at izu oshima volcano japan abstract volcanic tremor is frequently observed associated with activities of volcanoes in order to study the nature of the volcanic tremor in detail we installed a dense seismic network covering the area of the summit crater of izu oshima volcano we succeeded in observing volcanic tremors preceding a small explosive eruption on november 16 1987 the explosion was followed by a significant draining back of lava which had filled the summit crater of the 1986 eruption we analyzed the low frequency components of the tremor near 0 7 hz and found that it was generated by an isotropic volumetric source located 350 m above sea level assuming a cyclic pressure change at the source the amplitude of the seismic moment is 3 4 10 10 nm this corresponds to a pressure fluctuation of 3 10 3 pa 0 03 bar in the system containing the cavity of 10 7 m 3 which is comparable to the volume of subsidence at the 1987 eruption

excitation of a buried magmatic pipe a seismic source model for volcanic tremor recent observations of seismic events at various volcanoes suggest that harmonic tremor results from the sustained occurrence of so called long period or low frequency events accordingly we can view the long period volcanic event as the elementary process of tremor and interpret it as the impulse response of the tremor generating system we present a seismic model in which the source of tremor is the acoustic resonance of a fluid filled volcanic pipe triggered by excess gas pressure the model consists of three elements namely a triggering mechanism a resonator and a radiator for simplicity we assume a hemispherical trigger cylindrial resonator and circular radiator set in a vertical configuration with the trigger capping the top of the pipe and the disk shaped radiator shutting off its bottom considering the simple case of a source buried in a homogeneous half space we then apply the discrete wave number method to obtain a complete representation of the ground motion response at near and intermediate distances the results demonstrate that the displacement attributed to the pipe dominates the near field motion while that due to the disk is representative of the intermediate and far fields the trigger itself has a smaller contribution mainly limited to the field in the proximity of the source the characteristics displayed by the free surface response evolve from a strong impulsive signature in the immediate vicinity of the epicenter to a well developed harmonic wave train dominated by rayleigh waves at larger distances no clear shear arrival can be detected in the synthetic seismograms the displacement spectrum reflects the organ pipe modes of the conduit and the bandwidth associated with the dominant spectral peak of motion is controlled by the combined losses due to viscous attenuation in the fluid and elastic radiation into the solid in the case of the cylindrical magma column considered the radiation loss is proportional to the square of the pipe radius while the loss related to viscous damping is inversely proportional to the same factor indicating that the relative importance of the two loss mechanisms is critically dependent on the geometry of the magma reservoir the relative importance of the pipe and disk elements likewise is a function of the conduit cross section this suggests the possibility of determining the geometry of the source as well as the radiation loss and in situ magma viscosity from a comparison of near and far field observations

short term discharge and suspended sediment fluctuations in the proglacial skeldal river north east greenland the summer discharge pattern of the skeldal river which drains a 560 km2 partly glacierized catchment in north east greenland is dominated by diurnal oscillations reflecting variations in the melt rate of snow and ice in the basin superimposed on this diurnal pattern are numerous short lived discharge fluctuations of irregular periodicity and magnitude the larger fluctuations are described and attributed to both rainfall events and periodic collapse of the glacier margin damming flow from beneath the skelbrae glacier other minor fluctuations are less readily explained but are associated with changes in the channelized and distributed reservoirs and possibly temporary blockage of subglacial conduits caused by ice melt with subsequent damming fluctuations in suspended sediment concentration ssc are normally associated with discharge fluctuations although examples of transient flushes were observed where marked increases in ssc occurred in the absence of corresponding discharge variations a strong relationship between the event discharge increase and event ssc increase for rainfall induced events was established but no such relationship existed for non rainfall induced events there is some evidence for an exhaustion effect in the ssc patterns both at the event time scale and as the month proceeds a mean suspended sediment load of 1765 0 26 t day 1 was estimated for the study period which would be equivalent to a suspended sediment yield of 732 4 t km 2 year 1 copyright 2001 john wiley sons ltd

volcanic tremors at etna a model for hydraulic origin spectral analysis of volcanic tremors at etna reveals 2 characteristic features 1 the escaping vapor and gas streams drive the magma in random turbulent motion generating a seismic wave emission according toa f n exp m f k wherea amplitude f frequency n m k parameters 2 sudden transients in the boundary conditions of the piping system generate hydraulic pressure pulses which in the seismogram appear as wave packages their monochromatic frequency content suggests that their origin is associated with strong resonance phenomena it is conclusive that volcanic eruptions may be distinguished and classified according to different eigenfrequencies of the fissure and pipe network involved in the volcanic activity furthermore positions of active dykes may be traced following the amplitude distributions associated with pronounced spectral lines in the tremor spectra

assessing the transferability and robustness of an enhanced temperature index glacier melt model we investigate the transferability of an enhanced temperature index melt model that was developed and tested on haut glacier d arolla switzerland in the 2001 season the model s empirical parameters temperature factor tf and shortwave radiation factor srf are recalibrated for 1 other locations on haut glacier d arolla 2 subperiods of distinct meteorological conditions 3 different years on haut glacier d arolla and 4 other glaciers in different years the model parameters are optimized against simulations of an energy balance model validated against ablation observations results are compared with those obtained with the original parameters the model works very well when applied to other sites seasons and glaciers with the exception of overcast conditions differences are due to underestimation of high melt rates the parameter values are associated with the prevailing energy balance conditions showing that high srf are obtained on clear sky days whereas higher tf are typical of locations where glacier winds prevail and turbulent fluxes are high we also provide a range of parameters clearly associated with the site s location and its meteorological characteristics that could help to assign parameter values to sites where few data are available

glaciation of scotland during the younger dryas stadial a review considerable research effort in recent years has refined our knowledge of the landforms and sedimentary sequences generated by glaciers of the last cold episode in scotland the younger dryas and the chronology of their formation these recent advances have benefited from new technologies such as high resolution digital terrain models numerical ice sheet simulations and cosmic ray surface exposure dating this review presents a summary of recent younger dryas research in scotland particularly focusing on studies from the last two decades the findings of these investigations are compared with those of earlier workers and critically are appraised in terms of their relevance to five key aspects of younger dryas glaciation growth of the ice cap its geometry dynamics mode of decay and likely stadial climate minor differences aside the consensus view appears to favour a relatively thick but areally restricted central ice cap surrounded by satellite ice fields at its centre the ice cap was probably thickest during the first 500 years of the stadial after which it most likely thinned as the climate became increasingly arid outlet glaciers were dynamic only in marginal areas where they produced abundant moraine sequences during deglaciation but throughout much of the stadial central areas of the ice cap effected only minimal landscape modification deglaciation took place by active retreat under a rapidly warming climate overall the stadial was dominated by much colder and drier conditions than present and probably experi enced greater intra annual temperature ranges copyright 2009 john wiley sons ltd