More Like this

1. Stepping back to the source: The expression of eruption column collapse through submarine terraces

   Gilchrist, Johan; Jellinek, AM; Hooft, E; Wanket, S (December 2024, 2024 Fall Meeting, American Geophysical Union)

   Explosive eruption jets rising through relatively shallow water layers form eruption columns that can deliver volcanic ash, gases, and entrained water to the atmosphere and ocean in sequence or simultaneously, depending on eruption source parameters (Gilchrist et al. 2023). Despite the mesospheric eruption column height of the January 15, 2022 eruption of Hunga Tonga-Hunga Ha’apai (HTHH), the majority of erupted material was delivered to the surrounding seafloor via submarine pyroclastic density currents (PDCs). Deposits of HTHH show evidence of axisymmetric terraced deposits, which we show are linked to the mass eruption rate and dynamics of column collapse. We use scaled analog experiments on multiphase sand-water fountains injected into water layers of varying depth to model the collapse dynamics of shallow water eruption columns and to link fountain source conditions to deposit topography. The source strength of multiphase fountains predicts whether they collapse periodically or continuously via sedimentation waves with varying frequency and momentum. In turn, the frequency and momentum of sedimentation waves impacting the tank base determines whether ground-hugging gravity currents flowing out of the sedimentation wave impact zone are initially erosive or depositional. On the basis of experiments, we propose that syn-eruptive shallow submarine caldera deposits that show evidence of terracing and proximal scouring are linked to relatively strong eruption jets in the regime where the jet is in partial collapse or total collapse. In these regimes, the eruption jet collapses periodically as sedimentation waves that erode the deposit in the impact zone and transition into submarine PDCs that deposit the sedimentation wave mixture into regularly spaced terraces thereafter (Fig. 1, black boxes). In contrast, we expect weak eruption jets to occur in the total collapse regime where sedimentation waves descend in rapid succession and effectively supply submarine PDCs continuously which, in turn, build deposits lacking terraces (Fig. 1, blue box). For common values of caldera eruption source parameters, we link submarine PDC deposit morphology to eruption jet strength and plausible mass eruption rates. 

   more » « less
2. Runoff required to drive postimpact gully development on the walls of Meteor Crater (Arizona, USA)

   https://doi.org/10.1130/B36863.1  

   Palucis, Marisa C; Howard, Alan D; Kring, David A; Nishiizumi, Kuni; Caffee, Marc; Dietrich, William E (July 2023, Geological Society of America Bulletin)

   Since the impact ∼50,000 yr ago, surface runoff has entrained and transported sediment from the walls to the floor of Meteor Crater (Arizona, USA). Previous work interpreted this erosion and deposition to be due to predominantly fluvial (i.e., dilute water transport) processes. However, light detection and ranging (LiDAR)−derived topographic data and field observations indicate that debris flows dominated, which were likely generated by runoff that entrained the talus that borders bedrock cliffs high on the crater walls. The low gradient of the crater floor caused debris flows to stop, leaving lobate deposits, while fluvial processes delivered sediment toward the center of the crater. Cosmogenic radionuclide dating of levee deposits suggests that debris-flow activity ceased in the late Pleistocene, synchronous with regional drying. Assuming a rock-to-water ratio of 0.3 at the time of transport by mass flows, it would have taken ∼2 × 106 m3 of water to transport the estimated ∼6.8 × 106 m3 of debris-flow deposits found at the surface of the crater floor. This extensive erosion would require ∼6 m of total runoff over the 0.35 km2 upslope source area of the crater, or ∼18 mm of runoff per debris-flow event. Much more runoff did occur, as evidenced by crater lake deposits, Holocene fluvial activity (which produced little erosion), and contemporary rainfall rates. Rarely on Earth is the total amount of water that creates and runs through a landscape estimated, yet such calculations are commonly done on Mars. Our analysis suggests that erosional and depositional landforms may record only a small fraction of the total runoff. 

   more » « less
3. Magma mingling during the 1959 eruption of Kīlauea Iki, Hawaiʻi

   Marsh, Jennifer; Edmonds, Marie; Houghton, Bruce F; Buisman, I; Herd, R (May 2024, Bulletin of volcanology)

   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. The 1831 CE mystery eruption identified as Zavaritskii caldera, Simushir Island (Kurils)

   https://doi.org/10.1073/pnas.2416699122  

   Hutchison, William; Sugden, Patrick; Burke, Andrea; Abbott, Peter; Ponomareva, Vera V; Dirksen, Oleg; Portnyagin, Maxim V; MacInnes, Breanyn; Bourgeois, Joanne; Fitzhugh, Ben; et al (January 2025, Proceedings of the National Academy of Sciences)

   Polar ice cores and historical records evidence a large-magnitude volcanic eruption in 1831 CE. This event was estimated to have injected ~13 Tg of sulfur (S) into the stratosphere which produced various atmospheric optical phenomena and led to Northern Hemisphere climate cooling of ~1 °C. The source of this volcanic event remains enigmatic, though one hypothesis has linked it to a modest phreatomagmatic eruption of Ferdinandea in the Strait of Sicily, which may have emitted additional S through magma–crust interactions with evaporite rocks. Here, we undertake a high-resolution multiproxy geochemical analysis of ice-core archives spanning the 1831 CE volcanic event. S isotopes confirm a major Northern Hemisphere stratospheric eruption but, importantly, rule out significant contributions from external evaporite S. In multiple ice cores, we identify cryptotephra layers of low K andesite-dacite glass shards occurring in summer 1831 CE and immediately prior to the stratospheric S fallout. This tephra matches the chemistry of the youngest Plinian eruption of Zavaritskii, a remote nested caldera on Simushir Island (Kurils). Radiocarbon ages confirm a recent (<300 y) eruption of Zavaritskii, and erupted volume estimates are consistent with a magnitude 5 to 6 event. The reconstructed radiative forcing of Zavaritskii (−2 ± 1 W m⁻²) is comparable to the 1991 CE Pinatubo eruption and can readily account for the climate cooling in 1831–1833 CE. These data provide compelling evidence that Zavaritskii was the source of the 1831 CE mystery eruption and solve a confounding case of multiple closely spaced observed and unobserved volcanic eruptions. 

   more » « less
5. Turbulence, entrainment and low-order description of a transitional variable-density jet

   https://doi.org/10.1017/jfm.2017.822  

   Viggiano, B.; Dib, T.; Ali, N.; Mastin, L. G.; Cal, R. B.; Solovitz, S. A. (February 2018, Journal of Fluid Mechanics)

   Geophysical flows occur over a large range of scales, with Reynolds numbers and Richardson numbers varying over several orders of magnitude. For this study, jets of different densities were ejected vertically into a large ambient region, considering conditions relevant to some geophysical phenomena. Using particle image velocimetry, the velocity fields were measured for three different gases exhausting into air – specifically helium, air and argon. Measurements focused on both the jet core and the entrained ambient. Experiments considered relatively low Reynolds numbers from approximately 1500 to 10 000 with Richardson numbers near 0.001 in magnitude. These included a variety of flow responses, notably a nearly laminar jet, turbulent jets and a transitioning jet in between. Several features were studied, including the jet development, the local entrainment ratio, the turbulent Reynolds stresses and the eddy strength. Compared to a fully turbulent jet, the transitioning jet showed up to 50 % higher local entrainment and more significant turbulent fluctuations. For this condition, the eddies were non-axisymmetric and larger than the exit radius. For turbulent jets, the eddies were initially smaller and axisymmetric while growing with the shear layer. At lower turbulent Reynolds number, the turbulent stresses were more than 50 % higher than at higher turbulent Reynolds number. In either case, the low-density jet developed faster than a comparable non-buoyant jet. Quadrant analysis and proper orthogonal decomposition were also utilized for insight into the entrainment of the jet, as well as to assess the energy distribution with respect to the number of eigenmodes. Reynolds shear stresses were dominant in Q1 and Q3 and exhibited negligible contributions from the remaining two quadrants. Both analysis techniques showed that the development of stresses downstream was dependent on the Reynolds number while the spanwise location of the stresses depended on the Richardson number. 

   more » « less