The origin of accretionary lapilli

Experimental investigations in a recirculating wind tunnel of the mechanisms of formation of accretionary lapilli have demonstrated that growth is controlled by collision of liquid-coated particles, due to differences in fall velocities, and binding as a result of surface tension forces and secondary mineral growth. The liquids present on particle surfaces in eruption plumes are acid solutions stable at 100% relative humidity, from which secondary minerals, e.g. calcium sulphate and sodium chloride, precipitate prior to impact of accretionary lapilli with the ground. Concentric grain-size zones within accretionary lapilli build up due to differences in the supply of particular particle sizes during aggregate growth. Accretionary lapilli do not evolve by scavenging of particles by liquid drops followed by evaporation — a process which, in wind tunnel experiments, generates horizontally layered hemispherical aggregates. Size analysis of particles in the wind tunnel air stream and particles adhering to growing aggregates demonstrate that the aggregation coefficient is highly grain-size dependent. Theoretical simulation of accretionary lapilli growth in eruption plumes predicts maximum sizes in the range 0.7–20 mm for ash cloud thicknesses of 0.5–10 km respectively.     

Measurements of electric charge distribution in volcanic plumes at Sakurajima Volcano,Japan

Measurements of perturbations in the atmospheric potential gradient around volcanic plumes at multiple (from two to five) sites, and measurements of the charge-mass ratio of ash particles falling from volcanic plumes, were carried out at Sakurajima Volcano, Japan. Results from 28 and 29 October 1995, show that the nature of the perturbations depends on the intensity of plume activity. Although plume activity was vigorous on 28 October, negative perturbations were predominant. As plume activity peaked, the magnitude of negative perturbations decreased just below the plume and increased at an off-axis site. During the peak period, positively charged ash particles fell out from the plume. This suggests that the active plume dominantly contained negatively charged materials, and that positively and negatively charged materials were added to the lower and upper parts of the plume, respectively, during the peak period. On the other hand, as plume activity became less vigorous on 29 October, the perturbations were characterized by a positive anomaly followed by a negative anomaly. Because wind velocity increased with altitude that day, we infer that positive and negative charges were distributed in the upper and lower parts of the plume, respectively. The differences in perturbations observed on 28 and 29 October suggest that volcanic plumes are generally composed of three parts: an upper part with positively charged gas and aerosol, a middle part with negatively charged fine ash particles, and a lower part with positively charged coarse ash particles. The compilation of present and previous results from Sakurajima and other volcanoes indicates that the effect of the negative charge in the middle part was predominant in most cases, although positive perturbations caused by the upper part were observed around some weak plumes. The effect of positively charged particles in the lower part was observable only when plume activity was sufficiently strong because positively charged coarse particles tended to fall out near the vent.     

Sedimentation model and analysis for differential settling of two - particle - size suspensions in the Stokes region

A model is developed for predicting the settling velocity in suspensions of particles of two different sizes based on experimental data for the settling rate of two-size suspensions in various liquids using particles of equal density. In these experiments, the retarding effect of the smaller particles on the settling velocities of the larger ones is taken into account. The model considers Steinour’s fundamental equation and assumes a fixed arrangement of particles and constant velocity in a single-size susp...     

Ejecta velocities,magma chamber pressure and kinetic energy associated with the 1968 eruption of arenal volcano

During the initial explosive phase of the eruption of Arenal volcano small projectiles were thrown a maximum distance of 5 km. Considering the effect of atmospheric drag these projectiles must have had initial velocities of at least 600 m/sec. For this velocity, the gas pressure in the magma chamber must have reached at least 4700 bars and the kinetic energy of the initial explosion is estimated as 2.4 ± 1.2 × 10^a ergs. Had the effect of aerodynamic braking been ignored in making these calculations, as has always been done in the past, the calculated initial velocity would have been 220 m/sec; chamber pressure and kinetic energy estimates would thus be substantially lower. Clearly, velocities of ejecta, chamber pressures and kinetic energies for many explosive volcanic events have been seriously underestimated in the recent past, as has been the ability of overlying materials to contain, in certain cases, tremendous overpressures for short periods of time. A projectile with an initial velocity of 600 m/sec would have a maximum range of more than 200 km on the moon. Thus, the presence of far-reaching secondary crater fields on the moon cannot, at this time, be considered evidence for an impact origin of the parent crater. 600 m/sec is not the upper limit for initial velocities of volcanic ejecta. There is some indication that such velocities could reach values greater than 2 km/sec, suggesting that volcanic as well as impact mechanisms may be able to impart escape velocity to lunar materials.     

Continuous monitoring of volcanic eruption dynamics: a review of various techniques and new results from a frequency-modulated radar Doppler system

 In situ measurement of volcanic eruption velocities is one of the great challenges left in geophysical volcanology. In this paper we report on a new radar Doppler technique for monitoring volcanic eruption velocities. In comparison with techniques employed previously (e.g., photographic methods or acoustic Doppler measurements), this method allows continuous recordings of volcanic eruptions even during poor visibility. Also, radar Doppler instruments are usually light weight and energy efficient, which makes them superior to other Doppler techniques based on laser light or sound. The proposed new technique was successfully tested at Stromboli Volcano in late 1996 during a period of low activity. The recorded data allow a clear distinction between particles rising from the vent and particles falling back towards the vent. The mean eruption velocity was approximately 10 m/s. Most of the eruptions recorded by radar were correlated to seismic recordings. The correlation between the magnitude of the volcanic shocks and the eruption force index defined in the paper may provide new insights into magma transport in the conduit. Received: 15 May 1998 / Accepted: 15 December 1998     

Grain size,areal thickness distribution and controls on sedimentation of the 1991 Mount Pinatubo tephra layer in the South China Sea

The June 15, 1991 climactic eruption of Mt. Pinatubo produced an extensive, largely co-ignimbrite-derived airfall ash layer on Luzon Island and across the central South China Sea. The layer covers an area of ~4×10⁵ km² with a volume of 5.5 km³. Near the coast of Luzon, the deposit consists of two units: a normally graded basal ash bed, unimodal in grain size, and a finer-grained, internally structureless upper ash bed showing grain size bimodality. With increasing distance from the source, the coarse particle populations of the two units merge and migrate towards a near-constant fine population (~11 m); the distal region is covered by a fine-mode dominated, virtually ungraded single ash layer. The reversal of the winds from easterly directions at upper-tropospheric and stratospheric levels to westerly directions in the middle and lower troposphere indicates that both the coarse- and fine-mode components fell out from high-altitude eruption clouds. The high-velocity upper-level winds, however, would have transported fine-grained ash particles far beyond the South China Sea, which suggests that their settling was accelerated by aggregation. The boundary between the units thus marks a change from fallout of predominantly discrete pyroclasts to simultaneous fallout of aggregated fines and freely falling, coarse-grained particles. The particle populations composing the upper ash bed were almost completely removed from the proximal areas by the upper-level winds. At lower elevations, the counterclockwise circulation of a typhoon over the coastal area advected the ash south and eastward, producing a thickness maximum in the medial region (at about 160 km from source). The strong displacement of fines, possibly aided by wind turbulence, led to a break in bulk tephra thinning rates close to the coastline. In the distal region, outside the influence of the typhoon, southwest monsoonal winds caused a distinct lobe axis inflection and thickness asymmetry. Within this region, at about 420 km from source, fallout of particle aggregates created a second thickness maximum. Comparison of the field data with previous experimental observations and tephra flux records in the deep sea (Wiesner et al. 1995; Carey 1997; McCool 2002) implies that the transport of ash in the water column was largely determined by vertical density currents. Differences in the reaction of coarse and fine particles to turbulence in the descending plumes probably suppressed the segregation of fines but allowed the coarser pyroclasts to maintain their initial order of arrival at the sea surface. Considering typical fall rates of convective plumes, modifications of the initial fallout position of the particles by the South China Sea current system are on the order of only a few kilometers. The results suggest that convective sedimentation processes ensure the preservation of atmospheric particle transport directions, distances, and fallout modes in the deep sea.Editorial responsibility: R. CioniAn erratum to this article can be found at     

Controls on sonic velocity in carbonates

Compressional and shear-wave velocities (V _p andV _s) of 210 minicores of carbonates from different areas and ages were measured under variable confining and pore-fluid pressures. The lithologies of the samples range from unconsolidated carbonate mud to completely lithified limestones. The velocity measurements enable us to relate velocity variations in carbonates to factors such as mineralogy, porosity, pore types and density and to quantify the velocity effects of compaction and other diagenetic alterations.Pure carbonate rocks show, unlike siliciclastic or shaly sediments, little direct correlation between acoustic properties (V _p andV _s) with age or burial depth of the sediments so that velocity inversions with increasing depth are common. Rather, sonic velocity in carbonates is controlled by the combined effect of depositional lithology and several post-depositional processes, such as cementation or dissolution, which results in fabrics specific to carbonates. These diagenetic fabrics can be directly correlated to the sonic velocity of the rocks.At 8 MPa effective pressureV _p ranges from 1700 to 6500 m/s, andV _s ranges from 800 to 3400 m/s. This range is mainly caused by variations in the amount and type of porosity and not by variations in mineralogy. In general, the measured velocities show a positive correlation with density and an inverse correlation with porosity, but departures from the general trends of correlation can be as high as 2500 m/s. These deviations can be explained by the occurrence of different pore types that form during specific diagenetic phases. Our data set further suggests that commonly used correlations like Gardner's Law (V _p-density) or the time-average-equation (V _p-porosity) should be significantly modified towards higher velocities before being applied to carbonates.The velocity measurements of unconsolidated carbonate mud at different stages of experimental compaction show that the velocity increase due to compaction is lower than the observed velocity increase at decreasing porosities in natural rocks. This discrepancy shows that diagenetic changes that accompany compaction influence velocity more than solely compaction at increasing overburden pressure.The susceptibility of carbonates to diagenetic changes, that occur far more quickly than compaction, causes a special velocity distribution in carbonates and complicates velocity estimations. By assigning characteristic velocity patterns to the observed diagenetic processes, we are able to link sonic velocity to the diagenetic stage of the rock.     

Volcanic particle aggregation in explosive eruption columns. Part II: Numerical experiments

The goal of this paper is to determine the parameters that control the aggregation efficiency and the growth rate of volcanic particles within the eruption column. Numerical experiments are performed with the plume model ATHAM (Active Tracer High resolution Atmospheric Model). In this study we employ the parameterizations described in a companion paper (this issue). The presence of hydrometeors promotes the aggregation of ash particles, which strongly increases their fall velocities and thus their environmental impact. The tephra mass is about two orders of magnitude greater than that of hydrometeors during typical Plinian eruptions without interaction of external water. Ice is highly dominant in comparison to liquid water (> 99% by mass). This is caused by the fast column rise (> 100 m s^− 1 on average) to very cold altitudes. Most particles occur in the form of frozen aggregates with low ice content.The collection efficiency is governed by the availability of hydrometeors acting as adhesives at the particles’ surface in our study, and wet ash particles have a higher sticking capacity than icy ones. Therefore, aggregation is fastest during the eruption within the column when limited regions of liquid water exist and when particle concentrations are very high (of the order of 10⁵ cm^− 3). Increased humidity in the background atmosphere generally leads to enhanced ice formation, but shows only a weak influence on the aggregation process. First sensitivity studies showed, however, a significant increase of the liquid water fraction when considering salinity effects. The availability of water or ice at the particles' surfaces is also governed by the surface properties, the porosity and permeability of ash, which are not well established to date. Particle growth is significantly enhanced for greater differences in the sizes and fall velocities among particles, as gravitational capture becomes more efficient. Our experiments indicate a major influence of the erupted particle size distribution. First sensitivity studies show that electrostatic forces result in a significant enhancement of aggregated particles.The present exploratory study provides new insights into the sensitivity of the ash aggregation process to a number of key parameters. Our results indicate the need of further constraining particle composition, size, porosity, permeability, and surface properties at low temperatures by in situ observations in the laboratory and in the field. In addition further research on electrostatic aggregation would be desirable.     

Bifurcation of volcanic plumes in a crosswind

Bent-over buoyant jets distorted by a crosscurrent develop a vortex pair structure and can bifurcate to produce two distinct lobes which diverge from one another downwind. The region downwind of the source between the lobes has relatively low proportions of discharged fluid. Factors invoked by previous workers to cause or enhance bifurcation include buoyancy, release of latent heat at the plume edge by evaporating water droplets, geometry and orientation of the source, and the encounter with a density interface on the rising path of the plume. We suggest that the pressure distribution around the vortex pair of a rising plume may initially trigger bifurcation. We also report new experimental observations confirming that bifurcation becomes stronger for stronger bent-over plumes, identifying that bifurcation can also occur for straight-edged plumes but gradually disappears for stronger plumes which form a gravity current at their final level and spread for a significant distance against the current. Observations from satellites and the ground are reviewed and confirm that volcanic plumes can show bifurcation and a large range of bifurcation angles. Many of the bifurcating plumes spread out at the tropopause level and suggest the tropopause may act on the plumes as a density interface enhancing bifurcation. Even for quite moderate bifurcation angles, the two plume lobes become rapidly separated downwind by distances of tens of kilometers. Such bifurcating plumes drifting apart can only result in bilobate tephra fall deposits. The tephra fall deposit from the 16 km elevation, SE spreading, bifurcating volcanic plume erupted on 15 May 1981 from Mt Pagan was sampled by previous workers and clearly displayed bilobate characteristics. Examples of bilobate tephra fall deposits are reviewed and their origin briefly discussed. Bilobate deposits are common and may result from many causes. Plume bifurcation should be considered one of the possible mechanisms which can account for come examples of bilobate tephra fall deposits.     

Doppler radar sounding of volcanic eruption dynamics at Mount Etna

Besides their common use in atmospheric studies, Doppler radars are promising tools for the active remote sensing of volcanic eruptions but were little applied to this field. We present the observations made with a mid-power UHF Doppler radar (Voldorad) during a 7-h Strombolian eruption at the SE crater of Mount Etna on 11–12 October 1998. Main characteristics of radar echoes are retrieved from analysis of Doppler spectra recorded in the two range gates on either side of the jet axis. From the geometry of the sounding, the contribution of uprising and falling ejecta to each Doppler spectrum can be discriminated. The temporal evolution of total power backscattered by uprising targets is quite similar to the temporal evolution of the volcanic tremor and closely reproduces the overall evolution of the eruption before, during and after its paroxysm. Moreover, during the sharp decrease of eruptive activity following the paroxysm, detailed analysis of video (from camera recording), radar and seismic measurements reveals that radar and video signals start to decrease simultaneously, approximately 2.5 min after the tremor decline. This delay is interpreted as the ascent time through a magma conduit of large gas slugs from a shallow source roughly estimated at about 500 m beneath the SE crater. Detailed analysis of eruptive processes has been also made with Voldorad operating in a high sampling rate mode. Signature of individual outburst is clearly identified on the half part of Doppler spectra corresponding to rising ejecta: temporal variations of the backscattered power exhibit quasi periodic undulations, whereas the maximum velocity measured on each spectrum displays a sharp peak at the onset of each outburst followed by a slow decay with time. Periodicity of power variations (between 3.8 and 5.5 s) is in agreement with the occurrence of explosions visually observed at the SE vent. Maximum vertical velocities of over 160 m s^–1 were measured during the paraoxysmal stage and the renewed activity. Finally, by using a simplified model simulating the radar echoes characteristics, we show that when Voldorad is operating in high sampling rate mode, the power and maximum velocity variations are directly related to the difference in size and velocity of particles crossing the antenna beam.Editorial responsibility: A. Woods     

Unique phenomenon of the accumulation of <Emphasis Type="Italic">terrestrial</Emphasis> metal iron particles in lacustrine deposits: Zhombolok volcanic region,East Sayan

The native iron particles that were previously detected by thermomagnetic and microprobe analyses in the sediments of different age in many regions of the world are of extraterrestrial origin. The similarity in the compositions, grain shapes, and sizes observed in the extraterrestrial and terrestrial particles of native iron testifies to the common production conditions of iron particles during the formation of planets. In this paper, the single finding of terrestrial iron in the lacustrine sediments of the Zhombolok volcanic region, East Sayan, is discussed. The uniqueness of the results indicates that the spatial distribution of the particles of native iron is limited to a fairly narrow area around their source—volcanic eruption or/and the fall of a large meteorite.     

Sedimentation of tephra by volcanic plumes: I. Theory and its comparison with a study of the Fogo A plinian deposit,Sao Miguel (Azores)

Sedimentation of ejecta from volcanic plumes has been studied as a function of distance from the source in the Fogo A plinian deposit, Sao Miguel, Azores. The Fogo A trachytic pumice deposit is reversely graded and can be divided into two parts on the basis of pumice colour, abundance of syenite accessory lithic clasts and distribution. The lower syenite-poor part was dispersed to the south and was clearly influenced by wind. The upper syenite-rich part is coarsegrained and has a nearly symmetrical distribution around the vent. Elongation of isopachs to the east indicate a weak wind influence. The grain-size variations of lithic and crystal components in the upper coarse part were studied. Total accumulation and accumulation per unit area (expressed in kg/m²) show good fits to a gaussian function at distances greater than 7 km for grain diameters less than 2 cm. These results agree with a theoretical model for a radially spreading turbulent current moving over a quiescent fluid. The gaussian coefficient is shown to be a function of grain size and the flow rate of material into the umbrella region of the eruption column. The coefficient is therefore also a function of column height. The column height deduced from these data is 21 km, which is in broad agrrement with the column height of 27 km deduced from maximum clast dispersal using the method of Carey and Sparks (1986). The accumulation of clasts larger than 2 cm agrees with a theory for the fallout of clasts from the margins of the ascending eruption column, which treats the plume as a succession of large eddies that decrease their mass of particles as an exponential function of time. Calculations are also presented for the influence of the radial inflow of surrounding air into the column on the deposition of clasts. These calculations constrain the wind speed during the later part of the Fogo A eruption to be at most a few metres per second. The study has allowed four different dynamic categories of clast behaviour to be recognised in eruption columns.     

Do conservative solutes migrate at average pore-water velocity?

According to common understanding, the advective velocity of a conservative solute equals the average linear pore-water velocity. Yet direct monitoring indicates that the two velocities may be different in heterogeneous media. For example, at the Camp Dodge, Iowa, site the advective velocity of discrete Cl- plumes was less than one tenth of the average pore-water velocity calculated from Darcy's law using the measured hydraulic gradient, effective porosity, and hydraulic conductivity (K) from large-scale three-dimensional (3D) techniques, e.g., pumping tests. Possibly, this difference reflects the influence of different pore systems, if the K relevant to transient solute flux is influenced more by lower-K heterogeneity than a steady or quasi-steady water flux. To test this idea, tracer tests were conducted under controlled laboratory conditions. Under one-dimensional flow conditions, the advective velocity of discrete conservative solutes equaled the average pore-water velocity determined from volumetric flow rates and Darcy's law. In a larger 3D flow system, however, the same solutes migrated at approximately 65% of the average pore-water velocity. These results, coupled with direct observation of dye tracers and their velocities as they migrated through both homogeneous and heterogeneous sections of the same model, demonstrate that heterogeneity can slow the advective velocity of discrete solute plumes relative to the average pore-water velocity within heterogeneous 3D flow sytems.     

Gravity currents down canyons: effects of rotation

The flow of dense water in a V-shaped laboratory-scale canyon is investigated by using a non-hydrostatic numerical ocean model with focus on the effects of rotation. By using a high-resolution model, a more detailed analysis of plumes investigated in the laboratory (Deep-Sea Res I 55:1021–1034 2008) for laminar flow is facilitated. The inflow rates are also increased to investigate plume structure for higher Reynolds numbers. With rotation, the plumes will lean to the side of the canyon, and there will be cross-canyon geostrophic currents and Ekman transports. In the present study, it is found that the cross-canyon velocities are approximately 5 % of the down-canyon velocities over the main body of the plume for the rotational case. With rotation, the flow of dense water through the body of the plume and into the plume head is reduced. The plume head becomes less developed, and the speed of advance of the head is reduced. Fluid parcels near the top of the plume will to a larger extent be left behind the faster flowing dense core of the plume in a rotating system. Near the top of the plume, the cross-canyon velocities change direction. Inside the plume, the cross-flow is up the side of the canyon, and above the interface to the ambient there is a compensating cross-flow down the side of the canyon. This means that parcels of fluid around the interface become separated. Parcels of fluid around the interface with small down-canyon velocity components and relative large cross-canyon components will follow a long helix-like path down the canyon. It is found that the entrainment coefficients often are larger in the rotational experiments than in corresponding experiments without rotation. The effects of rotation and higher inflow rates on the areal patterns of entrainment velocities are demonstrated. In particular, there are bands of higher entrainment velocities along the lateral edges of the plumes in the rotational cases.     

Sedimentation of tephra by volcanic plumes. Part 2: controls on thickness and grain-size variations of tephra fall deposits

A model for sedimentation from turbulent suspensions predicts that tephra concentration decreases exponentially with time in an ascending volcanic column and in the overlying umbrella cloud. For grain-size distributions typical of plinian eruptions application of the model predicts for thickness variations in good agreement with the exponential thinning observed in tephra fall deposits. The model also predicts a proximal region where fallout from the plume margins results in a more rapid decrease in thickness so that the deposit shows two segments on a thickness versus distance plot. Several examples of deposits with two segments are known. The distance at which the two segments intersect is a measure of eruption column height. The thickness half-distance ( equivalent to the dispersal index of Walker) is strongly correlated with column height, but is also weakly dependent on grain-size distribution of the ejecta. For a dispersal index of 500 km² (the plinian/subplinian boundary of Walker) column heights between 14 and 18 km are calculated. For ultraplinian deposits with D>50000 km² column heights of at least 45 km are implied. Model grain-size distributions of the deposits have sorting values comparable to those observed in tephra fall deposits formed from eruption columns in a weak or negligible cross-wind. Median diameter decreases exponentially with distance as is observed. Sorting () improves with distance as is observed in plinian deposits in a weak wind. However, tephra fall deposits formed in strong winds do not show improved sorting with distance and proximal deposits are typically somewhat better sorted than the model calculations. Differences are attributed to the influence of wind which disperses particles further than predicted in our model and which has an increasing influence as particle size decreases.     

Surface area, porosity and water adsorption properties of fine volcanic ash particles

Our understanding on how ash particles in volcanic plumes react with coexisting gases and aerosols is still rudimentary, despite the importance of these reactions in influencing the chemistry and dynamics of a plume. In this study, six samples of fine ash (<100 m) from different volcanoes were measured for their specific surface area, a_s, porosity and water adsorption properties with the aim to provide insights into the capacity of silicate ash particles to react with gases, including water vapour. To do so, we performed high-resolution nitrogen and water vapour adsorption/desorption experiments at 77 K and 303 K, respectively. The nitrogen data indicated a_s values in the range 1.1–2.1 m²/g, except in one case where a a_s of 10 m²/g was measured. This high value is attributed to incorporation of hydrothermal phases, such as clay minerals, in the ash surface composition. The data also revealed that the ash samples are essentially non-porous, or have a porosity dominated by macropores with widths >500 Å. All the specimens had similar pore size distributions, with a small peak centered around 50 Å. These findings suggest that fine ash particles have relatively undifferentiated surface textures, irrespective of the chemical composition and eruption type. Adsorption isotherms for water vapour revealed that the capacity of the ash samples for water adsorption is systematically larger than predicted from the nitrogen adsorption a_s values. Enhanced reactivity of the ash surface towards water may result from (i) hydration of bulk ash constituents; (ii) hydration of surface compounds; and/or (iii) hydroxylation of the surface of the ash. The later mechanism may lead to irreversible retention of water. Based on these experiments, we predict that volcanic ash is covered by a complete monolayer of water under ambient atmospheric conditions. In addition, capillary condensation within ash pores should allow for deposition of condensed water on to ash particles before water reaches saturation in the plume. The total mass of water vapour retained by 1 g of fine ash at 0.95 relative water vapour pressure is calculated to be ~10^–2 g. Some volcanic implications of this study are discussed.Editorial responsibility: J. Gilbert     

Dilute gravity current and rain-flushed ash deposits in the 1.8 ka Hatepe Plinian deposit,Taupo, New Zealand

Two groups of poorly sorted ash-rich beds, previously interpreted as rain-flushed ashes, occur in the ca. AD 180 Hatepe Plinian pumice fall deposit at Taupo volcano, New Zealand. Two ash beds with similar dispersal patterns and an aggregate thickness of up to 13 cm make up the lowermost group (A). Group A beds extend 45 km north-east of the vent and cover 290 km². In the southern part of the group A distribution area, a coarse ash to lapilli-size Plinian pumice bed (deposit B) separates the two group A beds. The scarcity of lapilli (material seen elsewhere from the still-depositing pumice fall) in group A beds indicates that they were rapidly transported and deposited. However, this rapid transportation and deposition did not produce cross-bedding, nor did it erode the underlying deposits. It is proposed that thick (>600 m) but dilute gravity currents generated from the collapsing outer margin of the otherwise buoyant Hatepe Plinian eruption column deposited the group A beds. The upper ash beds (group C) consist of one to seven layers, attain an aggregate thickness of 35 cm, and vary considerably in thickness and number of beds with respect to distance from vent. Group C beds contain variable amounts of ash mixed with angular Plinian pumices and are genuine rain-flushed ashes. Several recent eruptions at other volcanoes (Ukinrek Maars, Vulcan, Rabaul, La Soufrère de Guadeloupe and Soufrière, St Vincent) have produced gravity currents similar in style, but much smaller than those envisaged for group A deposits. The overloaded margins of otherwise buoyant eruption plumes generated these gravity currents. Laboratory studies have produced experimental gravity current analogues. Hazards from dilute gravity currents are considerable but often overlooked, thus the recognition of gravity current deposits will contribute to more thorough volcanic hazard assessment of prehistoric eruption sequences.     

Base surge in recent volcanic eruptions

A base surge, first identified at the Bikini thermonuclear undersea explosion, is a ring-shaped basal cloud that sweeps outward as a density flow from the base of a vertical explosion column. Base surges are also common in shallow underground test explosions and are formed by expanding gases which first vent vertically and then with continued expansion rush over the crater lip (represented by a large solitary wave in an underwater explosion), tear ejecta from it, and feed a gas-charged density flow, which is the surge cloud. This horizontally moving cloud commonly has an initial velocity of more than 50 meters per second and can carry clastic material many kilometers. Base surges are a common feature of many recent shallow, submarine and phreatic volcanic eruptions. They transport ash, mud, lapilli, and blocks with great velocity and commonly sandblast and knock down trees and houses, coat the blast side with mud, and deposit ejecta at distances beyond the limits of throw-out trajectories. Close to the eruption center, the base surge can erode radial channels and deposit material with dune-type bedding.     

Compressional- and shear-wave velocities as a function of confining stress in dry sandstones

A laboratory study was carried out to investigate the influence of confining stress on compressional- and shear-wave velocities for a set of rock samples from gas-producing sandstone reservoirs in the Cooper Basin, South Australia. The suite of samples consists of 22 consolidated sublitharenites with helium porosity ranging from 2.6% to 16.6%. We used a pulse-echo technique to measure compressional- and shear-wave velocities on dry samples (cylindrical 4.6 × 2 cm) at room temperature and at elevated confining stress (≤ 60 MPa). Compressional- and shear-wave velocities in samples increase non-linearly with confining stress. A regression equation of the form V = A ? Be^?DP gives a good fit to the measured velocities with improved prediction of velocities at high confining stresses compared with equations suggested by other studies. The predicted microcrack-closure stresses of the samples show values ranging from 70 MPa to 95 MPa and insignificant correlation with porosity, permeability or clay content. There is a positive correlation between change in velocity with core porosity and permeability, but this association is weak and diminishes with increasing confining stress. Experimental results show that pore geometry, grain-contact type, and distribution and location of clay particles may be more significant than total porosity and clay content in describing the stress sensitivity of sandstones at in situ reservoir effective stress. The stress dependence of Cooper Basin sandstones is very large compared with data from other studies. The implication of our study for hydrocarbon exploration is that where the in situ reservoir effective stress is much less than the microcrack-closure stress of the reservoir rocks, the variation of reservoir effective stress could cause significant changes in velocity of the reservoir rocks. The velocity changes induced by effective stress in highly stress-sensitive rocks can be detected at sonic-log and probably surface-seismic frequencies.