Magma Plumbing System Research Articles

Constraints of texture and composition of clinopyroxene phenocrysts of Holocene volcanic rocks on a magmatic plumbing system beneath Tengchong, SW China

Understanding processes of magma replenishment in a magma plumbing system is essential to predict eruption potential of a dormant volcano. In this study, we present new petrologic and thermobarometric data for youngest lava flows from the Holocene Heikongshan volcano in the Tengchong area, SW China. Clinopyroxene phenocrysts from the trachytic lava flows display various textural/compositional zoning styles (i.e., normal, reverse and oscillatory). Such zoning patterns are indicative of an open magmatic plumbing system with multiphase magma replenishment and mixing, which were likely a key drive of the volcanic eruptions. Thermobarometric calculations of these zoned clinopyroxene phenocrysts yield crystallization pressures of 3.8–7.1 kbar (peak at 4.5–7.0 kbar), corresponding to a magma chamber at depths of 14–21 km. The calculated depths are consistent with the large low-resistivity body at 12–30 km beneath the Heikongshan volcano, implying that the magmatic plumbing system may still be active. Recent earthquakes in the Tengchong area suggest that the regional strike-slip faulting are still active, and may trigger future volcanic eruptions if the magma chamber(s) beneath the Tengchong volcanic field is disturbed, in spite of the volcanic quiescence since 1609 CE.

Phase-equilibrium geobarometers for silicic rocks based on rhyolite-MELTS. Part 4: Plagioclase, orthopyroxene, clinopyroxene, glass geobarometer, and application to Mt. Ruapehu, New Zealand

A new phase equilibria geobarometer determines magmatic storage and crystallization conditions, including pressure, temperature, oxygen fugacity ( $${f_{{{\text{o}}_2}}}$$ ), and the presence of a fluid phase for glass-bearing rocks containing the assemblage plagioclase + pyroxene(s). This newly developed geobarometer can better constrain crystallization conditions of shallow (< 500 MPa; <~ 20 km), glass-bearing andesites to dacites. The geobarometer utilizes rhyolite-MELTS to determine crystallization conditions in natural pumice and scoria samples. The validity of the geobarometer is tested by comparing it to results from experiments. Uncertainties are assessed using Monte Carlo simulations. We apply the geobarometer to the plag + opx + cpx-bearing system of Mt. Ruapehu, in the southern Taupo Volcanic Zone, New Zealand. The samples from Mt. Ruapehu are tested from ~ 5 to ~ 400 MPa and from super-liquidus to 90% crystalline (~ 1200 to ~ 700 °C). Mt. Ruapehu serves as a methodological testing ground for the geobarometer, and results from our geobarometer agree with recent Mt. Ruapehu studies. Results show a distribution of crystallization pressures ranging from 50 to 150 MPa (~ 2.0 to 5.9 km) for different eruptions, with modes of 110 MPa (~ 4.3 km) and 130 MPa (~ 5.1 km). These are consistent with field interpretations of different eruptive styles based on juvenile clast textures and previous knowledge of the magma plumbing system. Mt. Ruapehu magmas are fluid saturated, with $${f_{{{\text{o}}_2}}}$$ of ΔQFM ~ + 1 (NNO).

Nature of the magma storage system beneath the Damavand volcano (N. Iran): An integrated study

Damavand intraplate stratovolcano constructed upon a moderately thick crust (58–67km) over the last 2Ma. The erupted products are dominantly trachyandesite-trachyte (TT) lavas and pyroclasts, with minor mafic magmas including tephrite-basanite-trachybasalt and alkali olivine basalts emplaced as cinder cones at the base of the stratovolcano. The TT products are characterized by a mineral assemblage of clinopyroxene (diopside-augite), orthopyroxene (clinoenstatite), feldspar (An2–58, Ab6–69, Or2–56), high Ti phlogopite, F-apatite, FeTi oxides, and minor amounts of olivine (Fo73–80), amphibole and zircon, whereas olivine (Fo78–88), high Mg# (80–89) diopside, feldspar, apatite and FeTi oxide occur in the mafic magmas. The presence of hydrous and anhydrous minerals, normal zonings, mafic cumulates, and the composition of magmatic inclusions in the TT products suggest evolutionary processes in polybaric conditions. In the same way, disequilibrium textures - including orthopyroxene mantled with clinopyroxene, reaction rim of phlogopite and amphibole, the coexistence of olivine and orthopyroxene, reverse, oscillatory and complex zonings of pyroxene and feldspar crystals - suggest magmatic evolutions in open systems with a varying temperature, oxygen fugacity, water as well as pressure and, to a lesser extent, melt chemistry. Mineral assemblages are used to model the physicochemical conditions and assess default parameters for the thermodynamic simulation of crystallization using MELTS software to track the P-T-H2O-ƒO2 evolution of the magma plumbing system. Thermobarometry and MELTS models estimated the initial nucleation depth at 16–17kb (56–60km) for olivine (Fo89) and high Al diopside crystals occurring in the mafic primary magma; it then stopped and underwent fractionation between 8 and 10kb (28–35km), corresponding with Moho depth, and continued to differentiate in the lower crust, in agreement with the geophysical models. The mafic rocks were formed by crystal fractionation from the reconstructed primary magma (13wt% MgO) with the minor role of recharge and crustal assimilation. Phenocrysts in TT lavas recorded a wide range of temperature and pressure of crystallization; at least three main levels of magma storage can be recognized according to the statistical analysis of the models, at 6–8kb (22–28km), 4–6kb (15–22km), and 0.6–3kb (2–11km), respectively. The temperature of crystallization ranged from 1430 to 1180°C for primary mafic magma to alkali olivine basalts and 1180–800°C for TT suite. According to the current geophysical models, the present structure of the magma storage system in the crust beneath the Damavand volcano consists of three major accumulation zones located at 20km, 6–8km and 3–4.5km depth. Our data enlarge this scenario, suggesting a more complex magma storage system strongly controlled by the transpressional tectonic regime. Multi-depth magma reservoirs may account for the local thickening of crust below the volcano. The polybaric fractionation model, using the MELTS algorithm, reproduces mineralogy and chemical variations of minerals and whole rock of the Damavand TT lavas. However, some discrepancies between major elements of models and trends of data can be ascribed to the recharge of more mafic magma, minor crustal assimilation, disaggregation of crystal-rich mushes and uptake of magmatic inclusions, as well as crystals from different crustal levels. The polybaric differentiation as the fractionation and/or accumulation of crystals was the probable mechanism for explaining the scarcity of mafic volcanic rocks at the Damavand volcano.

The genetic relationship between andesites and dacites at Tungurahua volcano, Ecuador

Volcanic eruptions of intermediary and silica-rich magmas (andesites, dacites and rhyolites) in convergent arc settings generate voluminous and explosive eruptions that can strongly affect human activity and have significant environmental impacts. It is therefore crucial to understand how these magmas are generated in order to anticipate their potential impact. At convergent margins, primitive magmas (primitive basalts and/or andesites) are derived from the mantle wedge and they are progressively modified by physical and chemical processes operating between the melting zone and the surface to produce silica-rich magmas.In order to elucidate the relationship between andesites and dacites, we focus on Tungurahua volcano, located in the Ecuadorian Andes. We collected a set of samples comprising such lithologies that were erupted during the last 3000 year BP. This relatively short period of time allows us to assume that the geodynamic parameters remain constant. Petrology and major-trace element compositions of these lavas have already been examined, and so we performed a complementary Pb-Sr isotope study in order to determine the nature and origin of the components involved in andesite and dacite genesis. Sr isotopes range from 0.70417 to 0.70431, and Pb isotope compositions range from 18.889 to 19.154 for 206Pb/204Pb, from 15.658 to 15.696 for 207Pb/204Pb, and from 38.752 to 38.918 for 208Pb/204Pb. Dacites display a remarkably homogeneous Pb isotopic composition, with higher 206Pb/204Pb values for a given 207-208Pb/204Pb compared to andesites. Andesites show notable 207Pb/206Pb variations for a given SiO2 content, whereas dacites have lower and homogenous 207Pb/206Pb values. Andesite and dacite altogether plot in a roughly triangular distribution, with dacitic magmas systematically plotting at the high SiO2 and 87Sr/86Sr and low 207Pb/206Pb fields. Based on our new dataset, we show that at least 3 different components are required to explain the Tungurahua compositional and isotope variation: one corresponds to the mantle, the second has a deep origin (slab component or lower crust), and a mixture between these two components explains andesite heterogeneity. The third component is derived from the underlying upper continental crust. While andesites are derived from deep components, dacites are derived from the andesitic magmas that underwent an assimilation-fractional crystallization (AFC) process with incorporation of the local metamorphic basement. Finally, we used the geochemical and isotopic data to produce a model of the magmatic plumbing system beneath Tungurahua, consistent with geophysical and experimental petrology constraints. We conclude that melt migration and storage in the upper crust appears to be a key parameter for controlling volcanic behavior though time.

Understanding cyclic seismicity and ground deformation patterns at volcanoes: Intriguing lessons from Tungurahua volcano, Ecuador

Cyclic seismicity and ground deformation patterns are observed on many volcanoes worldwide where seismic swarms and the tilt of the volcanic flanks provide sensitive tools to assess the state of volcanic activity. Ground deformation at active volcanoes is often interpreted as pressure changes in a magmatic reservoir, and tilt is simply translated accordingly into inflation and deflation of such a reservoir. Tilt data recorded by an instrument in the summit area of Tungurahua volcano in Ecuador, however, show an intriguing and unexpected behaviour on several occasions: prior to a Vulcanian explosion when a pressurisation of the system would be expected, the tilt signal declines significantly, hence indicating depressurisation. At the same time, seismicity increases drastically. Envisaging that such a pattern could carry the potential to forecast Vulcanian explosions on Tungurahua, we use numerical modelling and reproduce the observed tilt patterns in both space and time. We demonstrate that the tilt signal can be more easily explained as caused by shear stress due to viscous flow resistance, rather than by pressurisation of the magmatic plumbing system. In general, our numerical models prove that if magma shear viscosity and ascent rate are high enough, the resulting shear stress is sufficient to generate a tilt signal as observed on Tungurahua. Furthermore, we address the interdependence of tilt and seismicity through shear stress partitioning and suggest that a joint interpretation of tilt and seismicity can shed new light on the eruption potential of silicic volcanoes.

High-resolution 40Ar/39Ar geochronology of volcanic rocks from the Siebengebirge (Central Germany)—Implications for eruption timescales and petrogenetic evolution of intraplate volcanic fields

A key parameter in understanding mantle dynamics beneath continents is the temporal evolution of intraplate volcanism in response to lithospheric thinning and asthenospheric uplift. To contribute to a better understanding of how intraplate volcanic fields evolve through time, we present a high precision 40Ar/39Ar age dataset for volcanic rocks from the Siebengebirge volcanic field (SVF) from central Germany, one of the best studied and compositionally most diverse intraplate volcanic fields of the Cenozoic Central European Volcanic Province (CEVP). Petrological and geochemical investigations suggest that the formation of the different rock types that occur in the SVF can be explained by a combination of assimilation and fractional crystallisation processes, starting from at least two different parental magmas with different levels of silica saturation (alkali basaltic and basanitic), and originating from different mantle sources. These evolved along two differentiation trends to latites and trachytes, and to tephrites and tephriphonolites, respectively. In contrast to their petrogenesis, the temporal evolution of the different SVF suites is poorly constrained. Previous K/Ar ages suggested a time of formation between about 28 and 19 Ma for the mafic rocks, and of about 27 to 24 Ma for the differentiated rocks. Our results confirm at high precision that the differentiated lithologies of both alkaline suites (40Ar/39Ar ages from 25.3 ± 0.2 Ma to 25.9 ± 0.3 Ma) erupted contemporaneously within a very short time period of ~ 0.6 Ma, whereas the eruption of mafic rocks (basanites) lasted at least 8 Ma (40Ar/39Ar ages from 22.2 ± 0.2 Ma to 29.5 ± 0.3 Ma). This implies that felsic magmatism in the central SVF was likely a single event, possibly triggered by an intense phase of rifting, and that ongoing melting and eruption of mostly undifferentiated mafic lavas dominate the > 8 Ma long magmatic history of this region. Among the mafic lavas, most basanites and tephrites predate the alkali basalts and hawaiites, suggesting an overall temporal evolution towards less SiO2-undersaturated primary melts and increasing degrees of melting over time. The peak in alkali basaltic to hawaiitic magmatism slightly post dates the flare-up of genetically related felsic magmatism, by no more than ~ 1 Ma. This is consistent with a model in which the magmatic plumbing system erupted successively from upper to lower levels, i.e. from more evolved to more primitive compositions. One young age for a basanitic sample suggests that silica saturation decreased again towards the end of volcanic activity. This chronology of volcanic events is in good agreement with previous models, suggesting continuous lithospheric thinning beneath the SVF as a response to an extensional regime and asthenospheric uplift in the northern alpine realm.

Geochemistry and Origin of the Neoproterozoic Natkusiak Flood Basalts and Related Franklin Sills, Victoria Island, Arctic Canada

The Natkusiak continental flood basalts and Franklin sills of Victoria Island preserve an exceptional record of the ca. 716–723 Ma Franklin large igneous province and are synchronous with major climatic variations and breakup of the supercontinent Rodinia. The Natkusiak Formation basalts record an early phase of discontinuous rubbly flows (<100 m, low-Ti Type 1 magmas) overlain by a thicker series of extensive tholeiitic sheet flows (~1100 m, high-Ti Type 2 magmas). Coeval intrusions hosted by underlying Shaler Supergroup sedimentary rocks are differentiated low-Ti Type 1 Franklin sills and doleritic high-Ti Type 2 sills, both of which show correlations in isotope plots with the northernmost basalts on Victoria Island. Whole-rock Pb-Sr-Nd-Hf isotopic compositions from 66 samples indicate that the earliest magmas (Type 1) had similar primary melt compositions (Fo90 olivine) to oceanic island basalts and incorporated up to 10% granitoid basement (initial eNd = –0.8 to –7, Nb/La = 0.42 to 0.67), a relatively weak continental signature compared to many other continental flood basalt provinces. Type 2 doleritic sills and the northern sheet flow basalts incorporated up to 5% granitoid (initial eNd = +0.9 to +5.5), consistent with a waning continental influence during maturation of the magmatic system. Radiogenic isotope ratios are not correlated with indices of fractional crystallisation, which indicates that the continental material was either dispersed within the melt source, or that magmas were heterogeneously contaminated prior to differentiation. In the southwestern part of Victoria Island, Type 1 basalts show negligible continental influence (Nb/La = 0.81 to 0.94) and have unusually high initial eNd ratios (+4.4 to +11.8) that are decoupled from initial eHf (+0.8 to +11.1). These radiogenic eNd compositions persist throughout the southern volcanic stratigraphy and indicate involvement of a component with high time-integrated Sm/Nd that lacked correspondingly high Lu/Hf. We suggest that the source region for the southwestern Natkusiak basalts and related sills included isotopically matured oceanic crust, which was recycled through the asthenospheric mantle into a laterally heterogeneous plume. The distinct trace element signatures of the southern and northern sources became attenuated with the onset of voluminous melting (corresponding to emplacement of the Type 2 doleritic sills and sheet flow basalts) and may reflect contributions from hydrous eclogitic material emplaced into the lithospheric mantle during the ca. 1.9 Ga Wopmay Orogeny. As both the northern and southern volcanic rocks exhibit contrasting isotopic signatures throughout the preserved stratigraphy, the magma plumbing system must have limited lateral mixing and homogenisation, which allowed for the expression of distinct mantle source signatures in the high-level sills and basaltic lavas.

Petrology of mafic enclaves in the 2006–2012 eruptive products of Bezymianny Volcano, Kamchatka

This paper reports the results of the first comprehensive petrological study of mafic enclaves widespread in the products of recent (2006–2012) eruptions of Bezymianny Volcano, Kamchatka. Four types of mafic enclaves were distinguished on the basis of the composition and morphology of minerals, P–T conditions of formation of mineral assemblages, and structural and textural characteristics of the rocks. Disequilibrium assemblages of mafic enclaves indicate a complex structure of the magmatic plumbing system of the volcano, including a shallow chamber with andesite–basaltic andesite magmas and a deep reservoir filled in part with plagioclase–hornblende cumulates and fed by basic magmas with mantle harzburgite xenoliths. The mafic enclaves were formed at different levels of the magmatic plumbing system of the volcano and correspond to different degrees of mixing of interacting magmas. The most abundant enclaves were formed during magma ascent from the deep reservoir (960–1040°C, 5–9 kbar) into the shallow andesitic chamber (940–980°C). Enclaves of plagioclase–hornblende cumulates from the basic magmas feeding the deep reservoir (T > 1090°C and P > 9 kbar) are much less common.

Timescales of magma ascent and degassing and the role of crustal assimilation at Merapi volcano (2006–2010), Indonesia: Constraints from uranium-series and radiogenic isotopic compositions

We present new 238U-230Th-226Ra-210Pb-210Po, 87Sr/86Sr and 143Nd/144Nd isotopic data of whole-rock samples and plagioclase separates from volcanic deposits of the 2006 and 2010 eruptions at Merapi volcano, Java, Indonesia. These data are combined with available eruption monitoring, petrographic, mineralogical and Pb isotopic data to assess current theories on the cause of a recent transition from effusive dome-building (2006) to explosive (2010) activity at the volcano, as well as to further investigate the petrogenetic components involved in magma genesis and evolution. Despite the significant difference in eruption style, the 2006 and 2010 volcanic rocks show no significant difference in (238U/232Th), (230Th/232Th) and (226Ra/230Th) activity ratios, with all samples displaying U and Ra excesses. The 226Ra and 210Pb excesses observed in plagioclase separates from the 2006 and 2010 eruptions indicate that a proportion of the plagioclase grew within the decades preceding eruption. The 2006 and 2010 samples were depleted in 210Po relative to 210Pb ((210Po/210Pb)i < 1) at the time of eruption but were variably degassed (69%–100%), with the degree of 210Pb degassing strongly related to sample texture and eruption phase. In good agreement with several activity monitoring parameters, 210Po ingrowth calculations suggest that initial intrusion into the shallow magma plumbing system occurred several weeks to a few months prior to the initial 2010 eruption. The 2006 and 2010 samples show a wide range in (210Pb/226Ra) activity ratio within a single eruption at Merapi and are largely characterised by 210Pb deficits ((210Pb/226Ra) < 1). Assuming a model of complete radon degassing, the 210Pb deficits in the 2006 volcanic rocks indicate relatively longer degassing timescales of ∼2–4 years than those given by the 2010 samples of ∼0–3 years. The uranium-series and radiogenic isotopic data do not support greater crustal assimilation of carbonate material as the explanation for the more explosive behaviour of Merapi in 2010 (as has been previously suggested) and instead indicate that relatively rapid ascent of a more undegassed magma was the primary difference responsible for the transition in explosive behaviour. This interpretation is in good agreement with gas monitoring data, previous petrological studies (mineral, microlite and melt inclusion work) and maximum calculated timescale estimates using Fe-Mg compositional gradients in clinopyroxene, that also suggest more rapid movement of relatively undegassed magma in 2010 relative to 2006.

Compositional and volumetric development of a monogenetic lava flow field: The historical case of Paricutin (Michoacán, Mexico)

Paricutin volcano is the youngest and most studied monogenetic volcano in the Michoacán-Guanajuato volcanic field (Mexico), with an excellent historical record of its nine years (February 1943 to March 1952) of eruptive activity. This eruption offered a unique opportunity to observe the birth of a new volcano and document its entire eruption. Geologists surveyed all of the eruptive phases in progress, providing maps depicting the volcano's sequential growth. We have combined all of those previous results and present a new methodological approach, which utilizes state of the art GIS mapping tools to outline and identify the 23 different eruptive phases as originally defined by Luhr and Simkin (1993). Using these detailed lava flow distribution maps, the volume of each of the flows was estimated with the aid of pre- and post-eruption digital elevation models. Our procedure yielded a total lava flow volume ranging between 1.59 and 1.68km3 DRE, which is larger than previous estimates based on simpler methods. In addition, compositional data allowed us to estimate magma effusion rates and to determine variations in the relative proportions of the different magma compositions issued during the eruption. These results represent the first comprehensive documentation of the combined chemical, temporal, and volumetric evolution of the Paricutin lava field and provide key constraints for petrological interpretations of the nature of the magmatic plumbing system that fed the eruption.

Paleomagnetically inferred ages of a cluster of Holocene monogenetic eruptions in the Tacámbaro-Puruarán area (Michoacán, México): Implications for volcanic hazards

The paleomagnetic dating procedure was applied to a cluster of four partly overlapping monogenetic Holocene volcanoes and associated lava flows, namely La Tinaja, La Palma, Mesa La Muerta, and Malpaís de Cutzaróndiro, located in the Tacámbaro-Puruarán area, at the southeastern margin of the Michoacán-Guanajuato volcanic field. For this purpose, 21 sites distributed as far apart as possible from each other were sampled to obtain a well-averaged mean paleomagnetic direction for each single lava flow. For intensity determinations, double-heating Thellier experiments using the IZZI protocol were conducted on 55 selected samples. La Tinaja is the oldest of these flows and was dated by the 14C method at ~5115±130years BP (cal 4184–3655BCE). It is stratigraphically underneath the other three flows with Malpaís de Cutzaróndiro lava flow being the youngest. The paleomagnetic dating procedure was applied using the Matlab archaeo-dating tool in couple with the geomagnetic field model SHA.DIF.14k. Accordingly, for La Tinaja several possible age ranges were obtained, of which the range 3650–3480BCE is closest to the 14C age. Paleomagnetic dating on La Palma produced a unique age range of 3220–2880BCE. Two ages ranges of 2240–2070BCE and 760–630BCE were obtained for Mesa La Muerta and a well-constrained age of 420–320BCE for Malpaís de Cutzaróndiro. Although systematic archaeological excavations have so far not been carried out in this area, it is possible that the younger eruptions were contemporary to local human occupation. Our paleomagnetic dates indicate that all four eruptions, although closely clustered in space, occurred separately in time with varying recurrence intervals ranging between ~300 and ~2300years. This finding should be considered when constraining the nature of the magmatic plumbing system and developing a strategy aimed at reducing risk in the volcanically active Michoacán-Guanajuato volcanic field, where several young monogenetic volcano clusters have been identified recently. These enigmatic small “flare-ups” (outbursts of small pods of magma in geologically short periods of time within a small area) have also been encountered in other subduction-related volcanic fields around the globe (e.g. Cascades arc in the western U.S.A.) and still require to be investigated by geophysical and petrological means in order to understand their origin.

Relationship Betweenb-Value Distribution and the Magma Plumbing System in and Around Mt. Tarumae, Japan

Mt. Tarumae is an active volcano located in the southeast of the Shikotsu caldera, Hokkaido, Japan. Recently, crustal expansion occurred in 1999–2000 and 2013 near the summit of Mt. Tarumae, with a M5.6 earthquake recorded west of the summit on July 8, 2014. In this study, we determined hypocenter distributions and performedb-value analysis for the period between August 1, 2014 and August 12, 2016 to improve our understanding of the geometry of the magma system beneath the summit of Mt. Tarumae. Hypocenters were mainly distributed in two regions: 3–5 km west of Mt. Tarumae, and beneath the volcano. We then determinedb-value distributions. Regions with relatively highb-values (∼1.3) were located at depths of –0.5 to 2.0 km beneath the summit and at depths greater than 6.0 km about 1.5–3.0 km northwest of the summit, whereas a region with relatively lowb-values (∼0.6) was located at depths of 2.0–6.0 km beneath the summit. Based on comparison of theb-value distributions with other geophysical observations, it was found that the highb-value region from –0.5 to 2.0 km in depth corresponded to regions of lower resistivity, positive self-potential anomaly, and an inflation source detected in 1999–2000. Therefore, it is inferred that this region was generated by crustal heterogeneity, a decrease in effective normal stress, and change of frictional properties caused by the development of faults and fissures and the circulation of hydrothermal fluids. On the other hand, the inflation source detected in 2013 was located near the boundary between the lowb-value region beneath the summit and the deeper highb-value region about 1.5–3.0 km northwest of the summit. Studies of other volcanoes have suggested that such highb-values likely correspond to the presence of a magma chamber. Based on the deeper highb-value region estimated in this study, the magma chamber is inferred to be located at depths greater than 6.0 km about 1.5–3.0 km northwest of the summit. Thus, these findings contribute to our understanding of the magma plumbing system beneath the summit of Mt. Tarumae.

Lithospheric magma dynamics beneath the El Hierro Volcano, Canary Islands: insights from fluid inclusions

At active volcanoes, petrological studies have been proven to be a reliable approach in defining the depth conditions of magma transport and storage in both the mantle and the crust. Based on fluid inclusion and mineral geothermobarometry in mantle xenoliths, we propose a model for the magma plumbing system of the Island of El Hierro (Canary Islands). The peridotites studied here were entrained in a lava flow exposed in the El Yulan Valley. These lavas are part of the rift volcanism that occurred on El Hierro at approximately 40–30 ka. The peridotites are spinel lherzolites, harzburgites, and dunites which equilibrated in the shallow mantle at pressures between 1.5 and 2 GPa and at temperatures between 800 and 950 °C (low-temperature peridotites; LT), as well as at higher equilibration temperatures of 900 to 1100 °C (high-temperature peridotites; HT). Microthermometry and Raman analyses of fluid inclusions reveal trapping of two distinct fluid phases: early type I metasomatic CO2-N2 fluids (X N2 = 0.01–0.18; fluid density (d) = 1.19 g/cm3), coexisting with silicate-carbonate melts in LT peridotites, and late type II pure CO2 fluids in both LT (d = 1.11–1.00 and 0.75–0.65 g/cm3) and HT (d = 1.04–1.11 and 0.75–0.65 g/cm3) peridotites. While type I fluids represent metasomatic phases in the deep oceanic lithosphere (at depths of 60–65 km) before the onset of magmatic activity, type II CO2 fluids testify to two fluid trapping episodes during the ascent of xenoliths in their host mafic magmas. Identification of magma accumulation zones through interpretation of type II CO2 fluid inclusions and mineral geothermobarometry indicate the presence of a vertically stacked system of interconnected small magma reservoirs in the shallow lithospheric mantle between a depth of 22 and 36 km (or 0.67 to 1 GPa). This magma accumulation region fed a short-lived magma storage region located in the lower oceanic crust at a depth of 10–12 km (or 0.26–0.34 GPa). Following our model, the 40–30-ka-old volcanic activity of El Hierro is related to this mantle-based magma system, a system that we propose fed the recent 2011–2012 eruption.

The shallow depth emplacement of mafic intrusions on a magma-poor rifted margin: An example from the Bight Basin, southern Australia

Three-dimensional (3D) seismic data has provided important insights into how magma is transported and stored in sedimentary basins. Much of our understanding is based on volcanic-rifted margins. In contrast, the magma plumbing systems of magma-poor rifted margins are less well studied. This study uses 3D seismic data to describe the magma plumbing system of the Bight Basin Igneous Complex; a volcanic province located along the magma-poor southern Australian margin. Here, magma emplacement occurred in the Mid Eocene, some 40 million years after the onset of seafloor spreading. Our data images a variety of 5–270 m thick, 2–23 km diameter mafic intrusions that are confined to within 1200 m of the paleoseabed. Intrusions emplaced at ≤200 m depth formed hybrid and compound sills and are typified by elongate, lava flow-like morphologies. Intrusions emplaced at depths of 200–1200 m formed saucer-shaped and compound sills, as well as laccoliths. Approximately 60% of the intrusions have lava flows and volcanogenic vents above their shallowest tips, as shown from previous studies. This suggests that the intrusions played a crucial role in transporting magma to the paleoseabed. Furthermore, many of the intrusions are overlain by forced folds, highlighting the role of magma emplacement in inducing overburden deformation. The sills and laccoliths rarely form interlinked complexes and instead most commonly occur as isolated bodies. This suggests that they were fed by dykes. We infer that high rates of magma ascent in the dykes prevented their transition into sills within sediments at >1.2 km depth. Our study highlights that the magma plumbing system of the Bight Basin Igneous Complex contains a diversity of magmatic intrusions, the morphology of which is linked to their emplacement depth, host sediment rheology and the physical properties of the magma. This plumbing system contrasts markedly with those found along better-studied volcanic rifted margins.

Insights into the petrogenesis of low- and high-Ti basalts: Stratigraphy and geochemistry of four lava sequences from the central Paraná basin

Lava flow sequences were sampled in the central part of the Paraná basin aiming to verify the time-related evolution of the Paraná basaltic magmatism. It is shown that low- and high-Ti basalts were erupted synchronously. In particular, Esmeralda and Pitanga flows are interlayered, with the former prevailing in the upper part of the sequence. Evidence for synchronously active magma plumbing systems is also supported by mineralogical data, showing signs of mixing between the two groups.Geochemical data, including Sr-Nd-Pb isotopic compositions are furthermore used to define the mantle source of various low- (Esmeralda and Gramado) and high-Ti (Pitanga and Urubici) magma types. Involvement of a carbonatitic component is proposed for the genesis of the basalts (particularly for the Urubici ones) as suggested by trace element enrichments unrelated to significant isotopic variations. This carbonatitic signature of the mantle source may be conveyed by CO2-rich metasomatic fluids or melts percolating upwards within the sub-continental lithospheric mantle (SCLM) leading to rapid and selective enrichment of incompatible trace elements. Metasomatism was probably localized at the outskirts of the basin, were Urubici tholeiites and contemporaneous carbonatites were erupted. Geochemical data also suggest the occurrence of significant amounts of crustal contamination in the LTi magmas (mainly in the Gramado and in the late Esmeralda lavas) while crustal assimilation seems negligible in the HTi samples. Globally, a very complex picture arises for the genesis of the Paraná tholeiites, with near-synchronous and geographically coincident flows undergoing significantly different extents of interaction with the crust and tapping different mantle sources.

Magmatic plumbing systems of the Koryakskii–Avacha Volcanic Cluster as inferred from observations of local seismicity and from the regime of adjacent thermal springs

An analysis of local seismicity within the Avacha–Koryakskii Volcanic Cluster during the 2000–2016 period revealed a sequence of plane-oriented earthquake clusters that we interpret as a process of dike and sill emplacement. The highest magmatic activity occurred in timing with the 2008–2009 steam–gas eruption of Koryakskii Volcano, with magma injection moving afterwards into the cone of Avacha Volcano (2010–2016). The geometry of the magma bodies reflects the NF geomechanical conditions (tension and normal faults, $$S_V>S_{H_{\text{max}}}>S_{h_{\text{min}}}$$ ) at the basement of Koryakskii Volcano dominated by vertical stresses S v , with the maximum horizontal stress $$S_{h_{\text{max}}}$$ pointing north. A CFRAC simulation of magma injection into a fissure under conditions that are typical of those in the basement of Koryakskii Volcano (the angle of dip is 60°, the size is 2 × 2 km2, and the depth is –4 km abs.) showed that when the magma discharge is maintained at the level of 20000 kg/s during 24 hours the fissure separation increases to reach 0.3 m and the magma injection is accompanied by shear movements that occur at a rate as high as 2 × 10–3 m/s, thus corresponding to the conditions of local seismic events with Mw below 4.5. We are thus able to conclude that the use of planeoriented clusters of earthquakes for identification of magma emplacement events is a physically sound procedure. The August 2, 2011 seismicity increase in the area of the Izotovskii hot spring (7 km from the summit of Koryakskii Volcano), which is interpreted as the emplacement of a dike, has been confirmed by an increase in the spring temperature by 10–12°C during the period from October 2011 to July 2012.

Magmatic architecture within a rift segment: Articulate axial magma storage at Erta Ale volcano, Ethiopia

Understanding the magmatic systems beneath rift volcanoes provides insights into the deeper processes associated with rift architecture and development. At the slow spreading Erta Ale segment (Afar, Ethiopia) transition from continental rifting to seafloor spreading is ongoing on land. A lava lake has been documented since the twentieth century at the summit of the Erta Ale volcano and acts as an indicator of the pressure of its magma reservoir. However, the structure of the plumbing system of the volcano feeding such persistent active lava lake and the mechanisms controlling the architecture of magma storage remain unclear. Here, we combine high-resolution satellite optical imagery and radar interferometry (InSAR) to infer the shape, location and orientation of the conduits feeding the 2017 Erta Ale eruption. We show that the lava lake was rooted in a vertical dike-shaped reservoir that had been inflating prior to the eruption. The magma was subsequently transferred into a shallower feeder dike. We also find a shallow, horizontal magma lens elongated along axis inflating beneath the volcano during the later period of the eruption. Edifice stress modeling suggests the hydraulically connected system of horizontal and vertical thin magmatic bodies able to open and close are arranged spatially according to stresses induced by loading and unloading due to topographic changes. Our combined approach may provide new constraints on the organization of magma plumbing systems beneath volcanoes in continental and marine settings.

伊豆大島火山のマグマ供給システムと噴火予測（特集「無人観測ロボット その2」）

伊豆大島火山のマグマ供給システムと噴火予測（特集「無人観測ロボット その2」）

Geochronological and geochemical constraints on sulfide mineralization in the Qingmingshan mafic intrusion in the western part of the Proterozoic Jiangnan orogenic belt along the southern margin of the Yangtze Craton

The Qingmingshan magmatic Fe-Ni-Cu sulfide deposit is located in the southern margin of the Yangtze Craton, in the western part of the Proterozoic Jiangnan orogenic belt. Sulfide mineralization occurs as disseminated zones in a dike-like mafic body intruding Early Neoproterozoic meta-sedimentary rocks. The host intrusion can be further divided into an upper gabbronorite unit, a middle dolerite unit and a lower gabbronorite unit. Each unit contains at least one disseminated sulfide zone. Sulfide mineral assemblages are composed of pyrrhotite, pentlandite and chalcopyrite. Zircon U–Pb dating yields a crystallization age of 847.8±3.8Ma (2σ) for the intrusion. The host rocks are characterized by enrichment in light rare earth elements (REE) relative to heavy REE and show pronounced negative Nb-Ta anomalies, similar to the characteristics of arc mafic igneous rocks worldwide. Mixing calculations indicate that the variations of (Th/Yb)N and (Th/Nb)N in the intrusion are consistent with between 10 and 20wt% crustal contamination of the parental magma. The sulfide separates are characterized by elevated δ34S values from 1.5 to 3.5‰, with a peak at 2.6‰, and by γOs values from 2 to 89, indicating addition of crustal sulfur and Os to the parental magma. The trace element data and S-Os isotopes support the premise that sulfide saturation in the Qingmingshan magma was triggered by crustal contamination. Variations of PGE tenors in different sulfide zones can be explained by variable R-factors from 300 to 1000 during sulfide segregation from the magma. Mass balance calculations indicate that the parental magma was severely depleted in PGE, which could be a primary signature of arc basaltic magma formed by low degrees of partial melting or due to previous sulfide segregation at depth. The lithological and mineralization zonation across the dike-like Qingmingshan mafic intrusion indicates that at least three separate pulses of sulfide-charged magma were involved in the formation of the Qingmingshan magmatic sulfide deposit. The Qingmingshan mafic intrusion is interpreted to be a part of a sub-arc magma plumbing system. Coeval ultramafic intrusions formed from more primitive magma in the region are good exploration targets for magmatic sulfide deposits with higher Ni/Cu ratios than the Qingmingshan deposit.

Volcano Geodesy: Recent developments and future challenges

Ascent of magma through Earth's crust is normally associated with, among other effects, ground deformation and gravity changes. Geodesy is thus a valuable tool for monitoring and hazards assessment during volcanic unrest, and it provides valuable data for exploring the geometry and volume of magma plumbing systems. Recent decades have seen an explosion in the quality and quantity of volcano geodetic data. New datasets (some made possible by regional and global scientific initiatives), as well as new analysis methods and modeling practices, have resulted in important changes to our understanding of the geodetic characteristics of active volcanism and magmatic processes, from the scale of individual eruptive vents to global compilations of volcano deformation. Here, we describe some of the recent developments in volcano geodesy, both in terms of data and interpretive tools, and discuss the role of international initiatives in meeting future challenges for the field.