Active Lava Lake Research Articles

Hot rings on Io observed by Juno/JIRAM

It has been thought that Io’s many paterae may contain lava lakes, but observations by NASA’s Galileo spacecraft at sufficiently high resolution were limited to a few locations, such as Loki Patera. Data acquired by NASA’s Juno mission in May 2023 reveal a common set of thermal characteristics for at least ten paterae on Io, with bright (hot) “thermal rings” around the perimeter of their floors. Loki, Surt, Fuchi, Amaterasu, Mulungu, Chors, and Dazhbog Patera, and four unnamed paterae, all show the same pattern of surface temperature. This new Juno/JIRAM data suggests that hot rings around paterae are a common phenomenon, and that they are indicative of active lava lakes. All the investigated paterae lack recent lava flows on their flanks, suggesting that at the time of observations, the level of the lake was not high enough to overflow the rim. These observations provide insight about the characteristics of paterae’s activity, which may involve either central upwelling of magma, or up-and-down “piston-like” vertical motion of the lake surface. Tidal forces, which are extreme at Io, could play a role as well. Future observations from Juno, particularly during the closest flybys, may indicate which mechanism is more plausible.

Magma recharge in persistently active basaltic–andesite systems and its geohazards implications: the case of Villarrica volcano, Chile

We report whole-rock chemistry, mineral chemistry, and volatile content from Villarrica volcano’s major recent paroxysms and background activity. Composition of the volcanic products are basalt to basaltic andesite with whole-rock SiO2 content between 50 and 56 wt%, and a mineralogy dominated by olivine (Fo71-80), clinopyroxene (Mg# ~ 50) and plagioclase (An60–80). Volatile contents in melt inclusions are up to 1.5 wt% H2O, 500 ppm CO2, 1230 ppm sulfur and 580 ppm chlorine. Regardless of the type of activity, there are no substantial variations in whole-rock composition or the volatile content when the activity switches from background activity to a major paroxysm, strongly suggesting that this shift does not just depend on the arrival of new magma in the shallow magmatic system. Geothermobarometry constrains crystallization of the major mineral phases at various depths between 3 and 12.7 km, suggesting that degassing of a volatile-rich recharge magma occurs deeper than 12 km, producing efficient mixing throughout the whole system, and sustaining the lava lake activity in Villarrica’s summit crater. The occurrence of a permanent lava lake also suggests that the magma recharge must be close to continuous and therefore sudden changes between background and paroxysmal volcanic activity are likely controlled by relatively small changes in the rate of recharge and/or the volatile release rate in the magmatic system. This has important implications for the understanding of eruption triggers and the forecasting of volcanic eruptions.Graphical abstract

Plumbing System Architecture and Differentiation Processes of the Nyiragongo Volcano, DR Congo

ABSTRACT The Nyiragongo volcano is one of the most alkali-rich volcanic centres on the planet (Na2O + K2O generally >10 wt.%, agpaitic index up to 1.34), characterized by a semi-permanently active lava lake which hosts silica-undersaturated (SiO2 < 40 wt.%), low viscosity lavas. To improve our understanding of this unique magmatic system, we present a set of 291 samples, acquired during new field excursions between 2017 and 2021. The major and trace element composition of all samples was measured, revealing a lithological range extending from primitive picrites (Mg# 82) erupted from parasitic cones to a variety of highly evolved nephelinites, leucitites, and melilitites erupted from the main edifice as recently as 2002, 2016, and 2021. We measured major and trace element compositions from the full spectrum of minerals present in all sampled lithologies. From these we calculated that the main magma reservoirs feeding Nyiragongo are at approximately 9–15 and 21–33 km depth, in agreement with recent seismic observations. Fractional crystallization modelling using observed mineral compositions and proportions was performed to quantitatively link the lithologies to specific residual liquid fractions assuming evolution from an olivine-melilite parental melt. Our modelling indicates that fractionation and cumulate formation in deep chambers reduces the remaining melt fraction to ~60%, after which melts are injected into upper, liquid dominated magma chambers where fractionation and accumulation of clinopyroxene, melilite, and feldspathoids dominates. Characterisation of mineral textures and geochemistry reveals high crystal mobility in a repeatedly recharging plumbing system split between liquid-dominated, evolved magma chambers and more solid-dominated, primitive mushes, decreasing in liquid fraction with depth.

Life-cycles of a lava lake: Ambrym volcano, Vanuatu

Eighteen years of heat flux data collected by the Moderate Resolution Imaging Spectroradiometer (MODIS) are used to investigate cycles of lava lake activity at Ambrym volcano in Vanuatu. Lava lakes are dynamic manifestations of open-vent volcanic behaviour and time-series analysis of this behaviour reveals the timing and magnitude of magmatic processes. Pulsatory fluctuations in heat flux of > 50 MW over periods of 30 to >100 days during lava lake activity are interpreted to represent episodes of magma replenishment. These are periodically interrupted by rapid falls in heat flux (50–100 MW over <5 days) or variations in the relative heat flux from lava lakes in different vents of the volcano, which are interpreted to represent changes to the shallow plumbing system feeding lava lake activity. Sudden lava lake drainage and flank eruptions can interrupt or terminate lava lake cycles. Heat flux patterns preceding such events do not show an obvious consistent pattern, suggesting that different magmatic processes may initiate these potentially hazardous surges in eruption intensity. Comparison with the dynamics of other lava lakes reveals similarities in long-term behaviour, despite differences in lava lake size, surface motion or degassing. Collectively this suggests fundamental similarities in the processes that sustain open-vent behaviour, including, the stability of magma supply, conduit-plumbing/geometry and local strain-state.

Investigating Phases of Thermal Unrest at Ambrym (Vanuatu) Volcano through the Normalized Hot Spot Indices Tool and the Integration with the MIROVA System

Ambrym is an active volcanic island, located in the Vanuatu archipelago, consisting of a 12 km-wide summit caldera. This open vent volcano is characterized by an almost persistent degassing activity which occurs in the Benbow and Marum craters, which were also the site of recent lava lakes. On 15 December 2018, about three years after an intense lava effusion, the first recorded since 1989, a small-scale intra-caldera fissure eruption occurred. On 16 December, the eruption stopped, and the lava lakes at the Benbow and Marum craters were drained. In this work, we investigated the thermal activity of the Ambrym volcano, before, during, and after the 15 December 2018 eruption, using daytime Sentinel-2 (S2) Multispectral Instruments (MSI) and Landsat-8 (L8) Operational Land Imager (OLI) data, at a mid-high spatial resolution. The results were integrated with Moderate Resolution Imaging Spectroradiometer (MODIS) observations. Outputs of the Normalized Hotspot Indices (NHI) tool, retrieved from S2-MSI and L8-OLI data, show that the thermal activity at the Ambrym craters increased about three weeks before the 15 December 2018 lava effusion. This information is consistent with the estimates of volcanic radiative power (VRP), which were performed by the Middle Infrared Observation of Volcanic Activity (MIROVA) system, by analyzing the nighttime MODIS data. The latter revealed a significant increase of VRP, with values above 700 MW at the end of the October–November 2018 period. Moreover, the drastic reduction of thermal emissions at the craters, marked by the NHI tool since the day of the fissure eruption, is consistent with the drop in the lava lake level that was independently suggested in a previous study. These results demonstrate that the S2-MSI and L8-OLI time series, combined with infrared MODIS observations, may contribute to detecting increasing trends in lava lake activity, which may precede effusive eruptions at the open vent volcanoes. This study addresses some challenging scenarios regarding the definition of possible threshold levels (e.g., in terms of VRP and total Short Wave Infrared radiance) from the NHI and MIROVA datasets, which could require special attention from local authorities in terms of the occurrence of possible future eruptions.

Spatio-temporal evolution of the magma plumbing system at Masaya Caldera, Nicaragua

Volcanic unrest in calderas can be exhibited through a variety of different mechanisms, such as changes in seismicity and ground deformation, as well as variations in thermal and/or gas emissions. However, not all caldera unrest results in explosive caldera-forming volcanic activity. Alternative activity may include periods of quiescence, passive degassing, effusive activity (e.g., lava flows lava lakes and dome formation), and/or magma injection into the shallow magma system. In this study, we perform a long-term study (spanning 2011–2019) of ground deformation at Masaya using six Interferometric Synthetic Aperture Radar (InSAR) datasets. Masaya exhibited bi-modal eruptive behavior between 2011 and 2019, dominated by open-vent lava lake activity and punctuated by short-lived summit explosions. The Multidimensional Small BAseline Subset time-series analysis approach was used to take advantage of the temporally dense SAR datasets. Between 2012 and early 2015, we observed degassing-induced pressurization of the Masaya Central Reservoir (MCR) at an estimated volume change rate of ~ 0.28 $$\times$$ 106 m3/year. In May 2015, magma was supplied into the MCR at a rate of ~ 5.6 $$\times$$ 106 m3/year, leading to the appearance of a summit lava lake in December 2015. Over the next 6 months, rapid magma supply continued to drive lava lake activity and was followed by a cessation of magma supply into the MCR for another 11 months. From mid-2017 to end-2019, we observed depressurization (~ − 0.67 $$\times$$ 106 m3/year) of the MCR due to a lack of magma supply and continued high rates of degassing in-conjunction with declining lava lake activity.

Top–down control on eruptive style at Masaya volcano inferred from melt composition

Highly explosive Plinian eruptions of basaltic magma are enigmatic because low melt viscosities should inhibit such eruptive style. Masaya volcano, Nicaragua, is a persistently active basaltic system capable of a wide range of eruptive styles, from open-conduit lava lake activity to voluminous Plinian eruptions; it is thus an ideal natural laboratory to constrain potential controls on basaltic eruption style. Here we report the major, trace, and volatile (CO2, H2O, S, Cl, F) element composition of olivine-, plagioclase- and clinopyroxene-hosted melt inclusions as well as matrix glasses from lava lake ejecta and two Plinian tephra deposits—the 2.1 ka Masaya Triple Layer and the 1.8 ka Ticuantepe Lapilli—to test whether pre-eruptive volatile contents and degassing history may be linked to eruptive style. All samples display a relatively narrow and largely overlapping basaltic–basaltic andesitic compositional range (51.7±1.0 wt.% SiO2, 4.8±0.4 wt.% MgO) with similar trace element signatures (e.g., Ce/Y=0.82±0.10, Ba/La=74±11). However, lava lake and Plinian samples show systematic differences in pre-eruptive volatile contents, forming distinct groups with mean H2O contents of 0.6±0.2 wt.% (lava lake), 1.1±0.2 wt.% (Masaya Triple Layer), and 1.9±0.3 wt.% (Ticuantepe Lapilli). Together, these groups generate broad positive correlations between S, Cl and H2O concentrations, with maximum values reaching 920 ppm, 1300 ppm and 2.3 wt.%, respectively, which are low compared to typical Central American arc magmas. Magma temperature estimates overlap and average at 1115±30°C, while volatile saturation pressures are low, mainly <100 MPa, although only lava lake samples record pressures <31 MPa. These observations reiterate the compositionally buffered state of the volcano's magmatic system highlighted by previous work and demonstrate that — regardless of eruption style — all Masaya magmas undergo variable, but extensive, pre-eruptive degassing at low pressure. Geohygrometry, gas emissions, and H2O/Ce–Ba/La systematics suggest initial, undegassed H2O contents on the order of 3.9–5.5 wt.%. Our results imply that pre-eruptive volatile contents are not the culprit for Plinian events at Masaya. Instead, we propose that the volcano's vigorous magma supply is modulated in a top–down manner to produce a wide range of eruptive styles, whereby temporary sealing of the conduit may instigate a transition to explosive behavior. In this model, rapid magma ascent is triggered when the seal eventually breaks from degassing-induced pressurization, yielding high degrees of undercooling, rapid microlite growth, and a dramatic increase in magma viscosity and explosive eruption potential. There may thus be a thin line between open-conduit conditions and Plinian eruptions at Masaya.

SO&amp;lt;sub&amp;gt;2&amp;lt;/sub&amp;gt; and BrO emissions of Masaya volcano from 2014 to 2020

Abstract. Masaya (Nicaragua, 12.0∘ N, 86.2∘ W; 635 m a.s.l.) is one of the few volcanoes hosting a lava lake, today. This study has two foci: (1) discussing the state of the art of long-term SO2 emission flux monitoring with the example of Masaya and (2) the provision and discussion of a continuous data set on volcanic gas data with a large temporal coverage, which is a major extension of the empirical database for studies in volcanology as well as atmospheric bromine chemistry. We present time series of SO2 emission fluxes and BrO/SO2 molar ratios in the gas plume of Masaya from March 2014 to March 2020 – covering the three time periods (1) before the lava lake appearance, (2) a period of high lava lake activity (November 2015 to May 2018), and (3) after the period of high lava lake activity. For these three time periods, we report average SO2 emission fluxes of (1000±200), (1000±300), and (700±200) t d−1 and average BrO/SO2 molar ratios of (2.9±1.5)×10-5, (4.8±1.9)×10-5, and (5.5±2.6)×10-5. Our SO2 emission flux retrieval is based on a comprehensive investigation of various aspects of spectroscopic retrievals, the wind conditions, and the plume height. We observed a correlation between the SO2 emission fluxes and the wind speed in the raw data. We present a partial correction of this artefact by applying dynamic estimates for the plume height as a function of the wind speed. Our retrieved SO2 emission fluxes are on average a factor of 1.4 larger than former estimates based on the same data. Further, we observed different patterns in the BrO/SO2 time series: (1) an annual cyclicity with amplitudes between 1.4 and 2.5×10-5 and a weak semi-annual modulation, (2) a step increase by 0.7×10-5 in late 2015, (3) a linear trend of 1.4×10-5 per year from November 2015 to March 2018, and (4) a linear trend of -0.8×10-5 per year from June 2018 to March 2020. The step increase in 2015 coincided with the lava lake appearance and was thus most likely caused by a change in the magmatic system. We suggest that the cyclicity might be a manifestation of meteorological cycles. We found an anti-correlation between the BrO/SO2 molar ratios and the atmospheric water concentration (correlation coefficient of −0.47) but, in contrast to that, neither a correlation with the ozone mixing ratio (+0.21) nor systematic dependencies between the BrO/SO2 molar ratios and the atmospheric plume age for an age range of 2–20 min after the release from the volcanic edifice. The two latter observations indicate an early stop of the autocatalytic transformation of bromide Br− solved in aerosol particles to atmospheric BrO.

The petrologic and degassing behavior of sulfur and other magmatic volatiles from the 2018 eruption of Kīlauea, Hawaiʻi: melt concentrations, magma storage depths, and magma recycling

Kīlauea Volcano’s 2018 lower East Rift Zone (LERZ) eruption produced exceptionally high lava effusion rates and record-setting SO2 emissions. The eruption involved a diverse range of magmas, including primitive basalts sourced from Kīlauea’s summit reservoirs. We analyzed LERZ matrix glasses, melt inclusions, and host minerals to identify melt volatile contents and magma storage depths. The LERZ glasses and melt inclusions span nearly the entire compositional range previously recognized at Kīlauea. Melt inclusions in Fo86-89 olivine from the main eruptive vent (fissure 8) underwent 70–170 °C cooling during transport in LERZ carrier melts, causing extensive post-entrapment crystallization and sulfide precipitation. Many of these melt inclusions have low sulfur (400–900 ppm) even after correction for sulfide formation. CO2 and H2O vapor saturation pressures indicate shallow melt inclusion trapping depths (1–5 km), consistent with formation within Kīlauea’s Halemaʻumaʻu and South Caldera reservoirs. Many of these inclusions also have degassed δ34S values (− 1.5 to − 0.5‰). Collectively, these results indicate that some primitive melts experienced near-surface degassing before being trapped into melt inclusions. We propose that decades-to-centuries of repeated lava lake activity and lava drain-back during eruptions (e.g., 1959 Kīlauea Iki) recycled substantial volumes of degassed magma into Kīlauea’s shallow reservoir system. Degassing and magma recycling from the 2008–2018 Halemaʻumaʻu lava lake likely reduced the volatile contents of LERZ fissure 8 magmas, resulting in lower fountain heights compared to many prior Kīlauea eruptions. The eruption’s extreme SO2 emissions were due to high lava effusion rates rather than particularly volatile-rich melts.

How to turn off a lava lake? A petrological investigation of the 2018 intra-caldera and submarine eruptions of Ambrym volcano

In December 2018, an unusually large intra- and extra-caldera eruption took place at Ambrym volcano (Vanuatu). The eruption drained the volcano’s five active lava lakes and terminated, at least momentarily, the surface activity that had been ongoing for decades to hundreds of years, sustaining the largest recorded persistent degassing on the planet. Here, we investigate the mechanisms and dynamics of this major eruption. We use major elements and volatiles in olivine and clinopyroxene-hosted melt inclusions, embayments, crystals and matrix glasses together with clinopyroxene geobarometry as well as olivine and clinopyroxene geothermometry and diffusion modelling in crystals and embayments to reconstruct the chronology and timing of the subsurface processes that accompanied the eruption. We find that the eruption began with the meeting, mingling and limited chemical mixing of mostly two magma bodies occupying similar vertical but different horizontal locations in the crust, one corresponding to the main plumbing system at Ambrym that fed the lava lakes and the other corresponding to an older, previously cutoff and more chemically evolved branch of the plumbing system. Within the primitive magma, two texturally distinct components—one microlite rich and one microlite poor—can further be identified. The 2018 eruption hence provides a detailed image of Ambrym’s complex plumbing system. Our diffusion timescales and geobarometric estimates coincide closely with geophysical observations. They point to a reconnection of the evolved magmatic branch with the main system occurring less than 10 h prior to the intra-caldera eruption and a period of 2 days for the subsequent > 30-km lateral magma transport along a deeper dike prior to submarine eruption just off the SE coast of the island with the more primitive magma reaching first followed by mingled magma containing both compositions. Magma ascent rates are estimated at 95 ± 24 m/s in the last ~ 2.5 km of ascent during the intra-caldera eruption and at 80 ± 6 m/s in the last ~ 4 km of ascent during the submarine eruption. Comparison with other lava lake draining eruptions reveals striking similarities both in terms of precursory activity, with lake level rising prior to the eruption in all cases, and in terms of plumbing system organization with the presence of peripheral magma pockets, isolated from the main magmatic system but that can be mobilized and erupted when met by dikes propagating laterally from the main system.

Nyiragongo and Nyamuragira: a review of volcanic activity in the Kivu rift, western branch of the East African Rift System

Nyiragongo and Nyamuragira are two active volcanoes of the Western branch of the East African Rift in the Virunga area. They were built at the Kivu rift axis ca. 12,000 years ago and set above two tectonic steps separated by the Kameronze Fault. Both volcanoes have displayed a succession of intra-crater and flank eruptions that have been observed and documented since the end of the nineteenth century. Here, we have collated and reviewed these publications and reports. Nyiragongo is famous for its semi-permanently active lava lake, which at the time of writing (2020) was the largest in the world. During the construction of the main stratovolcano, which ended a few centuries ago with a caldera collapse, the lava composition changed from melilitite to leucite and then to melilite-bearing nephelinite. The historically active lava lake is believed to be directly fed from an upper intra-volcano reservoir, a shallow reservoir situated a few kilometres below the volcano in the granite basement, and a deeper intra-crustal magma chamber. Historic activity has been documented since 1894 and can be divided into eight stages, on the basis of sudden changes between lava filling and draining, with cycles of rising lava lake activity and overflows, followed by sinking and complete or partial drainage. Twice in recent history (in 1977 and 2002), major flank eruptions were accompanied by complete drainage of the lava lake and the upper plumbing system. The lava that filled the crater since 1948, and then again after the 1977 and 2002 drainage events have been calculated at a cumulative volume of around 324 × 106 m3. In comparison, the 1977 and 2002 flank eruptions involved 47 × 106 m3 of lava. The average annual output rate associated with crater filling is thus estimated at between 4 and 13 × 106 m3. Lava lake behaviour changes from equilibrium, with alternation between gas pistoning and spattering regimes through disequilibrium with intermittent activity, to complete disappearance of the lava lake. These changes can be related to the conditions of the descent of dense degassed magma from the upper conduit into the shallow reservoir. However, since 1959, the chemical composition of the leucite and melilite-bearing nephelinite lavas has not significantly changed, which implies a magma supply from the same magma batch. Nyamuragira was characterised by a shield building phase of activity until a caldera collapse, ca. 300 to 500 years ago. The post-caldera phase has involved effusive activity in the caldera and at numerous flank fissures. The plumbing system consists of an upper reservoir roughly at the basement-volcano interface and averaging a volume of 400 × 106 m3, a shallow upper-crust stratified reservoir, and a middle-crust mafic magma chamber. Volcanic activity has involved a succession of filling and emptying events at the upper reservoir. Lava volumes of historic eruptions reveal annual output rates averaging 14 × 106 m3 between 1901 and 1976, and 40 × 106 m3 between 1976 and 2012. A drastic increase in activity occurred in December 1976. This event coincided with the January 1977 flank eruption of Nyiragongo and resulted from a main tectonic event in the rift basement that improved the efficiency of magma ascent at both volcanoes. The historic lava composition can be related to six cycles of magma accumulation in, and withdrawal of, the upper reservoir from the shallow stratified reservoir. Similar magma storage and transport systems are known in many effusive systems: the Nyiragongo lava lake shares behavioural characteristics similar to those observed at Kilauea, Erebus, and Erta’Ale. At Nyiragongo and Nyamuragira, magma supply and persistent activity with sudden changes of the magma output rates in relation to tectonic events are also comparable with those of Kilauea, Piton de la Fournaise of Reunion island, and Mount Etna.

Recent Activity of Nyiragongo (Democratic Republic of Congo): New Insights From Field Observations and Numerical Modeling

AbstractNyiragongo volcano is known for its active lava lake and for socioeconomic issues arising from future possible eruptive events having major impacts on the community living in the Virunga region. The 2020 field expedition inside the summit crater has allowed the collection of unprecedented field observations to state on the current activity. Since the February 2016 intracrater event, the crater floor level has been rising much faster than during the 2010–2016 period. The current activity is reminiscent of the 1970–1972 and 1994–1995 periods preceding the lava lake drainage events in 1977 and 2002. Numerical simulations, successfully validated with data over the past 30 years, show that the rising of the crater floor could slow down in the next months/years and reach a critical equilibrium. Based on the past eruptive history and on the current activity, a flank eruption in the March 2024 to November 2027 interval could be a possible scenario.

Magma reservoir failure and the onset of caldera collapse at Kīlauea Volcano in 2018.

Caldera-forming eruptions are among Earth's most hazardous natural phenomena, yet the architecture of subcaldera magma reservoirs and the conditions that trigger collapse are poorly understood. Observations from the formation of a 0.8-cubic kilometer basaltic caldera at Kīlauea Volcano in 2018 included the draining of an active lava lake, which provided a window into pressure decrease in the reservoir. We show that failure began after <4% of magma was withdrawn from a shallow reservoir beneath the volcano's summit, reducing its internal pressure by ~17 megapascals. Several cubic kilometers of magma were stored in the reservoir, and only a fraction was withdrawn before the end of the eruption. Thus, caldera formation may begin after withdrawal of only small amounts of magma and may end before source reservoirs are completely evacuated.

Periodic outgassing as a result of unsteady convection in Ray lava lake, Mount Erebus, Antarctica

Persistently active lava lakes show continuous outgassing and open convection over years to decades. Ray Lake, the lava lake at Mount Erebus, Ross Island, Antarctica, maintains long-term, near steady-state behavior in temperature, heat flux, gas flux, lake level, and composition. This activity is superposed by periodic small pulses of gas and hot magma every 5-18 min and disrupted by sporadic Strombolian eruptions. The periodic pulses have been attributed to a variety of potential processes including unstable bidirectional flow in the conduit feeding the lake. In contrast to hypotheses invoking a conduit source for the observed periodicity, we test the hypothesis that the behavior could be the result of dynamics within the lake itself, independent of periodic influx from the conduit. We perform numerical simulations of convection in Ray Lake driven by both constant and periodic inflow of gas-rich magma from the conduit to identify whether the two cases have different observational signatures at the surface. Our simulations show dripping diapirs or pulsing plumes leading to observable surface behavior with periodicities in the range of 5-20 min. We conclude that a convective speed faster than the inflow speed can result in periodic behavior without requiring periodicity in conduit dynamics. This finding suggests that the surface behavior of lava lakes might be less indicative of volcanic conduit processes in persistently outgassing volcanoes than previously thought, and that dynamics within the lava lake itself may modify or overprint patterns emerging from the conduit. Beyond its importance for Erebus, this result highlights the need to understand periodic processes emerging at all levels of the volcanic plumbing system in order to interpret surface observations.

Seismicity and outgassing dynamics of Nyiragongo volcano

Active lava lakes at volcanoes can be regarded as open windows to their magmatic systems. The dynamics of such lakes may therefore provide decisive insights into deeper magmatic processes, potentially leading to fundamental theoretical implications and volcano monitoring improvements. Among the rare volcanoes on Earth hosting a persistent lava lake, Nyiragongo in D.R. Congo directly threatens a massive population of roughly 1 million inhabitants. Here we analyze close-range (i.e., summit) and distant (around 17 km) seismic measurements acquired at this African volcano between 2011 and 2018 in order to better understand the seismic signature of the lava lake activity and how it relates to the deeper volcanic processes. Both summit and distant seismic records contain a similar continuous tremor pattern attributable to the lava lake activity. Combining this information with time-lapse camera images and lava lake level measurements confirms the mechanism of gas pistoning at Nyiragongo, which is characterized by short-duration (a few minutes long) and meter-scale level variations during the period of observation. We also characterize the dominant periodicity of this shallow tremor signature of about a few tens of minutes. Because this marked periodic pattern varies during a significant one-month fluctuation of SO2 emissions (estimated from space), we suggest that this particular seismic periodicity corresponds to the convective lake movement driven by the persistent degassing typical of active open-vent volcanoes. Finally, new seismic evidence reveals the effect of deep magmatic intrusion and consequent major pressure changes in the plumbing system, resulting in sudden and large drops of the lava lake level associated with strong degassing. Such transient episodes have similar characteristics to total lava lake drainage associated with flank eruptions already observed at this volcano in 1977 and 2002, or at Kīlauea volcano in 2018.

A global synthesis of lava lake dynamics

Active lava lakes represent a variety of open-vent volcanism in which a sizable body of lava accumulates at the top of the magma column, constrained by the vent and/or crater geometry. The longevity of lava lakes reflects a balancing of cooling and outgassing occurring at the surface and input of hot and gas-rich magma from below. Due to their longevity and relative accessibility, lava lakes provide a natural laboratory for studying fundamental volcanic processes such as degassing, convection and cooling. This article examines all seven lakes that existed at the time of writing in 2018, located in the Pacific, Antarctica, Africa, and South and Central America. These lakes span all tectonic environments, and a range of magma compositions. We focus on analysis of the lake surface motion using image velocimetry, which reveals both similarities and contrasts in outgassing and lake dynamics when comparing the different lakes. We identify two categories of lake behavior: Organized (Erta’Ale, Nyiragongo, Kīlauea after 2011, and Erebus) and Chaotic (Villarrica, Masaya, Marum). This division does not map directly to lake size, viscosity, gas emission rate, or temperature. Instead, when examined together, we find that the lakes follow a linear relationship between average surface speed and the ratio of total gas flux to lake surface area. This relationship points to the combined importance of both flux and lake size in addition to the total volume of gas emission, and suggests that a shared deep mechanism controls the supply of heat and gas to all lakes. On the other hand, the differences between Chaotic and Organized lakes highlight the important role of the geometry of the conduit-lake transition, which superimposes a shallow signal on that of the deep circulation. The spatial patterns of surface motion we document suggest that the release of gas bubbles at Chaotic lakes is more efficient (i.e., bubbles are less likely to be retained and recycled) compared with Organized lakes. In addition, the data presented here indicate that the solidified crust of Organized lakes plays a role in regulating convection and outgassing in lava lakes.

A Rapidly Convecting Lava Lake at Masaya Volcano, Nicaragua

Lava lakes provide a rare opportunity to study conduit flow processes through direct observation of the exposed magma surface. The recent lava lake activity at Masaya volcano (Nicaragua), which commenced in 2015, displayed several unusual phenomena. We report on the dynamics of this rapidly convecting lake, which, to the best of our knowledge manifested the highest lava flow velocities ever reported for a lava lake: 13.7–16.4 m s−1, in addition to unusual fluid dynamic behavior involving alteration in surface flow direction. We studied this system with multiparametric and high time resolution remote sensing measurements, performed during June 2017, including ultraviolet camera observations of SO2 emission rates, near infrared thermal camera measurements and video analyses of the lake surface. Median SO2 emission rates of 3.1 (±0.8) and 3.7 (±0.9) kg s−1 were found, which are lower than previously published estimates, and could represent challenging remote sensing conditions or a waning in lava lake activity. Video analyses enabled characterization of frequent bursts of approximately hemispherical spherical-cap bubbles on the surface with diameters ranging 0.6–8.5 m (median of 2.6 m), and calculation of individual bubble masses, which contribute to active bubble bursting values estimated at 1.9 to 3.9 kg s−1. We show that only a small fraction, 7–17%, of total emission volumes are contributed by these bubbles, based on estimated emission rates of 22.5 and 26.9 kg s−1. Furthermore, periodicity analysis reveals regular 200–300 s oscillations in SO2 emissions. These are not shared by any of our other datasets and particularly during the contemporaenously acquried thermal data, hence, we tentatively assign an atmospheric causal generation mechanism, driven by atmospheric transport and turbulence phenomena, such as eddying. Overall, we highlight the uniquely high velocity and fluid dynamic behavior of Masaya lava lake.

A review of controls on lava lake level: insights from Halema‘uma‘u Crater, Kīlauea Volcano

The height of the lava column is a fundamental measure of open-vent volcanic activity, but little continuous long-term data exist to understand this parameter. The recent (2008–2018) lava lake activity at the summit of Kīlauea Volcano provides an opportunity to track and understand the processes that control lava level over timescales ranging from minutes to years. We review recently published data as well as analyze the long-term record of lava level at Kīlauea’s summit during 2009–2017. Longer timescale fluctuations, over days to months, have a strong positive linear correlation with ground deformation, suggesting they reflect pressure changes in the summit magma reservoir. Shorter timescale fluctuations, over minutes to hours, have an inverse relationship with spattering (i.e., outgassing) intensity at the lake surface—increased outgassing lowers the lake level and vice versa. Our analysis of the long-term lava level data thus confirms that lava level at Halema‘uma‘u is controlled by these two types of processes: (1) those related to magma reservoir pressure changes (such as magma supply rate) and (2) shallow outgassing fluctuations (such as gas pistoning). Frequency filtering can isolate pressure- and outgassing-driven components of lava level change. Time series analysis indicates that there was no large, persistent periodicity in the lava level; a minor fortnightly oscillation might be related to Earth tides but was not continuously present. The correlation between lava level and deformation of the summit indicates that the lava lake acts as a reliable “piezometer”; tracking lava level over time may thus provide an indication of flank eruptive potential. We show that long-term interdisciplinary monitoring is necessary to discriminate the processes that control lava level.

Slug Stability in Flaring Geometries and Ramifications for Lava Lake Degassing

AbstractLong‐lived, persistently active lava lakes are rare but have been studied extensively because they may provide a natural laboratory for studying magmatic convection in basaltic and other low‐viscosity volcanoes. However, lava lakes differ from other volcanic systems in various ways, particularly geometry: the presence of the lava lake requires a sudden and significant flaring of the conduit. The goal of this paper is to advance our understanding of the effect that a sudden increase in the conduit width has on the stability of large gas slugs and bubbles ascending from the plumbing system. We investigate this question by linking analog laboratory experiments to direct numerical simulations in two dimensions. We find that the rapid change in the geometry associated with the transition from the conduit to the lava lake causes large gas slugs to break up. The nondimensional regime over which we observe the breakup depends sensitively on the flare angle of the lava lake and the importance of inertial effects in the flow, implying that degassing‐related eruptive activity at open‐system volcanoes is sensitive to the upper conduit and lake geometries. Applying this idea to the likely geometries and viscosities of actual lava lakes, our study suggests that all lava lakes except Erebus fall well within the regime where large gas slugs are prone to breakup.

Radar Altimetry as a Robust Tool for Monitoring the Active Lava Lake at Erebus Volcano, Antarctica

AbstractThe level of lava within a volcanic conduit reflects the overpressure within a connected magma reservoir. Continuous monitoring of lava level can therefore provide critical insights into volcanic processes and aid hazard assessment. However, accurate measurements of lava level are not easy to make, partly owing to the often dense fumes that hinder optical techniques. Here we present the first radar instrument designed for the purpose of monitoring lava level and report on its successful operation at Erebus volcano, Antarctica. We describe the hardware and data‐processing steps followed to extract a time series of lava lake level, demonstrating that we can readily resolve ∼1 m cyclic variations in lake level that have previously been recognized at Erebus volcano. The performance of the radar (continuous, automated data collection in temperatures of around −30 °C) indicates the suitability of this approach for sustained automated measurements at Erebus and other volcanoes with lava lakes.