Two-dimensionality of gravity water flows governed by the equatorial f-plane approximation

We show that gravity wave trains governed by the equatorial f-plane approximation propagate at the free surface of a rotational water flow of constant vorticity vector $$(\Omega _1, \Omega _2, \Omega _3)$$(Ω1,Ω2,Ω3) over a flat bed only if the flow is two-dimensional. Owing to the presence of Coriolis effects, our result is also true even if the vorticity vector vanishes. This represents a striking difference when compared with the cases without geophysical effects discussed in Constantin (Europhys Lett 86:29001, 2009, Eur J Mech 30:12–16; 2011) and Martin (J Math Fluid Mech 2016. doi:10.1007/s00021-016-0306-1), where the conclusion about the two-dimensionality of the flow was possible under the assumption of constant nonvanishing vorticity vector. Another upshot is that the only nonzero component of the vorticity that may not vanish is $$\Omega _2$$Ω2, that is, the one pointing in the horizontal direction orthogonal to the direction of wave propagation.


Introduction
There is an increased interest nowadays in the study of geophysical water flows. These are flows whose motion is influenced by the Earth's rotation reflected through the manifestation of the Coriolis force. The thorough understanding of geophysical flows is highly dependent upon inclusion of nonlinear aspects that considerably complicate the matter to the point of eluding exact treatment; this drawback is somehow mitigated by the discovery by Gerstner [17] who presented in 1802 the only known explicit exact solution (in Lagrangian variables) for the governing equations for periodic homogeneous two-dimensional traveling gravity B Calin Iulian Martin calin.martin@univie.ac.at 1 Faculty of Mathematics, University of Vienna, Oskar-Morgenstern-Platz 1, 1090 Vienna, Austria water waves. Gerstner solution was extended in 1932 to a heterogeneous fluid by Dubreil-Jacotin [13]; we refer the reader to Constantin [1] and Henry [18] for modern presentations of the Gerstner wave solution.
The specific general difficulties raised by the nonlinear character of equations governing the water wave propagation are enhanced by the inclusion of the before mentioned Coriolis effects; especially the complex dynamics of flows near the equatorial Pacific presents some peculiarities like a strong stratification (evidentiated by the presence of an interface, called thermocline, separating two layers of constant density), the presence strong depth-dependent underlying currents, cf. [9], with flow reversal at a depth of about 100-200 m and a wide range of other wave propagation phenomena. The softening of these intricacies becomes possible through the employment of suitable approximations of the geophysical governing equations; one of these is the f -plane approximation which is appropriate for oceanic flows within a restricted meridional range of about 2 • latitude from the Equator [9,12,21].
The rigorous mathematical study of geophysical water flows has resulted in the appearance of a bulk of papers commencing with the modeling of wave-current interactions in the fplane approximation for underlying currents of constant vorticity, by Constantin [4] and with the derivation of an exact solution (in the Lagrangian framework) of the nonlinear geophysical water waves in the β-plane near the Equator by Constantin in [5] by a suitable and substantial modification of the Gerstner wave solution to the case of three-dimensional flows. These studies were followed by papers presenting explicit exact solutions describing flows accommodating flow reversal and allowing for underlying currents in the β-plane setting, cf. [7,8,19]. An additional aspect of centrifugal forces in the β-plane approximation was addressed recently by Henry in [22] where an exact and explicit solution describing equatorially trapped waves was obtained.
Circumscribed to the same vein of solution flows exhibiting a prefered direction of propagation are recent works describing oceanic flows propagating in an azimuthal direction and exhibiting a complicated vertical structure [10,11,25]. Such geophysical flows were also analyzed from the perspective of stability/instability properties [6,16,20,23,24].
Our intention here is to show that gravity wave trains can propagate at the free surface of a rotational water flow of constant vorticity and governed by the equatorial f -plane approximation only if the flow has a two-dimensional character. From a historical perspective, the property that the flow below a wave train does not vary in the horizontal direction perpendicular to the direction of motion was validated by previous works of Constantin and Kartashova [2] proving the two-dimensionality of constant vorticity water flows below wave trains driven by capillarity and of Constantin [3] showing the two-dimensionality of constant vorticity gravity flows below wave trains; see also the works of Wahlén [30] and Stuhlmeier [29], the latter in the case of solitary waves.

The governing equations
We choose a rotating framework with the origin at a point on the Earth's surface, with the x axis pointing horizontally due East, the y-axis horizontally due North, while the z-axis is oriented upward. We will consider here regular wave trains of water waves propagating steadily in the direction of the horizontal x axis, periodic (of period L) in the variable x, and exhibiting no variation in the y direction. We will, in fact, show that the lack of the y-dependence of the free surface manifests throughout the entire flow.
The fluid domain is bounded below by the impermeable flat bed z = −d and above by the free surface z = η(x − ct), where the function η gives the wave profile and c > 0 is the wave speed.
If (u, v, w) denotes the velocity field, t is the time variable, P the pressure, ρ the constant density, g is the constant acceleration of gravity, and ω = 73 × 10 −6 rad/s is the constant rotational speed of Earth round the polar axis, the governing equations in the f -plane approximation near the Equator are the Euler's equations and the equation of mass conservation To single out the water wave problem from a vast range of hydrodynamical ones we impose the kinematic boundary conditions and together with the dynamic boundary condition To model observed phenomena, like underlying currents, it is essential to include the vorticity (x, y, z, t), obtained as the curl of the velocity field, that is, (2.6) The evolution of the vorticity is governed by the equation Remark 2.1 Throughout the paper we shall work under the assumption that 2 + 2ω = 0. This is a reasonable assumption, since ω ≈ 0.73 × 10 −4 rad/s, while a typical value for 2 in the equatorial Pacific is 25 × 10 −3 s −1 , cf. [4].
After these preliminary considerations we are able to formulate the main result of the paper. Proof We start the proof by noticing that, due to the assumption of constant vorticity, the evolution equation (2.7) takes the form We are going to use the third equation in the system above, which can be written as to prove that 3 = 0. We will resort to the method of proof by contradiction and assume that 3 = 0. Consequently, w is constant in the direction of the vector ( 1 , 2 + 2ω, 3 ), which, owing to 3 = 0, is not parallel to the flat bed z = −d. We then use the kinematic condition (2.4) and argue as in [3] to conclude that w = 0 within the fluid domain. Moreover, w = 0 implies that u z = 2 and v z = − 1 . The latter two relations yield the existence of two functionsũ =ũ(x, y, t) andṽ =ṽ(x, y, t) such that The advantage of working withũ,ṽ emerges as the possibility of deriving from (2.2) the existence of a function ψ(x, y, t) satisfying The latter equation can be used to rewrite 3 = v x − u y as Moreover, the first two equations of (2.9) can be written in terms of the function ψ as Note that the system (2.13) implies that 2 1 + ( 2 + 2ω) 2 > 0. Therefore, we obtain relations from which we infer that where Since η is not identically zero and the equality in (2.16) is true for all x it follows from above that F = 0. Since 2 + 2ω = 0, it necessarily follows from F = 0 that 2 = 0, because we are also in the hypothesis that 3 = 0. Moreover, we also have that E = 0 and D = − 1 2 3 . Consequently, from relation (2.16) we conclude that for all x it holds (e(t) − c)η (x − ct) = 0, which implies that e(t) ≡ c, and therefore From the third of Euler's equations (2.1) we conclude that which implies that for some function p(x, y, t). From (2.17) and the first of the Euler's equations we obtain that p x = P x = 0, while using also the second of the Euler's equations we get that p y = P y = ρ(d (t) − c 3 ). Consequently, for some functionf . The dynamic boundary condition (2.5) can now be written We infer from above that the coefficient of y has to vanish and thatf and η have to be constant functions. But η constant is equivalent with our flow having a flat surface, which is a contradiction.
The contradiction that we reached shows that 3 = 0. In the following, we will be concerned with proving that 1 also vanishes. Moreover, we will show that the pressure P and the velocity field (u, v, w) are independent of y.
Our next goal is to show that 1 = 0 and for this purpose we assume for the sake of contradiction that 1 = 0.
We note first that from (2.8) and using that 3 = 0 we obtain 1 u x + ( 2 + 2ω)u y = 0, On the other hand, from 0 = 3 = v x − u y we obtain a function ϕ = ϕ(x, y, z, t) which satisfies Up to addition to ϕ of a function that depends only on z and t we infer from the expressions of 1 and 2 that With the help of the above expressions of u, v, w in terms of ϕ we see that the equation of mass conservation (2.2) becomes ϕ x x + ϕ yy + ϕ zz = 0. (2.22) The harmonicity of w emerges now as a direct consequence of differentiating (2.22) with respect to z. We are going to show now that w = 0 on a plane parallel to the flat bed and contained entirely in the fluid domain. To be more precise, let be a trough line located at some position x 0 at some moment t 0 . Using the kinematic condition (2.3) and the fact that η x = 0 along l we conclude that w = 0 on the whole of l. (2.23) The third relation in (2.21) allows us to conclude that w is constant in the direction of the vector ( 1 , 2 + 2ω, 0) which, due to the assumption 1 = 0, is not collinear with l. Thus, w = 0 at all points of the plane containing l and spanned by l and the vector ( 1 , 2 +2ω, 0), which is in fact the plane z = η(x 0 , t) and which is parallel to the bed z = −d and contained entirely in the fluid domain.
Using that w is a bounded harmonic function that vanishes on z = −d and on z = η(x 0 , t) we conclude from the Phragmen-Lindelöf maximum principle, cf. [15], that w = 0 in the region of the fluid domain bounded by the flat bed and the plane z = η(x 0 , t). We can also infer from [14] that, owing to its harmonicity, w is real-analytic at any instant t. Thus, since w = 0 on an open set, it must vanish on the entire fluid domain. Then, as in the beginning of the proof we infer the existence of a function ψ(x, y, t) satisfying relations (2.10) and (2.11). In addition, since now 3 = 0, we have for the function ψ the relations for all x, y, t. We rely once more on the kinematic condition (2.3) to derive the relation From the latter we infer first that 2 = 0 and then, as a consequence, it follows that e(t) ≡ c.
We follow now the same route as in the beginning of the proof and reach that η has to be a constant. The latter contradiction shows that 1 necessarily has to vanish. We are going to prove now that the velocity field does not depend on the y variable. Moreover, we will prove that v is in fact a constant and the y-derivative of the pressure function vanishes. To this end note that the vorticity equation (2.8) takes the form ( 2 + 2ω)u y = ( 2 + 2ω)v y = ( 2 + 2ω)w y = 0, (2.26) which, due to 2 + 2ω = 0, implies u y = v y = w y = 0. Using also 0 = 1 = w y − v z and 0 = 3 = v x − u y we infer that v x = v z = 0. Thus, the second component of the velocity field depends only on t.
A consequence of the previous findings is that the second equation in (2.1) reduces to v t = − 1 ρ P y .
Integrating with respect to y in the above relation we find that P(x, y, z, t) = −ρv t y + F(x, z) + R(t), for some functions F and R whose variables are as indicated in the formula above. With the latter structure for P, the dynamic boundary condition (2.5) becomes −ρv t y + F(x, η(x)) + R(t) = P atm , for all x, y and t. But this implies that v t ≡ 0. Therefore, P y ≡ 0. We have thus proved that the flow has a two-dimensional character.
Remark 2.3 Note that, due to the assumption 2 + 2ω = 0-reasonable in the setting of equatorial waves-we can dispense of the constant nonzero vorticity requirement made in [2] in the context of capillary wave trains, in [3] for the case of gravity wave trains and in [26] regarding capillary-gravity water waves.