Magma plumbing system beneath Ontake Volcano, central Japan

Ontake Volcano in central Japan was last active from ~ 100–35 Ka. The eruptions contained rhyodacite pumice and lavas in the first stage (stage O1, > 33 km³), followed by eruptions of andesite lavas and pyroclastics (stages O2 and O3, > 16 km³). Modeling of major and incompatible elements with Sr isotope ratios suggests that the primary magma was a high-alumina basalt. One andesite magma type appears to have evolved from the basalt in a closed system magma chamber, in part by fractional crystallization, and its generation included crustal assimilation. The other andesite magma type is considered to have evolved in an open system magma chamber in which repeated input of primary magma occurred together with wall-rock assimilation and fractional crystallization. The rhyodacite is inferred to have evolved in a closed system magma chamber by fractional crystallization of the second type of andesite. These genetic relationships require that the magma chamber functioned alternately as an open and a closed system. Geobarometry indicates that there may have been multiple magma chambers, located in the upper crust for the rhyodacite, near the upper–lower crust interface for the andesite and in the mid-lower crust for the basalt. These chambers were stacked to form the magma plumbing system of Ontake. Incompatible element compositions of the basalt are considered to have changed during the eruptions, suggesting that two different plumbing systems for stage O1 magma and for stages O2, O3 magmas existed during the 65 Ka of activity. Evolutionary history of the systems implies that the primary magma was introduced into the magma plumbing system each for ~ 17 500 years and that the life span of a magma plumbing system was shorter than 40 Ka.     

Dissecting large earthquakes in Japan: Role of arc magma and fluids

We synthesized information from recent high-resolution tomographic studies of large crustal earthquakes which occurred in the Japanese Islands during 1995–2008. Prominent anomalies of low-velocity and high Poisson's ratio are revealed in the crust and uppermost mantle beneath the mainshock hypocenters, which may reflect arc magma and fluids that are produced by a combination of subducting slab dehydration and corner flow in the mantle wedge. Distribution of 164 crustal earthquakes ( M 5.7–8.0) that occurred in Japan during 1885–2008 also shows a correlation with the distribution of low-velocity zones in the crust and uppermost mantle. A qualitative model is proposed to explain the geophysical observations recorded so far in Japan. We consider that the nucleation of a large earthquake is not entirely a mechanical process, but is closely related to the subduction dynamics and physical and chemical properties of materials in the crust and upper mantle; in particular, the arc magma and fluids.     

Deep seismic sounding investigation into the deep structure of the magma system in Changbaishan-Tianchi volcanic region

IntroductionThe Changbaishan-Tianchi volcano, located in Jilin Province, northeast of China, is one of the well-known active volcanoes in China. In the year 1 215(15, the volcano underwent an eruption that is believed to be one of the greatest explosive eruptions on earth in the past 2000 years with regard to its scale and violence. The eruption amounted to over 100 km3, and the eruption cylinder reached a height of more than 30 km, the volcanic ashes floated and fell south of Hokkaido, Jap…     

岩浆驱动破裂的形态研究

轴地壳岩浆房是活动扩张中心海洋地壳结构的一个重要组成部分，轴地壳岩浆房通过深部岩浆源的补充，内部岩浆的同化熔融、结晶分异等岩浆房过程，其中的岩浆会破裂上覆的岩石层形成岩浆破裂，并沿岩浆破裂继续向上迁移。本文建立了岩浆迁移的层流模型，从理论上对岩浆迁移问题进行了探讨，并将遗传算法引入到该问题中来，用遗传算法求解了描述岩浆驱动破裂传播的积分方程。如果假设岩浆破裂在远离破裂末梢处的权限宽度为1M，则靠近末梢，破裂的宽度逐渐加大，在末梢处宽度为2m左右。并根据文献对岩浆流体的一些观测参数计算得出，岩浆破裂权限宽度不会很大，一般在1m左右的量级。     

Nishinoshima volcano in the Ogasawara Arc: New continent from the ocean?

Nishinoshima, a submarine volcano in the Ogasawara Arc, approximately 1 000 km south of Tokyo, Japan, suddenly erupted in November 2013, after 40 years of dormancy. Olivine‐bearing phenocryst‐poor andesites found in older submarine lavas from the flanks of the volcano have been used to develop a model for the genesis of andesitic lavas from Nishinoshima. In this model, primary andesite magmas originate directly from the mantle as a result of shallow and hydrous melting of plagioclase peridotites. Thus, it only operates beneath Nishinoshima and submarine volcanoes in the Ogasawara Arc and other oceanic arcs, where the crust is thin. The primary magma compositions have changed from basalt, produced at considerable depth, to andesite, produced beneath the existing thinner crust at this location in the arc. This reflects the thermal and mechanical evolution of the mantle wedge and the overlying lithosphere. It is suggested that continental crust‐like andesitic magma builds up beneath submarine volcanoes on thin arc lithosphere today, and has built up beneath such volcanoes in the past. Andesites produced by this shallow and hydrous melting of the mantle could accumulate through collisions of plates to generate continental crust.     

滇西地区壳幔解耦与腾冲火山区岩浆活动的深部构造研究

根据青藏东部边缘的深部地球物理资料,分析了滇西地区壳幔耦合和腾冲火山区岩浆活动的深部构造特征,确认了地幔各向异性与上地幔速度结构(包括P波速度和S波速度)的内在联系,指出产生这一结果的原因与以腾冲火山区为中心的地幔热物质上涌有关:上地幔顶部平均温度升高导致介质强度降低,在印支块体的侧向挤压或印缅块体的向东俯冲作用下发生韧性变形,造成滇西地区地幔各向异性的快波方向与青藏东部地壳块体的旋转方向不一致.此外,鉴于中下地壳低速层的横向非均匀性,估计韧性流动并非贯通青藏高原的东部边缘,而是被不同的构造块体和边界断裂限定在局部地区.总体而言,滇西地区下地壳的地震波速度和电阻率偏低,具备发生韧性变形的构造条件.作为地壳和上地幔之间的解耦层,它使得青藏东部地壳块体旋转产生的构造应力未能传输至上地幔.腾冲火山区的地壳结构与不同时期的岩浆活动有关,火山区东侧的高速结构代表了上新世时期火山通道内冷凝固结的岩浆侵入体或难以挥发的高密度残留物质,火山区西侧的低速结构反映了更新世以来持续至今的岩浆活动,壳内岩浆源主要分布在10～20km的深度范围内,横向尺度约为15～20km,有可能通过地壳深部的断裂与上地幔岩浆源区相连,估计腾冲火山区下方的岩浆活动将持续进行.     

Subsurface Temperature Changes Due to the Crustal Magmatic Activity – Numerical Simulation

Geothermal aspects of the hypothesis, relating the earthquake swarms in the West Bohemia/Vogtland seismoactive region to magmatic activity, are addressed. A simple 1-D geothermal model of the crust was used to assess the upper limit of the subsurface heating caused by magma intrusion at the assumed focal depth of 9 km. We simulated the process by solving the transient heat conduction equation numerically, considering the heat of magma crystallization to be gradually released in the temperature interval 1100°C to 900°C. The temperature field prior to the intrusion was in steady-state with a surface temperature of 10°C and heat flow of 80 mWm ^–2 , the temperature at the 9 km depth was 270°C. The results suggest that the temperature and heat flow in the uppermost 1 km of the crust begin to grow 100 ka after the intrusion emplacement only, and that the amplitudes of the changes for the realistic lateral extent (a few kilometres) of the intrusion are very small. It was also found that the rate of magma solidification depends strongly on the thickness of the intrusion. It takes about 100 years for a 50 m thick sill to cool down from 1100°C to 600°C, which value represents the lower limit of the solidus temperature. The same cooling takes only 60 days if the sill is 2 m thick. If the nature of the strongly reflected boundaries, interpreted from the January 1997 Nový Kostel seismograms, is connected with the fresh emplacement of magma, the calculated cooling rates have a predictive potential for the temporal changes of the waveforms.     

Petrochemistry of the Gjálp-1996 subglacial eruption, Vatnajökull, SE Iceland

In October 1996 a subglacial fissure to the north of the Grimsvötn caldera in W-Vatnajökull produced about 0.4 km³ of Fe-rich basaltic andesite–icelandite—in an area characterized mostly by tholeiitic basalt. In this paper the chemical composition of volcanic systems in the region is discussed with the help of six new analyses and others from the literature, and a tentative model for their evolution is proposed, in which magma produced by the partial melting of a two-component mantle mixes with hydrous, silicic melt in the crust. The Vatnajökull 1996 magma belongs to a separate volcanic system, intermediate between Bardarbunga and Grimsvötn.     

腾冲火山区水平形变初探

介绍了采用精密大地测量的方法,对腾冲火山区域的水平形变情况进行了探索。通过两期相距３８年的实测数据,对火山“休眠期”中岩浆活动与地壳局部变形的相关性进行研究,并提供了部分计算结果。作者认为,采用精密大地测量的方法对火山区域地壳变形进行研究是可行的;火山区域水平形变与构造形变的区别是存在的;腾冲火山区内近几十年内岩浆的动也是明显的,但目前总体为间歇性收缩期;从水平形变测量的结果来看,显示了南面活跃的迹象;且认为可能是岩浆囊的存在才会引起此大范围的地壳活动。     

Mechanical decoupling and thermal structure at the East Pacific Rise axis 9°N: Constraints from axial magma chamber geometry and seafloor structures

We study the relationships between the seafloor structures and the axial magma chamber geometry in the 9°N overlapping spreading center (OSC) area on the fast spreading East Pacific Rise (EPR). Our observations are based on a new high resolution bathymetric map of the 9°N OSC area derived from picks of the seafloor arrival on 3D seismic data, and on previously published data that constrain the presence and distribution of melt below the 9°N OSC. Differences in the orientation of structures between the seafloor and the magma chamber indicate a sharp change in principal stress directions with depth, suggesting that the brittle crust above the melt sill is decoupled from the melt sill itself and the ductile crust underlying it. The stress-field within the brittle upper crust results from a local interaction of the two overlapping spreading centers, whereas the stress-field in the crust below the melt sill corresponds to the regional stress-field imposed by plate separation. Given this mechanical structure of the crust, the melt sill shape and location appear to be controlled by the following factors: the location of the deep melt source below the melt sill, the ambient stress-field at the depth of the melt sill, and the stress-field in the brittle upper crust above the melt sill, which thermally shapes the roof of the melt sill through repeated eruptions.     

火山区地表形变及应力分布的有限元模拟——以长白山天池火山为例

地壳的非弹性变形以及岩浆与岩石变形的关系是近年构造物理研究的一个热点。文中主要研究火山区的地表形变以及应力场问题,设计了含有岩浆囊的上地壳有限元模型。该模型计算结果与Mogi模型的对比表明,有限元模型不仅可以实现Mogi模型的结果①,而且有很多优势:主要表现在地表对不同岩浆囊形状的响应;区域平面应力场和剖面应力场的体现以及应力场、位移随时间的变化趋势。利用该模型对比分析了岩浆囊半径差异以及岩浆囊深度差异导致的不同形变结果,认为在上地壳浅部,当岩浆囊半径小到一定程度(如0.5km),地表垂向位移的变化受岩浆囊压力变化的影响较小。同时,初步分析了区域最大应力场的分布与实测火山地震分布的关系,发现有较好的对应关系;模拟对比了只有构造应力场和带有岩浆上拱的构造应力场,发现后者导致了最大剪应力场的逐渐增大     

Reservoir-forming features of abiotic origin gas in Songliao Basin

The vertical structure of the crustal block of the Songliao Basin can be divided into upper, middle and low Earth’s crust according to density. There is an about 3-km-thick low density interval between the upper crust and the middle crust. This interval may be a magma chamber accumulated in crust by “fluid phase” which is precipitated and separated from upper mantle meltmass. The abiogenetic natural gas, other gaseous mass and hydrothermal fluids are provided to the Songliao rifted basin through crustal faults and natural earthquakes. This is a basic condition to form an abiogenetic gas reservoir in the Songliao Basin. On both flanks of the upper crust (or named basin basement) fault there are structural traps in and above the basement and unconformity surface or lateral extended sand, which contains communicated pores, as migration pathway and natural gas reservoir; up to gas reservoirs there is shale as enclosed cap rock, and the suitable arrangement of these conditions is the basic features of abiogenetic gas reservoir. Project supported by the National Natural Sc~ence Foundation of China.     

The 10 ka multiple vent pyroclastic eruption sequence at Tongariro Volcanic Centre, Taupo Volcanic Zone, New Zealand: Part 2. Petrological insights into magma storage and transport during regional extension

The six eruption episodes of the 10 ka Pahoka–Mangamate (PM) sequence (see companion paper) occurred over a ?200–400-year period from a 15-km-long zone of multiple vents within the Tongariro Volcanic Centre (TgVC), located at the southern end of the Taupo Volcanic Zone (TVZ). Most TgVC eruptives are plagioclase-dominant pyroxene andesites and dacites, with strongly porphyritic textures indicating their derivation from magmas that ascended slowly and stagnated at shallow depths. In contrast, the PM pyroclastic eruptives show petrographic features (presence of phenocrystic and groundmass hornblende, and the coexistence of olivine and augite without plagioclase during crystallisation of phenocrysts and microphenocrysts) which suggest that their crystallisation occurred at depth. Depths exceeding 8 km are indicated for the dacitic magmas, and >20 km for the andesitic and basaltic andesitic magmas. Other petrographic features (aphyric nature, lack of reaction rims around hornblende, and the common occurrence of skeletal microphenocrystic to groundmass olivine in the andesites and basaltic andesites) suggest the PM magmas ascended rapidly immediately prior to their eruption, without any significant stagnation at shallow depths in the crust. The PM eruptives show three distinct linear trends in many oxide–oxide diagrams, suggesting geochemical division of the six episodes into three chronologically-sequential groups, early, middle and late. Disequilibrium features on a variety of scales (banded pumice, heterogeneous glassy matrix and presence of reversely zoned phenocrysts) suggest that each group contains the mixing products of two end-member magmas. Both of these end-member magmas are clearly different in each of the three groups, showing that the PM magma system was completely renewed at least three times during the eruption sequence. Minor compositional diversity within the eruptives of each group also allows the PM magmas to be distinguished in terms of their source vents. Because petrography suggests that the PM magmas did not stagnate at shallow levels during their ascent, the minor diversity in magmas from different vents indicates that magmas ascended from depth through separate conduits/dikes to erupt at different vents either simultaneously or sequentially. These unique modes of magma transport and eruption support the inferred simultaneous or sequential tapping of small separate magma bodies by regional rifting in the southern Taupo Volcanic Zone during the PM eruption sequence (see companion paper).     

Magma transfer processes at persistently active volcanoes: insights from gravity observations

Magma transfer processes at persistently active volcanoes are distinguished by the large magma flux required to sustain the prodigious quantities of heat and gas emitted at the surface. Although the resulting degassed magma has been conjectured to accumulate either deep within the volcanic edifice or in the upper levels of the sub-edifice system, no direct evidence for such active accumulation has been reported. Temporal gravity data are unique in being able to quantify mass changes and have been successfully used to model shallow magma movements on different temporal scales, but have not generally been applied to the investigation of postulated long-term accumulation of magma at greater spatial scales within volcanic systems. Here, we model the critical data acquisition parameters required to detect mass flux at volcanoes, we review existing data from a number of volcanoes that exemplify the measurement of shallow mass changes and present new data from Poas and Telica volcanoes. We show that if a substantial proportion of degassed magma lodges within the sub-edifice region, it would result in measurable annual to decadal gravity increases occurring over spatial scales of tens of kilometres and propose that existing microgravity data from Sakurajima and, possibly, Etna volcanoes could be interpreted in these terms. Furthermore, such repeat microgravity data could be used to determine whether the accumulation rate is in equilibrium with the rate of production of degassed magma as calculated from the surface gas flux and hence identify the build-up of gas-rich magma at depth that may be significant in terms of eruption potential. We also argue that large magma bodies, both molten and frozen, modelled beneath volcanoes from seismic and gravity data, could represent endogenous or cryptic intrusions of degassed magma based on order of magnitude calculations using present-day emission rates and typical volcano lifetimes.     

2-D crustal Poisson's ratio from seismic travel time inversion in Changbaishan Tianchi volcanic region

Based on the inversion method of 2D velocity structure and interface, the crustal velocity structures of P-wave and S-wave along the profile L1 are determined simultaneously with deep seismic sounding data in Changbaishan Tianchi volcanic region, and then its Poisson's ratio is obtained. Calculated results show that this technique overcomes some defects of traditional forward calculation method, and it is also very effective to determine Poisson's ratio distribution of deep seismic sounding profile, especially useful for study on volcanic magma and crustal fault zone. Study result indicates that there is an abnormally high Poisson's ratio body that is about 30 km wide and 12 km high in the low velocity region under Tianchi crater. Its value of Poisson's ratio is 8% higher than that of surrounding medium and it should be the magma chamber formed from melted rock with high temperature. There is a high Poisson's ratio zone ranging from magma chamber to the top of crust, which may be the uprise passage of hot substance. The lower part with high Poisson's ratio, which stretches downward to Moho, is possibly the extrusion way of hot substance from the uppermost mantle. The conclusions above are consistent with the study results of both tomographic determination of 3D crustal structure and magnetotelluric survey in this region.     

A nonstationary model of the thermal regime of axial zones of mid-ocean ridges: Formation of crustal and mantle magma chambers

A nonstationary model of spreading with periodic intrusions of a molten material into an axial zone of a mid-ocean ridge (MOR) is applied to numerical analysis of the thermal state in MOR axial zones and the formation of crustal and mantle magma chambers in them. The model satisfactorily explains the positions, dimensions, and shapes of magma chambers, as well as variations in these parameters depending on the spreading rate, temperature, and composition of crustal and mantle rocks. The release and absorption of the latent heat of rock melting, hydrothermal heating of the crust, and variations in the solidus and liquidus temperatures of crustal and mantle rocks as a function of their composition are factors controlling the shape and position of crustal magma chambers.     

Petrogenesis of the Mesozoic intrusive rocks in the Tongling area,Anhui Province,China and their constraint on geodynamic process

TheTonglingarea,whichiscalledtheChineseCopperCapital,isoneofthemostimportantnon-ferrousmetalproducersinChina(e.g.Cu,AuandAg,especiallyCu).ManyresearchershavenotedthatthemetaldepositsarecloselyrelatedtotheMesozoicintrusiverocksinthisarea.Therefore,theTongl…     

云南腾冲火山区上部地壳三维地震速度层析成像

根据1999年腾冲火山区人工地震探测资料，用层析成像方法反演了腾冲火山区上地壳三维P波速度结构．层析成像结果显示出，在腾冲火山区下方的上地壳内存在明显的P波低速区，该低速区在地下7～8km深度，可能与上地壳内的岩浆囊或部分熔融体有关．层析结果还显示出低速区主体在热海地区的东北侧．地表附近的低速区可能是构造破裂发育地区形成了有利于深部热流体上升的通道，因此，在热海形成最活跃的地热显示．      

Replenishment rates of crustal magma and their bearing on potential sources of thermal energy

A general model of magma intrusion into the crust is developed which is based on a viscous-dissipation, forced-convection flow process driven by gravitational-buoyancy forces. Although some of the points in this general model have been studied before, it is possible with the present model to go further and calculate magma volumetric intrusion rates from fundamental properties and parameters. Equations for forced convection in a conduit with viscous dissipation are combined with results for the temperature dependence of magma viscosity. The volumetric intrusion rate is shown to be not a function of viscosity as might be expected, but rather a function primarily of the rate of change of viscosity with temperature. The model predictions for intrusion rate correlate well with field results for several sites where data exist for both intrusion or extrusion rate and for the temperature-dependent behavior of magma viscosity. The model predicts magma chamber replenishment rates equivalent to thermal energy rates on the order 10 GW (gigawatts) for a single active magma site. Assuming active magma sites on a 50-km spacing along volcanic lineaments leads to an estimate of a renewable magma intrusion rate into the crust of the western U.S. on the order of 2 TW (terawatts).     

火山活动状态及其机理的初步探讨

根据火山喷发实例总结了火山喷发在不同阶段的活动状态,并探讨了可能的物理机理。火山活动从岩浆补给到岩浆喷发的物理过程可分为3个阶段:1)岩浆补给阶段,岩浆囊压力差或过剩压力的大小决定了火山活动是否休眠或扰动,岩浆补给速率对压力差起了决定性的作用;2)通道形成阶段,当过剩压力超过围岩破裂强度时,围岩开始破裂,之后水热活动起了重要的作用;3)岩浆运移与失稳喷发阶段,主要是岩浆运移与地壳盖层的相互作用与失稳的过程。文中还讨论了火山活动状态与火山喷发危险性等级之间的关系,7个危险性等级分别对应于火山活动的7种状态,即休眠、平静、扰动、动荡、临界、活动、灾变