Attention:The NSF Public Access Repository (PAR) system and access will be unavailable from 5:00 PM ET until 8:00 PM ET on Friday, September 11 due to maintenance. We apologize for the inconvenience.


Title: A discontinuous Galerkin finite-element model for fast channelized lava flows v1.0
Abstract. Lava flows present a significant natural hazard to communities around volcanoes and are typically slow-moving (<1 to 5 cm s−1) and laminar. Recent lava flows during the 2018 eruption of Kīlauea volcano, Hawai'i, however, reached speeds as high as 11 m s−1 and were transitional to turbulent. The Kīlauea flows formed a complex network of braided channels departing from the classic rectangular channel geometry often employed by lava flow models. To investigate these extreme dynamics we develop a new lava flow model that incorporates nonlinear advection and a nonlinear expression for the fluid viscosity. The model makes use of novel discontinuous Galerkin (DG) finite-element methods and resolves complex channel geometry through the use of unstructured triangular meshes. We verify the model against an analytic test case and demonstrate convergence rates of P+1/2 for polynomials of degree 𝒫. Direct observations recorded by unoccupied aerial systems (UASs) during the Kīlauea eruption provide inlet conditions, constrain input parameters, and serve as a benchmark for model evaluation.  more » « less
Award ID(s):
1654588
PAR ID:
10321175
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Geoscientific Model Development
Volume:
14
Issue:
6
ISSN:
1991-9603
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Basaltic lava flows can be highly destructive. Forecasting the future path and/or behavior of an active lava flow is challenging because topography is often poorly constrained and lava has a complex rheology and emplacement history. Preserved lavas are an important source of information which, combined with observations of active flows, underpins conceptual models of lava flow emplacement. However, the value of preserved lavas is limited because pre-eruptive topography and, thus, syn-eruptive lava flow geometry are usually not known. Here, we use tree-mold data to constrain pre-eruptive topography and syn-eruptive lava flow geometry of the July 1974 flow of Kīlauea (USA). Tree molds, which are formed after advancing lava encloses standing trees, preserve the lava inundation height and the final preserved thickness of lava. We used data from 282 tree molds to reconstruct the temporal and spatial evolution of the ~ 2.1 km-long July 1974 flow. The tree mold dataset yields a detailed dynamic picture of staged emplacement, separated by intervals of ponding. In some ponded areas, flow depth during emplacement (~ 5 m) was twice the preserved thickness of the final lava (2–3 m). Drainage of the ponds led to episodic surges in flow advancement, decoupled from fluctuations in vent discharge rate. We infer that the final breakout occurred after the cessation of fountaining. Such complex emplacement histories may be common for pāhoehoe lavas at Kīlauea and elsewhere in situations where the terrain is of variable slope, and/or where lava is temporarily perched and stored. 
    more » « less
  2. Waight, Tod (Ed.)
    Abstract The 2018 eruption of Kīlauea volcano in its Lower East Rift Zone began with the discharge of evolved high-Ti basalt as weak lava fountains and short, slow-moving lava flows. The lavas were quickly geochemically recognized as being derived from magmas stored within the rift zone and remobilized by a new intrusion, a sequence that is common at Kīlauea. This initial phase of the 2018 eruption, referred to as phase 1a, lasted for 6 days and was followed by extrusion of mixed magma after a 3-day pause. Even though remobilization of older rift zone magmas is common within Kīlauea’s rift zones, it is difficult to determine which past intrusion(s) may have initially emplaced those stored magmas. This difficulty stems from the tendency for Kīlauea magmas to follow very similar differentiation paths without significant variations in major, minor, or even trace element chemistry. We investigate possible magma sources for the lavas erupted during phase 1a of the 2018 eruption using whole-rock, mineral, and glass major and trace element compositions from historical East Rift Zone eruptions with adjacent fissures. We consider two primary hypotheses for the phase 1a source: magmas associated with the 1955 Lower East Rift Zone eruption or the nine eruptions in the Middle and Upper East Rift Zone during the 1960s. Our results suggest that magma associated with the earliest phases of Kīlauea’s 1955 eruption was the most likely source of the 2018 phase 1a remobilized magma. We determine volatile saturation pressures from melt inclusion chemistry and find similar storage depths for the 2018 phase 1a and early 1955 magmas. The phase 1a and early 1955 lavas are nearly indistinguishable in all of the compositional criteria considered, implying that the leftover 1955 magma body barely cooled and differentiated in the 63 years between eruptions (cooling rates of ~0.1 °C/year). This study sheds light on the potential for protracted storage of eruptible magmas in rift zones at Kīlauea, and highlights some of the challenges and solutions to identifying genetic relationships between magmas at Kīlauea. 
    more » « less
  3. Magma mingling and mixing are common processes at basaltic volcanoes and play a fundamental role in magma petrogenesis and eruption dynamics. Mingling occurs most commonly when hot primitive magma is introduced into cooler magma. Here, we investigate a scenario whereby cool, partially degassed lava is drained back into a conduit, where it mingles with hotter, less degassed magma. The 1959 eruption of Kīlauea Iki, Hawaiʻi involved 16 high fountaining episodes. During each episode, fountains fed a lava lake in a pit crater, which then partially drained back into the conduit during and after each episode. We infer highly crystalline tachylite inclusions and streaks in the erupted crystal-poor scoria to be the result of the recycling of this drain-back lava. The crystal phases present are dendrites of plagioclase, augite and magnetite/ilmenite, at sizes of up to 10 μm. Host sideromelane glass contains 7–8 wt% MgO and the tachylite glass (up to 0.5% by area) contains 2.5–6 wt% MgO. The vesicle population in the tachylite is depleted in the smallest size classes (< 0.5 mm) and has overall lower vesicle number densities and a higher degree of vesicle coalescence than the sideromelane component. The tachylite exhibits increasingly complex ‘stretching and folding’ mingling textures through the episodes, with discrete blocky tachylite inclusions in episodes 1 and 3 giving way to complex, folded, thin filaments of tachylite in pyroclasts erupted in episodes 15 and 16. We calculate that a lava lake crust 8–35 cm thick may have formed in the repose times between episodes, and then foundered and been entrained into the conduit during drain-back. The recycled fragments of crust would have been reheated in the conduit, inducing glass devitrification and crystallisation of pyroxene, magnetite and plagioclase dendrites and eventually undergoing ductile flow as the temperature of the fragments approached the host magma temperature. We use simple models of magma mingling to establish that stretching and folding of recycled, ductile lava could involve thinning of the clasts by up to a factor of 10 during the timescale of the eruption, consistent with observations of streaks and filaments of tachylite erupted during episodes 15 and 16, which may have undergone multiple cycles of eruption, drain-back and reheating. 
    more » « less
  4. Changes in magma chemistry that affect eruptive behavior occur during many volcanic eruptions, but typical analytical techniques are too slow to contribute to hazard monitoring. We used rapid energy-dispersive x-ray fluorescence analysis to measure diagnostic elements in lava samples within a few hours of collection during the 2018 Kīlauea eruption. The geochemical data provided important information for field crews and civil authorities in advance of changing hazards during the eruption. The appearance of hotter magma was recognized several days before the onset of voluminous eruptions of fast-moving flows that destroyed hundreds of homes. We identified, in near real-time, interactions between older, colder, stored magma—including the unexpected eruption of andesite—and hotter magma delivered during dike emplacement. 
    more » « less
  5. The science of volcanology advances disproportionately during exceptionally large or well-observed eruptions. The 2018 eruption of Kīlauea Volcano (Hawai‘i) was its most impactful in centuries, involving an outpouring of more than one cubic kilometer of basalt, a magnitude 7 flank earthquake, and the volcano's largest summit collapse since at least the nineteenth century. Eruptive activity was documented in detail, yielding new insights into large caldera-rift eruptions; the geometry of a shallow magma storage-transport system and its interaction with rift zone tectonics; mechanisms of basaltic tephra-producing explosions; caldera collapse mechanics; and the dynamics of fissure eruptions and high-volume lava flows. Insights are broadly applicable to a range of volcanic systems and should reduce risk from future eruptions. Multidisciplinary collaboration will be required to fully leverage the diversity of monitoring data to address many of the most important outstanding questions. ▪ Unprecedented observations of a caldera collapse and coupled rift zone eruption yield new opportunities for advancing volcano science. ▪ Magma flow to a low-elevation rift zone vent triggered quasi-periodic step-like collapse of a summit caldera, which pressurized the magma system and sustained the eruption. ▪ Kīlauea's magmatic-tectonic system is tightly interconnected over tens of kilometers, with complex feedback mechanisms and interrelated hazards over widely varying time scales. ▪ The eruption revealed magma stored in diverse locations, volumes, and compositions, not only beneath the summit but also within the volcano's most active rift zone. Expected final online publication date for the Annual Review of Earth and Planetary Sciences, Volume 52 is May 2024. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates. 
    more » « less