Jie Wu

The 2022 Hunga volcano eruption generated one of the largest and most destructive volcanic tsunamis ever recorded. We interpret the timing, nature and relative hazard of volcanic-tsunamigenic processes by combining geophysical, oceanographic, atmospheric, and eye-witness observations of this event. Hydro-acoustic signal ( T waves) detected on seismic stations up to 2600 km away provided a key to discerning tsunami triggers from explosions, volcanic mass flows, and caldera subsidence. The largest tsunami was generated by the first of at least two sudden caldera-subsidence events. This occurred 1 hour after globally detected explosions that produced smaller tsunami. Our results provide rare insights into the rates of caldera subsidence and highlight the complex multiple sources of tsunamis resulting from eruptions of oceanic volcanoes.

We examine the connected history of dacite-dominant volcanic rocks of the Tschicoma Formation, erupted between 5.5 and 2 Ma from the Jemez Mountains volcanic field, western USA. Zircon samples from two separate eruptions have continuous SHRIMP U–Pb age spectra spanning 0.84–1.08 Myr duration (3.12–3.96 Ma and 3.50–4.58 Ma, respectively), following an episode of zircon crystallization 0.28–0.50 Myr earlier (at 4.46 Ma and 4.86 Ma, respectively). Zircon chemical variations, as well as ubiquitous resorption textures that commonly show large core-rim age differences (up to 720–740 kyr), suggest that they grew in separate melt lenses. Zircons were likely stored at near-solidus or even sub-solidus conditions after crystallization, but may have been reactivated in response to at least four major magma recharge events every 300–400 kyr and smaller events in between. A cycle of zircon dissolution (from heating), recrystallization (during cooling), and storage repeated in different locations in the Tschicoma mush system throughout its lifespan; each recharge-induced heating stage may last for several hundred to more than a thousand years based on calculations of zircon dissolution. We envisage the melt lenses to be distributed in a crystal mush zone, coalescing into a single magma batch as magma recharge occurs shortly before eruption. Once active, increasing magma supply rates may trigger large-scale partial melting of the pre-existing mush and caldera-forming eruptions.

The 15 Jan 2022 eruption of Hunga-Tonga-Hunga-Ha’apai was the largest explosive volcanic event in the last 30 years. These islands represent the subaerially exposed summit of the Hunga Volcano, merged into a single land mass during the most recent eruption in 2014-2015. The 2022 eruption likely represents a 1-in-1000-year event for the Hunga Volcano, with the previous large-magnitude eruption occurring in ~1100 CE during a series of caldera-forming events. The 2022 erupted magma is plagioclase-, orthopyroxene- and clinopyroxene-bearing basaltic andesite to andesite dominated by blocky, poorly vesicular glassy ash with lesser amounts of vesicular pumiceous ash and fine lapilli. Melt Inclusions (MIs) hosted in plagioclase, clinopyroxene and orthopyroxene are abundant and glassy, some displaying shrinkage bubbles, with no evidence of secondary crystallization along the walls or within the MI glass. The groundmass glass and MI in the three main phenocryst phases were analysed for major, trace and volatile element concentrations to enable identification of magmatic sources and to better constrain processes happening at depth. Preliminary data indicate that plagioclase phenocrysts range from An93 to An78, and MI range from 54.1 to 58.7 wt % SiO2, with MgO from 2.5 to 5.3 wt %. Clinopyroxene phenocrysts range from En42 to En50, and MI range from 51.6 to 65.1 wt % SiO2, with MgO from 1.1 to 5.7 wt %. Orthopyroxene phenocrysts range from En68 to En77, and MI range from 55.7 to 59.6 wt % SiO2, with MgO from 2.5 to 5.3 wt %. Clinopyroxene MI span the full range of SiO2 compositions observed from the Hunga Volcano, from the host 2022 event (SiO2: ~57.5 wt %), the 1100 CE event (SiO2: ~60 wt %), the 2014-2015 event (SiO2: ~60.5 wt %), and the most evolved 2009 event (SiO2: ~63 wt %) and extend an additional ~4 wt % SiO2 to more mafic compositions. Orthopyroxene MI most closely resemble the 1100 CE event and the average groundmass glass compositions of the 2022 event. Plagioclase MI overlap the least silicic compositions observed in the 2022 groundmass glass (58.6 wt% SiO2) and extend down to 54 wt % SiO2, overlapping the main field of clinopyroxene MI. Both plagioclase and clinopyroxene MI tend to show higher MgO as compared to the 2022 groundmass glass at the same SiO2 concentration, whereas orthopyroxene shows lower MgO than the groundmass glass. SO3 in MI ranges up to 1600 ppm, significantly higher than the 2022 groundmass glass which averages 200 ppm, with both plagioclase and clinopyroxene MI preserving the highest observed concentrations. In contrast, Cl concentrations in MI extend to 2000 ppm, with the highest values in orthopyroxene and clinopyroxene, and plagioclase MI are lower and generally overlie the main groundmass glass concentrations (~1300 ppm). F was below detection limits. We postulate that clinopyroxene crystals reflect a more primitive basaltic andesite magma, whereas orthopyroxene crystals were likely derived from the magmatic remnants of the 2009 and 2014/2015 events in the upper magma system, and plagioclase crystals were sourced from the full range of magma sources.