Bubble growth in rhyolitic melts: experimental and numerical investigation

 Bubble growth controlled by mass transfer of water from hydrated rhyolitic melts at high pressures and temperatures was studied experimentally and simulated numerically. Rhyolitic melts were hydrated at 150 MPa, 780–850  °C to uniform water content of 5.5–5.3 wt%. The pressure was then dropped and held constant at 15–145 MPa. Upon the drop bubbles nucleated and were allowed to grow for various periods of time before final, rapid quenching of the samples. The size and number density of bubbles in the quenched glasses were recorded. Where number densities were low and run duration short, bubble sizes were in accord with the growth model of Scriven (1959) for solitary bubbles. However, most results did not fit this simple model because of interaction between neighboring bubbles. Hence, the growth model of Proussevitch et al. (1993), which accounts for finite separation between bubbles, was further developed and used to simulate bubble growth. The good agreement between experimental data, numerical simulation, and analytical solutions enables accurate and reliable examination of bubble growth from a limited volume of supersaturated melt. At modest supersaturations bubble growth in hydrated silicic melts (3–6 wt% water, viscosity 10⁴–10⁶ Pa·s) is diffusion controlled. Water diffusion is fast enough to maintain steady-state concentration gradient in the melt. Viscous resistance is important only at the very early stage of growth (t<1 s). Under the above conditions growth is nearly parabolic, R²=2Dtρ_m(C₀–C_f)/ρ_g until the bubble approaches its final size. In melts with low water content, viscosity is higher and maintains pressure gradients in the melt. Growth may be delayed for longer times, comparable to time scales of melt ascent during eruptions. At high levels of supersaturation, advection of hydrated melt towards the growing bubble becomes significant. Our results indicate that equilibrium degassing is a good approximation for modeling vesiculation in melts with high water concentrations (C₀>3 wt%) in the region above the nucleation level. When the melt accelerates and water content decreases, equilibrium can no longer be maintained between bubbles and melt. Supersaturation develops in melt pockets away from bubbles and new bubbles may nucleate. Further acceleration and increase in viscosity cause buildup of internal pressure in the bubbles and may eventually lead to fragmentation of the melt. Received: 19 June 1995 / Accepted: 27 December 1995     

Physical properties of the 1980 Mount St. Helens cryptodome magma

 Physical properties of cryptodome and remelted samples of the Mount St. Helens grey dacite have been measured in the laboratory. The viscosity of cryptodome dacite measured by parallel–plate viscometry ranges from 10.82 to 9.94 log₁₀ η (Pa s) (T=900–982  °C), and shrinkage effects were dilatometrically observed at T>900  °C. The viscosity of remelted dacite samples measured by the micropenetration method is 10.60–9.25 log₁₀ η (Pa s) (T=736–802  °C) and viscosities measured by rotational viscometry are 3.22–1.66 log₁₀ η (Pa s) (T=1298–1594  °C). Comparison of the measured viscosity of cryptodome dacitic samples with the calculated viscosity of corresponding water-bearing melt demonstrates significant deviations between measured and calculated values. This difference reflects a combination of the effect of crystals and vesicles on the viscosity of dacite as well as the insufficient experimental basis for the calculation of crystal-bearing vesicular melt viscosities at low temperature. Assuming that the cryptodome magma of the 18 May 1980 Mount St. Helens eruption was residing at 900  °C with a phenocryst content of 30 vol.%, a vesicularity of 36 vol.% and a bulk water content of 0.6 wt.%, we estimate the magma viscosity to be 10^10.8 Pa s. Received: 25 August 1996 / Accepted: 19 July 1997     

Geochemistry of the thermal springs and fumaroles of Basse-Terre Island,Guadeloupe, Lesser Antilles

 The purpose of this work was to study jointly the volcanic-hydrothermal system of the high-risk volcano La Soufrière, in the southern part of Basse-Terre, and the geothermal area of Bouillante, on its western coast, to derive an all-embracing and coherent conceptual geochemical model that provides the necessary basis for adequate volcanic surveillance and further geothermal exploration. The active andesitic dome of La Soufrière has erupted eight times since 1660, most recently in 1976–1977. All these historic eruptions have been phreatic. High-salinity, Na–Cl geothermal liquids circulate in the Bouillante geothermal reservoir, at temperatures close to 250  °C. These Na–Cl solutions rise toward the surface, undergo boiling and mixing with groundwater and/or seawater, and feed most Na–Cl thermal springs in the central Bouillante area. The Na–Cl thermal springs are surrounded by Na–HCO₃ thermal springs and by the Na–Cl thermal spring of Anse à la Barque (a groundwater slightly mixed with seawater), which are all heated through conductive transfer. The two main fumarolic fields of La Soufrière area discharge vapors formed through boiling of hydrothermal aqueous solutions at temperatures of 190–215  °C below the "Ty" fault area and close to 260  °C below the dome summit. The boiling liquid producing the vapors of the Ty fault area has δD and δ¹⁸O values relatively similar to those of the Na–Cl liquids of the Bouillante geothermal reservoir, whereas the liquid originating the vapors of the summit fumaroles is strongly enriched in ¹⁸O, due to input of magmatic fluids from below. This process is also responsible for the paucity of CH₄ in the fumaroles. The thermal features around La Soufrière dome include: (a) Ca–SO₄ springs, produced through absorption of hydrothermal vapors in shallow groundwaters; (b) conductively heated, Ca–Na–HCO₃ springs; and (c) two Ca–Na–Cl springs produced through mixing of shallow Ca–SO₄ waters and deep Na–Cl hydrothermal liquids. The geographical distribution of the different thermal features of La Soufrière area indicates the presence of: (a) a central zone dominated by the ascent of steam, which either discharges at the surface in the fumarolic fields or is absorbed in shallow groundwaters; and (b) an outer zone, where the shallow groundwaters are heated through conduction or addition of Na–Cl liquids coming from hydrothermal aquifer(s). Received: 9 November 1998 / Accepted: 15 July 1999     

Effusive eruption of viscous silicic magma triggered and driven by recharge: a case study of the Cerro Chascon-Runtu Jarita Dome Complex in Southwest Bolivia

 The Cerro Chascon-Runtu Jarita Complex is a group of ten Late Pleistocene (∼85 ka) lava domes located in the Andean Central Volcanic Zone of Bolivia. These domes display considerable macroscopic and microscopic evidence of magma mixing. Two groups of domes are defined chemically and geographically. A northern group, the Chascon, consists of four lava bodies of dominantly rhyodacite composition. These bodies contain 43–48% phenocrysts of plagioclase, quartz, sanidine, biotite, and amphibole in a microlite-poor, rhyolitic glass. Rare mafic enclaves and selvages are present. Mineral equilibria yield temperatures from 640 to 750  °C and log ƒO₂ of –16. Geochemical data indicate that the pre-eruption magma chamber was zoned from a dominant volume of 68% to minor amounts of 76% SiO₂. This zonation is best explained by fractional crystallization and some mixing between rhyodacite and more evolved compositions. The mafic enclaves represent magma that intruded but did not chemically interact much with the evolved magmas. A southern group, the Runtu Jarita, is a linear chain of six small domes (<1 km³ total volume) that probably is the surface expression of a dike. The five most northerly domes are composites of dacitic and rhyolitic compositions. The southernmost dome is dominantly rhyolite with rare mafic enclaves. The composite domes have lower flanks of porphyritic dacite with ∼35 vol.% phenocrysts of plagioclase, orthopyroxene, and hornblende in a microlite-rich, rhyodacitic glass. Sieve-textured plagioclase, mixed populations of disequilibrium plagioclase compositions, xenocrystic quartz, and sanidine with ternary composition reaction rims indicate that the dacite is a hybrid. The central cores of the composite domes are rhyolitic and contain up to 48 vol.% phenocrysts of plagioclase, quartz, sanidine, biotite, and amphibole. This is separated from the dacitic flanks by a banded zone of mingled lava. Macroscopic, microscopic, and petrologic evidence suggest scavenging of phenocrysts from the silicic lava. Mineral equilibria yield temperatures of 625–727  °C and log ƒO₂ of –16 for the rhyolite and 926–1000  °C and log ƒO₂ of –9.5 for the dacite. The rhyolite is zoned from 73 to 76% SiO₂, and fractionation within the rhyolite composition produced this variation. Most of the 63–73% SiO₂ compositional range of the lava in this group is the result of mixing between the hybrid dacite and the rhyolite. Eruption of both groups of lavas apparently was triggered by mafic recharge. A paucity of explosive activity suggests that volatile and thermal exchanges between reservoir and recharge magmas were less important than volume increase and the lubricating effects of recharge by mafic magmas. For the Runtu Jarita group, the eruption is best explained by intrusion of a dike of dacite into a chamber of crystal-rich rhyolite close to its solidus. The rhyolite was encapsulated and transported to the surface by the less-viscous dacite magma, which also acted as a lubricant. Simultaneous effusion of the lavas produced the composite domes, and their zonation reflects the subsurface zonation. The role of recharge by hotter, more fluid mafic magma appears to be critical to the eruption of some highly viscous silicic magmas. Received: 23 August 1998 / Accepted: 10 March 1999     

Viscosity of microlite-bearing rhyolitic obsidians: an experimental study

 To investigate the influence of microlites on lava flow rheology, the viscosity of natural microlite-bearing rhyolitic obsidians of calc-alkaline and peralkaline compositions containing 0.1–0.4 wt.% water was measured at volcanologically relevant temperatures (650–950  °C), stresses (10³–10⁵ Pa) and strain rates (10^–5 to 10^–7 s^–1). The glass transition temperatures (T _g ) were determined from scanning calorimetric measurements on the melts for a range of cooling/heating rates. Based on the equivalence of enthalpic (calorimetric) and shear (viscosity) relaxation, we calculated the viscosity of the melt in crystal-bearing samples from the T _g data. The difference between the calculated viscosity of the melt phase and the measured viscosity for the crystal-bearing samples is interpreted to be the physical effect of microlites on the measured viscosity. The effect of <5 vol.% rod-like microlites on the melt rheology is negligible. Microlite-rich and microlite-poor samples from the same lava flow and with identical bulk chemistry show a difference of 0.6 log₁₀ units viscosity (Pa s), interpreted to be due to differences in melt chemistry caused by the presence of microlites. The only major differences between measured and calculated viscosities were for two samples: a calc-alkaline rhyolite with 1 vol.% branching crystals, and a peralkaline rhyolite containing crystal-rich bands with >45 vol.% crystals. For both of these samples a connectivity factor is apparent, with, for the latter, a close packing framework of crystals which is interpreted to influence the apparent viscosity. Received: 14 March 1996 / Accepted: 30 May 1996     

High-temperature emplacement of the Cerro Galán and Toconquis Group ignimbrites (Puna plateau,NW Argentina) determined by TRM analyses

Estimates of pyroclastic flow emplacement temperatures in the Cerro Galán ignimbrite and Toconquis Group ignimbrites were determined using thermal remanent magnetization of lithic clasts embedded within the deposits. These ignimbrites belong to the Cerro Galán volcanic system, one of the largest calderas in the world, in the Puna plateau, NW Argentina. Temperature estimates for the 2.08-Ma Cerro Galán ignimbrite are retrieved from 40 sites in 14 localities (176 measured clasts), distributed at different distances from the caldera and different stratigraphic heights. Additionally, temperature estimates were obtained from 27 sample sites (125 measured clasts) from seven ignimbrite units forming the older Toconquis Group (5.60–4.51 Ma), mainly outcropping along a type section at Rio Las Pitas, Vega Real Grande. The paleomagnetic data obtained by progressive thermal demagnetization show that the clasts of the Cerro Galán ignimbrite have one single magnetic component, oriented close to the expected geomagnetic field at the time of emplacement. Results show therefore that most of the clasts acquired a new magnetization oriented parallel to the magnetic field at the moment of the ignimbrite deposition, suggesting that the clasts were heated up to or above the highest blocking temperature (T _b) of the magnetic minerals (T _b = 580°C for magnetite; T _b = 600–630°C for hematite). We obtained similar emplacement temperature estimations for six out of the seven volcanic units belonging to the Toconquis Group, with the exception of one unit (Lower Merihuaca), where we found two distinct magnetic components. The estimation of emplacement temperatures in this latter case is constrained at 580–610°C, which are lower than the other ignimbrites. These estimations are also in agreement with the lowest pre-eruptive magma temperatures calculated for the same unit (i.e., 790°C; hornblende–plagioclase thermometer; Folkes et al. 2011b). We conclude that the Cerro Galán ignimbrite and Toconquis Group ignimbrites were emplaced at temperatures equal to or higher than 620°C, except for Lower Merihuaca unit emplaced at lower temperatures. The homogeneity of high temperatures from proximal to distal facies in the Cerro Galán ignimbrite provides constraints for the emplacement model, marked by a relatively low eruption column, low levels of turbulence, air entrainment, surface–water interaction, and a high level of topographic confinement, all ensuring minimal heat loss.     

Vesiculation path of ascending magma in the 1983 and the 2000 eruptions of Miyakejima volcano, Japan

The vesiculation of magma during the 1983 eruption of Miyakejima Volcano, Japan, is discussed based on systematic investigations of water content, vesicularity, and bubble size distribution for the products. The eruption is characterized by simultaneous lava effusion and explosive sub-plinian (‘dry’) eruptions with phreatomagmatic (‘wet’) explosions. The magmas are homogeneous in composition (basaltic andesite) and in initial water content (H₂O = 3.9±0.9 wt%), and residual groundmass water contents for all eruption styles are low (H₂O <0.4 wt%) suggestive of extensive dehydration of magma. For the scoria erupted during simultaneous ‘dry’ and ‘wet’ explosive eruptions, inverse correlation was observed between vesicularity and residual water content. This relation can be explained by equilibrium exsolution and expansion of ca. 0.3 wt% H₂O at shallow level with different times of quenching, and suggests that each scoria with different vesicularity, which was quenched at a different time, provides a snapshot of the vesiculation process near the point of fragmentation. The bubble size distribution (BSD) varies systematically with vesicularity, and total bubble number density reaches a maximum value at vesicularity Φ ∼ 0.5. At Φ  ∼ 0.5, a large number of bubbles are connected with each other, and the average thickness of bubble walls reaches the minimum value below which they would rupture. These facts suggest that vesiculation advanced by nucleation and growth of bubbles when Φ < 0.5, and then by expansion of large bubbles with coalescence of small ones for Φ > 0.5, when bubble connection becomes effective. Low vesicularity and low residual water content of lava and spatter (Φ  < 0.1, H₂O  < 0.1 wt%), and systematic decrease in bubble number density from scoria through spatter to lava with decrease in vesicularity suggest that effusive eruption is a consequence of complete degassing by bubble coalescence and separation from magma at shallow levels when magma ascent rate is slow.

Fumarolic activity of Avachinsky and Koryaksky volcanoes, Kamchatka, from 1993 to 1994

 Volcanic gas and condensate samples were collected in 1993–1994 from fumaroles of Koryaksky and Avachinsky, basaltic andesite volcanoes on the Kamchatka Peninsula near Petropavlovsk–Kamchatsky. The highest-temperature fumarolic discharges, 220  °C at Koryaksky and 473  °C at Avachinsky, are water-rich (940–985 mmol/mol of H₂O) and have chemical and isotopic characteristics typical of Kamchatka–Kurile, high- and medium-temperature volcanic gases. The temperature and chemical and water isotopic compositions of the Koryaksky gases have not changed during the past 11 years. They represent an approximate 2 : 1 mixture of magmatic and meteoric end members. Low-temperature, near-boiling-point discharges of Avachinsky Volcano are water poor (≈880 mmol/mol); Their compositions have not changed since the 1991 eruption, and are suggested to be derived from partially condensed magmatic gases at shallow depth. Based on a simple model involving mixing and single-step steam separation, low water and high CO₂ contents, as well as the observed Cl concentration and water isotopic composition in low-temperature discharges, are the result of near-surface boiling of a brine composed of the almost pure condensed magmatic gas. High methane content in low-temperature Avachinsky gases and the 220  °C Koryaksky fumarole, low C isotopic ratio in CO₂ at Koryaksky (–11.8‰), and water isotope data suggest that the "meteoric" end member contains considerable amounts of the regional methane-rich thermal water discovered in the vicinity of both volcanoes. Received: 2 May 1996 / Accepted: 5 November 1996     

The influence of trace amounts of water on the viscosity of rhyolites

 As a major volatile in volcanic systems, water has a significant influence on the rheological properties of silicic magmas. This is especially so at minor water contents relevant to the emplacement of silicic lavas. To investigate the influence of water on the viscosity of natural rhyolitic obsidians, a novel strategy has been adopted employing parallel-plate and micropenetration techniques. Viscosities have been determined on three types of material: (a) raw water-bearing obsidians; (b) remelted (1650  °C, 1 atm) degassed glasses of the obsidians; and (c) hydrothermally hydrated (1300  °C, 3 kbar) obsidians. Ten natural rhyolitic obsidians (peraluminous, calc-alkaline and peralkaline) were employed: seven originated from lava flows and contained <0.2 wt.% H₂O, two samples were F-rich from pyroclastic successions, and one was an obsidian cobble with 1.5 wt.% water also associated with pyroclastic units. Melt compositions and water contents were stable during viscometry. The measured decreases in activation energies of viscous flow and viscosity with small amounts of water are much greater than the Shaw calculation scheme predicts. In addition, a marked non-linear decrease in η exists with increasing water content. In contrast to the case for peralkaline rhyolites, 0.1–0.2 wt.% water decreases activation energies significantly (up to 30%) for calc-alkaline compositions. These results have important implications for the ease of near-surface degassing of silicic magmas during emplacement and permit the testing of calculational models for viscosity, largely based on synthetic systems. Received: 7 July 1997 / Accepted: 6 April 1998     

Experimental simulations of gas-driven eruptions: kinetics of bubble growth and effect of geometry

 Simulated gas-driven eruptions using CO₂–water-polymer systems are reported. Eruptions are initiated by rapidly decompressing CO₂–saturated water containing up to 1.0 wt.% CO₂. Both cylindrical test cells and a flask test cell were used to examine the effect of magma chamber/conduit geometry on eruption dynamics. Bubble-growth kinetics are examined quantitatively in experiments using cylindrical test cells. Uninhibited bubble growth can be roughly expressed as dr/dt≈λD(β-1)/(γt ^1/3) for a CO₂–water-polymer system at 0–22  °C and with viscosities up to 5 Pa·s, where r is the radius of bubbles, λ and D are the Ostwald solubility coefficient and diffusivity of the gas in the liquid, β is the degree of saturation (decompression ratio), and γ characterizes how the boundary layer thickness increases with time and is roughly 1.0×10^–5 m/s^1/3 in this system. Unlike the radius of cylindrical test cells, which does not affect the eruption threshold and dynamics, the shape of the test cells (flask vs cylindrical) affects the dynamics but not the threshold of eruptions. For cylindrical test cells, the front motion is characterized by constant acceleration with both Δh (the height increase) and ΔV (the volume increase) being proportional to t ²; for the flask test cell, however, neither Δh nor ΔV is proportional to t ² as the conduit radius varies. Test-cell geometry also affects foam stability. In the flask test cell, as it moves from the wider base chamber into the narrower conduit, the bubbly flow becomes fragmented, affecting the eruption dynamics. The fragmentation may be caused by a sudden increase in acceleration induced by conduit-shape change, or by the presence of obstacles to the bubbly flow. This result may help explain the range in vesicularities of pumice and reticulite. Received: 16 May 1997 / Accepted: 11 October 1997     

Cooling and crystallization of lava in open channels, and the transition of Pāhoehoe Lava to 'A'ā

 Samples collected from a lava channel active at Kīlauea Volcano during May 1997 are used to constrain rates of lava cooling and crystallization during early stages of flow. Lava erupted at near-liquidus temperatures (∼1150  °C) cooled and crystallized rapidly in upper parts of the channel. Glass geothermometry indicates cooling by 12–14  °C over the first 2 km of transport. At flow velocities of 1–2 m/s, this translates to cooling rates of 22–50  °C/h. Cooling rates this high can be explained by radiative cooling of a well-stirred flow, consistent with observations of non-steady flow in proximal regions of the channel. Crystallization of plagioclase and pyroxene microlites occurred in response to cooling, with crystallization rates of 20–50% per hour. Crystallization proceeded primarily by nucleation of new crystals, and nucleation rates of ∼10⁴/cm³s are similar to those measured in the 1984 open channel flow from Mauna Loa Volcano. There is no evidence for the large nucleation delays commonly assumed for plagioclase crystallization in basaltic melts, possibly a reflection of enhanced nucleation due to stirring of the flow. The transition of the flow surface morphology from pāhoehoe to 'a'ā occurred at a distance of 1.9 km from the vent. At this point, the flow was thermally stratified, with an interior temperature of ∼1137  °C and crystallinity of ∼15%, and a flow surface temperature of ∼1100  °C and crystallinity of ∼45%. 'A'ā formation initiated along channel margins, where crust was continuously disrupted, and involved tearing and clotting of the flow surface. Both observations suggest that the transition involved crossing of a rheological threshold. We suggest this threshold to be the development of a lava yield strength sufficient to prevent viscous flow of lava at the channel margin. We use this concept to propose that 'a'ā formation in open channels requires both sufficiently high strain rates for continued disruption of surface crusts and sufficient groundmass crystallinity to generate a yield strength equivalent to the imposed stress. In Hawai'i, where lava is typically microlite poor on eruption, these combined requirements help to explain two common observations on 'a'ā formation: (a) 'a'ā flow fields are generated when effusion rates are high (thus promoting crustal disruption); and (b) under most eruption conditions, lava issues from the vent as pāhoehoe and changes to 'a'ā only after flowing some distance, thus permitting sufficient crystallization. Received: 3 September 1998 / Accepted: 12 April 1999     

Sub-volcanic crystallization at Stromboli (Aeolian Islands, southern Italy) preceding the Sciara del Fuoco sector collapse: evidence from monzonite lithic suite

 We describe texture, mineralogy and whole-rock composition of cognate monzonite sub-volcanic clasts within debris flow deposits related to the 5000 years catastrophic phreatomagmatic eruption probably linked to the Sciara del Fuoco sector collapse. The debris flows are at the top of accretionary lapilli-rich ash deposits overlying potassic (KS, shoshonites) lavas of the Neostromboli period. The monzonites are inferred to be crystallized in situ, at low P, at the side walls and/or roof margins of a shallow magma chamber and to be cogenetic with the KS Neostromboli extrusives. They can be considered "ideal orthocumulates" since they approximately retain a bulk liquid composition and possibly represent "slowly cooled equivalents" of their KS shoshonite host rock. The "closure temperature" of final solidification of the monzonite lithic suite was estimated through ternary-feldspar geothermobarometry, plagioclase–K-feldspar and K-feldspar–biotite equilibria and is in the range of 750–790  °C with a maximum –logfO₂ around 15.1–15.3. The estimated pressure of crystallization is <0.5 kbar. Potassic lavas and dikes, previously emplaced during the Neostromboli period, also resemble the monzonites in both major trace elements and mineral chemistry. The cogenetic relationship between KS Neostromboli extrusives and the monzonite host-rock magma from which the sub-volcanic clasts were derived is clear evidence that a shallow magma chamber existed between the caldera collapse of the Vancori period and the Sciara del Fuoco sector collapse (i.e. between 13 000 and 5000 years). The monzonite clasts were derived from crystallization at very shallow depth (ca. 1 km) and strongly support the hypothesis of violent decompression of the shallow magmatic plumbing system during the Sciara del Fuoco sector collapse. Climax of the regressive landslide event, with maximum disruption of the chamber walls, took place during emplacement of the debris flows, i.e. during the late stage of the Neostromboli phreatomagmatic eruption. Received: 15 September 1996 / Accepted: 5 May 1997     

Geochemistry of the volcano-hydrothermal system of El Chichón Volcano, Chiapas, Mexico

 The 1982 eruption of El Chichón volcano ejected more than 1 km³ of anhydrite-bearing trachyandesite pyroclastic material to form a new 1-km-wide and 300-m-deep crater and uncovered the upper 500 m of an active volcano-hydrothermal system. Instead of the weak boiling-point temperature fumaroles of the former lava dome, a vigorously boiling crater spring now discharges  / 20 kg/s of Cl-rich (∼15 000 mg/kg) and sulphur-poor ( / 200 mg/kg of SO₄), almost neutral (pH up to 6.7) water with an isotopic composition close to that of subduction-type magmatic water (δD=–15‰, δ¹⁸O=+6.5‰). This spring, as well as numerous Cl-free boiling springs discharging a mixture of meteoric water with fumarolic condensates, feed the crater lake, which, compared with values in 1983, is now much more diluted (∼3000 mg/kg of Cl vs 24 030 mg/kg), less acidic (pH=2.6 vs 0.56) and contains much lower amounts of S ( / 200 mg/kg of SO₄, vs 3550 mg/kg) with δ³⁴S=0.5–4.2‰ (+17‰ in 1983). Agua Caliente thermal waters, on the southeast slope of the volcano, have an outflow rate of approximately 100 kg/s of 71  °C Na–Ca–Cl water and are five times more concentrated than before the eruption (B. R. Molina, unpublished data). Relative N₂, Ar and He gas concentrations suggest extensional tectonics for the El Chichón volcanic centre. The ³He/⁴He and ⁴He/²⁰Ne ratios in gases from the crater fumaroles (7.3R_a, 2560) and Agua Caliente hot springs (5.3R_a, 44) indicate a strong magmatic contribution. However, relative concentrations of reactive species are typical of equilibrium in a two-phase boiling aquifer. Sulphur and C isotopic data indicate highly reducing conditions within the system, probably associated with the presence of buried vegetation resulting from the 1982 eruption. All Cl-rich waters at El Chichón have a common source. This water has the appearence of a "partially matured" magmatic fluid: condensed magmatic vapour neutralized by interaction with fresh volcaniclastic deposits and depleted in S due to anhydrite precipitation. Shallow ground waters emerging around the volcano from the thick cover of fresh pumice deposits (Red waters) are Ca–SO₄–rich and have a negative oxygen isotopic shift, probably due to ongoing formation of clay at low temperatures. Received: 21 July 1997 / Accepted: 4 December 1997     

An evidence for earthquake occurrence time related with geologic structure — Relations between local mean lunar times of earthquake occurrence and their fault trends

The relation between the local mean lunar time τ of earthquake occurrence and their fault trends is studied in this paper. The local mean lunar times τ of 53 earthquakes in 24 groups are calculated. Because the tidal generation force arisen by the moon is a cyclic function of about 12 hours 25 minutes in the main, the two tidal generation forces anywhere in the earth arising by the moon are equal in general when the moon lies to the two sites of 180° interval of local mean lunar time. Based on this phenomenon the values Δτ of τ₁–τ₂ or τ₁–τ₂ ± 180° of two earthquakes occurring repetitiously in the same place are also calculated. The calculated results show that if the fault trends of the two earthquakes in the same place is near, the Δτ is usually smaller and if the fault trends of the two ones is not near, the Δτ is usually larger and the distribution of the local mean lunar time τ of earthquakes in different places is dispersive even if fault trends of these earthquakes are near, and the τ does not concentrate on the lower and upper transit of the moon. The above phenomena clear up that the triggering earthquake of earth solid tide arisen by the moon is relative with the fault trends of earthquakes and we ought to think over the difference of environmental conditions of earthquake preparation of each seismogenic zone and can not make statistics to earthquakes in different places when we study the relation between the solid earth tide arisen by the moon and earthquakes.     

Geochemistry of gases and waters discharged by the mud volcanoes at Paternò, Mt. Etna (Italy)

 Approximately 20 km south of Mt. Etna craters, at the contact between volcanic and sedimentary formations, three mud volcanoes discharge CO₂-rich gases and Na–Cl brines. The compositions of gas and liquid phases indicate that they are fed by a hydrothermal system for which temperatures of 100–150  °C were estimated by means of both gas and solute geothermometry. The hydrothermal system may be associated with CO₂-rich groundwaters over a large area extending from the central part of Etna to the mud volcanoes. Numerous data on the He, CH₄, CO₂ composition of the gases of the three manifestations, sampled over the past 5 years, indicate clearly that variations are due to separation processes of a CO₂-rich gas phase from the liquid. The effects of these processes have to be taken into account in the interpretation of the monitoring data collected for the geochemical surveillance of Etna volcano. Received: 4 September 1995 / Accepted: 14 February 1996     

Magnetic properties and emplacement of the Bishop tuff,California

 Anisotropy of magnetic susceptibility (AMS) and characteristic remanence were measured for 45 sites in the 0.76 Ma Bishop tuff, eastern California. Thirty-three sites were sampled in three stratigraphic sections, two in Owens gorge south of Long Valley caldera, and the third in the Adobe lobe north of Long Valley. The remaining 12 sites are widely distributed, but of limited stratigraphic extent. Weakly indurated, highly porous to dense, welded ash-flow tuffs were sampled. Saturation magnetization vs temperature experiments indicate two principal iron oxide phases: low Ti magnetites with 525–570  °C Curie temperatures, and maghemite with 610°–640  °C Curie temperatures. AF demagnetization spectra of isothermal remanent magnetizations are indicative of magnetite/maghemite predominantly in the multidomain to pseudo-single domain size ranges. Remeasurement of AMS after application of saturating direct fields indicates that randomly oriented single-domain grains are also present. The degree of anisotropy is only a few percent, typical of tuffs. The AMS ellipsoids are oblate with K_min axes normal to subhorizontal foliation and K_max axes regionally aligned with published source vents. For 12 of 16 locality means, K_max axes plunge sourceward, confirming previous observations regarding flow sense. Topographic control on flow emplacement is indicated by the distribution of tuff deposits and by flow directions inferred from K_max axes. Deposition east of the Benton range occurred by flow around the south end of the range and through two gaps (Benton notch and Chidago gap). Flow down Mammoth pass of the Sierra Nevada is also evident. At least some of the Adobe lobe in the northeast flowed around the west end of Glass mountain. Eastward flow directions in the upper Owens gorge and southeast directions in the lower Owens gorge are parallel to the present canyon, suggesting that the present drainage has been established along the pre-Bishop paleodrainage. Characteristic remanence directions from 45 sites (267 samples) yield an overall mean of D=348°, I=53° for the Bishop tuff. A correlation is found in two of the three profiles between density and remanence inclination. A mean remanence direction based on 13 localities together with data from uncompacted xenoliths and data from the ash-fall tuff at Lake Tecopa is: D=353°, I=54°, k=172, α₉₅=2.9°, N=15. Received: 11 July 1995 / Accepted: 29 February 1996     

Gas retention in fine-grained pyroclastic flow materials at high temperatures

The ability of a dense pyroclastic flow to maintain high gas pore pressure, and hence low friction, during runout is determined by (1) the strengths and longevities of gas sources, and (2) the ability of the material to retain residual gas once those sources become ineffective. The latter is termed the gas retention capacity. Gas retention capacity in a defluidizing granular material is governed by three timescales: one for the evacuation of bubbles (t _be ; brief and not considered in this paper), one for hindered settling from the expanded state (t _sett), and one for diffusive release of residual pore pressure from the non-expanded state (t _diff). The relative magnitides of t _sett and t _diff depend on bed thickness, t _sett dominating in thin systems and t _diff in thick ones. Three pyroclastic flow materials, two ignimbrites and a block-and-ash flow sample, were studied experimentally to investigate expansion behaviour under gas flow and to determine gas retention times. Effects of particle size were evaluated by using two size cuts (<4 mm and <250 μm) from each sample. Careful drying of the materials was necessary to avoid effects of humidity-related cohesion. Two sets of experiments were carried out: (1) expansion in the non-bubbling regime at 50–200°C, (2) bed collapse tests from the initially bubbling state at 50–550°C. Provided that gas channelling was avoided by gentle stirring, all the samples exhibited a regime of uniform expansion prior to the onset of bubbling. Fine particle size (in particular high fines content), low particle density and high temperature all favoured smoother fluidization by increasing the maximum expansion possible in the non-bubbling state. An empirical equation describing the uniform expansion of the materials was determined. High temperature also favoured greater gas partitioning into the dense phase of the bubbling bed, as well (in finer-grained samples) as higher voidage in the settled bed. Large values of t _sett and t _diff were favoured by fine particle size. Temperature had less influence, suggesting that experimental results at low temperatures (50–200°C) can be extrapolated to higher temperatures. Gas retention times provide insight into the ability of pyroclastic flows in expanded (t _sett) or non-expanded (t _diff) flow states to retain gas once air ingestion or gas production have become ineffective. Finer-grained pyroclastic flows are expected to retain gas longer, and hence to have higher apparent ‘mobilities’, than coarser-grained ones of comparable volume, as has been observed on Montserrat.     

Dependence of the energy released during earthquake on ambient shear stress

Starting with dislocation model, using the result of the fracture mechanics: the slip displacement at the crack tip is proportional to the length of the crack and the applied ambient shear stressτ ₀ ² , we consider the dislocation in the earthquake to be the slip displacement at the crack tip and have obtained the analysis expresses of displacement and velocity pulse for the circular crack and have calculated the seismic wave energy radiated by earthquake. The seismic wave energyE ∞M ₀ τ ₀ ² f(v) _r , i. e.E is proportional to the seismic momentM ₀ and the square of the ambient shear stressτ ₀ ² and increases with the rupture velocityv _r . In frequency domain, integrating the square of source velocity spectrum derived from our the scaling law model, we have also obtained the seismic wave energyE released by earthquake and earthquake radiated effficiencyη.E ∞M ₀ τ ₀ ² also. If takingτ ₀ = 10.0 MPa, E=4.79M ₀. This result is consistent with the estimate by Vassiliou and Kanamori (1982). Theη=5.26%. The distribution of the seismic wave energy is that most of the energy contains in the frequency range between the first corner frequencyf _c1 and thirdf _c3, amount to 92.3% the energy in the rangef<f _c1 is about 3.85% and 3.85% whenf>f _c3. Thef _c3 is about 8Hz forM ⩾ 6, thus most of radiated energy is below 2Hz. This phenomenon had been verified by Vassiliou Kanamori. Previous results show the energy radiated by earthquake to be strongly dependent on ambient shear stress. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 146–152, 1993. This work was supported by the Deutsche Forschungsgemeinschaft, Bonn, F. R. Germany. The support is grateful acknowledged. The authors are also grateful to Professor Klussmann and Mr. Hasthoff for their lots of help.     

The space and time distribution characteristics of the shear stress field for the sequence of the Wuding earthquake

Follow Chen and Duda's model of spectral fall-off of (3, the dependence of peak parameters of ground motion, peak displacement dm, peak velocity vm and peak acceleration am, upon the environment stress (0-values are studied using near source seismic digital recordings for the sequence of the Wuding, Yunnan, M = 6.5 earthquake, in which, as a new thought, the peak parameters are assumed to be related to the medium Q-value. Three formulae for estimating the environment stress (0-values by the peak parameters of three types of ground motions are derived. Using these formulae, the environment stress (0-values are calculated for the sequence of the Wuding earthquake. The result show that (0-values calculated by the three formulae are constant largely, the averages of (0 are in the range of 5.0～35 MPa for most earthquakes. It belongs to the high-stress earthquakes sequence: the high-stress values are restricted to the relatively small area closely near to the epicenter of the main shock. The fine distribution structure for the contours of the environment stress (0-values is related closely to the strong aftershocks. The analysis of spatial and temporal feature of (0-values suggests that the earthquakes sequence in a rupture process generated at the specific intersection zone of seismo-tectonics under high-stress background.     

Fluid geochemistry of Montserrat Island,West Indies

 Two geochemical surveys carried out in March 1991 and September 1992 revealed the existence of a hydrothermal system in the southern portion of Montserrat Island, below Soufrière Hills Volcano. This conclusion is supported by the presence of: (a) the thermal springs of Plymouth which are fed by deep Na–Cl waters (Cl concentration ∼25 000 mg/kg, temperature ca. 250  °C) mixed with shallow steam-heated waters; (b) the four fumarolic fields of Galway's Soufrière, Gages Upper Soufrière, Gages Lower Soufrière, and Tar River Soufrière, where acid to neutral, steam-heated waters are present together with several fumarolic vents, discharging vapors formed through boiling of hydrothermal aqueous solutions. Involvement of magmatic fluids in the recharge of the hydrothermal aquifers is suggested by: (a) the high ³He/⁴He ratios of fumarolic fluids, i.e., 8.2 R_A at Galway's Soufrière and 5.9 R_A at Gages Lower Soufrière; (b) the δD and δ¹⁸O values of Na–Cl thermal springs and steam condensates, indicating the involvement of arc-type magmatic water in the formation of deep geothermal liquids; and (c) the CH₄/CO₂ ratios of fumarolic fluids, which are lower than expected for equilibrium with the FeO–FeO_1.5 hydrothermal rock buffer, but being shifted towards the SO₂–H₂S magmatic gas buffer. Received: 26 March 1996 / Accepted: 19 July 1996