skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on December 1, 2025

Title: Stepping back to the source: The expression of eruption column collapse through submarine terraces
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
Award ID(s):
2023338
PAR ID:
10567335
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
2024 Fall Meeting, American Geophysical Union
Date Published:
Format(s):
Medium: X
Location:
https://agu.confex.com/agu/agu24/meetingapp.cgi/Paper/1519009
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract On 15 January 2022, Hunga volcano erupted, creating an extensive and high-reaching umbrella cloud over the open ocean, hindering traditional isopach mapping and fallout volume estimation. In MODIS satellite imagery, ocean surface water was discolored around Hunga following the eruption, which we attribute to ash fallout from the umbrella cloud. By relating intensity of ocean discoloration to fall deposit thicknesses in the Kingdom of Tonga, we develop a methodology for estimating airfall volume over the open ocean. Ash thickness measurements from 41 locations are used to fit a linear relationship between ash thickness and ocean reflectance. This produces a minimum airfall volume estimate of$${1.8}_{-0.4}^{+0.3}$$ 1.8 - 0.4 + 0.3 km3. The whole eruption produced > 6.3 km3of uncompacted pyroclastic material on the seafloor and a caldera volume change of 6 km3DRE. Our fall estimates are consistent with the interpretation that most of the seafloor deposits were emplaced by gravity currents rather than fall deposits. Our proposed method does not account for the largest grain sizes, so is thus a minimum estimate. However, this new ocean-discoloration method provides an airfall volume estimate consistent with other independent measures of the plume and is thus effective for rapidly estimating fallout volumes in future volcanic eruptions over oceans. 
    more » « less
  2. Esposti Ongaro, T. (Ed.)
    A multiphase fluid dynamic model is used to explore the effects of entrainment of granular debris into sustained volcanic jets such as those which produce sub-Plinian to Plinian eruption columns. The debris may be sourced from processes such as avalanches from crater walls or from recycling of previously erupted material. The results indicate that debris is not immediately, homogeneously mixed into a jet but instead forms a dense sheath that is dragged upward around the jet margin. While very small volumes of debris relative to the eruptive discharge rate mix progressively into the jet with increasing altitude, the dense sheath can inhibit entrainment of air into the lower portions of the jet, which may explain signs of column instability such as increased stratification in fallout deposits where lithic content increases. As debris volume increases, the dense sheath can collapse from a range of elevations to feed pyroclastic currents. The presence of the sheath of entrained debris contradicts some assumptions such as the top-hat profile for density and velocity that is commonly used in 1-D models. Transitions from fallout-producing buoyant column to collapsing behavior can be related to debris entrainment without any changes in primary eruption parameters such as vent size, exit velocity, or gas content. Boiling-over behavior can also be caused by debris entrainment, including recycling of previously erupted material such as might occur in a crater with restricted outlet. When entrained debris is relatively fine-grained such that it can couple well with the erupting mixture, complex, highly transient overpressured jet processes can occur due to the pinching effect of debris flowing into the base of the jet. Increasingly coarse debris causes collimation of the jet within the sheath of entrained material. The results suggest that accounting for the effects of debris entrainment is likely important for theoretical assessment of many natural eruption sequences and for assessment of hazard scenarios for potential sub-Plinian to Plinian activity. 
    more » « less
  3. Tsunamis from volcanic ‘explosive’ eruptions are rare, with the last catastrophic event being Krakatau in 1883 (Verbeek, 1885), during which, tsunamis were generated in the far-field by pressure shock-waves and in the nearfield of the volcano, in the Sunda Straits, by several potential geological mechanisms including pyroclastic flows, ash column, and/or caldera collapse. On 1/22/55, at about 4:15 UTC, a one in 1,000 year eruption of the Hunga Tonga-Hunga Ha’a-pai Volcano (HTHHV), that had started on12/20/21, reached its paroxysm with a series of large underwater explosions, releasing enormous energy (4-18 Mt of TNT), and ejecting a large ash plume 58 km into the stratosphere. We simulate both the near- and far-field tsunami generation from the eruption, but in this paper we focus on analyzing and validating the near-field impact against field data. 
    more » « less
  4. Abstract The April 2021 La Soufrière of St Vincent eruption generated several pyroclastic density currents (PDCs) during the 2 weeks of the crisis, from 9 to 22 April. To support the hazard assessment team during this eruption, numerical simulations were performed in real time and generated rapid scenario-based PDC invasion maps with the two-phase version of the code VolcFlow, which was able to simulate both the concentrated and dilute regime of PDCs. To generate the maps, only the source properties (shape and location) and the initial volume used to generate the PDCs were varied, all other input parameters were kept constant and estimated from previous simulations. New simulations were then performed based on the field-based deposit map to assess the code's ability to simulate such PDCs. Results show that the syn-crisis invasion maps satisfactorily mimic the observed valley-confined PDCs, while unconfined dilute PDCs were overestimated. The results also highlight that simulation results are greatly improved with additional field-based data, which help constrain the PDC sequence. Numerous lessons were learned, including (1) how to choose the most critical input parameters, (2) the importance of syn-eruptive radar imagery and (3) the potential of this two-phase model for rapid hazard assessment purposes. 
    more » « less
  5. Despite their global societal importance, the volumes of large-scale volcanic eruptions remain poorly constrained. Here, we integrate seismic reflection and P-wave tomography datasets with computed tomography-derived sedimentological analyses to estimate the volume of the iconic Minoan eruption. Our results reveal a total dense-rock equivalent eruption volume of 34.5 ± 6.8 km³, which encompasses 21.4 ± 3.6 km³ of tephra fall deposits, 6.9 ± 2 km³ of ignimbrites, and 6.1 ± 1.2 km³ of intra-caldera deposits. 2.8 ± 1.5 km³ of the total material consists of lithics. These volume estimates are in agreement with an independent caldera collapse reconstruction (33.1 ± 1.2 km³). Our results show that the Plinian phase contributed most to the distal tephra fall, and that the pyroclastic flow volume is significantly smaller than previously assumed. This benchmark reconstruction demonstrates that complementary geophysical and sedimentological datasets are required for reliable eruption volume estimates, which are necessary for regional and global volcanic hazard assessments. 
    more » « less