Open Conduit System Research Articles

Crystal mush interaction controls eruptive style during the 2018 Kīlauea fissure eruption

We use new geochemical, petrological, and rheological data to constrain the formation and emplacement of the highly compositionally unusual(andesitic basalt) Kīlauea 2018 Fissure 17 (F17) eruptive products. Despite the restricted spatial and temporal distribution, F17 samples are texturally and geochemically diverse. The western samples are enriched in SiO2 by up to 10 wt%, relative to their eastern equivalents; additionally, the western samples contain microcrystalline enclaves, absent from the homogenous eastern samples. The compositions erupted along F17 suggest interaction between the basaltic 2018 juvenile magma and a crystal mush at depth, likely a left-over from the nearby 1955 eruption. Magma mingling caused heating and local melting of remnant mush, leading to melt hybridization and volatile exsolution. Rapid water exsolution likely caused overpressurization of the reservoir underneath the western side of F17, leading to Strombolian explosions of viscous magma, in contrast to sustained Hawaiian fountaining on the eastern side. Remelting of remnant crystal mush and melt hybridization in open-conduit systems may hence be an effective mechanism in inducing volatile saturation.

A long time of rest at Llaima volcano following the 2010 MW 8.8 Maule earthquake, Chile

Llaima Volcano (38.692°S, 71.729°W) is one of the most active centers in South America, with ∼56 eruptions since 1852 CE. During the last eruptive cycle (2007–2009), the volcano displayed typical features of open conduit systems, with frequent long-period volcanic seismicity (LP). After the MW 8.8 Maule megathrust earthquake (27 February 2010), Llaima volcano located 305 km to SE of the epicenter, has experienced one of the longest periods of quiescence since 1852, with reduced activity observed since the beginning of its continuous instrumental monitoring in the year 2007. To better understand this behavior, we track the repose time between eruptions in the 20th century to depict the current period of quiescence and provide details of the pattern of volcano seismicity after the earthquake. The number of LP events decreased ∼90% shortly after the earthquake. In turn, crustal faults nearby (∼20 km south of the main crater) showed an increased level of seismicity lasting for approximately six months after the mainshock before returning to background levels. In order to propose driving mechanisms for such behavior, we computed the static changes of the stress tensor on chosen receiver structures. The acting receiver structure differs from those related to the long-term stress regime controlled by regional tectonics and was inferred from the pre-earthquake shallow seismicity (2007–2009). LP and volcanotectonic seismicity (VT) before the earthquake were located in a NW-SE pattern, in contrast to the NE-trending alignment of flank vents (and hence dike swarms) usually interpreted as an indicator of the maximum horizontal stress axis. From the static stress transfer imparted by the earthquake, we find that an opening event would have occurred at a NW-trending receiver structure coupled with closing of other potential structures, which promoted magma storage along this blind structure but precluding further propagation to the surface. Our results thus show a rare example of earthquake-induced suppression of volcanic activity, where the geometry of the fault-fracture network that shape the plumbing system likely plays a major role.

Pre-eruptive Conditions of the 3 March 2015 Lava Fountain of Villarrica Volcano (Southern Andes)

Villarrica or Rukapillan (35.9°S; 2,847 m a.s.l.) is one of the most active volcanoes in South America and is the highest-risk volcano in Chile. It has an open conduit with a persistent lava lake. On the 3 March 2015, Strombolian activity rapidly progressed into a 1.5-km-high lava fountain, erupting at least ∼ 2.4 × 106 m3 of tephra. Soon after, the activity returned to mild Strombolian “background” explosions, which lasted until early 2017. Understanding the pre-eruptive conditions of such paroxysmal events is fundamental for volcanic hazard assessment. We present major and trace element geochemistry for glass and crystalline phases of basaltic andesite paroxysm pyroclasts (52–56 wt.% SiO2), and for the subsequent Strombolian “background” activity through February 2017 (54–56 wt.% SiO2). The lava fountain source magma was initially stored in a deeper and hotter region (9.4–16.3 km; ca. 1140 °C) and was then resident in a shallow (≤ 0.8 km) storage zone pre-eruption. During storage, crystallising phases comprised plagioclase (An66–86), olivine (Fo75–78) and augite (En46–47). Equilibrium crystallisation occurred during upper-crustal magmatic ascent. During storage in the shallower region, magma reached H2O saturation, promoting volatile exsolution and over-pressurization, which triggered the eruption. In contrast, subsequent “background” explosions involving basaltic-andesite were sourced from a depth of ≤ 5.3 km (ca. 1110 °C). Pre-eruptive conditions for the 2015 lava fountain contrast with historical twentieth-century eruptions at Villarrica, which were likely driven by magma that underwent a longer period of mixing to feed both effusive and explosive activity. The rapid transition to lava-fountaining activity in 2015 represents a challenging condition in terms of volcano monitoring and eruption forecasting. However, our petrological study of the pyroclastic materials that erupted in 2015 offers significant insights into eruptive processes involving this type of eruption. This aids in deciphering the mechanisms behind sudden eruptions at open conduit systems.

Volcanic Activity of Sakurajima Monitored Using Global Navigation Satellite System

A dense Global Navigation Satellite System (GNSS) network has been deployed at Sakurajima volcano since 1995 and extends to the surrounding area of the Aira caldera. The ground deformation obtained by GNSS observation corresponds to transient eruptive activity of Sakurajima volcano, which has produced frequent vulcanian eruptions since 1955. Inflation of the volcano was detected prior to the increase in vulcanian eruptions in 1999, and resumption of the eruptions at the Showa crater. Magma intrusion events and an increase in eruptions in late 2009, late 2011, and early 2015 suggest the existence of an open-conduit system from the Aira caldera to the vents at the summit area of the Minamidake cone, through the sub reservoir beneath the older Kitadake cone. Ground deformation induced by sudden dike intrusion is different from that of previous intrusions, as revealed by the dense GNSS network. GNSS data are useful in evaluating and forecasting volcanic activity, and are available to grasp the advection and diffusion of volcanic ash.

Rapid mixing and short storage timescale in the magma dynamics of a steady-state volcano

Steady-state volcanic activity implies equilibrium between the rate of magma replenishment and eruption of compositionally homogeneous magmas, lasting for tens to thousands of years in an open conduit system. The Present-day activity of Stromboli volcano (Aeolian Islands, Southern Italy) has long been recognised as typical of a steady-state volcano, with a shallow magmatic reservoir (highly porphyritic or hp-magma) continuously refilled by more mafic magma (with low phenocryst content or lp-magma) at a constant rate and accompanied by mixing, crystallisation and eruption. Our aim is to clarify the timescale and dynamics of the plumbing system at the establishment of the Present-day steady-state activity (<1.2 ka) to pinpoint the onset of the steady-state regime. We investigated the Post-Pizzo (PP) pyroclastic sequence (∼1.7–1.5 ka) and one of the Early Paroxysms (EP) of the Present-day activity, focusing on the clinopyroxene population.Whole rock and clinopyroxene compositional variation among the PP and EP magmas is consistent with the time progression of the Stromboli system towards more mafic and lower 87Sr/86Sr compositions, pointing to the chemical and isotopic signature of the Present-day activity. Clinopyroxenes from both PP and EP record a complex history with compositional zoning that reflects growth in three different melt domains: a high-Mg# proto-lp recharging magma, a low-Mg# proto-hp resident magma, and a transient intermediate-Mg# magma. These are the result of complex turbulent flow fields and mixing regimes produced by repeated injections of the proto-lp magma in the shallow proto-hp magma reservoir. During the PP period the magmatic system was already able to regain the pre-input proto-hp composition, gradually changing toward a less evolved signature after the injection(s) of the more mafic proto-lp magma, owing to efficient (days to a few years) stirring and melt homogenisation (i.e., homogenisation time < residence time).Based upon Fe–Mg diffusion in clinopyroxene the total residence time during PP and EP periods, from the arrival of the mafic magma in the shallow system until the eruption, ranges from 1 to ∼50 years. Longer residence times (up to 150 years) have been recorded in the initial phase of the PP sequence, possibly testifying to the transition from a closed- to the open-conduit, steady-state regime of the Present-day activity.Some clinopyroxenes from the PP recorded the mafic triggering event of the feeding proto-lp magma occurring within few months to a few days before eruption. Remarkably, other clinopyroxene portions crystallised and captured the rapid timescales (a few days) of the on-going mixing and homogenisation process between the proto-lp and the proto-hp magmas leading to the eruption.The modelling of clinopyroxene zoning events at Stromboli provides evidence for growth and storage in three different melt domains, and sets robust constraints on their residence time from lp-magma recharge(s) to eruption, along with the timescales of melt homogenisation and triggering events. The lifetime history captured by Fe–Mg zoning of Stromboli clinopyroxenes suggests that the interplay between rapid mixing and short storage timescales can be a key parameter controlling the dynamics of the plumbing system of steady-state volcanoes.

Ultrafast syn-eruptive degassing and ascent trigger high-energy basic eruptions

Lithium gradients in plagioclase are capable of recording extremely short-lived processes associated with gas loss from magmas prior to extrusion at the surface. We present SIMS profiles of the 7Li/30Si ion ratio in plagioclase crystals from products of the paroxysmal sequence that occurred in the period 2011–2013 at Mt. Etna (Italy) in an attempt to constrain the final ascent and degassing processes leading to these powerful eruptions involving basic magma. The observed Li concentrations reflect cycles of Li addition to the melt through gas flushing, and a syn-eruptive stage of magma degassing driven by decompression that finally produce significant Li depletion from the melt. Modeling the decreases in Li concentration in plagioclase by diffusion allowed determination of magma ascent timescales that are on the order of minutes or less. Knowledge of the storage depth beneath the volcano has led to the quantification of a mean magma ascent velocity of ~43 m/s for paroxysmal eruptions at Etna. The importance of these results relies on the application of methods, recently used exclusively for closed-system volcanoes producing violent eruptions, to open-conduit systems that have generally quiet eruptive periods of activity sometimes interrupted by sudden re-awakening and the production of anomalously energetic eruptions.

Volatiles and energy released by Puracé volcano

Total CO2 output of Purace volcano (Colombia) was estimated on the basis of fluids discharged by fumaroles, soil gases, and dissolved carbon species in the aquifer. The soil CO2 emission was computed from a field survey of 512 points of CO2 soil flux measurements at the main degassing areas of Purace volcano. The CO2 flux from Purace’s plume was estimated using an indirect method, that used the SO2 plume flux and CO2/SO2 ratio of the main high temperature fumarole. The total output of CO2 was estimated at ≅ 1500 t/day. The main contribution of CO2 comes from the plume (summit degassing) and from soil degassing that emit 673 and 812 t/day, respectively. The contributions of summit and soil degassing areas are comparable, indicating an intermediate degassing style partitioned between closed and open conduit systems. The estimated water vapor discharge (as derived from the chemical composition of the fumaroles, the H2O/CO2 ratio, and the SO2 plume flux) allowed calculation of the total thermal energy (fumarolic, soil degassing, and aquifer) released from the Purace volcanic system. This was 360 MW.

Tracking dynamics of magma migration in open-conduit systems

Open-conduit volcanic systems are typically characterized by unsealed volcanic conduits feeding permanent or quasi-permanent volcanic activity. This persistent activity limits our ability to read changes in the monitored parameters, making the assessment of possible eruptive crises more difficult. We show how an integrated approach to monitoring can solve this problem, opening a new way to data interpretation. The increasing rate of explosive transients, tremor amplitude, thermal emissions of ejected tephra, and rise of the very-long-period (VLP) seismic source towards the surface are interpreted as indicating an upward migration of the magma column in response to an increased magma input rate. During the 2014 flank eruption of Stromboli, this magma input preceded the effusive eruption by several months. When the new lateral effusive vent opened on the Sciara del Fuoco slope, the effusion was accompanied by a large ground deflation, a deepening of the VLP seismic source, and the cessation of summit explosive activity. Such observations suggest the drainage of a superficial magma reservoir confined between the crater terrace and the effusive vent. We show how this model successfully reproduces the measured rate of effusion, the observed rate of ground deflation, and the deepening of the VLP seismic source. This study also demonstrates the ability of the geophysical network to detect superficial magma recharge within an open-conduit system and to track magma drainage during the effusive crisis, with a great impact on hazard assessment.

Deep fluid transfer evidenced by surface deformation during the 2014–2015 unrest at Piton de la Fournaise volcano

Identifying the onset of volcano unrest and providing an unequivocal identification of volcano reawakening remain challenging problems in volcanology. At Piton de la Fournaise, renewal of eruptive activity in 2014–2015, after 41months of quiescence and deflation, was associated with long-term continuous edifice inflation measured by GNSS. Inflation started on June 9, 2014, and its rate progressively increased through 2015. Inflation onset was rapidly followed by an eruption on June 20–21, 2014, showing that volcano reactivation can be extremely fast, even after long non-eruptive phases. This short-lived eruption involved a shallow source (1.3–1.9km depth below the summit). The inflation that followed, and eruptions in 2015, involved a larger depth range of fluid accumulation, constrained by inverse modeling at ca. 3.9 to 1.2–1.7km depth. This time evolution reveals that volcano reawakening was associated with continuous pressurization of the shallowest parts of its plumbing system, triggered by progressive upwards transfer of magma from greater depth. A deep magma pulse occurred in mid-April 2015 and was associated with deep seismicity (3 to 9.5km depth) and CO2 enrichment in fluids emitted by summit fumaroles. From this date, ground deformation accelerated and the output rates of eruptions increased, culminating in the long-lasting, large-volume, August–October eruption (~36Mm3). This evolution suggests that deep magma/fluid transfer through an open conduit system first provoked the expulsion of the top of the plumbing system in June 2014, and then induced the progressive vertical transfer of the entire plumbing system down to 9km (four eruptions in 2015). The new sustained feeding of the volcano was also at the origin of the hydrothermal system perturbation and the acceleration of the eastern flank motion, which favor lateral dike propagation and the occurrence of frequent and increasingly large eruptions. Our results highlight the fast and progressive way in which basaltic magmatic systems can wake up.

Shallow magma‐mingling‐driven Strombolian eruptions at Mt. Yasur volcano, Vanuatu

Mt. Yasur volcano (Vanuatu) has been increasingly recognized for its high‐frequency Strombolian eruptions. Strombolian activity is often regarded as a product of the rapid ascent of gas slugs originating from a deep magma, which mingle with a batch of shallow magma upon eruption. Heterogeneous crystal distribution as well as bimodal bubble‐size distributions, in the eruptive products, generally supports this view. Here, the Strombolian activity at Mt. Yasur is analyzed. A rheological investigation indicates that the basaltic‐andesitic eruptive products contain an apparently homogeneous glass phase and yet, exhibit evidence of a distinct range of glass transition temperatures with multiple peaks occurring in individual samples. Such anomalous behavior is proposed to result from the mingling of magmas with contrasting oxidation states. The unstable nature of the measured glass transition behavior leads us to the inference that mingling is located in the shallow parts of the eruptive conduits, partly driven by rejuvenation of material slumped from the crater walls into an open conduit system. The dynamics of this process may expose the periodicity of the eruptions themselves.

Regimes of magma recharge and their control on the eruptive behaviour during the period 2001–2005 at Mt. Etna volcano

Mount Etna volcano (Italy) during the period 2001–2005 has undergone a period of intense eruptive activity marked by three large eruptions (2001, 2002–2003 and 2004–2005). These eruptions encompassed diverse eruptive styles and regimes: from intensely explosive, during 2001 and 2002–2003 eruptions, to exclusively effusive in the 2004–2005 event. In this work, we put forward the idea that these three eruptions are the response of the progressive arrival into the uppermost segment of the open-conduit system of a new magma, which was geochemically distinct in terms of trace element and Sr–Nd–Pb isotope signature from the products previously emitted by the Etnean volcano. The magma migrated upwards mainly through a peripheral tectonic system, which can be considered as eccentric in spite of its relative proximity to the main system. The ingress of the new magma and its gradual displacement from the eccentric system into the uppermost sector of the open-conduit gave rise to different eruptive behaviours. At the beginning, the ascent of the undegassed magma, able to exsolve a gas phase at depth, and its interaction with closed-system magma reservoirs less than 10 km deep gave rise to the explosive events of 2001 and 2002–2003. Later, when the same magma entered into the open-conduit system, it took part in the steady-state degassing and partially lost its volatile load, leading to a totally effusive eruption during the 2004–2005 event. One further consideration highlighted here is that in 2001–2005, migration of the feeding axis from an eccentric and peripheral position towards the main open-conduit has led to the development of a new vent (South East Crater 2) located at the eastern base of the South East Crater through which most of the subsequent Etnean activity occurred.

Platinum group element and nickel sulphide ore tenors of the Mount Keith nickel deposit, Yilgarn Craton, Australia

A set of platinum group element (PGE) analyses of about 120 samples from a 250-m continuous drill core through the Mount Keith komatiite-hosted nickel orebody, combined with Ni, Cu, Co, S, and major elements, reveals a complex trend of covariance between the original cumulus components of a thick sequence of nearly pure olivine–sulphide liquid adcumulates. The intersection is divided into informal chemostratigraphic zones, defined primarily by combinations of fine-scale cyclicity in original olivine composition, defined by Mg#, and sulphide composition, defined by Pt/S and Ni/S. Contents of Ni and PGE in 100% sulphides (tenors) were determined from linear regressions of the Ni–S and PGE–S covariance for each zone. Inferred olivine compositions range from about Fo92 to Fo94.6 and show a broad decrease from bottom to top of the sequence complicated by numerous reversals, revealing crystallisation in an open conduit system. Ni and PGE tenors of Mount Keith sulphide ores have typical values similar to the type I deposits of the Kambalda Dome. Mobility of S, at least on the scale of 2-m sample composites, is evidently relatively minor. Tenors for the various zones range 12–22% Ni, 370–1540 ppb Pt, 970–3670 ppb Pd, 100–460 ppb Ir, 170–460 ppb Rh, and 710–1260 ppb Ru. Pt, Pd, and Rh tenors are very strongly correlated, but the iridium group of platinum group elements (IPGEs; Ir and Ru) less so. Tenor variations are predominantly controlled by variations in magma/sulphide ratio R (100–350), with a minor component of variance from equilibrium crystallisation trends in the parent magma. PGE depletion in the silicate melt due to sulphide liquid extraction is limited by entrainment of sulphide liquid droplets and continuous equilibration with the transporting silicate magma. Ratios of the PGEs to one another are similar to those in the host komatiite magma, with the exception of Pt, which is systematically depleted in ores, relative to Rh and Pd and relative to host magma, by a consistent factor of about 2 to 2.5. This anomalous Pt depletion relative to PGE element ratios in unmineralized komatiitic rocks matches that observed in bulk compositions of many komatiite-hosted orebodies. The highly consistent nature of this depletion, and particularly the very strong correlation between Pt, Pd, and Rh in the Mount Keith deposit, argue that this depletion is a primary magmatic signal and not an artefact of alteration. Differential diffusion rates between Pt and the other PGEs, giving rise to a low effective partition coefficient for Pt into sulphide liquid, is advanced as a possible but not definitive explanation.

Magma vesiculation and infrasonic activity at Stromboli open conduit volcano

Explosive activity at Stromboli is explained in terms of dynamics of large gas bubbles that ascend in the magma conduit and burst at the free surface generating acoustic pressure that propagates as infrasonic signals in the atmosphere. The rate and the amplitude of the infrasonic activity is directly linked to the rate and the overpressure of the bursting gas bubbles and thus reflects the rate at which magma column degasses under non-equilibrium pressure conditions. We investigate the link between explosive degassing and magma vesiculation by comparing the rate of infrasonic activity with the bubble size distributions (BSDs) of scoria clasts collected during several days of explosive activity at Stromboli. BSDs of scoria show a characteristic power law distribution, which reflect a gas bubble concentration mainly controlled by a combined process of bubble nucleation and coalescence. The cumulative distribution of the infrasonic pressure follows two power laws, indicating a clear separation between the frequent, but weak, bursting of small gas bubbles (puffing) and the more energetic explosions of large gas slugs. The exponents of power laws derived for puffing and explosive infrasonic activity show strongly correlated (0.96) changes with time indicating that when the puffing rate is high, the number of energetic explosions is also elevated. This correlation suggests that both puffing and explosive activity are driven by the same magma degassing dynamics. In addition, changes of both infrasonic power law exponents are very well correlated (0.92 with puffing and 0.87 with explosions) with variations of the BSD exponents of the scoria clasts, providing evidence of the strong interplay between scoria vesiculation and magma explosivity. Our analysis indicates that variable magma vesiculation regimes recorded in the scoria correlate with the event number and energy of the explosive activity. We propose that monitoring infrasound on active volcanoes may be an alternative way to look at the vesiculation process in open conduit systems.

The contemporaneous emission of low-K and high-K trachybasalts and the role of the NE Rift during the 2002 eruptive event, Mt. Etna, Italy

Mount Etna volcano erupted almost simultaneously on its northeastern and southern flanks between October 27 and November 3, 2002. The eruption on the northeastern flank lasted for 8 days, while on the southern flank it continued for 3 months. The northeastern flank eruption was characterized by the opening of a long eruptive fracture system between 2,900 and 1,900 m.a.s.l. A detailed survey indicates that the fractures’ direction shifted during the opening from N10W (at the NE Crater, 2,900 m) to N45E (at its lowest portion, 1,900 m) and that distinct magma groups were erupted at distinct fracture segments. Based on their petrological features, three distinct groups of rocks have been identified. The first group, high-potassium porphyritic (HKP), is made up of porphyritic lavas with a Porphyritic Index (P.I.) of 20–32 and K2O content higher than 2 wt%. The second group is represented by lavas and tephra with low modal phenocryst abundance (P.I. < 20) named here oligo-phyric (low-phyric), and K2O content higher than 2 wt% (HKO, high-potassium oligophyric). The third group, low-potassium oligophyric (LKO), consists of tephra with oligophyric texture (P.I. < 20) but K2O content < 2 wt%. K-rich magmas (HKP and HKO) are similar to the magma erupted on the southern flank, and geochemical variations within these groups can be accounted for by a variable degree of fractionation from a single parent magma. The K-poor magma (LKO), erupted only in the upper segment of the fracture, cannot be placed on the same liquid line of descent of the HK groups, and it is similar to the magmas that fed the activity of Etna volcano prior to the eruption of 1971. This is the first time since then that a magma of this composition has been documented at Mt. Etna, thus providing a strong indication for the existence of distinct batches of magma whose rise and differentiation are independent from the main conduit system. The evolution of this eruption provides evidence that the NE Rift plays a very active role in the activity of Mt. Etna volcano, and that its extensional tectonics allows the intrusion and residence of magma bodies at various depths, which can therefore differentiate independently from the main open conduit system.

Infrasonic monitoring at Stromboli volcano during the 2003 effusive eruption: Insights on the explosive and degassing process of an open conduit system

Stromboli is an open conduit volcano characterized by persistent conduit degassing and explosive dynamics. During the 2003 effusive eruption we monitored this activity with a five‐element array deployed at ∼450 m from the active craters. The array allows us to track in real time changes of explosive activity in terms of source position and excess pressure release. Infrasonic monitoring located the new active vents in the same preeruptive position, revealing that the collapse of the crater terrace has not changed the geometry of the shallow feeding system. Elevation of the infrasonic source is stable at ∼750 m above sea level, coincident with the elevation of the crater terrace. This suggests a shallow position of the magma column or, if the source is embedded in the conduit, diffraction at the crater rim. At the end of the eruption, back azimuth analysis indicates that explosive activity was progressively reestablished at all of the craters following a NE toward SW direction along the collapsed trench. The array detected significant small‐amplitude (∼1 Pa) infrasonic activity occurring at a rate of approximately one pulse every 2 s and generally located in only one vent at a time. This persistent activity is not related to Strombolian explosion, but we infer it is reflecting the overpressurized degassing of the magma column. Amplitude distribution of ∼7 million infrasonic signals shows that this degassing is 1 order of magnitude more energetic than the explosions. Besides, infrasonic amplitude distribution presents a decay with a double slope, indicating that explosions and degassing are driven by conduit dynamics acting at two different rates.

Magma mixing during the 2001 event at Mount Etna (Italy): Effects on the eruptive dynamics

During the 2001 eruptive episode three different magmas were erupted on the southern flank of Mount Etna volcano from distinct vent systems. Major and minor element chemistry of rocks and minerals shows that mixing occurred, and that the mixed magma was erupted during the last eruptive phases. The space–time integrated analysis of the eruption, supported by geophysical data, together with major and trace element bulk chemistry (XRF, ICP-MS) and major and trace mineral chemistry (EPMA, LAM ICP-MS), support the following model: 1) trachybasaltic magma rises through a NNW–SSE trending structure, connected to the main open conduit system; 2) ascent of an amphibole-bearing trachybasaltic magma from a 6 km deep eccentric reservoir through newly open N–S trending fractures; 3) just a few days following the eruption onset the two tectonic systems intersect at the Laghetto area; 4) at the Laghetto vent a mixed magma is erupted. Mixing occurred between the amphibole-bearing trachybasaltic magma and an inferred deep more basic end-member. The most relevant aspect in the eruptive dynamics is that the eruption of the mixed magma at the Laghetto vent was highly explosive due to volatile content in the magma. The gas phase formed, mainly because of the decreased volatile solubility due to rapid fractures opening and increased T, related to mixing, and partially because of the amphibole breakdown.

Research and Teaching in Forensic Science

Identifying the onset of volcano unrest and providing an unequivocal identification of volcano reawakening remain challenging problems in volcanology. At Piton de la Fournaise, renewal of eruptive activity in 2014–2015, after 41 months of quiescence and deflation, was associated with long-term continuous edifice inflation measured by GNSS. Inflation started on June 9, 2014, and its rate progressively increased through 2015. Inflation onset was rapidly followed by an eruption on June 20–21, 2014, showing that volcano reactivation can be extremely fast, even after long non-eruptive phases. This short-lived eruption involved a shallow source (1.3–1.9 km depth below the summit). The inflation that followed, and eruptions in 2015, involved a larger depth range of fluid accumulation, constrained by inverse modeling at ca. 3.9 to 1.2–1.7 km depth. This time evolution reveals that volcano reawakening was associated with continuous pressurization of the shallowest parts of its plumbing system, triggered by progressive upwards transfer of magma from greater depth. A deep magma pulse occurred in mid-April 2015 and was associated with deep seismicity (3 to 9.5 km depth) and CO2 enrichment in fluids emitted by summit fumaroles. From this date, ground deformation accelerated and the output rates of eruptions increased, culminating in the long-lasting, large-volume, August–October eruption (~ 36 Mm3). This evolution suggests that deep magma/fluid transfer through an open conduit system first provoked the expulsion of the top of the plumbing system in June 2014, and then induced the progressive vertical transfer of the entire plumbing system down to 9 km (four eruptions in 2015). The new sustained feeding of the volcano was also at the origin of the hydrothermal system perturbation and the acceleration of the eastern flank motion, which favor lateral dike propagation and the occurrence of frequent and increasingly large eruptions. Our results highlight the fast and progressive way in which basaltic magmatic systems can wake up.