Scenario&ndash;based Modelling of Waves Generated by Sublacustrine Explosive Eruptions at Lake Taupō, New Zealand

Hayward, Matthew W.; Lane, Emily M.; Whittaker, Colin N.; Leonard, Graham S.; Power, William L.

doi:https://doi.org/10.5194/egusphere-2022-1152

Preprints

https://doi.org/10.5194/egusphere-2022-1152

Preprints

21 Nov 2022

| 21 Nov 2022

Scenario–based Modelling of Waves Generated by Sublacustrine Explosive Eruptions at Lake Taupō, New Zealand

Matthew W. Hayward, Emily M. Lane, Colin N. Whittaker, Graham S. Leonard, and William L. Power

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 minutes, 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 corresponds to a mass eruption rate of 5.8x10⁷ kg s^-1, or VEI (Volcanic Explosivity Index) 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 run-up. 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.

Received: 25 Oct 2022 – Discussion started: 21 Nov 2022

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this preprint. The responsibility to include appropriate place names lies with the authors.

Download & links

Preprint (PDF, 51179 KB)

Notice on discussion status
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
Preprint (51179 KB)

Download & links

The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.

Journal article(s) based on this preprint

03 Mar 2023

Scenario-based modelling of waves generated by sublacustrine explosive eruptions at Lake Taupō, New Zealand

Matthew W. Hayward, Emily M. Lane, Colin N. Whittaker, Graham S. Leonard, and William L. Power

Nat. Hazards Earth Syst. Sci., 23, 955–971, https://doi.org/10.5194/nhess-23-955-2023,https://doi.org/10.5194/nhess-23-955-2023, 2023

Short summary

Matthew W. Hayward, Emily M. Lane, Colin N. Whittaker, Graham S. Leonard, and William L. Power

Interactive discussion

Status: closed

RC1:
'Comment on egusphere-2022-1152', Raphael Paris, 09 Dec 2022

Dear colleagues,

Congratulations for this original work on tsunami hazard from underwater explosions in Lake Taupo. The manuscript represents a substantial contribution to hazard assesment in the Laker Taupo area, and the recent tsunami observed on November 30, 2022 shows how this study is timely (although the source of the tsunami is not yet elucidated). The overall quality of the scientific and methodological approach is very good, and the results are presented in a clear and concise way. The manuscript is well-written and well-illustrated. I have only some minor comments that are listed below:

- The use of the Volcanic Explosivity Index VEI (e.g. at lines 8-9, 264, 266, 320-324, 371) is not fully appropriate here, because it is based only on the total erupted volume. Here you are dealing with short discrete explosions. You can keep the VEI values in Table 1, for information. However I would recommend refering to the energy or mass eruption rates in the text.

- I suggest adding a sentence on the tectonic setting, which is too briefly mentioned at line 24: "The 350 km-long TVZ forms the southern part of a back-arc basin behind the Tonga-Kermadec subduction zone. Most of the eruptive centres are aligned along the intra-arc Taupo rift."

- The TVZ cannot be considered as a purely silicic magmatic system, because the magmas erupted range from andesites to rhyolites.

- At line 169 you state that "the water and explosion depths are equivalent." This is valid only in the case of a new vent formed on flat lake bottom. If the explosion comes from a pre-existing submarine cone, then explosion depth and water depth are not equivalents. This statement should therefore be qualified.

- I suggest slightly reorganizing the paper as follows:

3. Results

3.1. Tsunami propagation and wave heights at the coast

3.2. Tsunami inundation and potential impact on infratructures

- In the discussion, we lack a short discussion (at line 324?) on the probabilities mentioned on Table 1, in the light of the new results obtained here.

- The scenario presented at lines 330-334 is very pragmatic and relevant. Indeed in the case of such an eruption, all people would be probably evacuated before a tsunami could impact them. Could you add a short sentence of the consequences in terms of evacuation policies?

- In the conclusion, maybe you could formulate some recommendations for a local tsunami warning system in Lake Taupo?

Line-by-line corrections/suggestions:

-line 7: "This minimum size corresponds to..."

-9: "slope" rather than "runu-up"?

-37: "(HTHH, January 2022)"

-37: You could cite other references on the HTHH tsunami (e.g. Carvajal et al., 2002; Omira et al., 2022).

-41: cite also Maeno & Imamura (2011) for the 1883 Krakatau tsunami debate.

-56: "and in the present work it is..."

-59-63: this sentence is perhaps too long.

-93 refined against what?

-127: 7% of what?

-128-129: not sure if these lines are really useful here.

-194: "'We selected the five..."

-230-232: again, this sentence is a little bit long.

-248: refer to Fig 6 (c-d).

-257: 1 m instead of 1 cm?

-340: "asteroid-ocean impact, megatsunami from ocean-island flank collapse, and the recent tsunami..."

-345: "over half a decade old " -> a reference is needed here.

- 350-352: I was quite surprised that you don't cite here the work of Shen et al. (2021a, 2021b).

-367: "of what any tsunami hazard"

-384: What is the address of the website where Basilisk can be found?

Suggestions on figures:

- change the color of the eruption sites on figure 2 (white or yellow?).

- Fig 4: add source location (eruption sites) on the maps (same for fig 5).

- On the left side of fig 4, could you please mention the location number (1 to 5) as in Table 2, for more clarity (same for fig 5).

- Fig 6: add a dotted line to indicate the 1 m threshold value?

- Fig 6 c) d): in the caption, indicate that the first numbers refer to the scenarios and the name to the source locations 1 to 5.

- Fig. 8: Indicate to which closest source it corresponds on the different maps (same for fig 9).

Suggested references:

- Carvajal, M., Sepúlveda, I., Gubler, A., Garreaud, R.(2022). Worldwide Signature of the 2022 Tonga Volcanic Tsunami. Geophysical Research Letters, 49 (6), art. no. e2022GL098153; doi: 10.1029/2022GL098153

- Omira, R., Ramalho, R.S., Kim, J., González, P.J., Kadri, U., Miranda, J.M., Carrilho, F., & Baptista, M.A. (2022). Global Tonga tsunami explained by a fast-moving atmospheric source. Nature; doi.org/10.1038/s41586-022-04926-4

- Maeno, F., Imamura, F., 2011. Tsunami generation by a rapid entrance of pyroclastic flow into the sea during the 1883 Krakatau eruption, Indonesia. Journal of Geophysical Research 116, B09205

- Shen, Y., Whittaker, C.N., Lane, E.M., White, J.D.L., Power, W., Nomikou, P., 2021a. Laboratory experiments on tsunamigenic discrete subaqueous volcanic eruptions. Part 1: Free surface disturbances. Journal of Geophysical Research, Oceans 126, e2020JC016588.

- Shen, Y., Whittaker, C.N., Lane, E.M., White, J.D.L., Power, W., Nomikou, P., 2021a. Laboratory experiments on tsunamigenic discrete subaqueous volcanic eruptions. Part 2: Properties of generated waves. Journal of Geophysical Research, Oceans 126, e2020JC016587.

Citation: https://doi.org/10.5194/egusphere-2022-1152-RC1
- AC1: 'Reply on RC1', Matthew Hayward, 08 Jan 2023
  
  Thank you for the very thorough review and generous comments on the work under review. Our response is in the attached document / supplement.
  
  Citation: https://doi.org/10.5194/egusphere-2022-1152-AC1
RC2:
'Comment on egusphere-2022-1152', Anonymous Referee #2, 14 Dec 2022

The article discussed the possible hazard caused by volcano eruption at Lake Taupo through numerical simulations. The manuscript is well-written, and the overall quality of the scientific approach is very good, and the results are presented clearly. As a numerical modeler, I have several comments on this article.

First, what is the motivation of using MLSW for this case?

Authors used multi-layer SWE for numerical simulations, and showed that the computational efficiency of the MLSW solver is poor. Recent studies on HTHH (Pakoksun et al. 2022) showed that numerical simulations with NLSW or Boussinesq-type equations are in good agreement with the near-field observation. It would be desirable for the authors to perform numerical simulations with NLSW and Boussinesq eqs, and compare the maximum inundation heights and computational efficiency.

Authors may argue that the waves generated by the volcanic eruption are very dispersive. However, I am not sure that Lake Taupo is large enough for the waves to develop dispersive characteristics since the waves reach the shore in 12 minutes or less.

Second, it is unclear how the authors handled the bottom friction.

One of the important factors for predicting current speed and Inundation height is the bottom friction coefficient. In the study of inundation by storm surges, Madli and Dawson (2014) claimed that variable friction can be important to take into account in storm surge simulations, and numerical models apply complex bottom friction fields. I would like to know how authors considered the bottom friction, and wet/dry interface for simulations.

Comments

L59-63 unclear

L77 “have been attempted using” -> used

Section 2.1 & 2.2 are similar to Hayward et al. 2022b. Authors may shorten these parts.

Inconsistent in the simulated time. Which is correct? L204 “All runs were executed for 24 minutes of simulated time” L219 “All scenarios were computed until a simulated time of 1400 s”

L211 “The number of these placed along a section is set to match the maximum horizontal grid resolution.” unclear

L240 “leading to a longer duration from the first arrival to the maximum amplitude wave at greater distances from the source” unclear

L309 “For these sources, this 310 does not happen and therefore the reflection, and the incidence of troughs between wave peaks at the inlet, produces negative (or reverse) discharges from the inlet towards the lake.” unclear

L312 72% of something?

Figure 4 & 5 The information on the wave height is hard to discern from the figures. It would be better to plot the maximum wave height along the perimeter of the lake (in 2-D)

References

Pakoksung, Kwanchai, Anawat Suppasri, and Fumihiko Imamura. "The near-field tsunami generated by the 15 January 2022 eruption of the Hunga Tonga-Hunga Ha’apai volcano and its impact on Tongatapu, Tonga." 12.1 (2022): 1-15.

Mandli, Kyle T., and Clint N. Dawson. "Adaptive mesh refinement for storm surge." 75 (2014): 36-50.

Citation: https://doi.org/10.5194/egusphere-2022-1152-RC2
- AC2: 'Reply on RC2', Matthew Hayward, 08 Jan 2023
  
  Thank you for the helpful and insightful review and the effort taken over the comments. Our response is in the attached document / supplement.
  
  Citation: https://doi.org/10.5194/egusphere-2022-1152-AC2

Interactive discussion

Status: closed

RC1:
'Comment on egusphere-2022-1152', Raphael Paris, 09 Dec 2022

Dear colleagues,

Congratulations for this original work on tsunami hazard from underwater explosions in Lake Taupo. The manuscript represents a substantial contribution to hazard assesment in the Laker Taupo area, and the recent tsunami observed on November 30, 2022 shows how this study is timely (although the source of the tsunami is not yet elucidated). The overall quality of the scientific and methodological approach is very good, and the results are presented in a clear and concise way. The manuscript is well-written and well-illustrated. I have only some minor comments that are listed below:

- The use of the Volcanic Explosivity Index VEI (e.g. at lines 8-9, 264, 266, 320-324, 371) is not fully appropriate here, because it is based only on the total erupted volume. Here you are dealing with short discrete explosions. You can keep the VEI values in Table 1, for information. However I would recommend refering to the energy or mass eruption rates in the text.

- I suggest adding a sentence on the tectonic setting, which is too briefly mentioned at line 24: "The 350 km-long TVZ forms the southern part of a back-arc basin behind the Tonga-Kermadec subduction zone. Most of the eruptive centres are aligned along the intra-arc Taupo rift."

- The TVZ cannot be considered as a purely silicic magmatic system, because the magmas erupted range from andesites to rhyolites.

- At line 169 you state that "the water and explosion depths are equivalent." This is valid only in the case of a new vent formed on flat lake bottom. If the explosion comes from a pre-existing submarine cone, then explosion depth and water depth are not equivalents. This statement should therefore be qualified.

- I suggest slightly reorganizing the paper as follows:

3. Results

3.1. Tsunami propagation and wave heights at the coast

3.2. Tsunami inundation and potential impact on infratructures

- In the discussion, we lack a short discussion (at line 324?) on the probabilities mentioned on Table 1, in the light of the new results obtained here.

- The scenario presented at lines 330-334 is very pragmatic and relevant. Indeed in the case of such an eruption, all people would be probably evacuated before a tsunami could impact them. Could you add a short sentence of the consequences in terms of evacuation policies?

- In the conclusion, maybe you could formulate some recommendations for a local tsunami warning system in Lake Taupo?

Line-by-line corrections/suggestions:

-line 7: "This minimum size corresponds to..."

-9: "slope" rather than "runu-up"?

-37: "(HTHH, January 2022)"

-37: You could cite other references on the HTHH tsunami (e.g. Carvajal et al., 2002; Omira et al., 2022).

-41: cite also Maeno & Imamura (2011) for the 1883 Krakatau tsunami debate.

-56: "and in the present work it is..."

-59-63: this sentence is perhaps too long.

-93 refined against what?

-127: 7% of what?

-128-129: not sure if these lines are really useful here.

-194: "'We selected the five..."

-230-232: again, this sentence is a little bit long.

-248: refer to Fig 6 (c-d).

-257: 1 m instead of 1 cm?

-340: "asteroid-ocean impact, megatsunami from ocean-island flank collapse, and the recent tsunami..."

-345: "over half a decade old " -> a reference is needed here.

- 350-352: I was quite surprised that you don't cite here the work of Shen et al. (2021a, 2021b).

-367: "of what any tsunami hazard"

-384: What is the address of the website where Basilisk can be found?

Suggestions on figures:

- change the color of the eruption sites on figure 2 (white or yellow?).

- Fig 4: add source location (eruption sites) on the maps (same for fig 5).

- On the left side of fig 4, could you please mention the location number (1 to 5) as in Table 2, for more clarity (same for fig 5).

- Fig 6: add a dotted line to indicate the 1 m threshold value?

- Fig 6 c) d): in the caption, indicate that the first numbers refer to the scenarios and the name to the source locations 1 to 5.

- Fig. 8: Indicate to which closest source it corresponds on the different maps (same for fig 9).

Suggested references:

- Carvajal, M., Sepúlveda, I., Gubler, A., Garreaud, R.(2022). Worldwide Signature of the 2022 Tonga Volcanic Tsunami. Geophysical Research Letters, 49 (6), art. no. e2022GL098153; doi: 10.1029/2022GL098153

- Omira, R., Ramalho, R.S., Kim, J., González, P.J., Kadri, U., Miranda, J.M., Carrilho, F., & Baptista, M.A. (2022). Global Tonga tsunami explained by a fast-moving atmospheric source. Nature; doi.org/10.1038/s41586-022-04926-4

- Maeno, F., Imamura, F., 2011. Tsunami generation by a rapid entrance of pyroclastic flow into the sea during the 1883 Krakatau eruption, Indonesia. Journal of Geophysical Research 116, B09205

- Shen, Y., Whittaker, C.N., Lane, E.M., White, J.D.L., Power, W., Nomikou, P., 2021a. Laboratory experiments on tsunamigenic discrete subaqueous volcanic eruptions. Part 1: Free surface disturbances. Journal of Geophysical Research, Oceans 126, e2020JC016588.

- Shen, Y., Whittaker, C.N., Lane, E.M., White, J.D.L., Power, W., Nomikou, P., 2021a. Laboratory experiments on tsunamigenic discrete subaqueous volcanic eruptions. Part 2: Properties of generated waves. Journal of Geophysical Research, Oceans 126, e2020JC016587.

Citation: https://doi.org/10.5194/egusphere-2022-1152-RC1
- AC1: 'Reply on RC1', Matthew Hayward, 08 Jan 2023
  
  Thank you for the very thorough review and generous comments on the work under review. Our response is in the attached document / supplement.
  
  Citation: https://doi.org/10.5194/egusphere-2022-1152-AC1
RC2:
'Comment on egusphere-2022-1152', Anonymous Referee #2, 14 Dec 2022

The article discussed the possible hazard caused by volcano eruption at Lake Taupo through numerical simulations. The manuscript is well-written, and the overall quality of the scientific approach is very good, and the results are presented clearly. As a numerical modeler, I have several comments on this article.

First, what is the motivation of using MLSW for this case?

Authors used multi-layer SWE for numerical simulations, and showed that the computational efficiency of the MLSW solver is poor. Recent studies on HTHH (Pakoksun et al. 2022) showed that numerical simulations with NLSW or Boussinesq-type equations are in good agreement with the near-field observation. It would be desirable for the authors to perform numerical simulations with NLSW and Boussinesq eqs, and compare the maximum inundation heights and computational efficiency.

Authors may argue that the waves generated by the volcanic eruption are very dispersive. However, I am not sure that Lake Taupo is large enough for the waves to develop dispersive characteristics since the waves reach the shore in 12 minutes or less.

Second, it is unclear how the authors handled the bottom friction.

One of the important factors for predicting current speed and Inundation height is the bottom friction coefficient. In the study of inundation by storm surges, Madli and Dawson (2014) claimed that variable friction can be important to take into account in storm surge simulations, and numerical models apply complex bottom friction fields. I would like to know how authors considered the bottom friction, and wet/dry interface for simulations.

Comments

L59-63 unclear

L77 “have been attempted using” -> used

Section 2.1 & 2.2 are similar to Hayward et al. 2022b. Authors may shorten these parts.

Inconsistent in the simulated time. Which is correct? L204 “All runs were executed for 24 minutes of simulated time” L219 “All scenarios were computed until a simulated time of 1400 s”

L211 “The number of these placed along a section is set to match the maximum horizontal grid resolution.” unclear

L240 “leading to a longer duration from the first arrival to the maximum amplitude wave at greater distances from the source” unclear

L309 “For these sources, this 310 does not happen and therefore the reflection, and the incidence of troughs between wave peaks at the inlet, produces negative (or reverse) discharges from the inlet towards the lake.” unclear

L312 72% of something?

Figure 4 & 5 The information on the wave height is hard to discern from the figures. It would be better to plot the maximum wave height along the perimeter of the lake (in 2-D)

References

Pakoksung, Kwanchai, Anawat Suppasri, and Fumihiko Imamura. "The near-field tsunami generated by the 15 January 2022 eruption of the Hunga Tonga-Hunga Ha’apai volcano and its impact on Tongatapu, Tonga." 12.1 (2022): 1-15.

Mandli, Kyle T., and Clint N. Dawson. "Adaptive mesh refinement for storm surge." 75 (2014): 36-50.

Citation: https://doi.org/10.5194/egusphere-2022-1152-RC2
- AC2: 'Reply on RC2', Matthew Hayward, 08 Jan 2023
  
  Thank you for the helpful and insightful review and the effort taken over the comments. Our response is in the attached document / supplement.
  
  Citation: https://doi.org/10.5194/egusphere-2022-1152-AC2

Peer review completion

AR: Author's response | RR: Referee report | ED: Editor decision | EF: Editorial file upload

ED: Publish subject to minor revisions (review by editor) (28 Jan 2023) by Rachid Omira

AR by Matthew Hayward on behalf of the Authors (30 Jan 2023) Author's response Author's tracked changes Manuscript

ED: Publish subject to technical corrections (19 Feb 2023) by Rachid Omira

AR by Matthew Hayward on behalf of the Authors (19 Feb 2023) Manuscript

Journal article(s) based on this preprint

03 Mar 2023

Scenario-based modelling of waves generated by sublacustrine explosive eruptions at Lake Taupō, New Zealand

Matthew W. Hayward, Emily M. Lane, Colin N. Whittaker, Graham S. Leonard, and William L. Power

Nat. Hazards Earth Syst. Sci., 23, 955–971, https://doi.org/10.5194/nhess-23-955-2023,https://doi.org/10.5194/nhess-23-955-2023, 2023

Short summary

Matthew W. Hayward, Emily M. Lane, Colin N. Whittaker, Graham S. Leonard, and William L. Power

Viewed

Total article views: 310 (including HTML, PDF, and XML)

HTML	PDF	XML	Total	BibTeX	EndNote
229	69	12	310	1	2

HTML: 229
PDF: 69
XML: 12
Total: 310
BibTeX: 1
EndNote: 2

Views and downloads (calculated since 21 Nov 2022)

Month	HTML	PDF	XML	Total
Nov 2022	63	21	3	87
Dec 2022	81	24	5	110
Jan 2023	53	13	4	70
Feb 2023	30	10	0	40
Mar 2023	2	1	0	3
Apr 2023	0
May 2023	0
Jun 2023	0
Jul 2023	0
Aug 2023	0
Sep 2023	0
Oct 2023	0
Nov 2023	0
Dec 2023	0
Jan 2024	0
Feb 2024	0
Mar 2024	0
Apr 2024	0
May 2024	0
Jun 2024	0
Jul 2024	0
Aug 2024	0
Sep 2024	0

Cumulative views and downloads (calculated since 21 Nov 2022)

Month	HTML	PDF	XML	Total
Nov 2022	63	21	3	87
Dec 2022	81	24	5	110
Jan 2023	53	13	4	70
Feb 2023	30	10	0	40
Mar 2023	2	1	0	3
Apr 2023	0
May 2023	0
Jun 2023	0
Jul 2023	0
Aug 2023	0
Sep 2023	0
Oct 2023	0
Nov 2023	0
Dec 2023	0
Jan 2024	0
Feb 2024	0
Mar 2024	0
Apr 2024	0
May 2024	0
Jun 2024	0
Jul 2024	0
Aug 2024	0
Sep 2024	0

Viewed (geographical distribution)

Total article views: 304 (including HTML, PDF, and XML) Thereof 304 with geography defined and 0 with unknown origin.

Country	#	Views	%

Latest update: 06 Sep 2024

Download

The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.

Preprint (51179 KB)
Metadata XML

Short summary

In this paper, 20 explosive volcanic eruption scenarios of differing location and magnitude are simulated to investigate tsunami generation in Lake Taupō, New Zealand. A non-hydrostatic multilayer numerical scheme is used to resolve the highly dispersive generated wavefield. Inundation, hydrographic and related hazard outputs are produced, indicating significant inundation begins around the lake shore above 5 on the Volcanic Explosivity Index.


Total:	0
HTML:	0
PDF:	0
XML:	0