the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
SEATANI: hazards from seamounts in SouthEast Asia, Taiwan, and Andaman and Nicobar Islands (eastern India)
Abstract. Submarine volcanism makes up approximately 85 % of volcanism taking place on Earth, and its eruptions can be particularly hazardous, with the potential to cause large-scale sector collapse of the volcanic edifice, tsunamis, and ash dispersal. Recent examples include the eruptions in Japan and in the Kingdom of Tonga in 2021 and 2022 respectively, but there has been little to no study of submarine volcanoes in Southeast Asia and its surroundings. Here we provide a compilation of 466 seamounts from the region, from different published sources, through the SEATANI dataset (Southeast Asia + Taiwan + Andaman & Nicobar Islands). We use this newly compiled dataset to assess on a regional basis the seamount hazard potential and exposure potential as a springboard for future more quantitative hazard studies for the region. The hazard potential was assessed through seamount morphological/structural analyses, to determine the seamount evolution stage and, grade of maturity. The exposure potential was evaluated through two different approaches: An areal analysis of the number of assets within a 100 km radius of each seamount; and the development of a hazard-weighted seamount density map to highlight potential areas of interest for future more-in-depth studies. Our results show that there are several potentially hazardous seamounts in this region, and Taiwan had the highest hazard and exposure potential, for all assets considered, while Philippines, Indonesia and Vietnam have relatively high exposure potential for submarine communication cables and ship traffic density. The results from this work serve as a first step for southeast Asian and neighbouring countries to become more resilient against and prepared for submarine volcanic eruptions in the region.
-
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
(0 KB)
-
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
Interactive discussion
Status: closed
-
RC1: 'Comment on egusphere-2023-2153', Edgardo Cañón-Tapia, 23 Nov 2023
In this paper the authors focus attention on the hazard potential of submerged volcanoes in Southeast Asia. Hazard estimates are based on the location and characteristics of 466 seamounts. As the authors point out, seamount hazard estimation has been largely overlooked, and despite the limitations acknowledged by the authors, this paper aims at providing fundamental elements for future studies. With that in mind, the submitted paper has the potential of being a good contribution to the understanding of natural hazards, although there are at least three aspects that deserve more elaboration, one concerning the database used for the analysis and two concerning part of the method of analysis, as described next.
- Database used for the analysis.
The base of the hazard estimation reported in the current submission is the dataset of Gevorgian et al. (2023), which is an updated version of the dataset of Kim and Wessel (2011). As noted recently by Cañón-Tapia (2023, Geoscience Frontiers), there are other databases that can be used for analysis like this one. In particular, the database of Yesson et al (2011) would have been useful in the present context to bracket hazard estimates. When the knolls and seamounts databases of Yesson are examined, there are 13696 features within the area of the current study, which is a much larger number than the 466 seamounts analyzed in the current submission. The results of the hazard analysis based on the location of such a larger number of possible seamounts would therefore be very different from the results reported in the work by Verolino et al.
As elements for the discussion on this subject, it can be noted that in the GeoscFront2023 paper, it was shown that both the Kim-Wessel, (and consequently the extension by Gevorgian, hereafter GKW) and Yesson et al (hereafter Yetal) datasets are based on the same gravity signal. The difference in number of the reported features is related with the filters applied to that signal. Perhaps the most equilibrated assessment of both datasets is something like this: GKW over-filtered the signal, leaving only a very small proportion of seamounts, whereas Yetal under-filtered the signal allowing the introduction of some noise. Within that context, it is extremely important to note that one of the imposed filters in GKW was precisely to eliminate signals in continental margins and their vicinity. Therefore, that database has BY DESIGN, only a small fraction of seamounts in the region studied by Vitalino et al. This aspect of the dataset used for the current analysis presented by Vitalino et al. therefore leads to an underestimation of the hazard potential in the area of study. It can be argued that the Yetal dataset would lead to an overestimation of that hazard potential because it includes many places that might not be true volcanoes. Ideally, the bracketing of hazard estimation using the Yetal dataset should be included in the revised version of the work by Vitalino et al., but if this is not practical for logistic/financial support reasons, I think that at the very least the work by Vitalino et al should mention the possibility of finding a different hazard estimate using a different database, leaving the door open for a future study that is made using the Yetal dataset and that offers an upper bracket of current hazard estimation in this region.
- Method of analysis
2a) KDE implementation
Although the current version of the paper by Vitalinom et al. states that details of the weighting process are given in the methods section, there is no factual information about that process that can be consulted by the reader, other than a statement indicating that KDE was performed on ESRI ArcMap 10.7.1 that uses a default bandwidth (no mention of how the weighting was managed within the calculations made by KDE). Although such ommision can be corrected very easily, a main subject of discussion should be the reliability of the default bandwidth on the ArcMap software. As discussed at length by Cañón-Tapia (2020 Geomorphology, 2021 ESR and especially 2022 FEART and the November 2023 issue of the BSGM), method selection may introduce unsuspected biases on the results of an analysis of spatial distribution. In particular, the automatized selection of a bandwidth parameter may not be adequate to gain a complete picture of the characteristics of the distribution. In the current case under discussion, the bandwidth selected by the software clearly includes an influence of the large areas in which there are not observations. Some of those areas are the result of a too restrictive database (as mentioned above), but other areas without data are the natural result of the presence of emerged lands. Thus, ArcMap calculates a blanket bandwidth that is too large for the small population of volcanoes, therefore biasing the results towards a too large scale distribution. As a result, the hazards associated with the groups of volcanoes on the southwest (actually all along the east Java and Sumatra coasts) are entirely sub-estimated, whereas those on the north and southeast are somewhat overestimated. Also, it must be noted that the automatically selected bandwidth (>600 km) is in contrast with the interest mentioned by the authors on the 100 km limit. This discrepancy on the scale of different stages of the analysis should be discussed in a revised version of the current submission.
To address these issues, it is recommended that the results of an exploration of the KDE are reported, in which the bandwidth is manually selected within the 50 to 600 km (increments of 50 km should suffice), to identify the scales at which hazards (concentrations of volcanoes and proximity with human-made infrastructure) may be changing in different sections of the area of study. Given the small number of volcanoes used in the currently used dataset, the general trends may not change much from those already reported, except from the possibility of increasing hazard level in the SW. Anyway, as mentioned above, this would be a lower bracket of hazard estimation, and this aspect of the analysis should be clearly highlighted in the revised version. The results of a similar analysis completed with a much larger population of seamounts is likely to produce entirely different results, which correspond to an upper bracket on the hazard estimation.
2b) Elimination of large emerged volcanoes.
Although this criteria eliminated only 16 volcanoes, it is strange that those volcanoes were eliminated on a work focusing on hazard estimation. Large eruptions from emerged volcanoes can lead to the entrance of pyroclastic deposits to the sea. Recent examples from Montserrat come to mind. Also, the potential for tsunamis generated by either pyroclastic deposits entering shallow waters and by (perhaps less likely) sector collapses of large edifices should be considered. The historic eruption of Krakatoa and some documented tsunamis on the islands of Hawaii come to mind. Indeed, this type of danger is different from that forming the bulk of the reported work, but I think it deserves at least some mention somewhere in the text.
In summary, although the work under discussion addresses a relevant scientific question within the scope of NHESS, the analysis reported is not thorough enough to be considered an accurate estimation of potential hazard. Technically, the paper is well written and concise, and figures are adequate. Consequently, once the above mentioned issues have been addressed, the paper will surely be a good contribution for this journal.
-
AC1: 'Reply on RC1', Andrea Verolino, 24 Nov 2023
We deeply thank the reviewer and appreciate the comments provided. The assessment of our manuscript highlights the reviewer’s knowledge in the field and will surely help to improve it substantially.
In response to the above-mentioned analysis, as an anticipation of what will be done in the revised version of our manuscript, we comment below:
Choice of dataset: We were aware of the other datasets mentioned in Cañón-Tapia (2023), however, we decided to use the one proposed by Gevorgian et al. (2023) for several reasons: (1) it’s the most up-to-date in terms of data quality of the Vertical Gravity Gradient (data noise reduction of 40% from previous VGG datasets) used for the seamount detection; (2) it only focuses on volcanic seamounts, while the other datasets may also include non-volcanic features, biasing our results, and also small features (knolls) that will be little relevant in terms of hazard potential for our study; (3) the Gevorgian et al. dataset might have been over-filtered, removing possible seamounts from areas such as continental margins, however, at this stage we prefer to maintain a conservative approach in terms of seamount detection, still highlighting the fact that higher resolution bathymetry is needed in areas such as the Sunda Shelf. Nevertheless, the above points will be properly reported in the respective sections of the manuscript, clarifying aspects that were not clear so far, and eventually add relevant discussions.
KDE implementation: We will conduct sensitivity tests with different bandwidths, to assess whether there are substantial changes in the resulting Kernel map, and will adjust the text accordingly. We will also look into the papers/reports suggested by the reviewer in order to make sure that our analysis is not biased, or we will add relevant discussions if needed.
Elimination of large emerged volcanoes: The threshold at which a volcanic island is still considered a seamount has not been previously defined. Here we provide guidelines of this possible threshold, with the intent of maintaining the focus on the unknown hazard potential of seamounts, which is indeed referring to their submerged parts. Large volcanic islands with more than 30% of their volume above the water level were deliberately considered out of the scope of this work, with many of them already included in previous studies (as highlighted already in the manuscript). Despite this, we still included seamounts partially emerged, in order to include those with similar characters to the Hunga Tonga Hunga-Ha’apai volcano (e.g. Krakatau), which recently showed the devastating potential of these types (stage3/stage4) of seamounts.
Citation: https://doi.org/10.5194/egusphere-2023-2153-AC1
-
RC2: 'Comment on egusphere-2023-2153', Eugenio Nicotra, 05 Dec 2023
The manuscript "SEATANI: hazards from seamounts in SouthEast Asia, Taiwan, and Andaman and Nicobar Islands (eastern India)" by Verolino et al. is an interesting paper about the potential hazard from seamounts in the SouthEast Asia.
The manuscript is well written, maybe too long and dispersive in some sections (see annotated pdf). Altough I'm not a mother tongue, I think that authors abused of bullet points all over the text, especially in the Introduction. A language refining would be appreciated.
Some minor points of discussion:
- Within the volcanic seamounts classification, which is a qualitative analysis of the shape of the seamount, I think that the importance of the ridge/fissural/linear shape is a little bit underestimated. South China Sea and Banda Sea are two geodynamic settings characterized by continental rifting and back-arc extension, two geodynamic settings which preferably generate volcanic ridge. A typical example has been recently studied in the Tyrrhenian Sea (back-arc setting), with the new bathymetric model for Marsili seamount (Nicotra et al., 2023). So, if data available to the authors are defined enough, it could be useful also to evaluate the elongation of the seamounts.
- On my opinion, section 5.3.1 about the geodynamic context is not useful at that point of the manuscript. So: 1) or it is moved in a background section after Introduction; 2) or it becomes a point of discussion in and background in the 5.1 section
- Authors well know the limitation of their work (and also dedicate the 5.4 section to this), mainly due to the association of geological objects having millions of years of difference. Although results are very interesting, maybe the introduction is a bit over overloaded about the importance of this paper in terms of hazard from seamounts. A lot of work is still needed, but this culd represent a first step in the costruction of a useful hazard database for seamounts of the South-East sea.
-
AC2: 'Reply on RC2', Andrea Verolino, 08 Dec 2023
We thank Dr E. Nicotra for the suggestions provided. We believe that they are of value, and we will adjust the text accordingly by applying the following:
- We will expand the text in the introduction and/or discussions on volcanic ridges, in reference to their geodynamic context.
- Regarding section 5.3.1 (“Geodynamic context”), we believe that it is more appropriate to have it in the discussion because it refers to some of the results obtained in this work. However, we will consider whether to keep it as 5.3.1, or move it to a subsection of 5.1, in the context of potential sources of interest.
- We will lighten the introduction, where suggested by the reviewer.
Citation: https://doi.org/10.5194/egusphere-2023-2153-AC2
-
AC2: 'Reply on RC2', Andrea Verolino, 08 Dec 2023
Interactive discussion
Status: closed
-
RC1: 'Comment on egusphere-2023-2153', Edgardo Cañón-Tapia, 23 Nov 2023
In this paper the authors focus attention on the hazard potential of submerged volcanoes in Southeast Asia. Hazard estimates are based on the location and characteristics of 466 seamounts. As the authors point out, seamount hazard estimation has been largely overlooked, and despite the limitations acknowledged by the authors, this paper aims at providing fundamental elements for future studies. With that in mind, the submitted paper has the potential of being a good contribution to the understanding of natural hazards, although there are at least three aspects that deserve more elaboration, one concerning the database used for the analysis and two concerning part of the method of analysis, as described next.
- Database used for the analysis.
The base of the hazard estimation reported in the current submission is the dataset of Gevorgian et al. (2023), which is an updated version of the dataset of Kim and Wessel (2011). As noted recently by Cañón-Tapia (2023, Geoscience Frontiers), there are other databases that can be used for analysis like this one. In particular, the database of Yesson et al (2011) would have been useful in the present context to bracket hazard estimates. When the knolls and seamounts databases of Yesson are examined, there are 13696 features within the area of the current study, which is a much larger number than the 466 seamounts analyzed in the current submission. The results of the hazard analysis based on the location of such a larger number of possible seamounts would therefore be very different from the results reported in the work by Verolino et al.
As elements for the discussion on this subject, it can be noted that in the GeoscFront2023 paper, it was shown that both the Kim-Wessel, (and consequently the extension by Gevorgian, hereafter GKW) and Yesson et al (hereafter Yetal) datasets are based on the same gravity signal. The difference in number of the reported features is related with the filters applied to that signal. Perhaps the most equilibrated assessment of both datasets is something like this: GKW over-filtered the signal, leaving only a very small proportion of seamounts, whereas Yetal under-filtered the signal allowing the introduction of some noise. Within that context, it is extremely important to note that one of the imposed filters in GKW was precisely to eliminate signals in continental margins and their vicinity. Therefore, that database has BY DESIGN, only a small fraction of seamounts in the region studied by Vitalino et al. This aspect of the dataset used for the current analysis presented by Vitalino et al. therefore leads to an underestimation of the hazard potential in the area of study. It can be argued that the Yetal dataset would lead to an overestimation of that hazard potential because it includes many places that might not be true volcanoes. Ideally, the bracketing of hazard estimation using the Yetal dataset should be included in the revised version of the work by Vitalino et al., but if this is not practical for logistic/financial support reasons, I think that at the very least the work by Vitalino et al should mention the possibility of finding a different hazard estimate using a different database, leaving the door open for a future study that is made using the Yetal dataset and that offers an upper bracket of current hazard estimation in this region.
- Method of analysis
2a) KDE implementation
Although the current version of the paper by Vitalinom et al. states that details of the weighting process are given in the methods section, there is no factual information about that process that can be consulted by the reader, other than a statement indicating that KDE was performed on ESRI ArcMap 10.7.1 that uses a default bandwidth (no mention of how the weighting was managed within the calculations made by KDE). Although such ommision can be corrected very easily, a main subject of discussion should be the reliability of the default bandwidth on the ArcMap software. As discussed at length by Cañón-Tapia (2020 Geomorphology, 2021 ESR and especially 2022 FEART and the November 2023 issue of the BSGM), method selection may introduce unsuspected biases on the results of an analysis of spatial distribution. In particular, the automatized selection of a bandwidth parameter may not be adequate to gain a complete picture of the characteristics of the distribution. In the current case under discussion, the bandwidth selected by the software clearly includes an influence of the large areas in which there are not observations. Some of those areas are the result of a too restrictive database (as mentioned above), but other areas without data are the natural result of the presence of emerged lands. Thus, ArcMap calculates a blanket bandwidth that is too large for the small population of volcanoes, therefore biasing the results towards a too large scale distribution. As a result, the hazards associated with the groups of volcanoes on the southwest (actually all along the east Java and Sumatra coasts) are entirely sub-estimated, whereas those on the north and southeast are somewhat overestimated. Also, it must be noted that the automatically selected bandwidth (>600 km) is in contrast with the interest mentioned by the authors on the 100 km limit. This discrepancy on the scale of different stages of the analysis should be discussed in a revised version of the current submission.
To address these issues, it is recommended that the results of an exploration of the KDE are reported, in which the bandwidth is manually selected within the 50 to 600 km (increments of 50 km should suffice), to identify the scales at which hazards (concentrations of volcanoes and proximity with human-made infrastructure) may be changing in different sections of the area of study. Given the small number of volcanoes used in the currently used dataset, the general trends may not change much from those already reported, except from the possibility of increasing hazard level in the SW. Anyway, as mentioned above, this would be a lower bracket of hazard estimation, and this aspect of the analysis should be clearly highlighted in the revised version. The results of a similar analysis completed with a much larger population of seamounts is likely to produce entirely different results, which correspond to an upper bracket on the hazard estimation.
2b) Elimination of large emerged volcanoes.
Although this criteria eliminated only 16 volcanoes, it is strange that those volcanoes were eliminated on a work focusing on hazard estimation. Large eruptions from emerged volcanoes can lead to the entrance of pyroclastic deposits to the sea. Recent examples from Montserrat come to mind. Also, the potential for tsunamis generated by either pyroclastic deposits entering shallow waters and by (perhaps less likely) sector collapses of large edifices should be considered. The historic eruption of Krakatoa and some documented tsunamis on the islands of Hawaii come to mind. Indeed, this type of danger is different from that forming the bulk of the reported work, but I think it deserves at least some mention somewhere in the text.
In summary, although the work under discussion addresses a relevant scientific question within the scope of NHESS, the analysis reported is not thorough enough to be considered an accurate estimation of potential hazard. Technically, the paper is well written and concise, and figures are adequate. Consequently, once the above mentioned issues have been addressed, the paper will surely be a good contribution for this journal.
-
AC1: 'Reply on RC1', Andrea Verolino, 24 Nov 2023
We deeply thank the reviewer and appreciate the comments provided. The assessment of our manuscript highlights the reviewer’s knowledge in the field and will surely help to improve it substantially.
In response to the above-mentioned analysis, as an anticipation of what will be done in the revised version of our manuscript, we comment below:
Choice of dataset: We were aware of the other datasets mentioned in Cañón-Tapia (2023), however, we decided to use the one proposed by Gevorgian et al. (2023) for several reasons: (1) it’s the most up-to-date in terms of data quality of the Vertical Gravity Gradient (data noise reduction of 40% from previous VGG datasets) used for the seamount detection; (2) it only focuses on volcanic seamounts, while the other datasets may also include non-volcanic features, biasing our results, and also small features (knolls) that will be little relevant in terms of hazard potential for our study; (3) the Gevorgian et al. dataset might have been over-filtered, removing possible seamounts from areas such as continental margins, however, at this stage we prefer to maintain a conservative approach in terms of seamount detection, still highlighting the fact that higher resolution bathymetry is needed in areas such as the Sunda Shelf. Nevertheless, the above points will be properly reported in the respective sections of the manuscript, clarifying aspects that were not clear so far, and eventually add relevant discussions.
KDE implementation: We will conduct sensitivity tests with different bandwidths, to assess whether there are substantial changes in the resulting Kernel map, and will adjust the text accordingly. We will also look into the papers/reports suggested by the reviewer in order to make sure that our analysis is not biased, or we will add relevant discussions if needed.
Elimination of large emerged volcanoes: The threshold at which a volcanic island is still considered a seamount has not been previously defined. Here we provide guidelines of this possible threshold, with the intent of maintaining the focus on the unknown hazard potential of seamounts, which is indeed referring to their submerged parts. Large volcanic islands with more than 30% of their volume above the water level were deliberately considered out of the scope of this work, with many of them already included in previous studies (as highlighted already in the manuscript). Despite this, we still included seamounts partially emerged, in order to include those with similar characters to the Hunga Tonga Hunga-Ha’apai volcano (e.g. Krakatau), which recently showed the devastating potential of these types (stage3/stage4) of seamounts.
Citation: https://doi.org/10.5194/egusphere-2023-2153-AC1
-
RC2: 'Comment on egusphere-2023-2153', Eugenio Nicotra, 05 Dec 2023
The manuscript "SEATANI: hazards from seamounts in SouthEast Asia, Taiwan, and Andaman and Nicobar Islands (eastern India)" by Verolino et al. is an interesting paper about the potential hazard from seamounts in the SouthEast Asia.
The manuscript is well written, maybe too long and dispersive in some sections (see annotated pdf). Altough I'm not a mother tongue, I think that authors abused of bullet points all over the text, especially in the Introduction. A language refining would be appreciated.
Some minor points of discussion:
- Within the volcanic seamounts classification, which is a qualitative analysis of the shape of the seamount, I think that the importance of the ridge/fissural/linear shape is a little bit underestimated. South China Sea and Banda Sea are two geodynamic settings characterized by continental rifting and back-arc extension, two geodynamic settings which preferably generate volcanic ridge. A typical example has been recently studied in the Tyrrhenian Sea (back-arc setting), with the new bathymetric model for Marsili seamount (Nicotra et al., 2023). So, if data available to the authors are defined enough, it could be useful also to evaluate the elongation of the seamounts.
- On my opinion, section 5.3.1 about the geodynamic context is not useful at that point of the manuscript. So: 1) or it is moved in a background section after Introduction; 2) or it becomes a point of discussion in and background in the 5.1 section
- Authors well know the limitation of their work (and also dedicate the 5.4 section to this), mainly due to the association of geological objects having millions of years of difference. Although results are very interesting, maybe the introduction is a bit over overloaded about the importance of this paper in terms of hazard from seamounts. A lot of work is still needed, but this culd represent a first step in the costruction of a useful hazard database for seamounts of the South-East sea.
-
AC2: 'Reply on RC2', Andrea Verolino, 08 Dec 2023
We thank Dr E. Nicotra for the suggestions provided. We believe that they are of value, and we will adjust the text accordingly by applying the following:
- We will expand the text in the introduction and/or discussions on volcanic ridges, in reference to their geodynamic context.
- Regarding section 5.3.1 (“Geodynamic context”), we believe that it is more appropriate to have it in the discussion because it refers to some of the results obtained in this work. However, we will consider whether to keep it as 5.3.1, or move it to a subsection of 5.1, in the context of potential sources of interest.
- We will lighten the introduction, where suggested by the reviewer.
Citation: https://doi.org/10.5194/egusphere-2023-2153-AC2
-
AC2: 'Reply on RC2', Andrea Verolino, 08 Dec 2023
Peer review completion
Journal article(s) based on this preprint
Viewed
Since the preprint corresponding to this journal article was posted outside of Copernicus Publications, the preprint-related metrics are limited to HTML views.
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
200 | 0 | 0 | 200 | 0 | 0 |
- HTML: 200
- PDF: 0
- XML: 0
- Total: 200
- BibTeX: 0
- EndNote: 0
Viewed (geographical distribution)
Since the preprint corresponding to this journal article was posted outside of Copernicus Publications, the preprint-related metrics are limited to HTML views.
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1
Andrea Verolino
Su Fen Wee
Susanna F. Jenkins
Fidel Costa
Adam D. Switzer
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.