Distribution of thermal areas on an active lava flow field: Landsat observations of Kilauea,Hawaii, July 1991

A Landsat Thematic Mapper (TM) image acquired on 23 July 1991 recorded widespread activity associated with the Episode 48 of the Pu'u 'O'o-Kupaianaha eruption of Kilauea Volcano, Hawaii. The scene contains a very large number (>3500) of thermally elevated near infrared (0.8–2.35 m) pixels (each 900 m²), which enable the spatial distribution of volcanic activity to be identified. This activity includes a lava lake within Pu'u 'O'o cone, an active lava tube system (7.9 km in length) with skylights between the Kupaianaha lava shield and several ocean entry points, and extensive active surface flows (total area of 1.3 km²) within a much larger area of cooling flows (total16 km²). The production of an average flux density map from the TM data of the flow field, wherein the average flux density is defined in units of Wm^-2, allows for the chronology of emplacement of active and cooling flows to be determined. The flux density map reveals that there were at least three breakouts (>5000 Wm^-2) feeding active flows, but on the day that the data were collected the TM recorded a waning phase of surface activity in this area, based on the relatively large amount of intermediate power-emitting (cooling) flows compared to high power-emitting (active) flows. The production of a comparable flux density map for future eruptions would aid in the assessment of volcanic hazards if the data were available in near-real time.     

Cooling history and differentiation of a thick North Mountain Basalt flow (Nova Scotia,Canada)

A thick (<175 m) North Mountain Basalt flow at McKay Head, Nova Scotia (Canada) shows 25-cm-thick differentiated layers separated by 130 cm of basalt in its upper 34m. Upper layers (5 m below the lava top) are highly vesicular whereas lower ones are pegmatitic and contain a thin (2 cm) rhyolite band. The layering of the flow closely resemble that of some Hawaiian lava lakes. The eesicular basalts and mafic pegmatites are inferred to be liquid-rich segregations which drained into horizontal cracks that formed within a crystalline mush. The cracks resulted from a thermal contraction associated with cooling and shrinkage of the mush. Rhyolites were formed by in situ differentiation. Gas overpressures fractured the pegmatites and gas effervescence filter pressing forced silicarich residual liquid from pegmatite interstices into the fractures creating bands. Chemical differences between the pegmatitic layers and early formed, highly differentiated upper vesicular layers may reflect a role for volatiles in the differentiation process along with crystal fractionation.     

Mixed-magma pyroclastic surge deposits associated with debris avalanche deposits at Colima volcanoes,Mexico

The transition between the terminal cones and the ancestral edifices of Nevado de Colima and Fuego de Colima volcanoes is marked by the deposits of gigantic volcanic debris avalanches of the Mount St. Helens (MSH) or Bezymianny type. Unusual mafic juvenile fragments and cauliflower bombs as well as juvenile fragments of mixed and more evolved composition are abundant in dune-bedded pyroclastic-surge deposits directly associated with these catastrophic events at both volcanoes. At Nevado, these mafic juvenile fragments represent the most primitive magma ever erupted by the volcano (SiO₂52.50%). The lavas directly preceding and following the debris-avalanche event are silicic andesites (SiO₂59%). At Fuego these juvenile fregments have 56% SiO₂. The lavas from the upper parts of the caldera wall are dacites (65% SiO₂), whereas the terminal cone is composed of andesites (57% to 62% SiO₂). At Nevado, petrologic evidence for interaction of mafic magma with andesitic or dacitic magma in a high-level magma chamber, just before the eruption of pyroclastic surge deposits, consists of: (1) banded juvenile bombs of intermediate composition; (2) the range of composition of these bombs from SiO₂52% to 58%; (3) the presence of highly magnesian olivine with reaction rims; (4) inverse zoning in clinopyroxene with strong Mg enrichment towards the rim; (5) resorption of plagioclase; and (6) significant compositional heterogeneity in the vitric phase. Volcanic debris-avalanche events at Nevado and Fuego de Colima may thus correspond with major breaks in the petrological evolution of the volcanoes and the start of a new magmatic cycle. Injection of mafic magma into the presently perched viscous surface dome of the active Fuego cone, as occurred in 1818 and 1913, could enhance the likelihood of southward collapse of the flank of an already unstable edifice, and it must be considered in future hazard assessment of this active volcano. Risk to life and property for the entire Colima region associated with such catastrophic phenomena would be immeasurably greater in comparison with hazards related to the last explosive outburst in 1913, which resulted in emplacement of pyroclastic flows over uninhabited areas of the upper flanks of the volcano.     

Lunar and solar barometric tides at Nandi,Fiji

Lunar and solar atmospheric tidal oscillations have been determined with satisfactory accuracy from 17 years 11 months of mean sea-level barometric pressure observations taken at Nandi, Fiji. In many respects, the results are consistent with previous tidal determinations in the south-west Pacific region, although these are few and widely scattered. However, the mean annual amplitude of the lunar tide at Nandi, as determined in this study, 88 b, is much greater than might have been expected from currently available global amplitude maps. Nevertheless, the probable correctness of this result has been confirmed by the analysis of nearly 6 years of similar data from Nausori (130 km E.S.E. of Nandi), which yielded a mean annual lunar amplitude of 88 b, compared with a Nandi amplitude of 83 b for a closely corresponding period.     

Digital topography of volcanoes from radar interferometry: an example from Mt Vesuvius,Italy

A new airborne radar technique can generate digital topographic data for volcanoes at a scale of 10 m spatial and 1–5 m vertical, with a swath width of 6.4 km. Called TOPSAR, the intrument is an interferometric radar flown on the NASA DC-8 aircraft. TOPSAR data permit the quantification of volcano slopes, volumes, and heights, and as such will be valuable for the analysis of lava flows, domes, and lahar channels. This instrument will be flown over several volcanoes in the near future, providing volcanologists with valuable data sets for the analysis of high-resolution topography. We briefly illustrate the potential use of TOPSAR data through examples from Mt Somma and Vesuvius, Italy.     

Ion distribution dynamics near the Earth’s bow shock: first measurements with the 2D ion energy spectrometer CORALL on the INTERBALL/Tail-probe satellite

The dynamics of the ion distribution function near the Earths bow shock is studied on the basis of quasi-3D measurements of ion energy spectra in the range of 30–24200 eV/q with the Russian-Cuban CORALL instrument on the INTERBALL/Tail-probe satellite. The instrument was designed for observations of magnetospheric plasma and measures ions, in an angular range of 36°-144° from the Earth-Sun direction. Ion populations generated by the Earth bow shock are often observed upstream from the bow shock. In the solarwind stream compressed and heated by the passing of very dense magnetic cloud (CME), two types of these ion populations were measured upstream and before the bow shock crossing on 25 August 1995 at 07:37 UT. Both populations were observed in the energy range above 2 keV. At 06:20 UT, when the angle between the direction of the interplanetary magnetic field and normal to the bow shock _Bn was w 43° the instrument observed a narrow, fast (800 km/s) field-aligned beam moving from the Earth. At 07:30, when _Bn 28°, the wide ion pitch-angle distribution was observed. A similar suprathermal ion population is observed in the magnetosheath simultaneously with the solar-wind ion population being heated and deflected from the Sun-Earth direction. The similarity of observations during the mentioned time-interval and under usual solar-wind conditions allows us to conclude that types of suprathermal ion populations upstream and downstream from the bow shock do not depend on the solar-wind disturbance generated by magnetic cloud.     

On the hydrodynamic behavior of supercooled water drops interacting with columnar ice crystals

A numerical evaluation of the complete Navier-Stokes equations of motion for steady-state, incompressible flow past an infinite circular cylinder is given in terms of the stream function, vorticity, and pressure distribution past such bodies. A method is described which allows use of these flow characteristics: (1) to approximate the characteristics of air flow past hexagonal columnar ice crystals falling under gravity at terminal velocity in air, (2) to compute the trajectory of supercooled cloud drops relative to such ice crystals, and (3) to determine the efficiency with which short columnar ice crystals and needle shaped ice crystals collide with supercooled cloud drops. It is found that for all columnar type ice crystals riming is negligible if the cloud drop size is less than 5 m, and that for riming to commence short columnar crystals must have diameters larger than 50 m, while needle crystals must have diameters larger than 40 m. It is further shown that the collision efficiency cut-offs at the small drop radius and at the large drop radius end of the collision efficiency diagram can be explained on the basis of the cloud drop trajectories for these drop size ranges.     

Factors influencing the height of Hawaiian lava fountains: implications for the use of fountain height as an indicator of magma gas content

The heights of lava fountains formed in Hawaiian-style eruptions are controlled by magma gas content, volume flux and the amounts of lava re-entrainment and gas bubble coalescence. Theoretical models of lava fountaining are used to analyse data on lava fountain height variations collected during the 1983–1986 Pu'u 'O'o vent of Kilauea volcano, Hawaii. The results show that the variable fountain heights can be largely explained by the impact of variations in volume flux and amount of lava re-entrainment on erupting magmas with a constant gas content of 0.32 wt.% H₂O. However, the gas content of the magma apparently declined by 0.05 wt.% during the last 10 episodes of the eruption series and this decline is attributed to more extensive pre-eruption degassing due to a shallowing of the sub-vent feeder dike. It is concluded that variations in lava fountain height cannot be simply interpreted as variations in gas content, as has previously been suggested, but that fountain height can still be a useful guide to minimum gas contents. Where sufficient data are available on eruptive volume fluxes and extent of lava entrainment, greatly improved estimates can be made of magma gas content from lava fountain height.     

Microseismicity, tectonics and seismic potential in southern Caribbean and northern Venezuela

Approximately one thousand microearthquakes with body-wave magnitude mb have been located in northern Venezuela and the southern Caribbean region (9–12° N; 64–70° W) since the installation in 1980 of the Venezuelan Seismological Array, together with forty events of mb 4, one of them with surface-wave magnitude Ms 6. Focal depths are in the range of 0 to <15 km. This geologically complex region is part of the boundary between the Caribbean and the South American Plates. Epicentral locations indicate that this E–W oriented portion of the boundary is formed by two 400 km long subparallel fault zones: San Sebastián fault zone (SSF), 20 km north of Caracas along the coast; and La Victoria fault zone (LVF), 25 km south of the city. They are clearly delineated by the microseismicity. New composite focal mechanism solutions (CFMS) along these faults show right-lateral strike-slip (RLSS) motion on nearly E–W oriented fault planes. NW-striking subsidiary active faults occur in the region and intercept the two main E–W fault zones. These interceptions show high levels of microearthquake activity and seismic moment release when compared to other portions of both, the main and subsidiary faults. New CFMS at those fault crossing sites show NW-striking RLSS motion and normal faulting, in an en-echelon-like structural behavior. Geological data and quantitative comparisons with other transcurrent plate boundaries in the world suggest that the rate of plate motion in this area is on the order of 20 mm/y. Several moderate and large shocks have occurred along the SSF and LVF since 1640, including an Ms 7.6 event in 1900 on SSF. Although the region may be relatively far from a repeat of this earthquake, seismicity data indicate that strong shocks could take place along segments of the seismically active faults identified in this study.     

Long-term variation of solar UVB (290–330 nm) observed at the earth's surface

During solar cycle 21 (1976–86), the primary solar irradiance at 300 nm was steady during 1980–82 and thereafter decreased until 1986 by only 2–3%. The stratospheric ozone in middle latitudes had a QBO of 3–4% in this interval but the long-term ozone trend was less than 3% per decade, which could result in a UVB increase of only 5–6% per decade. Thus, the combined effect of changes in primary solar irradiance and ozone changes could be an increase of 5–6% in UVB, observed at ground during 1977–81 and a steady level during 1981–86. During 1976–86, the average cloudiness changed by less than 5% indicating UVB changes of 5% or less on this count. The aerosol level was almost constant during 1976–82 and increased abruptly in 1982 due to the E1 Chichon eruption and decayed slowly unitl 1986. Thus, due to aerosols only, the UVB was expected to be constant during 1976–82, to decrease sharply in 1982 and to recoup slowly thereafter.Measurements of clear-sky solar UVB at ground made at Jungfraujoch (Swiss Alps, 47°N, 8°E) during 1981–89 and at Rockville, USA (39°N, 77°W) were not comparable between themselves and did not follow the above expected patterns. Neither did the all-day R-B meter UVB measurements at Philadelphia, USA (40°N, 75°W) and Minneapolis, USA (45°N, 93°W). We suspect that some of these measurements are erroneous. This needs further detailed scrutiny.     

Multi-instrument observations of the electric and magnetic field structure of omega bands

High time resolution data from the CUTLASS Finland radar during the interval 01:30–03:30 UT on 11 May, 1998, are employed to characterise the ionospheric electric field due to a series of omega bands extending 5° in latitude at a resolution of 45 km in the meridional direction and 50 km in the azimuthal direction. E-region observations from the STARE Norway VHF radar operating at a resolution of 15 km over a comparable region are also incorporated. These data are combined with ground magnetometer observations from several stations. This allows the study of the ionospheric equivalent current signatures and height integrated ionospheric conductances associated with omega bands as they propagate through the field-of-view of the CUTLASS and STARE radars. The high-time resolution and multi-point nature of the observations leads to a refinement of the previous models of omega band structure. The omega bands observed during this interval have scale sizes 500 km and an eastward propagation velocity 0.75 km s^–1. They occur in the morning sector (05 MLT), simultaneously with the onset/intensification of a substorm to the west during the recovery phase of a previous substorm in the Scandinavian sector. A possible mechanism for omega band formation and their relationship to the substorm phase is discussed.     

Coseismic slip in the 1964 Prince William Sound earthquake: A new geodetic inversion

The 1964 Prince William Sound earthquake (March 28, 1964;M _w =9.2) caused crustal deformation over an area of approximately 140,000 km² in south central Alaska. In this study geodetic and geologic measurements of this surface deformation were inverted for the slip distribution on the 1964 rupture surface. Previous seismologic, geologic, and geodetic studies of this region were used to constrain the geometry of the fault surface. In the Kodiak Island region, 28 rectangular planes (50 by 50 km each) oriented 218°N, with a dip varying from 8^o nearest the Aleutian trench to 9^o below Kodiak Island, define the rupture surface. In the Prince William Sound region 39 planes with variable dimensions (40 by 50 km near the trench, 64 by 50 km inland) and orientation (218°N in the west and 270°N in the east) were used to approximate the complex faulting. Prior information was introduced to constrain offshore dip-slip values, the strike-slip component, and slip variation between adjacent planes. Our results suggest a variable dip-slip component with local slip maximums occurring near Montague Island (up to 30 m), further to the east near Kayak Island (up to 14 m), and trenchward of the northeast segment of Kodiak Island (up to 17m). A single fault plane dipping 30°NW, corresponding to the Patton Bay fault, with a slip value of 8 m modeled the localized but large uplift on Montague Island. The moment calculated on the basis of our geodetically derived slip model of 5.0×10²⁹ dyne cm is 30% less than the seismic moment of 7.5×10²⁹ dyne cm calculated from long-period surface waves (Kanamori, 1970) but is close to the seismic moment of 5.9×10²⁹ dyne cm obtained byKikuchi andFukao (1987).     

A simple formula for interpolation inn-dimensional space

n¶rt;m u u¶rt; umnu n-a nmam, ma ¶rt;mm n¶rt; u (auum m um ¶rt;uam, aum ua¶rt;m umnua uu). ¶rt;a nu¶rt;um nmu u¶rt; ma umnu u m ma.     

Extension of Antarctic ozone hole to lower latitudes in the South-American region

A comparison of monthly mean values of total ozone at South Pole, Buenos Aires (Argentina), Cachoeira Paulista and Natal (Brazil), and Huancayo (Peru) revealed that whereas South Pole showed an ozone depletion of 45% in October 1987 (as compared to October, 1977), Buenos Aires showed a small decrease (10%) while the other locations showed very small decreases (1–2%). When daily values are considered, the Antarctic ozone hole of October 1987 seems to have caused 10% depletion at Buenos Aires and 5% at Natal and Huancayo in December 1987. However, a large part of this is normal seasonal variation, except at Huancayo, where a residual effect of 5% depletion in December 1987 remains. The QBO effects (5–8% changes in the ozone level in 2–3 years) could cause 10–15% fluctuations in solar UVB on the ground on clear-sky days and could be a possible health hazard unless factors like cloudiness reduce the UVB intensities.     

Volcán Quizapu,Chilean Andes

Quizapu is a flank vent of the basalt-to-rhyodacite Holocene stratocone, Cerro Azul, and lies at the focus of a complex Quaternary volcanic field on the Andean volcanic front. The Quizapu vent originated in 1846 when 5 km³ of hornblende-dacite magma erupted effusively with little accompanying tephra. Between 1907 and 1932, phreatic and strombolian activity reamed out a deep crater, from which 4 km³ of dacite magma identical to that of 1846 fed the great plinian event of 10–11 April 1932. Although a total of >9 km³ of magma was thus released in 86 years, there is no discernible subsidence. As the pre-plinian crater was lined by massive lavas, 1932 enlargement was limited and the total plinian deposit contains only 0.4 wt % lithics. Areas of 5-cm and 1-cm isopachs for compacted 1932 fallout are about half of those estimated in the 1930's, yielding a revised ejecta volume of 9.5 km³. A strong inflection near the 10-cm isopach (downwind 110 km) on a plot of log Thickness vs Area^1/2 reflects slow settling of fine plinian ash — not of coignimbrite ash, as the volume of pyroclastic flows was trivial (<0.01 km³). About 17 vol.% of the fallout lies beyond the 1-cm isopach, and 82 wt% of the ejecta are finer than 1 mm. A least 18 hours of steady plinian activity produced an exceptionally uniform fall deposit. Observed column height (27–30 km) and average mass eruption rate (1.5x10⁸ kg/s) compare well with values for height and peak intensity calculated from published eruption models. The progressive aeolian fractionation of downwind ash (for which Quizapu is widely cited) is complicated by the large compositional range of 1932 juvenile pumice (52–70% SiO₂). The eruption began with andesitic scoria and ended with basaltic scoria, but >95% of the ejecta are dacitic pumice (67–68% SiO₂); minor andesitic scoria and frothier rhyodacite pumice (70% SiO₂) accompanied the dominant dacite. Phenocrysts (pl>hbopx>mt>ilmcpx) are similar in both abundance and composition in the 1846 (effusive) and 1932 (plinian) dacites. Despite the contrast in mode of eruption, bulk compositions are also indistinguishable. The only difference so far identified is a lower range of D values for 1846 hornblende, consistent with pre-eruptive degassing of the effusive batch.     

Evolution of the Olympus Mons Caldera,Mars

Synoptic images of the Martian volcano Olympus Mons are of a quality and quantity that are unique for mars and, somewhat surprisingly, are appreciably better than image data that exist for many volcanoes on Earth. Useful information about the evolution of shield volcanoes on Earth can thus be derived from the investigation of this extraterrestrial example. We have used shadow-length measurements and photoclinometrically derived profiles to supplement and refine the topographic map of the Olympus Mons caldera. As much as 2.5 km of collapse took place within the 80×65 km diameter caldera and the elevation of the caldera rim varies by almost 2.0 km (low around the oldest collapse events, high around the youngest). An eight-stage evolutionary sequence for the caldera of Olympus Mons is identified which shows that caldera subsidence was a longterm process rather than the near-instantaneous event that has been interpreted from comparable terrestrial examples. Tectonic features on the caldera floor indicate a transition from an extensional environment (graben formation) around the perimeter of the caldera to compression (ridge formation) towards the caldera center. This transition from a compressional to extensional environment is surprisingly sudden, occurs at a radial distance of 17 km from the caldera center, and is import because it can be used to infer that the magma chamber was relatively shallow (thought to be at a depth of <16 km beneath the caldera floor; Zuber and Mouginis-Mark 1990). Ample evidence is also found within the Olympus Mons caldera for solidified lava lakes more than 30 km in width, and for the localzed overturning and/or withdrawal of lava within these lakes.     

Coastal lava flows from Mauna Loa and Hualalai volcanoes,Kona, Hawaii

A major carbonate reef which drowned 13 ka is now submerged 150 m below sea level on the west coast of the island of Hawaii. A 25-km span of this reef was investigated using the submersibleMakali'i. The reef occurs on the flanks of two active volcanoes, Mauna Loa and Hualalai, and the lavas from both volcanoes both underlie and overlie the submerged reef. Most of the basaltic lava flows that crossed the reef did so when the water was much shallower, and when they had to flow a shorter distance from shoreline to reef face. Lava flows on top of the reef have protected it from erosion and solution and now occur at seaward-projecting salients on the reef face. These relations suggest that the reef has retreated shoreward as much as 50 m since it formed. A 7-km-wide shadow zone occurs where no Hualalai lava flows cross the reef south of Kailua. These lava flows were probably diverted around a large summit cone complex. A similar shadow zone on the flank of Mauna Loa volcano in the Kealakekua Bay region is downslope from the present Mauna Loa caldera, which ponds Mauna Loa lava and prevents it from reaching the coastline. South of the Mauna Loa shadow zone the - 150 m reef has been totally covered and obscured by Mauna Loa lava. The boundary between Hualalai and Mauna Loa lava on land occurs over a 6-km-wide zone, whereas flows crossing the - 150 m reef show a sharper boundary offshore from the north side of the subaerial transition zone. This indicates that since the formation of the reef, Hualalai lava has migrated south, mantling Mauna Loa lava. More recently, Mauna Loa lava is again encroaching north on Hualalai lava.     

Evolution of Icelandic central volcanoes: evidence from the Austurhorn intrusion,southeastern Iceland

The Austurhorn intrusive complex in southeastern Iceland represents an exhumed Tertiary central volcano. The geometry of the intrusion and geochemistry of the mafic and felsic rocks indicate Austurhorn was a volcanic center analogous to Eyjafjallajökull and Torfajökull in Iceland's eastern neovolcanic zone (EVZ). Early transitional tholeiitic basalt magmatism at Austurhorn formed a shallow crustal chamber 5 km in diameter. Apparent rhythmic modal layering of, and intrusive contacts within, the gabbro indicate the mafic chamber was replenished frequently as it cooled and crystallized. Felsic activity postdated near-solidification of the gabbro; numerous granitic magmas intruded along gabbro margins and within the adjacent crust. Field relations indicate that infrequent felsic replenishment prevented convective mixing of the Austurhorn chamber during this time, although commingled mafic and felsic magmas are observed in an extensive net veined complex. Late stage mafic dikes intrude the entire complex, suggesting that magmatic heat was abundantly available throughout the evolution of the Austurhorn system. Plagioclase and clinopyroxene compositions in mafic through felsic rocks, including gabbros, support a model of progressive differentiation. Field relations constrain the felsic magmas to originate at P1 kbar, presumably by fractional crystallization. The structure and geochemistry of the Austurhorn intrusive complex suggest formation in an immature rift environment similar to the modern EVZ. The proposed rift segment was parallel to the western and eastern neovolcanic zones, and probably resulted from a reorganization of plate boundaries 7 Ma (Saemundsson 1979; Helgason 1985; Jancin et al. 1985) triggered by activity of the Iceland mantle plume.     

Lengths and hazards from channel-fed lava flows on Mauna Loa, Hawai‘i, determined from thermal and downslope modeling with FLOWGO

Using the FLOWGO thermo-rheological model we have determined cooling-limited lengths of channel-fed (i.e. a) lava flows from Mauna Loa. We set up the program to run autonomously, starting lava flows from every 4th line and sample in a 30-m spatial-resolution SRTM DEM within regions corresponding to the NE and SW rift zones and the N flank of the volcano. We consider that each model run represents an effective effusion rate, which for an actual flow coincides with it reaching 90% of its total length. We ran the model at effective effusion rates ranging from 1 to 1,000 m³ s^–1, and determined the cooling-limited channel length for each. Keeping in mind that most flows extend 1–2 km beyond the end of their well-developed channels and that our results are non-probabilistic in that they give all potential vent sites an equal likelihood to erupt, lava coverage results include the following: SW rift zone flows threaten almost all of Mauna Loas SW flanks, even at effective effusion rates as low as 50 m³ s^–1 (the average effective effusion rate for SW rift zone eruptions since 1843 is close to 400 m³ s^–1). N flank eruptions, although rare in the recent geologic record, have the potential to threaten much of the coastline S of Keauhou with effective effusion rates of 50–100 m³ s^–1, and the coast near Anaehoomalu if effective effusion rates are 400–500 m³ s^–1 (the 1859 a flow reached this coast with an effective effusion rate of 400 m³ s^–1). If the NE rift zone continues to be active only at elevations >2,500 m, in order for a channel-fed flow to reach Hilo the effective effusion rate needs to be 400 m³ s^–1 (the 1984 flow by comparison, had an effective effusion rate of 200 m³ s^–1). Hilo could be threatened by NE rift zone channel-fed flows with lower effective effusion rates but only if they issue from vents at 2,000 m or lower. Populated areas on Mauna Loas SE flanks (e.g. Phala), could be threatened by SW rift zone eruptions with effective effusion rates of 100 m³ s^–1.Editorial responsibility: J Donnelly-Nolan     

Thermodynamics of gas and steam-blast eruptions

Eruptions of gas or steam and non-juvenile debris are common in volcanic and hydrothermal areas. From reports of non-juvenile eruptions or eruptive sequences world-wide, at least three types (or end-members) can be identified: (1) those involving rock and liquid water initially at boiling-point temperatures (boiling-point eruptions); (2) those powered by gas (primarily water vapor) at initial temperatures approaching magmatic (gas eruptions); and (3) those caused by rapid mixing of hot rock and ground- or surface water (mixing eruptions). For these eruption types, the mechanical energy released, final temperatures, liquid water contents and maximum theoretical velocities are compared by assuming that the erupting mixtures of rock and fluid thermally equilibrate, then decompress isentropically from initial, near-surface pressure (10 MPa) to atmospheric pressure. Maximum mechanical energy release is by far greatest for gas eruptions (1.3 MJ/kg of fluid-rock mixture)-about one-half that of an equivalent mass of gunpowder and one-fourth that of TNT. It is somewhat less for mixing eruptions (0.4 MJ/kg), and least for boiling-point eruptions (0.25 MJ/kg). The final water contents of crupted boiling-point mixtures are usually high, producing wet, sloppy deposits. Final erupted mixtures from gas eruptions are nearly always dry, whereas those from mixing eruptions vary from wet to dry. If all the enthalpy released in the eruptions were converted to kinetic energy, the final velocity (v _max) of these mixtures could range up to 670 m/s for boiling-point eruptions and 1820 m/s for gas eruptions (highest for high initial pressure and mass fractions of rock (m _r) near zero). For mixing eruptions, v _max ranges up to 1150 m/s. All observed eruption velocities are less than 400 m/s, largely because (1) most solid material is expelled when m _r is high, hence v _max is low; (2) observations are made of large blocks the velocities of which may be less than the average for the mixture; (3) heat from solid particles is not efficiently transferred to the fluid during the eruptions; and (4) maximum velocities are reduced by choked flow or friction in the conduit.