Long wave generation and coastal amplification due to propagating atmospheric pressure disturbances

Meteotsunamis are long waves generated by displacement of a water body due to atmospheric pressure disturbances that have similar spatial and temporal characteristics to landslide tsunamis. NAMI DANCE that solves the nonlinear shallow water equations is a widely used numerical model to simulate tsunami waves generated by seismic origin. Several validation studies showed that it is highly capable of representing the generation, propagation and nearshore amplification processes of tsunami waves, including inundation at complex topography and basin resonance. The new module of NAMI DANCE that uses the atmospheric pressure and wind forcing as the other inputs to simulate meteotsunami events is developed. In this paper, the analytical solution for the generation of ocean waves due to the propagating atmospheric pressure disturbance is obtained. The new version of the code called NAMI DANCE SUITE is validated by comparing its results with those from analytical solutions on the flat bathymetry. It is also shown that the governing equations for long wave generation by atmospheric pressure disturbances in narrow bays and channels can be written similar to the 1D case studied for tsunami generation and how it is integrated into the numerical model. The analytical solution of the linear shallow water model is defined, and results are compared with numerical solutions. A rectangular shaped flat bathymetry is used as the test domain to model the generation and propagation of ocean waves and the development of Proudman resonance due to moving atmospheric pressure disturbances. The simulation results with different ratios of pressure speed to ocean wave speed (Froude numbers) considering sub-critical, critical and super-critical conditions are presented. Fairly well agreements between analytical solutions and numerical solutions are obtained. Additionally, basins with triangular (lateral) and stepwise shelf (longitudinal) cross sections on different slopes are tested. The amplitudes of generated waves at different time steps in each simulation are presented with discussions considering the channel characteristics. These simulations present the capability of NAMI DANCE SUITE to model the effects of bathymetric conditions such as shelf slope and local bathymetry on wave amplification due to moving atmospheric pressure disturbances.


Introduction
Tsunamis are the destructive long waves generated by various sources under the sea surface (e.g., earthquakes, submarine landslides and/or volcanic eruptions) or above the sea surface (e.g., meteorological events, asteroid impacts). While tsunamis generated by earthquakes and landslides have been studied immensely, there is less focus on research related to meteotsunamis (tsunamis of meteorological origin) as they are considered as rare and underrated phenomenon (Pattiaratchi and Wijeratne 2015;Rabinovich 2020). Meteotsunamis are long waves generated by a displacement of a water body due to atmospheric pressure disturbances and often have local names such as "Rissaga" (Balearic Islands), "Marubbio" (Sicily), "Milghuba" (Malta), "Abiki" (Nagasaki Bay) and "Seebär" (Baltic Sea) (Monserrat and Rabinovich 2006).
Long ocean waves of meteorological origin have been observed and described in the scientific literature at the following locations around the world; USA East Coast (Churchill et al. 1995;Sallenger et al. 1995;Pasquet and Vilibić 2013;Lipa et al. 2014), Gulf of Mexico (Paxton and Sobien 1998), The Great Lakes (Ewing et al. 1954;Donn 1959;Bechle and Wu 2014), Atlantic Ocean (Mercer 2002;Candella 2009;Dragani et al. 2009), Adriatic Sea (Vučetić et al. 2009;Šepić et al. 2012), Mediterranean (Airy 1878; Monserrat et al. 1991;Rabinovich and Monserrat 1996;Vilibić et al. 2008), the Aegean Sea (Papadopoulos 1993), Black Sea (Vilibić et al. 2010), East China Sea (Hibiya and Kajiura 1982;Tanaka 2010), Sea of Japan (Park 1986) and Yellow Sea (Wang et al. 1987;Cho et al. 2013). Moreover, Pattiaratchi and Wijeratne (2015) listed meteotsunami events around the world according to maximum wave heights, and recently Rabinovich (2020) overviewed 49 known events to illustrate the origin of these waves and their correlation to atmospheric disturbances. He also classified these meteotsunami events based on good or bad weather conditions as well as their location of impact. One or a few solitary waves are observed on open coasts such as along straight beaches, whereas more intense harbor oscillations are observed in embayments.
General characteristics of meteotsunamis are similar to those of tsunamis of tectonic or landslide origin. They have similar spatial scales with landslide tsunamis due to being local events, whereas temporal scales of several minutes to hours are common for all of them. These tsunami waves may cause devastating destructions in the coastal zone (Monserrat et al. 2006). The difference is in the mechanisms of their generation. For meteotsunamis, sources are spatial and temporal pressure distributions, atmospheric gravity waves and squall lines (Rabinovich 2020). However, in order to be noticeable, meteotsunamis require amplification as the equilibrium response to atmospheric pressure variations (inverted barometer principle) can only generate a few centimeters of water level change (Romero et al. 2019).
There have been several studies focused on the mechanisms of meteotsunami generation and propagation. Some researchers also used analytical and numerical models to analyze the mechanism of ocean wave generation by atmospheric pressure disturbances under idealized conditions. Vilibić (2008) simulated the Proudman resonance caused by sinusoidal air pressure disturbances. He also discussed the similarities between meteotsunamis and landslide tsunamis as both mechanisms generate forced waves and free waves, which then propagate along the ocean surface. A series of analytical studies on the trap and refraction of ocean waves induced by moving atmospheric pressure over a step bottom and continental shelf were conducted by Vennell (2007Vennell ( , 2010, Thiebaut and Vennell (2011) and Melinand (2015). Niu and Zhou (2015) studied the shape of water surface oscillations under a moving low-pressure system with different speeds by adopting a nonlinear shallow water model. They used a moving low pressure with Gaussian distribution in an unbounded flat-bottom water area and analyzed the effects of speed, central pressure drop and spatial scales on the water elevation. Choi and Seo (2017) investigated the ocean waves induced by a low-pressure system moving over a slope towards shore based on an idealized problem from the landing process of a typhoon or a hurricane. They also investigated the runup and development of edge wave patterns generated by a moving a circular shape atmospheric pressure disturbance across a straight shoreline on a sloping beach. A Gaussian shape pressure drop of 20 hPa is used at the center in the circular area of a 200 km radius. Various speeds of pressure disturbance are selected up to 50 m/s and tested. The change in the maximum runup is presented according to different moving speeds for the given atmospheric pressure disturbance. While wave runup and rundown are repeated at the shoreline, the generation and propagation of edge waves in the alongshore direction due to wave refraction are also observed.
In this study, the atmospheric pressure disturbance acts on the water surface and generates free and forced waves, which are amplified due to the Proudman resonance and shoaling. Similarly, Chen and Niu (2018) focused on the occurrence of Proudman resonance when the speed of moving pressure disturbance is close to the local shallow water wave velocity (celerity) as the disturbance moves to the shore over a slope. They performed a numerical investigation on the generation and propagation of forced waves induced by a circular shape atmospheric pressure drop. The wave pattern evolves from a concentric-circle type into a triangular type with the increase in the Froude number on the slope, which can be defined as Fr = V/c; the ratio of the atmospheric pressure speed, V, to the celerity of long ocean waves c = √ gh where g is the gravitational acceleration and h is the water depth. The effects of the characteristic parameters of the circular pressure disturbances and slope gradient are discussed. They found that it is not always possible to observe the significant peak of the maximum water elevation before the landing of pressure disturbances. The extremely high runup is observed when the forced wave hits the shore. In another recent study, Niu and Chen (2020) state that energy is transferred from the atmosphere to the ocean during the response of the sea to a disturbance, while the forced wave induced by the disturbance gradually grows and finally reaches a quasi-steady state. The time required to excite the forced wave induced by pressure disturbance as it reaches a steady state is quantified based on the nonlinear shallow water wave model by Niu and Chen (2020).
Many meteotsunami events are reproduced by numerical modeling with different degrees of accuracy (Rabinovich et al. 1999;Liu et al. 2003a;Vilibić et al. 2004Vilibić et al. ,2008Dragani, 2007;Orlić et al. 2010;Tanaka, 2010;An et al. 2012;Bechle and Wu, 2014;Šepić et al. 2018;Heidarzadeh et al. 2020). Romero et al. (2019) state that a simple shallow water model may suffice to represent the dynamics and amplification of long ocean waves provided the atmospheric coupling and resonant mechanisms are accounted for. Much of the research on meteotsunami modeling uses available numerical ocean models that utilize Navier-Stokes equations with hydrostatic and Boussinesq approximations. Romero et al. (2019) use shallow water equations after the inclusion of the appropriate atmospheric forcing term, and they include a drag force in the momentum equation to integrate bottom roughness and seabed frictional resistance. Šepić et al. (2018) use an ocean model code developed by Dr. Isaac Fine (Institute of Ocean Sciences, Sidney, BC) to simulate Odessa meteotsunami based on depth-integrated shallow water equations with Coriolis term and all nonlinear terms, except bottom friction, omitted. In this code, the time derivative of air pressure is a principal input parameter and it is independent of other model outputs and parameters. Bubalo et al. (2018) simulate Chrystal and Proudman resonances in a closed rectangular basin with two finite element (ADCIRC and SCHISM) models and one finite difference (ROMS) ocean model. They state that even though all three models are solving the same equations with similar approximations, there is a substantial difference in timestepping integration technique between them, and the Courant-Friedrichs-Levy (CFL) condition becomes important when these models are used to describe Proudman resonance as accurately as possible. They conclude that the errors that may arise from choosing an inadequate integration time step for all areas of the grid may be very significant and affect the quality of the results of Proudman resonance in the open sea, which will, therefore, affect the reproduction in coastal areas. Denamiel et al. (2019) discuss the challenges in meteotsunami modeling, stating that resolution and bathymetry used in ocean models, as well as meteorological forcing, are not accurate enough to truly recreate the generation, propagation, and amplification of these events. Following the experiences of many researchers (Monserrat et al. 2006;Tanaka, 2010;Horvath and Vilibić, 2014), they conclude that while the main challenge is reproducing the atmospheric physics responsible for creating the atmospheric disturbances, the development of reliable atmospheric-ocean coupled models for the short-term forecast is also important. They recommend that for an ocean model to capture the offshore meteotsunami dynamics, the atmospheric forcing should be imposed every minute to capture the speed and amplitude of pressure disturbances, and the spatial resolution should allow for an accurate representation of the continental shelf. However, Denamiel et al. (2019) conclude that spatial discretization and lack of precise bathymetric data limits the accuracy of the present state-of-the-art ocean numerical models to represent the ocean and complex geomorphology of the coast responsible for the meteotsunami wave generation and amplification (Proudman and topographic resonances). Their discussion focuses on the operation of a regional meteotsunami system developed for the Adriatic region and aims to work in real-time for early warning and forecasting purposes.
For the ocean, Proudman and topographic resonances are the main generation and amplification processes of long oceanic waves (Monserrat et al. 2006;Vennell, 2010). Due to spatial discretization and lack of precise bathymetric data, it is unlikely that present state-of-the-art ocean numerical models can accurately represent the ocean and coast complex geomorphology responsible for the meteotsunami wave generation and amplification, which is particularly true for the middle and southern parts of the eastern Adriatic. However, given these overwhelming limitations, the AdriSC Meteotsunami Forecast component is showing skills in fair reproduction of meteotsunami events, which might be used for operational forecasting. Bubalo et al. (2019) attempted to improve the numerical simulations by including the flooding and drying option of a finite element model (ADCIRC) for Adriatic since many of the previous studies either had cut-off depths set at several meters or did not utilize flooding and drying option. They also state that flooding and drying are not frequently used in meteotsunami studies, contrary to storm surge and tsunami research. They show that the inclusion of flooding and drying leads to a much more accurate reproduction of the event and grid refinement is very important as well. They recommend the use of a flooding and drying algorithm for modeling meteotsunami events.
This paper presents the new module of NAMI DANCE tsunami numerical code that solves the nonlinear shallow water equations with the air pressure forcing term to simulate meteotsunamis. An ocean model should capture both the offshore generation and resonant amplification of meteotsunamis, and the nearshore amplification of the tsunami-like wave by topography and harbor resonance (Rabinovich 2009). Based on the vast validation of NAMI DANCE model, the nearshore amplification process of tsunami waves has already been satisfied, including harbor resonance. ; Lynett et al. (2017) and Dogan et al. (2021) are some of the works that present the capability of NAMI DANCE model in properly representing the nearshore amplification process due to topography applied in various tsunami events. The performance of NAMI DANCE model to represent harbor resonance and the amplification process of long waves are presented in Sozdinler et al. (2015); Yalciner et al. (2015) and Sogut and Yalciner (2019). The previous applications of NAMI DANCE model include the tsunami waves generated by earthquakes and landslides. The additional module to solve generation, propagation and coastal amplification of the ocean waves due to spatial and temporal change of pressure and wind fields is introduced in this paper. The cases simulated in this paper provides validation of NAMI DANCE by comparing the numerical results with analytical solutions and show that it is capable of modeling the generation, propagation, and coastal amplification of ocean waves and Proudman resonance due to moving atmospheric disturbances.

Resonance by coupling of ocean wave motion and atmospheric pressure forcing
When the air pressure decreases gradually, the water level in the sea region rises in the order of 1 cm with respect to 1 hPa pressure decrease. The water level changes can be amplified by the resonance when the ocean wave speed coincides with the speed of the atmospheric pressure disturbance. The three most important resonance mechanisms are (i) Proudman resonance (Proudman, 1929), (ii) Greenspan resonance (Greenspan, 1956) and (iii) shelf resonance (Monserrat et al. 2006). Pattiaratchi and Wijeratne (2015) specified that Proudman resonance is the main cause of meteotsunamis happening all over the world. Proudman resonance occurs when the velocity (V) of the atmospheric pressure movement coincides with the celerity (c) of the ocean wave, which denotes V = c. The changes in the water surface level due to air pressure disturbances depend on the Froude number (Fr). Šepić et al. (2015) stated that when Fr ≈ 1.0, during movement of the ocean waves, atmospheric energy is absorbed and, therefore, strong amplification occurs.
Greenspan resonance occurs where the air pressure velocity in alongshore direction ( V l ) is close to the propagation phase speed of j-th mode of edge waves ( c j ). The generated edge waves amplify due to Greenspan resonance in relation to the amplitude and the propagation speed of the pressure disturbance and a beach slope.
The atmospheric pressure disturbance causes the initial change in water level, and water oscillations develop due to the shelf geometry such as length, water depth, topographic funneling, and bed slope. Shelf resonance occurs when the period or wavelength of atmospherically induced ocean waves and waves in the continental shelf region coincide.
The internal resonance mechanism has another important role in wave amplification inside enclosed basins and harbors. When the natural resonance frequency of a semienclosed basin coincides with the frequency of incoming wave or forcing to the basin, then strong amplification can be observed inside the basin. A short cut numerical method for the determination of periods of free oscillations for basins with irregular geometry and bathymetry is presented in Yalciner and Pelinovsky (2007). Bailey et al. (2014) stated that in long, narrow, and shallow bays, stronger oscillations occur while the system tends to have a low rate of energy dissipation.

Wave generation and propagation due to atmospheric pressure disturbances
The energy transfer in the oceans by ocean waves is generally categorized according to the energy sources such as wind, fault break, submarine or subaerial landslides, pressure disturbance, or solar and planetary attractions. Air pressure disturbances cause water level changes and hence generate long-period ocean waves. High atmospheric pressure causes a static water level drop and contrariwise, low atmospheric pressure causes a static water level rise. Due to the different water level changes at high pressure and low pressure regions, the water level differences cause the wave generation in the sea, and consequently, propagation and coastal amplification in the ocean area.
On the other hand, bathymetry is an important factor that may cause amplification of the long waves generated by atmospheric pressure disturbance. The speed of pressure disturbance (V) and wave celerity (c) are the two important parameters that govern the characteristics of the generated long waves. The speed of the ocean waves is controlled by the water depth and it may be faster or slower or equal to the speed of pressure disturbance. Therefore, the shape and propagation of the generated ocean waves differ according to the different speeds of pressure disturbances as well as the magnitude and gradient of the pressure disturbance. In this study, the effects of the speed of pressure disturbance on the generation, propagation and amplification of the ocean waves in the basins with different bathymetries are investigated.
Simple analytical solutions for long wave generation in a basin of constant depth (Proudman resonance) and a basin with a linear sloping bottom (Greenspan resonance) are given in their original papers (Proudman 1929;Greenspan 1956) and then reproduced in many papers. This paper aims to define the analytical solution of the linear shallow water model and compare the results with the numerical solution of the linear and nonlinear shallow water model. Thereafter, analytical and numerical models are applied to certain test cases. This helps to interpret the verification of the numerical model and the wave generation and propagation mechanism due to the speed of pressure disturbances over different types of basins.

Analytical solution
The governing equations for tsunami wave generation by atmospheric disturbances in narrow bays and channels can be written similar to the 1D case studied for tsunami landslide generation (Tinti et al. 2001;Tinti and Bortolucci 2000;Liu et al. 2003b;Pelinovsky, 2006;Didenkulova et al. 2010;Didenkulova and Pelinovsky, 2013): where H(x, t) is the water depth along the channel axis, u(x, t) is the flow velocity averaged over a cross section, g is the gravity acceleration, S(H, x, t) is the variable cross section of the bay filled by water, p atm (x,t) is atmospheric pressure, ρ is the density of the water, x is coordinate, and t is time. Resonant amplification of tsunami waves induced by moving landslides is studied in (Didenkulova et al. 2011;Didenkulova and Pelinovsky 2012).
Generated ocean waves and the atmospheric disturbance are assumed to be uniform over the cross section.
In the important case of the quasi-parabolic channel, its cross section is described by where m is an arbitrary positive constant. In particular, m = 1 describes the triangular channel, m = 2 describes the parabolic channel, and m → ∞ corresponds to the rectangular channel. As a result, we may compute a function S(H) For such channels, we have Eqs. (1) and (2) for wave components: The set of analytical solutions can be obtained if the atmospheric disturbance is weak. In this case, we may use the linear version of these equations where (x, t) = H(x, t) − h is the vertical displacement of the water surface, and h is an unperturbed water depth along the longitudinal axis.
This system can be reduced to the wave equation where where c is the wave celerity in the channel. It is important to note that Eq. (9) has a universal form for waves in channels of any quasi-parabolic shape; the channel shape is described by the coefficient m in Eq. (10). Namely, this equation is studied for describing the Proudman resonance (Melinand, 2015). By introducing Equation (9) can be rewritten in the following form In particular, the response on the static atmospheric disturbance is described by the well-known "inverse barometer law." It means that low atmospheric pressure can cause the sea level rise, and high atmospheric pressure leads to depression in the mean sea level. If the pressure is moving, the resonance can occur when the pressure velocity is close to the wave celerity. Then, replacing we may transform Eq. (12) to In this form, Eq. (15) coincides with those used for tsunamis generated by landslides (Tinti et al. 2001;Liu et al. 2003b;Pelinovsky, 2006;Didenkulova et al. 2010;Didenkulova and Pelinovsky, 2013).
Initial conditions for the shallow water system (7-8) are applied to both (x, t) and u (x, t) . If at the initial moment the ocean rests, they can be transformed into conditions for the function (x, t) The linear Eq. (15) with initial conditions (16) is the Froude number. The wave celerity varies with the water depth, and Fr varies with the speeds of pressure disturbance and generated water waves.
The response to the moving disturbance represents three waves. The first term in Eq. (19) describes the forced wave, which propagates together with the atmospheric disturbance. The forced wave can have different polarities depending on the regime of the disturbance motion. Assuming ζ a is positive, the forced wave is positive in the super-critical regime ( Fr > 1 ) and negative in the sub-critical regime ( Fr < 1 ). The second term describes a free wave, which propagates with its own speed c. As the forced wave, it can also be of different polarities: negative in the super-critical regime and positive in the subcritical regime. The last term corresponds to the free wave of positive polarity. It is important to mention that the solution (19) has the same form for waves generated by a landslide moving with constant speed in a basin of constant depth (Tinti et al. 2001;Pelinovsky, 2006), for 1D landslide moving along the non-reflecting bottom profile h ∼ x 4∕3 (Didenkulova et al. 2010) and along narrow bays of variable cross sections (Didenkulova and Pelinovsky, 2013).
For practice, it is useful to compare amplitudes of waves moving in the same direction as an atmospheric disturbance. The ratio of the amplitude (Q) of the free wave (second term) to the amplitude (P) of the forced wave (first term) is The general sectional view of long waves generated by moving constant atmospheric pressure is shown in Fig. 1 for the case Fr < 1.
Formula (21) is correct for any Fr. But for Fr < 1, it is the ratio of the front wave to the rear wave (forced wave is behind the free wave), as shown in Fig. 1. For Fr > 1, the forced wave is in front of the free wave, and therefore, Eq. (21) gives the ratio of the rear wave to the front wave.

Numerical model
Atmospheric pressure disturbances on the water body may induce long waves. The three distinctions are important for the long waves: i) generation of the wave, ii) propagation of the wave in deep water, and iii) propagation and coastal amplification in shallow water and on the shore. Wave propagation in deep water can be explained by the linear shallow water equations. In contrast, in shallow water, wave dynamics is better described by the nonlinear shallow water equations, which also cover the surface and bottom stresses, Coriolis force, and shoaling (Eze et al. 2009). A full potential model or Boussinesq-type equations are even better in shallower regions. Several numerical models were tested and compared in Liu et al. (2008) and Lynett et al. (2017). Some of those numerical codes are COMCOT (Liu et al. 1994;2008), TUNAMI-N2 (Imamura 1996;Imamura et al. 2006;Liu et al. 2008), MOST (Titov and Synolakis, 1998;Liu et al. 2008;Lynett et al. 2017), NAMI DANCE (Liu et al. 2008;Sozdinler et al. 2015;Metin 2016;Lynett et al. 2017;Sogut and Yalciner 2019;Dogan et al. 2021), which use shallow water wave equations, while FUNWAVE (Kirby et al. 1998) and GEOWAVE (Watts et al. 2003) use Boussinesq-type equations. The aforementioned numerical models are developed for the numerical solution of co-seismic and/or landslide generated tsunamis without considering the atmospheric pressure disturbance.
In the present study, the nonlinear shallow water equations containing the air pressure forcing and wind field forcing terms are solved numerically to simulate generation, propagation and coastal amplification of long waves generated by the atmospheric pressure disturbances and wind field forcing. In this direction, NAMI DANCE is upgraded to the new version called NAMI DANCE SUITE with the capability of using the Graphical Processing Unit of the graphic card. It computes the main hydrodynamic parameters of long waves, which are (i) maximum positive amplitude, (ii) maximum flow depth, (iii) maximum current velocity, (iv) maximum momentum flux, (v) maximums hydrodynamic forces, (vi) maximum negative amplitude, (vii) maximum wave arrival time, (viii) initial wave arrival time and (ix) durations of inundation and shoreline withdrawal Sozdinler et al. 2015;.
The set of two-dimensional equations with atmospheric pressure and wind field terms in Cartesian coordinates are given in Eqs. (23)-(25): where t is time; x and y are spatial coordinates in West-East and South-North directions, respectively; is the water surface elevation; w is the water density; air is the air density; D is the water depth, P atm is the atmospheric pressure in Pascal; x and y are the bottom shear stresses; U w10 and V w10 are wind velocities at 10 m elevation in x and y directions; C D is the wind drag coefficient computed by the following equations (Garrat 1977;Powell et al. 2003): where M and N are the discharge fluxes in the x and y directions: In Eq. (28), u and v are the water current velocities in x and y directions (horizontal plane), respectively. The bottom shear stresses are computed using Eqs. (29) and (30). where f b denotes the bottom friction, which is 0.015, and Manning's coefficient n is calculated by Since Coriolis force and friction are not included in the analytical solution, those are not included in the numerical tests for proper comparisons with the analytical solution.
In the numerical model, in addition to standard inputs (bathymetry, initial water surface displacement, and fluxes if needed), the spatial distribution of the barometric pressure at the sea level in Pascal at one-minute intervals during the simulations is employed. As a set of output parameters, NAMI DANCE SUITE computes water surface elevations, magnitude and directions of water current velocities, momentum fluxes, overland flow depths throughout the domain at specified time intervals, and extremums of all these parameters. Additionally, the inundation and withdrawal limits of the shoreline motion; and time histories of these parameters at specified numerical wave gauge locations are computed.

Verification of the model
In order to verify the model to correctly calculate the long wave motion due to the effects of atmospheric pressure disturbances, the test simulations are conducted with a regular shaped flat bathymetry (200 km of distance between 0.9 W and 0.9E longitudes, 195 km of distance between 0.885S and 0.885 N latitudes and 200 m of constant depth). The North and South boundaries of the basin are used as solid impermeable boundaries in order to prevent the escape of waves through these borders. The East and West borders are open (outgoing) boundaries. The 16 km wide high pressure disturbance (102,000 Pa) propagating with a constant velocity from West to East is used in all simulations. The standard pressure over the basin is set as 100,000 Pa. The borders of the pressure disturbance are smoothed by the 2 km linear increase from and decrease to 100,000 Pa at the front and rear sides of the disturbance (Fig. 2). The selected pressure disturbance band is an exceptional case. The main reason to select this disturbance is to obtain higher water fluctuations for a better comparison of the simulation results.  Fig. 3.
To verify the accuracy of the numerical model, the simulations are performed for six different Froude numbers. The ratios of the amplitude of the free wave (determined by the wave celerity) to the amplitude of the forced wave (determined by the speed of pressure disturbance) are computed and compared with the analytical predictions (see Table 1). In the numerical simulations for comparison and validation with the analytical solution, the linear wave equations are used. The comparisons show that the majority of the % Error (the ratio of the absolute difference between numerical results and the corresponding analytical results) for each Fr is less than 1%. Furthermore, the same simulations are performed using different grid sizes of 50 m, 100 m, 177 m, 354 m and 708 m for the sensitivity analysis due to the grid size. The numerical results for Fr = 1 indicate that no significant difference occurred with changing grid size in the range of used sizes with a % Error less than 1.5%.

3 4 Numerical tests for different types of basins
The investigation of the long wave generation and propagation due to the moving pressure disturbance with a constant speed is important to understand considering the possible wave amplification at shallow regions in relation to bathymetric and morphologic conditions in basins and/or marine environment. In this direction, three different basins are selected for applications with different cross sections in lateral or longitudinal directions with respect to directions of pressure disturbance and ocean wave; (i) triangular lateral cross section, (ii) stepwise shelf from 200 to 20 m of depth on 1:10 slope in the longitudinal direction, (iii) stepwise shelf from 200 to 20 m of depth on 1:700 slope in the longitudinal direction. Additionally, similar stepwise shelf basins on slopes varying from 1:100 to 1:600 are simulated to expand the numerical tests over different slopes.

Simulations in the basin with triangular cross section
The basin with a triangular cross section has a 200 m water depth at the center along the direction of pressure disturbance with 1:444 side slopes. The top and sectional views of the basin are given in Fig. 4. According to Eq. 10, the representative wave celerity of this For three different speeds of pressure disturbance, the maximum wave amplification is observed for the case of 44 m/s speed of pressure disturbance at 200 m depth, the deepest part (along the central axis) of the basin. A similar wave amplification was observed for the case of 31 m/s speed of the pressure disturbance, which coincides with the corresponding speed of the free wave at about 94 m water depth. In the case of V = 22 m/s, the extrema of positive and negative amplitudes are less compared to the other two cases observed at a depth of ~ 49 m.  The results for the two different basins are given in the following sections. In addition, six different slopes varying from 1:100 to 1:600 are also numerically tested for 31 m/s speed of pressure disturbance to elaborate on the effect of slope in section (ii) on wave amplification. The results of these simulations show that as the slope becomes milder, there is an increase in the extrema of the positive and negative wave amplitudes during the 120 min simulation (the row indicated with a star symbol in Table 2). The amplification in the amplitudes of the depression (front part) and the elevation (rear part) waves is also higher for a milder slope but observed at different time steps of the simulation depending on the location of the slope (bold amplitude values in Table 2).

Stepwise shelf bathymetry with 1:10 shelf slope
The top and sectional views of the stepwise shelf bathymetry with 1:10 shelf slope are shown in Fig. 6. The simulations are performed using three different speeds of pressure disturbance (22 m/s, 31 m/s and 44 m/s) propagating over this basin from the left (West) to the right (East).
The distributions of water surface (sea states) and sectional views at zero latitude at selected time steps (t = 10 min, 20 min, 40 min, 60 min and 90 min) for three different speeds of pressure disturbance moving over the stepwise shelf bathymetry with 1:10 shelf slope are presented in Fig. 7.
Negative and positive amplitudes of the front and rear parts of the generated wave at different time steps for three different speeds of pressure disturbance are given in Table 3.
The simulations for the case of V = 22 m/s speed of pressure disturbance ( Fig. 7 and Table 2

Stepwise shelf bathymetry with 1:700 shelf slope
The top and sectional views of the stepwise shelf bathymetry with 1:700 shelf slope are shown in Fig. 8. The simulations are performed using three different speeds of the pressure disturbance (22 m/s, 31 m/s and 44 m/s) propagating over this basin from the left (West) to the right (East). The distributions of water surface (sea states) and sectional views at zero latitude at selected time steps (t = 10 min, 20 min, 40 min, 60 min and 90 min) for three different speeds of the pressure disturbance, moving over the stepwise shelf bathymetry with 1:700 shelf slope from 200 to 20 m water depths are presented in Fig. 9. The negative and positive amplitudes of the front and rear parts of the generated waves at different time steps for three different speeds of the pressure disturbance are given in Table 4.
It is seen from Fig. 9 and Table 4 709 m), whereas the positive amplitude shows the highest maximum at around 90 min (3.552 m). In the stepwise shelf bathymetry with a 1:700 slope basin, the maximum wave amplification is observed for the case of V = 44 m/s, which satisfies the Proudman resonance in the deep (left) part of the basin while approaching the shelf slope. In other words, the same speeds of pressure disturbance and the wave celerity in the deep (left) part of the basin caused larger amplification of the wave before arriving at the shelf. Furthermore, a higher wave amplification is observed as the wave reaches the top (end) of the When we compare the results given in Tables 3 and 4, it is observed that the milder shelf slopes can cause higher wave amplification.

Discussions and conclusions
The effects related to the speed of the atmospheric pressure disturbance on the generation and amplification of a meteotsunami in different basins are investigated by using analytical and numerical approaches. In this study, the analytical solution of the wave generation and propagation due to the propagation of the atmospheric pressure disturbance over a flat bottom basin is developed. The numerical model NAMI DANCE solving tsunamis of seismic and landslide origins is upgraded to solving the long-period ocean wave generation and motion due to the spatial and temporal change of the atmospheric pressure and velocity fields. The model is verified by comparing the results with those from analytical solutions on a flat bottom basin. The numerical model is applied to investigate the generation and coastal amplification of the ocean waves in the basins with a triangular lateral cross section and stepwise sloping shelf with different shelf slopes towards shore in the propagation direction. Analytically, the new solution for meteotsunami generation in the water channel of 'power' cross section (including the channel of triangular cross section) is developed.
For verification of the numerical model NAMI DANCE SUITE, a flat bottom basin, which corresponds to the constant wave celerity (c), is selected. Different speeds of pressure disturbance (V = 11 m/s, 22 m/s, 33 m/s, 44 m/s, 55 m/s, 66 m/s) propagating over the 200 m deep flat bottom basin corresponding to different flow conditions (sub-critical, critical and super-critical) are used as the forcing mechanism. Satisfactory agreements between the numerical and analytical results are obtained.
Furthermore, numerical simulations are performed for different water basins: (i) a channel with a triangular cross section, (ii) a stepwise shelf bathymetry with 1:10 shelf slope from 200 to 20 m water depths, and (iii) a stepwise shelf bathymetry with 1:700 shelf slope from 200 to 20 m water depths. For each basin, three different speeds of the pressure disturbance (22 m/s, 31 m/s and 44 m/s) are considered. Additionally, numerical tests for stepwise shelf bathymetry are extended to use different shelf slopes from 1:10 to 1:700.
One of the governing parameters which cause the sea level fluctuations and wave amplification is the relation between the speed of pressure disturbance and the speed of the generated wave in the basin. When the speed of pressure disturbance exceeds the speed of the generated water wave, the generated wave cannot evolve because the forcing mechanism (pressure disturbance) goes faster and leaves the basin much earlier than the generated wave. On the other hand, when the speed of pressure disturbance is slower than the speed of the generated water wave, the water wave evolves to a certain amplitude (higher amplitudes than in the previous case) and propagates faster than the forcing mechanism (pressure disturbance). Finally, when the pressure disturbance propagates with the same speed of the generated water wave, the wave evolves with the pressure disturbance and propagates with the highest amplitude, as the forcing mechanism propagates together with the generated wave and continuously inputs the energy directly to the generated wave. Another important factor that leads to the amplification of meteotsunami generated by the propagating pressure disturbance is the bathymetric conditions of the basins. In order to investigate the effects of bathymetric conditions on wave amplification due to the propagating pressure disturbance, different basins with different bathymetric conditions are considered.
Shown for the basin of triangular cross section, the maximum wave amplification is observed at the water depth, where the speed of the pressure disturbance coincides with the wave celerity.
The shelf slope is another parameter affecting the amplification of the wave. Eight different basins of stepwise shelf bathymetry with slopes ranging from 1:10 to 1:700 (sloping section from 200 to 20 m water depth) are used in the simulations. It is shown that for the same speed of the atmospheric pressure disturbance, the amplification over the stepwise shelf bathymetry with milder (1:700) shelf slope is greater than for the steeper (1:10) one, which indicates that the milder shelf slopes can cause higher wave amplification. This is also confirmed with the computed amplitudes for the other slopes (1:100-1:600).