Determining magma flow in sills,dykes and laccoliths and their implications for sill emplacement mechanisms

Three-dimensional seismic data from the Faeroe-Shetland Basin provides detailed information on the relationships between sills, dykes, laccoliths and contemporaneous volcanic activity. The data shows that sills are predominantly concave upwards, being complete or partial versions of radially or bilaterally symmetrical forms that possess flat inner saucers connected to a flat outer rim by a steeply inclined sheet. Such morphologies are only partially modified by pre-existing faults. Sills can be sourced from dykes or the steep climbing portions of deeper sills. Both sills and dykes can provide magma to overlying volcanic fissures and sills can be shown to feed shallow laccoliths. Magma flow patterns, as revealed by opacity rendering, suggest that sills propagate upwards and outwards away from the magma feeder. As an individual sill can consist of several leaves emplaced at different stratigraphic levels, and as a sill or dyke can provide magma to volcanic fissures, other sills and laccoliths, the data suggests that neutral buoyancy concepts may not provide a complete explanation for the mechanism and level of sill emplacement. Instead, the data suggests that the presence of lithological contrasts, particularly ductile horizons such as overpressured shales may permit sill formation at any level below the neutrally buoyant level. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Ken Thomson–deceased, April 2007     

Magmatic gas scrubbing: implications for volcano monitoring

Despite the abundance of SO_2(g) in magmatic gases, precursory increases in magmatic SO_2(g) are not always observed prior to volcanic eruption, probably because many terrestrial volcanoes contain abundant groundwater or surface water that scrubs magmatic gases until a dry pathway to the atmosphere is established. To better understand scrubbing and its implications for volcano monitoring, we model thermochemically the reaction of magmatic gases with water. First, we inject a 915°C magmatic gas from Merapi volcano into 25°C air-saturated water (ASW) over a wide range of gas/water mass ratios from 0.0002 to 100 and at a total pressure of 0.1 MPa. Then we model closed-system cooling of the magmatic gas, magmatic gas-ASW mixing at 5.0 MPa, runs with varied temperature and composition of the ASW, a case with a wide range of magmatic–gas compositions, and a reaction of a magmatic gas–ASW mixture with rock. The modeling predicts gas and water compositions, and, in one case, alteration assemblages for a wide range of scrubbing conditions; these results can be compared directly with samples from degassing volcanoes. The modeling suggests that CO_2(g) is the main species to monitor when scrubbing exists; another candidate is H₂S_(g), but it can be affected by reactions with aqueous ferrous iron. In contrast, scrubbing by water will prevent significant SO_2(g) and most HCl_(g) emissions until dry pathways are established, except for moderate HCl_(g) degassing from pH<0.5 hydrothermal waters. Furthermore, it appears that scrubbing will prevent much, if any, SO_2(g) degassing from long-resident boiling hydrothermal systems. Several processes can also decrease or increase H_2(g) emissions during scrubbing making H_2(g) a poor choice to detect changes in magma degassing.We applied the model results to interpret field observations and emission rate data from four eruptions: (1) Crater Peak on Mount Spurr (1992) where, except for a short post-eruptive period, scrubbing appears to have drastically diminished pre-, inter-, and post-eruptive SO_2(g) emissions, but had much less impact on CO_2(g) emissions. (2) Mount St. Helens where scrubbing of SO_2(g) was important prior to and three weeks after the 18 May 1980 eruption. Scrubbing was also active during a period of unrest in the summer of 1998. (3) Mount Pinatubo where early drying out prevented SO_2(g) scrubbing before the climactic 15 June 1991 eruption. (4) The ongoing eruption at Popocatépetl in an arid region of Mexico where there is little evidence of scrubbing.In most eruptive cycles, the impact of scrubbing will be greater during pre- and post-eruptive periods than during the main eruptive and intense passive degassing stages. Therefore, we recommend monitoring the following gases: CO_2(g) and H₂S_(g) in precursory stages; CO_2(g), H₂S_(g), SO_2(g), HCl_(g), and HF_(g) in eruptive and intense passive degassing stages; and CO_2(g) and H₂S_(g) again in the declining stages. CO_2(g) is clearly the main candidate for early emission rate monitoring, although significant early increases in the intensity and geographic distribution of H₂S_(g) emissions should be taken as an important sign of volcanic unrest and a potential precursor. Owing to the difficulty of extracting SO_2(g) from hydrothermal waters, the emergence of >100 t/d (tons per day) of SO_2(g) in addition to CO_2(g) and H₂S_(g) should be taken as a criterion of magma intrusion. Finally, the modeling suggests that the interpretation of gas-ratio data requires a case-by-case evaluation since ratio changes can often be produced by several mechanisms; nevertheless, several gas ratios may provide useful indices for monitoring the drying out of gas pathways.     

Development and sensitivity analysis of a model for assessing stratification and safety of Lake Nyos during artificial degassing

To prevent the recurrence of a disastrous eruption of carbon dioxide (CO₂) from Lake Nyos, a degassing plan has been set up for the lake. Since there are concerns that the degassing of the lake may reduce the stability of the density stratification, there is an urgent need for a simulation tool to predict the evolution of the lake stratification in different scenarios. This paper describes the development of a numerical model to predict the CO₂ and dissolved solids concentrations, and the temperature structure as well as the stability of the water column of Lake Nyos. The model is tested with profiles of CO₂ concentrations and temperature taken in the years 1986 to 1996. It reproduces well the general mixing patterns observed in the lake. However, the intensity of the mixing tends to be overestimated in the epilimnion and underestimated in the monimolimnion. The overestimation of the mixing depth in the epilimnion is caused either by the parameterization of the k-epsilon model, or by the uncertainty in the calculation of the surface heat fluxes. The simulated mixing depth is highly sensitive to the surface heat fluxes, and errors in the mixing depth propagate from one year to the following. A precise simulation of the mixolimnion deepening therefore requires high accuracy in the meteorological forcing and the parameterization of the heat fluxes. Neither the meteorological data nor the formulae for the calculation of the heat fluxes are available with the necessary precision. Consequently, it will be indispensable to consider different forcing scenarios in the safety analysis in order to obtain robust boundary conditions for safe degassing. The input of temperature and CO₂ to the lake bottom can be adequately simulated for the years 1986 to 1996 with a constant sublacustrine source of 18 l s^–1 with a CO₂ concentration of 0.395 mol l^–1 and a temperature of 26 °C. The results of this study indicate that the model needs to be calibrated with more detailed field data before using it for its final purpose: the prediction of the stability and the safety of Lake Nyos during the degassing process.Responsible Editor: Hans Burchard     

Pit crater structure and processes governing persistent activity at Masaya Volcano, Nicaragua

 Persistent activity at Masaya Volcano, Nicaragua, is characterised by cycles of intense degassing, lava lake development and pit crater formation. It provides a useful site to study the processes which govern such activity, because of its easy accessibility and relatively short cycles (years to decades). An understanding of the present activity is important because Masaya is visited by large numbers of tourists, is located close to major cities and has produced voluminous lavas, plinian eruptions and ignimbrites in the recent past. We provide structural and geophysical data that characterise the "normal" present state of activity. These indicate that the ongoing degassing phase (1993 to present) was not caused by fresh magma intrusion. It was associated with shallow density changes within the active Santiago pit crater. The activity appears to be associated predominantly with shallow changes in the pit crater structure. More hazardous activity will occur only if there are significant departures from the present gravity, deformation and seismic signatures. Received: 16 May 1997 / Accepted: 29 October 1997     

Modelling the dynamics and thermodynamics of volcanic degassing

 The rates of passive degassing from volcanoes are investigated by modelling the convective overturn of dense degassed and less dense gas-rich magmas in a vertical conduit linking a shallow degassing zone with a deep magma chamber. Laboratory experiments are used to constrain our theoretical model of the overturn rate and to elaborate on the model of this process presented by Kazahaya et al. (1994). We also introduce the effects of a CO₂–saturated deep chamber and adiabatic cooling of ascending magma. We find that overturn occurs by concentric flow of the magmas along the conduit, although the details of the flow depend on the magmas' viscosity ratio. Where convective overturn limits the supply of gas-rich magma, then the gas emission rate is proportional to the flow rate of the overturning magmas (proportional to the density difference driving convection, the conduit radius to the fourth power, and inversely proportional to the degassed magma viscosity) and the mass fraction of water that is degassed. Efficient degassing enhances the density difference but increases the magma viscosity, and this dampens convection. Two degassing volcanoes were modelled. At Stromboli, assuming a 2 km deep, 30% crystalline basaltic chamber, containing 0.5 wt.% dissolved water, the ∼700 kg s^–1 magmatic water flux can be modelled with a 4–10 m radius conduit, degassing 20–100% of the available water and all of the 1 to 4 vol.% CO₂ chamber gas. At Mount St. Helens in June 1980, assuming a 7 km deep, 39% crystalline dacitic chamber, containing 4.6 wt.% dissolved water, the ∼500 kg s^–1 magmatic water flux can be modelled with a 22–60 m radius conduit, degassing ∼2–90% of the available water and all of the 0.1 to 3 vol.% CO₂ chamber gas. The range of these results is consistent with previous models and observations. Convection driven by degassing provides a plausible mechanism for transferring volatiles from deep magma chambers to the atmosphere, and it can explain the gas fluxes measured at many persistently active volcanoes. Received: 26 September 1997 / Accepted: 11 July 1998     

Magma fragmentation dynamics: experiments with analogue porous low-strength material

 Experiments were conducted on the fragmentation of analogue low-strength porous material (plastiprin) by rapid decompression in a shock-tube-type apparatus. The porous samples (length=365 mm, cross-section dimensions 40×40 mm) pressurized by air to pressures up to 0.9 MPa, were rapidly decompressed to 0.1 MPa. Rapid decompression of samples caused fragmentation and ejection of the fragmentation products into a large volume tank. The process of analogue material fragmentation was documented using high-speed cinematography and dynamic pressure measurements. The duration of the fragmentation event is significantly shorter than that of the ejection event. The fragmentation of material precedes the acceleration of fragments. As a result of fragmentation, sub-parallel fractures are generated. The characteristic fragment size decreases as the initial pressure differential increases. The ejected fragments obtain velocities of 60 m/s. The mechanisms of material fragmentation during unloading and fragmentation wave propagation are discussed. The experimental results provide insight into the fragmentation dynamics of highly viscous magmas in which brittle failure at high strain rate is possible. Received: 23 July 1997 / Accepted: 23 November 1997     

新疆希勒库都克铜钼矿区岩浆混合作用: 来自锆石U-Pb年代学的证据

新疆北部富蕴县内希勒库都克铜钼矿区的花岗闪长岩及其包体岩相学、矿物化学和岩石地球化学特征及野外地质特征显示其为岩浆混合作用的结果。本文获得花岗闪长岩及其中暗色微粒包体玄武安山玢岩、细粒辉长闪长岩锆石U-Pb年龄为326.8±2.1Ma、327.6±2.4Ma、329.3±2.3Ma,年龄值基本一致,这一结果从年代学角度为花岗闪长岩及其中暗色微粒包体的岩浆混合作用成因提供证据。偏酸性的花岗闪长质与偏基性玄武安山质岩浆混合作用形成了331.9±2.1Ma的安山玢岩脉。     

冈底斯明则-程巴斑岩-夕卡岩型Mo-Cu矿床成矿时代与含矿岩石成因

冈底斯东南缘克鲁-冲木达矿带以发育渐新世斑岩-夕卡岩型Mo-W-Cu成矿作用而有别于冈底斯中新世斑岩-夕卡岩型Cu-Mo-Au矿带,但是对渐新世成矿岩石类型及成因以及渐新世与中新世斑岩成矿作用存在差异的原因尚不清楚,严重制约了该带矿产资源潜力评价.为此,本文选择矿带内程巴斑岩型Mo矿床及明则夕卡岩型Cu矿床开展了成岩与成矿年代学及含矿岩石成因研究.明则矿区与夕卡岩矿化有关的二长岩的结晶年龄为30.4±0.6Ma,表明夕卡岩矿化时代约30Ma,这与程巴斑岩型Mo矿床矿石中辉钼矿Re-Os同位素分析所揭示的约30Ma矿化时代一致.明则二长岩的SiO2含量为55％ ～57％,K2O含量为2.7％ ～5.0％,属于钾玄岩,具有较高的MgO含量(3.5％～6.9％)、Mg#值(57.6～67.2)以及相容元素含量(Cr =34×10-6 ～379×10-6;Ni为48×10-6～116×10-6).而程巴花岗闪长岩的SiO2含量为65％ ～ 67％,K2O含量为3.2％～4.1％,属于高钾钙碱性岩,具有较低的MgO含量(1.7％～2.1％)、Mg#值(49.5～51.1)以及相容元素含量(Cr =20×10-6～39×10-6;Ni为16×10-6 ～25×10-6).二长岩和花岗闪长岩均显示富集轻稀土元素与大离子亲石元素,亏损高场强元素,具有弱的Eu负异常,但二长岩的REE含量明显高于花岗闪长岩,且具有较低的Sr/Y比值(24 ～49),后者Sr/Y比值为54～68.两种岩石均具有类似的Hf同位素组成,二长岩的εHf(t)值为+2.8～+6.8,花岗闪长岩εHf(t)值为+4.2～ +6.1.综合分析表明,二长岩是受大洋板片流体交代的富集岩石圈地幔部分熔融的产物,而花岗闪长岩是增厚的新生下地壳部分熔融的产物,二者在形成过程中可能发生了幔源镁铁质岩浆与壳源长英质岩浆的混合作用,这种壳幔相互作用导致区域斑岩-夕卡岩型Cu-Mo矿化,但由于该区隆升剥蚀强烈,上部的斑岩型Cu矿体基本被剥蚀而没有保存下来,因此现今表现为斑岩型Mo矿床.     

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.     

长白山天池火山区电性结构和岩浆系统

在长白山火山区及其附近 ,沿 5条伞状测线和 1条北西向测线进行了 6 1个大地电磁测深点的观测。用Robust等方法对观测资料进行了处理 ;应用阻抗张量分解方法消除了局部不均匀体的影响 ;计算了反映地下横向不均匀性的地磁感应矢量。利用RRI二维自动反演等技术对资料进行了解释。结果表明 :在近地表存在低阻体 ,并分别与长白山天池水体、聚龙泉温泉、锦江温泉及长白山山门附近的地下水相对应。在长白山天池及其以北和以东地区 ,约 12km深处存在电阻率很低的地质体 ,电阻率为 10到几十Ω·m ,可能是地壳岩浆囊。在测区范围内 ,地壳电性结构在南北方向存在着明显的横向变化 ,而东西方向相对均匀 ,变化比较平缓