Explicit two-dimensional solutions for the ocean flow in arctic gyres

We investigate the ocean flow in an arctic gyre. For given linear or constant oceanic vorticites, we prove explicit solutions which depend on both the polar and the azimuthal angles in the spherical coordinates associated with a rotating Earth. This enables us to identify resonant modes.


Introduction
In an overly-simplified model of our planet, the major ocean currents would be solely driven by the winds. In reality however, the Coriolis force, due to the rotation of the Earth, deflects these practically steady currents so that they produce in the Northern Hemisphere a clockwise circular motion and in the Southern Hemisphere counterclockwise circular paths. Furthermore, continental boundaries interrupt and divert their flow as well. These large-scale linked-up ocean currents rimmed by land masses are known as "gyres" [2]. There are three types of gyres: tropical, which form near the Equator (confined to either the Northern or the Southern hemisphere since gyres do not cross the Equator, which acts like a wave-guide due to the vanishing of the meridional component of the Coriolis force, thereby facilitating east-west flow propagation-Communicated by A. Constantin. B Susanna V. Haziot susanna.haziot@univie.ac.at 1 Fakultät für Mathematik, Universität Wien, Oskar-Morgenstern-Platz 1, 1090 Wien, Austria see the discussion in [7,8]), subtropical (between polar and equatorial regions) and subpolar (in the polar regions). In this paper we will focus on the Arctic basin. It is worth mentioning that the Arctic is very different from the Antarctic: whereas the Arctic is a semi-enclosed ocean almost completely surrounded by landmasses, the Antarctic is almost its complete geographical opposite, a single landmass encircled by a very powerful current known as the Antarctic Circumpolar Current (ACC) [20]; we refer to [9,18,19,22] for studies of the gyre-like flow of the ACC.
To study gyres, one can reduce their complexity by using some of their properties: as their vertical velocity is negligeable with respect to their horizontal one, they can be considered as being two-dimensional flows. This feature was exploited recently in [10] to derive a model for ocean gyres in spherical coordinates. In [3], the author uses the stereographical projection to transform this model into a planar elliptic boundary-value problem which he then reduces to a one-dimensional ODE by neglecting azimuthal variations. This permits him to provide explicit solutions for given oceanic vorticities.
In this paper we replace the stereographic projection by the Mercator projection, which reduces the model in [10] to a semilinear elliptic PDE that is simpler than the equation obtained recently in [3]. The advantage of this new transformation is that for linear oceanic vorticity functions, thanks to the periodicity of the azimuthal velocity, we can use Fourier series to separate our x and y variables (in a planar geometry context, as we can take advantage of the Mercator projection). This in turn provides us with explicit solutions which take into account possible azimuthal variations. The availability of explicit linear solutions opens up possibilities for the investigation of weakly nonlinear effects by means of perturbation theory, along the lines of the developments in [4,5] of the considerations made in [3].

The Mercator projection
We begin by recalling the recently derived model for gyres in spherical coordinates. [10] Let θ ∈ [0, π) be the polar angle, where θ = 0 corresponds to the South Pole. Our latitude angle is therefore θ − π 2 . Let ϕ ∈ [0, 2π) be the azimuthal (or longitude) angle. The polar and azimuthal velocity components of the horizontal flow on the Earth are given by respectively, where ψ(θ, ϕ) represents the stream function in spherical coordinates. The governing equation for gyres is given by where (θ, ϕ) = ψ(θ, ϕ) + ω cos(θ ) is associated with the vorticity of motion of the ocean relative to the Earth's surface, 2ω cos(θ ) is the spin vorticity due to the rotation of the Earth and F( − ω cos(θ )) is the oceanic vorticity, due to the motion of the ocean and specific to a particular gyre. Given ω and F, we have to solve the governing equation in a given spherical region whose boundary is ψ 0 (ϕ). For this, we turn to the Mercator projection. The Mercator projection maps the sphere onto the plane such that the north-south direction is the horizontal direction, the east-west direction is the vertical direction with the length of the Equator preserved, and all paths of equal compass bearing on the sphere are straight lines. It is conformal but distorts areas [12]. The change of variables satisfies the above mentioned properties. With this change of variables, the North Pole (θ = π) corresponds to x = −∞ and the Equator (θ = π 2 ) to x = 0. Therefore, x < 0 in the Northern Hemisphere with We can then rewrite the governing equation as the following semilinear elliptic PDE with boundary condition if the boundary of our spherical region corresponds to the parallel θ 0 = 2 arctan(e −x 0 ) ∈ π 2 , π , situated in the Northern Hemisphere ( Fig. 1).
Since u(x, y) is periodic in the y-variable, we can reformulate u as the Fourier series and therefore, if F is linear, the linear elliptic PDE reduces to an infinite number of boundary-value problems for second-order ODEs: when k ∈ Z and k = 0, and When u(x 0 , y) is a constant the boundary is defined as a streamline (i.e. a level set of ψ). For our general solution to be physically admissible, we require that our stream function be equal to some constant at the North Pole, and that the North Pole be a stagnation point. Therefore, u(x, y) needs to satisfy the following asymptotic conditions: uniformly in y ∈ R since the horizontal velocity is 1 sin(θ) ψ ϕ , −ψ θ = (u y , u x ) cosh(x).

Constant oceanic vorticity
Let us consider the case when the oceanic vorticity is constant, F ≡ γ with γ ∈ R. Therefore, (2)-(3) takes the form: with boundary conditions The general solution to (6) is Since (5) yields we must have Theorem 1 If u(x 0 , y) = ψ 0 for all y ∈ [0, 2π), then u k ≡ 0 for k = 0.
Proof We note that u(x 0 , y) = k∈Z u k (x 0 )e iky for k = 0 is only equal to some constant for all y ∈ [0, 2π) when u k (x 0 ) = 0 for all k = 0. Because of (11), this means u k ≡ 0 for all k = 0.
We are now ready to state and prove our main result. Proof We begin by looking at the case k = 0. We have: It can be easily verified that is a particular solution of (12) and therefore, the general solution can be written in the form: We claim that for δ = γ and β = −ω+γ ln(2), (13) satisfies the asymptotic conditions (5). Indeed, setting δ = γ , we get: For the first condition in (5) From the boundary condition (4), we get that We now look at the case k = 1 and k = −1. In order for u 1 (x) = c 1 e x + c 2 e −x to satisfy (5), we need to set c 1 = c 2 = 0. Therefore, for k = ±1, u k ≡ 0. Consequently, from the boundary condition (4) we can deduce that for α 1 = 0 or α −1 = 0, we have no physically relevant solution. From (11), we can see that (5) is satisfied for all |k| ≥ 2. In addition, using the boundary condition (4), we have u(x 0 , y) = k∈Z α k e iky from which we get

Remark 1
The problem with no y-dependence was considered in [3]. The result in [3] corresponds to the setting u k ≡ 0 for k = 0.
Concerning the physical relevance of the above considerations, the basic sources of oceanic vorticity are wind force and the gravitational forces due to the relative motions of the Moon, the Sun and the Earth in the form of the flood and ebb tidal currents (see the discussions in [6,16]). Both oceanic vorticities can be realisticallly regarded as non-zero constants (see [11,14]), with the sign depending on the prevalent wind direction, and, respectively, on whether the flow is of ebb or flood type. Let us note that wave-current interactions in flows with vorticity are a topic of great current interest (see the discussions in [6,11,13,15]), at the large scales that are relevant for gyre flows the presence of surface waves is not relevant.

Some special linear oceanic vorticities
Let us now consider the case of linear vorticity, setting F(u k ) = au k . (2) then takes the form: for k = 0 and for k = 0. The case for k = 0 has been shown in [3].
To begin with, we make the change of variable and thus, using the fact that 1 − z 2 = cosh −2 (x), (14) becomes We rewrite Setting (17) into (16) yields which, if l is a positive integer, is known as the associated Legendre Differential equation, see [1,17]. Therefore, the general solution to (18) is where P k l (z) denotes the associated Legendre function.
, where l is any positive integer, and for any b ∈ R, is the general solution to (2) with F(u) = au + b satisfying the asymptotic conditions (5) and the boundary condition (4), Proof To begin with, the associated Legendre polynomial for 0 ≤ k ≤ l is defined as: and we can therefore write where n ≥ 2, n ∈ N. Therefore, for (20) to satisfy the asymptotic conditions (5), we need to set c 2 = 0. From the boundary condition (4), we have thus concluding the proof for k > 0.
For negative k, the associated Legendre polynomial P k l (z) is defined as and therefore the proof follows directly from the one for k > 0. Finally, from the previous section, we know that the asymptotic conditions are satisfied for . From the boundary condition (4), we have

Some considerations about general linear oceanic vorticity
It is interesting to observe that (14) can be reduced to a Riccati's equation, which in turn is equivalent to solving an Ermakov equation. Indeed, assuming that the solution to (14) is of the form u k (x) = e f (x) , where f is an arbitrary function, solving (14) is equivalent to solving the following Riccati equation (see [17]) with g(x) = f (x). Using the following substitution and choosing s(x) such that we can write (22) as where We assume for now that such an s exists. A condition for its existence will be given later on.
We can now rewrite (25) as and using the substitution The solution to (28) is and therefore a particular solution to (14) is It remains to find s(x). From (24), we get the following second-order ODE and using the substitution γ (x) = s −1/2 (x) we can reduce (30) to the following Ermakov equation ( see [21]) The general solution to an Ermakov equation (see [21]) is of the form where p 0 (x) is a particular solution tõ Remark 2 If a = −l(l + 1), for some integer l ≥ 1, then from (19) we already know that a particular solution to (33) is the associated Legendre polynomial P k l (z) with z = tanh(x) and l = −1+ √ 1−4a 2 . We can therefore conclude that