Vincent Duchêne

Tatsuo Iguchi

A mathematical analysis of the Kakinuma model for interfacial gravity waves. Part II: Justification as a shallow water approximation

. In this present paper, we show that the Kakinuma model is a higher order shallow water approximation to the full model for interfacial gravity waves with an error of order O(δ 4N +2 1

) in the sense of consistency, where δ 1 and δ 2 are shallowness parameters, which are the ratios of the mean depths of the upper and the lower layers to the typical horizontal wavelength, respectively, and N is, roughly speaking, the size of the Kakinuma model and can be taken an arbitrarily large number. Moreover, under a hypothesis of the existence of the solution to the full model with a uniform bound, a rigorous justification of the Kakinuma model is proved by giving an error estimate between the solution to the Kakinuma model and that of the full model. An error estimate between the Hamiltonian of the Kakinuma model and that of the full model is also provided.

.

, we analyzed the initial value problem of the Kakinuma model when the approximated velocity potentials are defined by (1.2). We found that the Kakinuma model has a stability regime which can be expressed as

0 on the interface, where H 1 := h 1 -ζ and H 2 := h 2 + ζ -b are the depths of the upper and the lower layers, P app 1 and P app 2

Introduction

We will consider the motion of the interfacial gravity waves at the interface between two layers of immiscible fluids in (n + 1)-dimensional Euclidean space. Let t be the time, x = (x 1 , . . . , x n ) the horizontal spatial coordinates, and z the vertical spatial coordinate. We assume that the layers are infinite in the horizontal directions, bounded from above by a flat rigid-lid, and from below by a time-independent variable topography. The interface, the rigid-lid, and the bottom are represented as z = ζ(x, t), z = h 1 , and z = -h 2 + b(x), respectively, where ζ = ζ(x, t) is the elevation of the interface, h 1 and h 2 are mean depths of the upper and lower layers, and b = b(x) represents the bottom topography. See Figure 1.1. We assume that the fluids in the upper and the lower layers are both incompressible and inviscid fluids with constant densities ρ 1 and ρ 2 , respectively, and that the flows are both irrotational. Then, the motion of the fluids is described by the velocity potentials Φ 1 (x, z, t) and Φ 2 (x, z, t) and the pressures P 1 (x, z, t) and P 2 (x, z, t) in the upper and the lower layers. We recall the governing equations, referred as the full model for interfacial gravity waves, in Section 2 below. Generalizing the work of J. C. Luke [START_REF] Luke | A variational principle for a fluid with a free surface[END_REF], these equations can be obtained as the Euler-Lagrange equations associated with the Lagrangian density L (Φ 1 , Φ 2 , ζ) given by the vertical integral of the pressure in both water regions. Building on this variational structure, T. Kakinuma [START_REF] Kakinuma | [END_REF][START_REF] Kakinuma | A set of fully nonlinear equations for surface and internal gravity waves[END_REF][START_REF] Kakinuma | A nonlinear numerical model for surface and internal waves shoaling on a permeable beach[END_REF] proposed and studied numerically the model obtained as the Euler-Lagrange equations for an approximated Lagrangian density, L (Φ app 1 , Φ app 2 , ζ), where

(1.1) Φ app ℓ (x, z, t) = N ℓ i=0
Z ℓ,i (z; hℓ (x))ϕ ℓ,i (x, t)

for ℓ = 1, 2, and {Z 1,i } and {Z 2,i } are appropriate function systems in the vertical coordinate z and may depend on h1 (x) and h2 (x), respectively, which are the depths of the upper and

x z ζ(x, t) h 1 h 2 ρ 1 ρ 2 Figure 1
.1: Internal gravity waves the lower layers in the rest state, whereas ϕ ℓ = (ϕ ℓ,0 , ϕ ℓ,1 , . . . , ϕ ℓ,N ℓ ) T , ℓ = 1, 2, are unknown variables. This yields a coupled system of equations for ϕ 1 , ϕ 2 , and ζ, depending on the function systems {Z 1,i } and {Z 2,i }, which we named Kakinuma model. Note that in our setting of the problem we have h1 (x) = h 1 and h2 (x) = h 2 -b(x). In this work we study the Kakinuma model obtained when the approximate velocity potentials are defined by (1.2)

             Φ app 1 (x, z, t) := N i=0
(-z + h 1 ) 2i ϕ 1,i (x, t), Φ app 2 (x, z, t) :=

N * i=0 (z + h 2 -b(x)) p i ϕ 2,i (x, t),
where N, N * , and p 0 , p 1 , . . . , p N * are nonnegative integers satisfying 0 = p 0 < p 1 < • • • < p N * . Specifically, we show that the Kakinuma model obtained through the approximated potentials (1.2) with (H1) N * = N and p i = 2i (i = 0, 1, . . . , N ) in the case of the flat bottom b(x) ≡ 0, (H2) N * = 2N and p i = i (i = 0, 1, . . . , 2N ) in the case with general bottom topographies, provides a higher order shallow water approximation to the full model for interfacial gravity waves in the strongly nonlinear regime. The choice of the function systems as well as N, N * , and p 0 , p 1 , . . . , p N * is discussed and motivated later on.

Comparison with surface gravity waves The Kakinuma model is an extension to interfacial gravity waves of the so-called Isobe-Kakinuma model for surface gravity waves, that is, water waves, in which Luke's Lagrangian density L Luke (Φ, ζ), where ζ is the surface elevation and Φ is the velocity potential of the water, is approximated by a density L app (ϕ, ζ) = L Luke (Φ app , ζ), where

(1.3) Φ app (x, z, t) = N i=0 Z i (z; b(x))ϕ i (x, t).
The Isobe-Kakinuma model was first proposed by M. Isobe [START_REF] Isobe | A proposal on a nonlinear gentle slope wave equation[END_REF][START_REF] Isobe | Time-dependent mild-slope equations for random waves[END_REF] and then applied by T. Kakinuma to simulate numerically the water waves. Recently, this model was analyzed from a mathematical point of view when the function system {Z i } is a set of polynomials in z: Z i (z; b(x)) = (z + h -b(x)) p i with integers p i satisfying 0 = p 0 < p 1 < • • • < p N . The initial value problem was analyzed by Y. Murakami and T. Iguchi [START_REF] Murakami | Solvability of the initial value problem to a model system for water waves[END_REF] in a special case and by R. Nemoto and T. Iguchi [START_REF] Nemoto | Solvability of the initial value problem to the Isobe-Kakinuma model for water waves[END_REF] in the general case. The hypersurface t = 0 in the space-time R n × R is characteristic for the Isobe-Kakinuma model in the sense that the operator acting on time derivatives of the unknowns has a non-trivial kernel. As a consequence, one needs to impose some compatibility conditions on the initial data for the existence of the solution. Under these compatibility conditions, the non-cavitation condition, and a Rayleigh-Taylor type condition -∂ z P app ≥ c 0 > 0 on the water surface, where P app is an approximate pressure in the Isobe-Kakinuma model calculated from Bernoulli's equation, they showed the well-posedness of the initial value problem in Sobolev spaces locally in time. Moreover, T. Iguchi [START_REF] Iguchi | Isobe-Kakinuma model for water waves as a higher order shallow water approximation[END_REF][START_REF] Iguchi | A mathematical justification of the Isobe-Kakinuma model for water waves with and without bottom topography[END_REF] showed that under the choice of the function system

(1.4) Z i (z; b(x)) = (z + h) 2i
in the case of the flat bottom, (z + h -b(x)) i in the case of a variable bottom, the Isobe-Kakinuma model is a higher order shallow water approximation for the water wave problem in the strongly nonlinear regime. Furthermore, V. Duchêne and T. Iguchi [START_REF] Duchêne | A Hamiltonian structure of the Isobe-Kakinuma model for water waves[END_REF] showed that the Isobe-Kakinuma model also enjoys a Hamiltonian structure analogous to the one exhibited by V. E. Zakharov [START_REF] Zakharov | Stability of periodic waves of finite amplitude on the surface of a deep fluid[END_REF] on the full water wave problem and that the Hamiltonian of the Isobe-Kakinuma model is a higher order shallow water approximation to the one of the full water wave problem.

Our aim in the present paper and the companion paper [START_REF] Duchêne | A mathematical analysis of the Kakinuma model for interfacial gravity waves. Part I: Structures and well-posedness[END_REF] is to extend these results on surface gravity waves to the framework of interfacial gravity waves. With respect to surface gravity waves, our interfacial gravity waves framework brings two additional difficulties. The first one is that, due to the rigid-lid assumption, the full system for interfacial gravity waves described in Section 2 features only one evolution equation for the two velocity potentials, and a constraint associated with the fixed fluid domain. From a physical perspective, the unknown velocity potential at the interface may be interpreted as a Lagrange multiplier associated with the constraint. A second important difference between water waves and interfacial gravity waves is that the latter suffer from Kelvin-Helmholtz instabilities. As a consequence the initial value problem of the full model for interfacial gravity waves is ill-posed in Sobolev spaces; see T. Iguchi, N. Tanaka, and A. Tani [START_REF] Iguchi | On the two-phase free boundary problem for twodimensional water waves[END_REF], V. Kamotski and G. Lebeau [START_REF] Kamotski | On 2D Rayleigh-Taylor instabilities[END_REF]. This raises the question of the validity of any model for interfacial gravity waves. A partial answer is offered by the work of D. Lannes [START_REF] Lannes | A stability criterion for two-fluid interfaces and applications[END_REF], which proves the existence and uniqueness of solutions over large time intervals in the presence of interfacial tension. While interfacial tension effects is not expected to be the relevant regularization mechanism for the propagation of waves between, for instance, fresh and salted water, the key observation is that physical systems allow the propagation of waves with large amplitude and long wavelengths provided that some mechanism tames Kelvin-Helmholtz instabilities acting on the high-frequency component of the flow. This description is consistent with the fact that the initial value problem of the bi-layer shallow water system for the propagation of interfacial gravity waves in the hydrostatic framework is well-posed in Sobolev spaces under some hyperbolicity condition describing the absence of lowfrequency Kelvin-Helmholtz instabilities, as proved by D. Bresch and M. Renardy [START_REF] Bresch | Well-posedness of two-layer shallow water flow between two horizontal rigid plates[END_REF]. Let us mention however that such a property is not automatic for higher order shallow water models. Specifically, we note that the Miyata-Choi-Camassa model derived by M. Miyata [START_REF] Miyata | An internal solitary wave of large amplitude[END_REF] and W. Choi and R. Camassa [START_REF] Choi | Fully nonlinear internal waves in a two-fluid system[END_REF] and which can be regarded as a two-layer generalization of the are approximate pressures of the fluids in the upper and the lower layers, α 1 and α 2 are positive constants depending only on N and on p 0 , p 1 , . . . , p N * , respectively. This is a generalization of the aforementioned Rayleigh-Taylor type condition for the Isobe-Kakinuma model. It is worth noticing that, consistently with the expectation that the Kakinuma model is a higher order model for the full system for interfacial gravity waves and that the latter suffers from Kelvin-Helmholtz instabilities, the constants α 1 and α 2 converge to 0 as N and N * go to infinity so that the stability condition becomes more and more stringent as N and N * grow. When N = N * = 0, the Kakinuma model coincides with the aforementioned bi-layer shallow water system, and the stability regime coincides with the hyperbolic domain exhibited in [START_REF] Bresch | Well-posedness of two-layer shallow water flow between two horizontal rigid plates[END_REF]. Moreover, when the motion of the fluids together with the motion of the interface is in the rest state, the above stability condition is reduced to the well-known stable stratification condition (1.6) (ρ 2 -ρ 1 )g > 0.

In [START_REF] Duchêne | A mathematical analysis of the Kakinuma model for interfacial gravity waves. Part I: Structures and well-posedness[END_REF], we showed that under the stability condition (1.5), the non-cavitation assumptions (1.7)

H 1 ≥ c 0 > 0, H 2 ≥ c 0 > 0,
and intrinsic compatibility conditions on the initial data, the initial value problem for the Kakinuma model is well-posed in Sobolev spaces locally in time. We also showed in [START_REF] Duchêne | A mathematical analysis of the Kakinuma model for interfacial gravity waves. Part I: Structures and well-posedness[END_REF] that the Kakinuma model enjoys a Hamiltonian structure analogous to the one exhibited by T. B. Benjamin and T. J. Bridges [START_REF] Benjamin | Reappraisal of the Kelvin-Helmholtz problem. Part 1. Hamiltonian structure[END_REF] on the full model for interfacial gravity waves.

Comparison with other higher order models The Isobe-Kakinuma and the Kakinuma models belong to higher order models for the water waves and for the full interfacial gravity waves, respectively. By this we mean a family of systems of equations parametrized by nonnegative integers describing the order of the system within the family, that is N for the Isobe-Kakinuma model,and whose solutions are expected to approach solutions to the full system as the order increases. Several such models have been introduced in the literature, mostly in the water waves framework, and we will restrict the discussion to water waves in this paragraph.

Based on a Taylor expansion of the Dirichlet-to-Neumann operator at stake in the water waves system with respect to the shape of the domain D. G. Dommermuth and D. K. Yue [START_REF] Dommermuth | A high-order spectral method for the study of nonlinear gravity waves[END_REF], B. J. West et al. [START_REF] West | A new numerical method for surface hydrodynamics[END_REF] and W. Craig and C. Sulem [START_REF] Craig | Numerical simulation of gravity waves[END_REF] have proposed the so-called high order spectral (HOS) models. While these models have been successfully employed in efficient numerical schemes (see recent accounts by J. Wilkening and V. Vasan [START_REF] Wilkening | Comparison of five methods of computing the Dirichlet-Neumann operator for the water wave problem[END_REF], D. P. Nicholls [START_REF] Nicholls | High-order perturbation of surfaces short course: boundary value problems[END_REF] and P. Guyenne [START_REF] Guyenne | HOS simulations of nonlinear water waves in complex media[END_REF]), the equations feature Fourier multipliers which prevent their direct use in situations involving non-trivial geometries such as horizontal boundaries. Moreover, the rigorous justification of HOS models is challenged by well-posedness issues; see the discussion in D. M. Ambrose, J. L. Bona, and D. P. Nicholls [START_REF] Ambrose | On ill-posedness of truncated series models for water waves[END_REF], and V. Duchêne and B. Melinand [START_REF] Duchêne | Rectification of a deep water model for surface gravity waves[END_REF].

A second class of higher order models originate from formal shallow water expansions put forward by J. Boussinesq [START_REF] Boussinesq | Addition au mémoire sur la théorie des ondes et des remous qui se propagent le long d'un canal rectangulaire, etc[END_REF] and J. W. S. Rayleigh [START_REF] Rayleigh | On waves[END_REF]. A systematic derivation procedure has been described by K. O. Friedrichs in the appendix to [START_REF] Stoker | The formation of breakers and bores. The theory of nonlinear wave propagation in shallow water and open channels[END_REF]. Recently, these higher order shallow water models have been described and discussed by Y. Matsuno in [START_REF] Matsuno | Hamiltonian formulation of the extended Green-Naghdi equations[END_REF][START_REF] Matsuno | Hamiltonian structure for two-dimensional extended Green-Naghdi equations[END_REF] and W. Choi in [START_REF] Choi | On Rayleigh expansion for nonlinear long water waves[END_REF][START_REF] Choi | High-order strongly nonlinear long wave approximation and solitary wave solution[END_REF]. The derivation procedure displays formula for approximate velocity potentials under the form (1.3)-(1.4) (in particular, only even powers appear in the flat bottom case), with the important difference that the functions ϕ i (i = 0, . . . , N ) are prescribed through explicit recursion relations. An important consequence of this derivation is that the resulting systems of equations involve only standard differential operators. However the order of the differential operators at stake augments with the order of the system, which renders such models impractical for numerical simulations.

By contrast, the Isobe-Kakinuma model features only differential operators of order at most two acting on the variables ϕ i (i = 0, . . . , N ) which are unknowns of the system. Notice that the size of the system augments with its order, N . However the degrees of freedom do not augment with the order since, as mentioned above, some compatibility conditions must be satisfied. In fact all quantities are uniquely determined by two scalar functions which represent the canonical variables in the Hamiltonian formulation of the water waves system. Let us mention that function systems different from (1.4) have been considered by G. A. Athanassoulis and K. A. Belibassakis [START_REF] Athanassoulis | A consistent coupled-mode theory for the propagation of small-amplitude water waves over variable bathymetry regions[END_REF], G. Klopman, B. van Groesen, and M. W. Dingemans [START_REF] Klopman | A variational approach to Boussinesq modelling of fully nonlinear water waves[END_REF] and C. E. Papoutsellis and G. A. Athanassoulis [START_REF] Papoutsellis | A new efficient hamiltonian approach to the nonlinear water-wave problem over arbitrary bathymetry[END_REF] (see also references therein). While the systems obtained in these works have a similar nature, they are all different. We let the reader refer to V. Duchêne [START_REF] Duchêne | Many Models for Water Waves[END_REF]Chapter D] for an extended discussion and comparison of these models.

The choice of the function systems in (1.4) is motivated by the aforementioned Friedrichs expansion and is essential in the analysis of T. Iguchi [START_REF] Iguchi | Isobe-Kakinuma model for water waves as a higher order shallow water approximation[END_REF][START_REF] Iguchi | A mathematical justification of the Isobe-Kakinuma model for water waves with and without bottom topography[END_REF] proving that the Isobe-Kakinuma model is a higher order shallow water approximation for the water wave problem in the strongly nonlinear regime. We note that one may modify (1.4) by putting all odd and even terms (z + h) i for i = 0, 1, . . . in the case of the flat bottom. However, in that case one needs to use the terms up to order 2N to keep the same precision of the approximation. Therefore, such a choice increases the number of unkonwns and equations by N so that it is undesirable for practical application.

In other words, one can save memories in numerical simulations by using only even terms in the case of the flat bottom. On the contrary, if we put only odd terms (z + h -b(x)) 2i for i = 0, 1, 2, . . . in the case of a non-flat bottom, then the corresponding Isobe-Kakinuma model does not give any good approximation even if we take N a sufficiently large number, because the corresponding approximate velocity potential Φ app cannot approximate the boundary condition on the bottom so well due to the lack of odd order terms (z + h -b(x)) 2i+1 for i = 0, 1, 2, . . .. Following this discussion, the choice of the function systems (1.2) with (H1) or (H2) in our interfacial waves framework is very natural. In particular, the rigid-lid is assumed to be flat so that we do not need to use odd order terms (-z + h 1 ) 2i+1 for i = 0, 1, 2, . . ., in the approximate velocity potential Φ app 1 to obtain a good approximation, because Φ app 1 can approximate the boundary condition on the rigid-lid without such terms.

Description of the results

In the present paper we show that the Kakinuma model obtained through the approximated potentials (1.2) with (H1) N * = N and p i = 2i (i = 0, 1, . . . , N ) in the case of the flat bottom b(x) ≡ 0, (H2) N * = 2N and p i = i (i = 0, 1, . . . , 2N ) in the case with general bottom topographies, provides a higher order shallow water approximation to the full model for interfacial gravity waves in the strongly nonlinear regime. Our results apply to the dimensionless Kakinuma model obtained after suitable rescaling. The system of equations then depend on the positive dimensionless parameters δ 1 and δ 2 which are shallowness parameters related to the upper and the lower layers, respectively, that is, δ ℓ = h ℓ λ (ℓ = 1, 2) with the typical horizontal wavelength λ. The shallow water regime is described through the smallness of the parameters δ 1 and δ 2 . What is more, our results are uniform with respect to parameters satisfying either ρ 2 ≲ ρ 1 < ρ 2 , or ρ 1 ≪ ρ 2 and h 2 ≲ h 1 . We notice that the rigid-lid framework is expected to be invalid in the regime ρ 1 ≪ ρ 2 and h 1 ≪ h 2 which is excluded in this paper; see V. Duchêne [START_REF] Duchêne | Boussinesq/Boussinesq systems for internal waves with a free surface, and the KdV approximation[END_REF].

Our first result extends the result of [START_REF] Duchêne | A mathematical analysis of the Kakinuma model for interfacial gravity waves. Part I: Structures and well-posedness[END_REF] on the well-posedness of the initial value problem by showing that solutions to the dimensionless Kakinuma model are defined on a time interval which does not vanish for arbitrarily small values of δ 1 and δ 2 .

Theorem 1.1 (Long-time well-posedness). Under the (dimensionless) stability condition (1.5), the (dimensionless) non-cavitation assumptions (1.7), and intrinsic compatibility conditions on the initial data, the initial value problem for the Kakinuma model is well-posed in Sobolev spaces on a time interval which is independent of δ 1 ∈ (0, 1] and δ 2 ∈ (0, 1].

While the non-cavitation assumption and the stability condition are automatically satisfied for small initial data and small bottom topography b, an arrangement of nontrivial initial data satisfying the compatibility conditions with suitable bounds is a non-trivial issue, and demands a specific analysis. Proposition 1.2. Initial data satisfying the compatibility conditions and necessary bounds in Theorem 1.1 are uniquely determined (up to an additive constant) by sufficiently regular initial data for the canonical variables of the Hamiltonian structure.

Then, we show that under the special choice of the indices p 0 , p 1 , . . . , p N * as in (H1) or (H2), the dimensionless Kakinuma model is consistent with the full model for interfacial gravity waves with an error of order O(δ 4N +2 1

+ δ 4N +2 2 ).
Theorem 1.3 (Consistency). Assume (H1) or (H2). The solutions to the dimensionless Kakinuma model constructed in Theorem 1.1 produce functions that satisfy approximately the dimensionless full interfacial gravity waves system up to error terms of size O(δ 4N +2 1

+ δ 4N +2 2 ).
Conversely, solutions to the dimensionless full interfacial gravity waves system satisfying suitable uniform bounds produce through Proposition 1.2 functions that satisfy approximately the dimensionless Kakinuma model up to error terms of size O(δ 4N +2 1

+ δ 4N +2 2 ).
In the last result we assume the existence of a solution to the full model with a uniform bound since for general initial data in Sobolev spaces, one cannot expect to construct a solution to the initial value problem, due to the ill-posedness of the problem discussed previously. The same issue arises for the full justification of the Kakinuma model. Theorem 1.4 (Full justification). Assuming the existence of a solution to the dimensionless full interfacial gravity waves system with a uniform bound and satisfying initially the (dimensionless) stability condition (1.5) and (dimensionless) non-cavitation assumptions (1.7), then the Kakinuma model with (H1) or (H2) and appropriate initial data produces an approximate solution with the error estimate

|ζ K (x, t) -ζ IW (x, t)| ≲ δ 4N +2 1 + δ 4N +2 2
on some time interval independent of δ 1 ∈ (0, 1] and δ 2 ∈ (0, 1], where ζ K and ζ IW are solutions to the dimensionless Kakinuma model and to the full model, respectively.

In our last main result, we show that the Hamiltonian structure of the Kakinuma model is a shallow water approximation of the Hamiltonian structure of the full interfacial gravity waves model.

Theorem 1.5 (Hamiltonians). Assume (H1) or (H2). Under appropriate assumptions on the canonical variables (ζ, ϕ), we have

|H K (ζ, ϕ) -H IW (ζ, ϕ)| ≲ δ 4N +2 1 + δ 4N +2 2 ,
where H K and H IW are the Hamiltonians of the dimensionless Kakinuma model and of the dimensionless full interfacial gravity waves model, respectively.

Remark 1.6. The precise statements of our main results are displayed in Section 3. Specifically, Theorem 1.1 corresponds to Theorem 3.1, Proposition 1.2 corresponds to Proposition 3.4, Theorem 1.3 corresponds to Theorem 3.5 and Theorem 3.6 (see also Remark 3.8), Theorem 1.4 corresponds to Theorem 3.9, and Theorem 1.5 corresponds to Theorem 3.10.

Structures of the Kakinuma model

In order to obtain our main results, we exploit several structures of the Kakinuma model. The Kakinuma model can be written compactly as

(1.8)            l 1 (H 1 )∂ t ζ + L 1 (H 1 )ϕ 1 = 0, l 2 (H 2 )∂ t ζ -L 2 (H 2 , b)ϕ 2 = 0, ρ 1 l 1 (H 1 ) • ∂ t ϕ 1 + 1 2 |u 1 | 2 + w 2 1 -ρ 2 l 2 (H 2 ) • ∂ t ϕ 2 + 1 2 |u 2 | 2 + w 2 2 + (ρ 1 -ρ 2 )gζ = 0,
where we denote ) T , and the linear operators L ℓ , and functions u ℓ and w ℓ for ℓ = 1, 2 are defined (after non-dimensionalization) in Section 3. Here we recognize the fact that the hypersurface t = 0 in the space-time R n × R is characteristic for the Kakinuma model, since the system of evolution equations is overdetermined for the variable ζ, and underdetermined for the variables ϕ 1 and ϕ 2 . As a consequence, solutions to the Kakinuma model must satisfy some compatibility conditions. Introducing linear operators L 1,i (i = 0, . . . , N ) acting on φ 1 = (φ 1,0 , . . . , φ 1,N ) T and L 2,i (i = 0, . . . , N * ) acting on φ 2 = (φ 2,0 , . . . , φ 2,N * ) T by

ϕ 1 := (ϕ 1,0 , ϕ 1,1 , . . . , ϕ 1,N ) T , ϕ 2 := (ϕ 2,0 , ϕ 2,1 , . . . , ϕ 2,N * ) T , put l 1 (H 1 ) := (1, H 2 1 , H
                                   L 1,0 (H 1 )φ 1 := N j=0 L 1,0j (H 1 )φ 1,j , L 1,i (H 1 )φ 1 := N j=0 (L 1,ij (H 1 )φ 1,j -H 2i 1 L 1,0j (H 1 )φ 1,j ) for i = 1, 2, . . . , N, L 2,0 (H 2 , b)φ 2 := N * j=0 L 2,0j (H 2 , b)φ 2,j , L 2,i (H 2 , b)φ 2 := N * j=0 (L 2,ij (H 2 , b)φ 2,j -H p i 2 L 2,0j (H 2 , b)φ 2,j ) for i = 1, 2, . . . , N * ,
the necessary conditions can be written simply as (1.9)

     L 1,i (H 1 )ϕ 1 = 0 for i = 1, 2, . . . , N, L 2,i (H 2 , b)ϕ 2 = 0 for i = 1, 2, . . . , N * , L 1,0 (H 1 )ϕ 1 + L 2,0 (H 2 , b)ϕ 2 = 0.
The first two vectorial identities are analogous to the compatibility conditions of the Isobe-Kakinuma model for water waves, while the last identity is specific to the bi-layer framework and is related to the continuity of the normal component of the velocity at the interface.

A first key ingredient of the analysis is the fact that for sufficiently regular functions ζ, b and ϕ 1 (respectively ϕ 2 ), there exists a unique solution ϕ 1 (respectively ϕ 2 ) to the problems (1.10)

l 1 (H 1 ) • ϕ 1 = ϕ 1 , L 1,i (H 1 )ϕ 1 = 0 for i = 1, 2, . . . , N, l 2 (H 2 ) • ϕ 2 = ϕ 2 , L 2,i (H 2 , b)ϕ 2 = 0 for i = 1, 2, . . . , N *
satisfying suitable elliptic estimates. What is more, the well-defined linear operators

Λ (N ) 1 (ζ) : ϕ 1 → L 1,0 (H 1 )ϕ 1 , Λ (N * ) 2 (ζ, b) : ϕ 2 → L 2,0 (H 2 , b)ϕ 2 , are found to approximate the corresponding Dirichlet-to-Neumann maps Λ 1 (ζ) and Λ 2 (ζ, b) defined by Λ 1 (ζ)ϕ 1 := -∂ z Φ 1 + ∇Φ 1 • ∇ζ z=ζ(x,t) , Λ 2 (ζ, b)ϕ 2 := ∂ z Φ 2 -∇Φ 2 • ∇ζ z=ζ(x,t) ,
where Φ 1 and Φ 2 are the unique solutions to Laplace's equations

     ∆Φ 1 + ∂ 2 z Φ 1 = 0 in Ω 1 (t), Φ 1 = ϕ 1 on Γ(t), ∂ z Φ 1 = 0 on Σ 1 , and
     ∆Φ 2 + ∂ 2 z Φ 2 = 0 in Ω 2 (t), Φ 2 = ϕ 2 on Γ(t), ∇Φ 2 • ∇b -∂ z Φ 2 = 0 on Σ 2 ,
where we denote the upper layer, the lower layer, the interface, the rigid-lid, and the bottom at time t by Ω 1 (t), Ω 2 (t), Γ(t), Σ 1 , and Σ 2 , respectively. Specifically it is proved that, under the special choice of the indices p 0 , p 1 , . . . , p N * in (H1) or (H2) and after suitable rescaling, the difference between the dimensionless operators is of size

O(δ 4N +2 1 + δ 4N +2 2
). This analysis, which follows directly from the corresponding analysis for surface waves developed in [START_REF] Iguchi | A mathematical justification of the Isobe-Kakinuma model for water waves with and without bottom topography[END_REF] and scaling arguments, provides the key argument in the proof of the consistency result described in Theorem 1. [START_REF] Benjamin | Reappraisal of the Kelvin-Helmholtz problem. Part 1. Hamiltonian structure[END_REF].

In order to study the Kakinuma model, we also need to analyze the full elliptic system (1.11)

           L 1,i (H 1 )ϕ 1 = f 1,i for i = 1, 2, . . . , N, L 2,i (H 2 , b)ϕ 2 = f 2,i for i = 1, 2, . . . , N * , L 1,0 (H 1 )ϕ 1 + L 2,0 (H 2 , b)ϕ 2 = ∇ • f 3 , -l 1 (H 1 ) • ϕ 1 + l 2 (H 2 ) • ϕ 2 = f 4 ,
for sufficiently regular functions ζ, b and

f 1 = (f 1,1 , . . . f 1,N ) T , f 2 = (f 2,1 , . . . , f 2,N * ) T , f 3 , f 4 .
The ellipticity of the problem relies on the coercivity of the corresponding operators L 1 (H 1 ) and L 2 (H 2 ). The solvability of (1.11) is essential in several directions. Firstly, it provides an alternative consistency result, where solutions to the full interfacial gravity waves system produce approximate solutions to the Kakinuma model but satisfying exactly and not approximately the necessary conditions (1.9). In turn, this provides a crucial ingredient to the full justification of the Kakinuma model described in Theorem 1.4. Furthermore, the arrangement of initial data satisfying the compatibility conditions as stated in Proposition 1.2 amounts to solving (1.11) with f 1 = 0, f 2 = 0, f 3 = 0 and f 4 = ϕ. Similarly, our result on the Hamiltonians H K and H IW described in Theorem 1.5 relies on a comparison of solutions to (1.11) with f 1 = 0, f 2 = 0, f 3 = 0 and f 4 = ϕ and solutions to

                     ∆Φ 1 + ∂ 2 z Φ 1 = 0 in Ω 1 (t), ∆Φ 2 + ∂ 2 z Φ 2 = 0 in Ω 2 (t), ∂ z Φ 1 = 0 on Σ 1 , ∇Φ 2 • ∇b -∂ z Φ 2 = 0 on Σ 2 , (∇Φ 1 • ∇ζ -∂ z Φ 1 ) -(∇Φ 2 • ∇ζ -∂ z Φ 2 ) = 0 on Γ(t), ρ 2 Φ 2 -ρ 1 Φ 1 = ϕ on Γ(t),
thus extending to the interfacial gravity waves framework the analysis in [START_REF] Duchêne | A Hamiltonian structure of the Isobe-Kakinuma model for water waves[END_REF]. Finally, the solvability of (1.11) allows to determine and control time derivatives ∂ t ϕ 1 and ∂ t ϕ 2 of sufficiently regular solutions to the Kakinuma model (1.8) by using the equations obtained when differentiating with respect to time the compatibility conditions (1.9) combined with the last equation of (1.8). This is a crucial ingredient for the analysis of the initial value problem. Another crucial ingredient for the analysis of the initial value problem concerns uniform energy estimates on the linearized Kakinuma system. To this end, we write the linearized system under the form (1.12)

A 1 (∂ t + u • ∇) U + A mod 0 U = Ḟ ,
where U := ( ζ, φ1 , φ2 ) T is the deviation from the reference state U := (ζ, ϕ 1 , ϕ 2 ) T , u is a suitable velocity which is a convex combination of u 1 and u 2 whose weights depend on ρ ℓ , H ℓ as well as α ℓ (ℓ = 1, 2) the positive constants mentioned previously, Ḟ represent lower order terms and A 1 := A 1 (U ) is a skew-symmetric matrix and A mod 0 := A mod 0 (U ) is a linear operator symmetric in L 2 . The energy function associated to (1.12) is given by (A mod 0 U , U ) L 2 , and we prove that

(A mod 0 U , U ) L 2 ≃ E ( U ) := ∥ ζ∥ 2 L 2 + ℓ=1,2 ρ ℓ (∥∇ φℓ ∥ 2 L 2 + ∥ φ′ ℓ ∥ 2 L 2 )
under the non-cavitation assumption (1.7) and the stability condition (1.5). Because the structure of (1.12) is not standard, the control of the energy function is obtained by testing (1.12) with the time derivatives, ∂ t U . This, together with suitable product and commutator estimates in Sobolev spaces, provides the a priori control of the energy function for solutions to the Kakinuma model and their derivatives, and we show that this control is uniform in the shallow water regime after suitable rescaling. Since the construction and uniqueness of a solution was obtained in the companion paper [START_REF] Duchêne | A mathematical analysis of the Kakinuma model for interfacial gravity waves. Part I: Structures and well-posedness[END_REF], the uniform estimates provide the proof of the long-time well-posedness of the initial value problem for the Kakinuma model result stated in Theorem 1.1. Furthermore, using the aforementioned consistency result, we prove that the difference between solutions to the full interfacial gravity waves system and corresponding solutions to the Kakinuma model satisfy an identity analogous to (1.12), and hence infer a control of the energy function of the difference and its derivatives, which yields the full justification of the Kakinuma model stated in Theorem 1.4.

Outline The contents of this paper are as follows. In Section 2 we first recall the basic equations governing the interfacial gravity waves and write down the [P, Q] = P Q-QP denotes the commutator and [P ; u, v] = P (uv)-(P u)v -u(P v) denotes the symmetric commutator. For a matrix A we denote by A T the transpose of A. O denotes a zero matrix. For a vector ϕ = (ϕ 0 , ϕ 1 , . . . , ϕ N ) T we denote the last N components by ϕ ′ = (ϕ 1 , . . . , ϕ N ) T . f ≲ g means that there exists a non-essential positive constant C such that f ≤ Cg holds. f ≃ g means that f ≲ g and g ≲ f hold. and the bottom at time t by Ω 1 (t), Ω 2 (t), Γ(t), Σ 1 , and Σ 2 , respectively. The velocity potentials Φ 1 (x, z, t) and Φ 2 (x, z, t) in the upper and lower layers, respectively, satisfy Laplace's equations

∆Φ 1 + ∂ 2 z Φ 1 = 0 in Ω 1 (t), (2.1) ∆Φ 2 + ∂ 2 z Φ 2 = 0 in Ω 2 (t), (2.2) where ∆ = ∂ 2 1 + • • • + ∂ 2 n
is the Laplacian with respect to the horizontal space variables x = (x 1 , . . . , x n ). Bernoulli's laws of each layers have the form

ρ 1 ∂ t Φ 1 + 1 2 (|∇Φ 1 | 2 + (∂ z Φ 1 ) 2 ) + gz + P 1 = 0 in Ω 1 (t), (2.3) ρ 2 ∂ t Φ 2 + 1 2 (|∇Φ 2 | 2 + (∂ z Φ 2 ) 2 ) + gz + P 2 = 0 in Ω 2 (t), (2.4)
where ∇ = (∂ 1 , . . . , ∂ n ), the positive constant g is the acceleration due to gravity, and P 1 (x, z, t) and P 2 (x, z, t) are pressures in the upper and lower layers, respectively. The dynamical boundary condition on the interface is given by (2.5)

P 1 = P 2 on Γ(t).
The kinematic boundary conditions on the interface, the rigid-lid, and the bottom are given by

∂ t ζ + ∇Φ 1 • ∇ζ -∂ z Φ 1 = 0 on Γ(t), (2.6) ∂ t ζ + ∇Φ 2 • ∇ζ -∂ z Φ 2 = 0 on Γ(t), (2.7) ∂ z Φ 1 = 0 on Σ 1 , (2.8) ∇Φ 2 • ∇b -∂ z Φ 2 = 0 on Σ 2 . (2.9)
These are the basic equations for interfacial gravity waves. It follows from Bernoulli's laws (2.3)-(2.4) and the dynamical boundary condition (2.5) that

ρ 1 ∂ t Φ 1 + 1 2 (|∇Φ 1 | 2 + (∂ z Φ 1 ) 2 ) (2.10) -ρ 2 ∂ t Φ 2 + 1 2 (|∇Φ 2 | 2 + (∂ z Φ 2 ) 2 ) = (ρ 2 -ρ 1 )gζ on Γ(t).
We will always assume the stable stratification condition (ρ 2 -ρ 1 )g > 0. As in the case of surface water waves, the basic equations have a variational structure and the corresponding Luke's Lagrangian is given, up to terms which do not contribute to the variation of the Lagrangian, by the vertical integral of the pressure in the water regions. After using Bernoulli's laws (2.3)-(2.4) we can find the Lagrangian density

L (Φ 1 , Φ 2 , ζ) = -ρ 1 h 1 ζ ∂ t Φ 1 + 1 2 (|∇Φ 1 | 2 + (∂ z Φ 1 ) 2 ) dz (2.11) -ρ 2 ζ -h 2 +b ∂ t Φ 2 + 1 2 (|∇Φ 2 | 2 + (∂ z Φ 2 ) 2 ) dz - 1 2 (ρ 2 -ρ 1 )gζ 2 .
In fact, one checks readily that (2.1)-(2.2) and (2.6)-(2.10) are Euler-Lagrange equations associated with the action function

J (Φ 1 , Φ 2 , ζ) := t 1 t 0 R n L (Φ 1 , Φ 2 , ζ) dx dt.
We proceed to the Kakinuma model. Let N and N * be nonnegative integers. In view of the analysis for the Isobe-Kakinuma model for surface water waves, we approximate the velocity potentials Φ 1 and Φ 2 in the Lagrangian by (2.12)

             Φ app 1 (x, z, t) = N i=0 (-z + h 1 ) 2i ϕ 1,i (x, t), Φ app 2 (x, z, t) = N * i=0 (z + h 2 -b(x)) p i ϕ 2,i (x, t),
where p 0 , p 1 , . . . , p N * are nonnegative integers satisfying 0 = p 0 < p 1 < • • • < p N * . Plugging (2.12) into the Lagrangian density (2.11), we obtain an approximate Lagrangian density

L app (ϕ 1 , ϕ 2 , ζ) := L (Φ app 1 , Φ app 2 , ζ),
where ϕ 1 := (ϕ 1,0 , ϕ 1,1 , . . . , ϕ 1,N ) T and ϕ 2 := (ϕ 2,0 , ϕ 2,1 , . . . , ϕ 2,N * ) T . The corresponding Euler-Lagrange equation is the Kakinuma model, which has the form (2.13)

                                                                   H 2i 1 ∂ t ζ - N j=0 ∇• 1 2(i + j) + 1 H 2(i+j)+1 1 ∇ϕ 1,j - 4ij 2(i + j) -1 H 2(i+j)-1 1 ϕ 1,j = 0 for i = 0, 1, . . . , N, H p i 2 ∂ t ζ + N * j=0 ∇ • 1 p i + p j + 1 H p i +p j +1 2 ∇ϕ 2,j - p j p i + p j H p i +p j 2 ϕ 2,j ∇b + p i p i + p j H p i +p j 2 ∇b • ∇ϕ 2,j - p i p j p i + p j -1 H p i +p j -1 2 (1 + |∇b| 2 )ϕ 2,j = 0 for i = 0, 1, . . . , N * , ρ 1 N j=0 H 2j 1 ∂ t ϕ 1,j + gζ + 1 2 N j=0 H 2j 1 ∇ϕ 1,j 2 + N j=0 2jH 2j-1 1 ϕ 1,j 2 -ρ 2 N * j=0 H p j 2 ∂ t ϕ 2,j + gζ + 1 2 N * j=0 (H p j 2 ∇ϕ 2,j -p j H p j -1 2 ϕ 2,j ∇b) 2 + N * j=0 p j H p j -1 2 ϕ 2,j 2 = 0,
where H 1 and H 2 are depths of the upper and the lower layers, that is,

H 1 (t, x) := h 1 -ζ(x, t), H 2 (x, t) := h 2 + ζ(x, t) -b(x).
In (2.13), we used the notational convention 0/0 = 0. More precisely, this convention was used so as to dictate p 0 /(p 0 + p 0 ) = 0 and p 0 p 1 /(p 0 + p 1 -1) = 0 in the case p 1 = 1. We recall also that p 0 = 0 is always assumed.

The dimensionless equations

In order to rigorously validate the Kakinuma model (2.13) as a higher order shallow water approximation of the full model for interfacial gravity waves (2.1)-(2.9), we first introduce nondimensional parameters and then non-dimensionalize the equations, through a convenient rescaling of variables. Let λ be a typical horizontal wavelength. Following D. Lannes [START_REF] Lannes | A stability criterion for two-fluid interfaces and applications[END_REF], we introduce a nondimensional parameter δ by

δ := h λ with h := h 1 h 2 ρ 1 h 2 + ρ 2 h 1 ,
where ρ 1 and ρ 2 are relative densities. We also need to use relative depths h 1 and h 2 of the layers. These nondimensional parameters are defined by

ρ ℓ := ρ ℓ ρ 1 + ρ 2 , h ℓ := h ℓ h (ℓ = 1, 2),
which satisfy the relations (2.14)

ρ 1 + ρ 2 = 1, ρ 1 h 1 + ρ 2 h 2 = 1. Note also that min{h 1 , h 2 } ≤ h ≤ max{h 1 , h 2 }. It follows from the second relation in (2.14) that (2.15) 1 < min h 1 ρ 1 , h 2 ρ 2 ≤ 2.
Here, we note that the standard shallowness parameters δ 1 := h 1 λ and δ 2 := h 2 λ relative to the upper and the lower layers, respectively, are related to the above parameters by δ ℓ = h ℓ δ for ℓ = 1, 2. In many results of this paper, we restrict our consideration to the parameter regime

(2.16) h -1 1 + h -1 2 ≲ 1. To understand this restriction, it is convenient to use nondimensional parameters γ := ρ 1 ρ 2 and θ := h 1 h 2 .
In terms of these parameters, h -1 ℓ (ℓ = 1, 2) can be represented as

h -1 1 = γ + 1 γ + θ , h -1 2 = γ -1 + 1 γ -1 + θ -1 .
Therefore, the only cases that (2.16) excludes are the case γ, θ ≪ 1 and the case γ, θ ≫ 1. Since we shall also assume the stable stratification condition (ρ 2 -ρ 1 )g > 0, we can describe the two regimes considered in this paper as

(i) γ ≃ 1, i.e., ρ 1 ≃ ρ 2 , (ii) γ ≪ 1 and θ ≳ 1, i.e., ρ 1 ≪ ρ 2 and h 2 ≲ h 1 .
Introducing c SW := (ρ 2 -ρ 1 )gh the speed of infinitely long and small interfacial gravity waves, we rescale the independent and the dependent variables by

x = λ x, z = hz, t = λ c SW t, ζ = h ζ, b = h b, Φ ℓ = λc SW Φℓ (ℓ = 1, 2).
Plugging these into the full model (2.1)-(2.2) and (2.6)-(2.10) and dropping the tilde sign in the notation we obtain

                                 ∆Φ 1 + δ -2 ∂ 2 z Φ 1 = 0 in Ω 1 (t), ∆Φ 2 + δ -2 ∂ 2 z Φ 2 = 0 in Ω 2 (t), ∂ t ζ + ∇Φ 1 • ∇ζ -δ -2 ∂ z Φ 1 = 0 on Γ(t), ∂ t ζ + ∇Φ 2 • ∇ζ -δ -2 ∂ z Φ 2 = 0 on Γ(t), ∂ z Φ 1 = 0 on Σ 1 , ∇Φ 2 • ∇b -δ -2 ∂ z Φ 2 = 0 on Σ 2 , ρ 1 ∂ t Φ 1 + 1 2 |∇Φ 1 | 2 + 1 2 δ -2 (∂ z Φ 1 ) 2 -ρ 2 ∂ t Φ 2 + 1 2 |∇Φ 2 | 2 + 1 2 δ -2 (∂ z Φ 2 ) 2 -ζ = 0 on Γ(t),
where in this scaling the upper layer Ω 1 (t), the lower layer Ω 2 (t), the interface Γ(t), the rigid-lid Σ 1 , and the bottom Σ 2 are written as

Ω 1 (t) = {X = (x, z) ∈ R n+1 ; ζ(x, t) < z < h 1 }, Ω 2 (t) = {X = (x, z) ∈ R n+1 ; -h 2 + b(x) < z < ζ(x, t)}, Γ(t) = {X = (x, z) ∈ R n+1 ; z = ζ(x, t)}, Σ 1 = {X = (x, z) ∈ R n+1 ; z = h 1 }, Σ 2 = {X = (x, z) ∈ R n+1 ; z = -h 2 + b(x)}. Denoting ϕ ℓ (x, t) := Φ ℓ (x, ζ(x, t), t) (ℓ = 1, 2)
and using the chain rule, the above system can be written in a more compact and closed form as (2.17)

                 ∂ t ζ + Λ 1 (ζ, δ, h 1 )ϕ 1 = 0, ∂ t ζ -Λ 2 (ζ, b, δ, h 2 )ϕ 2 = 0, ρ 1 ∂ t ϕ 1 + 1 2 |∇ϕ 1 | 2 - 1 2 δ 2 (Λ 1 (ζ, δ, h 1 )ϕ 1 -∇ζ • ∇ϕ 1 ) 2 1 + δ 2 |∇ζ| 2 -ρ 2 ∂ t ϕ 2 + 1 2 |∇ϕ 2 | 2 - 1 2 δ 2 (Λ 2 (ζ, b, δ, h 2 )ϕ 2 + ∇ζ • ∇ϕ 2 ) 2 1 + δ 2 |∇ζ| 2 -ζ = 0,
where Λ 1 (ζ, δ, h 1 ) and Λ 2 (ζ, b, δ, h 2 ) are the Dirichlet-to-Neumann maps for Laplace's equations. More precisely, these are defined by

Λ 1 (ζ, δ, h 1 )ϕ 1 := -δ -2 ∂ z Φ 1 + ∇Φ 1 • ∇ζ z=ζ(x,t) , Λ 2 (ζ, b, δ, h 2 )ϕ 2 := δ -2 ∂ z Φ 2 -∇Φ 2 • ∇ζ z=ζ(x,t) ,
where Φ 1 and Φ 2 are unique solutions to the boundary value problems

     ∆Φ 1 + δ -2 ∂ 2 z Φ 1 = 0 in Ω 1 (t), Φ 1 = ϕ 1 on Γ(t), ∂ z Φ 1 = 0 on Σ 1 , and
     ∆Φ 2 + δ -2 ∂ 2 z Φ 2 = 0 in Ω 2 (t), Φ 2 = ϕ 2 on Γ(t), ∇Φ 2 • ∇b -δ -2 ∂ z Φ 2 = 0 on Σ 2 .
As for the Kakinuma model, we introduce additionally the rescaled variables

ϕ 1,i := λc SW h 2i 1 φ1,i , ϕ 2,i := λc SW h p i 2 φ2,i ,
where we recall that p 0 , p 1 , . . . , p N * are nonnegative integers satisfying 0 = p 0 < p 1 < • • • < p N * appearing in the approximation (2.12). Plugging these and the previous scaling into the Kakinuma model (2.13) and dropping the tilde sign in the notation we obtain the Kakinuma model in the nondimensional form, which is written as

(2.18)                                                                H 2i 1 ∂ t ζ -h 1 N j=0 ∇• 1 2(i + j) + 1 H 2(i+j)+1 1 ∇ϕ 1,j - 4ij 2(i + j) -1 H 2(i+j)-1 1 (h 1 δ) -2 ϕ 1,j = 0 for i = 0, 1, . . . , N, H p i 2 ∂ t ζ + h 2 N * j=0 ∇ • 1 p i + p j + 1 H p i +p j +1 2 ∇ϕ 2,j - p j p i + p j H p i +p j 2 ϕ 2,j h -1 2 ∇b + p i p i + p j H p i +p j 2 h -1 2 ∇b • ∇ϕ 2,j - p i p j p i + p j -1 H p i +p j -1 2
((h 2 δ) -2 + h -2 2 |∇b| 2 )ϕ 2,j = 0 for i = 0, 1, . . . , N * , ρ 1 N j=0 H 2j 1 ∂ t ϕ 1,j + 1 2 N j=0 H 2j 1 ∇ϕ 1,j 2 + (h 1 δ) -2 N j=0 2jH 2j-1 1 ϕ 1,j 2 -ρ 2 N * j=0 H 2j 2 ∂ t ϕ 2,j + 1 2 N * j=0 (H p j 2 ∇ϕ 2,j -p j H p j -1 2 ϕ 2,j h -1 2 ∇b) 2 + (h 2 δ) -2 N * j=0 p j H p j -1 2 ϕ 2,j 2 -ζ = 0,
where we used the notational convention 0/0 = 0, and

(2.19) H 1 (x, t) := 1 -h -1 1 ζ(x, t), H 2 (x, t) := 1 + h -1 2 ζ(x, t) -h -1 2 b(x).
We impose the initial conditions to the Kakinuma model of the form

(2.20) (ζ, ϕ 1 , ϕ 2 ) = (ζ (0) , ϕ 1(0) , ϕ 2(0) ) at t = 0.

Hamiltonian structures

T. B. Benjamin and T. J. Bridges [START_REF] Benjamin | Reappraisal of the Kelvin-Helmholtz problem. Part 1. Hamiltonian structure[END_REF] found that the full model for interfacial gravity waves can be written in Hamilton's canonical form

∂ t ζ = δH IW δϕ , ∂ t ϕ = - δH IW δζ ,
where the canonical variable ϕ is defined by

(2.21) ϕ = ρ 2 ϕ 2 -ρ 1 ϕ 1
and the Hamiltonian H IW is the total energy E written in terms of the canonical variables (ζ, ϕ). Specifically, E is the sum of the kinetic energies of the fluids in the upper and the lower layers and the potential energy due to the gravity defined as

E := ℓ=1,2 Ω ℓ (t) 1 2 ρ ℓ |∇Φ ℓ (x, z, t)| 2 + δ -2 (∂ z Φ ℓ (x, z, t)) 2 dxdz + R n 1 2 ζ(x, t) 2 dx = ℓ=1,2 1 2 ρ ℓ (Λ ℓ (ζ)ϕ ℓ (t), ϕ ℓ (t)) L 2 + 1 2 ∥ζ(t)∥ 2 L 2 .
Here and in what follows, we denote simply Λ

1 (ζ) = Λ 1 (ζ, δ, h 1 ) and Λ 2 (ζ) = Λ 2 (ζ, b, δ, h 2 )
. It follows from the kinematic boundary conditions on the interface that Λ 1 (ζ)ϕ 1 + Λ 2 (ζ)ϕ 2 = 0, so that ϕ 1 and ϕ 2 can be written in terms of the canonical variables (ζ, ϕ) as

ϕ 1 = -(ρ 1 Λ 2 (ζ) + ρ 2 Λ 1 (ζ)) -1 Λ 2 (ζ)ϕ, ϕ 2 = (ρ 1 Λ 2 (ζ) + ρ 2 Λ 1 (ζ)) -1 Λ 1 (ζ)ϕ.
Therefore, the Hamiltonian H IW (ζ, ϕ) of the full model for interfacial gravity waves is given explicitly by

(2.22) H IW (ζ, ϕ) = 1 2 ((ρ 1 Λ 2 (ζ) + ρ 2 Λ 1 (ζ)) -1 Λ 1 (ζ)ϕ, Λ 2 (ζ)ϕ) L 2 + 1 2 ∥ζ∥ 2 L 2 .
As was shown in the companion paper [START_REF] Duchêne | A mathematical analysis of the Kakinuma model for interfacial gravity waves. Part I: Structures and well-posedness[END_REF], the Kakinuma model (2.18) also enjoys a Hamiltonian structure analogous to that of the full model for interfacial gravity waves. The canonical variables are the elevation of the interface ζ and ϕ defined by

ϕ(x, t) := ρ 2 Φ app 2 (x, ζ(x, t), t) -ρ 1 Φ app 1 (x, ζ(x, t), t) (2.23) = ρ 2 N * i=0 H 2 (x, t) p i ϕ 2,i (x, t) -ρ 1 N i=0 H 1 (x, t) 2i ϕ 1,i (x, t),
where Φ app ℓ (ℓ = 1, 2) are nondimensional versions of the approximate velocity potentials, which are defined by

(2.24)              Φ app 1 (x, z, t) := N i=0 (1 -h -1 1 z) 2i ϕ 1,i (x, t), Φ app 2 (x, z, t) := N * i=0 (1 + h -1 2 (z -b(x))) p i ϕ 2,i (x,
H K (ζ, ϕ) := ℓ=1,2 Ω ℓ 1 2 ρ ℓ |∇Φ app ℓ (x, z, t)| 2 + δ -2 (∂ z Φ app ℓ (x, z, t)) 2 dxdz + R n 1 2 ζ(x, t) 2 dx.

Statements of the main results

Before stating the main results in this paper, let us introduce some notations which allow in particular to rewrite (2.18) in a compact form. We introduce second order differential operators

L 1,ij = L 1,ij (H 1 , δ, h 1 ) (i, j = 0, 1, . . . , N ) and L 2,ij = L 2,ij (H 2 , b, δ, h 2 ) (i, j = 0, 1, . . . , N * ) by L 1,ij φ 1,j := -∇ • 1 2(i + j) + 1 H 2(i+j)+1 1 ∇φ 1,j + 4ij 2(i + j) -1 H 2(i+j)-1 1 (h 1 δ) -2 φ 1,j , (3.1) L 2,ij φ 2,j := -∇ • 1 p i + p j + 1 H p i +p j +1 2 ∇φ 2,j - p j p i + p j H p i +p j 2 φ 2,j h -1 2 ∇b (3.2) - p i p i + p j H p i +p j 2 h -1 2 ∇b • ∇φ 2,j + p i p j p i + p j -1 H p i +p j -1 2 ((h 2 δ) -2 + h -2 2 |∇b| 2 )φ 2,j ,
where we use the notational convention 0/0 = 0. Notice that we have (L ℓ,ij

) * = L ℓ,ji for ℓ = 1, 2, where (L ℓ,ij ) * is the adjoint operator of L ℓ,ij in L 2 (R n ). We put ϕ 1 := (ϕ 1,0 , ϕ 1,1 , . . . , ϕ 1,N ) T , ϕ 2 := (ϕ 2,0 , ϕ 2,1 , . . . , ϕ 2,N * ) T , and (3.3)
                     l 1 (H 1 ) := (1, H 2 1 , H 4 1 , . . . , H 2N 1 ) T , l ′ 1 (H 1 ) := (0, 2H 1 , . . . , 2N H 2N -1 1 ) T , l ′′ 1 (H 1 ) := (0, 2, . . . , 2N (2N -1)H 2N -2 1 ) T , l 2 (H 2 ) := (1, H p 1 2 , H p 2 2 , . . . , H p N * 2 ) T , l ′ 2 (H 2 ) := (0, p 1 H p 1 -1 2 , . . . , p N * H p N * 2 ) T , l ′′ 2 (H 2 ) := (0, p 1 (p 1 -1)H p 1 -2 2 , . . . , p N * (p N * -1)H p N * 2
) T , and define u ℓ and w ℓ for ℓ = 1, 2, which represent approximately the horizontal and the vertical components of the velocity field on the interface from the water region Ω ℓ (t), by (3.4)

u 1 := (l 1 (H 1 ) ⊗ ∇) T ϕ 1 , w 1 := -l ′ 1 (H 1 ) • ϕ 1 , u 2 := (l 2 (H 2 ) ⊗ ∇) T ϕ 2 -(l ′ 2 (H 2 ) • ϕ 2 )h -1 2 ∇b, w 2 := l ′ 2 (H 2 ) • ϕ 2 .
Then, denoting L 1 := (L 1,ij ) 0≤i,j≤N and L 2 := (L 2,ij ) 0≤i,j≤N * we can write the Kakinuma model (2.18) more compactly as

(3.5)            l 1 (H 1 )∂ t ζ + h 1 L 1 (H 1 , δ, h 1 )ϕ 1 = 0, l 2 (H 2 )∂ t ζ -h 2 L 2 (H 2 , b, δ, h 2 )ϕ 2 = 0, ρ 1 l 1 (H 1 ) • ∂ t ϕ 1 + 1 2 |u 1 | 2 + (h 1 δ) -2 w 2 1 -ρ 2 l 2 (H 2 ) • ∂ t ϕ 2 + 1 2 |u 2 | 2 + (h 2 δ) -2 w 2 2 -ζ = 0.
By eliminating ∂ t ζ from the first two vectorial identities in (3.5), we obtain N + N * + 1 scalar relations which are necessary conditions for the existence of solutions to the Kakinuma model, as stated below. Introducing linear operators

L 1,i := L 1,i (H 1 , δ, h 1 ) (i = 0, . . . , N ) acting on φ 1 = (φ 1,0 , . . . , φ 1,N ) T and L 2,i := L 2,i (H 2 , b, δ, h 2 ) (i = 0, . . . , N * ) acting on φ 2 = (φ 2,0 , . . . , φ 2,N * ) T by (3.6)                                    L 1,0 φ 1 := N j=0 L 1,0j φ 1,j , L 1,i φ 1 := N j=0 (L 1,ij φ 1,j -H 2i 1 L 1,0j φ 1,j ) for i = 1, 2, . . . , N, L 2,0 φ 2 := N * j=0 L 2,0j φ 2,j , L 2,i φ 2 := N * j=0 (L 2,ij φ 2,j -H p i 2 L 2,0j φ 2,j ) for i = 1, 2, . . . , N * ,
the necessary conditions can be written simply as

(3.7)      L 1,i (H 1 , δ, h 1 )ϕ 1 = 0 for i = 1, 2, . . . , N, L 2,i (H 2 , b, δ, h 2 )ϕ 2 = 0 for i = 1, 2, . . . , N * , h 1 L 1,0 (H 1 , δ, h 1 )ϕ 1 + h 2 L 2,0 (H 2 , b, δ, h 2 )ϕ 2 = 0.
Hereafter, these necessary conditions will be referred to as the compatibility conditions. Notice that under these compatibility conditions we have for ℓ = 1, 2

(3.8) L ℓ ϕ ℓ = l ℓ L ℓ,0 ϕ ℓ ,
where l ℓ = l ℓ (H ℓ ) and similar simplifications of notations will be used in the following without any comments. In connection with the stability condition (1.5), we introduce a function

a := 1 + ρ 1 h -1 1 {l ′ 1 (H 1 ) • (∂ t + u 1 • ∇)ϕ 1 -(h 1 δ) -2 w 1 l ′′ 1 (H 1 ) • ϕ 1 } (3.9) + ρ 2 h -1 2 {l ′ 2 (H 2 ) • (∂ t + u 2 • ∇)ϕ 2 + (h 2 δ) -2 w 2 -h -1 2 ∇b • u 2 l ′′ 2 (H 2 ) • ϕ 2 }, which corresponds to -(∂ z (P app 2 -P app 1 ))| Γ(t)
in the stability condition. Our first main result in this paper is the existence of the solution to the initial value problem (2.18)-(2.20) for the Kakinuma model on a time interval independent of parameters, especially, the shallowness parameters δ 1 = h 1 δ and δ 2 = h 2 δ together with a uniform bound of the solution. For simplicity, we denote H ℓ(0) := H ℓ | t=0 , u ℓ(0) := u ℓ | t=0 for ℓ = 1, 2, and a (0) := a| t=0 , which can be written in terms of the initial data according to the initial condition (2.20). Although the function a includes the terms (∂ t ϕ ′ ℓ )| t=0 for ℓ = 1, 2, where ϕ ′ 1 = (ϕ 1,1 , . . . , ϕ 1,N ) T and ϕ ′ 2 = (ϕ 2,1 , . . . , ϕ 2,N * ) T , and the hypersurface t = 0 is characteristic for the Kakinuma model, we can uniquely determine them in terms of the initial data. For details, we refer to Remark 5.3. Theorem 3.1. Let c 0 , M 0 , h min be positive constants and m an integer such that m > n 2 + 1. There exist a time T > 0 and a constant M > 0 such that for any positive parameters ρ 1 , ρ 2 , h 1 , h 2 , δ satisfying the natural restrictions (2.14), h 1 δ, h 2 δ ≤ 1, as well as the condition h min ≤ h 1 , h 2 , if the initial data (ζ (0) , ϕ 1(0) , ϕ 2(0) ) and the bottom topography b satisfy

(3.10)      ∥ζ (0) ∥ 2 H m + ℓ=1,2 ρ ℓ h ℓ ∥∇ϕ ℓ(0) ∥ 2 H m + (h ℓ δ) -2 ∥ϕ ′ ℓ(0) ∥ 2 H m ≤ M 0 , h -1 2 ∥b∥ W m+1,∞ + (h 2 δ)∥b∥ W m+2,∞ ≤ M 0 ,
the non-cavitation assumption

(3.11) H 1(0) (x) ≥ c 0 , H 2(0) (x) ≥ c 0 for x ∈ R n ,
the stability condition

(3.12) a (0) (x) - ρ 1 ρ 2 ρ 1 h 2 H 2(0) (x)α 2 + ρ 2 h 1 H 1(0) (x)α 1 |u 1(0) (x) -u 2(0) (x)| 2 ≥ c 0 for x ∈ R n ,
with positive constants α 1 and α 2 defined by (3.16), and the compatibility conditions

(3.13)      L 1,i (H 1(0) , δ, h 1 )ϕ 1(0) = 0 for i = 1, 2, . . . , N, L 2,i (H 2(0) , b, δ, h 2 )ϕ 2(0) = 0 for i = 1, 2, . . . , N * , h 1 L 1,0 (H 1(0) , δ, h 1 )ϕ 1(0) + h 2 L 2,0 (H 2(0) , b, δ, h 2 )ϕ 2(0) = 0, then the initial value problem (2.18)-(2.20) has a unique solution (ζ, ϕ 1 , ϕ 2 ) on the time interval [0, T ] satisfying ζ, ∇ϕ 1,0 , ∇ϕ 2,0 ∈ C([0, T ]; H m ) ∩ C 1 ([0, T ]; H m-1 ), ϕ ′ 1 , ϕ ′ 2 ∈ C([0, T ]; H m+1 ) ∩ C 1 ([0, T ]; H m
), where we recall the notation ϕ ′ 1 = (ϕ 1,1 , ϕ 1,2 , . . . , ϕ 1,N ) T and ϕ ′ 2 = (ϕ 2,1 , ϕ 2,2 , . . . , ϕ 2,N * ) T . Moreover, the solution satisfies the uniform bound

(3.14) ∥ζ(t)∥ 2 H m + ℓ=1,2 ρ ℓ h ℓ ∥∇ϕ ℓ (t)∥ 2 H m + (h ℓ δ) -2 ∥ϕ ′ ℓ (t)∥ 2 H m ≤ M for t ∈ [0, T ] together with (3.15)    a(x, t) - ρ 1 ρ 2 ρ 1 h 2 H 2 (x, t)α 2 + ρ 2 h 1 H 1 (x, t)α 1 |u 1 (x, t) -u 2 (x, t)| 2 ≥ c 0 /2, H 1 (x, t) ≥ c 0 /2, H 2 (x, t) ≥ c 0 /2 for x ∈ R n , t ∈ [0, T ].
Remark 3.2. The constants α 1 and α 2 are defined by

(3.16) α ℓ := det A ℓ,0 det Ãℓ,0 , Ãℓ,0 := 0 1 T -1 A ℓ,0 , for ℓ = 1, 2
, where 1 := (1, . . . , 1) T and the matrices A 1,0 and A 2,0 are defined by

       A 1,0 := 1 2(i + j) + 1 0≤i,j≤N , A 2,0 := 1 p i + p j + 1 0≤i,j≤N * .
Hence, α 1 and α 2 are positive constants depending only on N and the nonnegative integers 0 = p 0 < p 1 < . . . < p N * , respectively, and go to 0 as N, N * → ∞.

Remark 3.3. It is easy to check that the non-cavitation assumption (3.11) and the stability condition (3.12) are automatically satisfied for small initial data (ζ (0) , ϕ 1(0) , ϕ 2(0) ) and small bottom topography b, whereas an arrangement of nontrivial initial data satisfying the compatibility conditions (3.13) together with the uniform bound (3.10) is a non-trivial issue. To this end, we use the canonical variable ϕ defined by (2.23), which can be written as

(3.17) ϕ = ρ 2 l 2 (H 2 ) • ϕ 2 -ρ 1 l 1 (H 1 ) • ϕ 1 .
Given the initial data (ζ (0) , ϕ (0) ) for the canonical variables (ζ, ϕ), and the bottom topography b, the necessary conditions (3.7) and the above relation (3.17) determine the initial data (ϕ 1(0) , ϕ 2(0) ) for the Kakinuma model (2.18)-(2.20) satisfying the compatibility conditions (3.13) and the uniform bound (3.10), which is unique up to an additive constant of the form (Cρ 2 , Cρ 1 ) to (ϕ 1,0(0) , ϕ 2,0(0) ). In fact, we have the following proposition, which is a simple corollary of Lemma 5.1 given in Section 5.

Proposition 3.4. Let c 0 , M 0 be positive constants and m an integer such that m > n 2 +1. There exists a positive constant C such that for any positive parameters ρ 1 , ρ 2 , h 1 , h 2 , δ satisfying the natural restrictions (2.14)

and h 1 δ, h 2 δ ≤ 1, if the initial data (ζ (0) , ϕ (0) ) ∈ H m × Hm of the canonical variables, the bottom topography b ∈ W m,∞ , and initial depths H 1(0) := 1 -h -1 1 ζ (0) and H 2(0) := 1 + h -1 2 ζ (0) -h -1 2 b satisfy h -1 1 ∥ζ (0) ∥ H m + h -1 2 ∥ζ (0) ∥ H m + h -1 2 ∥b∥ W m,∞ ≤ M 0 , H 1(0) (x) ≥ c 0 , H 2(0) (x) ≥ c 0 for x ∈ R n ,
then there exist initial data (ϕ 1(0) , ϕ 2(0) ) satisfying the compatibility conditions (3.13) as well as

ϕ (0) = ρ 2 l 2 (H 2(0) ) • ϕ 2(0) -ρ 1 l 1 (H 1(0) ) • ϕ 1(0) . Moreover, we have ℓ=1,2 ρ ℓ h ℓ ∥∇ϕ ℓ(0) ∥ 2 H m-1 + (h ℓ δ) -2 ∥ϕ ′ ℓ(0) ∥ 2 H m-1 ≤ C∥∇ϕ (0) ∥ 2 H m-1 .
The next theorem shows that the Kakinuma model (2.18) is consistent with the full model for interfacial gravity waves (2.17) at order O((h 1 δ) 4N +2 + (h 2 δ) 4N +2 ) under the special choice of the indices p 0 , p 1 , . . . , p N * as (H1) N * = N and p i = 2i (i = 0, 1, . . . , N ) in the case of the flat bottom b(x) ≡ 0, (H2) N * = 2N and p i = i (i = 0, 1, . . . , 2N ) in the case with general bottom topographies. Theorem 3.5. Let c, M be positive constants and m an integer such that m ≥ 4(N + 1) and m > n 2 + 1. We assume (H1) or (H2). There exists a positive constant C such that for any positive parameters ρ 1 , ρ 2 , h 1 , h 2 , δ satisfying h 1 δ, h 2 δ ≤ 1 and for any solution (ζ, ϕ 1 , ϕ 2 ) to the Kakinuma model (2.18) on a time interval [0, T ] with a bottom topography b ∈ W m+1,∞ satisfying

(3.18) h -1 1 ∥ζ(t)∥ H m + h -1 2 ∥ζ(t)∥ H m + h -1 2 ∥b∥ W m+1,∞ ≤ M, H 1 (x, t) ≥ c, H 2 (x, t) ≥ c for x ∈ R n , t ∈ [0, T ], if we define ϕ ℓ := l ℓ (H ℓ ) • ϕ ℓ for ℓ = 1, 2, then (ζ, ϕ 1 , ϕ 2 )
satisfy approximately the full model for interfacial gravity waves as

                 ∂ t ζ + Λ 1 (ζ, δ, h 1 )ϕ 1 = r 1 , ∂ t ζ -Λ 2 (ζ, b, δ, h 2 )ϕ 2 = r 2 , ρ 1 ∂ t ϕ 1 + 1 2 |∇ϕ 1 | 2 - 1 2 δ 2 (Λ 1 (ζ, δ, h 1 )ϕ 1 -∇ζ • ∇ϕ 1 ) 2 1 + δ 2 |∇ζ| 2 -ρ 2 ∂ t ϕ 2 + 1 2 |∇ϕ 2 | 2 - 1 2 δ 2 (Λ 2 (ζ, b, δ, h 2 )ϕ 2 + ∇ζ • ∇ϕ 2 ) 2 1 + δ 2 |∇ζ| 2 -ζ = r 0 .
Here, the errors (r 1 , r 2 , r 0 ) satisfy

     ∥r ℓ (t)∥ H m-4(N +1) ≤ Ch ℓ (h ℓ δ) 4N +2 ∥∇ϕ ℓ (t)∥ H m-1 (ℓ = 1, 2), ∥r 0 (t)∥ H m-4(N +1) ≤ C ℓ=1,2 ρ ℓ (h ℓ δ) 4N +2 ∥∇ϕ ℓ (t)∥ 2 H m-1 for t ∈ [0, T ].
Particularly, we see that under the special choice of indices (H1) or (H2) the solutions to the Kakinuma model (2.18)-(2.20) constructed in Theorem 3.1 satisfy approximately the full model for interfacial gravity waves (2.17) with the choice

ϕ ℓ = l ℓ (H ℓ ) • ϕ ℓ (ℓ = 1, 2) and that the error is of order O((h 1 δ) 4N +2 + (h 2 δ) 4N +2 ).
Conversely, the next theorem shows that the full model for interfacial gravity waves is consistent with the Kakinuma model at order O((h 1 δ) 4N +2 + (h 2 δ) 4N +2 ) under the special choice of indices (H1) or (H2). Theorem 3.6. Let c, M be positive constants and m an integer such that m ≥ 4(N + 1) and m > n 2 + 1. We assume (H1) or (H2). There exists a positive constant C such that for any positive parameters ρ 1 , ρ 2 , h 1 , h 2 , δ satisfying h 1 δ, h 2 δ ≤ 1 and for any solution (ζ, ϕ 1 , ϕ 2 ) to the full model for interfacial gravity waves (2.17) on a time interval [0, T ] with a bottom topography b ∈ W m+1,∞ satisfying (3.18), if we define H 1 and H 2 as in (2. [START_REF] Iguchi | Isobe-Kakinuma model for water waves as a higher order shallow water approximation[END_REF]) and ϕ 1 and ϕ 2 as the unique solutions to the problems

(3.19) l 1 (H 1 ) • ϕ 1 = ϕ 1 , L 1,i (H 1 , δ, h 1 )ϕ 1 = 0 for i = 1, 2, . . . , N, l 2 (H 2 ) • ϕ 2 = ϕ 2 , L 2,i (H 2 , b, δ, h 2 )ϕ 2 = 0 for i = 1, 2, . . . , N * , then (ζ, ϕ 1 , ϕ 2 )
satisfy approximately the Kakinuma model as

           l 1 (H 1 )h -1 1 ∂ t ζ + L 1 (H 1 , δ, h 1 )ϕ 1 = r1 , l 2 (H 2 )h -1 2 ∂ t ζ -L 2 (H 2 , b, δ, h 2 )ϕ 2 = r2 , ρ 1 l 1 (H 1 ) • ∂ t ϕ 1 + 1 2 |u 1 | 2 + (h 1 δ) -2 w 2 1 -ρ 2 l 2 (H 2 ) • ∂ t ϕ 2 + 1 2 |u 2 | 2 + (h 2 δ) -2 w 2 2 -ζ = r0 .
Here, the errors (r 1 , r2 , r0 ) satisfy

(3.20)      ∥r ℓ (t)∥ H m-4(N +1) ≤ C(h ℓ δ) 4N +2 ∥∇ϕ ℓ (t)∥ H m-1 (ℓ = 1, 2), ∥r 0 (t)∥ H m-4(N +1) ≤ C ℓ=1,2 ρ ℓ (h ℓ δ) 4N +2 ∥∇ϕ ℓ (t)∥
           L 1,i (H 1 , δ, h 1 )ϕ 1 = 0 for i = 1, 2, . . . , N, L 2,i (H 1 , b, δ, h 2 )ϕ 2 = 0 for i = 1, 2, . . . , N * , h 1 L 1,0 (H 1 , δ, h 1 )ϕ 1 + h 2 L 2,0 (H 2 , b, δ, h 2 )ϕ 2 = 0, ρ 2 l 2 (H 2 ) • ϕ 2 -ρ 1 l 1 (H 1 ) • ϕ 1 = ϕ,
where ϕ = ρ 2 ϕ 2 -ρ 1 ϕ 1 is the canonical variable for the full model for interfacial gravity waves.

The above system is nothing but the compatibility conditions (3.7) together with the definition (3.17) of the canonical variable for the Kakinuma model. The existence of the approximate solution (ϕ 1 , ϕ 2 ) is guaranteed by Lemma 5.1 given in Section 5. Then, we have similar error estimates to (3.20). For details, we refer to Proposition 7.6.

The above Theorems 3.5 and 3.6 concern essentially the approximation of the equations. To give a rigorous justification of the Kakinuma model (2.18) as a higher order shallow water approximation to the full model for interfacial gravity waves (2.17), one needs to give an error estimate between solutions to the Kakinuma model and that to the full model. However, we cannot expect to construct general solutions to the initial value problem for the full model for interfacial gravity waves because the initial value problem is ill-posed. Nevertheless, if we assume the existence of a solution to the full model with a uniform bound with respect to the shallowness parameters δ 1 = h 1 δ and δ 2 = h 2 δ, then we can give an error estimate with respect to a solution to the Kakinuma model by making use of the well-posedness of the initial value problem for the Kakinuma model as we can see in the following theorem. Theorem 3.9. Let c, M, h min be positive constants and m an integer such that m > n 2 + 4(N + 1). We assume (H1) or (H2). Then, there exist a time T > 0 and a constant C > 0 such that the following holds true. Let ρ 1 , ρ 2 , h 1 , h 2 , δ be positive parameters satisfying the natural restrictions (2.14), h 1 δ, h 2 δ ≤ 1, and the condition h min ≤ h 1 , h 2 , and let b ∈ W m+2,∞ such that h -1 2 ∥b∥ W m+2,∞ ≤ M . Suppose that the full model for interfacial gravity waves (2.17) possesses a solution

(ζ IW , ϕ IW 1 , ϕ IW 2 ) ∈ C([0, T IW ]; H m+1 × Hm+1 × Hm+1 ) satisfying a uniform bound      ∥ζ IW (t)∥ 2 H m+1 + ℓ=1,2 ρ ℓ h ℓ ∥∇ϕ IW ℓ (t)∥ 2 H m ≤ M, H IW 1 (x, t) ≥ c, H IW 2 (x, t) ≥ c for x ∈ R n , t ∈ [0, T IW ],
where we denote

H IW 1 := 1 -h -1 1 ζ IW and H IW 2 := 1 + h -1 2 ζ IW -h -1 2 b. Let ζ (0) := ζ IW | t=0 and ϕ (0) := (ρ 2 ϕ IW 2 -ρ 1 ϕ IW 1 )
| t=0 be the initial data for the canonical variables, and let (ϕ 1(0) , ϕ 2(0) ) be the initial data to the Kakinuma model constructed from (ζ (0) , ϕ (0) ) by Proposition 3.4. Assume moreover that the initial data (ζ (0) , ϕ 1(0) , ϕ 2(0) ) satisfy the stability condition (3.12),

let (ζ K , ϕ K 1 , ϕ K
2 ) be the solution to the initial value problem for the Kakinuma model (2.18)-(2.20) on the time interval [0, T ] whose unique existence is guaranteed by Theorem 3.1, and put

ϕ K ℓ = l ℓ (H ℓ ) • ϕ K ℓ for ℓ = 1, 2.
Then, we have the error bound

∥ζ K (t) -ζ IW (t)∥ H m-4(N +1) + ℓ=1,2 ρ ℓ h ℓ ∥∇ϕ K ℓ (t) -∇ϕ IW ℓ (t)∥ H m-(4N +5) ≤ C((h 1 δ) 4N +2 + (h 2 δ) 4N +2 ) for 0 ≤ t ≤ min{T, T IW }.
The next theorem is the final main result in this paper and states the consistency of the Hamiltonian H K (ζ, ϕ) of the Kakinuma model with respect to the Hamiltonian H IW (ζ, ϕ) of the full model for interfacial gravity waves. We recall that these Hamiltonians are defined in (2.25) and (2.22), respectively. Theorem 3.10. Let c, M, h min be positive constants and m an integer such that m > n 2 + 1 and m ≥ 4(N + 1). We assume (H1) or (H2). There exists a positive constant C such that for any positive parameters ρ 1 , ρ 2 , h 1 , h 2 , δ satisfying the natural restrictions (2.14), h 1 δ, h 2 δ ≤ 1, and the condition h min ≤ h 1 , h 2 , and for any

(ζ, ϕ) ∈ H m × H4(N+1) and b ∈ W m+1,∞ satisfying h -1 1 ∥ζ∥ H m + h -1 2 ∥ζ∥ H m + h -1 2 ∥b∥ W m+1,∞ ≤ M, H 1 (x) ≥ c, H 2 (x) ≥ c for x ∈ R n ,
with H 1 and H 2 defined by (2.19), we have

|H K (ζ, ϕ) -H IW (ζ, ϕ)| ≤ C∥∇ϕ∥ H 4N +3 ∥∇ϕ∥ L 2 ((h 1 δ) 4N +2 + (h 2 δ) 4N +2 ).
4 Consistency of the Kakinuma model; proof of Theorems 3.5 and 3.6

In this section we show that under the special choice of the indices p 0 , p 1 , . . . , p N * as (H1) N * = N and p i = 2i (i = 0, 1, . . . , N ) in the case of the flat bottom b(x) ≡ 0, (H2) N * = 2N and p i = i (i = 0, 1, . . . , 2N ) in the case with general bottom topographies, the Kakinuma model (2.18)is a higher order model to the full model for interfacial gravity waves (2.17) in the limit δ 1 = h 1 δ → 0, δ 2 = h 2 δ → 0, in the sense of consistency. Specifically, we prove Theorems 3.5 and 3.6. Our proof relies essentially on results obtained in the framework of surface waves in [START_REF] Iguchi | A mathematical justification of the Isobe-Kakinuma model for water waves with and without bottom topography[END_REF], which are recalled in Subsection 4.1. The extension to the framework of interfacial waves and the completion of the proof are provided in Subsection 4.2.

Results in the framework of surface waves

In this subsection, we consider the case of surface waves where the water surface and the bottom of the water are represented as z = ζ(x) and z = -1 + b(x), respectively. Here, the time t is fixed arbitrarily, so that we omit the dependence of t in notations. Let H(x) = 1 + ζ(x) -b(x) be the water depth. For a nonnegative integer N , let N * and p 0 , p 1 , . . . , p N * be nonnegative integers satisfying the condition (H1) or (H2). Put

L ij φ j := -∇ • 1 p i + p j + 1 H p i +p j +1 ∇φ j - p j p i + p j H p i +p j φ j ∇b (4.2) - p i p i + p j H p i +p j ∇b • ∇φ j + p i p j p i + p j -1 H p i +p j -1 (δ -2 + |∇b| 2 )φ j ,
where we use the notational convention 0/0 = 0. Introduce linear operators L i = L i (H, b, δ) (i = 0, 1, . . . , N * ) acting on φ = (φ 0 , . . . , φ N * ) T by (4.3)

             L 0 φ := N * j=0 L 0j φ j , L i φ := N * j=0 (L ij φ j -H p i L 0j φ j ) for i = 1,
(4.4) ∥ζ∥ H m + ∥b∥ W m,∞ ≤ M, H(x) ≥ c for x ∈ R n ,
then for any k = ±0, . . . , ±(m -1), any δ ∈ (0, 1], and any ϕ ∈ Hk+1 there exists a unique solution ϕ = (ϕ 0 , ϕ 1 , . . . , ϕ N * ) = (ϕ 0 , ϕ ′ ) ∈ Hk+1 × (H k+1 ) N * to the problem

(4.5) L i (H, b, δ)ϕ = 0 for i = 1, 2, . . . , N * , l(H) • ϕ = ϕ.
Moreover, the solution satisfies ∥∇ϕ∥

H k + δ -1 ∥ϕ ′ ∥ H k ≤ C∥∇ϕ∥ H k .
As a corollary of this lemma, under the assumptions of Lemma 4.1

Λ (N ) (ζ, b, δ) : ϕ → L 0 (H, b, δ)ϕ,
where ϕ is the unique solution to (4.5), is defined as a bounded linear operator from Hk+1 to H k-1 for any k = ±0, . . . , ±(m -1). A key result is that the operator Λ (N ) (ζ, b, δ) provide good approximations in the shallow water regime δ ≪ 1 to the corresponding Dirichlet-to-Neumann map Λ(ζ, b, δ), which is defined by

(4.6) Λ(ζ, b, δ)ϕ := δ -2 ∂ z Φ -∇ζ • ∇Φ z=ζ ,
where Φ is the unique solution to the boundary value problem (4.7)

     ∆Φ + δ -2 ∂ 2 z Φ = 0 in -1 + b(x) < z < ζ(x), Φ = ϕ on z = ζ(x), ∇b • ∇Φ -δ -2 ∂ z Φ = 0 on z = -1 + b(x).
More precisely, we have the following Lemma. Lemma 4.2. Let c, M be positive constants and m, j integers such that m > n 2 +1, m ≥ 2(j +1), and 1 ≤ j ≤ 2N + 1. We assume (H1) or (H2). There exists a positive constant

C such that if ζ ∈ H m , b ∈ W m+1,∞ , and H = 1 + ζ -b satisfy (4.8) ∥ζ∥ H m + ∥b∥ W m+1,∞ ≤ M, H(x) ≥ c for x ∈ R n ,
then for any ϕ ∈ Hk+2(j+1) with 0 ≤ k ≤ m -2(j + 1) and any δ ∈ (0, 1] we have

∥Λ (N ) (ζ, b, δ)ϕ -Λ(ζ, b, δ)ϕ∥ H k ≤ Cδ 2j ∥∇ϕ∥ H k+2j+1 .
Proof. We observe that the bound on The proof is also valid in the case 1 ≤ j ≤ 2N + 1 and 0 ≤ k ≤ m -2(j + 1).

r 1 := Λ (N ) (ζ, b, δ)ϕ -Λ(ζ
The above estimate allows us to obtain the desired consistency result on the equations describing the conservation of mass. We need a similar estimate for the contributions of Bernoulli's equation. To this end, we denote (4.9)

B(ϕ; ζ, b, δ) := 1 2 |∇ϕ| 2 - 1 2 δ 2 (Λ(ζ, b, δ)ϕ + ∇ζ • ∇ϕ) 2 1 + δ 2 |∇ζ| 2 and (4.10) B (N ) (ϕ; ζ, b, δ) := 1 2 |u| 2 + δ -2 w 2 -wΛ (N ) (ζ, b, δ)ϕ with u := (l(H) ⊗ ∇) T ϕ -(l ′ (H) • ϕ)∇b, w := l ′ (H) • ϕ,
where l ′ (H) := (0, p 1 H p 1 -1 , . . . , p N * H p N * -1 ) T and ϕ := (ϕ 0 , ϕ 1 , . . . , ϕ N * ) T is the solution to (4.5), whose unique existence is guaranteed by Lemma 4.1. Then, the following lemma shows that B (N ) (ϕ; ζ, b, δ) is a higher order approximation to B(ϕ; ζ, b, δ) in the shallow water regime δ ≪ 1.

∥B (N ) (ϕ; ζ, b, δ) -B(ϕ; ζ, b, δ)∥ H m-4(N +1) ≤ Cδ 4N +2 ∥∇ϕ∥ 2 H m-1 .
Proof. Notice first that differentiating ϕ = l(H) • ϕ we have ∇ϕ = u + w∇ζ, so that

B (N ) (ϕ; ζ, b, δ) = 1 2 |∇ϕ| 2 + δ -2 w 2 (1 + δ 2 |∇ζ| 2 ) -w ∇ζ • ∇ϕ + Λ (N ) (ζ, b, δ)ϕ = 1 2 |∇ϕ| 2 + δ -2 w 2 (1 + δ 2 |∇ζ| 2 ) -w Λ(ζ, b, δ)ϕ + ∇ζ • ∇ϕ + w Λ(ζ, b, δ)ϕ -Λ (N ) (ζ, b, δ)ϕ .
If we introduce a residual r by

r = (δ -2 ∂ z Φ app -∇ζ • ∇Φ app )| z=ζ -(δ -2 ∂ z Φ -∇ζ • ∇Φ)| z=ζ ,
where Φ is the solution to the boundary value problem (4.7) and Φ app is an approximate velocity potential defined by

Φ app (x, z) = N * i=0 (z + 1 -b(x)) p i ϕ i (x), then we have r = δ -2 w -∇ζ • u -Λ(ζ, b, δ)ϕ = δ -2 w(1 + δ 2 |∇ζ| 2 ) -∇ζ • ∇ϕ -Λ(ζ, b, δ)ϕ.
Therefore, we obtain

B (N ) (ϕ; ζ, b, δ) -B(ϕ; ζ, b, δ) = 1 2 δ 2 r 2 1 + δ 2 |∇ζ| 2 + w Λ(ζ, b, δ)ϕ -Λ (N ) (ζ, b, δ)ϕ .
The desired estimate for the second term readily follows from Lemma

Results in the framework of interfacial waves

In this section, we prove Theorems 3.5 and 3.6. To this end, we first rewrite the Kakinuma model (2.18) using a formulation which allows a direct comparison with the full model for interfacial gravity waves (2.17), thanks to the following Lemma. Lemma 4.4. Let c, M be positive constants and m an integer such that m > n 2 + 1. There exists a positive constant C such that for any positive parameters

h 1 , h 2 , δ satisfying h 1 δ, h 2 δ ≤ 1, if ζ ∈ H m , b ∈ W m,∞ , H 1 = 1 -h -1 1 ζ, and H 2 = 1 + h -1 2 ζ -h -1 2 b satisfy (4.11) h -1 1 ∥ζ∥ H m + h -1 2 ∥ζ∥ H m + h -1 2 ∥b∥ W m,∞ ≤ M, H 1 (x) ≥ c, H 2 (x) ≥ c for x ∈ R n ,
then for any k = 0, ±1, . . . , ±(m -1) and any ϕ 1 , ϕ 2 ∈ Hk+1 there exists a unique solution

ϕ 1 = (ϕ 1,0 , ϕ ′ 1 ) ∈ Hk+1 × (H k+1 ) N , ϕ 2 = (ϕ 2,0 , ϕ ′ 2 ) ∈ Hk+1 × (H k+1 ) N * to the problem (4.12) l 1 (H 1 ) • ϕ 1 = ϕ 1 , L 1,i (H 1 , δ, h 1 )ϕ 1 = 0 for i = 1, 2, . . . , N, l 2 (H 2 ) • ϕ 2 = ϕ 2 , L 2,i (H 2 , b, δ, h 2 )ϕ 2 = 0 for i = 1, 2, . . . , N * .
Moreover, the solution satisfies

∥∇ϕ ℓ ∥ H k + (h ℓ δ) -1 ∥ϕ ′ ℓ ∥ H k ≤ C∥∇ϕ ℓ ∥ H k for ℓ = 1, 2.
Proof. Notice that we have identities

L 1,ij (H 1 , δ, h 1 ) = L ij (H 1 , 0, h 1 δ), L 2,ij (H 2 , b, δ, h 2 ) = L ij (H 2 , h -1 2 b, h 2 δ)
with suitable choices of indices {p i }. Hence, Lemma 4.1 gives the desired result.

As a corollary of this lemma, under the assumptions of Lemma 4.4

Λ (N ) 1 (ζ, δ, h 1 ) : ϕ 1 → L 1,0 (H 1 , h 1 , δ)ϕ 1 , Λ (N ) 2 (ζ, b, δ, h 2 ) : ϕ 2 → L 2,0 (H 2 , b, h 2 , δ)ϕ 2 ,
where (ϕ 1 , ϕ 2 ) is the unique solution to (4.12), are defined as bounded linear operators from Hk+1 to H k-1 for any k = ±0, . . . , ±(m -1). Using these definitions and noting the relations (3.8) and l ℓ (H ℓ )

• ∂ t ϕ ℓ = ∂ t (l ℓ (H ℓ ) • ϕ ℓ ) -w ℓ h -1 ℓ ∂ t ζ
           ∂ t ζ + h 1 Λ (N ) 1 (ζ, δ, h 1 )ϕ 1 = 0, ∂ t ζ -h 2 Λ (N ) 1 (ζ, b, δ, h 2 )ϕ 2 = 0, ρ 1 ∂ t ϕ 1 + 1 2 |u 1 | 2 + (h 1 δ) -2 w 2 1 + w 1 Λ (N ) 1 (ζ, δ, h 1 )ϕ 1 -ρ 2 ∂ t ϕ 2 + 1 2 |u 2 | 2 + (h 2 δ) -2 w 2 2 -w 2 Λ (N ) 2 (ζ, b, δ, h 2 )ϕ 2 -ζ = 0,
where we recall that u 1 , u 2 , w 1 , and w 2 are uniquely determined from ϕ 1 and ϕ 2 by (3.4), wherein ϕ 1 and ϕ 2 are defined as the solutions to (4.12). We further introduce notations, which are contributions of Bernoulli's equation and interfacial versions of B and B (N ) defined by (4.9) and (4.10). We denote

       B 1 (ϕ 1 ; ζ, δ, h 1 ) := 1 2 |∇ϕ 1 | 2 - 1 2 δ 2 (Λ 1 (ζ, δ, h 1 )ϕ 1 -∇ζ • ∇ϕ 1 ) 2 1 + δ 2 |∇ζ| 2 , B 2 (ϕ 2 ; ζ, b, δ, h 2 ) := 1 2 |∇ϕ 2 | 2 - 1 2 δ 2 (Λ 2 (ζ, b, δ, h 2 )ϕ 2 + ∇ζ • ∇ϕ 2 ) 2 1 + δ 2 |∇ζ| 2 , and      B (N ) 1 (ϕ 1 ; ζ, δ, h 1 ) := 1 2 |u 1 | 2 + (h 1 δ) -2 w 2 1 + w 1 Λ (N ) 1 (ζ, δ, h 1 )ϕ 1 , B (N ) 2 (ϕ 2 ; ζ, b, δ, h 2 ) := 1 2 |u 2 | 2 + (h 2 δ) -2 w 2 2 -w 2 Λ (N ) 2 (ζ, b, δ, h 2 )ϕ 2 .
Then, the full model for interfacial gravity waves (2.17) and the Kakinuma model (4.13) can be written simply as

     ∂ t ζ + Λ 1 (ζ, δ, h 1 )ϕ 1 = 0, ∂ t ζ -Λ 2 (ζ, b, δ, h 2 )ϕ 2 = 0, ρ 1 ∂ t ϕ 1 + B 1 (ϕ 1 ; ζ, δ, h 1 ) -ρ 2 ∂ t ϕ 2 + B 2 (ϕ 2 ; ζ, b, δ, h 2 ) -ζ = 0, and      ∂ t ζ + h 1 Λ (N ) 1 (ζ, δ, h 1 )ϕ 1 = 0, ∂ t ζ -h 2 Λ (N ) 1 (ζ, b, δ, h 2 )ϕ 2 = 0, ρ 1 ∂ t ϕ 1 + B (N ) 1 (ϕ 1 ; ζ, δ, h 1 ) -ρ 2 ∂ t ϕ 2 + B (N ) 2 (ϕ 2 ; ζ, b, δ, h 2 ) -ζ = 0, respectively. The following lemmas show that h 1 Λ (N ) 1 , h 2 Λ (N ) 2 , B (N )
1 , and B (N ) 2 are higher order approximations in the shallow water regime δ 1 = h 1 δ ≪ 1 and δ 2 = h 2 δ ≪ 1 to Λ 1 , Λ 2 , B 1 , and B 2 , respectively. Lemma 4.5. Let c, M be positive constants and m, j integers such that m > n 2 +1, m ≥ 2(j +1), and 1 ≤ j ≤ 2N + 1. We assume (H1) or (H2). There exists a positive constant C such that for any positive parameters

h 1 , h 2 , δ satisfying h 1 δ, h 2 δ ≤ 1, if ζ ∈ H m , b ∈ W m+1,∞ , H 1 = 1-h -1 1 ζ, and H 2 = 1 + h -1 2 ζ -h -1 2 b satisfy (4.14) h -1 1 ∥ζ∥ H m + h -1 2 ∥ζ∥ H m + h -1 2 ∥b∥ W m+1,∞ ≤ M, H 1 (x) ≥ c, H 2 (x) ≥ c for x ∈ R n , Lemma 5.1.
Let c, M be positive constants and m an integer such that m > n 2 +1. There exists a positive constant C such that for any positive parameters ρ

1 , ρ 2 , h 1 , h 2 , δ satisfying h 1 δ, h 2 δ ≤ 1, if ζ ∈ H m , b ∈ W m,∞ , H 1 = 1 -h -1 1 ζ, and H 2 = 1 + h -1 2 ζ -h -1 2 b satisfy (4.11), then for any f ′ 1 = (f 1,1 , . . . , f 1,N ) T ∈ (H k ) N , f ′ 2 = (f 2,1 , . . . , f 2,N * ) T ∈ (H k ) N * , f 3 ∈ (H k ) n
, and f 4 ∈ Hk+1 with k ∈ {0, 1, . . . , m -1}, there exists a solution (φ 1 , φ 2 ) to (5.1)

           L 1,i (H 1 , δ, h 1 )φ 1 = f 1,i for i = 1, 2, . . . , N, L 2,i (H 2 , b, δ, h 2 )φ 2 = f 2,i for i = 1, 2, . . . , N * , h 1 L 1,0 (H 1 δ, h 1 )φ 1 + h 2 L 2,0 (H 2 , b, δ, h 2 )φ 2 = ∇ • f 3 , -ρ 1 l 1 (H 1 ) • φ 1 + ρ 2 l 2 (H 2 ) • φ 2 = f 4 , satisfying ℓ=1,2 ρ ℓ h ℓ ∥∇φ ℓ ∥ 2 H k + (h ℓ δ) -2 ∥φ ′ ℓ ∥ 2 H k ≤ C ℓ=1,2 ρ ℓ h ℓ min ∥f ′ ℓ ∥ 2 H k-1 , (h ℓ δ) 2 ∥f ′ ℓ ∥ 2 H k + min ρ 1 h 1 , ρ 2 h 2 ∥f 3 ∥ 2 H k + min h 1 ρ 1 , h 2 ρ 2 ∥∇f 4 ∥ 2 H k .
Moreover, the solution is unique up to an additive constant of the form (Cρ 2 , Cρ 1 ) to (φ 1,0 , φ 2,0 ).

Proof. The existence and uniqueness up to an additive constant of the solution has been given in the companion paper [START_REF] Duchêne | A mathematical analysis of the Kakinuma model for interfacial gravity waves. Part I: Structures and well-posedness[END_REF]Lemma 6.4]. We focus here on the derivation of uniform estimates. By direct rescaling within the proof of [14, Lemma 6.1], we infer that

(L ℓ φ ℓ , φ ℓ ) L 2 ≃ ∥∇φ ℓ ∥ 2 L 2 + (h ℓ δ) -2 ∥φ ′ ℓ ∥ 2 L 2
for ℓ = 1, 2. We note the identities

L 1 φ 1 = l 1 L 1,0 φ 1 + (0, L 1,1 φ 1 , . . . , L 1,N φ 1 ) T , L 2 φ 2 = l 2 L 2,0 φ 2 + (0, L 2,1 φ 2 , . . . , L 2,N * φ 2 ) T ,
so that for the solution (φ 1 , φ 2 ) to (5.1) we have ℓ=1,2

ρ ℓ h ℓ (L ℓ φ ℓ , φ ℓ ) L 2 = ℓ=1,2 ρ ℓ h ℓ (L ℓ,0 φ ℓ , l ℓ • φ ℓ ) L 2 + ℓ=1,2 ρ ℓ h ℓ (f ′ ℓ , φ ′ ℓ ) L 2 (5.2) =: I 1 + I 2 .
Therefore, it is sufficient to evaluate I 1 and I 2 . As for the term I 2 we have

|(f ′ ℓ , φ ′ ℓ ) L 2 | ≤ min{∥f ′ ℓ ∥ H -1 ∥φ ′ ℓ ∥ H 1 , ∥f ′ ℓ ∥ L 2 ∥φ ′ ℓ ∥ L 2 } ≤ min{∥f ′ ℓ ∥ H -1 , (h ℓ δ)∥f ′ ℓ ∥ L 2 }(∥∇φ ℓ ∥ L 2 + (h ℓ δ) -1 ∥φ ′ ℓ ∥ L 2 ).
As for the term I 1 , we note the trivial identities ℓ=1,2

ρ ℓ h ℓ (L ℓ,0 φ ℓ , l ℓ • φ ℓ ) L 2 = (h 1 L 1,0 φ 1 + h 2 L 2,0 φ 2 , ρ 1 l 1 • φ 1 ) L 2 + (h 2 L 2,0 φ 2 , ρ 2 l 2 • φ 2 -ρ 1 l 1 • φ 1 ) L 2 , (h 1 L 1,0 φ 1 + h 2 L 2,0 φ 2 , ρ 2 l 2 • φ 2 ) L 2 + (h 1 L 1,0 φ 1 , ρ 1 l 1 • φ 1 -ρ 2 l 2 • φ 2 ) L 2 .
Therefore, the term I 1 in (5.2) can be expressed in two ways as

I 1 = ρ 1 (∇ • f 3 , l 1 • φ 1 ) L 2 + h 2 (L 2,0 φ 2 , f 4 ) L 2 , ρ 2 (∇ • f 3 , l 2 • φ 2 ) L 2 -h 1 (L 1,0 φ 1 , f 4 ) L 2 .
By the linearity of (5.1) it is sufficient to evaluate it in the case f 4 = 0 and in the case f 3 = 0, separately. In the case f 4 = 0, we evaluate it as

|I 1 | ≤ min{ρ 1 ∥f 3 ∥ L 2 ∥∇(l 1 • φ 1 )∥ L 2 , ρ 2 ∥f 3 ∥ L 2 ∥∇(l 2 • φ 2 )∥ L 2 } = min ρ 1 h 1 ∥f 3 ∥ L 2 ρ 1 h 1 ∥∇(l 1 • φ 1 )∥ L 2 , ρ 2 h 2 ∥f 3 ∥ L 2 ρ 2 h 2 ∥∇(l 2 • φ 2 )∥ L 2 ≲ min ρ 1 h 1 , ρ 2 h 2 ∥f 3 ∥ L 2 ℓ=1,2 ρ ℓ h ℓ (∥∇φ ℓ ∥ L 2 + ∥φ ′ ℓ ∥ L 2 ).
In the case f 3 = 0 we evaluate it as

|I 1 | ≲ min{h 1 ∥∇φ 1 ∥ L 2 ∥∇f 4 ∥ L 2 , h 2 (∥∇φ 2 ∥ L 2 + ∥φ ′ 2 ∥ L 2 )∥∇f 4 ∥ L 2 } = min h 1 ρ 1 ∥∇f 4 ∥ L 2 ρ 1 h 1 ∥∇φ 1 ∥ L 2 , h 2 ρ 2 ∥∇f 4 ∥ L 2 ρ 2 h 2 (∥∇φ 1 ∥ L 2 + ∥φ ′ 2 ∥ L 2 ) ≤ min h 1 ρ 1 , h 2 ρ 2 ∥∇f 4 ∥ L 2 ℓ=1,2 ρ ℓ h ℓ (∥∇φ ℓ ∥ L 2 + ∥φ ′ ℓ ∥ L 2 ).
From the above estimates we deduce immediately the desired inequality for k = 0. In order to obtain the desired inequality on derivatives, we let k ∈ {1, 2, . . . , m -1} and β be a multi-index such that 1 ≤ |β| ≤ k. Applying the differential operator ∂ β to (5.1), we have

           L 1,i ∂ β φ 1 = ∂ β f 1,i + f 1,i,β for i = 1, 2, . . . , N, L 2,i ∂ β φ 2 = ∂ β f 2,i + f 2,i,β for i = 1, 2, . . . , N * , h 1 L 1,0 ∂ β φ 1 + h 2 L 2,0 ∂ β φ 2 = ∇ • (∂ β f 3 + h 1 f 3,1,β + h 2 f 3,2,β ), -ρ 1 l 1 • ∂ β φ 1 + ρ 2 l 2 • ∂ β φ 2 = ∂ β f 4 + ρ 1 f 4,1,β + ρ 2 f 4,2,β , where                      f 1,i,β := -[∂ β , L 1,i (H 1 , δ, h 1 )]φ 1 for i = 1, 2, . . . , N, f 2,i,β := -[∂ β , L 2,i (H 2 , b, δ, h 2 )]φ 2 for i = 1, 2, . . . , N * , ∇ • f 3,1,β := -[∂ β , L 1,0 (H 1 , δ, h 1 )]φ 1 , ∇ • f 3,2,β := -[∂ β , L 2,0 (H 2 , b, δ, h 2 )]φ 2 , f 4,1,β := [∂ β , l 1 (H 1 )] • φ 1 , f 4,2,β := -[∂ β , l 2 (H 2 )] • φ 2 .
We put f 1,β = (0, f 1,1,β , . . . , f 1,N,β ) and f 2,β = (0, f 2,1,β , . . . , f 2,N * ,β ). Then, with a suitable decomposition f ℓ,β = f high ℓ,β + f low ℓ,β for ℓ = 1, 2, we see that

∥f high ℓ,β ∥ H -1 + (h ℓ δ)∥f low ℓ,β ∥ L 2 + ∥f 3,ℓ,β ∥ L 2 + ∥∇f 4,ℓ ∥ L 2 ≲ ∥∇φ ℓ ∥ H k-1 + (h ℓ δ) -1 ∥φ ′ ℓ ∥ H k-1
for ℓ = 1, 2. Therefore, in view of the linearity of (5.1) the desired inequality for k ≥ 1 follows by induction on k.

From the above elliptic estimates we deduce the following bounds on time derivatives of regular solutions to the Kakinuma model (2.18). We introduce a mathematical energy E m (t) for a solution (ζ, ϕ 1 , ϕ 2 ) to the Kakinuma model by

(5.3) E m (t) := ∥ζ(t)∥ 2 H m + ℓ=1,2 ρ ℓ h ℓ (∥∇ϕ ℓ (t)∥ 2 H m + (h ℓ δ) -2 ∥ϕ ′ ℓ (t)∥ 2 H m ),
where ϕ ′ 1 = (ϕ 1,1 , . . . , ϕ 1,N ) T and ϕ ′ 2 = (ϕ 2,1 , . . . , ϕ 2,N * ) T .

Lemma 5.2. Let c, M 1 , h min be positive constants and m an integer such that m > n 2 + 1. There exists a positive constant C 1 such that for any positive parameters ρ 1 , ρ 2 , h 1 , h 2 , δ satisfying the natural restrictions (2.14), h 1 δ, h 2 δ ≤ 1, and the condition

h min ≤ h 1 , h 2 , if a regular solution (ζ, ϕ 1 , ϕ 2 ) to the Kakinuma model (2.18) with bottom topography b ∈ W m+1,∞ satisfy E m (t) + h -1 2 ∥b∥ W m+1,∞ ≤ M 1 , H 1 (x, t) ≥ c, H 2 (x, t) ≥ c for x ∈ R n , 0 ≤ t ≤ T, then we have ∥∂ t ζ(t)∥ 2 H m-1 + ℓ=1,2 ρ ℓ h ℓ (∥∇∂ t ϕ ℓ (t)∥ 2 H m-1 + (h ℓ δ) -2 ∥∂ t ϕ ′ ℓ (t)∥ 2 H m-1 ) (5.4) + ∥∂ 2 t ζ(t)∥ 2 H m-2 + ℓ=1,2 ρ ℓ h ℓ (∥∇∂ 2 t ϕ ℓ (t)∥ 2 H m-2 + (h ℓ δ) -2 ∥∂ 2 t ϕ ′ ℓ (t)∥ 2 H m-2 ) ≤ C E m (t)
for 0 ≤ t ≤ T .

Proof. First, we recall that the Kakinuma model (2.18) can be written compactly as (3.5). It follows from the first component of the first two equations in (3.5) that ∂ t ζ can be written in two ways as

∂ t ζ = -h 1 L 1,0 ϕ 1 = h 2 L 2,0 ϕ 2 , so that ∥∂ t ζ∥ 2 H m-1 = min{h 2 1 ∥L 1,0 ϕ 1 ∥ 2 H m-1 , h 2 2 ∥L 2,0 ϕ 2 ∥ 2 H m-1 } ≲ min{h 2 1 ∥∇ϕ 1 ∥ 2 H m , h 2 2 (∥∇ϕ 2 ∥ 2 H m + ∥ϕ ′ 2 ∥ 2 H m )} ≤ min h 1 ρ 1 , h 2 ρ 2 E m ≤ 2E m ,
where we used (2.15). As for the estimate of (∂ t ϕ 1 , ∂ t ϕ 2 ), we differentiate the compatibility conditions (3.7) with respect to time and use the last equation in (3.5). Then, we have (5.5)

           L 1,i ∂ t ϕ 1 = f 1,i for i = 1, 2, . . . , N, L 2,i ∂ t ϕ 2 = f 2,i for i = 1, 2, . . . , N * , h 1 L 1,0 ∂ t ϕ 1 + h 2 L 2,0 ∂ t ϕ 2 = ∇ • f 3 , -ρ 1 l 1 • ∂ t ϕ 1 + ρ 2 l 2 • ∂ t ϕ 2 = f 4 , where (5.6)            f 1,i := -[∂ t , L 1,i (H 1 , δ, h 1 )]ϕ 1 for i = 1, 2, . . . , N, f 2,i := -[∂ t , L 2,i (H 2 , b, δ, h 2 )]ϕ 2 for i = 1, 2, . . . , N * , f 3 := (u 2 -u 1 )∂ t ζ, f 4 := 1 2 ρ 1 |u 1 | 2 + (h 1 δ) -2 w 2 1 -1 2 ρ 2 |u 2 | 2 + (h 2 δ) -2 w 2 2 -ζ.
Therefore, by Lemma 5.1 we have ℓ=1,2

ρ ℓ h ℓ (∥∇∂ t ϕ ℓ ∥ 2 H m-1 + (h ℓ δ) -2 ∥∂ t ϕ ′ ℓ ∥ 2 H m-1 ) (5.7) ≲ ℓ=1,2 ρ ℓ h ℓ (h ℓ δ) 2 ∥f ′ ℓ ∥ 2 H m-1 + min ρ 1 h 1 , ρ 2 h 2 ∥f 3 ∥ 2 H m-1 + ∥f 4 ∥ 2 H m , where f ′ 1 = (f 1,1 , . . . , f 1,N ) T , f ′ 2 = (f 2,1 , . . . , f 2,N * )
T , and we used (2.15). We proceed to evaluate the right-hand side. By writing down the operators L ℓ,i explicitly, we see that the operators do not include any derivatives of H ℓ . Therefore, we can write f ℓ,i as

f 1,i = ∂ ∂H 1 L 1,i ϕ 1 h -1 1 ∂ t ζ, f 2,i = - ∂ ∂H 2 L 2,i ϕ 2 h -1 2 ∂ t ζ.
We note also that the operators ∂ ∂H ℓ L ℓ,i have a similar structure as L ℓ,i . Therefore,

ρ ℓ h ℓ (h ℓ δ) 2 ∥f ′ ℓ ∥ 2 H m-1 ≲ ρ ℓ h ℓ (h ℓ δ) 2 (∥∇ϕ ℓ ∥ 2 H m + (h ℓ δ) -4 ∥ϕ ′ ℓ ∥ 2 H m-1 )∥h -1 ℓ ∂ t ζ∥ 2 H m-1 ≲ E 2 m for ℓ = 1, 2,
where, here and henceforth, we utilize fully our restriction

h -1 1 , h -1 2 ≲ 1.
In view of the definition (3.4) of u 1 , u 2 , w 1 , and w 2 , we see easily that

(5.8) ℓ=1,2 ρ ℓ h ℓ (∥u ℓ ∥ 2 H m + (h ℓ δ) -2 ∥w ℓ ∥ 2 H m ) ≲ E m .
We evaluate the term on f 3 as

min ρ 1 h 1 , ρ 2 h 2 ∥f 3 ∥ 2 H m-1 ≲ ℓ=1,2 ρ ℓ h ℓ ∥u ℓ ∂ t ζ∥ 2 H m-1 ≲ ℓ=1,2 ρ ℓ h ℓ ∥u ℓ ∥ 2 H m-1 ∥h -1 ℓ ∂ t ζ∥ 2 H m-1 ≲ E 2 m .
Similarly, we have

∥f 4 ∥ 2 H m ≲ ℓ=1,2 ρ 2 ℓ (∥u ℓ ∥ 2 H m + (h 1 δ) -2 ∥w ℓ ∥ 2 H m ) 2 + ∥ζ∥ 2 H m ≲ ℓ=1,2 h -2 ℓ {ρ ℓ h ℓ (∥u ℓ ∥ 2 H m + (h 1 δ) -2 ∥w ℓ ∥ 2 H m )} 2 + ∥ζ∥ 2 H m ≲ E 2 m + E m .
Plugging in (5.7) the above estimates, we obtain the desired estimate for (∂ t ϕ 1 , ∂ t ϕ 2 ). Finally, the estimate of ∂ 2 t ζ can be obtained by differentiating

∂ t ζ = -h 1 L 1,0 ϕ 1 = h 2 L 2,
0 ϕ 2 with respect to time. Then, the estimate of (∂ 2 t ϕ 1 , ∂ 2 t ϕ 2 ) can be obtained by differentiating (5.5) with respect to time once more and applying Lemma 5.1.

Remark 5.3. In view of the above arguments, we see easily that for the Kakinuma model (2.18), (∂ t ϕ 1 , ∂ t ϕ 2 )| t=0 can be determined from the initial data (ζ (0) , ϕ 1(0) , ϕ 2(0) ) and the bottom topography b, although the hypersurface t = 0 is characteristic for the model. They are unique up to an additive constant of the form (Cρ 2 , Cρ 1 ) to

(∂ t ϕ 1,0 , ∂ t ϕ 2,0 )| t=0 . Particularly, (∂ t ϕ ′ 1 , ∂ t ϕ ′ 2 )
| t=0 and hence a| t=0 with the function a given in (3.9) can be uniquely determined from the data.

Uniform energy estimates; proof of Theorem 3.1

In this section we provide uniform energy estimates for solutions to the Kakinuma model. Consequently, we prove Theorem 3.1. We recall that the Kakinuma model (2.18) can be written compactly as (6.1)

           l 1 (H 1 )∂ t ζ + h 1 L 1 (H 1 , δ, h 1 )ϕ 1 = 0, l 2 (H 2 )∂ t ζ -h 2 L 2 (H 2 , b, δ, h 2 )ϕ 2 = 0, ρ 1 l 1 (H 1 ) • ∂ t ϕ 1 + 1 2 |u 1 | 2 + (h 1 δ) -2 w 2 1 -ρ 2 l 2 (H 2 ) • ∂ t ϕ 2 + 1 2 |u 2 | 2 + (h 2 δ) -2 w 2 2 -ζ = 0,
where we recall that

H 1 := 1 -h -1 1 ζ, H 2 := 1 + h -1 2 ζ -h -1 2 b, ϕ 1 := (ϕ 1,0 , ϕ 1,1 , . . . , ϕ 1,N ) T , ϕ 2 := (ϕ 2,0 , ϕ 2,1 , . . . , ϕ 2,N * ) T , and l 1 , l 2 , L 1 , L 2 , u 1 , u 2 , w 1 , w 2 are defined in Section 3.

Analysis of linearized equations

Before deriving linearized equations to the Kakinuma model (6.1), we introduce some more notations. For ℓ = 1, 2, the coefficient matrices of the principal part and the singular part with respect to the small parameter δ ℓ = h ℓ δ of the operator L ℓ are denoted by A ℓ (H ℓ ) and C ℓ (H ℓ ), respectively, that is, (6.2)

       A 1 (H 1 ) := 1 2(i + j) + 1 H 2(i+j)+1 1 0≤i,j≤N , A 2 (H 2 ) := 1 p i + p j + 1 H p i +p j +1 2 0≤i,j≤N *
, and (6.3)

       C 1 (H 1 ) := 4ij 2(i + j) -1 H 2(i+j)-1 1 0≤i,j≤N , C 2 (H 2 ) := p i p j p i + p j -1 H p i +p j -1 2 0≤i.j≤N * .
We put also (6.4)

         B 2 (H 2 ) := p j p i + p j H p i +p j 2 0≤i,j≤N * , B2 (H 2 ) := B 2 (H 2 ) -B 2 (H 2 ) T , C2 (H 2 , h -1 2 b) := |h -1 2 ∇b| 2 C 2 (H 2 ) + h -1 2 (∆b)B 2 (H 2 ).
In the above expressions, we used the notational convention 0/0 = 0. Then, the operators L 1 and L 2 can also be written as (6.5)

L 1 ϕ 1 = -A 1 ∆ϕ 1 -l 1 (u 1 • ∇H 1 ) + (h 1 δ) -2 C 1 ϕ 1 , L 2 ϕ 2 = -A 2 ∆ϕ 2 -l 2 (u 2 • ∇H 2 ) + (h 2 δ) -2 C 2 ϕ 2 + B2 (h -1 2 ∇b • ∇)ϕ 2 + C2 ϕ 2 .
For ℓ = 1, 2, we decompose the operator L ℓ as L ℓ = L pr ℓ + L low ℓ , where (6.6)

L pr ℓ (H ℓ )φ ℓ := - n l=1 ∂ l (A ℓ (H ℓ )∂ l φ ℓ ) + (h ℓ δ) -2 C ℓ (H ℓ )φ ℓ .
We now linearize the Kakinuma model (6.1) around an arbitrary flow (ζ, ϕ 1 , ϕ 2 ) and denote the variation by ( ζ, φ1 , φ2 ). After neglecting lower order terms, the linearized equations have the form (6.7)

     l 1 (H 1 )(∂ t + u 1 • ∇) ζ + h 1 L pr 1 (H 1 , δ, h 1 ) φ1 = ḟ1 , l 2 (H 2 )(∂ t + u 2 • ∇) ζ -h 2 L pr 2 (H 2 , δ, h 2 ) φ2 = ḟ2 , ρ 1 l 1 (H 1 ) • (∂ t + u 1 • ∇) φ1 -ρ 2 l 2 (H 2 ) • (∂ t + u 1 • ∇) φ2 -a ζ = ḟ0 ,
where the function a is defined by (3.9). In order to derive a good symmetric structure of the equations, following the companion paper [START_REF] Duchêne | A mathematical analysis of the Kakinuma model for interfacial gravity waves. Part I: Structures and well-posedness[END_REF] we introduce (6.8)

θ 1 := ρ 2 h 1 H 1 α 1 ρ 1 h 2 H 2 α 2 + ρ 2 h 1 H 1 α 1 , θ 2 := ρ 1 h 2 H 2 α 2 ρ 1 h 2 H 2 α 2 + ρ 2 h 1 H 1 α 1 ,
where (6.9)

α ℓ := det A ℓ,0 det Ãℓ,0 , Ãℓ,0 := 0 1 T -1 A ℓ,0 , A ℓ,0 := A ℓ (1)
for ℓ = 1, 2 and 1 := (1, . . . , 1) T . Then, we have θ 1 + θ 2 = 1. We recall that α 1 and α 2 are positive constants depending only on N and the nonnegative integers 0 = p 0 < p 1 < . . . < p N * , respectively, and go to 0 as N, N * → ∞. We also introduce

u := θ 2 u 1 + θ 1 u 2 , v := u 2 -u 1 .
Then, we have u 1 = u-θ 1 v and u 2 = u+θ 2 v. Plugging these into the linearized equations (6.7), we can write them in a matrix form as (6.10)

A 1 (∂ t + u • ∇) U + A mod 0 U = Ḟ ,
where

U :=   ζ φ1 φ2   , Ḟ :=    ḟ0 ρ 1 ( ḟ1 -(∇ • (θ 1 l 1 ⊗ v)) ζ ρ 2 ( ḟ2 -(∇ • (θ 2 l 2 ⊗ v)) ζ   ,
and

A 1 :=   0 -ρ 1 l T 1 ρ 2 l T 2 ρ 1 l 1 O O -ρ 2 l 2 O O   , A mod 0 :=   a ρ 1 θ 1 l T 1 (v • ∇) ρ 2 θ 2 l T 2 (v • ∇) (v • ∇) * (ρ 1 θ 1 l 1 • ) ρ 1 h 1 L pr 1 O (v • ∇) * (ρ 2 θ 1 l 2 • ) O ρ 2 h 2 L pr 2   .
Here, (v

• ∇) * denotes the adjoint operator of v • ∇ in L 2 , that is, (v • ∇) * f = -∇ • (f v).
We note that A 1 is a skew-symmetric matrix and A mod 0 is symmetric in L 2 . Therefore, the corresponding energy function is given by (A mod 0 U , U ) L 2 . We put (6.11)

E ( U ) := ∥ ζ∥ 2 L 2 + ℓ=1,2 ρ ℓ h ℓ (∥∇ φℓ ∥ 2 L 2 + (h ℓ δ) -2 ∥ φ′ ℓ ∥ 2 L 2 ).
The following lemma shows that (A mod 0 U , U ) L 2 ≃ E ( U ) under the non-cavitation assumption and the stability condition, stated respectively as (3.11) and (3.12) in Theorem 3.1. Lemma 6.1. Let c, M, h min be positive constants. There exists a positive constant C such that for any positive parameters ρ 1 , ρ 2 , h 1 , h 2 , δ satisfying the condition

h min ≤ h 1 , h 2 , if H 1 , H 2 , u 1 , u 2 ,
and the function a satisfy (6.12)

             ℓ=1,2 ∥H ℓ ∥ L ∞ + ρ ℓ h ℓ ∥u ℓ ∥ L ∞ + ∥a∥ L ∞ ≤ M, a(x) - ρ 1 ρ 2 ρ 1 h 2 H 2 (x)α 2 + ρ 2 h 1 H 1 (x)α 1 |u 2 (x) -u 1 (x)| 2 ≥ c, H 1 (x) ≥ c, H 2 (x) ≥ c for x ∈ R n , then for any U = ( ζ, φ1 , φ2 ) T ∈ L 2 × ( H1 × (H 1 ) N ) × ( H1 × (H 1 ) N * ) we have C -1 E ( U ) ≤ (A mod 0 U , U ) L 2 ≤ CE ( U ).
Proof. This lemma can be shown along with the proof of [START_REF] Duchêne | A mathematical analysis of the Kakinuma model for interfacial gravity waves. Part I: Structures and well-posedness[END_REF]Lemma 7.4]. For the sake of completeness, we sketch the proof. We first note that

(A mod 0 U , U ) L 2 = (a ζ, ζ) L 2 + ℓ=1,2 {ρ ℓ h ℓ (L pr ℓ φℓ , φℓ ) L 2 + 2ρ ℓ (θ ℓ l ℓ • (v • ∇) φℓ , ζ) L 2 } = (a ζ, ζ) L 2 + ℓ=1,2 ρ ℓ h ℓ n l=1 (A ℓ ∂ l φℓ , ∂ l φℓ ) L 2 + (h ℓ δ) -2 (C ℓ φℓ , φℓ ) L 2 + 2ρ ℓ (θ ℓ v • (l ℓ ⊗ ∇) T φℓ , ζ) L 2 ,
where we used the identity a

• (v • ∇)φ = v • (a ⊗ ∇) T φ.
On the other hand, we can put

q ℓ (H ℓ ) q ℓ (H ℓ ) T -q ℓ (H ℓ ) Q ℓ (H ℓ ) := 0 l ℓ (H ℓ ) T -l ℓ (H ℓ ) A ℓ (H ℓ ) -1
for ℓ = 1, 2. Then, we see that q ℓ (H ℓ ) = H ℓ α ℓ and that Q ℓ (H ℓ ) is nonnegative. Moreover, the identity (6.13)

A ℓ (H ℓ )φ ℓ • φ ℓ = q ℓ (H ℓ )(l ℓ (H ℓ ) • φ ℓ ) 2 + Q ℓ (H ℓ )A ℓ (H ℓ )φ ℓ • A ℓ (H ℓ )φ ℓ holds for any φ ℓ . Therefore, n l=1 (A ℓ ∂ l φℓ , ∂ l φℓ ) L 2 = n l=1 {(q ℓ l ℓ • ∂ l φℓ , l ℓ • ∂ l φℓ ) L 2 + (Q ℓ A ℓ ∂ l φℓ , A ℓ ∂ l φℓ ) L 2 } = (H ℓ α ℓ (l ℓ ⊗ ∇) T φℓ , (l ℓ ⊗ ∇) T φℓ ) L 2 + n l=1 (Q ℓ A ℓ ∂ l φℓ , A ℓ ∂ l φℓ ) L 2 , so that (A mod 0 U , U ) L 2 = (a ζ, ζ) L 2 + ℓ=1,2 {ρ ℓ h ℓ (H ℓ α ℓ (l ℓ ⊗ ∇) T φℓ , (l ℓ ⊗ ∇) T φℓ ) L 2 + 2ρ ℓ (θ ℓ v • (l ℓ ⊗ ∇) T φℓ , ζ) L 2 } + ℓ=1,2 ρ ℓ h ℓ n l=1 (Q ℓ A ℓ ∂ l φℓ , A ℓ ∂ l φℓ ) L 2 + (h ℓ δ) -2 (C ℓ φℓ , φℓ ) L 2 =: I 1 + I 2 .
We proceed to evaluate I 1 .

I 1 ≥ R n a ζ2 + ℓ=1,2 ρ ℓ h ℓ H ℓ α ℓ |(l ℓ ⊗ ∇) T φℓ | 2 -2ρ ℓ θ ℓ |v||(l ℓ ⊗ ∇) T φℓ || ζ| dx = R n A 0     ζ ρ 1 h 1 |(l 1 ⊗ ∇) T φ1 | ρ 2 h 2 |(l 2 ⊗ ∇) T φ2 |     •     ζ ρ 1 h 1 |(l 1 ⊗ ∇) T φ1 | ρ 2 h 2 |(l 2 ⊗ ∇) T φ2 |     dx,
where the matrix A 0 is given by

A 0 =      a -ρ 1 /h 1 θ 1 |v| -ρ 2 /h 2 θ 2 |v| -ρ 1 /h 1 θ 1 |v| H 1 α 1 0 -ρ 2 /h 2 θ 2 |v| 0 H 2 α 2     
.

Here, we see that

det A 0 = H 1 H 2 α 1 α 2 a - ρ 1 ρ 2 ρ 1 h 2 H 2 α 2 + ρ 2 h 1 H 1 α 1 |v| 2 ≥ c 3 α 1 α 2 > 0,
so that A 0 is positive definite by Sylvester's criterion. Moreover, tr A 0 ≤ max{1, α 1 , α 2 } M ≲ 1 and the minimal eigenvalue of the matrix A 0 is bounded from below by 4 det A 0 /(tr A 0 ) 2 ≳ 1. Therefore, we obtain

I 1 ≳ R n ζ2 + ℓ=1,2 ρ ℓ h ℓ H ℓ α ℓ |(l ℓ ⊗ ∇) T φℓ | 2 dx.
As for I 2 , it is easy to see that (C ℓ φℓ , φℓ ) L 2 ≃ ∥ φ′ ℓ ∥ 2 L 2 for ℓ = 1, 2. Summarizing the above estimates and using the decomposition (6.13) again, we obtain (A mod 0 U , U ) L 2 ≳ E ( U ). In order to obtain the estimate of (A mod 0 U , U ) L 2 from above, it is sufficient to show that each element of the matrix A 0 is uniformly bounded. Since θ 1 + θ 2 = 1, we have

   ρ 1 /h 1 θ 1 |v| ≤ h -1 1 ρ 1 h 1 |u 1 | + ρ 1 /h 1 θ 1 |u 2 |, ρ 2 /h 2 θ 2 |v| ≤ ρ 2 /h 2 θ 2 |u 1 | + h -1 2 ρ 2 h 2 |u 2 |.
Here, we see that

ρ 1 /h 1 θ 1 |u 2 | = 1 h 2 H 1 α 1 H 2 α 2 (ρ 1 h 2 H 2 α 2 )(ρ 2 h 1 H 1 α 1 ) ρ 1 h 2 H 2 α 2 + ρ 2 h 1 H 1 α 1 ρ 2 h 2 |u 2 | ≤ 1 2h 2 H 1 α 1 H 2 α 2 ρ 2 h 2 |u 2 | ≤ 1 2h min M α 1 cα 2 M ≲ 1.
Similarly, we have

ρ 2 /h 2 θ 2 |u 1 | ≲ 1. Therefore, we obtain (A mod 0 U , U ) L 2 ≲ E ( U ).
In the following Lemma we provide uniform energy estimates for regular solutions to the linearized Kakinuma model (6.7). Proposition 6.2. Let c, M, M 1 , h min be positive constants. There exist positive constants C = C(c, M, h min ) and C 1 = C 1 (c, M, M 1 , h min ) such that for any positive parameters ρ 1 , ρ 2 , h 1 , h 2 , δ satisfying the natural restrictions (2.14) and the condition

h min ≤ h 1 , h 2 , if H 1 , H 2 , u 1 , u 2 ,
and the function a satisfy (6.12) and ℓ=1,2

∥∂ t H ℓ ∥ L ∞ + ∥∇H ℓ ∥ L ∞ + ρ ℓ h ℓ (∥∂ t u ℓ ∥ 2 L ∞ + ∥∇u ℓ ∥ 2 L ∞ ) + ∥∂ t a∥ L ∞ + ∥∇a∥ L ∞ ≤ M 1 ,
then for any regular solution U = ( ζ, φ1 , φ2 ) T to the linearized Kakinuma model (6.7) we have

E ( U (t)) ≤ Ce C 1 t E ( U (0)) + C 1 t 0 e C 1 (t-τ ) ∥ ḟ0 (τ )∥ H 1 (∥∂ t ζ(τ )∥ H -1 + ∥ ζ(τ )∥ L 2 ) + ℓ=1,2 ρ ℓ (∥ ḟℓ (τ )∥ L 2 + ∥ ζ(τ )∥ L 2 )∥(∂ t φℓ (τ ), ∇ φℓ (τ ))∥ L 2 dτ.
Proof. We deduce from (6.10) that

d dt (A mod 0 U , U ) L 2 = ([∂ t , A mod 0 ] U , U ) L 2 + 2(A mod 0 ∂ t U , U ) L 2 = ([∂ t , A mod 0 ] U , U ) L 2 + 2((∂ t + u • ∇) U , A mod 0 U ) L 2 -2((u • ∇) U , A mod 0 U ) L 2 = ([∂ t , A mod 0 ] U , U ) L 2 -2((u • ∇) U , A mod 0 U ) L 2 + 2((∂ t + u • ∇) U , Ḟ ) L 2 =: I 1 + I 2 + I 3 ,
where we used the fact that A mod 0 is a symmetric operator in L 2 and that A 1 is a skew-symmetric matrix. As for I 1 , we have

I 1 = ((∂ t a) ζ, ζ) L 2 + ℓ=1,2 ρ ℓ h ℓ n l=1 ((∂ t A ℓ )∂ l φℓ , ∂ l φℓ ) L 2 + (h ℓ δ) -2 ((∂ t C ℓ ) φℓ , φℓ ) L 2 +2ρ ℓ ([∂ t , θ ℓ l T ℓ (v • ∇)] φℓ , ζ) L 2 .
Here, as in the proof of Lemma 6.1 we have

ρ ℓ /h ℓ θ ℓ (|v| + |∂ t v|) ≲ 1 for ℓ = 1, 2. In view of the relations ∂ t θ 1 = -∂ t θ 2 = θ 1 θ 2 (H -1 1 ∂ t H 1 -H -1 2 ∂ t H 2 ), we have |∂ t θ ℓ | ≲ θ 1 θ 2 for ℓ = 1, 2. Therefore, we obtain |I 1 | ≲ E ( U ).
As for I 2 , by integration by parts we have

I 2 = ((∇ • (au)) ζ, ζ) L 2 - ℓ=1,2 ρ ℓ h ℓ n l=1 2(A ℓ ∂ l φℓ , ((∂ l u) • ∇) φℓ ) L 2 + (((u • ∇) * A ℓ )∂ l φℓ , ∂ ℓ φℓ ) L 2 + (h ℓ δ) -2 (((u • ∇) * C ℓ ) φℓ , φℓ ) L 2 + 2 ℓ=1,2 ρ ℓ ((∇ • u) ζ, θ ℓ l T ℓ (v • ∇) φℓ ) L 2 + ( ζ, [u • ∇, θ ℓ l T ℓ (v • ∇)] φℓ ) L 2 .
By using (2.14), we see that

θ 1 ≃ ρ 2 h 1 ρ 1 h 2 + ρ 2 h 1 = ρ 2 h 2 , θ 2 ≃ ρ 1 h 2 ρ 1 h 2 + ρ 2 h 1 = ρ 1 h 1 . Therefore, we have |u| ≤ θ 2 |u 1 | + θ 1 |u 2 | ≲ 1. In view of |∇θ ℓ | ≲ θ 1 θ 2 for ℓ = 1
, 2, we have also |∇u| ≲ 1 and ρ ℓ /h ℓ θ l |∇v| ≲ 1 for ℓ = 1, 2. Hence, we obtain |I 2 | ≲ E ( U ). Finally, as for I 3 , we have

I 3 = 2(∂ t ζ, ḟ0 ) L 2 -2( ζ, ∇ • (u ḟ0 )) L 2 + 2 ℓ=1,2 ρ ℓ ((∂ t + u • ∇) φℓ , ḟℓ -(∇ • (θ ℓ l ℓ ⊗ v)) ζ) L 2 ≲ ∥ ḟ0 ∥ H 1 (∥∂ t ζ∥ H -1 + ∥ ζ∥ L 2 ) + ℓ=1,2 ρ ℓ (∥ ḟℓ ∥ L 2 + ∥ ζ∥ L 2 )∥(∂ t φℓ , ∇ φℓ )∥ L 2 .
Summarizing the above estimates we obtain

d dt (A mod 0 U , U ) L 2 ≲ E ( U ) + ∥ ḟ0 ∥ H 1 (∥∂ t ζ∥ H -1 + ∥ ζ∥ L 2 ) + ℓ=1,2 ρ ℓ (∥ ḟℓ ∥ L 2 + ∥ ζ∥ L 2 )∥(∂ t φℓ , ∇ φℓ )∥ L 2 .
This together with Lemma 6.1 and Gronwall's inequality gives the desired estimate.

Energy estimates

In this subsection, we will complete the proof of Theorem 3.1. The existence and the uniqueness of the solution to the initial value problem for the Kakinuma model (6.1) has already been established in the companion paper [START_REF] Duchêne | A mathematical analysis of the Kakinuma model for interfacial gravity waves. Part I: Structures and well-posedness[END_REF], so that it is sufficient to derive the uniform bound (3.14) of the solution for some time interval [0, T ] independent of parameters. The following lemma can be shown in the same way as the proof of [20, Lemma 4.2].

Lemma 6.3. Let c, M be positive constants and m an integer such that m > n 2 + 1. There exists a positive constant C such that for any positive parameters

h 1 , h 2 , δ satisfying h 1 δ, h 2 δ ≤ 1, if ζ ∈ H m-1 , b ∈ W m,∞ , H 1 = 1 -h -1 1 ζ, and H 2 = 1 + h -1 2 ζ -h -1 2 b satisfy h -1 1 ∥ζ∥ H m-1 + h -1 2 ∥ζ∥ H m-1 + h -1 2 ∥b∥ W m,∞ ≤ M, H 1 (x) ≥ c, H 2 (x) ≥ c for x ∈ R n ,
and if φ 1 and φ 2 satisfy

L 1,i (H 1 , δ, h 1 )φ 1 = f 1,i for i = 1, 2, . . . , N, L 2,i (H 2 , b, δ, h 2 )φ 2 = f 2,i for i = 1, 2, . . . , N * ,
then for any k = 0, ±1, . . . , ±(m -1) we have

(h ℓ δ) -2 ∥φ ′ ℓ ∥ H k ≤ C(∥∇φ ℓ ∥ H k+1 + ∥φ ′ ℓ ∥ H k+1 + ∥f ′ ℓ ∥ H k ) (ℓ = 1, 2).
The next lemma gives an energy estimate of the solution to the Kakinuma model (6.1) under appropriate assumptions on the solution. We recall that the mathematical energy function E m (t) is defined by (5.3). Lemma 6.4. Let c, M, M 1 , h min be positive constants. There exist two positive constants C = C(c, M, h min ) and C 1 = C 1 (c, M, M 1 , h min ) such that for any positive parameters ρ 1 , ρ 2 , h 1 , h 2 , δ satisfying the natural restrictions (2.14), h 1 δ, h 2 δ ≤ 1, and the condition h min ≤ h 1 , h 2 , if a regular solution (ζ, ϕ 1 , ϕ 2 ) to the Kakinuma model (6.1) with a bottom topography b satisfies (6.12), h -1

2 ∥b∥ W m+1,∞ + (h 2 δ)∥b∥ W m+2,∞ ≤ M 1 , and E m (t) ≤ M 1 for some time interval [0, T ], then we have E m (t) ≤ Ce C 1 t E m (0) for 0 ≤ t ≤ T .
Proof. Let β be a multi-index such that 1 ≤ |β| ≤ m. Applying ∂ β to the Kakinuma model (6.1), after a tedious but straightforward calculation, we obtain (6.14)

     l 1 (H 1 )(∂ t + u 1 • ∇)∂ β ζ + h 1 L pr 1 (H 1 , δ, h 1 )∂ β ϕ 1 = f 1,β , l 2 (H 2 )(∂ t + u 2 • ∇)∂ β ζ -h 2 L pr 2 (H 2 , δ, h 2 )∂ β ϕ 2 = f 2,β , ρ 1 l 1 (H 1 ) • (∂ t + u 1 • ∇)∂ β ϕ 1 -ρ 2 l 2 (H 2 ) • (∂ t + u 2 • ∇)∂ β ϕ 2 -a∂ β ζ = f 0,β ,
where L pr 1 and L pr 2 are defined by (6.6), the function a by (3.9), and

f 1,β := -[∂ β , l 1 (H 1 )]∂ t ζ + h 1 [∂ β , A 1 (H 1 )]∆ϕ 1 -(l 1 (H 1 ) ⊗ l 1 (H 1 ))(∇H 1 • ∇)∂ β ϕ 1 (6.15) + [∂ β , l 1 (H 1 ) ⊗ u 1 ]∇H 1 -(h 1 δ) -2 [∂ β , C 1 (H 1 )]ϕ 1 , f 2,β := -[∂ β , l 2 (H 2 )]∂ t ζ -h 2 [∂ β , A 2 (H 2 )]∆ϕ 2 -(l 2 (H 2 ) ⊗ l 2 (H 2 ))(∇H 2 • ∇)∂ β ϕ 2 (6.16) + [∂ β , l 2 (H 2 ) ⊗ u 2 ]∇H 2 -(h 2 δ) -2 [∂ β , C 2 (H 2 )]ϕ 2 -l 2 (H 2 )(u 2 • ∂ β (h -1 2 ∇b)) -∂ β B2 (H 2 )(h -1 2 ∇b • ∇)ϕ 2 + C2 (H 2 , h -1 2 b)ϕ 2 , f 0,β := -ρ 1 [∂ β , l 1 (H 1 )] -l ′ 1 (H 1 )(∂ β H 1 ) T ∂ t ϕ 1 (6.17) + 1 2 [∂ β ; u 1 , u 1 ] + 1 2 (h 1 δ) -2 [∂ β ; w 1 , w 1 ] + u 1 • [∂ β , l 1 (H 1 )] -l ′ 1 (H 1 )(∂ β H 1 ) ⊗ ∇ T ϕ 1 -(h 1 δ) -2 w 1 [∂ β , l ′ 1 (H 1 )] -l ′′ 1 (H 1 )(∂ β H 1 ) T ϕ 1 + l ′ 1 (H 1 ) • ∂ β ϕ 1 + ρ 2 [∂ β , l 2 (H 2 )] -l ′ 2 (H 2 )(∂ β H 2 ) -l ′ 2 (H 2 )(∂ β (h -1 2 b)) T ∂ t ϕ 2 + 1 2 [∂ β ; u 2 , u 2 ] + 1 2 (h 2 δ) -2 [∂ β ; w 2 , w 2 ] + u 2 • [∂ β , l 2 (H 2 )] -l ′ 2 (H 2 )(∂ β H 2 ) -l ′ 2 (H 2 )(∂ β (h -1 2 b)) ⊗ ∇ T ϕ 2 -u 2 • [∂ β , h -1 2 ∇b ⊗ ϕ 2 ]l ′ 2 (H 2 ) -(u 2 • h -1 2 ∇b)ϕ 2 • ∂ β l ′ 2 (H 2 ) -l ′′ 2 (H 2 )(∂ β H 2 ) -l ′′ 2 (H 2 )(∂ β (h -1 2 b)) + (h 2 δ) -2 w 2 [∂ β , l ′ 2 (H 2 )] -l ′′ 2 (H 2 )(∂ β H 2 ) -l ′′ 2 (H 2 )(∂ β (h -1 2 b)) T ϕ 2 + l ′ 2 (H 2 ) • ∂ β ϕ 2 .
Here,

[∂ β ; u, v] = ∂ β (uv) -(∂ β u)v -u(∂ β v) is the symmetric commutator. For vector valued functions, it is defined by [∂ β ; u, v] = ∂ β (u • v) -(∂ β u) • v -u • (∂ β v).
On the other hand, by Lemma 5.2 we have the estimate (5.4) for time derivatives of the solution. Particularly, we have (6.18) ℓ=1,2

ρ ℓ h ℓ ∥∂ t u ℓ ∥ 2 H m-1 + (h ℓ δ) -2 ∥∂ t w ℓ ∥ 2 H m-1 + ∥∂ t ϕ ′ ℓ ∥ 2 H m + ∥∂ 2 t ϕ ′ ℓ ∥ 2 H m-1 ≲ E m .
Note that we have also the estimate (5.8) for the velocities (u ℓ , w ℓ ) (ℓ = 1, 2). Moreover, it follows from Lemma 6.3 that

ρ ℓ h ℓ (h ℓ δ) -4 ∥ϕ ′ ℓ ∥ 2 H m-1 ≲ E m for ℓ = 1, 2.
In view of the definition (3.9) of the function a, it is not difficult to check the estimate ∥a-1∥ 2

H m +∥∂ t a∥ 2 H m-1 ≲ E m . Therefore, by
the Sobolev imbedding theorem we see that all the assumptions in Proposition 6.2 are satisfied, so that for the solution U = (ζ, ϕ 1 , ϕ 2 ) T we have

E (∂ β U (t)) ≤ Ce C 1 t E (∂ β U (0)) + C 1 t 0 e C 1 (t-τ ) F β (τ )dτ,
where

F β = ∥f 0,β ∥ H 1 (∥∂ t ∂ β ζ∥ H -1 + ∥∂ β ζ∥ L 2 ) + ℓ=1,2 ρ ℓ (∥f ℓ,β ∥ L 2 + ∥∂ β ζ∥ L 2 )∥(∂ t ∂ β ϕ ℓ , ∇∂ β ϕ ℓ )∥ L 2 .
In view of the estimates (5.4), (5.8), and (6.18) together with

∥([∂ β , l ℓ (H ℓ )] -l ′ ℓ (H 1 )(∂ β H ℓ )) T φ ℓ ∥ H 1 ≲ ∥φ ′ ℓ ∥ H m for ℓ = 1, 2, we obtain F β ≲ E m .
We note that the multi-index β is assumed to satisfy 1 ≤ |β| ≤ m. As for the case

β = 0, in view of d dt E (U (t)) ≲ E m (t) we infer the inequality E (U (t)) ≤ E (U (0)) + C 1 t 0 E m (τ )dτ .
Summarizing the above estimates we obtain

E m (t) ≤ Ce C 1 t E m (0) + C 1 t 0 e C 1 (t-τ ) E m (τ )dτ with constants C = C(c, M, h min ) and C 1 = C 1 (c, M, M 1 , h min ).
Therefore, Gronwall's inequality gives the desired estimate. Now, we are ready to prove Theorem 3.1. Suppose that the initial data (ζ (0) , ϕ 1(0) , ϕ 2(0) ) and the bottom topography b satisfy (3.10)- (3.13). Let C 0 be a positive constant such that ℓ=1,2

(∥H ℓ(0) ∥ L ∞ + ρ ℓ h ℓ ∥u ℓ(0) ∥ 2 L ∞ ) + ∥a (0) ∥ L ∞ ≤ C 0 .
Such a constant C 0 exists as a constant depending on c 0 , M 0 , h min , and m. We will show that the solution (ζ, ϕ 1 , ϕ 2 ) satisfies (3.14), (3.15), and

(6.19) ℓ=1,2 (∥H ℓ (t)∥ L ∞ + ρ ℓ h ℓ ∥u ℓ (t)∥ 2 L ∞ ) + ∥a(t)∥ L ∞ ≤ 2C 0
for 0 ≤ t ≤ T with a constant M and a time T which will be determined below. We note that (3.14) is equivalent to E m (t) ≤ M . To this end, we assume that the solution satisfies (3.14), (3.15), and (6.19) for 0 ≤ t ≤ T . In the following, the constant depending on c 0 , C 0 , h min , m but not on M is denoted by C and the constant depending also on M by C 1 . These constants may change from line to line. Then, it follows from Lemma 6.4 that E m (t) ≤ Ce C 1 t M 0 for 0 ≤ t ≤ T . Therefore, if we chose M = 2CM 0 and if T is so small that T ≤ C -1 1 log 2, then (3.14) holds in fact for 0 ≤ t ≤ T . It remains to show (3.15) and (6.19). As before, we can check

     ℓ=1,2 ∥∂ t H ℓ (t)∥ L ∞ + ρ ℓ h ℓ ∥∂ t u ℓ (t)∥ L ∞ + ∥∂ t a(t)∥ L ∞ ≤ C 1 , ∥∂ t a(t) - ρ 1 ρ 2 ρ 1 H 2 (t)α 2 + ρ 2 H 1 (t)α 1 |u 1 (t) -u 2 (t)| 2 ∥ L ∞ ≤ C 1 .
Therefore, if T is so small that T ≤ (2C 1 ) -1 c 0 and T ≤ ((2C 1/2 0 + 1)C 1 ) -1 C 0 , then the lower bound (3.15) and the upper bound (6.19) hold in fact for 0 ≤ t ≤ T . This completes the proof of Theorem 3.1. respectively. Put (7.1)

           Φapp (x, z) := 2N * +2 i=0 (z + 1 -b(x)) p i φi (x), Ψ(x, z) := 2N * +2 i=0 (z + 1 -b(x)) p i ψ i (x),
and Φ res := Φ-Φapp . We note that Φapp is a higher order approximation of the velocity potential Φ and that it satisfies the boundary value problem (4.7) approximately in the sense that

     ∆ Φapp + δ -2 ∂ 2 z Φapp = R in -1 + b(x) < z < ζ(x), Φapp = ϕ on z = ζ(x), ∇b • ∇ Φapp -δ -2 ∂ z Φapp = r B on z = -1 + b(x),
where the residual R can be written in the form

R(x, z) = 2N * +2 i=0 (z + 1 -b(x)) p i r i (x).
Estimates for the residuals (r 0 , r 1 , . . . , r 2N * +2 ) and r B were given in [20, Lemmas 6.4 and 6.9].

In fact, we have ∥(r 0 , r 1 , . . . , r

2N * +2 )∥ H k + ∥r B ∥ H k ≲ δ 2j ∥∇ϕ∥ H k+2j+1 for -m ≤ k ≤ m -2(j + 1)
and 0 ≤ j ≤ 2N + 1. Now, with a slight modification from the strategy in [START_REF] Iguchi | A mathematical justification of the Isobe-Kakinuma model for water waves with and without bottom topography[END_REF], we use the identity

I 1 = Ω I δ ∇ X Φ res • I δ ∇ X ΨdX,
where we denote Ω := {X = (x, z) ; -1 + b(x) < z < ζ(x)}, I δ := diag(1, . . . , 1, δ -1 ), and

∇ X := (∇, ∂ z ) = (∂ 1 , . . . , ∂ n , ∂ z )
. Indeed, we have on one hand

(Λϕ, ψ) L 2 = Ω I δ ∇ X Φ • I δ ∇ X ΨdX
as a consequence of (4.7), Ψ(x, ζ(x)) = ψ(x), and Green's identity, and on the other hand

(Λ (2N +2) ϕ, ψ) L 2 = (L (2N +2) 0 φ, l (2N +2) • ψ) L 2 = 2N * +2 i=0 (H p i L (2N +2) 0 φ, ψ i ) L 2 = 2N * +2 i,j=0 (L ij φj , ψ i ) L 2 = Ω I δ ∇ X Φapp • I δ ∇ X ΨdX,
where the last identity follows from the expressions (4.2) and (7.1).

To evaluate I 1 , it is convenient to transform the water region Ω into a simple flat domain Ω 0 = R n × (-1, 0) by using a diffeomorphism which simply stretches the vertical direction Θ(x, z) = (x, θ(x, z)) : Ω 0 → Ω, where θ(x, z) = ζ(x)(z + 1) + (1 -b(x))z. Put Φres = Φ res • Θ and Ψ = Ψ • Θ. Then, the above integral is transformed into

I 1 = Ω 0 PI δ ∇ X Φres • I δ ∇ X ΨdX,
where

P = det ∂Θ ∂X I -1 δ ∂Θ ∂X -1 I 2 δ ∂Θ ∂X -1 T I -1 δ .
Therefore, under the restriction |k| ≤ m -1 and using the hypothesis (4.8), we have

|I 1 | ≲ ∥J k I δ ∇ X Φres ∥ L 2 (Ω 0 ) ∥J -k I δ ∇ X Ψ∥ L 2 (Ω 0 ) , where J = (1 -∆) 1 2
. Moreover, Φres satisfies the boundary value problem

     ∇ X • I δ PI δ ∇ X Φres = -R in Ω 0 , Φres = 0 on z = 0, e z • I δ PI δ ∇ X Φres = -r B on z = -1, where R = R • Θ = 2N * +2 i=0
(z + 1) p i H p i r j and e z = (0, . . . , 0, 1) T . By applying the standard theory of elliptic partial differential equations to the above problem, for 0 ≤ k ≤ m -1 we have

∥J k I δ ∇ X Φres ∥ L 2 (Ω 0 ) ≲ δ(∥J k R∥ L 2 (Ω 0 ) + ∥r B ∥ H k ) ≲ δ(∥(r 0 , r 1 , . . . , r 2N * +2 )∥ H k + ∥r B ∥ H k ).
Moreover, in view of Ψ = 2N * +2 i=0 (z + 1) p i H p i ψ j and by Lemma 4.1, we have

∥J -k I δ ∇ X Ψ∥ L 2 (Ω 0 ) ≲ ∥∇ψ∥ H -k + δ -1 ∥ψ ′ ∥ H -k ≲ ∥∇ψ∥ H -k
for |k| ≤ m -1. Summarizing the above estimates we have |I 1 | ≲ δ 2j+1 ∥∇ϕ∥ H k+2j+1 ∥∇ψ∥ H -k for 0 ≤ k ≤ m -2(j + 1) and 0 ≤ j ≤ 2N + 1.

As for the term I 2 , the evaluation is exactly the same as in [START_REF] Iguchi | A mathematical justification of the Isobe-Kakinuma model for water waves with and without bottom topography[END_REF]. In fact, the identities

I 2 = 2N * +2 i,j=0 (L ij φj , ψ i ) L 2 - N * j=0 (L 0j ϕ j , ψ) L 2 = N * j=0 2N * +2 i=N * +1 ((L ij -H p i L 0j )φ j , ψ i ) L 2 - 2N * +2 i,j=N * +1 ((L ij -H p i L 0j ) φj , ψ i ) L 2
were shown in [START_REF] Iguchi | A mathematical justification of the Isobe-Kakinuma model for water waves with and without bottom topography[END_REF]Equation (7.7)], where φ := (φ 0 , φ 1 , . . . , φ N * ) was defined by φ i := ϕ i -φi for i = 0, 1, . . . , N * . Now, we decompose j = j 1 + j 2 such that 1 ≤ j 1 ≤ N + 1 and 0 ≤ j 2 ≤ N . Then, by [20, Lemmas 5.2, 5.4, 6.2 and 6.7] we see that

|I 2 | ≲ {∥φ∥ H k+2j 1 +1 + ∥( φN * +1 , . . . , φ2N * +2 )∥ H k+2j 1 +1 + δ -2 (∥φ∥ H k+2j 1 -1 + ∥( φN * +1 , . . . , φ2N * +2 )∥ H k+2j 1 -1 )}∥(ψ N * +1 , . . . , ψ 2N * +2 )∥ H -(k+2j 1 -1) ≲ δ 2(j 1 +j 2 ) ∥∇ϕ∥ H k+2(j 1 +j 2 ) ∥∇ψ∥ H -k if max{|k|, |k+2j 1 -2|, |k+2j 1 +1|, |k+2(j 1 +j 2 )|} ≤ m-1 and max{|k|, |k+1|, |k+2j 1 -1|} ≤ m.
These conditions are satisfied under the restriction -m + 1 ≤ k ≤ m -2(j + 1).

To summarize, we obtain as desired |((Λ -Λ (N ) )ϕ, ψ) L 2 | ≲ δ 2j ∥∇ϕ∥ H k+2j+1 ∥∇ψ∥ H -k for 0 ≤ k ≤ m -2(j + 1) and 1 ≤ j ≤ 2N + 1. The proof is complete. This lemma and the scaling relations (4.15) imply immediately the following lemma. 0 ≤ k ≤ m -4(N + 1) satisfying the compatibility condition Λ 1 (ζ, δ, h 1 )ϕ 1 + Λ 2 (ζ, b, δ, h 2 )ϕ 2 = 0 the solution (ϕ 1 , ϕ 2 ) to (7.3) and the solution ( φ1 , φ2 ) to (7.2) satisfy ℓ=1,2

ρ ℓ h ℓ (∥∇( φℓ -ϕ ℓ )∥ 2 H k + (h ℓ δ) -2 ∥ φ′ ℓ -ϕ ′ ℓ ∥ 2 H k + (h ℓ δ) -4 ∥ φ′ ℓ -ϕ ′ ℓ ∥ 2 H k-1 ) ≤ C ℓ=1,2 ρ ℓ h ℓ (h ℓ δ) 2(4N +2) ∥∇ϕ ℓ ∥ 2 H k+4N +3 .
Proof. For simplicity, we write L 1,i = L 1,i (H 1 , δ, h 1 ), l 1 = l 1 (H 1 ), and so on. We recall that Λ

(N )

1 : ϕ 1 → L 1,0 ϕ 1 and Λ (N ) 2 : ϕ 2 → L 2,0 ϕ 2 . Notice that φℓ -ϕ ℓ for ℓ = 1, 2 satisfy            L 1,i ( φ1 -ϕ 1 ) = 0 for i = 1, 2, . . . , N, L 2,i ( φ2 -ϕ 2 ) = 0 for i = 1, 2, . . . , N * , h 1 L 1,0 ( φ1 -ϕ 1 ) + h 2 L 2,0 ( φ2 -ϕ 2 ) = (Λ 1 -h 1 Λ (N ) 1 )ϕ 1 + (Λ 2 -h 2 Λ (N ) 2 )ϕ 2 , ρ 2 l 2 • ( φ2 -ϕ 2 ) -ρ 1 l 1 • ( φ1 -ϕ 1 ) = 0.
Since the right-hand side of the third equation can be written as ∇ • f 3 with

f 3 = -∇(-∆) -1 (Λ 1 -h 1 Λ (N ) 1 )ϕ 1 -(Λ 2 -h 2 Λ (N )
2 )ϕ 2 , by Lemmas 5.1 and 7.2 we obtain ℓ=1,2

ρ ℓ h ℓ (∥∇( φℓ -ϕ ℓ )∥ 2 H k + (h ℓ δ) -2 ∥ φ′ ℓ -ϕ ′ ℓ ∥ 2 H k ) ≲ min ρ 1 h 1 , ρ 2 h 2 ∥f 3 ∥ 2 H k ≲ l=1,2 ρ ℓ h ℓ ∥(-∆) -1 2 (Λ ℓ -h ℓ Λ (N ) ℓ )ϕ ℓ ∥ 2 H k ≲ l=1,2 ρ ℓ h ℓ (h ℓ δ) 2(4N +2) ∥∇ϕ ℓ ∥ 2 H k+4N +3 .
Moreover, it follows from Lemma 6.3 that

(h ℓ δ) -2 ∥ φ′ ℓ -ϕ ′ ℓ ∥ H k-1 ≲ ∥∇( φℓ -ϕ ℓ )∥ H k + (h ℓ δ) -1 ∥ φ′ ℓ -ϕ ′ ℓ ∥ H k for ℓ = 1, 2.
This completes the proof.

The following proposition gives another version of Theorem 3.6 for the consistency of the Kakinuma model. Proposition 7.6. Let c, M be positive constants and m an integer such that m ≥ 4N + 4 and m > n 2 + 2. We assume (H1) or (H2). There exists a positive constant C such that for any positive parameters ρ 1 , ρ 2 , h 1 , h 2 , δ satisfying h 1 δ, h 2 δ ≤ 1, and for any solution (ζ, ϕ 1 , ϕ 2 ) to the full model for interfacial gravity waves (2.17) on a time interval [0, T ] satisfying (3.18), if we define H 1 and H 2 as in (2. [START_REF] Iguchi | Isobe-Kakinuma model for water waves as a higher order shallow water approximation[END_REF]) and ( φ1 , φ2 ) as a solution to (7.2), then (ζ, φ1 , φ2 ) satisfy approximately the Kakinuma model as

(7.4)            l 1 (H 1 )h -1 1 ∂ t ζ + L 1 (H 1 , δ, h 1 ) φ1 = r 1 , l 2 (H 2 )h -1 2 ∂ t ζ -L 2 (H 2 , b, δ, h 2 ) φ2 = r 2 , ρ 1 l 1 (H 1 ) • ∂ t φ1 + 1 2 | ũ1 | 2 + (h 1 δ) -2 w2 1 -ρ 2 l 2 (H 2 ) • ∂ t φ2 + 1 2 | ũ2 | 2 + (h 2 δ) -2 w2
where ũ1 , ũ2 , w1 , w2 are defined by (3.4) with (ϕ 1 , ϕ 2 ) replaced by ( φ1 , φ2 ), and the errors (r 1 , r 2 , r 0 ) satisfy

(7.5)        ℓ=1,2 ρ ℓ h ℓ ∥r ℓ (t)∥ 2 H m-(4N +5) ≤ C ℓ=1,2 ρ ℓ h ℓ (h ℓ δ) 2(4N +2) ∥∇ϕ ℓ (t)∥ 2 H m-1 , ∥r 0 (t)∥ H m-4(N +1) ≤ C (h 1 δ) 4N +2 + (h 2 δ) 4N +2 (h -1 1 + h -1 2 ) ℓ=1,2 ρ ℓ h ℓ ∥∇ϕ ℓ (t)∥ 2 H m-1 , for t ∈ [0, T ].
Proof. Let ϕ 1 and ϕ 2 be the unique solutions to (7.3), and (r 1 , r2 , r0 ) the errors in Theorem 3.6.

Then, the errors (r 1 , r 2 , r 0 ) in the proposition can be written as

           r 1 = r1 -L 1 (H 1 , δ, h 1 )( φ1 -ϕ 1 ), r 2 = r2 + L 2 (H 2 , b, δ, h 2 )( φ2 -ϕ 2 ), r 0 = r0 + ρ 1 {h -1 1 (∂ t ζ)( w1 -w 1 ) -1 2 ( ũ1 + u 1 ) • ( ũ1 -u 1 ) + (h 1 δ) -2 ( w1 + w 1 )( w1 -w 1 ) } -ρ 2 {h -1 2 (∂ t ζ)( w2 -w 2 ) -1 2 ( ũ2 + u 2 ) • ( ũ2 -u 2 ) + (h 2 δ) -2 ( w2 + w 2 )( w2 -w 2 ) }. Therefore, we have ∥r ℓ -rℓ ∥ H k ≲ φℓ -ϕ ℓ )∥ H k+1 + ∥ φ′ ℓ -ϕ ′ ℓ ∥ H k+1 + (h ℓ δ) -2 ∥ φ′ ℓ -ϕ ′ ℓ ∥ H k for -m ≤ k ≤ m -1 and ℓ = 1, 2.
Applying this estimate with k = m -(4N + 5) and the estimate in Lemma 7.5 with k = m -4(N + 1) and using the result in Theorem 3.6, we obtain the first estimate in (7.5). Since m -2 > n 2 , we have

∥r 0 -r0 ∥ H k ≲ ℓ=1,2 ρ ℓ {(∥ ũℓ ∥ H m-2 + ∥u ℓ ∥ H m-2 )∥ ũℓ -u ℓ ∥ H k + h -1 ℓ ∥∂ t ζ∥ H m-2 + (h ℓ δ) -2 (∥ wℓ ∥ H m-2 + ∥w ℓ ∥ H m-2 ) ∥ wℓ -w ℓ ∥ H k } for |k| ≤ m -2.
Here, it follows from Lemmas 4.4, 5.1, and 7.5 that ℓ=1,2

ρ ℓ h ℓ (∥u ℓ ∥ 2 H m-1 + (h ℓ δ) -2 ∥w ℓ ∥ 2 H m-1 ) ≲ ℓ=1,2 ρ ℓ h ℓ (∥∇ϕ ℓ ∥ 2 H m-1 + (h ℓ δ) -2 ∥ϕ ′ ℓ ∥ 2 H m-1 ) ≲ ℓ=1,2 ρ ℓ h ℓ ∥∇ϕ ℓ ∥ 2 H m-1 , ℓ=1,2 ρ ℓ h ℓ (∥ ũℓ ∥ 2 H m-1 + (h ℓ δ) -2 ∥ wℓ ∥ 2 H m-1 ) ≲ ℓ=1,2 ρ ℓ h ℓ (∥∇ φℓ ∥ 2 H m-1 + (h ℓ δ) -2 ∥ φ′ ℓ ∥ 2 H m-1 ) ≲ min h 1 ρ 1 , h 2 ρ 2 ∥∇(ρ 2 ϕ 2 -ρ 1 ϕ 1 )∥ 2 H m-1 ≲ ℓ=1,2 ρ ℓ h ℓ ∥∇ϕ ℓ ∥ 2 H m-1 , and
ℓ=1,2 ρ ℓ h ℓ (∥ ũℓ -u ℓ ∥ 2 H k + (h ℓ δ) -2 ∥ wℓ -w ℓ ∥ 2 H k ) ≲ ℓ=1,2 ρ ℓ h ℓ (∥∇( φℓ -ϕ ℓ )∥ 2 H k + (h ℓ δ) -2 ∥ φ′ ℓ -ϕ ′ ℓ ∥ 2 H k ) ≲ ℓ=1,2 ρ ℓ h ℓ (h ℓ δ) 2(4N +2) ∥∇ϕ ℓ ∥ 2 H k+4N +3
for 0 ≤ k ≤ m-4(N +1). Moreover, it follows from Lemma 7.4 that

∥∂ t ζ∥ H m-2 = ∥Λ ℓ ϕ ℓ ∥ H m-2 ≲ h ℓ ∥∇ϕ ℓ ∥ H m-1 for ℓ = 1, 2.
Summarizing the above estimates and using the result in Theorem 3.6, we easily obtain the second estimate in (7.5). The proof is complete.

7.3 Completion of the proof of Theorem 3.9

Now we are ready to prove Theorem 3.9. Let (ζ IW , ϕ IW 1 , ϕ IW 2 ) be the solution to the full model for interfacial gravity waves (2.17) with uniform bound stated in the theorem, and define ϕ IW := ρ 2 ϕ IW 2 -ρ 1 ϕ IW 1 , which is a canonical variable of the full model. We first ensure a uniform bound on the time derivative of the canonical variables (ζ IW , ϕ IW ). It follows from the first and the second equations in (2.17) that

∂ t ζ IW = -Λ IW 1 ϕ IW 1 = Λ IW 2 ϕ IW 2 , where Λ IW 1 = Λ 1 (ζ IW , δ, h 1 ) and Λ IW 2 = Λ 2 (ζ IW , b, δ, h 2 )
. Similar notations will be used in the following without any comment. Therefore, by Lemma 7.4 we have

∥∂ t ζ IW ∥ 2 H m-1 = min{∥Λ IW 1 ϕ IW 1 ∥ 2 H m-1 , ∥Λ IW 2 ϕ IW 2 ∥ 2 H m-1 } ≲ min{h 2 1 ∥∇ϕ IW 1 ∥ 2 H m , h 2 2 ∥∇ϕ IW 2 ∥ 2 H m } ≲ min h 1 ρ 1 , h 2 ρ 2 ℓ=1,2 ρ ℓ h ℓ ∥∇ϕ IW ℓ ∥ 2 H m ≤ 2 ℓ=1,2 ρ ℓ h ℓ ∥∇ϕ IW ℓ ∥ 2 H m ,
where we used (2.15). It follows from the third equation in (2.17) that

∂ t ϕ IW = ρ 2 ∂ t ϕ IW 2 -ρ 1 ∂ t ϕ IW 1 = 1 2 ρ 1 |∇ϕ IW 1 | 2 -δ 2 (Λ IW 1 ϕ IW 1 -∇ζ IW • ∇ϕ IW 1 ) 2 1 + δ 2 |∇ζ IW | 2 - 1 2 ρ 2 |∇ϕ IW 2 | 2 -δ 2 (Λ IW 2 ϕ IW 2 + ∇ζ IW • ∇ϕ IW 2 ) 2 1 + δ 2 |∇ζ IW | 2 -ζ IW .
Here, we note that in view of the conditions h 1 δ, h 2 δ ≤ 1 and h -1 1 , h -1 2 ≲ 1 we have δ ≲ 1. Therefore, by Lemma 7.4 we have

∥∂ t ϕ IW ∥ H m-1 ≲ ∥ζ IW ∥ H m-1 + ℓ=1,2 ρ ℓ {∥∇ϕ IW ℓ ∥ 2 H m-1 + δ 2 (h 2 ℓ ∥∇ϕ IW ℓ ∥ 2 H m + ∥∇ϕ IW ℓ ∥ 2 H m-1 )} ≲ ∥ζ IW ∥ H m-1 + ℓ=1,2 ρ ℓ h ℓ ∥∇ϕ IW ℓ ∥ 2 H m .
Hence, we obtain

∥∂ t ζ IW ∥ H m-1 + ∥∂ t ϕ IW ∥ H m-1 ≲ 1. Let (
ρ ℓ h ℓ (∥∇ φIW ℓ ∥ 2 H m + (h ℓ δ) -2 ∥ φIW ′ ℓ ∥ 2 H m ) ≲ min h 1 ρ 1 , h 2 ρ 2 ∥∇ϕ IW ∥ 2 H m ≲ ℓ=1,2 ρ ℓ h ℓ ∥∇ϕ IW ℓ ∥ 2 H m ≲ 1, which yields ℓ=1,2 ρ ℓ h ℓ (∥ ũIW ℓ ∥ 2 H m + (h ℓ δ) -2 ∥ wIW ℓ ∥ 2 H m + (h ℓ δ) -4 ∥ φIW ′ ℓ ∥ 2 H m-1 ) ≲ 1,
where ũIW ). To this end, we derive equations for these time derivatives by differentiating (7.2) with respect to t. The procedure is almost the same as in the proof of Lemma 5.2. The only difference is the last equation in (5.5), especially, the expression of f 4 . In this case, f 4 has the form

f 4 = ∂ t ϕ IW + ρ 1 wIW 1 h -1 1 ∂ t ζ IW -ρ 2 wIW 2 h -1 2 ∂ t ζ IW , so that ∥f 4 ∥ H m-1 ≲ 1. Therefore, we obtain ℓ=1,2 ρ ℓ h ℓ (∥∇∂ t φIW ℓ ∥ 2 H m-2 + (h ℓ δ) -2 ∥∂ t φIW ′ ℓ ∥ 2 H m-2 ) ≲ 1. Let (ζ K , ϕ K 1 , ϕ K
2 ) be the solution to the initial value problem for the Kakinuma model (2.18)-(2.20) stated in the theorem, whose unique existence is guaranteed by Theorem 3.1 and Proposition 3.4. Note also that the solution satisfies the uniform bound (3.14) together with the stability and non-cavitation conditions (3.15). It follows from Lemma 6.3 that

ρ ℓ h ℓ (h ℓ δ) -4 ∥ϕ K ′ ℓ ∥ 2 H m-1 ≲ 1 for ℓ = 1, 2. Moreover, the time derivatives (∂ t ζ K , ∂ t ϕ K 1 , ∂ t ϕ K
2 ) satisfy (5.4) and (u K ℓ , w K ℓ ) (ℓ = 1, 2), which are defined by (3.4) with (ϕ 1 , ϕ 2 ) replaced by (ϕ K 1 , ϕ K 2 ), satisfy (5.8). Putting

ζ res := ζ K -ζ IW , ϕ res ℓ := ϕ K ℓ -φIW ℓ (ℓ = 1, 2),
we will show that (ζ res , ϕ res 1 , ϕ res 2 ) can be estimated by the errors (r 1 , r 2 , r 0 ). To this end, we are going to evaluate

E res k (t) := ∥ζ res (t)∥ 2 H k + ℓ=1,2 ρ ℓ h ℓ (∥∇ϕ res ℓ (t)∥ 2 H k + (h ℓ δ) -2 ∥ϕ res ′ ℓ (t)∥ 2 H k )
for an appropriate integer k by making use of energy estimates similar to the ones obtained in Sections 5 and 6 for the proof of the well-posedness of the initial value problem for the Kakinuma model (2.18)-(2.20). Here, we note that E res k (0) = 0. As in the case of the energy estimate for the Kakinuma model, we first need to evaluate times derivatives (∂ t ζ res , ∂ t ϕ res 1 , ∂ t ϕ res 2 ) in terms of E res k . By taking difference between the first components of the first two equations in (3.5) and (7.4), ∂ t ζ res can be written in two ways as

∂ t ζ res = -h 1 {L K 1,0 ϕ res 1 + (L K 1,0 -L IW 1,0 ) φIW 1 + r 1,0 } = h 2 {L K 2,0 ϕ res 2 + (L K 2,0 -L IW 2,0 ) φIW 2 + r 2,0 }, where L K 1,0 = L 1,0 (H K 1 , δ, h 1 ), H K 1 = 1 -h -1 1 ζ K ,
and similar simplifications are used, and r ℓ,0 is the 0th component of the error r ℓ for ℓ = 1, 2. Therefore, we have

∥∂ t ζ res ∥ H k-1 ≲ h ℓ {∥∇ϕ res ℓ ∥ H k + ∥ϕ res ′ ℓ ∥ H k + ∥ζ res ∥ H k (∥∇ φIW ℓ ∥ H m + ∥ φIW ′ ℓ ∥ H m ) + ∥r ℓ,0 ∥ H k-1 }
Here, f K 1,β , f K 2,β , and f K 0,β are those in (6.15)-(6.17

) with (ζ, ϕ 1 , ϕ 2 ) = (ζ K , ϕ K 1 , ϕ K 2 )
, and so on. As we saw, all the assumptions in Proposition 6.2 are satisfied, so that we have

E (∂ β U res (t)) ≲ t 0 F res β (τ )dτ,
where U res := (ζ res , ϕ res 1 , ϕ res 2 ) T , E is defined in (6.11), and

F res β = ∥f res 0,β ∥ H 1 (∥∂ t ζ res ∥ H k-1 + ∥ζ res ∥ H k ) + ℓ=1,2 ρ ℓ (∥f res ℓ,β ∥ L 2 + ∥ζ res ∥ H k )(∥∇∂ t ϕ res ℓ ∥ H k-1 + ∥∇ϕ res ℓ ∥ H k ). In view of ∥(ζ IW , ζ K )∥ H m ≲ 1, straightforward calculations yield ∥f res ℓ,β ∥ L 2 ≲ (∥∂ t ζ IW ∥ H m-1 + ∥ ũIW ℓ ∥ H m )∥ζ res ∥ H k + h ℓ (∥∇ φIW ℓ ∥ H m + ∥ φIW ′ ℓ ∥ H m + (h ℓ δ) -2 ∥ φIW ′ ℓ ∥ H m-1 )∥ζ res ∥ H k + h ℓ (∥∇ φres ℓ ∥ H k + ∥ φres ′ ℓ ∥ H k + (h ℓ δ) -2 ∥ φres ′ ℓ ∥ H k-1 ) + ∥∂ t ζ res ∥ H k-1 + ∥ ũIW ℓ -u K ℓ ∥ H k + h ℓ ∥r ℓ ∥ H k for ℓ = 1, 2 and n 2 < k ≤ m -1.
As for f res 0,β , we note the relation

[∂ β , l K 2 ] -l ′ 2 (H K 2 )(h -1 2 ∂ β ζ K ) -[∂ β , l IW 2 ] -l ′ 2 (H IW 2 )(h -1 2 ∂ β ζ IW ) T ∂ t φIW 2 = 1 0 [∂ β , l ′ 2 (sH IW + (1 -s)H K )] -l ′′ 2 (sH IW + (1 -s)H K )h -1 2 ∂ β (sζ IW + (1 -s)ζ K ) T (h -1 2 ζ res )∂ t φIW 2 + l ′ 2 (sH IW + (1 -s)H K ) [∂ β , h -1 2 ζ res ] -(∂ β (h -1 2 ζ res )) T ∂ t φIW 2 ds.
Therefore, straightforward calculations yield

∥f res 0,β ∥ H 1 ≲ l=1,2 ρ ℓ (∥∇∂ t φIW ℓ ∥ H m-2 + ∥∂ t φIW ′ ℓ ∥ H m-2 )∥ζ res ∥ H k + (∥ ũIW ℓ ∥ H m + ∥u K ℓ ∥ H m )(∥∇ φIW ℓ ∥ H m + ∥ φIW ′ ℓ ∥ H m )∥ζ res ∥ H k + (h ℓ δ) -2 ∥w K ℓ ∥ H m ∥ φIW ′ ℓ ∥ H m ∥ζ res ∥ H k + ∥∇∂ t φres ℓ ∥ H k-1 + ∥∂ t φres ′ ℓ ∥ H k-1 + ∥u K ℓ ∥ H m (∥∇ φres ℓ ∥ H k + ∥ φres ′ ℓ ∥ H k ) + (h ℓ δ) -2 ∥w K ℓ ∥ H m ∥ φres ′ ℓ ∥ H k + (∥u K ℓ ∥ H m + ∥ ũIW ℓ ∥ H m + ∥∇ φIW ℓ ∥ H m + ∥ φIW ′ ℓ ∥ H m )∥ ũIW ℓ -u K ℓ ∥ H k + (h ℓ δ) -2 (∥w K ℓ ∥ H m + ∥ wIW ℓ ∥ H m + ∥ φIW ′ ℓ ∥ H m )∥ wIW ℓ -w K ℓ ∥ H k + ∥r 0 ∥ H k+1 for n 2 < k ≤ m -2.
In view of the above estimates and (7.6)-(7.8) we obtain

F res β ≲ E res k + R k with R k := ∥r 0 ∥ 2 H k+1 + l=1,2 ρ ℓ h ℓ ∥r ℓ ∥ 2 H k .
We note that the multi-index β is assumed to satisfy 1 ≤ |β| ≤ k. As for the case β = 0, we have

(t) ≲ (h 1 δ) 4N +2 + (h 2 δ) 4N +2 for 0 ≤ t ≤ min{T, T IW }. It remains to evaluate ϕ IW ℓ -ϕ K ℓ for ℓ = 1, 2. Let (ϕ IW 1 , ϕ IW 2 ) be the solution to (3.19) with (ζ, ϕ 1 , ϕ 2 ) = (ζ IW , ϕ IW 1 , ϕ IW 2 ). Then, we have ϕ K ℓ -ϕ IW ℓ = l K ℓ •ϕ res ℓ +(l K ℓ -l IW ℓ )• φIW ℓ +l IW ℓ •( φIW ℓ -ϕ IW ℓ ), so that for any 0 ≤ k ≤ m -1 ∥∇ϕ K ℓ -∇ϕ IW ℓ ∥ H k ≲ ∥∇ϕ res ℓ ∥ H k + ∥ϕ res ′ ℓ ∥ H k + h -1 ℓ ∥ζ res ∥ H k+1 ∥ φIW ′ ℓ ∥ H m + ∥∇( φIW ℓ -ϕ IW ℓ )∥ H k + ∥ φIW ′ ℓ -ϕ IW ′ ℓ ∥ H k .
Therefore, the previous result together with Lemma 7.5 implies ℓ=1,2

ρ ℓ h ℓ ∥∇ϕ K ℓ -∇ϕ IW ℓ ∥ 2 H m-(4N +5) ≲ (h 1 δ) 4N +2 + (h 2 δ) 4N +2 .
This completes the proof of Theorem 3.9.

((h 1 δ) 4N +2 + (h 2 δ) 4N +2 ).

Preliminary elliptic estimates

We consider the following transmission problem (8.1)

           ∇ X • I 2 δ ∇ X Φ ℓ = 0 in Ω ℓ (ℓ = 1, 2), n • I 2 δ ∇ X Φ ℓ = 0 on Σ ℓ (ℓ = 1, 2), n • I 2 δ ∇ X Φ 2 -n • I 2 δ ∇ X Φ 1 = r S on Γ, ρ 2 Φ 2 -ρ 1 Φ 1 = ϕ on Γ,
where the rigid-lid Σ 1 of the upper layer Ω 1 , the bottom Σ 2 of the lower layer Ω 2 , and the interface Γ are defined by z = h 1 , z = -h 2 + b(x), and z = ζ(x), respectively, I δ := diag(1, . . . , 1, δ -1 ), ∇ X := (∇, ∂ z ) T = (∂ 1 , . . . , ∂ n , ∂ z ), and n is an upward normal vector, specifically, n = e z on Σ 1 , n = (-∇b, 1) T on Σ 2 , and n = (-ζ, 1) T on Γ.

Lemma 8.1. Let c, M be positive constants. There exists a positive constant C such that for any positive parameters

ρ 1 , ρ 2 , h 1 , h 2 , δ satisfying h 1 δ, h 2 δ ≤ 1, if ζ, b ∈ W 1,∞ , H 1 = 1 -h -1 1 ζ, and H 2 = 1 + h -1 2 ζ -h -1 2 b satisfy h -1 1 ∥ζ∥ W 1,∞ + h -1 2 ∥ζ∥ W 1,∞ + h -1 2 ∥b∥ W 1,∞ ≤ M, H 1 (x) ≥ c, H 2 (x) ≥ c for x ∈ R n , then for any (r S , ϕ) satisfying ∇ϕ ∈ H -1 2 and (-∆) -1 2 r S ∈ H 1 2
there exists a solution (Φ 1 , Φ 2 ) to the transmission problem (8.1). The solution is unique up to an additive constant of the form (ρ 2 C, ρ 1 C) and satisfies

ℓ=1,2 ρ ℓ ∥I δ ∇ X Φ ℓ ∥ 2 L 2 (Ω ℓ ) (8.2) ≤ C ∥((ρ 1 Λ 2,0 + ρ 2 Λ 1,0 ) -1 Λ 1,0 Λ 2,0 ) 1 2 ϕ∥ 2 L 2 + ρ 1 ρ 2 ∥(ρ 1 Λ 2,0 + ρ 2 Λ 1,0 ) -1 2 r S ∥ 2 L 2 ,
where Λ 1,0 = Λ 1 (0, δ, h 1 ) and Λ 2,0 = Λ 2 (0, 0, δ, h 2 ) are Dirichlet-to-Neumann maps in the case ζ(x) ≡ b(x) ≡ 0. Particularly, if we further impose ϕ ∈ H1 , (-∆) -1 2 r S ∈ H 1 , the natural restrictions (2.14), and h min ≤ h 1 , h 2 with a positive constant h min , then we have

(8.3) ℓ=1,2 ρ ℓ ∥I δ ∇ X Φ ℓ ∥ 2 L 2 (Ω ℓ ) ≤ C∥∇ϕ∥ 2 L 2 + C min ℓ=1,2 ρ ℓ h ℓ ∥((-∆) -1 2 + h ℓ δ)r S ∥ 2 L 2 ,
where the constant C depends also on h min .

Proof. The existence and the uniqueness of the solution is standard, so that we focus on deriving the uniform estimate of the solution. To this end, it is convenient to transform the water regions Ω 1 and Ω 2 into simple domains Ω

1,0 = R n × (0, h 1 ) and Ω 2,0 = R n × (-h 2 , 0) by using diffeomorphisms Θ ℓ (x, z) = (x, θ ℓ (x, z)) : Ω ℓ,0 → Ω ℓ (ℓ = 1, 2), respectively, where θ 1 (x, z) = (1 -h -1 1 ζ(x))z + ζ(x) and θ 2 (x, z) = (1 + h -1 2 (ζ(x) -b(x)))z + ζ(x). Put Φℓ = Φ ℓ • Θ ℓ (ℓ = 1, 2). Then, the transmission problem (8.1) is transformed into            ∇ X • I δ P ℓ I δ ∇ X Φℓ = 0 in Ω ℓ,0 (ℓ = 1, 2), e z • I δ P ℓ I δ ∇ X Φℓ = 0 on Σ ℓ,0 (ℓ = 1, 2), e z • I δ P 2 I δ ∇ X Φ2 -e z • I δ P 1 I δ ∇ X Φ1 = r S on Γ 0 , ρ 2 Φ2 -ρ 1 Φ1 = ϕ on Γ 0 ,
where Σ 1,0 , Σ 2,0 , and Γ 0 are represented as z = h 1 , z = -h 2 , and z = 0, respectively, and

P ℓ := det ∂Θ ℓ ∂X I -1 δ ∂Θ ℓ ∂X -1 I 2 δ ∂Θ ℓ ∂X -1 T I -1 δ (ℓ = 1, 2).
We note that ∥I δ ∇ X Φ ℓ ∥ L 2 (Ω ℓ ) ≃ ∥I δ ∇ X Φℓ ∥ L 2 (Ω ℓ,0 ) (ℓ = 1, 2). Let (Ψ 1 , Ψ 2 ) be a solution to the transmission problem

           ∇ X • I 2 δ ∇ X Ψ ℓ = 0 in Ω ℓ,0 (ℓ = 1, 2), e z • I 2 δ ∇ X Ψ ℓ = 0 on Σ ℓ,0 (ℓ = 1, 2), e z • I 2 δ ∇ X Ψ 2 -e z • I 2 δ ∇ X Ψ 1 = r S on Γ 0 , ρ 2 Ψ 2 -ρ 1 Ψ 1 = ϕ on Γ 0 ,
and put Φ res ℓ = Φℓ -Ψ ℓ (ℓ = 1, 2). Then, we can decompose

|I δ ∇ X Φ res ℓ | 2 -I δ ∇ X Φ res ℓ • (I -P ℓ )I δ ∇ X Φℓ = ∇ X Φ res ℓ • {(I δ P ℓ I δ ∇ X Φℓ -I 2 δ ∇ X Ψ ℓ )}
for ℓ = 1, 2 and ρ 1 Φ res 1 = ρ 2 Φ res 2 on z = 0. Therefore, denoting the unit outward normal vector to ∂Ω ℓ,0 by N ℓ (ℓ = 1, 2) we have Similarly, in view of the decomposition

I δ ∇ X Φ res ℓ • P ℓ I δ ∇ X Φ res ℓ -I δ ∇ X Φ res ℓ • (I -P ℓ )I δ ∇ X Ψ ℓ = ∇ X Φ res ℓ • {(I δ P ℓ I δ ∇ X Φℓ -I 2 δ ∇ X Ψ ℓ )} for ℓ = 1,
ρ ℓ ∥I δ ∇ X Φℓ ∥ 2 L 2 (Ω ℓ,0 ) ≃ ℓ=1,2 ρ ℓ ∥I δ ∇ X Ψ ℓ ∥ 2 L 2 (Ω ℓ,0 ) .
Therefore, it is sufficient to evaluate the right-hand side of the above equation. In other words, the evaluation is reduced to the simple case ζ(x) ≡ b(x) ≡ 0. Putting ψ ℓ = Ψ ℓ | z=0 (ℓ = 1, 2), we see that

ℓ=1,2 ρ ℓ ∥I δ ∇ X Ψ ℓ ∥ 2 L 2 (Ω ℓ,0 ) = ℓ=1,2
ρ ℓ (Λ ℓ,0 ψ ℓ , ψ ℓ ) L 2 and that Λ 1,0 ψ 1 + Λ 2,0 ψ 2 = r S , ρ 2 ψ 2 -ρ 1 ψ 1 = ϕ.

Particularly, we have

ψ 1 ψ 2 = (ρ 1 Λ 2,0 + ρ 2 Λ 1,0 ) -1 -Λ 2,0 ϕ + ρ 2 r S Λ 1,0 ϕ + ρ 1 r S .
Therefore,

ℓ=1,2 ρ ℓ ∥I δ ∇ X Ψ ℓ ∥ 2 L 2 (Ω ℓ,0 ) = ∥((ρ 1 Λ 2,0 + ρ 2 Λ 1,0 ) -1 Λ 1,0 Λ 2,0 ) 1 2 ϕ∥ 2 L 2 if r S = 0, ρ 1 ρ 2 ∥(ρ 1 Λ 2,0 + ρ 2 Λ 1,0 ) -1 2 r S ∥ 2 L 2 if ϕ = 0.
Hence, by the linearity of the problem we obtain (8.2). Finally, in order to show (8.3) it is sufficient to evaluate the symbols of the Fourier multipliers (ρ 1 Λ 2,0 + ρ 2 Λ 1,0 ) -1 Λ 1,0 Λ 2,0 and ρ 1 ρ 2 (ρ 1 Λ 2,0 + ρ 2 Λ 1,0 ) -1 . We recall that the symbol of the Dirichlet-to-Neumann map Λ ℓ,0 is given by σ(Λ ℓ,0 ) = δ -1 |ξ| tanh(h ℓ δ|ξ|) for ℓ = 1, 2. In view of 0 ≤ tanh ξ ≤ ξ for ξ ≥ 0, we have σ((ρ 1 Λ 2,0 + ρ 2 Λ 1,0 ) -1 Λ 1,0 Λ 2,0 ) ≤ min σ(Λ where we used 1 ≲ h 1 , h 2 . These estimates imply (8.3). The proof is complete.

Completion of the proof of Theorem 3.10

Now we are ready to prove Theorem 3.10. We recall the definitions (3.3) of l 1 (H 1 ), l 2 (H 2 ) and (3.6) of the operators L 1,i (H 1 , δ, h 1 ) and L 2,i (H 2 , b, δ, h 2 ). These depend on N , so that we denote them by l We proceed to evaluate I 2 , which can be written as

I 2 = ℓ=1,2 ρ ℓ h ℓ (L (2N +2) ℓ φℓ , φℓ ) L 2 - ℓ=1,2 ρ ℓ h ℓ (L (N ) ℓ ϕ ℓ , ϕ ℓ ) L 2 =: I 2,1 + I 2,2 .
In view of (3.8), we see that • ϕ 2 -ρ 1 l

I 2,1 = ρ 1 h 1 (L (2N +2) 1,0 φ1 , l (2N +2) 1 • φ1 ) L 2 + ρ 2 h 2 (L ( 2N
(N ) 1

• ϕ 1 ) L 2 = ρ 1 h 1 (L

(2N +2) 1,0 φ1 , l (N ) 1 • ϕ 1 ) L 2 + ρ 2 h 2 (L (2N +2) 2,0 φ2 , l (N ) 2 • ϕ 2 ) L 2 = ρ 1 h 1 N i=0 2N +2 j=0 (L 1,ij φ1,j , ϕ 1,i ) L 2 + ρ 2 h 2 N * i=0 2N * +2 j=0 (L 2,ij φ2,j , ϕ 2,i ) L 2 = ρ 1 h 1 N i=0 2N +2 j=0 (L 1,ji ϕ 1,i , φ1,j ) L 2 + ρ 2 h 2 N * i=0 2N * +2 j=0 (L 2,ji ϕ 2,i , φ2,j ) L 2 ,
where we used L * ℓ,ij = L ℓ,ji . Similarly, we see also that 2,0 ϕ 1 = N * i=0 L 2,ji ϕ 2,i hold only for j = 0, 1, . . . , N and for j = 0, 1, . . . , N * , respectively. Therefore, we have Hence, denoting by φ 1 = (φ 1,0 , φ 1,1 , . . . , φ 1,N ) T and φ 2 = (φ 2,0 , φ 2,1 , . . . , φ 2,N * ) T with φ ℓ,i = ϕ ℓ,i -φℓ,i we obtain Here, we note that (φ 1 , φ 2 ) satisfy

I 2,2 = ρ 1 h 1 (L (N ) 1,0 ϕ 1 , l (N ) 1 • ϕ 1 ) L 2 + ρ 2 h 2 (L (N ) 2,0 ϕ 2 , l (N ) 2 • ϕ 2 ) L 2 = (h 2 L (N ) 2,0 ϕ 2 , ρ 2 l (N ) 2 • ϕ 2 -ρ 1 l (N ) 1 • ϕ 1 ) L 2 = (h 2 L (N ) 2,0 ϕ 2 , ϕ) L 2 = (h 2 L (N ) 2,0 ϕ 2 , ρ 2 l (2N +2) 2 • φ2 -ρ 1 l (2N +2) 1 • φ1 ) L 2 = ρ 1 h 1 (L (N ) 1,0 ϕ 1 , l (2N +2) 1 • φ1 ) L 2 + ρ 2 h 2 (L (N ) 2,0 ϕ 2 , l (2N +2) 2 • φ2 ) L 2 = ρ 1 h 1 2N +2 j=0 (H 2j 1 L (N )
I 2,2 = ρ 1 h 1 N i=0 N j=0 (L 1,
|I 2 | ≲ ℓ=1,2 ρ ℓ h ℓ (∥∇φ ℓ ∥ 2 L 2 + (h ℓ δ) -2 ∥φ ′ ℓ ∥ 2 L 2 ) + ρ 1 h 1 ∥( φ1,
           L (N )
1,i φ 1 = r 1,i for i = 0, 1, . . . , N, L j=N +1 (L 1,ij -H 2i L 1,0j ) φ1,j for i = 0, 1, . . . , N, r 2,i = -2N * +2 j=N * +1 (L 2,ij -H p i L 2,0j ) φ2,j for i = 0, 1, . . . , N * , ∇ • r 3,1 = 2N +2 j=N +1 L 1,0j φ1,j , ∇ • r 3,2 = 2N * +2 j=N * +1 L 2,0j φ2,j , r 4,1 = 2N +2 j=N +1 H 2j 1 φ1,j , r 4,2 = -2N * +2 j=N * +1 H p j 2 φ2,j . We put r ′ 1 = (0, r 1,1 , . . . , r 1,N ) T and r ′ 2 = (0, r 2,1 , . . . , r 2,N ) T . Then, with a suitable decomposition r ℓ = r high ℓ + (h ℓ δ) -2 r low ℓ for ℓ = 1, 2, and using the linearity of (5.1), we see by Lemma

8 Figure 1 . 2 :

 812 Figure 1.2: Articulation of the proofs.

(4. 1 )

 1 l(H) := (1, H p 1 , . . . , H p N * ) T and define L ij = L ij (H, b, δ) (i, j = 0, 1, . . . , N * ) by

2 ,

 2 . . . , N * . The following Lemma has been proved in [20, Lemmas 3.2 and 3.4]. Lemma 4.1. Let c, M be positive constants and m an integer such that m > n 2 + 1. There exists a positive constant C such that if ζ ∈ H m , b ∈ W m,∞ , and H = 1 + ζ -b satisfy

  , b, δ)ϕ in the case j = 2N + 1 and k = m -4(N + 1) is given in [20, Theorem 2.2] and proved in [20, Sections 8.1 and 8.2].

Lemma 4 . 3 .n 2 + 1 .

 4321 Let c, M be positive constants and m an integer such that m ≥ 4(N + 1) and m > We assume (H1) or (H2). There exists a positive constant C such that if ζ ∈ H m , b ∈ W m+1,∞ , and H = 1 + ζ -b satisfy (4.8), then for any ϕ ∈ Hm and any δ ∈ (0, 1] we have

ℓ=1, 2 ρ 2 ρ

 22 ℓ Ω ℓ,0 |I δ ∇ X Φ res ℓ | 2 -I δ ∇ X Φ res ℓ • (I -P ℓ )I δ ∇ X Φℓ dX = ℓ=1,2 ∂Ω ℓ,0 ρ ℓ Φ res ℓ (N ℓ • I δ P ℓ I δ ∇ X Φℓ -N ℓ • I 2 δ ∇ X Ψ ℓ )dS = ℓ=1,2 R n ρ 1 Φ res 1 {(e z • I δ P 2 I δ ∇ X Φ2 -e z • I 2 δ ∇ X Ψ 2 ) -(e z • I δ P 1 I δ ∇ X Φ1 -e z • I 2 δ ∇ X Ψ 1 )} z=0 dx = 0, so that we obtain ℓ=1,2 ρ ℓ Ω ℓ,0 |I δ ∇ X Φ res ℓ | 2 dX = ℓ=1,ℓ Ω ℓ,0 I δ ∇ X Φ res ℓ • (I -P ℓ )I δ ∇ X Φℓ dX.

(H 1 (H 2

 12 (H 1 , δ, h 1 ) and L (N ) 2,i (H 2 , b, δ, h 2 ), respectively, in the following argument. For given (ζ, ϕ), let Φ be the solution to the transmission problem (8.1) with r S = 0 and let (ϕ 1 , ϕ 2 ) and ( φ1 , φ2 ) be the solutions to the problems (H 1 , δ, h 1 )ϕ 1 = 0 for i = 1, 2, . . . , N,L (N ) 2,i (H 2 , b, δ, h 2 )ϕ 2 = 0 for i = 1, 2, . . . , N * , h 1 L (N ) 1,0 (H 1 , δ, h 1 )ϕ 1 + h 2 L (N ) 2,0 (H 2 , b, δ, h 2 )ϕ 2 = 0, ρ 2 l (N ) 2 (H 2 ) • ϕ 2 -ρ 1 l (N ) 1 (H 1 ) • ϕ 1 = ϕ , δ, h 1 ) φ1 = 0 for i = 1,2, . . . , 2N + 2, , b, δ, h 2 ) φ2 = 0 for i = 1, 2, . . .

  φ 2 = r 2,i for i = 0, 1, . . . , N * , h 1 L (N ) 1,0 φ 1 + h 2 L (N ) 2,0 φ 2 = ∇ • (h 1 r 3,1 + h 2 r 3,2 ), ρ 2 l (N ) 2 • φ 2 -ρ 1 l (N ) 1 • φ 1 = ρ 1 r 4,1 + ρ 2 r 4,2 , i = -2N +2

  Kakinuma model that we are going to analyze in this paper, and then rewrite them in a nondimensional form by introducing several nondimensional parameters. Hamiltonians of the full model and of the Kakinuma model in the nondimensional variables are also provided. In Section 3 we first introduce some differential operators, which enable us to write the Kakinuma model simply in the form(1.8), and then we present the precise statements of our main results in this paper. In Section 4 we first recall results in the framework of surface waves related to the consistency of the Isobe-Kakinuma model, and then prove Theorems 3.5 and 3.6 concerning the consistency of the Kakinuma model by a simple scaling argument. In Section 5 we first derive an elliptic estimate related to the compatibility conditions for the Kakinuma model, which explains how to prepare the initial data, as stated in Proposition 3.4. Then we give uniform a priori bounds on regular solutions to the Kakinuma model, especially, a priori bounds of time derivatives. In Section 6 we provide uniform energy estimates for the solution to the Kakinuma model and prove Theorem 3.1, which ensures the existence of the solution to the initial value problem for the Kakinuma model on a time interval independent of parameters, especially, δ 1 and δ 2 , under the stability condition, the non-cavitation assumptions, and intrinsic compatibility conditions on the initial data, together with a uniform bound of the solution. In Section 7 we first give a supplementary estimate on an approximation of the Dirichlet-to-Neumann map, and then revisit the consistency of the Kakinuma model. We prove Proposition 7.6, which is another version of the consistency given in Theorem 3.6, where we adopt a different construction of an approximate solution to the Kakinuma model from the solution to the full model. Then, by making use of the well-posedness of the initial value problem for the Kakinuma model we prove Theorem 3.9 which provides a conditional rigorous justification of the Kakinuma model, that is, assuming the existence of a solution to the full model with a uniform bound we derive an error estimate between a corresponding solution to the Kakinuma model and that of the full model. Finally, in Section 8 we prove Theorem 3.10 which gives an error estimate between the Hamiltonian of the Kakinuma model and that of the full model. For the convenience of the reader, the structure of the paper and proofs dependencies are sketched in Figure1.2.Notation We denote by W m,p the L p Sobolev space of order m on R n and H m = W m,2 . We put Hm = {ϕ ; ∇ϕ ∈ H m-1 }. The norm of a Banach space B is denoted by ∥ • ∥ B . The L 2 -inner product is denoted by (•, •) L 2 . We put ∂ t = ∂ ∂t , ∂ j = ∂ x j = ∂ ∂x j , and ∂ z = ∂ ∂z .

  t), and H ℓ (ℓ = 1, 2) are depths of the upper and lower layers defined by(2.19). We note that if the canonical variables (ζ, ϕ) are given, then the Kakinuma model (2.18) determines ϕ 1 = (ϕ 1,0 , ϕ 1,1 , . . . , ϕ 1,N ) T and ϕ 2 = (ϕ 2,0 , ϕ 2,1 , . . . , ϕ 2,N * ) T , which are unique up to an additive constant of the form (Cρ 1 , Cρ 2 ) to (ϕ 1,0 , ϕ 2,0 ). For details, we refer to [14, Subsection 8.1] and Lemma 5.1 in Section 5. Then, the Hamiltonian H K (ζ, ϕ) of the Kakinuma model is given by(2.25)

  Remark 3.7. The unique existence of the solutions ϕ 1 and ϕ 2 to the problems (3.19) is guaranteed by Lemma 4.4 below under an additional assumption ϕ 1 (•, t), ϕ 2 (•, t) ∈ Hm . Lemma 4.4 is essentially a simple corollary of [20, Lemma 3.4]. Remark 3.8. In order to define the approximate solution (ϕ 1 , ϕ 2 ) to the Kakinuma model (2.18) from the solution (ζ, ϕ 1 , ϕ 2 ) to the full model, we can use, in place of (3.19), the following system

	of equations
	(3.21)

2 H m-1 for t ∈ [0, T ].

  4.1 and Lemma 4.2. As for the first term, in view of m > n 2 we can use a calculus inequality ∥r 2 ∥ H k ≲ ∥r∥ 2 H (m+k)/2 for k ∈ {0, 1, . . . , m}. Particularly, we have ∥r 2 ∥ H m-4(N +1) ≲ ∥r∥ 2 H m-2(N +1) . The last term can be evaluated by estimates in [20, Sections 8.1 and 8.2].

8

  Approximation of Hamiltonians; proof of Theorem 3.10As was shown in the companion paper[START_REF] Duchêne | A mathematical analysis of the Kakinuma model for interfacial gravity waves. Part I: Structures and well-posedness[END_REF] Theorem 8.4], the Kakinuma model (2.18) enjoys a Hamiltonian structure analogous to the one exhibited on the full model for interfacial gravity waves by T. B. Benjamin and T. J. Bridges in[START_REF] Benjamin | Reappraisal of the Kelvin-Helmholtz problem. Part 1. Hamiltonian structure[END_REF]. In this section, we will prove Theorem 3.10, which states that the Hamiltonian H K (ζ, ϕ) of the Kakinuma model defined in (2.25) approximates the Hamiltonian H IW (ζ, ϕ) of the full model defined in(2.22) with an error of order O

  2, we obtain ℓ=1,2ρ ℓ Ω ℓ,0 I δ ∇ X Φ res ℓ • P ℓ I δ ∇ X Φ res ℓ = ℓ=1,2 ρ ℓ Ω ℓ,0 I δ ∇ X Φ res ℓ • (I -P ℓ )I δ ∇ X Ψ ℓ dX.It follows from these two identities thatℓ=1,2 ρ ℓ ∥I δ ∇ X Φ res ℓ ∥ 2 L 2 (Ω ℓ,0 ) ≲ min ℓ=1,2 ρ ℓ ∥I δ ∇ X Φℓ ∥ 2 L 2 (Ω ℓ,0 ) , ℓ=1,2 ρ ℓ ∥I δ ∇ X Ψ ℓ ∥ 2 L 2 (Ω ℓ,0 ) ,

	which yields the equivalence
	ℓ=1,2

  where we used(2.15). In view of tanh ξ ≃ (1 + ξ) -1 ξ for ξ ≥ 0 and the relation (2.14), we haveσ(ρ 1 ρ 2 (ρ 1 Λ 2,0 + ρ 2 Λ 1,0 ) -1 ) ≃ ρ 1 ρ 2 h 1 h 2 (1 + h 1 δ|ξ|)(1 + h 2 δ|ξ|) (1 + δ|ξ|)|ξ| 2

						1,0 ) ρ 1	,	σ(Λ 2,0 ) ρ 2
			≤ min	h 1 ρ 1	,	h 2 ρ 2	|ξ| 2
			≤ 2|ξ| 2 ,
	≲ min	ρ 1 h 1	ρ 2	1 + h 1 δ|ξ| |ξ| 2	,	ρ 2 h 2	ρ 1	1 + h 2 δ|ξ| |ξ| 2
	≲ min	ρ 1 h 1	(|ξ| -1 + h 1 δ) 2 ,	ρ 2 h 2	(|ξ| -1 + h 2 δ) 2 ,

  , 2N * + 2, On the other hand, it follows from Lemmas 8.1, 4.5, 7.2, and 5.1 that ℓ=1,2ρ ℓ ∥I δ ∇ X Φ res,2 ℓ h ℓ (h ℓ δ) 2(4N +2) ∥∇( Φapp ℓ | z=ζ )∥ 2 ℓ h ℓ (h ℓ δ) 2(4N +2) (∥∇ φℓ ∥ 2 H 4N +3 + ∥ φ′ ℓ ∥ 2 H 4N +3 ) ≲ ((h 1 δ) 4N +2 + (h 2 δ) 4N +2 ) 2 ∥∇ϕ∥ 2 H 4N +3 .Summarizing the above estimates, we obtain|I 1 | ≲ ((h 1 δ) 4N +2 + (h 2 δ) 4N +2 )∥∇ϕ∥ H 4N +3 ∥∇ϕ∥ L 2 .

	h 1 L (2N +2) 1,0 ρ 2 l (2N +2) 2 (H 2 ) • φ2 -ρ 1 l (H 1 , δ, h 1 ) φ1 + h 2 L (2N +2) 2,0 (2N +2) 1 (H 1 ) • φ1 = ϕ, (H 2 , b, δ, h 2 ) φ2 = 0, respectively, and define (Φ app 1 , Φ app 2 ) and ( Φapp 1 , Φapp 2 ) by (2.24) and        Φapp 1 (x, z) = 2N +2 i=0 (1 -h -1 1 z) 2i φ1,i (x), ∥ 2 L 2 (Ω ℓ ) ≲ min ℓ=1,2 ρ ℓ h ℓ ∥((-∆) -1 2 + h ℓ δ)r S ∥ 2 L 2 ≲ ℓ=1,2 ρ ℓ h ℓ ∥((-∆) -1 2 + h ℓ δ)(h ℓ Λ (2N +2) ℓ -Λ ℓ )[ Φapp ℓ | z=ζ ]∥ 2 L 2 ≲ ℓ=1,2 H 4N +3 ≲ 2N ℓ ℓ=1,2
	     	Φapp 2 (x, z) =

* +2 i=0 (1 + h -1 2 (z -b(x))) p i φ2,i (x), ρ ρ

  1,0 ϕ 1 , φ1,j ) L 2 + ρ 2 h 2 = N i=0 L 1,ji ϕ 1,i and H

	Here, it follows from (3.8) that H 2j 1 L (N ) 1,0 ϕ 1 p j 2 L (N )
	2N * +2
	(H 2 L p j (N ) 2,0 ϕ 2 , φ2,j ) L 2 .
	j=0
	58

  ji ϕ 1,i , φ1,j ) L 2 + ρ 2 h 2 L 1,0i ϕ 1,i , φ1,j ) L 2 + ρ 2 h 2 (L 1,ji -H 2j 1 L 1,0i )(ϕ 1,i -φ1,i ), φ1,j ) L 2

							N *	N *
								(L 2,ji ϕ 2,i , φ2,j ) L 2
							i=0	j=0
	+ ρ 1 h 1	N	2N +2	(H 2j 1 N *	2N * +2	(H	p j 2 L 2,0i ϕ 2,i , φ2,j ) L 2 ,
		i=0	j=N +1				i=0	j=N * +1
	so that
			N		2N +2
	I 2 = ρ 1 h 1				((L 1,ji -H 2j 1 L 1,0i )ϕ 1,i , φ1,j ) L 2
			i=0	j=N +1
				N *	2N * +2
		+ ρ 2 h 2				((L 2,ji -H 2 L 2,0i )ϕ 2,i , φ2,j ) L 2 p j
				i=0	j=N * +1
			N		2N +2
		= ρ 1 h 1
			i=0	j=N +1
					2N +2	2N +2
		-ρ 1 h 1				((L 1,ji -H 2j 1 L 1,0i ) φ1,i , φ1,j ) L 2
				i=N +1	j=N +1
					2N * +2	2N
		-ρ 2 h 2
				i=N * +1

(+ ρ 2 h 2 N * i=0 2N * +2 j=N * +1 ((L 2,ji -H p j 2 L 2,0i )(ϕ 2,i -φ2,i ), φ2,j ) L 2 * +2 j=N * +1 ((L 2,ji -H p j 2 L 2,0i ) φ2,i , φ2,j ) L 2 .

  N+1 , φ1,N+2 , . . . , φ1,2N+2 )∥ 2H 1 + ρ 2 h 2 ∥( φ2,N * +1 , φ2,N * +2 , . . . , φ2,2N * +2 )∥ 2 H 1 + ρ 1 h 1 (h 1 δ) -2 ∥( φ1,N+1 , φ1,N+2 , . . . , φ1,2N+2 )∥ 2 L 2 + ρ 2 h 2 (h 2 δ) -2 ∥( φ2,N * +1 , φ2,N * +2 , . . . , φ2,2N * +2 )∥ 2 L 2 .

  5.1 that ℓ=1,2ρ ℓ h ℓ (∥∇φ ℓ ∥ 2 L 2 + (h ℓ δ) -2 ∥φ ′ ℓ ∥ 2 L 2 ) + (h ℓ δ) -2 ∥r low ℓ ∥ 2 L 2 + ∥r 3,ℓ ∥ 2 L 2 + ∥r 4,ℓ ∥ 2 H 1 ) ≲ ρ 1 h 1 ∥( φ1,N+1 , φ1,N+2 , . . . , φ1,2N+2 )∥ 2 H 1 + ρ 2 h 2 ∥( φ2,N * +1 , φ2,N * +2 , . . . , φ2,2N * +2 )∥ 2 H 1 + ρ 1 h 1 (h 1 δ) -2 ∥( φ1,N+1 , φ1,N+2 , . . . , φ1,2N+2 )∥ 2 L 2 + ρ 2 h 2 (h 2 δ) -2 ∥( φ2,N * +1 , φ2,N * +2 , . . . , φ2,2N * +2 )∥ 2 L 2 .Moreover, it follows from [20, Lemmas 5.2 and 5.4] that∥( φ1,N+1 , φ1,N+2 , . . . , φ1,2N+2 )∥ H k ≲ (h 1 δ) 2N +2-k ∥∇ φ1 ∥ H 2N +1 ∥( φ2,N * +1 , φ2,N * +2 , . . . , φ2,2N * +2 )∥ H k ≲ (h 2 δ) 2N +2-k (∥∇ φ2 ∥ H 2N +1 + ∥ φ′ 2 ∥ H 2N +1 ),for k = 0, 2, and hence also for k = 1 by interpolation, so that|I 2 | ≲ ℓ=1,2 ρ ℓ h ℓ (h ℓ δ) 4N +2 (∥∇ φℓ ∥ 2 H 2N +1 + ∥ φ′ ℓ ∥ 2 H 2N +1 ) ≲ ((h 1 δ) 4N +2 + (h 2 δ) 4N +2 )∥∇ϕ∥ 2 4N +2 + (h 2 δ) 4N +2 )∥∇ϕ∥ H 4N +2 ∥∇ϕ∥ L 2 ,where we used Lemma 5.1 with (2.15), and interpolation. This completes the proof of Theorem 3.10.

	≲	ℓ=1,2	ρ ℓ h ℓ (∥r high ℓ	∥ 2 H -1 H 2N +1

≲ ((h 1 δ)

The basic equations and the Kakinuma model2.1 Equations with physical variablesWe first recall the equations governing potential flows for two layers of immiscible, incompressible, homogeneous, and inviscid fluids, and then write down the Kakinuma model at stake in this work. In the following, we denote the upper layer, the lower layer, the interface, the rigid-lid,

-ζ = r 0 ,

Acknowledgement T. I. was partially supported by JSPS KAKENHI Grant Number JP17K18742 and JP22H01133. V. D. thanks the Centre Henri Lebesgue ANR-11-LABX-0020-01 for creating an attractive mathematical environment.

then for any ϕ 1 , ϕ 2 ∈ Hk+2(j+1) with 0 ≤ k ≤ m -2(j + 1) we have

Proof. By simple scaling arguments, we have (4.15)

). Therefore, the results follow from Lemma 4.2.

Lemma 4.6. Let c, M be positive constants and m an integer such that m ≥ 4(N + 1) and m > n 2 + 1. We assume (H1) or (H2). There exists a positive constant C such that for any positive parameters h 1 , h 2 , δ satisfying

2 b satisfy (4.14), then for any ϕ 1 , ϕ 2 ∈ Hm we have

H m-1 (h 2 δ) 4N +2 . Proof. By simple scaling arguments, we have

δ). Therefore, the results follow from Lemma 4.3.

We can now prove Theorems 3.5 and 3.6. In view of (3.8) the errors (r 1 , r 2 , r 0 ) and (r 1 , r2 , r0 ) can be written explicitly as

Therefore, the theorems are simple corollaries of the above Lemmas 4.5 and 4.6.

Elliptic estimates and time derivatives

In this section we derive useful uniform a priori bounds on regular solutions to the Kakinuma model (2.18). Firstly, due to the fact that the hypersurface t = 0 in the space-time R n × R is characteristic for the Kakinuma model, we need the following key elliptic estimate in order to be able to estimate time derivatives of the solution. Let us recall that the operators L 1,i for i = 0, 1, . . . , N and L 2,i for i = 0, 1, . . . , N * are defined by (3.6), and the vectors l 1 (H 1 ) and l 2 (H 2 ) are defined by (3.3). We recall the convention that for a vector ϕ = (ϕ 0 , ϕ 1 , . . . , ϕ N ) T we denote the last N components by ϕ ′ = (ϕ 1 , . . . , ϕ N ) T . 7 Approximation of solutions; proof of Theorem 3.9

In this section we prove Theorem 3.9, which gives a rigorous justification of the Kakinuma model (2.18) as a higher order shallow water approximation to the full model for interfacial gravity waves (2.17) under the hypothesis of the existence of the solution to the full model with uniform bounds.

Supplementary estimate for the Dirichlet-to-Neumann map

In this subsection, we give a supplementary estimate to Lemma 4.2 for the Dirichlet-to-Neumann map Λ(ζ, b, δ) defined by (4.6) appearing in the framework of surface waves. We recall the map Λ (N ) (ζ, b, δ) : ϕ → L 0 (H, b, δ)ϕ, where L 0 (H, b, δ) is defined by (4.3) and ϕ is the unique solution to (4.5). In this section we omit the dependence of t in notations.

Lemma 7.1. Let c, M be positive constants and m, j integers such that m > n 2 +1, m ≥ 2(j +1), and 1 ≤ j ≤ 2N + 1. We assume (H1) or (H2). There exists a positive constant

, and H = 1 + ζ -b satisfy (4.8), then for any ϕ ∈ Hk+2(j+1) with 0 ≤ k ≤ m -2(j + 1) and any δ ∈ (0, 1] we have

Proof. This lemma can be proved in a similar way to the proof of Lemma 4.2 with a slight modification. For the completeness, we sketch the proof. By the duality (H k ) * = H -k and the symmetry of the operator (-∆) -1 2 , it is sufficient to show the estimate

for any ϕ ∈ Hk+2(j+1) and any ψ ∈ H 1-k . We decompose it as

and evaluate the two components of the right-hand side separately. We recall the definitions (4.1) of the (N * + 1) vector-valued function l(H) and (4.3) of the operator L i (H, b, δ), which acts on (N * + 1) vector-valued functions. These depend on N , so that we denote them by l (N ) (H) and L (N ) i (H, b, δ), respectively, in the following argument. Let Φ be the solution to the boundary value problem (4.7) and let ϕ = (ϕ 0 , ϕ 1 , . . . , ϕ N * ), φ = ( φ0 , φ1 , . . . , φ2N * +2 ), and ψ = (ψ 0 , ψ 1 , . . . , ψ 2N * +2 ) be the solutions to the problems

Let c, M be positive constants and m, j integers such that m > n 2 +1, m ≥ 2(j +1), and 1 ≤ j ≤ 2N + 1. We assume (H1) or (H2). There exists a positive constant C such that for any positive parameters

We recall also the estimate for the Dirichlet-to-Neumann map Λ(ζ, b, δ) itself. The following lemma is now standard. For sharper estimates, we refer to T. Iguchi and D. Lannes [START_REF] Lannes | The water waves problem: mathematical analysis and asymptotics[END_REF].

This lemma and the scaling relations (4.15) imply immediately the following lemma.

Lemma 7.4. Let c, M be positive constants and m an integer such that m > n 2 + 2. There exists a positive constant C such that for any positive parameters

Consistency of the Kakinuma model revisited

As we mentioned in Remark 3.8, the approximate solution to the Kakinuma model (2.18) made from the solution (ζ, ϕ 1 , ϕ 2 ) to the full model can be constructed as a solution to (3.21), that is,

To show this fact, we need to guarantee that the difference between these two solutions is of order O((h 1 δ) 4N +2 + (h 2 δ) 4N +2 ). The following lemma gives such an estimate.

Lemma 7.5. Let c, M be positive constants and m an integer such that m > n 2 + 1 and m ≥ 4(N + 1). We assume (H1) or (H2). There exists a positive constant C such that for any positive parameters

for |k| ≤ m. We proceed to evaluate (∂ t ϕ res 1 , ∂ t ϕ res 2 ). We recall that (

2 ) and note that, differentiating the first three equations of (7.2) with respect to t and using the last equation in (7.4) ) and f 4 added with the error term -r 0 . By taking the difference between these equations, we have therefore

, where a ℓ,0 (H ℓ ) and b 2,0 (H 2 ) are the 0th columns of the matrixes A ℓ (H ℓ ) and B 2 (H 2 ) defined by (6.2) and (6.4), respectively, and so on. Note the relations

2 ∇b). Therefore, by Lemma 5.1 we have, for

We will evaluate each term in the right-hand side. For 1 ≤ k ≤ m -1, we see that

and

Moreover, for any 0 ≤ k ≤ m we have also (7.6) ℓ=1,2

Summarizing the above estimates and using

)

We need also to evaluate

Now, by deriving equations for spatial derivatives of (ζ res , ϕ res 1 , ϕ res 2 ) and applying the energy estimate obtained in Subsection 6.1 we will evaluate E res k . Let β be a multi-index such that 1 ≤ |β| ≤ k. Applying ∂ β to the Kakinuma model (3.5) for (ζ K , ϕ K 1 , ϕ K 2 ) and to (7.4) for (ζ IW , φIW 1 , φIW 2 ) and taking the difference between the resulting equations, we obtain

respectively. Then, by the definitions of the Hamiltonian functionals H IW (ζ, ϕ) and H K (ζ, ϕ) given in Section 2.3, we have

We will evaluate I 1 and I 2 , separately.

In order to evaluate I 1 , we put Φ res ℓ = Φ ℓ -Φapp ℓ (ℓ = 1, 2), so that

≲ ∥∇ϕ∥ 2 L 2 , where we used Lemma 5.1 and order to evaluate

Here, we note that R ℓ (ℓ = 1, 2) can be written the form

56

Estimates for the residuals (r 1,0 , r 1,1 , . . . , r 1,2N +2 ), (r 2,0 , r 2,0 , . . . , r 2,2N * +2 ), and r B were given in [20, Lemmas 6.4 and 6.9] and their proofs. In fact, we have

We decompose Φ res ℓ = Φ res,1 ℓ + Φ res,2 ℓ

, where (Φ res,1 1 , Φ res,1

2

) is a unique solution to the problem

≲ h 1 (h 1 δ) 2 ∥(r 1,0 , r 1,1 , . . . , r 1,2N +2 )∥ L 2 ≲ h 1 (h 1 δ) 2(4N +3) ∥∇ φ1 ∥ 2

and that