Search

Variable eruptive style and explosivity is common in basaltic to basaltic andesite volcanoes but can have uncertain origins. Veniaminof volcano in the Alaska-Aleutian arc is a frequently active open-vent center, regularly producing Strombolian eruptions and small lava flows from an intracaldera cone within an intracaldera ice cap. The September–December 2018 eruption of Veniaminof evolved in explosivity over time. The eruption was documented with frequent satellite observations, syn- and post-eruption structure-from-motion photo surveys, and post-eruption sampling of lava flows and tephra preserved in the syn-eruption snowpack. Lava flows with a total volume of ~ 6 × 106 m3 flowed down the cone flanks into the ice cap, overthickening at ice marginal flow fronts. Smaller tephra deposits were estimated at ~ 1–2 × 106 m3 dense rock equivalent, with almost half of this volume deposited directly on the eruptive cone itself. Erupted products were basaltic andesite, and composition (54 ± 0.7 wt% bulk rock, 58 ± 0.7 wt% glass SiO2), sideromelane microlite crystallinity (20–30%), and microlite number density (plagioclase 6.4 ± 2.6 × 105 n/mm3) did not change significantly over the eruption suggesting a similar magma source and ascent rate. We defined tephra componentry with groundmass microcrystalline textures using backscatter electron images. The componentry of tephra groundmass showed significant increases in tachylite grains, defined here by the presence of dendritic interstitial nanolites, corresponding to increasing seismic tremor and periods of increased ash emissions. We suggest that these componentry changes reflect increasing undercooled zones on the conduit margins that increased brittle shearing, fragmentation, and ultimately ash emissions.

The 2018 lower East Rift Zone (LERZ) eruption of Kīlauea, Hawai’i, provides an excellent natural laboratory with which to test models of lava flow propagation. During early stages of eruption crises, the most useful lava flow propagation equations utilize readily determined parameters and require fewer a priori assumptions about future behavior of the flow. Here, we leverage the numerous observations of lava flows collected over the duration of the eruption crisis at Kīlauea in 2018 to test simple lava flow propagation models. These models track the one-dimensional propagation of the flows according to three main rheological restraining forces: bulk viscosity, yield strength, and growth of a surface crust. We calculate the predicted changes in length through time of three flows that vary in bulk composition, crystal content, and total flow length. Cooler flows that are more crystal-rich tend to be more dominated by crust growth, though early stages of propagation can be controlled by bulk viscosity. We find that variations in effusion rate significantly impact flows that are short-lived; flows that are produced during steady-state effusion are readily approximated by average values for the entire flow. Thus, accurate knowledge of variations in effusion rate are critical to accurate lava flow propagation forecasting.

The 2018 eruption of Kīlauea Volcano produced large and destructive lava flows from the fissure 8 (Ahu ‘aila ‘au) vent with flow velocities up to 17 m s−1, highly variable effusion rates over both short (minutes) and long (hours) time scales, and a proximal channel or spillway that displayed flow features similar to open channel flow in river systems. Monitoring such dynamic vent and lava flow systems is a challenge. Our results demonstrate that infrasound, combined with ground-based observations and imagery from unoccupied aircraft systems (UAS), can be used to distinguish vent degassing activity from high-speed lava flow activity. We use spectral characteristics and the infrasound frequency index (FI) to distinguish spillway infrasound from vent infrasound. Comparing FI with flow speeds derived from UAS videos reveals that spillway infrasound only occurs when flow speeds were sufficiently high to cause a supercritical flow state and breaking waves (Froude values > 1.7), and we propose that the spillway signals are produced primarily through the interaction of the turbulent lava-free surface with the atmosphere. We show that FI can also provide a means to track bulk effusion rate. Our results indicate that infrasound offers a new way to characterize lava flow channel hydraulics and is a powerful tool for monitoring effusive eruptions when high-speed flows are possible.

Volcano monitoring, eruption response, and hazard assessment at volcanoes in the United States of America (US) fall under the mandate of five regional volcano observatories covering 161 active volcanoes. Working in a wide range of volcanic and geographic settings, US observatories must learn from and apply new knowledge and techniques to a great variety of scientific and hazard communication problems in volcanology. Over the past decade, experience during volcanic crises, such as the landmark 2018 eruption of Kīlauea, Hawai'i, has combined with investments and advances in research and technology, and the changing needs of partner agencies and the public, to transform the operations, science, and communication programs of US volcano observatories. Scientific and operational lessons from the past decade now guide new research and growing inter-observatory and external communication networks to meet new challenges and improve detection, forecasting, and response to volcanic eruptions in the US and around the world.

Harrat Rahat, one of several large, basalt-dominated volcanic fields in western Saudi Arabia, is a prime example of continental, intraplate volcanism. Excellent exposure makes this an outstanding site to investigate changing volcanic flux and composition through time. We present 93 40Ar/39Ar ages and six 36Cl surface-exposure ages for volcanic deposits throughout northern Harrat Rahat that, when integrated with a new geologic map, define 12 eruptive stages. Exposed volcanic deposits in the study area erupted 570 ka. Over the past 570 k.y., the average eruption rate was 0.14 km3/k.y., but volcanism was episodic with periods alternating between low (0.04–0.06 km3/k.y.) and high (0.1–0.3 km3/k.y.) effusion rates. Before 180 ka, eruptions vented from the volcanic field’s dominant eastern vent axis and from a subsidiary, diffuse, western vent axis. After 180 ka, volcanism focused along the eastern vent axis, and the composition of volcanism varied systematically along its length from basalt dominated in the north to trachyte dominated in the south. We hypothesize that these compositional variations

The 2018 eruption on the lower East Rift Zone of Kīlauea Volcano produced one of the largest and most destructive lava flows in Hawai’i during the past 200 years. Over the course of more than 3 months, twenty-four fissures erupted, and the rate of lava effusion varied by two orders of magnitude, with significant implications for evolving flow behavior and hazards. Syn-eruptive data were collected to quantify these changes in lava effusion rate, including video of flow through channels and digital elevation models acquired using small unoccupied aircraft systems, airborne lidar, and airborne single-pass interferometric synthetic aperture radar. Topographic data through time allowed calculation of subaerial lava flow volume and time-averaged discharge rate over the course of the eruption, which we integrated with pre- and post-eruption bathymetric surveys. Repeat videos of the near-vent channel were analyzed with particle velocimetry to extract flow velocities, and these were combined with open channel flow theory to calculate a time series of instantaneous effusion rates. Results show a general increase in dense rock equivalent (DRE) effusion rate from ~7 to ~100 m3/s from early to late May for the whole flow field and ≥ 200 m3/s by mid-June after the eruption had focused at a primary vent. By the end of the eruption in August, 0.9–1.4 km3 DRE of lava had erupted, with 0.4 km3 deposited on land and at least 0.5 km3 offshore. The trends in effusion rate through time reflect magmatic processes in the connected summit and rift zone system that controlled eruption rate, with resulting implications for lava flow dynamics and hazards.

The 2018 eruption of Kīlauea volcano produced the largest and most destructive lava flows in the lower East Rift Zone (LERZ) in the past 200 years. Average effusion rates exceeded 100 m3 s-1 (DRE) ...

Connecting caldera collapse The Kīlauea Volcano on the island of Hawai‘i erupted for 3 months in 2018. Neal et al. present a summary of the eruption sequence along with a variety of geophysical observations collected by the Hawaiian Volcano Observatory. The cyclic inflation, deflation, and eventual collapse of the summit was tied to lava eruption from lower East Rift Zone fissures. A total volume of 0.8 cubic kilometers of magma erupted, roughly the equivalent of 320,000 swimming pools, which matched the change in volume at the summit. Science , this issue p. 367

Unoccupied aircraft systems (UAS) are developing into fundamental tools for tackling the grand challenges in volcanology; here, we review the systems used and their diverse applications. UAS can typically provide image and topographic data at two orders of magnitude better spatial resolution than space-based remote sensing, and close-range observations at temporal resolutions down to those of video frame rates. Responsive deployments facilitate dense time-series measurements, unique opportunities for geophysical surveys, sample collection from hostile environments such as volcanic plumes and crater lakes, and emergency deployment of ground-based sensors (and robots) into hazardous regions. UAS have already been used to support hazard management and decision-makers during eruptive crises. As technologies advance, increasing system capabilities, autonomy and availability, supported by more diverse and lighter-weight sensors, will offer unparalleled potential for hazard monitoring. UAS will provide opportunities for pivotal advances in our understanding of complex physical and chemical volcanic processes.

INTRODUCTION The 2018 flank eruption and summit collapse of Kīlauea Volcano was one of the largest and most destructive volcanic events in Hawai‘i in the past 200 years. The eruption occurred on the volcano’s lower East Rift Zone (ERZ), draining magma at a high rate from the summit reservoir and triggering incremental collapse of the overlying caldera floor. Lava flows erupted for 3 months, destroying several residential subdivisions and burying miles of roads. The eruption rate exhibited cyclic behavior on multiple time scales, resulting in repeated lava breakouts and overflows. Multidisciplinary observations provide insight into the nature of these variations, driving forces in the magmatic system, and implications for hazard. RATIONALE Volumetric eruption rate is a primary control on the vigor and hazard of lava flows, but the processes that control its temporal variations are not well understood because of limited observational data. We integrated field observations, photos and video from time-lapse cameras and unmanned aircraft systems, seismic tremor, and infrasound to track the time scales and magnitude of fluctuations in eruption rate at the primary vent for the 2018 eruption on Kīlauea’s LERZ. We combined these data with documentation of summit caldera collapses to investigate the origins and impacts of these fluctuations. RESULTS Cyclic variations in eruption rate occurred on two disparate time scales. First, short-term fluctuations (“pulses”) in eruption vigor had periods of 5 to 10 min, but had no major implications for lava flow hazard. Flow rate in the lava channel was inversely related to fountaining and outgassing intensity at the vent. Second, long-term fluctuations (“surges”) had periods of 1 to 2 days and began within minutes of episodic caldera collapse events at the summit, 40 km upslope. These surges triggered overflows from the channel that produced hazardous enlargement of the lava flow field, which could be forecast several hours in advance. We also show that seismic tremor and infrasound were correlated with lava flow eruption rates. CONCLUSION We conclude that the two types of eruption rate cycles were controlled by two distinct processes. Short-term pulses were driven by changes in outgassing efficiency of the lava at shallow depths. Long-term surges in eruption rate were driven by pressure transients induced by the summit collapses and transmitted through the magma conduit over a distance of 40 km. The pressure-driven surges in eruption rate demonstrate that the episodic rhythm of summit caldera collapse sequences may be imparted on the accompanying flank eruption. The surges also help to constrain the efficient hydraulic connection between Kīlauea’s summit magma system and rift zones and demonstrate that pressure communication over distances of 40 km can occur on a time scale of minutes. Seismic tremor and infrasound may be effective proxies for lava flow eruption rates, allowing for improved tracking of lava flow hazards. Our multidisciplinary data provide a clear link between eruption rate fluctuations and their driving processes in the magmatic system.

Lava flows pose a significant hazard to infrastructure and property located close to volcanoes, and understanding how flows advance is necessary to manage volcanic hazard during eruptions. Compared to low-silica basaltic flows, flows of andesite composition are infrequently erupted and so relatively few studies of their characteristics and behaviour exist. We use El Reventador, Ecuador as a target to investigate andesitic lava flow properties during a 4.5 year period of extrusive eruption between February 2012 and August 2016. We use satellite radar to map the dimensions of 43 lava flows and look at variations in their emplacement behaviour over time. We find that flows descend the north and south flanks of El Reventador, and were mostly emplaced during durations shorter than the satellite repeat interval of 24 days.Flows ranged in length from 0.3 to 1.7 km, and the length of these flows decreased over the observation period. We measure a decrease in flow volume with time that is correlated with a long-term exponential decrease in eruption rate, and propose that this behaviour is caused by temporary magma storage in the conduit acting as a melt capacitor between the magma reservoir and the surface. We use the dimensions of the flow levees and widths to estimate the flow yield strengths, which were of the order of 10-100 kPa. We observe that some flows were diverted by topographic obstacles, and compare measurements of decreased channel width and increased flow thickness at the obstacles with observations from laboratory experiments. Radar observations, such as those presented here, could be used to map and measure properties of evolving lava flow fields at other remote or difficult to monitor volcanoes.

Laboratory experiments that simulate lava flows have been in use by volcanologists for many years. The behavior of flows in the lab, where “eruption” parameters, material properties, and environmental settings are tightly controlled, provides insight into the influence of various factors on flow evolution. A second benefit of laboratory lava flows is to provide a set of observations with which numerical models of flow emplacement can be tested. Models of lava flow emplacement vary in mathematical approach, physical assumptions, and computational cost. Nonetheless, all models require thorough testing and evaluation, and laboratory experiments produce an excellent test for models. This paper provides a primer on modern analog laboratory lava flow experiments. It reviews scaling con- siderations and provides quantitative information meant to guide future experimentalists in designing their experiments to be relevant to natural processes. Traditional and novel laboratory techniques are described, including a discussion of current limitations. New insights from recent experiments highlight the impact of topographic conditions and highlight the importance of considering bed roughness, major obstacles, and slope breaks. The influence of episodic or non-uniform effusion rate is demonstrated through recent experi- mental works. Lastly, the paper discusses several open questions about lava flow emplacement and the ways in which future improvements in experimental methods, such as the ability to utilize three-phase suspensions and materials with complex rheologies and to image the interior of flows could help answer these.

Mafic volcanic fields are widespread, but few have erupted in historic times, providing limited observations of the magnitudes, dynamics, and timescales of lava flow emplacement in these settings. To expand our knowledge of effusive mafic eruptions, we must evaluate solidified flows to discern syn-eruptive conditions. The Harrat Rahat volcanic field in western Saudi Arabia offers a good opportunity for this, with a historical eruption in 1256 CE and many well-preserved prehistoric flows. We combine historical observations and rheological and morphological analyses of the youngest flows with analytical models to reconstruct eruptive histories and lava flow emplacement conditions in Harrat Rahat. Petrologic analysis of samples for emplacement temperatures and crystallinities shows cooling trends from vent to toe of ~ 1140 to ~ 1090 °C at rates of 2–7 °C km−1, crystallinities increasing from 0.5 to 60%, and apparent viscosities increasing from 102 to 109 Pa s. High-resolution topographic data facilitates quantitative analysis of morphology and interpolation of pre-eruptive surfaces to measure flow thicknesses, channels, and levees, and enables calculation of eruptive volumes. Analytical models relating flow morphology to emplacement conditions are applied to estimate effusion rates. Within the suite of studied flows, volume estimates range from 0.07 to 0.42 km3 dense rock equivalent, with effusion rates on the order of 10 to 100 s of m3 s−1 and durations from 1 to 15 weeks. These integrated analyses quantify past lava flow emplacement conditions and dynamics in Harrat Rahat, improving our understanding and observations of fundamental parameters and controls of effusive eruptions in mafic volcanic fields.

A fundamental goal of igneous petrology is to quantify the duration of time required to produce evolved magmas following influx of basalt into the crust. However, in many cases, complex field relations and/or the presence of a long-lived magmatic system make it difficult to assess how basaltic inputs relate to more evolved magmas, therefore, precluding calculation of meaningful timescales. Here, we present field relations, geochemistry, 40Ar/39Ar ages, and 36Cl ages for volcanic rocks from the Harrat Rahat volcanic field, Saudi Arabia. These data document a systematic and repeated temporal progression from alkali basalt to trachyte for the youngest eruptives. From ~ 150 to ~ 17 ka the following eruptive sequence occurred four times: (1) alkali basalt, (2) hawaiite, mugearite, or benmoreite, and (3) trachyte. We interpret each eruptive sequence to result from injection of basalt into the crust, and its subsequent differentiation and eruption of progressively evolved magmas. We use the interval time between successive eruptions within a given sequence to calculate the duration of time required to produce trachyte from alkali basalt. Differentiation from alkali basalt to intermediate compositions (hawaiite, mugearite, benmoreite) took ≤ 2 kyr on average. Differentiation from intermediate compositions to trachyte took a maximum of 6.6 ± 3.5 to 22.5 ± 1.6 kyr. Thus, the total duration of differentiation was ~ 9 to ~ 25 kyr. Timescales presented here are insensitive to processes evoked to drive differentiation because they are based solely on the ages and compositions of eruptive products from a system characterized by a simple, repeated differentiation sequence.

Numerical simulations of lava flow emplacement are valuable for assessing lava flow hazards, forecasting active flows, designing flow mitigation measures, interpreting past eruptions, and understanding the controls on lava flow behavior. Existing lava flow models vary in simplifying assumptions, physics, dimensionality, and the degree to which they have been validated against analytical solutions, experiments, and natural observations. In order to assess existing models and guide the development of new codes, we conduct a benchmarking study of computational fluid dynamics (CFD) models for lava flow emplacement, including VolcFlow, OpenFOAM, FLOW-3D, COMSOL, and MOLASSES. We model viscous, cooling, and solidifying flows over horizontal planes, sloping surfaces, and into topographic obstacles. We compare model results to physical observations made during well-controlled analogue and molten basalt experiments, and to analytical theory when available. Overall, the models accurately simulate viscous flow with some variability in flow thickness where flows intersect obstacles. OpenFOAM, COMSOL, and FLOW-3D can each reproduce experimental measurements of cooling viscous flows, and OpenFOAM and FLOW-3D simulations with temperature-dependent rheology match results from molten basalt experiments. We assess the goodness-of-fit of the simulation results and the computational cost. Our results guide the selection of numerical simulation codes for different applications, including inferring emplacement conditions of past lava flows, modeling the temporal evolution of ongoing flows during eruption, and probabilistic assessment of lava flow hazard prior to eruption. Finally, we outline potential experiments and desired key observational data from future flows that would extend existing benchmarking data sets.

New high-resolution maps of Hawaiian lava flows highlight complex topographically controlled channel networks. Network geometries range from distributary systems dominated by branching around local obstacles, to tributary systems constricted by topography. We combine 2-D network analysis tools developed for river systems and neural networks with 3-D lidar morphologic analysis and historical records of flow emplacement to investigate both the origins of channel networks and their influence on flow morphology and behavior. We find that network complexity is a function of underlying slope and that the degree of flow branching, network connectivity, and longevity of individual channels all influence the final flow morphology (flow and channel widths and levee heights). Furthermore, because channel networks govern the distribution of lava supply within a flow, changes in the channel topology can dramatically alter the effective volumetric flux in any one branch, which affects both flow length and advance rate. Specifically, branching will slow and shorten flows, while merging can accelerate and lengthen them. Consideration of channel networks is thus important for predicting lava flow behavior and mitigating flow hazards with diversion barriers. Observed relationships between network geometry, flow parameters, and morphology also offer insight into the interpretation of these features elsewhere on Earth and other terrestrial planets.

While major explosive eruptions of anhydrite-bearing silicic magma have been known to influence earth's climate through the release of sulfate aerosols, quantifying the atmospheric loading of these eruptions remains an elusive goal. The eruptions that have had the most impact on climate typically erupted oxidized and anhydrite-bearing magma. The > 9 km3 of dacite magma erupted during the VEI 6 eruption of Huaynaputina (Peru) in 1600 is oxidized (ƒO2 NNO + 1.4) but does not contain magmatic anhydrite. Yet it is thought to have a significant impact on climate. Herein we follow the approach that S contents in magmatic apatite grains can be used to provide an independent estimate of the magnitude and source of the atmospheric sulfur loading. Apatite SO3 concentrations of 0.09 to 0.17 wt.% yield an estimate of ~0.8 Mt of total erupted sulfur derived from degassing of the melt. This is considerably less than the stratospheric sulfur input of 16–55 Mt of S estimated from ice core data and requires that much of the stratospheric sulfur load from Huaynaputina was from other sources than the melt. While a fluid phase in equilibrium with the melt is a likely source, a significant contribution from combustion of sulfur from a fossil hydrothermal system excavated during eruption is evidenced by comminuted hydrothermal debris and lithics that make up to 2.5 to 10 wt.% of erupted material. We calculate that up to 22% of the total S budget of the eruption could reasonably be attributed to combusted hydrothermal sulfur. Thus, three sources of sulfur: the melt, a coexisting fluid phase, and external hydrothermal sulfur, appear to explain the atmospheric input from this eruption. Our work confirms that unraveling the S budget of large silicic eruptions is extremely complex and requires a multi-faceted approach based on a detailed understanding of the eruption.

[1] Lava flow mapping is both an essential component of volcano monitoring and a valuable tool for investigating lava flow behavior. Although maps are traditionally created through field surveys, remote sensing allows an extraordinary view of active lava flows while avoiding the difficulties of mapping on location. Synthetic aperture radar (SAR) imagery, in particular, can detect changes in a flow field by comparing two images collected at different times with SAR coherence. New lava flows radically alter the scattering properties of the surface, making the radar signal decorrelated in SAR coherence images. We describe a new technique, SAR Coherence Mapping (SCM), to map lava flows automatically from coherence images independent of look angle or satellite path. We use this approach to map lava flow emplacement during the Pu‘u ‘Ō‘ō-Kupaianaha eruption at Kīlauea, Hawai‘i. The resulting flow maps correspond well with field mapping and better resolve the internal structure of surface flows, as well as the locations of active flow paths. However, the SCM technique is only moderately successful at mapping flows that enter vegetation, which is also often decorrelated between successive SAR images. Along with measurements of planform morphology, we are able to show that the length of time a flow stays decorrelated after initial emplacement is linearly related to the flow thickness. Finally, we use interferograms obtained after flow surfaces become correlated to show that persistent decorrelation is caused by post-emplacement flow subsidence.