2. General Equations and Asymptotic Formula for Wave Profile in the Uniform Depth Basin

2.1. The General Equations

[4]  Let us give first very rough arguments showing the possibility to use (1, 2) for a description of tsunami waves in a frame of the piston model over slow varying bottom. The equations (1) are obtained from the following linear equations of the piston model without the assumption of long wave approximation for the velocity potential F(x,z,t) and the elevation of the free surface h(x,t) (see e.g. [Whitmore and Sokolowski, 1996; Pelinovski, 1996; Kowalik and Murty, 1993; Lewis and Adams, 1983; Shokin et al., 1989; Tinti, 1993])

eq003.gif

eq004.gif(3)

eq005.gif

where d is the delta function, and it is assumed that Fequiv 0, hequiv 0 when t< 0.

[5]  The task is to evaluate h(x,t) at large distances |x|>>l from the source, near the wavefronts (where |h| has its maximum values), assuming that l>>H(x) and that the depth H(x) has a small variation at distances of the order l.

2.2. Asymptotic Formula

[6]  In the case of the basin of a uniform depth H0, the wavefronts are the circles |x|=C0 t, where C0=(gH0)1/2 is the velocity of the long waves. We seek the solution for the problem (3) in the form of the Fourier transform with respect to x=(x1,x2) and with Fourier parameter p=(p1,p2). Using in the integrand polar coordinates (r,f) defined by the formulas p1=r/l cosf, p2=r/l sinf we obtain

eq006.gif(4)

where

eq007.gif(5)

and the angle y is defined by the equations

eq008.gif

Below, we consider in detail the case, when the initial vertical displacement of the bottom is like a ridge (or a valley), though the approach developed here is applicable for any other localized disturbance.

[7]  To make final formulas more explicit let us model the initial displacement using the function of the type of the Gaussian exponent with oscillations:

eq009.gif

eq010.gif(6)

where a1, a2, 1 < b1, b2 < 3, q, c are dimensionless parameters and V0 is the parameter of the dimension of length.

[8]  Substituting for V(y) from (6) in (5) we obtain

eq011.gif(7)

where

eq012.gif(8)

Now we derive the equation (9) under the main assumption |x|>>1. We estimate the integral with respect to q by the stationary phase method. After some calculations we have the following asymptotic formula for h(x,t) (see [Borovikov and Kelbert, 1996; Dobrokhotov et al., 1993; Berry, 2005]):

eq013.gif(9)

[9]  From (7) and (8) one can conclude, reminding the above assumption on the values of the parameter b, that the main contribution to the integral (9) corresponds to the values of r in the interval 0 < r< 3. Moreover, since |x|/l>>1 integral (9) is not small only in the case when the expression in the big parentheses in the exponent of the integrand is small. Further, since H0/l<<1, the functionseq014.gif and eq015.gif in the integrand can be expanded in Taylor series in the powers of rH0 / l:

eq016.gif(10)

One can see that the integral (9) gets its largest values near the circle |x|=C0 t, i.e. in a neighborhood of the wavefront. Then, if we take in series (10) the first terms only, we obtain from (9) the formulae for h(x,t) in the long waves approximation. From (9) and (10) we can derive roughly the following requirement for the correctness of the long wave approximation

eq017.gif(11)

Taking H0=4 km, l=50 km we find from (11) that for the typical conditions of the ocean, equation (1) can be used instead of the general equations (3) in the piston model in the case of a basin of uniform depth. We note that the same conclusion can be done for initial disturbances of more general types than those described by (6) as well as in the case of a basin with a non-uniform slowly varying depth. According to this reason we will consider here the problem in the long wave approximation (1, 2) instead of the general equations (3).

[10]  Thus, from (9) an asymptotic formula for h(x,t) at any instant t and observation point x located near wavefront and satisfying requirement |x|>>l and (11) becomes

eq018.gif(12)

where y is the angle between the vector x and the axis Ox1,

eq019.gif(13)

In (13) the function

V(r,y)

is given by (5) for an initial disturbance of a general type and by (7) for the disturbance (6) considered here. In the last case the integral in (13) is evaluated analytically and F(z,y) becomes

eq020.gif(14)

where

eq021.gif

and 1F1(...) are the Kummer's confluent hypergeometric functions. The function F can be expressed also in terms of parabolic cylinder functions D<font face="Symbol">n (see [Dobrokhotov et al., 1991]), we use here the form given by Mathematica.

2006ES000215-fig01
Figure 1
[11]  One can see from (12) that the function F(z,y) defines the structure of wave profiles near the wave front (i.e. for values of z that are small or are of the order of several units). In turn, the function F(z,y) depends on the form of the initial bottom disturbance (through function V in (13)) and the variable z. The dependence Re F(z,y) on z
2006ES000215-fig02
Figure 2
is presented in Figures 1 and 2 for different values of the angle  y.

[12]  Figure 1 shows that in the case b1 neq b2 (non axially symmetric initial bottom disturbance) the wave profile has different forms at different values of the angle y, while Figure 2 shows that the difference becomes very strong if a1 neq a2. So the wave profile crucially depend on the form of initial disturbance determined by (6).

[13]  In the case of the initial disturbance of general type, the function F(z,y) can be calculated from (13) and (5) numerically. The calculation is simplified by the fact that the important part of the integral (13) corresponds to the values of the variable of integration r of the order of several units.

[14]  Keeping in mind that the structures of the asymptotic formulae for h(x,t) are similar in the cases of the basins of uniform and non uniform depth (see Sec 3), we give the following comments on the formulae (12) using terminology of the WKB theory.

[15]  The function F(x,t) in (12) can be called the phase because of two facts. First, from (12) and (13) it is clear that |h(x,t)| has a maximum when F(x,t)=0 and decays rapidly with increasing |F(x,t)|. Second, outside of the neighborhood of origin x=0 the function F(x,t) is the action or a solution of 2-D Hamilton-Jacobi equation

eq022.gif(15)

corresponding to (1) in the case of the basin of the uniform depth H0.

[16]  The integration of the last equation is connected with the Hamiltonian system (see e.g. [Maslov and Fedoiuk, 1981]); for (15) it is determined by the following 2-D Hamiltonian H=C0|p|, where p=(p1,p2) is a momentum. Thus the corresponding Hamiltonian system has the form

eq023.gif(16)

Consider the one-parametric family of trajectories of this system satisfying the following initial conditions

eq024.gif(17)

where

eq025.gif(18)

and the angle yin [0,2p] is a parameter. These trajectories are (vector-columns) X(t,y)=(C0 t cos y, C0 t sin y)T, P(t,y)= n(y). The family of trajectories X(t, y), P(t,y) go out from the point x=0 with a unit momentum p= n(y). The projection of the trajectory x=X(t, y) on the plane mathbbrx.gif is called "ray". Here the rays are the straight lines starting from the point x=0 and arriving at the instant t to the point x1=C0 t cos y, x2=C0 t sin y. At every instant t>0 the curve formed by the end of the rays is called the "wavefront". Here the wavefronts are the circles |x|=C0 t. One can see from (12) that at the wavefronts the phase F=0. It should be noted also that the points on the front are parameterized by the parameter yin (0,2p] which was introduced in (17) as the initial condition for the impulse p(t). So, y in (12) can be considered with the same point of view.

[17]  In the case of a basin with a non-uniform depth, asymptotic formulas for h(x,t) have a structure similar to that in (12). But the rays are no more straight lines, the Hamiltonian is more complicated, the factor (l / |x|)1/2 has to be replaced by other one and sometimes the power in the exponent e-ip/4 has to be multiplied by an integer number m which has a deep topological meaning since it is possible, under certain condition, to evaluate m as the Morse index of the trajectory.

[18]  Anyway, it is clear that the algorithm for calculating h(x,t) from the (12) is faster than the finite difference method since for finding the value of the amplitude of the tsunami front passing in the point x at time t it is not necessary to integrate over all the space but one need just to identify the trajectory starting from the point where the initial perturbation has taken place and arriving to the point x at time t. In other words this method of integration is "local" in the sense that it is concentrated on the trajectories of an Hamiltonian system. This Hamiltonian system is not difficult to solve numerically in the case of a basin of non-uniform depth.


RJES

Powered by TeXWeb (Win32, v.2.0).