Continuous chemical monitoring of volcanic gas in Izu-Oshima volcano, Japan

A new continuous monitoring system has been developed for the measurement of volcanic gas from the steam well located 3 km north from the summit of Izu-Oshima volcano, Japan. After removing the water vapor using three sequential dehydration methods, CO₂ and SO₂ contents are measured using IR sensors, and O₂ and H₂ using a zirconia sensor and a semiconductor sensor, respectively. This system has been in operation without any significant trouble for 3 years.The dehydrated volcanic gas from the well consists of a mixture of CO₂, O₂ and N₂. A decreasing trend of the CO₂ content was observed from 1995 to 1998 together with a decrease of volcanic activity. Seasonal changes have also been observed in CO₂ and O₂ contents, CO₂ being higher and O₂ lower in summer, which suggests larger contribution of magmatic components in summer. While changes in short-term variation in CO₂ and O₂ are influenced by atmospheric pressure changes; the CO₂ content correlates inversely with atmospheric pressure unlike O₂ with some hours delay. In contrast, the H₂ content increased intermittently up to 1200 ppm one to several hours after a sudden drop in the atmospheric pressure and without any apparent correlation with seasonal changes.This system allows us to study temporal variation in chemical composition of volcanic gas during quiescent periods of volcanic activity of Izu-Oshima volcano, and might help us detect anomalous changes before future eruptive events.     

Thermal, seismic and infrasonic evidences of variable degassing rates at Stromboli volcano

During June 1999, we measured the amplitude and rate (number of events per second) at which gas exited the vent at Stromboli volcano as discrete gas bursts or puffs. This allowed us to identify two styles of gas burst (puffing) activity. The first is characterized by frequent, rapidly rising puffs, the second by less frequent, slowly rising puffs. Each style persisted over 5–40-min-long durations and was associated with a high and low number of strombolian explosions per hour, respectively. Each period was also associated with characteristic delay times between the arrival of the infrasonic and thermal signals during strombolian explosions; the delays were longer during vigorous puffing periods. To explain our observations, we propose a model in which the degassing process cycles between vigorous and weak degassing phases. During vigorous degassing phases, bubble layers ascend the conduit at a frequency of 0.5–1.0 s⁻¹. This high degassing level reflects a gas-rich magma column and leads to an increased rate in the formation of shallow foams and, hence, an increase in puffing and explosive activity, as well as a higher free surface level and/or gas jet velocity. During weak phases, bubble layers ascend the conduit at a reduced frequency of 0.2–0.3 s⁻¹. During such times the magma column is poor in gas. This leads to a decreased rate of foam layer formation and hence a reduction in puffing and explosive activity, as well as a lower free surface level and/or gas jet velocity. Variations in puffing activity can thus be used to track changes in the rate at which the shallow system is supplied by fresh, gas-rich magma. Our observations indicate that the two degassing styles last from 5 to 40 min and that the switch from one to the other occurs over a matter of minutes.     

The magma evolution of Tianchi volcano, Changbaishan

The Changbaishan Tianchi volcano is composed of the basaltic rocks at the shield-forming stage, the trachyte and pantellerite at the cone-forming stage and modern eruption. Studies on their REE, incompatible elements and Sr, Nd, Pb isotopes suggest that rocks at different stages have a common magma genesis and close evolution relationship with differentiation crystallization playing the key role. The co-eruption of basaltic trachyandesite magma and pantellerite magma indicates that there exist both crustal magma chamber and mantle magma reservoir beneath the Tianchi volcano. Project supported by the National Natural Science Foundation of China (Grant No. 49672109).     

Sulfur dioxide emissions from Popocatépetl volcano (Mexico): case study of a high-emission rate, passively degassing erupting volcano

Popocatépetl volcano in central Mexico has been erupting explosively and effusively for almost 4 years. SO₂ emission rates from this volcano have been the largest ever measured using a COSPEC. Pre-eruptive average SO₂ emission rates (2–3 kt/d) were similar to the emission rates measured during the first part of the eruption (up to August 1995) in contrast with the effusive–explosive periods (March 1996–January 1998) during which SO₂ emission rates were higher by a factor of four (9–13 kt/d). Based on a chronology of the eruption and the average SO₂ emission rates per period, the total SO₂ emissions (up to 1 January 1998) are estimated to be about 9 Mt, roughly half as much as the SO₂ emissions from Mount Pinatubo in a shorter period. Popocatépetl volcano is thus considered as a high-emission rate, passively degassing eruptive volcano. SO₂ emission rates and SO₂ emissions are used here to make a mass balance of the erupted magma and related gases. Identified excess SO₂ is explained in terms of continuous degassing of unerupted magma and magma mixing. Fluctuations in SO₂ emission rate may be a result of convection and crystallization in the chamber or the conduits, cleaning and sealing of the plumbing system, and/or SO₂ scrubbing by the hydrothermal system.     

Simple stochastic modelling of the eruption history of a basaltic volcano: Nyamuragira,Zaire

Three simple models of the behaviour of a series of basaltic eruptions have been tested against the eruptive history of Nyamuragira. The data set contains the repose periods and the volumes of lava emitted in 22 eruptions since 1901. Model 1 is fully stochastic and eruptions of any volume with random repose intervals are possible. Models 2 and 3 are constrained by deterministic limits on the maximum capacity of the magma reservoir and on the lowest drainage level of the reservoir respectively. The method of testing these models involves (1) seeking change points in the time series to determine regimes of uniform magma supply rate, and (2) applying linear regression to these regimes, which for models 2 and 3 are the determinsstic limits to those models. Two change points in the time series for Nyamuragira, in 1958 and 1980, were determined using a Kolmogorov-Smirnov technique. The latter change involved an increase in the magma supply rate by a factor of 2.5, from 0.55 to 1.37 m³s^-1. Model 2 provides the best fit to the behavior of Nyamuragira with the ratio of variation explained by the model to total variation. R², being greater than 0.9 for all three regimes. This fit can be interpreted to mean that there is a determinstic limit to the elastic strength of the magma reservoir 4–8 km below the summit of the volcano.     

ACTIVE system for monitoring volcanic activity: A case study of the Izu-Oshima Volcano,Central Japan

A system is proposed for the monitoring of changes in the underground structure of an active volcano over time by applying a transient electromagnetic method. The monitoring system is named ACTIVE, which stands for Array of Controlled Transient-electromagnetics for Imaging Volcano Edifice. The system consists of a transmitter dipole used to generate a controlled transient electromagnetic (EM) field and an array of receivers used to measure the vertical component of the transient magnetic field at various distances, with automatic operation of both units. In order to verify the performance of the proposed system, numerical and field experiments were carried out by application of the system to the Izu-Oshima volcano, where a remarkable change in the apparent DC resistivity over time had been detected in association with the eruption in 1986.     

Variation of volatile concentration in a magma system of Satsuma-Iwojima volcano deduced from melt inclusion analyses

Chemical analyses of 30 melt inclusions from Satsuma-Iwojima volcano, Japan, were carried out to investigate volatile evolution in a magma chamber beneath the volcano from about 6300 yr BP to the present. Large variations in volatile concentrations of melts were observed. (1) Water concentration of rhyolitic melts decreases with time; 3–4.6 wt.% at the time of latest caldera-forming eruption of Takeshima pyroclastic flow deposit (ca. 6300 yr BP), 3 wt.% for small pyroclastic flow (ca. 1300 yr BP) of Iwodake, post-caldera rhyolitic dome, and 0.7–1.4 wt.% for submarine lava eruption (Showa-Iwojima) in 1934. (2) Rhyolitic melts of the Takeshima and Iwodake eruptions contained CO₂ of less than 40 ppm, while the Showa-Iwojima melt has higher CO₂ concentration of up to 140 ppm. (3) Water and CO₂ concentrations of basaltic to andesitic melt of Inamuradake, a post-caldera basaltic scoria cone, are 1.2–2.8 wt.% and ≤290 ppm, respectively.Volatile evolution in the magma chamber is interpreted as follows: (1) the rhyolitic magma at the time of the latest caldera-forming eruption (ca. 6300 yr BP) was gas-saturated due to pressure variation in the magma chamber because the large variation in water concentration of the melt was attributed to exsolution of volatile in the magma prior to the eruption. Iwodake eruption (ca. 1300 yr BP) was caused by a remnant of the caldera-forming rhyolitic magma, suggested from the similarity of major element composition between these magmas. (2) Volatile composition of the Showa-Iwojima rhyolitic melt agrees with that of magmatic gases presently discharging from a summit of Iwodake, indicating the low pressure degassing condition. (3) The degassing of the magma chamber by magma convection in a conduit of Iwodake during non-eruptive but active degassing period for longer than 800 years decreased water concentration of the rhyolitic magma. (4) Geological and petrological observations indicate that a stratified magma chamber, which consists of a lower basaltic layer and an upper rhyolitic layer, might have existed during the post-caldera stage. Addition of CO₂ from the underlying basaltic magma to the upper gas-undersaturated (degassed) rhyolitic magma increased CO₂ concentration of the rhyolitic magma.     

Groundmass pargasite in the 1991–1995 dacite of Unzen volcano: phase stability experiments and volcanological implications

Pargasite commonly occurs in the dacitic groundmass of the 1991–1995 eruption products of Unzen volcano. We described the occurrence and chemical compositions of amphibole in the dacite, and also carried out melting experiments to determine the low-pressure stability limit of amphibole in the dacite. The 1991–1995 ejecta of the Unzen volcano show petrographic evidence of magma mixing, such as reverse compositional zoning of plagioclase and amphibole phenocrysts, and we used a groundmass separate as a starting material for the experiments. Reversed experiments show that the maximum temperature for the crystallization of amphibole is 930°C at 196 MPa, 900°C at 98 MPa, and 820°C at 49 MPa. Compared with the experimental results on the Mount St. Helens dacite, present experiments on the Unzen dacitic groundmass show that amphibole is stable to pressures ca. 50 MPa lower at 850°C. Available Fe–Ti oxide thermometry indicates the crystallization temperature of the groundmass of the Unzen dacite to be 880±30°C, suggesting that the groundmass pargasite crystallized at >70 MPa, corresponding to a depth of more than 3 km in the conduit. The chlorine content of the groundmass pargasite is much lower than that of phenocrystic magnesiohornblende in the 1991–1995 dacite of Unzen volcano, indicating that vesiculation/degassing of magma took place before the crystallization of the groundmass pargasite. The present study shows that the magma was water oversaturated and that the degassing of magma along with magma mixing caused crystallization of the groundmass amphibole at depths of more than 3 km in the conduit.     

Anomalous Emissions of SO2 During the Recent Eruption of Santa Ana Volcano, El Salvador, Central America

Santa Ana volcano in western El Salvador, Central America, had a phreatic eruption at 8:05 am (local time) on October 1, 2005, 101 years after its last eruption. However, during the last one hundred years this volcano has presented periods of quiet degassing with fumarolic activity and an acidic lake within its crater. This paper presents results of frequent measurements of SO₂ degassing using the MiniDOAS (Differential Optical Absorption Spectroscopy) system and a comparison with the volcanic seismicity prior to the eruption, during, and after the eruption. Vehicle measurements of SO₂ flux were taken every hour during the first nine days of the eruption and daily after that. The period of time reported here is from August to December, 2005. Three periods of degassing are distinguished: pre-eruptive, eruptive, and post-eruptive periods. The intense activity at Santa Ana volcano started in July 2005. During the pre-eruptive period up to 4306 and 5154 ton/day of SO₂ flux were recorded on October 24 and September 9, 2005, respectively. These values were of the same order of magnitude as the recorded values just after the October 1^st eruption (2925 ton/day at 10:01 am). Hourly measurements of SO₂ flux taken during the first nine days after the main eruptive event indicate that explosions are preceded by an increase in SO₂ flux and that this parameter reaches a peak after the explosion took place. This behavior suggests that increasing accumulation of exsolved magmatic gases occurs within the magmatic chamber before the explosions, increasing the pressure until the point of explosion. A correlation between SO₂ fluxes and RSAM (Real Time Seismic Amplitude Measurements) is observed during the complete sampling period. Periodic fluctuations in the SO₂ and RSAM values during the entire study period are observed. One possible mechanism explaining these fluctuations it that convective circulation within the magmatic chamber can bring fresh magma periodically to shallow levels, allowing increasing degasification and then decreasing degasification as the batch of magma lowers its gas content, becomes denser, and sinks to give space to a new magma pulse. These results illustrate that the measurements of SO₂ flux can provide important warning signals for incoming explosive activity in active volcanoes.     

A 3D magnetic structure of Izu-Oshima Volcano and their changes after the eruption in 1986 as estimated from repeated airborne magnetic surveys

A 3D magnetic inversion method using a conjugate gradient method (CG method) was developed for constructing 3D magnetization models of a volcanic edifice and applied to aeromagnetic anomalies of Izu-Oshima Volcano surveyed in 1986 and in 1997. The calculated results of the 1986 data show that the volcanic edifice of Izu-Oshima Volcano has a mean magnetization intensity ranging from 10.4 to 12.1 A/m. The derived 3D magnetic structure shows low magnetization zones beneath the west-northwest of the western caldera rim, beneath the west-southwest of Mt. Mihara and beneath Mt. Shiroishi. These features may be related with demagnetizations, reflecting a high thermal state due to magma activities in the 1986 eruption. The comparison between 3D magnetization models in 1986 and in 1997, indicates meaningful changes beneath the C-craters erupted in 1986, suggesting a recovery process of demagnetizations and a considerable decrease of magnetization intensities in the foot of Mt. Futago, indicative of demagnetizations. A derived magnetization model including Izu-Oshima Volcano and its surrounding sea areas clarifies the submerged volcanic edifices around Izu-Oshima Island, and suggests that the old volcanic edifices of Fudeshima, Gyojyanoiwaya, and Okata Volcanoes have been affected by eastward migrations due to massive intrusions of a dike-like structure inferred at the base of Izu-Oshima Volcano.     

Equilibrium and disequilibrium degassing of a phonolitic melt (Vesuvius AD 79 “white pumice”) simulated by decompression experiments

Equilibrium and disequilibrium degassing of a volatile phase from a magma of K-phonolitic composition was investigated to assess its behavior upon ascent. Decompression experiments were conducted in Ar-pressurized externally heated pressure vessels at superliquidus temperature (1050 °C), in the pressure range 10–200 MPa using pure water as fluid phase. All experiments were equilibrated at 200 MPa and then decompressed to lower pressures with rates varying from 0.0028 to 4.8 MPa/s. Isobaric saturation experiments were performed at the same temperature and at 900–950 °C to determine the equilibrium water solubility in the pressure range 30–250 MPa. The glasses obtained from decompression experiments were analyzed for their dissolved water content, vesicularity and bubble size distribution. All decompressed samples presented a first event of bubble nucleation at the capsule–melt interface. Homogeneous bubble nucleation in the melt only occurred in fast-decompressed experiments (4.8 and 1.7 MPa/s), for ΔP ≅ 100 MPa. For these decompression rates high water over-saturations were maintained until a rapid exsolution was triggered at ΔP > 150 MPa. For slower rates (0.0028, 0.024, 0.17 MPa/s) the degassing of the melt took place by diffusive growth of the bubbles nucleating at the capsule–melt interface. This process sensibly reduced water over-saturation in the melt, preventing homogeneous nucleation to occur. For decompression rates of 0.024 and 0.17 MPa/s low water over-saturations were attained in the melt, gradually declining toward equilibrium concentrations at low pressures. A near-equilibrium degassing path was observed for a decompression rate of 0.0028 MPa/s. Experimental data combined with natural pumice textures suggest that both homogeneous and heterogeneous bubble nucleations occurred in the phonolitic magma during the AD 79 Vesuvius plinian event. Homogeneous bubble nucleation probably occurred at a depth of ∼ 3 km, in response to a fast decompression of the magma during the ascent.     

Magma degassing and basaltic eruption styles: a case study of 2000 year BP Xitle volcano in central Mexico

To investigate the relationship between volatile abundances and eruption style, we have analyzed major element and volatile (H₂O, CO₂, S) concentrations in olivine-hosted melt inclusions in tephra from the 2000 yr BP eruption of Xitle volcano in the central Trans-Mexican Volcanic Belt. The Xitle eruption was dominantly effusive, with fluid lava flows accounting for 95% of the total dense rock erupted material (1.1 km³). However, in addition to the initial, Strombolian, cinder cone-building phase, there was a later explosive phase that interrupted effusive activity and deposited three widespread ash fall layers. Major element compositions of olivine-hosted melt inclusions from these ash layers range from 52 to 58 wt.% SiO₂, and olivine host compositions are Fo_84–86. Water concentrations in the melt inclusions are variable (0.2–1.3 wt.% H₂O), with an average of 0.45±0.3 (1σ) wt.% H₂O. Sulfur concentrations vary from below detection (50 ppm) to 1000 ppm but are mostly ≤200 ppm and show little correlation with H₂O. Only the two inclusions with the highest H₂O have detectable CO₂ (310–340 ppm), indicating inclusion entrapment at higher pressures (700–900 bars) than for the other inclusions (≤80 bars). The low and variable H₂O and S contents of melt inclusions combined with the absence of less soluble CO₂ indicates shallow-level degassing before olivine crystallization and melt inclusion formation. Olivine morphologies are consistent with the interpretation that most crystallization occurred rapidly during near-surface H₂O loss. During cinder cone eruptions, the switch from initial explosive activity to effusive eruption probably occurs when the ascent velocity of magma becomes slow enough to allow near-complete degassing of magma at shallow depths within the cone as a result of buoyantly rising gas bubbles. This allows degassed lavas to flow laterally and exit near the base of the cone while gas escapes through bubbly magma in the uppermost part of the conduit just below the crater. The major element compositions of melt inclusions at Xitle show that the short-lived phase of renewed explosive activity was triggered by a magma recharge event, which could have increased overpressure in the storage reservoir beneath Xitle, leading to increased ascent velocities and decreased time available for degassing during ascent.     

Petrologic investigations of the 1995 and 1996 eruptions of Ruapehu volcano, New Zealand: formation of discrete and small magma pockets and their intermittent discharge

 Ruapehu volcano erupted intermittently between September and November 1995, and June and July 1996, producing juvenile andesitic scoria and bombs. The volcanic activity was characterized by small, sequential phreatomagmatic and strombolian eruptions. The petrography and geochemistry of dated samples from 1995 (initial magmatic eruption of 18 September 1995, and two larger events on 23 September and 11 October), and from 1996 (initial and larger eruptions on 17–18 June) suggest that episodes of magma mixing occurred in separate magma pockets within the upper part of the magma plumbing system, producing juvenile andesitic magma by mixing between relatively high (1000–1200  °C)- and low (∼1000  °C)- temperature (T) end members. Oscillatory zoning in pyroxene phenocrysts suggests that repeated mixing events occurred prior to and during the 1995 and 1996 eruptions. Although the 1995 and 1996 andesitic magmas are products of similar mixing processes, they display chronological variations in phenocryst clinopyroxene, matrix glass, and whole-rock compositions. A comparison of the chemistry of magnesian clinopyroxene in the four tephras indicates that, from 18 September through June 1996, the tephras were derived from at least two discrete high-temperature (high-T) batches of magma. Crystals of magnesian clinopyroxene in the 23 September and 11 October tephras appear to be derived from different high-T magma batches. Whole-rock and matrix-glass compositions of all tephras are consistent with their derivation from distinct mixed melts. We propose that, prior to 1995 there was a shallow low-temperature (low-T) magma storage system comprising crystal-rich mush and remnant magma from preceding eruptive episodes. Crystal clots and gabbroic inclusions in the tephras attest to the existence of relict crystal mush. At least two discrete high-T magmas were then repeatedly injected into the mush zone, forming discrete and mixed magma pockets within the shallow system. The intermittent 1995 and 1996 eruptions sequentially tapped these magma pockets. Received: 1 April 1998 / Accepted: 22 December 1998     

Numerical model of convection in sills

The numerical model of convection in magma sills is developed. The model is based on a full system of equations of fluid dynamics and includes heat transfer, buoyancy effects and diffusion of some minor component (marker). Solidification is treated as a phase transition. The results indicate that there are some qualitative differences between very thin sills with Rayleigh number Ra = 10⁵ and thin sills with Ra = 10⁶. For a basaltic magma the first case corresponds to the thickness of the sills of approximately 30 cm and the second case corresponds to the thickness of 60 cm. In the first case mixing is inefficient and conduction is the dominant form of heat transfer. In the second case mixing is efficient and convection is the dominant form of heat transfer. Some of the results can be scaled for the more viscous magmas in thicker sills.     

采用有限元等效体力方法数值分析火山区岩浆压力变形源—以长白山火山为例

火山区岩浆压力变形源的反演计算采用解析方法存在难以考虑地形的限制,采用传统有限元方法则存在网格依赖和计算量大的问题,反演过程中每一次正演由于岩浆房位置和大小变化都需要重新生成一次网格,耗费巨大的计算量和网格生成时间.为了克服上述问题,首次在长白山火山区使用"有限元等效体力"方法考虑地形影响反演地下岩浆压力变形源,计算岩浆应力扰动对周边断层稳定性的影响.在火山区地下压力变形源引起的地表形变计算中,地表地形影响不可忽略.埋深越浅,地表最大径向位移ur所在的位置越靠近岩浆囊中心.当坡度达到30°时,最大垂向位移uz所在位置不再位于岩浆囊正上方.椭球状岩浆囊压力源可以较好地模拟长白山火山地区2002—2003年间的GPS和水准测量.岩浆房扰动应力场和区域构造应力场的叠加有可能造成天池西部近EW向,天池北部以NW-NNW向为主的现今应力方向.岩浆房压力源引起的库仑应力变化有利于天池火山口NW向震群在空间上主要分布于火山口的西南和东北部.     

采用有限元等效体力方法数值分析火山区岩浆压力变形源——以长白山火山为例

火山区岩浆压力变形源的反演计算采用解析方法存在难以考虑地形的限制,采用传统有限元方法则存在网格依赖和计算量大的问题,反演过程中每一次正演由于岩浆房位置和大小变化都需要重新生成一次网格,耗费巨大的计算量和网格生成时间.为了克服上述问题,首次在长白山火山区使用"有限元等效体力"方法考虑地形影响反演地下岩浆压力变形源,计算岩浆应力扰动对周边断层稳定性的影响.在火山区地下压力变形源引起的地表形变计算中,地表地形影响不可忽略.埋深越浅,地表最大径向位移u_r所在的位置越靠近岩浆囊中心.当坡度达到30°时,最大垂向位移u_z所在位置不再位于岩浆囊正上方.椭球状岩浆囊压力源可以较好地模拟长白山火山地区2002—2003年间的GPS和水准测量.岩浆房扰动应力场和区域构造应力场的叠加有可能造成天池西部近EW向,天池北部以NW-NNW向为主的现今应力方向.岩浆房压力源引起的库仑应力变化有利于天池火山口NW向震群在空间上主要分布于火山口的西南和东北部.     

The characteristics of magma reservoir failure beneath a volcanic edifice

Eruptions fed from subsurface reservoirs commonly construct volcanic edifices at the surface, and the growth of an edifice will in turn modify the subsurface stress state that dictates the conditions under which subsequent rupture of the inflating reservoir can occur. We re-examine this problem using axisymmetric finite element models of ellipsoidal reservoirs beneath conical edifices, explicitly incorporating factors (e.g., full gravitational loading conditions, an elastic edifice instead of a surface load, reservoir pressures sufficient to induce tensile rupture) that compromise previous solutions to illustrate why variations in rupture behavior can occur. Relative to half-space model results, the presence of an edifice generally rotates rupture toward the crest of a spherical reservoir, with increasing flank slope (for an edifice of constant volume) and larger edifices (or greater reservoir scaled depths) normally serving to enhance this trend. When non-spherical reservoirs are considered, the presence of an edifice amplifies previously identified half-space failure characteristics, shifting rupture to the crest more rapidly for prolate reservoirs while forcing rupture closer to the midpoint of oblate reservoirs. Rupture is always observed to occur in the σ_t orientation, and depending on where initial failure occurs rupture favors the initial emplacement of either lateral sills, circumferential intrusions or vertically ascending dikes. Ultimately, integration of our numerical model results with other information, for instance the sequence of intrusion/eruption events observed at a given volcano, can provide useful new insight into how a volcano's subsurface magma plumbing system evolved. We demonstrate this process through application of our model to Summer Coon, a well-studied stratocone on Earth, and Ilithyia Mons, a large conical shield volcano on Venus.     

Erratum to Compositional variations and magma mixing in the 1991 eruptions of Hudson volcano,Chile

The August 1991 eruptions of Hudson volcano produced ~2.7 km³ (dense rock equivalent, DRE) of basaltic to trachyandesitic pyroclastic deposits, making it one of the largest historical eruptions in South America. Phase 1 of the eruption (P1, April 8) involved both lava flows and a phreatomagmatic eruption from a fissure located in the NW corner of the caldera. The paroxysmal phase (P2) began several days later (April 12) with a Plinian-style eruption from a different vent 4 km to the south-southeast. Tephra from the 1991 eruption ranges in composition from basalt (phase 1) to trachyandesite (phase 2), with a distinct gap between the two erupted phases from 54–60 wt% SiO₂. A trend of decreasing SiO₂ is evident from the earliest part of the phase 2 eruption (unit A, 63–65 wt% SiO₂) to the end (unit D, 60–63 wt% SiO₂). Melt inclusion data and textures suggest that mixing occurred in magmas from both eruptive phases. The basaltic and trachyandesitic magmas can be genetically related through both magma mixing and fractional crystallization processes. A combination of observed phase assemblages, inferred water content, crystallinity, and geothermometry estimates suggest pre-eruptive storage of the phase 2 trachyandesite at pressures between ~50–100 megapascal (MPa) at 972 ± 26°C under water-saturated conditions (log fO₂ –10.33 (±0.2)). It is proposed that rising P1 basaltic magma intersected the lower part of the P2 magma storage region between 2 and 3 km depth. Subsequent mixing between the two magmas preferentially hybridized the lower part of the chamber. Basaltic magma continued advancing towards the surface as a dyke to eventually be erupted in the northwestern part of the Hudson caldera. The presence of tachylite in the P1 products suggests that some of the magma was stalled close to the surface (<0.5 km) prior to eruption. Seismicity related to magma movement and the P1 eruption, combined with chamber overpressure associated with basalt injection, may have created a pathway to the surface for the trachyandesite magma and subsequent P2 eruption at a different vent 4 km to the south-southeast. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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.

Permeability measurements of natural and experimental volcanic materials with a simple permeameter: Toward an understanding of magmatic degassing processes

Permeability measurement of quenched volcanic porous materials is an important approach to understand permeability development and degassing of vesicular silicic magmas. In this study, we developed a gas permeameter to measure permeability of natural samples and experimental products. The permeameter has broad measurement ranges of pressure difference (10¹–10⁵ Pa) and gas-flow rate (10^− 9–10^− 5 m³/s). These ranges enable us to measure viscous permeability in the range of 10^− 17–10^− 9 m² for 1 centimeter-scale samples, using the Forchheimer equation, which includes the inertial effect of gas flow permeating through the samples. In addition, we improved the procedure for performing permeability measurements of mm-sized products of decompression experiments. Although a previous study reported the first permeability data for vesicular silicic glass products of decompression experiments, we found an overestimation in their permeability data due to problems in sample preparation, especially for very low permeability samples. Our improved measurements give lower permeability values than those of Takeuchi et al. (2005)(Takeuchi, S., Nakashima, S., Tomiya, A., Shinohara, H., 2005. Experimental constraints on the low gas permeability of vesicular magma during decompression. Geophys. Res. Lett., 32, L10312 doi:10.1029/2005GL022491).