Abstract
The massive eruption of the Hunga Tonga-Hunga Haâapai (HTHH) volcano in Tonga on 15 January 2022 at 04:15 UTC had a global impact and triggered an atmospheric wave and a tsunami. We first analyzed observation data from meteorological stations and tide gauges at 12 locations. Low-frequency trends in the observation data were removed by using a high-pass filter. Fourier and wavelet spectral analyses were applied to determine the frequency characteristics of the filtered data. Modal analysis was developed and used to investigate natural oscillation periods based on bathymetry. The results showed that the Lamb wave generated by the atmospheric pressure wave arrivedâ~â7 andâ~â44Â h after the eruption. The tsunami arrivedâ~â11 andâ~â45Â h after the eruption, which corresponded to the arrival time of the Lamb wave. The dominant periods of the Lamb waves wereâ~â7.7 andâ~â7.5Â min, and for the tsunamis they wereâ~â9.9 andâ~â28.7Â min. The periods derived from the spectral analysis matched the natural oscillation of the eigenperiod derived from the modal analysis, in eight out of the twelve stations. This study provides valuable insight and information regarding the nonseismic and far-field effects of tsunamis generated by volcanic eruptions.
Similar content being viewed by others
Introduction
The 2022 HTHH eruption produced an atmospheric pressure wave similar to a Lamb wave, which is characterized by a pressure pulse. Lamb waves transfer energy over a long distance to generate tsunami waves1. The nondispersive Lamb wave propagated along the horizontal direction at a speed of approximately 315Â m/s2,3. The tsunami generated by the Lamb wave was induced by the massive volcanic eruption, which was similar to the tsunami that occurred following the 1883 Krakatau eruption4. Similarly, the tsunami generated by meteorological pressure in Lake Michigan in 1954 was also reported and modeled5.
The tsunami was produced by the atmospheric pressure wave released from the 2022 HTHH eruption and propagated across the Pacific Ocean6,7. The maximum measured tsunami height was approximately 2.5 m in Central America8 and higher than 1 m in Chile, New Caledonia, Vanuatu, USA, and Japan9,10,11,12,13,14. The tsunami caused damage along the coastal area of the Pacific Ocean12,15,16,17,18. In Japan, the damage caused by this tsunami was located along the coastal areas on the eastern side of Japan and extended from the north to the south but was concentrated in the bays and ports. Damage was reported on marine vessels and aquaculture rafts12, as shown in Fig. 1, which indicates the location of damage at the prefecture scale. At present, the cause of the damage remains unclear, and may be solely related to the stiffness of the structure or perhaps to the physical tsunami mechanisms, such as resonance.
Study area. Observed station locations (atmospheric pressure and sea level) are presented by a blue triangle. These are arranged from north to south along the eastern coastline of Japan. The reported damage on the prefecture scale was obtained from Imamura et al.12 The figure was generated using Python version 3.866 with the Matplotlib library67.
Tsunami resonance occurs when the tsunami wave period is similar to the natural oscillation of the water surface19 and is produced by the reflection and interference of tsunami waves at the edge of a bay or port. The amplification of the natural oscillation of the bay can generate a higher water level in the bay and increase run-up along the coast19,20,21,22. The energy generated by natural high-frequency tsunami oscillations can generate strong currents inside bays that may lead to damage on boats, ships, port infrastructure, and aquaculture rafts23. Therefore, understanding resonance characteristics is important to help identify tsunami hazards for coastal planning and emergency management24,25,26. Some previous studies have been conducted to investigate tsunami resonance around the Pacific Ocean. Cortés et al.27 calculated the longwave resonance to investigate the regions that had an increased potential for tsunamis, specifically within the northern coastal area of Chile. Abe28 investigated the dominant period of tsunami resonance in the bays of northern Japan and compared them with the spectra of tsunamis from the near-field and trans-Pacific location. Wang et al.19 investigated the tsunami dominant period from several trans-Pacific tsunami events in the Japan, Hokkaido, and Sanriku regions and compared them with the natural dominant period.
In this study, we first investigated the resonance characteristics of the 2022 HTHH trans-Pacific tsunami event in Japan. The investigation of resonance characteristics is based on the observations of atmospheric pressure and tsunami height in the bay along the eastern coast area of Japan (see Fig. 1) made by 12 stations, and some stations were in the areas damaged by the 2022 HTHH event. First, the collected data were quality controlled and processed by removing a low-frequency trend with a high-pass filter. We analyzed the filtered atmospheric pressure and tsunami height data from the 12 stations to investigate the potential relationship between the atmospheric pressure wave and tsunami wave. The frequency characteristics (dominant period) of atmospheric pressure and tsunami waves were calculated using both wavelet and Fourier spectral analyses. We also calculated the eigenmodes of the natural oscillation of the bay/port where the 12 stations were located to compare the dominant period from the HTHH 2022 event and the natural oscillation (on regional and local domains). This is the first study to investigate the resonance characteristics from a real a nonseismic tsunami event in Japan in which the tsunami was related to the atmospheric pressure waves generated by a volcanic eruption. We also provide theoretical information on the tsunami generated by the 2022 HTHH volcanic eruption.
Results
Data processing results
Atmospheric pressure data were collected from Weathernews Inc.29. We selected 12 meteorological stations along the eastern coast area of Japan that were located in relation to the selected tsunami gauge station (see Fig. 1), which had at a sampling rate of 1 min. The atmospheric pressure data were selected from January 14, 2022, at 00:00 UTC to January 19, 2022, at 00:00 UTC, which covers more 90 h following the HTHH eruption, as shown in Fig. S1a. The raw data were quality controlled to remove spikes, gaps, or redundant values. Next, a high-pass filter was applied to remove diurnal pressure oscillations30 with a cutoff of 6.67âÃâ10â4 Hz (1500 s), and the low-frequency trend is shown in Fig. S1a. The raw data subtracted from the low-frequency data become the filtered dataset (see Fig. S1b), which is summarized in Fig. S2. Tsunami height data were collected from 12 tide gauges along the eastern Japan coast (see Fig. 1). All gauge stations are in bays or ports, and their sampling rate is 1 min. Data from the 11 stations were obtained from the Intergovernmental Oceanographic Commission31: Hansaki, Kushiro, Hakodate, Ofunato, Mera, Omaezaki, Kushimoto, Tosashimizu, Aburatsu, Naha, and Ishigakijima. However, data from 1 station was obtained from the Geospatial Information Authority of Japan32, Soma. The selected time series of tsunami data covered the time between January 14, 2022, 00:00 UTC, and January 19, 2022, 00:00 UTC (see Fig. S3a). After the raw data were subjected to basic quality control, a high-pass filter was applied to remove the tidal signals with a cutoff of 6.25âÃâ10â3 Hz (160 s)33,34,35, and the low-frequency trend is shown in Fig. S3a. The dataset was filtered by calculating the differences between the raw data and the low-frequency data (see Fig. S3b) and are summarized in Fig. S4.
Waveform analysis results
Atmospheric pressure and tsunami height from the filtered low-frequency data derived from the 12 stations are presented in Fig. 2. From the atmospheric pressure, two persistent waves (Lamb waves) were visible, with intervals of approximately 36 h. The duration of the first wave was from January 15, 2022, 11:00 UTC to January 17, 2022, 00:00 UTC, and the second wave occurred from January 17, 2022, 00:00 UTC to January 18, 2022, 12:00 UTC. The tsunami waves were also classified to be consistent with the duration of the Lamb wave, and two tsunami waves also arrived after the Lamb wave. The arrival times of both the Lamb wave and tsunami wave were generally consistent.
Time series of atmospheric pressure (green line) and sea level (blue line) records. The eruption of the HTHH volcano on 04:14:45 UTC, 15 January 2022, is indicated by the dashed black line. The first (12:00 UTC, 15 January 2022 to 00:00 UTC, 17 January 2022) and second (00:00 UTC, 17 January 2022 to 12:00 UTC, 18 January 2022) wave components were observed over approximately 36Â h.
The data from the 12 stations indicate that the first Lamb wave arrived in Japan at approximately 11:00 (UTC) on 15 January 2022, approximately 7 h after the HTHH eruption, as shown in Fig. 2. The amplitude at the crest of the first Lamb wave was approximately 2 hPa with a wavelength of approximately 30 min. The shape of the first Lamb wave was represented by an N-wave shape and a small wave of approximately 0.5 hPa was visible following the crest.
The sea surface gradually varied and rose for approximately 2 h after the arrival of the first Lamb, as indicated by the tide gauge observations of the first tsunami wave. The height of the leading tsunami wave varied positively from 5 to 10 cm in the early phase and was apparent after 13:00 (UTC) on 15 January 2022. The maximum amplitude of the tsunami wave overall occurred in the later phase, as shown in Fig. 2. Hanasaki station, located at the north end of Hokkaido Island, reveals that the maximum height of the first tsunami was approximately 0.5 m and occurred at 15:30 UTC on 15 January 2022 (00:30 JST on 16 January 2022). Several peaks were continually generated over a period of 7 h until the largest peak of approximately 0.6 m occurred at 20:00 UTC on 15 January 2022 (5:00 JST on 16 January 2022). The Kushiro and Hakodate stations had a similar pattern to the Hanasaki station, with the first maximum peaks at also have heights approximately equal to 0.5 m and 0.3 m, respectively. Ofunato, Soma, and Mera stations, located on the eastern side of Japan, showed that the small first peak occurred at 15:00 UTC on 15 January 2022 (00:00 JST on 16 January 2022). Several peaks were continually generated until the largest peak was observed at 18:00 UTC on 15 January 2022 (03:00 JST on 16 January 2022). The largest tsunami peaks at Ofunato, Soma, and Mera were 0.25 m, 0.55 m, and 0.5 m, respectively. The stations on the western side of Japan are Omaezaki, Kushimoto, Tosashimizu, and Aburatsu, which revealed that the peaks on this coast were shorter in duration when compared to the peaks observed by the station in the northern area. The four stations recorded the largest peak at 15:00 UTC on 15 January 2022 (00:00 JST on 16 January 2022). The largest peaks at Omaesaki, Kushimoto, Tosashimizu, and Aburatsu were 0.75 m, 0.9 m, 1.0 m, and 0.55 m, respectively. Tosashimizu recorded the largest peak value among the 12 stations. Naha and Ishigakijima, located on Okinawa Island, which is southwestern Japan, recorded peals with a smaller amplitude than recorded in the northern area. Both stations observed that the peak began small, and the tsunami wave continued until the largest peak was recorded at 14:00 UTC on 15 January 2022 (23:00 JST). The largest peaks at Naha and Ishigakijima were 0.2 m and 0.1 m, respectively.
The second Lamb wave was measured by atmospheric pressure data and arrived at 00:00 UTC on 17 January 2022, as shown in Fig. 2. The arrival time of the second wave was approximately 44 h after the HTHH eruption and 36 h after the first Lamb wave arrived. The amplitude of the second Lamb was approximately 0.5 hPa, which was smaller than half of the first Lamb wave. Overall, as shown in Fig. 2, tsunami waves gradually increased after the second Lamb wave arrived at approximately 1â2 h. The tsunami height during the early phase was approximately 0.03â0.05 m (see Fig. S3b), and the largest peak occurred in the later phase. The maximum tsunami heights recorded by all stations varied from 0.05 to 0.1 m.
Wavelet spectral analysis results
Wavelet spectra were used to present the duration of the Lamb and tsunami waves generated by the 2022 HTHH eruption. Determining wave duration was not easily accomplished by analyzing the waveform11. Figure 3 presents the clear pattern of the duration of two waves (the first and second in the waveform analysis), which are from atmospheric pressure (1st column) and tsunami (2nd column) waves.
The energy signal of the first Lamb wave in the wavelet spectra analysis revealed that the arrival time occurred at 11:00 UTC on 15 January 2022. This result was similar to the findings from the waveform analysis across all stations for the duration of the first wave. The oscillation time of this wave, as recorded by six stations (Hanasaki, Kushiro, Hakodate, Ofunato, Soma, and Ishigakijima), was approximately 3Â h, and persisted until 14:00 UTC on 15 January 2022. The main energy was concentrated in the period from 8 to 128Â min. The oscillation time was approximately 4Â h until 15:00 UTC on 15 January 2022, with the main energy in the period from 4 to 128Â min, as recorded by the other six stations, Mera, Omaezaki, Kushimoto, Tosashimizu, Aburatsu, and Naha.
The wavelet analysis of the first tsunami wave reveals that the energy signal had an arrival time of at 11:00 UTC on 15 January 2022, which was similar to the results of the waveform analysis for all stations. The main energy was concentrated in the period from 8 to 128Â min for five stations, Hanasaki, Kushiro, Ofunato, Omaezaki, and Naha. The five stations revealed differences in oscillation times of 36Â h, 36Â h, 18Â h, 36Â h, and 36Â h at Hanasaki, Kushiro, Ofunato, Omaezaki, and Naha, respectively. Two stations (Hakodate and Soma) recorded the main energy in the period range of 8 to 256Â min with an oscillating time of approximately 36Â h. The main energy of the period ranged from 4 to 64Â min at Mera, Kushimoto, and Tosashimizu with oscillating times of 12Â h, 36Â h, and 24Â h, respectively. Aburatsu station had the main energy in the period of 8 to 64Â min with an oscillation time longer than 36Â h. The main energy in a period of 4 to 256Â min was recorded at Ishigakijima station with an oscillating time of 15Â h.
For the duration of the second wave, the energy signal of the Lamb wave was identified based on an arrival time at 00:00 UTC on 17 January 2022, which was similar to the results of the waveform analysis. The main energy in the period was in the range from 8 to 128Â min, and the energy was lower than the energy in the first wave. The oscillation time of this wave was approximately 2Â h. For this second tsunami wave, the energy signal was also lower than the energy signal of the first wave. Eight stations recorded wave energy that continued from the first wave: Hanasaki, Kushiro, Hakodate, Soma, Omaesaki, Kushimoto, Aburatsu, and Naha. The eight stations presented an oscillation time approximately longer than 12Â h following the arrival of the second wave. The other four stations provided a lower response signal after the arrival time of the second wave, which was discontinuous with the signal energy of the first wave, Ofunato, Mera, Tosashimizu, and Ishigakijima.
Fourier spectral analysis results
Fourier spectra are used to present the energy spectra related to different periods (or frequencies) in this study. The Fourier spectra of the durations of the first and second waves were plotted and compared with the Fourier spectra of the background waves for atmospheric pressure and tsunami height, as shown in Fig. 4. In Fig. 4, the first column shows the Fourier spectra of atmospheric pressure (Lamb wave), and the second column shows the Fourier spectra of the tsunami height. The gap between wave spectra (first and second waves) and background spectra is the wave energy spectra for the 2022 HTHH eruption. The shapes of the wave spectra are normally followed by the background spectra to reveal the effect of local air pressure and bathymetry on the Lamb and tsunami waves, respectively11,36,37.
Fourier spectra of the atmospheric pressure wave and the tsunami wave from the 12 stations. Red and blue lines are the spectra of the first and second waves, respectively. The green line is the spectra of the background signals. For the energy of the atmospheric pressure wave, the dominant periods of the first tsunami wave are marked by the vertical line (red dashed line). For the tsunami wave energy, the vertical red dashed line marks the dominant periods. For the natural oscillation within the area around the observation points, the vertical black dashed line represents the period of the eigenmode (regional domain), and the vertical black line is based on the fundamental mode period (local domain). The first column shows the spectra of the atmospheric pressure wave, and the second column shows the spectra of the tsunami wave.
The Fourier spectra of the Lamb waves were plotted by varying the periods (frequency) from 400 to 2 min (0.0025 to 0.5 Hz) and energies from 1âÃâ10â3 to 1âÃâ10â8 m2/min for both wave durations (first and second), as shown in the first column of Fig. 4. The long period (low frequency) energy spectra for both wave durations were greater than 200 min (0.005 Hz) and were close to the energy spectra of the background. The gap in the energy spectra between the Lamb wave spectra and background spectra gradually increased after periods below 200 min for both wave durations. The first wave duration had a gap higher than the second wave duration. Overall, the maximum gaps were located between periods of 5 and 90 min, and the average dominant periods during the gap were approximately 8 to 30 min.
The second column of Fig. 4 shows the Fourier spectra of the tsunami waves that were plotted by the same scale as the Lamb wave (shown in the first column). The energy spectra of both durations for periods greater than 200 min (0.005 Hz) were close to the energy spectra of the background spectra, which was similar to the spectral behavior of the Lamb wave. Additionally, the energy spectra were close to the short period (high frequency) background spectra, which was approximately Below 5 min (0.2 Hz). A period of maximum energy spectra greater than 30 min was observed at five stations in the east, Hansaki, Kushiro, Hakodate, Ofunato, and Soma. The seven stations (Mera, Omaezaki, Kushimoto, Tosashimizu, Aburatsu, Naha, and Ishigakijima) reported a period of maximum energy spectra smaller than 30 min.
The source spectrum was assumed by the ratio between the background spectra before and during the tsunami event. The spectral ratios are independent of the local effect on the background condition36,38. We estimated the spectral ratio for all stations, as shown in Fig. S5. Figure 5 shows the summary of the spectral ratios of all stations and the mean of all stations with the same sampling rate. Figure 5a1 shows clear energy spectra for the first duration of the atmospheric pressure wave (Lamb wave) in periods ranging from 2 to 150 min. The clear peak in the energy spectra was reported at periods of 7.7 and 30.1 min. The second duration of the Lamb wave showed a clear period in the energy spectra ranging from 2 to 150 min, as shown in Fig. 5a2. The reported peaks of the energy spectra were 7.5 and 35.4 min and were approximately 10 times smaller than peaks for the first duration. The spectral ratio of the first duration of the tsunami wave is presented in Fig. 5b1. The clear energy spectra for the first tsunami duration was a period ranging from 2 to 200 min, with peak energy spectra at periods of 9.9, 23.4, and 34.6 min. Figure 5b2 presents clear energy spectra for the second tsunami duration in the periods of 2 to 200 min. The clear periods of the peak energy spectra were 28.7, 45.6, and 130.7 min. It was concluded by considering the atmospheric pressure data that the energy spectra in periods shorter than 150 min were generated by the Lamb wave. For the tsunami, the energy spectra in a period shorter than 200 min were generated by the tsunami that was triggered by the 2022 HTHH volcanic eruption.
Spectral ratios for the atmospheric pressure wave and tsunami wave. Thin lines are the spectral ratio (ratio between wave spectra and background spectra) for each observation station, and thick lines represent the mean of the 12 stations. The top row of the panel presents the spectral ratio of the atmospheric pressure, and the bottom row panel presents the spectral ratio of the tsunami wave. The first column shows the spectral ratio of the first wave via the red line, and the second column shows the spectral ratio of the second wave in blue.
Modal analysis results
Eigenmodes of regional natural oscillations were calculated by considering the tsunami alone, as shown in Fig. 6. An amplification map of the eigenmode was used to highlight where the zone of increased energy would be concentrated19,25. The amplification map is ordered by period. In this study, we plotted the first four eigenmodes that were matched to the tsunami spectral analysis from all selected stations. The periods of the eigenmodes are marked by the vertical lines in Fig. 4 and shown in Table 1. We considered the difference between the periods derived from eigenmodes and from the spectral analysis of temporal data to be less than 10% for increasing amplification. The periods of the first four eigenmodes in the area around the Hanasaki station are 41.19, 23.23, 21.60, and 18.13 min, as shown by the amplification map in Fig. 6a1âa4. Figure 6b1âb4 shows that the periods of the first four eigenmodes in the area around the Kushiro station are 48.03, 34.59, 27.60, and 21.74 min. The mode of the 3rd period from Hanasaki station is close to the dominant periods from the spectral analysis of this station by less than 10%. The mode of the 1st period that is within 10% of the dominant periods from the spectral analysis was reported by six stations, Kushiro, Hakodate, Ofunato, Omaezaki, Kushimoto, and Aburatsu. Additionally, the mode of the 1st period exceeded the dominant period (as determined by spectral analysis) by 10% in Soma, and Tosashimizu stations. However, the 3rd mode from Soma station was within 10% of the dominant period of the spectral analysis. At three stations, Mera, Naha, and Ishigakijima, the eigenmode period did not match the dominant periods of the spectral analysis. The fundamental local mode of oscillation was also calculated for a specific tsunami, as shown in Table 1. All stations are in ports that open to the sea. The periods of the local oscillation modes were almost less than 30 min. There are only two stations where the local oscillation model was longer than 30 min, Hakodate and Naha. At the Hanasaki station, the local oscillation period was close to the dominant period from the spectral analysis, and the 3rd eigenmodes were within 10%. The local oscillation period of Kushimoto station was within 10% of the dominant period from the spectral analysis. At the Aburatsu and Naha stations, the local oscillation period was within 10% of the dominant period but differed the period of the eigenmode by more than 10%. It was concluded that increasing amplification might have been caused by the natural oscillation mode at eight stations: Hanasaki, Kushiro, Hakodate, Ofunato, Omaezaki, Kushimoto, Aburatsu, and Naha.
Modal analysis shows representative eigenmodes in the area around each station for the first four modes. The black triangle represents the observation station. The maps were generated using Python version 3.866 with the Matplotlib library67, and the basemap was downloaded from the QuickMapServices plugin68 through QGIS60.
Conclusion and discussion
The eruption of the HTHH volcano generated atmospheric pressure as the Lamb waves that caused disturbances around the earth39. Meteorological stations along the Japanese coastline have measured the Lamb wave multiple times, as shown in Figs. 2 and S2. After the HTHH volcanic eruption, the first Lamb wave arrived after approximately 7 h, as observed from the collected data. The distance from the HTHH volcano to Japan is approximately 8000 km, and, using a sound velocity of approximately 315 m/s2,40 reveals that the theoretical arrival time needed for the Lamb wave to reach Japan was approximately 7.1 h. This theoretical arrival time agrees well with the observation data. The shape of the pulse wave from the eruption is an N-wave shape15,41,42. N-wave shapes were also present in the observation data for the atmospheric pressure propagation in North America, Australia, New Zealand, and China1,3,8,11,15,43. The Lamb waves were also observed at approximately 00:00 UTC on 17 January 2022 approximately 36 h after the arrival time of the first Lamb wave. The second Lamb wave was the reflection of the first Lamb wave considering that the circumference of the Earth is approximately 40,075 km. The amplitude of the second Lamb wave was approximately 4 times smaller than the amplitude of the first Lamb wave. The wavelet spectra analysis also indicated that the first and second arrival times of the Lamb wave were 11:00 UTC on 15 January and 00:00 UTC on 17 January, respectively. There is no evidence that the second Lamb wave was amplified. The source spectra, taken as a ratio between the Lamb wave and the background of the two waves, revealed similar shapes and peak spectra that can be identified as originating from the same source.
The tsunami that impacted the coastline of Japan, as observed from the 12 selected stations, revealed that the arrival time was strongly related to the Lamb wave. The first tsunami wave was generated approximately 2â3Â h (approximately 10Â h after the HTHH eruption) after the first Lamb waveâs arrival, which was notable for all stations. The second tsunami wave occurred 1 to 2Â h after the second Lamb wave arrived, which can be seen in the data. Based on the gravity wave mechanism, the propagation of free waves that emanated from the HTHH volcano toward Japan took approximately 16Â h. The observed arrival time from the stations was faster than the times predicted in accordance with the theory of gravity wave, and this was also noted was by Wang et al.11 who noted that the wave impacted Linging Bay in China earlier than expected. The wavelet analysis also revealed the tsunami signal approximately 10Â h after the HTHH eruption. The wavelet analysis showed the amplification of tsunami height that occurred during the later phase. The source spectra of the tsunami were investigated by taking the spectral ratio of the durations of the two tsunami waves and the background. This revealed that both durations were related to a similar mechanism with similar shapes and peak spectra.
Based on the amplification of tsunami height in the later phase, it is assumed that the dominant period of the tsunami is close to the period of natural oscillation11,27. Then, the eigenmodes were used to determine the natural oscillation period. The dominant period of the tsunami was close to the period of natural oscillation at eight stations: Hanasaki, Kushiro, Hakodate, Ofunato, Omaezaki, Kushimoto, Aburatsu, and Naha. The four stations (Ofunato, Kushimoto, Aburatsu, and Naha) that showed that the dominant period was close to the natural period were located in the damaged area, as mentioned by Imamura et al.12. This may indicate that the amplification of tsunami resonance was one of the causes of damage in these four areas. The reported damage to marine vessels and aquaculture rafts during the event was a result of resonance12. The studyâs analysis supports the notion that these four areas align with the reported damage, and the dominant period of the tsunami event was close to the natural oscillation period. Therefore, it may be concluded that the amplification of tsunami resonance played a significant part in the damage incurred within the bay.
The causes of the volcanic tsunami have been explained by several factors. This studyâs results provided tsunami characteristics in some observed bays along Japanâs coastline. The tsunami was generated by atmospheric pressure via the Lamb wave that was triggered by the eruption of the HTHH volcano39. The results were limited to the affected region within the selected Bays. These results also support the assumption that the tsunami in this event was generated from the Lamb wave. The lead Lamb wave caused the water surface to rise10. The rise in water level, which represents the first tsunami wave, occurred concurrently over the deep sea due to the Lamb wave. The first tsunami wave was converted to a gravity wave over the continental slope10,44. Then, the maximum height of the first tsunami wave arrived approximately 3â6Â h after the Lamb wave impacted Japanâs coastal area. Water subsidence occurred and was converted into total energy after the Lamb wave had passed7. The inverse effect of water subsidence generated the gravity wave, which amplified the tsunami energy entering the bay.
This study focuses on using recorded data to study a wave resonance analysis of a tsunami generated by a volcanic explosionâs atmospheric pressure wave (Lamb wave). The study was limited to understanding the tsunami resonance characteristics due to the limited number of observation points in the bay. Using a simplified bay model assumed constant water depth for the local domain, which might not completely represent the complex bathymetry in the bay. Numerical simulations are necessary to improve a comprehensive understanding of the tsunami resonance characteristics in the bay entire. However, simulating this non-seismic tsunami event requires incorporating several unique effects into the numerical tsunami model, including the Earthâs rotation and the atmospheric pressure wave triggered by the volcanic explosion. The interaction between atmospheric pressure and tsunami waves might remain uncertain. Therefore, this study lacks numerical simulation results that could be used for comparison with the analysis results of this study. These limitations highlight the need for the following future study to understand tsunami wave characteristics in a deeper and more accurate representation of this specific event to contribute to tsunami science.
Methods
Spectral analysis
We performed two types of spectral analysis on the observed atmospheric pressure and tsunami height data recorded at each station: Wavelet and Fourier analysis. Wavelet analysis is used to determine variations in the dominant peak periods within the time series. The method has been widely applied in wave analysis to understand the temporal variation in wave energy11,19,45,46,47. The adopted wavelet is based on the PyWavelets library in the Python system48. The mother function for the applied wavelet was modeled by the Morlet function11,49,50. We conducted wavelet analysis to determine the real waveform (atmospheric pressure and tsunami height) of this event, in which the waveform time series was applied for 96 h, to present the spectral energy related to the period and time. Fourier analysis is used to reveal oscillation patterns and frequency behavior51,52. In this study, we conducted Fourier analysis to analyze both waveforms (the atmospheric pressure and tsunami height) and to identify their wave energy. The wave energy from the Fourier analysis was proceeded by two intervals, i.e., the wave that impacted during the tsunami and the background conditions, following the visible of atmospheric pressure pulses shown in Fig. 2. The Fourier analysis was conducted using the Numpy library in the Python programming language53. The tsunami wave was divided into 2 components, the first wave and the second wave. The first wave component occurred from January 15, 2022, 11:00 UTC to January 17, 2022, 00:00 UTC, and the second wave occurred from January 17, 2022, 00:00 UTC to January 18, 2022, 12:00 UTC. The time series of background component, as derived from the Fourier analysis, occurred from January 14, 2022, 00:00 UTC to January 15, 2022, 11:00 UTC. The background wave spectrum was used to reveal the natural oscillation periods that are a local effect54. To remove the local effect, the Lamb wave source spectrum and tsunami source spectrum were calculated by taking the ratio between the wave spectrum during the tsunami and its background spectrum19,54,55. In this study, we used this method to present the frequency characteristics of the Lamb wave-generating mechanism and the tsunami-generating mechanism, which is consistent with several previous studies38,56,57,58.
Modal analysis
The tsunami resonance characteristic is required to determine the dominant periods and their spatial characteristics when investigating the eigenmodes of natural oscillations59. The eigenmodes of natural oscillations are based on bathymetry data. The bathymetry data used were resampled with a resolution of 1.6 arc-second by using the cubic spline method in QGIS60. The bathymetry data were obtained from the Association for Promotion of Infrastructure Geospatial Information Distribution of Japan61. We investigated the tsunami resonance characteristics based on 2 domains, regional and local. The regional domains were used to provide information on large-scale modes in the area (10â30Â km) that covers the bay where the selected observed stations were located, as shown in detail in Fig. S6. The local domains were used to provide information on smaller-scale modes, such as in the area (1â3Â km) that covers the port in the bay for the regional domains, as shown in detail in Fig. S7. We calculated the eigenmodes of the natural oscillation in each selected bay (regional domain) using our developed numerical model that a governing equation is originally from Loomis62. The numerical model required only the bathymetry data. This model was solved using the linear solution to the shallow water equation (SWE) in a staggered grid system on the finite difference method. The governing equation of the model is as follows:
where \({\varphi }_{x,y,t}\) is the water level, \({h}_{x,y}\) is the water depth, and g is the gravitational acceleration (9.81 m/s2). The coastal boundary was assumed to be fully reflective and was without flux, and the open sea boundary was assumed to have no reflection conditions. The coordinate system of the original model is based on the Cartesian system that is unsuitable for depicting the actual surface of the earth. Therefore, we modified the governing equation of the original model to a spherical coordinate system based on Wu and Satake63. The differential equation given in Eq. (1), in the form of a matrix-eigenvalue problem, was modified to spherical coordinates in this study, as shown by Eq. (2):
where \(\theta \) is the latitude, \(\varphi \) is the longitude, R is the Earthâs radius, and λâ=âÏ2/g is the eigenvalue. The detail of the numerical method for this model was presented in supplementary S1.
We applied the method proposed by Rabinovich64 and Ivanov et al.65 to estimate the natural oscillation modes of the port considered for the local domain, which is where the selected stations were located. All stations were located in ports that face the open sea. In the method used for studying tsunami resonance in the port, the period of oscillation mode is determined for a rectangular basin with a constant water depth, which can be expressed as:
where L is the length of the port measured in the direction perpendicular to the port entrance and W is the width of the port parallel to the entrance. The details of both dimensions are shown in the left panel of Figs. S7 and 1. h is the average water depth in the port selected by the maximum probability based on the histogram, and the details used to estimate this variable are shown in the right panel of Fig. S7. In this study, we considered only the fundamental mode with kâ=â0 and mâ=â0. We compared the calculated eigenmode of regional and local domains with the dominant period of the spectral analysis. It is possible that damage was due to the increased amplification caused by periods that deviate by less than 10% from the period of eigenmodes, as inferred from spectral analysis of the time series19,27.
Data availability
The atmospheric pressure data used in this study are from Weathernews Inc. (https://global.weathernews.com/news/16551/). The sea water level data is available from Intergovernmental Oceanographic Commission (http://www.ioc-sealevelmonitoring.org) and Geospatial Information Authority of Japan (https://www.gsi.go.jp/kanshi/tide_furnish.html). The bathymetry data was provided from Association for Promotion of Infrastructure Geospatial Information Distribution of Japan (https://front.geospatial.jp/). The code for Wavelet analysis can be downloaded from PyWavelets (https://pywavelets.readthedocs.io/). The code for Fourier analysis is availed from Numpy (https://numpy.org/doc/stable/reference/routines.fft.html).
Code availability
The code that supports the findings in this study are available from the corresponding author upon resonable request.
References
Amores, A. et al. Numerical simulation of atmospheric lamb waves generated by the 2022 Hunga-Tonga volcanic eruption. Geophys. Res. Lett. 49, e2022GL098240. https://doi.org/10.1029/2022GL098240 (2022).
Nishida, K., Kobayashi, N. & Fukao, Y. Background lamb waves in the earthâs atmosphere. Geophys. J. Int. 196, 312â316. https://doi.org/10.1093/gji/ggt413 (2014).
Gusman, A. R. et al. The 2022 Hunga Tonga-Hunga Haâapai volcano air-wave generated tsunami. Pure Appl. Geophys. 179, 3511â3525. https://doi.org/10.1007/s00024-022-03154-1 (2022).
Monserrat, S., VilibiÄ, I. & Rabinovich, A. B. Meteotsunamis: Atmospherically induced destructive ocean waves in the tsunami frequency band. Nat. Hazard. Earth Syst. Sci. 6, 1035â1051. https://doi.org/10.5194/nhess-6-1035-2006 (2006).
Bechle, A. J. & Wu, C. H. The lake Michigan meteotsunamis of 1954 revisited. Nat. Hazards 74, 155â177. https://doi.org/10.1007/s11069-014-1193-5 (2014).
RamÃrez-Herrera, M., Coca, O. & Vargas-Espinosa, V. Tsunami effects on the coast of mexico by the Hunga Tonga-Hunga Haâapai volcano eruption, Tonga. Pure Appl. Geophys. https://doi.org/10.1007/s00024-022-03017-9 (2022).
Kubota, T., Saito, T. & Nishida, K. Global fast-traveling tsunamis by atmospheric pressure waves on the 2022 Tonga eruption. Science 377, 91â94. https://doi.org/10.1126/science.abo4364 (2022).
Ortiz-Huerta, L. G. & Ortiz, M. On the Hunga-Tonga complex tsunami as observed along the pacific coast of Mexico on January 15, 2022. Pure Appl. Geophys. 179, 1139â1145. https://doi.org/10.1007/s00024-022-03027-7 (2022).
Manneela, S. & Kumar, S. Overview of the Hunga Tonga-Hunga Haâapai volcanic eruption and tsunami. J. Geol. Soc. India. 98, 299â304. https://doi.org/10.1007/s12594-022-1980-7 (2022).
Tanioka, Y., Yamanaka, Y. & Nakagaki, T. Characteristics of the deep sea tsunami excited offshore Japan due to the air wave from the 2022 Tonga eruption. Earth. Planet Space. 74, 61. https://doi.org/10.1186/s40623-022-01614-5 (2022).
Wang, Y., Wang, P., Kong, H. & Wong, C. S. Tsunamis in Lingding Bay, China, caused by the 2022 Tonga volcanic eruption. Geophys. J. Int. 232, 2175â2185. https://doi.org/10.1093/gji/ggac291 (2022).
Imamura, F. et al. Preliminary observations and impact in Japan of the tsunami caused by the Tonga volcanic eruption on January 15. Pure Appl. Geophys. 179, 1549â1560. https://doi.org/10.1007/s00024-022-03058-0 (2022).
Wang, Y. et al. Data assimilation using high-frequency radar for tsunami early warning: A case study of the 2022 Tonga volcanic tsunami. J. Geophys. Res. Solid Earth 128, e2022JB025153. https://doi.org/10.1029/2022JB025153 (2023).
Wang, Y., Imai, K., Mulia, I. E., Kusumoto, S. & Takahashi, N. Tsunami early warning of the Hunga volcanic eruption using an ocean floor observation network off the Japanese islands. Seismol. Res. Lett. 94(2A), 567â577. https://doi.org/10.1785/0220220098 (2022).
Lynett, P. et al. Diverse tsunamigenesis triggered by the Hunga Tonga-Hunga Haâapai eruption. Nature 609, 728â733. https://doi.org/10.1038/s41586-022-05170-6 (2022).
La Selle, S. M. et al. Observations of tsunami and runup heights in santa cruz harbor and surrounding beaches from the 2022 Hunga Tonga-Hunga Haâapai tsunami, Doi: https://doi.org/10.5066/P9ZVAB8D (2022).
INDECI. Inician acciones de respuesta luego de oleajes en el litoral. retrieved. Retrieved https://www.gob.pe/institucion/indeci/noticias/576687-inician-acciones-de-respuesta-luego-de-oleajes-en-el-litoral (2022).
UNOCHA. Peru: Oil spill, flash update no. 02. retrieved. Retrieved https://reliefweb.int/report/peru/peru-oil-spill-flash-update-no-02-27-january-2022 (2022).
Wang, Y., Zamora, N., Quiroz, M., Satake, K. & Cienfuegos, R. Tsunami resonance characterization in Japan due to trans-pacific sources: Response on the bay and continental shelf. J. Geophys. Res. Ocean. 126, 2020JC017037. https://doi.org/10.1029/2020JC017037 (2021).
Allgeyer, S., Hébert, H. & Madariaga, R. Modeling the tsunami free oscillations in the marquesas (French Polynesia). Geophys. J. Int. 193, 1447â1459. https://doi.org/10.1093/gji/ggt064 (2013).
Heidarzadeh, M. & Satake, K. Excitation of basin-side modes of the pacific ocean following the March 2011 Tohoku tsunami. Pure Appl. Geophys. 171, 3405â3419. https://doi.org/10.1007/s00024-013-0731-5 (2014).
Wilson, B. W. Modeling the tsunami free oscillations in the marquesas (French polynesia). Adv. Hydrosci. 8, 1â94. https://doi.org/10.1016/B978-0-12-021808-0.50006-1 (1972).
Okal, E. A. et al. Oman field survey after the December 2004 Indian ocean tsunami. Earthq. Spectra 22, 203â218. https://doi.org/10.1193/1.2202647 (2006).
Munger, S. & Cheung, K. F. Resonance in hawaii waters from the 2006 kuril islands tsunami. Geophys. Res. Lett. 35, L07605. https://doi.org/10.1029/2007GL032843 (2008).
Yamazaki, Y. & Cheung, K. Shelf resonance and impact of near filed tsunami generated by the 2010 chile earthquake. Geophys. Res. Lett. 38, L12605. https://doi.org/10.1029/2011GL047508 (2011).
Zhao, X. & Liu, H. Resonance of long wave runup on south China seaâs piecewise topographies. J. Earthq. Tsunami. 11, 1740003. https://doi.org/10.1142/S1793431117400036 (2017).
Cortés, P., Catalán, P., Aránguiz, R. & Bellotti, G. Tsunami and shelf resonance on the northern chile coast. J. Geophys. Res. 122, 7364â7379. https://doi.org/10.1002/2017JC012922 (2017).
Abe, K. Tsunami resonance curve from dominant periods observed in bays of northeastern japan. Adv. Nat. Technol. Hazards Res. 23, 97â113. https://doi.org/10.1007/1-4020-3331-1_6 (2005).
Inc., W. Atmospheric pressure data from soratena weather sensors offered to researhers free of charge. retrieved. Retrieved https://global.weathernews.com/news/16551/ (2022).
Haurwitz, B. & Cowley, A. D. The diurnal and semidiurnal barometric oscillations global distribution and annual variation. Pure Appl. Geophys. 102, 193â222. https://doi.org/10.1007/BF00876607 (1973).
VLIZ. Flanders marine institute, intergovernmental oceanographic commission (IOC) (2022): Sea level station monitoring facility. Accessed http://www.ioc-sealevelmonitoring.org (2022).
Geospatial Information Authority of Japan, G. Geospatial information authority of japan. Accessed https://www.gsi.go.jp/kanshi/tide_furnish.html (2022).
Pakoksung, K., Suppasri, A. & Imamura, F. 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. Sci Rep. 12, 15187. https://doi.org/10.1038/s41598-022-19486-w (2022).
Heidarzadeh, M. & Satake, K. Waveform and spectral analyses of the 2011 Japan tsunami records on tide gauge and dart stations across the Pacific ocean. Pure Appl. Geophys. 170, 1275â1293. https://doi.org/10.1007/s00024-012-0558-5 (2013).
Ren, Z., Hou, J., Wang, P. & Wang, Y. Tsunami resonance and standing waves in Hangzhou bay. Phys. Fluids 33, 081702. https://doi.org/10.1063/5.0059383 (2021).
Rabinovich, A. B., Titov, V. V., Moore, C. W. & Eblé, M. C. The 2004 Sumatra tsunami in the southeastern Pacific ocean: New global insight from observations and modeling. J. Geophys. Res. Ocean. 122, 7992â8019. https://doi.org/10.1002/2017JC013078 (2017).
Heidarzadeh, M. & Satake, K. The El Salvador and Philippines tsunamis of August 2012: Insights from sea level data analysis and numerical modeling. Pure Appl. Geophys. 171, 3437â3455. https://doi.org/10.1007/s00024-014-0790-2 (2014).
Toguchi, Y., Fujii, S. & Hinata, H. Tsunami waves and tsunami-induced natural oscillations determined by HF radar in Ise bay, Japan. J. Geophys. Res. Ocean. 123, 2965â2980. https://doi.org/10.1029/2017JC013626 (2018).
Nishikawa, Y. et al. Observation and simulation of atmospheric gravity waves exciting subsequent tsunami along the coastline of Japan after Tonga explosion event. Sci. Rep. 12, 22354. https://doi.org/10.1038/s41598-022-25854-3 (2022).
Abbrescia, M. et al. Observation of Rayleigh-Lamb waves generated by the 2022 Hunga-Tonga volcanic eruption with the pola detectors at ny-Ã lesund. Sci. Rep. 12, 19978. https://doi.org/10.1038/s41598-022-23984-2 (2022).
Fitzgerald, T. J. Observations of total electron content perturbations on GPS signals caused by a ground level explosion. J. Atmos. Sol. Terr. Phys. 59, 829â834. https://doi.org/10.1016/S1364-6826(96)00105-8 (1997).
Matoza, R. S. et al. Atmospheric waves and global seismoacoustic observations of the January 2022 Hunga eruption, Tonga. Science 337, 95â100. https://doi.org/10.1126/science.abo7063 (2022).
Hu, G., Li, L., Ren, Z. & Zhang, K. The characteristics of the 2022 Tonga volcanic tsunami in the pacific ocean. Nat. Hazards Earth Syst. Sci. 23, 675â691. https://doi.org/10.5194/nhess-23-675-2023 (2022).
Tonegawa, T. & Fukao, Y. Wave propagation of meteotsunamis and generation of free tsunamis in the sloping area of the japan trench for the 2022 HungaâTonga volcanic eruption. Earth Planets Space. 74, 159. https://doi.org/10.1186/s40623-022-01727-x (2022).
Nagarajan, B. et al. The great tsunami of 26 December 2004: A description based on tide-gauge data from the Indian subcontinent and surrounding areas. Earth Planets Space. 58, 211â215. https://doi.org/10.1186/BF03353380 (2006).
Heidarzadeh, M. & Satake, K. The 21 May 2003 tsunami in the western mediterranean sea: Statistical and wavelet analyses. Pure Appl. Geophys. 170, 1449â1462. https://doi.org/10.1007/s00024-012-0509-1 (2013).
VilibiÄ, I. & Å epiÄ, J. Global mapping of nonseismic sea level oscillations at tsunami timescales. Sci. Rep. 7, 40818. https://doi.org/10.1038/srep40818 (2017).
Lee, G. R., Gommers, R., Wasilewski, F., Wohlfahrt, K. & OâLeary, A. Pywavelets: A python package for wavelet analysis. J. Open Source Softw. 4, 1237. https://doi.org/10.21105/joss.01237 (2019).
Goupillaud, P., Grossmann, A. & Morlet, J. Cycle-octave and related transforms in seismic signal analysis. Geoexploration 23, 85â102. https://doi.org/10.1016/0016-7142(84)90025-5 (1984).
Torrence, C. & Compo, G. P. A practical guide to wavelet analysis. Bull. Am. Meteorol. Soc. 79, 61â78. https://doi.org/10.1175/1520-0477 (1998).
Ren, Z., Liu, H., Jimenez, C. & Wang, Y. Tsunami resonance and standing waves in the south China sea. Ocean. Eng. 262, 112323. https://doi.org/10.1016/j.oceaneng.2022.112323 (2022).
Catalán, P. A. et al. The 1 April 2014 Pisagua tsunami: Observations and modeling. Geophys. Res. Lett. 42, 2918â2925. https://doi.org/10.1002/2015GL063333 (2015).
Harris, C. R. et al. Array programming with numpy. Nature 585, 357â362. https://doi.org/10.1038/s41586-020-2649-2 (2020).
Rabinovich, B. Spectral analysis of tsunami waves: Separation of source and topography effects. J. Geophys. Res. 102, 12663â12676. https://doi.org/10.1029/97JC00479 (1997).
Heidarzadeh, M., Harada, T., Satake, K., Ishibe, T. & Gusman, A. R. Comparative study of two tsunamigenic earthquakes in the solomon islands: 2015 mw 7.0 normal-fault and 2013 santa cruz mw 8.0 megathrust earthquakes. Geophys. Res. Lett. 43, 4340â4349. https://doi.org/10.1002/2016GL068601 (2016).
Vich, M. & Monserrat, S. Source spectrum for the Algerian tsunami of 21 May 2003 estimated from coastal tide gauge data. Geophys. Res. Lett. https://doi.org/10.1029/2009GL039970.10.1029/2009GL039970 (2003).
Zaytsev, O., Rabinovich, A. B. & Thomson, R. E. A comparative analysis of coastal and open-ocean records of the great chilean tsunamis of 2010, 2014 and 2015 off the coast of mexico. Pure Appl. Geophys. 173, 4139â4178. https://doi.org/10.1007/s00024-016-1407-8 (2016).
Rabinovich, A. B., Fritz, H. M., Tanioka, Y. & Geist, E. L. Global tsunami science: Past and future. Pure Appl. Geophys. 174, 2883â2889. https://doi.org/10.1007/s00024-017-1638-3 (2017).
Bellotti, G., Briganti, R. & Beltrami, G. M. The combined role of bay and shelf modes in tsunami amplification along the coast. J. Geophys. Res. Ocean. 117, C08027. https://doi.org/10.1029/2012JC008061 (2012).
QGIS. Qgis geographic information system. QGIS Association http://www.qgis.org (2022).
Geospatial. Association for promotion of infrastructure geospatial information distribution of japan. Geospatial https://front.geospatial.jp/ (2023).
G., L. H. Normal modes of oscillation of honokohau harbor. (noaa-jtre-142, honolulu, hi: Hawaii institute of geophysics, university of hawaii. Loomis https://books.google.co.jp/books?id=3_F2OwAACAAJ (1975).
Wu, Y. & Satake, K. Synthesis and source characteristics of tsunamis in the sea of Japan based on normal-mode method. J. Geophys. Res. Solid Earth. 123, 5760â5773. https://doi.org/10.1029/2018JB015707 (2018).
Rabinovich, A. B. Seiches and harbor oscillations. In Handb. Coast. Ocean. Eng. (ed. Kim, Y. C.) 193â236 (World Sci Publ. Co. Pte. Ltd, 2010). https://doi.org/10.1142/9789812819307_0009.
Ivanov, V. A., Manilyuk, Y. V. & Sannikov, V. F. Seiches in a basin with an open entrance. J. Appl. Mech. Tech. Phys. 59, 594â600. https://doi.org/10.1134/S0021894418040041 (2018).
van Rossum, G. & Drake, F. L. Python 3 reference manual. Loomis http://www.python.org (Createspace, 2009).
Hunter, J. D. A 2D graphics environment. Comput. Sci. Eng. 9, 90â95. https://doi.org/10.1109/MCSE.2007.55 (2018).
NextGIS. Quickmapservices: Easy basemaps in qgis. Loomis https://nextgis.com/blog/quickmapservices/ (2019).
Acknowledgements
The authors gratefully acknowledge Association for Promotion of Infrastructure Geospatial Information Distribution of Japan for the bathymetry data, Weathernews Inc. for the atmospheric data, and Intergovernmental Oceanographic Commission and Geospatial Information Authority of Japan for the sea level data. This research was funded by the Arthur J. Gallagher & Co. global group of companies, and by financial support from Pacific Consultants Co., Ltd. through the International Research Institute of Disaster Science (IRIDeS) at Tohoku University. The authors would like to thank the editor of this journal and two reviewers.
Author information
Authors and Affiliations
Contributions
All authors read, reviewed, and approved the manuscript. Conceptualization, methodology, investigation, visualization, data collection, and writing: K.P. Supervision and funding: A.S. and F.I.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Pakoksung, K., Suppasri, A. & Imamura, F. Resonance characteristics of tsunami in bay of Japan by the Hunga Tonga-Hunga Haâapai volcano eruption on 15th January 2022. Sci Rep 13, 18385 (2023). https://doi.org/10.1038/s41598-023-45601-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-45601-6
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.