Combined hazard of typhoon-generated meteorological tsunamis and storm surges along the coast of Japan

Two hazardous typhoons, Lionrock (August 2016) and Jebi (September 2018), destructively affected the coast of Japan and produced extreme sea level variations. The results of field surveys in the impacted regions showed that multiple deaths and extensive floods were caused by the combined effect of low-frequency sea level raise (storm surges) and intensive high-frequency (HF) tsunami-like waves (meteotsunamis). The data from ten tide gauges for the 2016 event and eight gauges for the 2018 event were used to examine the properties of the observed sea levels, to estimate the relative contribution of the two sea level components and to evaluate their statistical characteristics (maximum wave heights, amplitudes and periods of individual components, etc.). For the 2016 event, we found that the surge heights were from 12 to 35 cm and that the mean contribution of surges into the total observed sea level heights was ~ 39%; the meteotsunami amplitudes were from 22 to 92 cm, and they contributed 61% of the total height. For the 2018 event, storm surges were significantly stronger, from 46 to 170 cm, while HF amplitudes were from 38 to 130 cm; their relative inputs were 67% and 33%, respectively. Combined, they formed total flood heights of up to 120 cm (2016 event) and 288 cm (2018 event). Previously, the contribution of storm seiches (meteotsunamis) in coastal floods had been underestimated, but results of the present study demonstrate that they can play the principal role. What is even more important, they produce devastating currents: according to our estimates, current speeds were up to 3 knots (1.5 m/s) during the Lionrock event and more than 5 knots (2.6 m/s) during Jebi; these strong currents appear to be the main reason for the resulting damage of coastal infrastructure. The most important characteristic of the recorded meteotsunamis is their trough-to-crest maximum height. During the 2016 event, these heights at three stations were > 1 m: 171 cm at Erimo, 109 cm at Hachijojima and 102 cm at Ayukawa. The 2018 event was stronger; maximum meteotsunami wave heights were 257 cm at Gobo, 138 cm at Kushimoto, 137 cm at Kumano and 128 cm at Murotomisaki. The 2018 Gobo height of 257 cm is much larger than historical non-seismic seiche maxima for the Pacific coast of Japan (140–169 cm) estimated by Nakano and Unoki (1962) for the period of 1930–1956.


Introduction
Meteorological tsunamis (meteotsunamis) are strong tsunami-like ocean waves generated by atmospheric processes. They have approximately the same periods and spatial scales as ordinary tsunami waves and can affect coastal areas in a similar destructive way (cf. Monserrat et al. 2006;Vilibić et al. 2014). In recent years, meteotsunamis have attracted significant attention from scientists and specialists. A number of devastating meteotsunami events occurred in various parts of world oceans, including the coasts of North and South America, SE Asia, South Africa and Australia, as well as the Great Lakes (Pattiaratchi and Wijeratne 2015;Rabinovich 2020). The 3.5-m Dayyer meteotsunami of 19 March 2017 killed five people and severely destructed local infrastructure on the coast of the Persian Gulf (Salaree et al. 2018;Heidarzadeh et al. 2020a). In bays, inlets, harbours and lakes, meteotsunamis, as well as seismically generated tsunamis, have the same character of eigen oscillations of the corresponding basins, i.e. "extreme seiches" (Rabinovich 2009(Rabinovich , 2020. Various types of atmospheric disturbances can produce meteotsunamis. In general, all meteotsunamis may be divided into two principal groups: "good-weather" and "badweather" (Rabinovich 2020). The events of the first group are mostly generated by smallscale air-pressure disturbances, or by propagating atmospheric gravity waves (cf. Monserrat et al. 1991;Monserrat and Thorpe 1996); the weather during the corresponding events is commonly reported as "nice" and "quiet". "Good-weather" meteotsunamis are very typical for the Mediterranean region (cf. Rabinovich and Monserrat 1996;Vilibić and Šepić 2009;Šepić et al. 2015). In contrast, "bad-weather" meteotsunamis are normally generated by incoming hurricanes, typhoons and deep cyclones and are associated with stormy weather conditions. Extreme storm seiches with trough-to-crest wave heights > 30-40 cm are generally identified as meteotsunamis (cf. Pattiaratchi and Wijeratne 2015;Olabarrieta et al. 2017;Dusek et al. 2019); they are common on the Atlantic coasts of the USA and western Europe, on the oceanic coast of Japan, and in the Great Lakes (de Jong et al. 2003;Šepić and Rabinovich 2014;Olabarrieta et al. 2017;Dusek et al. 2019;Tanaka 2019).
There are important differences between these two types of meteotsunami: 1. "Good-weather" meteotsunamis are typically local events. They have a resonant nature and normally occur in a specific bay/harbour or in a few neighbouring bays (cf. Vilibić and Šepić 2009). There are some "hot spots" in the world oceans, where such events take place regularly and achieve extreme heights: in particular, Ciutadella Harbour, Menorca Island, Spain (cf. Gomis et al. 1993;Rabinovich and Monserrat 1996); Vela Luka Bay, Korčula Island, Croatia (cf. Orlić et al. 2010;Vilibić et al. 2016); Nagasaki Bay, Japan (Honda et al. 1908;Hibiya and Kajiura 1982;Tanaka 2010). 2. Storm seiches caused by propagating hurricanes or typhoons are "Bad-weather meteotsunamis" and considered as large-scale events that can be observed along large segments of the mainland coastline, for example along the entire East Coast of USA (cf. Thomson et al. 2007;Pasquet et al. 2013;Dusek et al. 2019).
While good-weather meteotsunamis are individual hazardous events, "bad-weather" meteotsunamis often coincide with other dangerous types of long-wave oscillations, in particular, with storm surges (cf. Thomson et al. 2007;Roeber and Bricker 2015). The superposition of waves from these two natural phenomena strongly increases the damaging effect of this hazard.
Typhoons are usually associated with a number of hazards in a cascading order, including severe storm waves, wind-blown debris, storm surges, extreme seiches, intensive rainfall, landslides and debris flows (cf. Heidarzadeh et al. 2018Heidarzadeh et al. , 2020b. Among these hazards, the two common and coinciding ones are seiches (meteotsunamis) and storm surge. Both meteotsunamis and storm surges are atmospherically induced events. Therefore, for scientists and authorities working on hazard mitigation and coastal engineering, it appears to be difficult to distinguish the differences between these two phenomena (Rabinovich 2020). However, meteotsunamis and storm surges have quite dissimilar temporal and spatial scales and principally have different generation mechanisms. A storm surge is a lowfrequency process, with typical periods/durations from several hours to approximately one week; they can flood extensive coastal areas of several hundred kilometres (Pugh and Woodworth 2014). Meteotsunamis have the same periods and other properties as ordinary tsunami waves; in harbours and bays ("harbour meteotsunamis") they have a scale of the corresponding local basins and periods determined exactly by natural (eigen) periods of the basins. This means that their spatial scales are from a few hundred metres to a few tens of kilometres and their periods from 1 to 2 min to approximately two hours (Rabinovich 2020).
Storm surges are formed by the combined effect of low air pressure and wind, where the wind generally plays the determining role in producing destructive surges (cf. Pugh and Woodworth 2014). In contrast, meteotsunamis are mainly induced by direct air pressure forcing and are usually related to propagating atmospheric disturbances (Monserrat et al. 2006;Šepić et al. 2015). The principal differences between meteotsunamis and surges are discussed in detail by Rabinovich (2020). What is crucial, storm surges and seiches are typically combined together. During strong storms, normally associated with incoming hurricanes, typhoons and deep cyclones, background coastal seiches are considerably intensified and can be identified as "meteotsunamis" (cf. Olabarrieta et al. 2017;Dusek et al. 2019;Heidarzadeh et al. 2018). Such extreme storm seiches and associated reversing currents can lead to the breaking of mooring lines, fenders and piles and can seriously damage anchored boats (Rabinovich 2009). The damaging effect of storm surges is strongly amplified by their superposition with storm seiches/meteotsunamis (cf. Thomson et al. 2007); however, this effect is still underestimated or, at least, is poorly investigated. The main challenge for a combined examination is the very different temporal/spatial scales of these two phenomena; in particular, the sampling interval of 1 h for sea level data is good enough to operatively analyse storm surges (Pugh and Woodworth 2014), while for seiches/meteotsunamis the optimum sampling should be two orders finer (0.5-2.0 min) (Vilibić et al. 2016). Similarly, the numerical modelling of a storm surge can be effectively done based on a computational grid with 0.5 arc min resolution (~ 0.9 km) (Pugh and Woodworth 2014), whereas to reproduce actual meteotsunami oscillations we need a 0.5-1.0 arc sec grid (15-30 m) (Orlić et al. 2010;Rabinovich et al. 2020). That is why, there are only a few papers considering these two processes together (cf. Roeber and Bricker 2015;Bajo et al. 2019). In the present study, this exact effect, the combination of storm-generated meteotsunamis and storm surges, is examined for two devastating typhoons affecting the coast of Japan in 2016 and 2018.
Japan is the country most impacted by hazardous tsunami waves. At the same time, pronounced atmospherically induced tsunami-like seiches are also common in multiple bays and inlets of the Japanese coast (cf. Honda et al. 1908). Especially intense seiches, i.e. meteotsunamis, are generated by severe typhoons (with air pressure as low as 920-960 hPa), which destructively impact Japan in the summer and fall seasons. In fact, the term "meteorological tsunamis" originally came from Japan and was specifically related to 1 3 ruinous tsunami-like waves produced by typhoons (Nomitsu 1935). Severe typhoons are characterized by very low air pressure and are frequently accompanied by sustained hurricane onshore winds. These two factors cause powerful surges on the oceanic coasts of Japan. At the same time, this extreme atmospheric forcing produces strong storm seiches superimposed on surges. Heidarzadeh et al. (2020b) and Le et al. (2019) thoroughly examined the impacts of two recent catastrophic typhoons on the Pacific coast of Japan: the 2016 Lionrock and 2018 Jebi typhoons. The primary focuses of their studies were field surveys of the affected territories and numerical modelling of generated storm surges. The results of their studies showed that multiple deaths and extensive coastal floods were partly caused by intensive high-frequency tsunami-like waves observed at several sites. Here, we further expand those studies by extracting storm surge and seiche signals from tide gauge records, estimating the effect of their superposition and the potential risk to coastal communities.

The August 2016 Typhoon Lionrock
Typhoon Lionrock was a destructive powerful tropical cyclone that on 25-31 August 2016 strongly affected the coasts of the Philippines, Japan, North Korea and Russia. The total damage in these four countries was estimated to be approximately US $4.0 billion; the typhoon killed a total of about 550 people. The typhoon caused intense floods in coastal regions associated with extreme rainfalls, storm surges and seiches.

Description of the event on the east coast of Japan
The 2016 Typhoon Lionrock, also known as Typhoon No. 1610 in Japan (for being 10th in the 2016 typhoon season), was active in the Pacific Ocean region during the period of 16 − 31 August 2016. Lionrock made its main landfall on the coast of Japan on 30 August 2016 (Fig. 1) causing a series of hazardous effects, including severe flooding, hurricane winds and landslides; altogether 26 people were killed (FDMA 2017).
Lionrock was classified as a Category-5 hurricane in the Saffir-Simpson Hurricane Wind Scale (SSHWS) by reaching a maximum one-minute sustained wind speed of 60 m/s. As seen in Fig. 1, Lionrock travelled along a complicated track-it originated from the SE of Japan, moved towards the NW, then changed its direction to SW and after a full loop finally again began to propagate north-westward to make landfall in east Japan. The minimum typhoon pressure was 940 hPa (Nayak and Takemi 2019). Approaching the coast of Japan, Lionrock began to slightly weaken; it made landfall near Ōfunato (Iwate Prefecture) on 30 August; maximum winds at that time were up to 70 knots (35 m/s). As an example, Fig. 2a shows sea level pressure (SLP) records at Ofunato and nearby station Ishinomaki. Unfortunately, we did not have high-resolution air-pressure (AP) records at these or other Japanese stations that could probably uncover small-scale AP disturbances associated with this typhoon; normally such disturbances are the main source of meteotsunamis (cf. Monserrat et al. 2006;Rabinovich et al. 2020). Therefore, we used standard JMA 1 AP measurements with 1-h sampling (Fig. 2a), which revealed the marked V-shape of the AP records during the typhoon passage over the respective 1 3 sites, with sharp SLP minima of 970 hPa (Ofunato) and 972 hPa (Ishinomaki) at 16:00 UTC on 30 August. A satellite image of the Lionrock typhoon ( Fig. 2b) indicates that on 28 August 2016 the diameter of the Lionrock system was approximately 900 km.
Field survey of the typhoon-affected areas was conducted in a few days after the event (Heidarzadeh et al. 2020b). Lionrock caused widespread damage through its various cascading effects as shown in Fig. 3. Typhoon-triggered severe flooding brought extensive mud into the city of Miyako ( Fig. 3a; the location of Miyako is indicated in Fig. 1). In Ryouri Bay, a part of the caisson breakwater was damaged (Fig. 3b). An under-construction bridge in Miyako was partially damaged by extensive flood-driven debris (Fig. 3c), while wave-driven debris piled up at the coast of Miyako (Fig. 3d). The preliminary examination of the flooding and destruction demonstrated that they were associated with both long-period sea level rise and intensive short-period oscillations with periods from a few minutes to approximately 15-20 min (cf. FDMA 2017; Le et al 2019). Thus, the purpose of our further study was to examine both types of water motions.

The combined effect of storm surges and meteotsunamis
Sea level oscillations caused by typhoon Lionrock were recorded by a great number of tide gauges. For our analysis, we selected ten stations: four on the east and south-east coast of Hokkaido Island, three on the east coast of Honshu Island and three located on isolated Pacific islands (Fig. 1). The data were provided by the Japan Meteorological Agency (JMA); all records had a sampling interval of 15 s, except Ofunato which had 1-min sampling. The height resolution was 1 cm, except Erimo, Misaki-gyoko and Hachijojima, which had a resolution of several cm. The sea level data underwent quality control to remove gaps, spikes and outliers. The tides were calculated using the least-squares method (cf. Pugh and Woodworth 2014) and subtracted from the original records based on MAT-LAB T_Tide software (cf. Pawlowicz et al. 2002;Heidarzadeh et al. 2017); the resulting residual time series were used in all subsequent analyses.
As an example, Fig. 4a shows residual (de-tided) records at two stations: Kushiro and Ayukawa. The intensification of sea level variations during the typhoon passage is clearly seen in the records (Compare Fig. 4 with the Lionrock track shown in Fig. 1 and with AP records in Fig. 2a). At Ayukawa, the maximum variations were observed at about 15:00 UTC on 30 August 2016, at Kushiro approximately 9 h later. Two types of sea level changes are evident in the records: (1) low-frequency sea level rise with a duration of ~ 2 days and (2) high-frequency oscillations superimposed on the low-frequency rise. The first type of change is a storm surge (Pugh and Woodworth 2014), which appears to be produced by the combined effect of very low air pressure and surge wind. The second type are storm seiches (Rabinovich 2009), which are so strong that they can be considered as meteotsunamis.
To separate these two types of sea level oscillations, we used a Kaiser-Bessel filter (cf. Thomson and Emery 2014) with a 3-h window. This window is commonly used to isolate tsunami waves and to separate storm surges and seiches (cf. Thomson et al. 2007;Rabinovich 2020). The black curve in Fig. 4a contours the storm surge that is also isolated by shading. Seiches (meteotsunamis) for the segments denoted in Fig. 4a by grey bands are shown in more detail in Fig. 4b. Both the storm surge and the storm seiches are atmospherically generated processes specifically related to the typhoon's passage, but as seen in Fig. 4, their character and the main features are in strong contrast. The storm surge at Kushiro and Ayukawa, as well as at other examined stations, for this event was not substantial, but combined with storm seiches it became hazardous.

Observed sea level oscillations
The de-tided unfiltered and filtered records were used to estimate the main parameters of the observed oscillations ( Fig. 4): • The total observed height above mean sea level (MSL), i.e. the height of the surge together with seiches. • The surge height relative to MSL. • The maximum seiche (meteotsunamis) amplitude relative to MSL. • The observed dominant period of seiches. Residual (de-tided) sea level records at stations Kushiro and Ayukawa on the east coast of Japan during passage of the 2016 Lionrock typhoon. a Partitioning of the records into low-and high-frequency oscillations using a 3-h Kaiser-Bessel window. The thick black lines denote the low-frequency components of the records; the shaded transparent region denotes the storm surge. The maximum observed (total) amplitude and the amplitude of the surge are indicated; and b expanded versions of the high-frequency (stormgenerated) seiches at the two stations for the record segments shaded grey in a. The maximum amplitudes and trough-to-crest wave heights are indicated • The maximum trough-to-crest wave seiche (meteotsunami) height.
The latter parameter is important for meteotsunami research, since the maximum wave height is commonly used to identify and characterize meteotsunami events (cf. Rabinovich and Monserrat 1996;Pasquet et al. 2013;Dusek et al. 2019).
The parameters listed above have been estimated for all ten stations; the corresponding results are shown in Figs. 5 and 6 and presented in Table 1. The maximum observed (total) heights were from 33 cm (Hanasaki) to 120 cm (Erimo); the mean height (averaged over all 10 stations) was 66 cm. The maximum surge heights were significantly smaller, 12 − 35 cm, and the mean surge height was ~ 25 cm (Fig. 5). The relatively amplified surge heights of 30-35 cm were at four stations located on the typhoon track (Ayukawa and Ofunato) or on the right-hand side of this track (Hakodate and Erimo); the smallest heights of 12-14 cm were at two northernmost stations, Hanasaki and Kushiro, and on the remote easternmost island Chichijima (see Fig. 1 for the station locations). The seiche amplitudes were considerably larger, from 22 cm (Hanasaki) to 92 cm (Erimo) (Fig. 6); the mean seiche amplitude was ~ 44 cm. The remarkable fact is that at nine from ten tide gauges, the seiche amplitudes were larger than the surge height; the only exception is Hakodate where the seiche amplitude of 26 cm was slightly smaller than the surge height of 35 cm. The seiche amplitudes vary significantly from one station to another, strongly controlling the hazard level posed by the typhoon floods at each location. For instance, while both Erimo and Hakodate recorded the same surge height of 35 cm (Fig. 5), the maximum observed height at Erimo was 120 cm, while at Hakodate it was only 57 cm because Erimo experienced much stronger seiching than Hakodate. It is noteworthy that Erimo was not located directly on the track of typhoon Lionrock, but was distanced approximately 300 km northward from the track (Fig. 1a). We can assume that Erimo is a "hot spot" (similar to Atlantic City on the US East Coast, cf. Šepić and Rabinovich 2014) and that some local topographic properties (e.g. the specific resonant properties of the adjacent shelf, funnelling effect of the coast, and high Q-factor of the respective inner basin) of this site promote strong sea level response to the propagating typhoon and fiercely amplified seiches. In general, the results of our analysis indicate that for typhoon Lionrock seiches were a greater hazard than the surge, while surges are commonly assumed to be the main reason for coastal floods when typhoons pass over the coast of Japan, or when hurricanes affect the US East Coast.
At all stations, the observed total heights were smaller than the sum of the maximum surge height and seiche amplitude ( Table 1). The reason for this is the time lag between the occurrence of the corresponding maxima (Fig. 5). Nevertheless, in most cases the difference was less than 1-3 cm. At two stations, Erimo and Misaki-gyoko, the difference was 7 cm, but the relative difference was rather small: 5.5% and 8.7%, respectively. This means that the cumulative effect of the storm surge and seiches can roughly be counted as the sum of each of these phenomena. The data in Table 1 demonstrate that seiche heights at most stations were so large that they can be considered as meteotsunamis. In particular, at Erimo they were 171 cm, which is comparable with the greatest meteotsunamis recorded during the last 27 years in various regions of the world oceans (Rabinovich 2020) At two other stations, Ayukawa and Hachijojima, they were higher than 1 m and at four more, Kushiro, Ofunato, Misaki-gyoko and Miyakejima, higher than 65 cm. These heights are similar to those for the strongest hurricane-generated meteotsunamis on the coasts of the Gulf of Mexico and the East Coast of the USA (Olabarrieta et al. 2017;Dusek et al. 2019).
Based on visual analysis of the records themselves and on time-frequency (wavelet) analysis (see next section) we estimated the dominant periods of the recorded seiches were 45 min at Hakodate and 16 and 20 min at Hanasaki; at other stations, the observed seiche periods were short, mostly shorter than 8 min (Table 1). These are shorter than typical periods of meteotsunamis (10-40 min) associated with propagating small-scale atmospheric disturbances (Monserrat et al. 2006;Rabinovich 2020). High-frequency oscillations  The high-pass filtered (with a 3-h Kaiser-Bessel window) records of sea level oscillations at ten stations on the coast of Japan showing storm seiches (meteotsunamis) caused by the 2016 typhoon Lionrock. The purple and red numbers denote the maximum amplitudes and trough-to-crest wave heights, respectively

Time-frequency analysis
To examine temporal variations of the recorded seiche oscillations associated with the 2016 Lionrock typhoon in the frequency domain, we used a multiple-filter method, which is similar to wavelet analysis (cf. Thomson and Emery 2014;Thomson et al 2007;Heidarzadeh et al. 2020cHeidarzadeh et al. , 2017 and is based on narrow-band filters with a Gaussian window that isolates a specific central frequency, n = 2 f n , where f n is the frequency. This method has been effectively used previously to investigate various meteotsunami events (cf. Šepić and Rabinovich 2014;Šepić et al. 2015). The method enables us to determine the evolution of meteotsunami waves as a function of frequency, f, and time, t, and to build so called f-t diagrams that display possible non-stationary character of the observed waves.
We selected four stations, Kushiro, Ayukawa, Miyakejima and Chichijima, and for these stations constructed the respective "f-t diagrams" for the frequency band of 1-60 cph (cycles per hour), i.e. for the periods from 1 h to 1 min. We used 3-day data segments: 29 August-1 September 2016 for the first two stations (the more northerly stations) and 28-31 August 2016 for the other two stations. The most prominent feature of the resulting plots, shown in Fig. 7, are well defined relatively narrow frequency bands of significantly amplified energy, which are likely associated with the eigen frequencies of the respective harbours. These "bands" in the f-t diagrams approximately relate to the following periods: 4, 5 and 17 min at Kushiro; 7.5 min at Ayukawa; 4 min at Miyakejima; 1, 4 and 17 min at Chichijima. The same "bands" are clearly seen in background oscillations (before the typhoon arrival); however, during the event the respective oscillations strongly intensified. High-frequency stability and the persistent character of the corresponding "band" oscillations during the entire examined periods support the assumption of their resonant ("eigen") nature. For a few hours around the core time of the typhoon attack we can see dramatic amplification of these specific resonant oscillations and, at the same time, active generation of broad-band oscillations occupying almost the entire frequency band of 1-60 cph (Fig. 7). After the typhoon's passage, the broad-band oscillations decayed very fast, while the narrow-band resonant oscillations kept ringing. In general, the f-t analysis demonstrates evident sea level response to the propagating typhoon and the high importance of the local resonant properties of individual harbours in the formation of hazardous floods. Storm-related infragravity (IG) waves with periods < 5 min are another important factor amplifying high-frequency seiches at specific sites (for example, at Chichijima).

The September 2018 Typhoon Jebi
Typhoon Jebi was the most powerful typhoon to strike Japan in 25 years. It was formed as a tropical depression on 26 August 2018 and then rapidly intensified into a typhoon. This typhoon destructively affected the Northern Mariana Islands and then on 2 September caused huge waves and five fatalities on the coast of Taiwan. Propagating north-westward and then northward on 4 September, Jebi reached the south-eastern coast of Japan producing severe damage and intensive floods along the entire coasts of Shikoku Island and the Kansai region of Honshu Island. The typhoon existed for 15 days and ended on 9 September 2018.

Description of the event on the south-east coast of Japan
Typhoon Jebi, one of the largest typhoons to ever affect Japan, left 13 dead and hundreds injured (Le et al. 2019;Mori et al. 2019). Jebi, which made landfall in west Japan ( Fig. 1), caused serious disruption and damage in areas around Osaka, which is the second economic centre in Japan after Tokyo. As a result of Jebi, the normal activities at the Kansai international airport were halted for a few days (FDMA 2018) resulting in serious economic loss. The total economic losses in Japan, associated with Jebi, were estimated to be US 12.6 billion dollars, making Jebi the third costliest typhoon in the Japanese history. With a one-minute maximum sustained wind speed of 79 m/s (154 knots), Jebi was classified as a Category-5 hurricane in the SSHWS system. The minimum central pressure reached by Jebi was 915 hPa on 27 August 2016 (Takabatake et al. 2018). However, Jebi then weakened to a Category-3 hurricane during landfall in east Japan (Takabatake et al. 2018). Figure 8a presents sea level pressure (SLP) time series at Kobe and Osaka (Kansai region) indicating minimum SLPs of 959 and 962 hPa, respectively. In comparison with the Lionrock AP records (Fig. 3a), the records of Jebi are sharper and deeper (Fig. 8a). A satellite image of Jebi approaching Japan is shown in Fig. 8b; the size of the typhoon system at that time (3 September 2018) was approximately 450 km.
Significant damage on the east coast of Japan caused by Jebi was strongly associated with coastal floods. A detailed field survey description of the damage was given by Takabatake et al. (2018), Mori et al. (2019) and Le et al. (2019). Some examples of hazardous consequences of the typhoon are shown in Fig. 9: damaged road guardrails Nadakuroiwa (Fig. 9a), sea wall railings broken on the south coast of Hyogo Prefecture (Fig. 9b). Disrupted riprap revetment and wave-driven debris at the coast are presented in Figs. 9c and 9d, respectively.

Observed sea level oscillations
For analysis of sea level oscillations on the coast of Japan caused by the 2018 typhoon Jebi, we used eight stations, five on the south-eastern coast of Honshu Island (Kansai region)-Tannowa, Kumano, Sumoto, Kushimoto and Gobo-and three on the coast of Shikoku Island-Komatsujima, Murotomisaki and Tosashimizu (Fig. 1). The data with 15-s sampling were provided by the JMA. The entire procedure for data preparation and examination was the same as for the Lionrock data (Sect. 2.2). We de-tided the original series and used a 3-h Kaiser-Bessel window (cf. Thomson and Emery 2014) to separate the storm surge from high-frequency storm seiches (meteotsunamis). The results of analyses are shown in Figs. 10 and 11, and the estimated parameters are presented in Table 2. In general, sea level oscillations produced by typhoon Jebi were much stronger than those generated by typhoon Lionrock: the total observed sea levels were 99-288 cm (mean 151 cm), in comparison with 33-120 cm (mean 66 cm) for Lionrock; the Jebi surge heights were 46-173 cm (mean 102 cm), while for Lionrock they were only 12-35 cm (mean 25 cm). Both the maximum total height of 288 cm and the maximum surge height of 173 cm were recorded at Gobo, a small city located in the Wakayama Prefecture. This station is located exactly along the track of typhoon Jebi (Fig. 1). Substantial surges of more than 1 m also occurred at four other stations ( Table 2).
The 2018 typhoon induced intense seiches with maximum amplitudes of 38-131 cm (mean 66 cm) and maximum wave heights of 69-257 cm (mean 122 cm). However, in contrast to the 2016 event, at most stations (at six of the eight) the surge heights were larger Nevertheless, the seiche activity on the coasts of the Kansai region (Honshu) and Shikoku Island during the 2018 typhoon Jebi passage was quite high; seiche amplitudes varied from 38 to 131 cm (mean 66 cm), the maximum trough-to-crest seiche heights were from 69 to 257 cm (mean 122 cm). The latter value-a mean seiche height for eight stations of 122 cm-makes this event extraordinary. The maximum seiche height of 257 cm was observed at Gobo, i.e. at the same station where the maximum surge was recorded. This height of 257 cm is one of the largest meteotsunami heights ever instrumentally measured anywhere in the world (cf. Rabinovich 2020; Gusiakov 2020). The total height above MSL at this station was 288 cm; this means that to the maximum surge height of 173 cm, seiches added 115 cm more, which is an enormous contribution! Similarly to what was observed for the Lionrock event, the total observed height was approximately 11% smaller than the sum of the maximum surge height and the maximum seiche amplitude. The largest The dominant periods of the observed seiches for the 2018 event, based on inspection of the records and f-t diagrams (Sect. 3.3), were from 1 to 60 min (Table 2). In comparison with those observed in 2016 (Table 1), they were more polychromatic: oscillations with 2-3 periods were observed at each station. The exact periods of the recorded oscillations vary from one site to another, but in general periods of 3-15 min prevail.

Time-frequency analysis
The frequency properties/exact periods of the recorded 2018 seiche oscillations and their evolution during the typhoon's passage can be effectively examined based on the corresponding f--t diagrams. The entire procedure is the same as for the 2016 oscillations (Sect. 2-4). The quality of the 2018 data was significantly higher than of the 2016 data, and this enabled us to construct respective diagrams for all examined records (Fig. 12a,  b). For seven records, we used the same 3-day time segments: 3-6 September 2018. The exception was Kumano: the available original record at this station started on 4 September Fig. 11 The high-pass filtered (with a 3-h Kaiser-Bessel window) records of sea level oscillations at eight stations on the coast of Japan showing storm seiches (meteotsunamis) caused by the 2018 typhoon Jebi. The purple and red numbers denote the maximum amplitudes and trough-to-crest wave heights, respectively (see Figs. 10 and 11); therefore, the respective f-t diagram started 25.5 h later than for the other stations (Fig. 12b).
Despite the evident differences between oscillations (in particular, in wave heights and periods) observed at each site, there was a remarkable common feature: the strong broadband amplification of these oscillations during the typhoon's passage. A 7-h band of highly energetic oscillations in the f-t diagrams of all stations corresponds roughly to the mid day of 4 September 2018 and includes the entire frequency band of 1-60 cph (periods of 60 to 1 min). This evidence of a broad-band external (atmospheric) forcing, producing a broadband sea level response. At the same time, there are specific frequencies with pronounced energy peaks: 12 cph (5 min) at Gobo, 20 and 7 cph (3 and 9 min) at Kushimoto, 24 and 7.5 cph (2.5 and 8 min) at Murotomisaki, 15 and 2.4 cph (4 and 25 min) at Tosashimizu (Fig. 12a), 2 cph (30 min) at Komatsujima, 4.5 and 1.2 cph (13 and 50 min) at Tannowa, 12 cph (5 min) at Sumoto and 20 and 6 cph (3 and 10 min) at Kumano. As can be seen, almost all these peaks correspond to the definite frequency bands of increased energy that appear to be associated with the resonant (eigen) frequencies of the respective sites (they are evident in both background and event segments of the f-t diagrams). This means that the typhoon had generated sea level oscillations in a wide range of frequencies, but the strongest response occurs at the resonant frequencies. The entire character of seiche generation by the propagating 2018 typhoon was comparable to that for the 2016 typhoon Lionrock (Fig. 7). Similar to the 2016 event, the broad-band oscillations decay quickly after the 2018 typhoon's passage, while narrow-band oscillations at resonant frequencies keep reverberating. In general, we may conclude that local topography (eigen properties of specific harbours and bays) plays a critical role in the generation of the extreme seiches, which can be considered as meteotsunamis. In this regard, these conclusions are in good agreement with those by Honda et al. (1908) and Nomitsu (1935), who were the first to show that the resonant properties of specific sites are central to the formation of tsunamis and tsunami-like waves on the coast of Japan.

Combined effect of surges and seiches and associated currents
Both storm surges and meteotsunamis are damaging natural hazards, but combined together they are becoming much more destructive. This is a typical situation for typhoons and hurricanes impacting the coastal zones of Japan, USA, Mexico, Caribbean countries, India and some other regions. However, there are very few studies trying to examine the contribution of this cumulative effect (cf. Roeber and Bricker 2015;Bajo et al. 2019). In the present study, we used the sea level data from two intense typhoons, Lionrock and Jebi, that struck Japan in August 2016 and September 2018, respectively. These typhoons induced both storm surges and intense seiches; the two types of sea level oscillations merged together, roughly doubled the hazardous effect, and produced devastating floods along the east coast of Japan. Figure 13 illustrates the relative contribution of storm surges and seiches (meteotsunamis) into the absolute sea level heights during the 2016 and 2018 events. The area of inner (magenta) circles is proportional to the storm surge heights, the blue rings indicate the seiche supplement, and the areas of the entire circles are proportional to the total observed heights shown in Figs. 5 and 10 and indicated in Tables 1 and 2. The relative impact of surges and seiches is different for the two events; in particular, it is seen that for the 2016 event the seiches prevailed, while for the 2018 event the surge contribution was larger. But in general, both components were important, strongly amplifying the resulting flood. It can be seen that for Lionrock the surge heights remain approximately the same for all stations along the track of the typhoon (e.g. Ayukawa and Ofunato), as well as for stations far from the track (e.g. Erimo and Misaki-gyoko). However, for typhoon Jebi, the stations located directly along the track (e.g. Gobo and Murotomisaki) had much higher surges than those located off the track (e.g. Kumano). This could be attributed to the spans of the corresponding typhoons and associated distributions of atmospheric pressure: the respective cloud systems indicate that Lionrock was approximately 900 km in diameter (Fig. 3b), whereas Jebi's diameter was around 450 km (Fig. 8b). Therefore, Lionrock affected a much wider region in comparison to Jebi. However, Jebi had lower air pressures (compare Figs. 3a and 8a) and, consequently, created a substantially stronger surge.
Storm-generated seiches/meteotsunamis, however, are perilous not only because they supplement significant additional height to storm surges and enlarge the flood. What is even more important are the extreme currents that accompany seiche oscillations. Current speeds arising from wave-like motions entering into a bay or harbour are directly proportional to the wave height and inversely proportional to the wave period. The typical meteotsunami period of ~ 15 min for the fundamental (Helmholtz) eigen-mode is approximately 1 3 50 times smaller than the predominant 12.5-h semidiurnal (SD) tidal period. This means that seiche/meteotsunami currents will be fifty times stronger than tidal currents caused by a tide of the same height. For example, if tidal currents on the shelf or in narrows are 10-15 cm/s, then the meteotsunami currents can be ~ 5-7.5 m/s (10-15 knots). Typical durations of storm surges are even longer than tidal periods; therefore, storm generated seiche currents are approximately two orders faster than those associated with the respective storm surge.
Recently, tsunami generated currents began to attract much attention. Dengler et al. (2008), Borrero et al. (2015), Wilson et al. (2013), and Admire et al. (2014) reported significant damage at maritime facilities and harbour infrastructure produced by tsunami currents in ports of California and New Zealand caused by the 2006 Kuril, 2010 Chile and 2011 Tohoku trans-oceanic tsunamis. In particular, in November 2006, the harbour at Crescent City (California) was strongly affected by tsunami currents resulting from the M w 8.3 earthquake in the Kuril Islands; the total damage was almost $20 million (Dengler et al. 2008;Wilson et al. 2013). After the 2010 Chile earthquake tsunami currents of up to 16 knots (8 m/s) caused severe damage to docks in San Diego, Catalina Island, Ventura, and Santa Cruz (California) (Wilson et al. 2013). Lynett et al. (2014) and Borrero et al. (2015) presented many other historical examples of destructive tsunami currents in various regions of the world oceans. Quite often intense tsunami currents are observed in far-field ports, where no flood inducing water levels was recorded; these currents can appear many hours after the initial tsunami arrival and last for two-three days.
Strong tsunami-induced currents and associated eddies and swirls repeatedly occur along the coast of Japan. Thus, Horikawa (1961), following the 1960 Chilean tsunami, reported such currents in Susaki Port on Shikoku Island and at some other sites in Japan. At Hachinohe (the north-eastern coast of Hokkaido Island) 10-knot currents broke moorings, while at Miyako, located on the same coast (Fig. 1), currents of 5-8 knots impeded vessel traffic. The 2011 Tohoku tsunami provided a large set of observational and instrumental data, including the most vivid example of rotational (eddy) currents in the Port of Oarai, located on the east coast of Honshu about 100 km north from Tokyo (Lynett et al. 2012) Inazu et al. (2018. presented a wide collection of ship navigation records indicating that tsunami currents during the 2011 tsunami were observed in numerous bays and harbours along the Pacific coast of Japan. Hence, highly intense currents are a common feature of tsunami events in Japan, and we may assume the same situation for meteotsunami currents. There are several known events when relatively small-height meteotsunamis were accompanied by hazardous currents. In particular, on 27 June 2003 a meteotsunami on the west coast of Croatia affected Stari Grad Bay (Hvar Island) and Mali Ston Bay (near the isthmus of the Pelješac Peninsula). The meteotsunami, which are known locally as "šćiga", caused flooding in the seafront of Stari Grad and swift currents at Mali Ston, resulting in strong damage to shellfish farms in the corresponding bay, despite the fact that that the sea level oscillations at this site were only 40 cm (Vilibić et al. 2004). According to Pattiaratchi and Wijeratne (2015), the meteotsunami of 17 August 2014 in Perth, Australia, produced strong currents of > 1 m/s (2 knots), which resulted in the breakage of moorings and damage to ships inside the Fremantle Port. A tragic event occurred on the Warren Dunes State Park beach (Lake Michigan) on 4 July 2003: a moderate-height (~ 0.3 m) meteotsunami generated strong offshore currents that drowned seven people (Linares et al. 2019). Similar disastrous accidents were reported in this region in the past (Great Lakes Current Incident Database; http://www. misea grant .umich .edu/dcd/dcdse arch.php). Several meteotsunami events in British Columbia (Canada) have been identified not because of significant sea level oscillations, but because of anomalous extreme flows in coastal zones reported by eyewitnesses; such currents can cause damage to yachts, boats and port infrastructure (Thomson et al. 2009;Rabinovich et al. 2020).
Thus, storm seiches superimposed on storm surges not only increase the flooding effect, but can also produce extreme currents and major associated damage. Lynett et al. (2014) provided a special study to correlate tsunami current speed to the damage introduced: There were no current measurements at the examined sites during the 2016 and 2018 events or, at least, they are unknown to the authors of the present study. However, the corresponding speeds may be roughly estimated based on seiche amplitudes and water depths within the harbours. According to Sorensen (2010), the maximum current speed, V max , can be evaluated as: where g is the gravity acceleration, H is the water depth, and A is the zero-to-crest maximum amplitude of the seiche. The values of A were taken from Tables 1 and 2, while the water depth (H) at the respective harbours was taken from the following website: https :// webap p.navio nics.com. The estimated current speeds for the corresponding sites are presented in Table 3 and shown in Fig. 14. The 2018 seiche amplitudes (A) were higher than those in 2016; consequently, the estimated "Jebi currents" were stronger than the "Lionrock currents". We note that Eq. (1) implies that the current speed is directly correlated with the seiche amplitude and is inversely related to the water depth in the harbour.
The speed values in Table 3 should be considered as very approximate: in certain port/ harbour narrows the currents can be even faster. Also, as indicated by Borrero et al. (2015), large eddies and whirlpools are often generated in ports and harbours during tsunamis and these features lead to intense, dynamic and persistent currents, which can be much stronger than estimated by Eq. (1). However, even these preliminary estimates demonstrate that during the 2016 event currents at Erimo, Hachijojima and probably Miyakejima, were strong enough to create damage to anchored boats and cargo vessels and to substantially affect port facilities. The meteotsunami currents caused by typhoon Jebi were substantially stronger; it appears that at Gobo, Kushimoto and Murotomisaki they were of ~ 4-6 knots, while at other examined sites they were 2-3 knots. The reported destruction of port facilities and boats during the 2018 event was probably produced by these currents.

Discussion and perspectives
The seiche oscillations and associated currents generated by 2018 typhoon Jebi were much more intense than those generated by 2016 typhoon Lionrock. Three factors could potentially be the reason for this difference: (a) the stronger forcing of the 2018 event; (b) the more effective parameters of the 2018 typhoon for seiche (meteotsunami) generation; (c) specific resonant properties of the sites affected by the 2018 typhoon. In the following text, these factors are discussed one by one: (a) As evident in a comparison of Figs. 8a and 3a, typhoon Jebi was stronger and deeper than typhoon Lionrock; the minimum air pressure at the east coast of Japan was < 959 hPa for Jebi and ~ 970 hPa for Lionrock. What is especially important: Jebi was more compact and "sharper" than Lionrock (Figs. 8a and 3a); this means that the spatial AP gradients for Jebi were significantly larger. Consequently, Jebi was characterized by extremely fast winds; up to 46.5 m/s (gusts of more than 55 m/s) in the Kansai region, while the maximum Lionrock winds over Japan were up to 35 m/s. It is natural that a stronger atmospheric event could produce a stronger sea level response. (b) The effectivity of meteotsunami/seiche generation strongly depends on parameters of the propagating atmospheric disturbances (Rabinovich and Monserrat 1998). The most important tsunamigenic air pressure (AP) parameters are the speed and direction of the disturbance's propagation (cf. Orlić et al. 2010;Rabinovich et al. 2020). The most effective generation occurs during the Proudman resonance (cf. Hibiya and Kajiura 1982;Monserrat et al 2006;Šepić et al 2015), when the speed of the AP disturbances, U, is approximately equal to the longwave speed of ocean waves, c = √ gH , i.e. when U ≈ c and the Froude number, Fr = U/c ≈ 1.0 (Šepić et al. 2015), and the AP directed  towards the bay/harbour entrance (Rabinovich and Monserrat 1998). However, some other parameters are also important, in particular, the frequency composition of the corresponding disturbance relative to the eigen frequencies of the respective site (Rabinovich and Monserrat 1998). A remarkable example of resonant generation of extreme seiches was described by Mitsuta and Yoshizumi (1968). An extraordinarily strong typhoon (No. 6618, Cora) with a minimum AP of 918 hPa propagated over Miyakojima Island (the Ryukyu Islands, Japan). An analogue barograph at the Miyakojima Weather Station recorded regular AP oscillations with a monochromatic period of ~ 50 min and a double amplitude of up to 15 hPa overlapping the general V-shaped pressure change. These oscillations induced simultaneous sea level oscillations of the same period with a maximum trough-to-crest wave height of 115 cm, i.e. almost 8 times stronger than the AP oscillations (according to the inverse barometer response). We do not know the exact reason for anomalously strong seiche oscillations during 2018 typhoon Jebi, but we can assume that certain parameters of this typhoon, including its speed, direction and frequency content, were favourable for intense seiche generation. (c) As was mentioned above, there are specific sites around the world where meteotsunamis/extreme seiches occur regularly and achieve large heights: Ciutadella Harbour (Menorca Island, Spain), Vela Luka (Korčula Island, Croatia), Mazara del Vallo (Sicily, Italy), Longkou Harbour (China) and Nagasaki Bay (Japan) (Rabinovich and Monserrat 1996;Monserrat et al. 2006;Rabinovich 2020). At these sites, even moderate atmospheric disturbances can produce a strong sea level response. Šepić and Rabinovich (2014) introduced the "generation coefficient" parameter to characterize the efficiency of meteotsunami generation at various sites: where h j max is the maximum recorded trough-to-crest height of generated sea level oscillations at jth station and ΔP j max is the height of the respective air pressure disturbance. This enabled Šepić and Rabinovich (2014) to identify "hot spots" in the area of the Great Lakes and on the northeastern Atlantic coast of the USA, i.e. sites with enlarged sea level response and expected meteotsunamis: Port Calumet (Illinois); Atlantic City, Sandy Hook and Bergen Point (all in New Jersey). Similarly, Dusek et al. (2019Dusek et al. ( ) examined 10-20-year data series (1996Dusek et al. ( -2017 from tide gauge stations along the US East Coast, Caribbean and Bermuda and found that large meteotsunamis regularly occur at a few sites: Atlantic City, Cape Hatteras, Providence, and Port Canaveral. Because of certain combinations of shelf topography and coastal geometry these sites on the Atlantic coast of the USA are hot spots. The question is: Could the sites on the coast of Japan, where the strongest seiches were recorded during the 2016 and 2018 events, also be considered as "hot spots"? Two events are not enough to answer this question, especially taking into account that the 2016 and 2018 typhoons affected different parts of the Japanese coast and the corresponding seiche oscillations were recorded at different stations. Regrettably, we do not have high-resolution air-pressure data, similar to those used by Šepić and Rabinovich (2014) to estimate the "generation coefficient" for various sites on the East Coast of the USA. Therefore, we cannot estimate this coefficient for 18 sites examined in our paper for the 2016 and 2018 events, but we can assume that the sites with maximum recorded high-frequency oscillations (meteotsunamis) during these two events are "hot". Specific topographic features, in particular the funnelling effect of coastal topography (see Figs. 1b and 1c), can significantly amplify incoming waves. An additional important factor responsible for generation of destructive meteotsunamis is the Proudman resonance, U ≈ c. However, both typhoons were approaching the coast of Japan with relatively slow speed, U = 15-20 km/h (Fig. 1a), while the shelf depths in the areas of the typhoon landfalls were rather large (Figs. 1b,c); therefore, this effect could not be responsible for extreme seiches during the 2016 and 2018 events.
A natural approach would be to compare these two events with some historical events, described in the literature. Unfortunately, there are too few papers on analysis of recorded seiche heights in Japan. The most comprehensive study of this question was done by Nakano and Unoki (1962), i.e. almost 60 years ago. They examined multiyear series of seiche oscillations (known in Japan as "secondary undulations of tides") at 45 tide gauge stations located in various harbours and bays of the Japanese coast, including five stations that were examined in the present paper: Kushiro, Ayukawa and Chichijima (2016), Kushimoto and Tannowa (2018). The maximum wave heights recorded at these five stations are given in Table 4.
The data presented in Table 4 enable us to make certain conclusions: • Stations Kushiro, Chichijima and Tannowa cannot definitely be considered as "hot spots"; maximum historical wave heights of atmospherically induced seiches at these stations were significantly smaller than those measured during the 2016 and 2018 typhoons and, in fact, less than 54 cm. • Stations Kushimoto and Ayukawa potentially can be considered as hot spots: seiche heights of more than 1 m had been recorded several times at these stations in the past and also during the 2016 or 2018 events.
(2) Typhoons Lionrock (2016) and Jebi (2018) and associated extreme seiches/meteotsunamis should be regarded as exceptional events; maximum seiche heights at three stations were substantially higher than 10-20-year historical maxima, and two more stations they were of the same order.
It is interesting to compare maximum wave heights recorded during the 2016 and 2018 events with historical maxima for all stations and records examined by Nakano and Unoki (1962). Our data are related to only two specific events, while the historical data presented in that study are based on analysis of 45 stations with continuous sea level observations from 8 to 27 years (between 1930 and 1956). The pivot Table 5 includes all events of atmospheric origin with wave heights more than 120 cm. The maximum wave heights listed by Nakano and Unoki (1962), 169 and 156 cm, were recorded at station Mera located at the southernmost end of the Boso Peninsula approximately 90 km south from Tokyo. Both events (1943 and 1939) were associated with propagating typhoons. Two strong events (145 and 140 cm), initiated by a travelling anticyclone and low-pressure system over the Sea of Japan, respectively, were observed at Miyako (north-eastern coast of Honshu Island, see Fig. 1). Intense seiches with a wave height of 143 cm at Kushimoto were generated by typhoon Agnes (1952).
These were the five strongest events listed in Nakano and Unoki's table. However, both 257-cm seiches at Gobo induced by typhoon Jebi (2018) and 171-cm seiches at Erimo caused by typhoon Lionrock were stronger than these five historical events and all other events described by Nakano and Unoki (1962). Moreover, wave heights of 138 cm at Kushimoto and 137 cm at Kumano, both initiated by typhoon Jebi, are in Lines 8 and 9 in our table (Table 5). All these clearly demonstrate that typhoons Lionrock and Jebi, Table 4 Maximum recorded historical non-seismic trough-to-crest seiche heights at five stations examined in the present study, according to the analyses of Nakano and Unoki (1962) The observation periods at various sites are indicated. Maximum seiche heights recorded at these stations during 2016 typhoon Lionrock or 2018 typhoon Jebi are given in brackets following the historical maximum heights in column "Wave height (cm)" *The only data indicated by Nakano and Unoki (1962) are for the 1952 Kamchatka earthquake (95 cm) and the 1952 Tokachi earthquake (54 cm). This means that the strongest non-seismic seiches at this station for the observation period (1947)(1948)(1949)(1950)(1951)(1952)(1953)(1954)(1955)(1956)  Stationary front Stationary front Typhoon Babs Table 5 Maximum recorded seiche heights (ranked) on the coast of Japan based on historical data from Nakano and Unoki (1962)  Typhoon Tess especially the second one, were extraordinary events. Also we assume that sites listed in Table 5, in particular, Gobo, Erimo, Mera, Miayko, Kushimoto, Kumano, Sakai, etc., may be considered as "hot spots" for the coast of Japan; it appears that strong, atmospherically generated seiches, i.e. meteotsunamis, are a specific feature of these sites. The data presented in Table 5 show that seiches/meteotsunamis generated by typhoons affecting Japan are much stronger than those caused by hurricanes propagating along or across the East USA Coast and the Gulf of Mexico (Olabarrieta et al. 2017;Dusek et al. 2019 Providence; at all other sites the predicted 50-yr heights are significantly smaller. The reason for this difference is not clear: Either typhoons affecting the main coast of Japan are stronger than hurricanes impacting the Atlantic coast of the USA, or the local topographic conditions on the Japanese coast are more favourable to produce strong highfrequency sea level oscillations than those on the East Coast of the USA. It is interesting that some of the events listed in Table 5, in particular Miyako, 1947 (145 cm); Tomie, 1942 (130 cm) and 1940 (126 cm), occurred during anticyclones, i.e. when there were a high pressure and good weather. It is evident that the corresponding seiches had a principally different generation mechanism than those generated by hurricanes, but probably similar to the extreme seiches typical for the Mediterranean Sea (cf. Šepić et al. 2015) and may be considered as "good weather meteotsunamis" (Rabinovich 2020).
It should be emphasized that the information from Nakano and Unoki (1962) presented in Table 5 stops on 1956. Between 1956 and 2016, i.e. in 60 years, there were probably about 150 typhoons affecting the coasts of Japan including such catastrophic typhoons as Typhoon Vera, also known as "Isewan Typhoon" (September 1959), that killed over 5000 people and produced a 3.9 m surge in the Port of Nagoya (https ://wiki2 .org/en/Typho on_Vera). So, these statistics have to be considered as very preliminary. Nevertheless, it gives some important information about the character of atmospherically induced seiches in this region and helps to specify some hot spots.
The Gobo event on 4 September 2018, with a maximum recorded wave height of 257 cm (Table 5), appears to be one of the strongest recorded meteotsunamis generated by a typhoon or hurricane. Moreover, recently Gusiakov (2020) prepared a list of the 20 largest validated meteotsunamis observed in the world in 1978-2020 and that list includes the Gobo event (as No. 20). Certainly, such events as Gobo can produce severe damage to boats and port infrastructure, as well as a serious threat to human lives.

Conclusions
Two devastating typhoons, Lionrock in August 2016 and Jebi in September 2018, destructively affected the eastern coast of Japan and produced extreme sea level variations and coastal floods. The field surveys in the impacted regions showed that the multiple deaths and extensive floods in these regions were caused by the combined effect of low-frequency sea level rise (storm surges) and intensive high-frequency tsunami-like seiches (meteotsunamis). In the present study, we used the data from ten coastal tide gauges for the 2016 event and eight tide gauges for the 2018 event to examine the properties of the sea level oscillations in the affected regions, to estimate the relative contribution of these two sea level components and to evaluate their statistical characteristics. Our main findings are the following: • For the 2016 event, we found that the surge heights were from 12 cm (Hanasaki) to 35 cm (Erimo and Hakodate); the mean contribution of the storm surge into the total observed sea level heights was ~ 39%; the seiche/meteotsunami amplitudes were from 22 cm (Hanasaki) to 92 cm (Erimo), and they contributed ~ 61% of the total height. • For the 2018 event, storm surges were significantly stronger, from 46 cm (Kumano) to 170 cm (Gobo), while seiche amplitudes were from 38 cm (Komatsujima) to 130 cm (Gobo); their mean relative inputs were 67% and 33%, respectively. • The contribution of storm seiches in coastal floods had previously been underestimated, but results of the present study demonstrate that they can play the principal role. Combined together, the surge and tsunami-like seiches create prominent damaging effects. During the 2016 event, the maximum total flood height of 120 cm was at Erimo. The cumulative effect of surge and seiches (meteotsunamis) was exceptionally strong during the 2018 event: the maximum recorded surge + seiche wave heights were 288 cm at Gobo and 178 cm at Murotomisaki; at the other six stations examined they were also high-from 99 to 135 cm. • What is especially important, the extreme currents that produce severe damage to port facilities, anchored boats and coastal infrastructure, are specifically associated with the storm seiches/meteotsunamis. Our rough estimates indicate that the current speeds during the Lionrock event were up to 3 knots and during Jebi were more than 5 knots. • The most important characteristic of recorded meteotsunamis is their trough-to-crest heights. During the 2016 event these heights at three stations were > 1 m: 171 cm at Erimo, 109 cm at Hachijojima and 102 cm at Ayukawa. The 2018 event was stronger; maximum meteotsunami wave heights were 257 cm at Gobo, 138 cm at Kushimoto, 137 cm at Kumano and 128 cm at Murotomisaki. The 2018 Gobo height of 257 cm is much larger than historical non-seismic seiche maxima of 140-169 cm for the Pacific coast of Japan estimated by Nakano and Unoki (1962)