Tritium and stable isotopes of magmatic waters

To investigate the isotopic composition and age of water in volcanic gases and magmas, we analyzed samples from 11 active volcanoes ranging in composition from tholeiitic basalt to rhyolite: Mount St. Helens (USA), Kilauea (USA), Pacaya (Guatemala), Galeras (Colombia), Satsuma Iwo-Jima (Japan), Sierra Negra and Alcedo (Ecuador), Vulcano (Italy), Parícutin (Mexico), Kudryavy (Russia), and White Island (New Zealand). Tritium at relatively low levels (0.1–5 T.U.) is found in most emissions from high-temperature volcanic fumaroles sampled, even at discharge temperatures >700°C. Although magmatic fluids sampled from these emissions usually contain high CO₂, S_total, HCl, HF, B, Br, ³He R/R_A, and low contents of air components, stable isotope and tritium relations of nearly all such fluids show mixing of magmatic volatiles with relatively young meteoric water (model ages≤75 y). Linear δD/δ¹⁸O and ³H/δ¹⁸O mixing trends of these two end-members are invariably detected at arc volcanoes. Tritium is also detected in fumarole condensates at hot spot basalt volcanoes, but collecting samples approaching the composition of end-member magmatic fluid is exceedingly difficult. In situ production of ³H, mostly from spontaneous fission of ²³⁸U in magmas is calculated to be <0.001 T.U., except for the most evolved compositions (high U, Th, and Li and low H₂O contents). These values are below the detection limit of ³H by conventional analytical techniques (about 0.01 T.U. at best). We found no conclusive evidence that natural fusion in the Earth produces anomalous amounts of detectable ³H (>0.05 T.U.).     

Presence of native gold and tellurium in the active high-sulfidation hydrothermal system of the La Fossa volcano (Vulcano, Italy)

Grains of native gold and tellurium were found in siliceous hydrothermally altered rocks in the high-temperature (170–540°C) fumarolic field of the La Fossa volcano (Island of Vulcano). In addition to Au and Te, Pb–Bi sulfides (cannizzarite) and Tl-bromide chloride were found as sublimates in the hottest fumarolic vents of the crater rim. The chemical composition of altered rocks associated with sublimate deposition indicate the presence of a significant concentration of Te (up to 75 ppm), while gold concentrations are very low (<9 ppb). Pb, Bi and Tl are strongly enriched in the hottest and less oxidized fumarolic vents, reaching concentrations of 2186, 146 and 282 ppm, respectively. These elements are transported (generally as chloride complexes) to the surface by volcanic gases, and several of these (Bi, Te, Tl) are originated from magma degassing. The silicic alteration is produced by the flow of fluids with pH<2. High acidity results from introduction of magmatic gases such as SO₂, HCl and HF released by the shallow magmatic reservoir of La Fossa volcano. The silicic alteration found at Vulcano may represent an early stage of the `vuggy silica' facies which characterizes the high-sulfidation epithermal ore deposits, confirming the analogies existing between this type of ore deposit and magmatic-hydrothermal systems associated with island-arc volcanoes.     

Fumaroles in ice caves on the summit of Mount Rainier—preliminary stable isotope, gas, and geochemical studies

The edifice of Mount Rainier, an active stratovolcano, has episodically collapsed leading to major debris flows. The largest debris flows are related to argillically altered rock which leave areas of the edifice prone to failure. The argillic alteration results from the neutralization of acidic magmatic gases that condense in a meteoric water hydrothermal system fed by the melting of a thick mantle of glacial ice. Two craters atop a 2000-year-old cone on the summit of the volcano contain the world's largest volcanic ice-cave system. In the spring of 1997 two active fumaroles (T=62°C) in the caves were sampled for stable isotopic, gas, and geochemical studies.Stable isotope data on fumarole condensates show significant excess deuterium with calculated δD and δ¹⁸O values (−234 and −33.2‰, respectively) for the vapor that are consistent with an origin as secondary steam from a shallow water table which has been heated by underlying magmatic–hydrothermal steam. Between 1982 and 1997, δD of the fumarole vapor may have decreased by 30‰.The compositions of fumarole gases vary in time and space but typically consist of air components slightly modified by their solubilities in water and additions of CO₂ and CH₄. The elevated CO₂ contents (δ¹³C_CO2=−11.8±0.7‰), with spikes of over 10,000 ppm, require the episodic addition of magmatic components into the underlying hydrothermal system. Although only traces of H₂S were detected in the fumaroles, most notably in a sample which had an air δ¹³C_CO2 signature (−8.8‰), incrustations around a dormant vent containing small amounts of acid sulfate minerals (natroalunite, minamiite, and woodhouseite) indicate higher H₂S (or possibly SO₂) concentrations in past fumarolic gases.Condensate samples from fumaroles are very dilute, slightly acidic, and enriched in elements observed in the much higher temperature fumaroles at Mount St. Helens (K and Na up to the ppm level; metals such as Al, Pb, Zn Fe and Mn up to the ppb level and volatiles such as Cl, S, and F up to the ppb level).The data indicate that the hydrothermal system in the edifice at Mount Rainier consists of meteoric water reservoirs, which receive gas and steam from an underlying magmatic system. At present the magmatic system is largely flooded by the meteoric water system. However, magmatic components have episodically vented at the surface as witnessed by the mineralogy of incrustations around inactive vents and gas compositions in the active fumaroles. The composition of fumarole gases during magmatic degassing is distinct and, if sustained, could be lethal. The extent to which hydrothermal alteration is currently occurring at depth, and its possible influence on future edifice collapse, may be determined with the aid of on site analyses of fumarole gases and seismic monitoring in the ice caves.     

Evaluation of gases,condensates, and SO2 emissions from Augustine volcano,Alaska: the degassing of a Cl-rich volcanic system

After the March–April 1986 explosive eruption a comprehensive gas study at Augustine was undertaken in the summers of 1986 and 1987. Airborne COSPEC measurements indicate that passive SO₂ emission rates declined exponentially during this period from 380±45 metric tons/day (T/D) on 7/24/86 to 27±6 T/D on 8/24/87. These data are consistent with the hypothesis that the Augustine magma reservoir has become more degassed as volcanic activity decreased after the spring 1986 eruption. Gas samples collected in 1987 from an 870°C fumarole on the andesitic lava dome show various degrees of disequilibrium due to oxidation of reduced gas species and condensation (and loss) of H₂O in the intake tube of the sampling apparatus. Thermochemical restoration of the data permits removal of these effects to infer an equilibrium composition of the gases. Although not conclusive, this restoration is consistent with the idea that the gases were in equilibrium at 870°C with an oxygen fugacity near the Ni–NiO buffer. These restored gas compositions show that, relative to other convergent plate volcanoes, the Augustine gases are very HCl rich (5.3–6.0 mol% HCl), S rich (7.1 mol% total S), and H₂O poor (83.9–84.8 mol% H₂O). Values of D and ¹⁸O suggest that the H₂O in the dome gases is a mixture of primary magmatic water (PMW) and local seawater. Part of the Cl in the Augustine volcanic gases probably comes from this shallow seawater source. Additional Cl may come from subducted oceanic crust because data by Johnston (1978) show that Cl-rich glass inclusions in olivine crystals contain hornblende, which is evidence for a deep source (>25km) for part of the Cl. Gas samples collected in 1986 from 390°–642°C fumaroles on a ramp surrounding the inner summit crater have been oxidized so severely that restoration to an equilibrium composition is not possible. H and O isotope data suggest that these gases are variable mixtures of seawater, FMW, and meteoric steam. These samples are much more H₂O-rich (92%–97% H₂O) than the dome gases, possibly due to a larger meteoric steam component. The 1986 samples also have higher Cl/S, S/C, and F/Cl ratios, which imply that the magmatic component in these gases is from the more degassed 1976 magma. Thus, the 1987 samples from the lava dome are better indicators than the 1986 samples of degassing within the Augustine magma reservoir, even though they were collected a year later and contain a significant seawater component. Future gas studies at Augustine should emphasize fumaroles on active lava domes. Condensates collected from the same lava-dome fumarole have enrichments ot 10⁷–10² in Cl, Br, F, B, Cd, As, S, Bi, Pb, Sb, Mo, Zn, Cu, K, Li, Na, Si, and Ni. Lower-temperature (200°–650°C) fumaroles around the volcano are generally less enriched in highly volatile elements. However, these lower-termperature fumaroles have higher concentration of rock-forming elements, probably derived from the wall rock.     

Caldera formation at Volcán Colima, Mexico, by a large large holocene volcanic debris avalanche

About 4,300 years ago, 10 km³ of the upper cone of ancestral Volcán Colima collapsed to the southwest leaving a horseshoe-shaped caldera 4 km in diameter. The collapse produced a massive volcanic debris avalanche deposit covering over 1550 km² on the southern flanks of the volcano and extending at least 70 km from the former summit. The avalanche followed a steep topographic gradient unobstructed by barriers, resulting in an unusually high area/volume ratio for the Colima deposit. The apparent coefficient of friction (fall height/distance traveled) for the Colima avalanche is 0.06, a low value similar to those of other large-volume deposits. The debris avalanche deposit contains 40–75% angular volcanic clasts from the ancestral cone, a small proportion of vesicular blocks that may be juvenile, and in distal exposures, rare carbonate clasts plucked from the underlying surface by the moving avalanche. Clasts range in size to over 20 m in diameter and are brecciated to different degrees, pulverized, and surrounded by a rock-flour matrix. The upper surface of the deposit shows prominent hummocky topography with closed depressions and surface boulders. A thick, coarse-grained, compositionally zoned scoria-fall layer on the upper northeastern slope of the volcano may have erupted at the time of collapse. A fine-grained surge layer is present beneath the avalanche deposit at one locality, apparently representing an initial blast event. Most of the missing volume of the ancestral volcano has since been restored at an average rate of 0.002 km³/yr through repeated eruptions from the post-caldera cone. As a result, the southern slope of Volcán Colima may again be susceptible to collapse. Over 200,000 people are now living on primary or secondary deposits of the debris avalanche, and a repetition of this event would constitute a volcanic disaster of great magnitude.Ancestral Volcán Colima grew on the southern, trenchward flank of the earlier and larger volcano Nevado de Colima. Trenchward collapse was favored by the buttressing effect of Nevado, the rapid elevation drop to the south, and the intrusion of magma into the southern flank of the ancestral volcano. Other such trenchward-younging, paired volcanoes are known from Mexico, Guatemala, El Salvador, Chile, and Japan. The trenchward slopes of the younger cones are common sites for cone collapse to form avalanche deposits, as occurred at Colima and Popocatepetl in Mexico and at San Pedro Volcano in Chile.     

Fumarolic emissions from Mount St. Augustine,Alaska: 1979–1984 degassing trends,volatile sources and their possible role in eruptive style

Gas samples were collected from high-temperature, rooted summit vents at Mount St. Augustine in 1979, 1982, and 1984. All of the gas samples exhibit various degrees of disequilibrium. Thermodynamic restoration of the analyzed gases permits partial or complete removal of these disequilibrium effects and allows inference of equilibrium gas compositions. Long-term (1979–1984) degassing trends within resampled or adjacent vents are characterized by increases (from 97.4 to 99.8 mole%) in the H₂O fraction and major decreases in the residual gases. Over this same period total gas HCl contents decreased by a factor of 3 to 10 while dry gas (H₂O-free recalculated) HCl contents increased by a factor of 1.6 to 3. Dry gas mole proportions at these sites changed from being CO₂-dominated (46% CO₂, 24% H₂ in 1979) to H₂-dominated (49% H₂, 22% CO₂ in 1984). The overall trends in gas chemistry and the stable isotope patterns in gases and condensates from the summit fumaroles can be explained by progressive magmatic outgassing coupled with increasing proportions of seawater in the fumarole emissions.Studies of the gaseous emissions following the 1976 and 1986 Mount St. Augustine eruptions confirmed the Cl- and S-rich nature of the Mount St. Augustine emanations. Seawater, possibly derived from magmatic assimilation or dehydration of near-surface seawater-bearing sediments, could supply a portion of the outgassed Cl and S. Continued seawater influx through subvolcanic fractures or permeable sediments would recharge the seawater-depleted zone and provide a near-surface Cl and S source for the next eruptive cycle,Various lines of evidence support a phreatomagmatic component in the 1976 and 1986 Mount St. Augustine eruptions. We suggest that seawater may interact with magma or volcanic gases during the early explosive phase of Mount St. Augustine eruptions and that it continues to influence high-temperature fumarole emissions as the volcanic system cools.     

Temporal geochemical variations in volatile emissions from Mount St. Helens, USA, 1980–1994

Fumarole discharges (95–560°C) collected from the dacite dome inside Mount St. Helens crater show temporal changes in their isotopic and chemical compositions. A δD vs. δ¹⁸O plot shows that condensed waters from the gases are mixtures of meteoric and magmatic components, but that the apparent magmatic end-member in 1994 was depleted by about 7‰ in δD relative to the apparent end-member in 1980. Based on δD modeling, approximately 63% of shallow, post-1980 magma has yet to degas. Surprisingly, Cl and F contents in the 1994 samples were only 0.47 and 3.8%, respectively, of the concentrations determined for end-member magmatic fluid in 1980. The data indicate that Cl (and F and B) is degassed from magma relatively quickly compared to water and/or that most of the Cl degassed in later years is dissolved into the shallow Mount St. Helens hydrothermal system. Because metals are often transported in magmatic and hydrothermal fluids as Cl complexes, rapid changes in surface volatile compositions may have implications for the timing and location of metals transport and deposition in some volcanoes.     

Physical-chemical conditions of the Teide volcanic system (Tenerife, Canary Islands)

The Teide volcano (3717 m) is the central structure of the island of Tenerife and at present its morphology is that of a stratovolcano which has grown on a large caldera with a collapse 17 km in diameter, which was generated some 0.6 million years ago.The different studies that have been carried out seem to indicate that, in a oversimplified model, there is an intermediate magma chamber with an approximate volume of 30 km³ and located 2–3 km below the actual base of the caldera, i.e., almost at sea level, with a temperature of 430 ± 50°C, and a pressure of 400 ± 100 bar.The summit fumarole emissions are 85°C and are formed mainly of CO₂ with small amounts of sulphur species, H₂, CH₄ and He. The water vapor (68–82%) emitted with the gases comes from the vaporization of a perched aquifer in the upper cone, as shown by the isotopic analyses.     

The stable isotope geochemistry of volcanic lakes, with examples from Indonesia

Stable isotope compositions (δD, δ¹⁸O and δ³⁴S) of volcanic lake waters, gas condensates and spring waters from Indonesia, Italy, Japan, and Russia were measured. The spring fluids and gas samples plot in a broad array between meteoric waters and local high-temperature volcanic gas compositions. The δD and δ¹⁸O data from volcanic lakes in East Indonesia plot in a concave band ranging from local meteoric waters to evaporated fluids to waters heavier than local high-temperature volcanic gases. We investigated isotopic fractionation processes in volcanic lakes at elevated temperatures with simultaneous mixing of meteoric waters and volcanic gases. An elevated lake water temperature gives enhanced kinetic isotope fractionation and changes in equilibrium fractionation factors, providing relatively flat isotope evolution curves in δ¹⁸O–δD diagrams. A numerical simulation model is used to derive the timescales of isotopic evolution of crater lakes as a function of atmospheric parameters, lake water temperature and fluxes of meteoric water, volcanic gas input, evaporation, and seepage losses. The same model is used to derive the flux magnitude of the Keli Mutu lakes in Indonesia. The calculated volcanic gas fluxes are of the same order as those derived from energy budget models or direct gas flux measurements in open craters (several 100 m³ volcanic water/day) and indicate a water residence time of 1–2 decades. The δ³⁴S data from the Keli Mutu lakes show a much wider range than those from gases and springs, which is probably related to the precipitation of sulfur in these acid brine lakes. The isotopic mass balance and S/Cl values suggest that about half of the sulfur input in the hottest Keli Mutu lake is converted into native sulfur.     

Seismicity and magma supply rate of the 1998 failed eruption at Iwate volcano,Japan

Iwate volcano, Japan, showed significant volcanic activity including earthquake swarms and volcano inflation from the beginning of 1998. A large earthquake of magnitude 6.1 hit the south-west of the volcano on September 3. Although a 1 km² fumarole field formed, blighting plants on the ridge in the western part of the volcano in the spring of 1999, no magmatic eruptions occurred. We reconcile the spatio-temporal distributions of volcanic pressure sources determined by previously reported studies in which GPS, strain and tilt data from dense geodetic station networks are analyzed (Miura et al. Earth Planet Space 52:1003–1008, 2000; Sato and Hamaguchi J Volcanol Geotherm Res 155:244–262, 2006). We calculate the magma supply rates from their results and compare them with the occurrence rates of volcanic earthquakes. The results show that the magma supply rates are almost constant or even decrease with time while the earthquake occurrence rate increases with time. This contrast in their temporal changes is interpreted to result from stress accumulation in the volcanic edifice caused by constant magma supply without effusion of magma to the surface. We further show that data showing slight acceleration in strain can be best explained by magma ascent at a constant velocity, and that there is no evidence for increased magma buoyancy resulting from gas bubble growth. This consideration supports the interpretation that the magma stayed at 2 km depth and horizontally migrated. These findings relating magma supply rate and seismicity to magma ascent process are clues to understanding why no magmatic eruption occurred at Iwate volcano in 1998.     

Hydrogen and oxygen isotopic compositions of waters from fumaroles at Kilauea summit,Hawaii

Condensate samples were collected in 1992 from a high-temperature (300° C) fumarole on the floor of the Halemaumau Pit Crater at Kilauea. The emergence about two years earlier of such a hot fumarole was unprecedented at such a central location at Kilauea. The condensates have hydrogen and oxygen isotopic compositions which indicate that the waters emitted by the fumarole are composed largely of meteoric water, that any magmatic water component must be minor, and that the precipitation that was the original source to the fumarole fell on a recharge area on the slopes of Mauna Loa Volcano to the west. However, the fumarole has no tritium, indicating that it taps a source of water that has been isolated from atmospheric water for at least 40 years. It is noteworthy, considering the unstable tectonic environment and abundant local rainfall of the Kilauea and Mauna Loa regions, that waters which are sources to the hot fumarole remain uncontaminated from atmospheric sources over such long times and long transport distances. As for the common, boiling point fumaroles of the Kilauea summit region, their ¹⁸O, D and tritium concentrations indicate that they are dominated by recycling of present day meteoric water. Though the waters of both hot and boiling point fumaroles have dominantly meteoric sources, they seem to be from separate hydrological regimes. Large concentrations of halogens and sulfur species in the condensates, together with the location at the center of the Kilauea summit region and the high temperature, initially suggested that much of the total mass of the emissions of the hot fumarole, including the H₂O, might have come directly from a magma body. The results of the present study indicate that it is unreliable to infer a magmatic origin of volcanic waters based solely on halogen or sulfur contents, or other aspects of chemical composition of total condensates.     

Geochemical mapping of magmatic gas–water–rock interactions in the aquifer of Mount Etna volcano

Systematic analysis of major and minor elements in groundwaters from springs and wells on the slopes of Mt. Etna in 1995–1998 provides a detailed geochemical mapping of the aquifer of the volcano and of the interactions between magmatic gas, water bodies and their host rocks. Strong spatial correlations between the largest anomalies in pCO₂ (pH and alkalinity) K, Rb, Mg, Ca and Sr suggest a dominating control by magmatic gas (CO₂) and consequent basalt leaching by acidified waters of the shallow (meteoric) Etnean aquifer. Most groundwaters displaying this magmatic-type interaction discharge within active faulted zones on the S–SW and E lower flanks of the volcanic pile, but also in a newly recognised area on the northern flank, possibly tracking a main N–S volcano-tectonic structure. In the same time, the spatial distribution of T°C, TDS, Na, Li, Cl and B allows us to identify the existence of a deeper thermal brine with high salinity, high content of B, Cl and gases (CO₂, H₂S, CH₄) and low K/Na ratio, which is likely hosted in the sedimentary basement. This hot brine reaches the surface only at the periphery of the volcano near the Village of Paternò, where it gives rise to mud volcanoes called “Salinelle di Paternò”. However, the contribution of similar brines to shallower groundwaters is also detected in other sectors to the W (Bronte, Maletto), SW (Adrano) and SE (Acireale), suggesting its possible widespread occurrence beneath Etna. This thermal brine is also closely associated with hydrocarbon fields all around the volcano and its rise, generally masked by the high outflow of the shallow aquifer, may be driven by the ascent of mixed sedimentary–magmatic gases through the main faults cutting the sedimentary basement.     

Correlation of seismic activity and fumarole temperature at the Mt. Merapi volcano (Indonesia) in 2000

In this paper we present densely sampled fumarole temperature data, recorded continuously at a high-temperature fumarole of Mt. Merapi volcano (Indonesia). These temperature time series are correlated with continuous records of rainfall and seismic waveform data collected at the Indonesian–German multi-parameter monitoring network. The correlation analysis of fumarole temperature and precipitation data shows a clear influence of tropical rain events on fumarole temperature. In addition, there is some evidence that rainfall may influence seismicity rates, indicating interaction of meteoric water with the volcanic system. Knowledge about such interactions is important, as lava dome instabilities caused by heavy-precipitation events may result in pyroclastic flows. Apart from the strong external influences on fumarole temperature and seismicity rate, which may conceal smaller signals caused by volcanic degassing processes, the analysis of fumarole temperature and seismic data indicates a statistically significant correlation between a certain type of seismic activity and an increase in fumarole temperature. This certain type of seismic activity consists of a seismic cluster of several high-frequency transients and an ultra-long-period signal (<0.002 Hz), which are best observed using a broadband seismometer deployed at a distance of 600 m from the active lava dome. The corresponding change in fumarole temperature starts a few minutes after the ultra-long-period signal and simultaneously with the high-frequency seismic cluster. The change in fumarole temperature, an increase of 5 °C on average, resembles a smoothed step. Fifty-four occurrences of simultaneous high-frequency seismic cluster, ultra-long period signal and increase of fumarole temperature have been identified in the data set from August 2000 to January 2001. The observed signals appear to correspond to degassing processes in the summit region of Mt. Merapi.     

Magnetotelluric study of the Teide (Tenerife) and Timanfaya (Tanzarote) volcanic areas

The application of a magnetotelluric surveying method, developed specifically for detection of strong conductivity contrasts, has permitted interpretations of the deep structure of two important volcanic formations in the Canary Archipelago: the Timanfaya (Lanzarote) and Teide (Tenerife) volcanoes. Beneath the Timanfaya volcano, a highly conductive body (p ≤ 1 Ω m) has been detected, which is interpreted as a shallow magma chamber associated with the eruption which took place between the years 1730 and 1736. The core of this intrusion, still largely molten, would be located at approximately 4 km depth and would still maintain temperatures in excess of 900 °C.In the Telde caldera the asymmetry of the apparent resistivity components indicates a double depression; furthermore, a highly conductive level encountered in the centre of the eastern part of the volcanic complex has been attributed to the persistance of relatively high temperatures in the root of the last salic dome, extruded inside the caldera. In both zones, the deep limiting horizon is a highly conductive layer corresponding to the mantle (15 km in Lanzarote and 13 km in Tenerife), whereas other, more superficial, discontinuities have been related to local structures.     

Emplacement and arrest of sheets and dykes in central volcanoes

Sheet intrusions are of two main types: local inclined (cone) sheets and regional dykes. In Iceland, the inclined sheets form dense swarms of (mostly) basaltic, 0.5–1 m thick sheets, dipping either at 20–50° or at 75–90° towards the central volcano to which they belong. The regional dykes are (mostly) basaltic, 4–6 m thick, subvertical, subparallel and form swarms, less dense than those of the sheets but tens of kilometres long, in the parts of the volcanic systems that are outside the central volcanoes. In both types of swarms, the intrusion intensity decreases with altitude in the lava pile. Theoretical models generally indicate very high crack-tip stresses for propagating dykes and sheets. Nevertheless, most of these intrusions become arrested at various crustal depths and never reach the surface to supply magma to volcanic eruptions. Two principal mechanisms are proposed to explain arrest of dykes and sheets. One is the generation of stress barriers, that is, layers with local stresses unfavourable for the intrusion propagation. The other is mechanical anisotropy whereby sheet intrusions become arrested at discontinuities. Stress barriers may develop in several ways. First, analytical solutions for a homogeneous and isotropic crust show that the intensity of the tensile stress associated with a pressured magma chamber falls off rapidly with distance from the chamber. Thus, while dyke and sheet injection in the vicinity of a chamber may be favoured, dyke and sheet arrest is encouraged in layers (stress barriers) at a certain distance from the chamber. Second, boundary-element models for magma chambers in a mechanically layered crust indicate abrupt changes in tensile stresses between layers of contrasting Young’s moduli (stiffnesses). Thus, where soft pyroclastic layers alternate with stiff lava flows, as in many volcanoes, sheet and dyke arrest is encouraged. Abrupt changes in stiffness between layers are commonly associated with weak and partly open contacts and other discontinuities. It follows that stress barriers and discontinuities commonly operate together as mechanisms of dyke and sheet arrest in central volcanoes.     

Microgravity and height changes caused by volcanic activity: Four

Microgravity measurements and levelling surveys on volcanoes are not always easy to make, but are useful for studying volcanic processes quantitatively. Gravity changes associated with volcanic activity are not always significant. Precision of microgravity measurements depend critically on the procedures adopted, and those applied in the present paper are described. Levelling technique is now orthodox, and some empirical laws relating ground deformation to volcanic activity are deduced from the accumulated data. Gravity changes occur at the same time and places as ground deformations. The relationship between microgravity and height changes are discussed from the standpoint of analyzing the data obtained on volcanoes. The observational results obtained on four volcanoes in Japan are separately analyzed because each volcano exhibits different patterns of gravity changes and deformations. During the 1977–1982 activity of Usu volcano, deformation was accompanied by microgravity changes frequently observed at a particular benchmark at the base of the volcano for about five years. The gravity changes prove to be not a direct effect of magma movements but to be caused by the deformations of ground strata and aquifers around the benchmark. The 1983 eruption of Miyakejima volcano was associated with local gravity changes around the eruptive fissures due to magma intrusion which was approximately modelled. Similarly the 1986 eruption of Ooshima volcano caused gravity changes on the volcano, but these were poorly correlated with elevation changes and their origins were not uniquely interpreted. To detect gravity changes associated with the activity of Sakurajima volcano, an equigravity point was selected at the north of the volcano besides the gravity points on and around the volcano itself. The probable gradual accumulation of magmas beneath the volcano for eight years is substantiated by observed microgravity and elevation changes.     

Geology and eruptive history of Unzen volcano, Shimabara Peninsula, Kyushu, SW Japan

During the past 500 thousand years, Unzen volcano, an active composite volcano in the Southwest Japan Arc, has erupted lavas and pyroclastic materials of andesite to dacite composition and has developed a volcanotectonic graben. The volcano can be divided into the Older and the Younger Unzen volcanoes. The exposed rocks of the Older Unzen volcano are composed of thick lava flows and pyroclastic deposits dated around 200–300 ka. Drill cores recovered from the basal part of the Older Unzen volcano are dated at 400–500 ka. The volcanic rocks of the Older Unzen exceed 120 km³ in volume. The Younger Unzen volcano is composed of lava domes and pyroclastic deposits, mostly younger than 100 ka. This younger volcanic edifice comprises Nodake, Myokendake, Fugendake, and Mayuyama volcanoes. Nodake, Myokendake and Fugendake volcanoes are 100–70 ka, 30–20 ka, and <20 ka, respectively. Mayuyama volcano formed huge lava domes on the eastern flank of the Unzen composite volcano about 4000 years ago. Total eruptive volume of the Younger Unzen volcano is about 8 km³, and the eruptive production rate is one order of magnitude smaller than that of the Older Unzen volcano.     

Cladistic analysis applied to the classification of volcanoes

Cladistics is a systematic method of classification that groups entities on the basis of sharing similar characteristics in the most parsimonious manner. Here cladistics is applied to the classification of volcanoes using a dataset of 59 Quaternary volcanoes and 129 volcanic edifices of the Tohoku region, Northeast Japan. Volcano and edifice characteristics recorded in the database include attributes of volcano size, chemical composition, dominant eruptive products, volcano morphology, dominant landforms, volcano age and eruptive history. Without characteristics related to time the volcanic edifices divide into two groups, with characters related to volcano size, dominant composition and edifice morphology being the most diagnostic. Analysis including time based characteristics yields four groups with a good correlation between these groups and the two groups from the analysis without time for 108 out of 129 volcanic edifices. Thus when characters are slightly changed the volcanoes still form similar groupings. Analysis of the volcanoes both with and without time yields three groups based on compositional, eruptive products and morphological characters. Spatial clusters of volcanic centres have been recognised in the Tohoku region by Tamura et al. (Earth Planet Sci Lett 197:105–106, 2002). The groups identified by cladistic analysis are distributed unevenly between the clusters, indicating a tendency for individual clusters to form similar kinds of volcanoes with distinctive but coherent styles of volcanism. Uneven distribution of volcano types between clusters can be explained by variations in dominant magma compositions through time, which are reflected in eruption products and volcanic landforms. Cladistic analysis can be a useful tool for elucidating dynamic igneous processes that could be applied to other regions and globally. Our exploratory study indicates that cladistics has promise as a method for classifying volcanoes and potentially elucidating dynamic and evolutionary volcanic processes. Cladistics may also have utility in hazards assessment where spatial distributions and robust definitions of a volcano are important, as in locating sensitive facilities such as nuclear reactors and repositories.     

Cl,F, and SO2 in Central American volcanic gases

The Cl, F, SO₂ contents and Cl/F and Cl/SO₄ ratios in Central American volcanic gases are examined. 10³?10⁵ tons SO₂ per day are given off during eruptions and 10² when a prominent vapor cloud persists between eruptive periods. Data regarding Cl and F and SO₄ from leachates, condensates, and incrustations are compared. Our data suggest circumpacific volcanoes are SO₂ poor relative to Cl and may be F poor although F is higher in basaltic Central American volcanoes than others.     

Stress Field Estimations for Colima Volcano, Mexico, Based on Seismic Data

 For first time, during 1991, seismic activity was recorded during an eruption at Colima volcano. We analyze these data to obtain a stress pattern using a composite focal mechanism technique. From the analysis of regional seismicity, the Tamazula Fault and the Armeria River appear as active features and the dip of the slab east of the Jalisco Block is approximately 12°. Southwest of Colima volcano a vertical alignment of seismic events was observed. We estimate five different composite focal mechanism solutions from our data set, which indicate a change of the stress field at the volcano after the 1991 eruption. These solutions suggest that the stress field in the volcanic edifice was controlled by stresses related to the emplacement of magma superimposed on the regional stress field. No evidence of active local faults in the volcanic edifice was found. We propose a model for the eruptive process that involves tilting of the volcanic edifice. Received: 15 October 1995 / Accepted: 26 October 1998