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
Complex staged emplacement of a basaltic lava: The example of the July 1974 flow of Kīlauea
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
- Award ID(s):
- 2119838
- PAR ID:
- 10580480
- Publisher / Repository:
- Springer Science + Business Media
- Date Published:
- Journal Name:
- Bulletin of Volcanology
- Volume:
- 87
- Issue:
- 4
- ISSN:
- 1432-0819
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract On 3 May 2018, Kīlauea Volcano, one of the most active volcanoes in the world, entered a new eruptive phase because of a dike intrusion in the East Rift zone. One day later, an Mw 6.9 earthquake, which was likely trigged by the dike intrusion, occurred in the submarine south flank of Kīlauea Volcano. In mid-July, an ocean-bottom seismometer (OBS) array consisting of 12 stations was deployed on the submarine south flank of Kīlauea Volcano to monitor the aftershocks and lava–water interaction near the ocean entry. Eleven OBSs were recovered in mid-September. Preliminary evaluation of the data reveals a large number of seismic and acoustic events, which provide a valuable dataset for understanding flank deformation and stability as well as lava–water interaction. Here, we introduce this dataset and document notable instrument malfunctions along with some initial seismic and acoustic observations.more » « less
-
Abstract Preserved rocks in the Jurassic Ferrar Large Igneous Province consist mainly of intrusions, and extrusive rocks, the topic of this chapter, comprise the remaining small component. They crop out in a limited number of areas in the Transantarctic Mountains and southeastern Australia. They consist of thick sequences of lavas and sporadic occurrences of volcaniclastic rocks. The latter occur mainly beneath the lavas and represent the initial eruptive activity, but also are present within the lava sequence. The majority are basaltic phreatomagmatic deposits and in at least two locations form immense phreatocauldrons filled with structureless tuff breccias and lapilli tuffs with thicknesses of as much as 400 m. Stratified sequences of tuff breccias, lapilli tuffs and tuffs are up to 200 m thick. Thin tuff beds are sparsely distributed in the lava sequences. Lava successions are mainly 400–500 m thick, and comprise individual lavas ranging from 1 to 230 m thick, although most are in the range of 10–100 m. Well-defined colonnade and entablature are seldom displayed. Lava sequences were confined topographically and locally ponded. Water played a prominent role in eruptive activity, as exhibited by phreatomagmatism, hyaloclastites, pillow lava and quenching of lavas. Vents for lavas have yet to be identified.more » « less
-
Abstract Seafloor volcanic eruptions are difficult to directly observe due to lengthy eruption cycles and the remote location of mid‐ocean ridges. Volcanic eruptions in 2005–2006 at 9°50′N on the East Pacific Rise have been well documented, but the lava volume and flow extent remain uncertain because of the limited near‐bottom bathymetric data. We present near‐bottom data collected during 19 autonomous underwater vehicle (AUV)Sentrydives at 9°50′N in 2018, 2019, and 2021. The resulting 1 m‐resolution bathymetric grid and 20 cm‐resolution sidescan sonar images cover 115 km2, and span the entire area of the 2005–2006 eruptions, including an 8 km2pre‐eruption survey collected with AUVABEin 2001. Pre‐ and post‐eruption surveys, combined with sidescan sonar images and seismo‐acoustic impulsive events recorded during the eruptions, are used to quantify the lava flow extent and to estimate changes in seafloor depth caused by lava emplacement. During the 2005–2006 eruptions, lava flowed up to ∼3 km away from the axial summit trough, covering an area of ∼20.8 km2; ∼50% larger than previously thought. Where pre‐ and post‐eruption surveys overlap, individual flow lobes can be resolved, confirming that lava thickness varies from ∼1 to 10 m, and increases with distance from eruptive fissures. The resulting lava volume estimate indicates that ∼57% of the melt extracted from the axial melt lens probably remained in the subsurface as dikes. These observations provide insights into recharge cycles in the subsurface magma system, and are a baseline for studying future eruptions at the 9°50′N area.more » « less
-
Mocho-Choshuenco volcano (39.9°S, 72.0°W) produced ∼75 explosive eruptions following retreat of the >1.5-km-thick Patagonian Ice Sheet associated with the local Last Glacial Maximum (LGM, from 35 to 18 ka). Here, we extend this record of volcanic evolution to include pre- and syn-LGM lavas that erupted during the Pleistocene. We establish a long-term chronology of magmatic and volcanic evolution and evaluate the relationship between volcanism and loading/unloading of the Patagonian Ice Sheet via twenty-four 40Ar/39Ar and two 3He age determinations integrated with stratigraphy and whole-rock compositions of lava flows and glass compositions of tephra. Our findings reveal that the edifice is much younger than previously thought and preserves 106 km3 of eruptive products, of which 50% were emplaced immediately following the end of the penultimate glaciation and 20% after the end of the LGM. A period of volcanic inactivity between 37 and 26 ka, when glaciers expanded, was followed by the eruption of incompatible element-rich basaltic andesites. Several of these syn-LGM lavas dated between 26 and 16 ka, which crop out at 1500−1700 m above sea level, show ice contact features that are consistent with emplacement against a 1400- to 1600-m-thick Patagonian Ice Sheet. Small volume dacitic eruptions and two explosive rhyolitic eruptions dominate the volcanic output from 18 to 8 ka, when the Patagonian Ice Sheet began to retreat rapidly. We hypothesize that increased lithostatic loading as the Patagonian Ice Sheet grew prohibited dike propagation, thus stalling the ascent of magma, promoting growth of at least three discrete magma reservoirs, and enhancing minor crustal assimilation to generate incompatible element-rich basaltic andesitic to dacitic magmas that erupted between 26 and 17 ka. From an adjacent reservoir, incompatible element-poor dacites erupted from 17 to 12 ka. These lava flows were followed by the caldera-forming eruption at 11.5 ka of 5.3 km3 of rhyolite from a deeper reservoir atop which a silicic melt lens had formed and expanded. Subsequent eruptions of oxidized dacitic magmas from the Choshuenco cone from 11.5 to 8 ka were followed by andesitic to dacitic eruptions at the more southerly Mocho cone, as well as small flank vent eruptions of basaltic andesite at 2.5 and 0.5 ka. This complex history reflects a multi-reservoir plumbing system beneath Mocho-Choshuenco, which is characterized by depths of magma storage, oxidation states, and trace element compositions that vary over short periods of time (<2 k.y.).more » « less
An official website of the United States government
