Abstract. Theoretical source models of underwater explosions are often applied in studying tsunami hazards associated with subaqueous volcanism; however, their use in numerical codes based on the shallow water equations can neglect the significant dispersion of the generated wavefield. A non-hydrostatic multilayer method is validated against a laboratory-scale experiment of wave generation from instantaneous disturbances and at field-scale subaqueous explosions at Mono Lake, California, utilising the relevant theoretical models. The numerical method accurately reproduces the range of observed wave characteristics for positive disturbances and suggests a relationship of extended initial troughs for negative disturbances at low-dispersivity and high-non-linearity parameters. Satisfactory amplitudes and phase velocities within the initial wave group are found using underwater explosion models at Mono Lake. The scheme is then applied to modelling tsunamis generated by volcanic explosions at Lake Taupō, New Zealand, for a magnitude representing an ejecta volume of 0.1 km3. Waves reach all shores within 15 min with maximum incident crest amplitudes around 0.2 m at shores near the source. This work shows that the multilayer scheme used is computationally efficient and able to capture a wide range of wave characteristics, including dispersive effects, which is necessary when investigating subaqueous explosions. This research therefore provides the foundation for future studies involving a rigorous probabilistic hazard assessment to quantify the risks and relative significance of this tsunami source mechanism.

Abstract. Theoretical source models of underwater explosions are often applied in studying tsunami hazards associated with submarine volcanism; however, their use in numerical codes based on the shallow water equations can neglect the significant dispersion of the generated wavefield. A non-hydrostatic multilayer method is validated against a laboratory-scale experiment of wave generation from instantaneous disturbances and at field-scale submarine explosions at Mono Lake, California, utilising the relevant theoretical models. The numerical method accurately reproduces the range of observed wave characteristics for positive disturbances and suggests a previously unreported relationship of extended initial troughs for negative disturbances at low dispersivity and high nonlinearity parameters. Satisfactory amplitudes and phase velocities within the initial wave group are found using underwater explosion models at Mono Lake. The scheme is then applied to modelling tsunamis generated by volcanic explosions at Lake Taupō, New Zealand, for a magnitude range representing ejecta volumes between 0.04–0.4 km3. Waves reach all shores within 15 minutes with maximum incident crest amplitudes around 4 m at shores near the source. This work shows that the multilayer scheme used is computationally efficient and able to capture a wide range of wave characteristics, including dispersive effects, which is necessary when investigating submarine explosions. This research therefore provides the foundation for future studies involving a rigorous probabilistic hazard assessment to quantify the risks and relative significance of this tsunami source mechanism.

Abstract. Volcanogenic tsunami and wave hazard remains less understood than that of other tsunami sources. Volcanoes can generate waves in a multitude of ways, including subaqueous explosions. Recent events, including a highly explosive eruption at Hunga Tonga-Hunga Ha'apai and subsequent tsunami in January 2022, have reinforced the necessity to explore and quantify volcanic tsunami sources. We utilise a non-hydrostatic multilayer numerical method to simulate 20 scenarios of sublacustrine explosive eruptions under Lake Taupō, New Zealand, across five locations and four eruption sizes. Waves propagate around the entire lake within 15 min, and there is a minimum explosive size required to generate significant waves (positive amplitudes incident on foreshore of > 1 m) from the impulsive displacement of water from the eruption itself. This minimum size corresponds to a mass eruption rate of 5.8×107 kg s−1, or VEI 5 equivalent. Inundation is mapped across five built areas and becomes significant near shore when considering only the two largest sizes, above VEI 5, which preferentially impact areas of low-gradient slope. In addition, novel hydrographic output is produced showing the impact of incident waves on the Waikato River inlet draining the lake and is potentially useful for future structural impact analysis. Waves generated from these explosive source types are highly dispersive, resulting in hazard rapidly diminishing with distance from the source. With improved computational efficiency, a probabilistic study could be formulated and other, potentially more significant, volcanic source mechanisms should be investigated.

No abstract