A steady, purely azimuthal flow model for the Antarctic Circumpolar Current

We consider a steady, purely azimuthal eddy viscosity flow model for the Antarctic Circumpolar Current and derive an exact formula for the solution.

the reader to [10] for the analysis of an ACC model, which establishes a connection between the free surface and the pressure.
In this paper we consider a purely azimuthal divergence free flow, which is governed by the incompressible Euler equations in the f -plane approximation near the 45th parallel south. The action of mesoscale eddies that transfer the surface stress to the bottom (see [12]) is taken into account by adding an eddy-viscosity term to the system. Furthermore we assume the presence of a constant wind stress and a constant atmospheric pressure on the ocean surface. On the flat bottom we impose a no-slip boundary condition. This system can be solved exactly and the solution is written in terms of the viscosity function, cf. Sect. 3.

The model
As a starting point we consider the f -plane approximation of Euler's equations for three-dimensional divergence-free ocean flows in the region of the ACC. Generally, for a fixed specified polar angle θ , where θ takes values between −π/2 and π/2 (e.g. we have that θ = 0 along the Equator, θ = −π/2 at the South Pole and θ = π/2 at the North Pole) the f -plane approximation is given by where Here t denotes the time and x, y, z denote the directions of increasing azimuth, latitude and vertical elevation respectively. Similarly we denote by u = u(x, y, z, t), v = v(x, y, z, t) and w = w(x, y, z, t) the velocity components of the flow field in direction of increasing azimuth, latitude and elevation. Furthermore P = P(x, y, z, t) denotes the pressure field. The constant Ω ≈ 7.29·10 −5 rad s −1 denotes the rotational speed of the Earth around the polar axis and g ≈ 9.81 m s −2 is the gravitational acceleration. The fluid density ρ is taken to be constant, since there is low stratification in the Polar region [9,18]. For modeling the ACC we take θ = −π/4, which corresponds to the 45th parallel south. Furthermore we restrict our considerations to the special case of purely azimuthal flows, i.e. we assume that v = 0 and w = 0.
(2.2) Thus (2.1) reduces to the linear system In order to include the transfer of the wind-generated surface stress to the bottom, which is due to the presence of mesoscale eddies [12], we add a viscosity term of the form (νu z ) z to (2.3) similarly as in [14] in the context of equatorial flows. The classical model of uniform eddy viscosity is due to [17]. We follow the more realistic approach with a depth dependent viscosity function as it was introduced in [8]. For solutions of the equatorial current system and their qualitative properties under the assumption of depth-dependent eddy viscosity we refer to [3,13,16] and the references therein. Furthermore we add an a priori unknown forcing term F to (2.4), cf. Remark 1, which will be determined in the analysis in Sect. 3. The viscosity function ν = ν(z), which is assumed to be smooth with |ν| > δ for some δ > 0, only depends on the depths of water [12], whereas F may depend on both y and z. (It is clear that F is independent of x.) Thus in summary we end up with the following system: Before introducing the boundary conditions for u and P, we set which is referred to as wind stress.
Let d denote the depth of water, which is taken to be constant, i.e. we assume a flat ocean bed. The water surface is considered to be flat, being located at z = 0. Thus the fluid domain consists of the layer D = R 2 × [−d, 0] ⊆ R 3 . On the boundary of D, which consists of the bottom {z = −d} and the surface {z = 0}, we impose the following conditions: (2.14) Here P atm denotes the constant atmospheric pressure at sea level and τ 0 is a constant related to the wind stress at the surface of the ocean.

Remark 1
We note that the usage of the forcing term F is necessary to ensure the existence of a nontrivial flow: otherwise we would end up with the relation u y = u z holding within the entire fluid domain. To see this we first differentiate (2.8) with respect to z. Next we differentiate (2.9) with respect to y and combine the respective outcomes. If we assume that F = 0, we get that u y = u z and the no-slip boundary condition (2.14) implies that u = 0 throughout D.
Our assumption that the wind-stress at the surface τ 0 is independent of both x and y simplifies the model and can not be justified globally, since the wind system in the region of the ACC is rather complex. In fact, the ACC is-in contrast to most other ocean currents-not a single flow, but a fragmented system of strong, regionally bounded jets [15]. Locally, within such a jet, it is reasonable to take τ 0 to be constant; cf. the real-time animations of ocean currents and winds on https://earth.nullschool. net/. To get a rough global picture, we may take an average value for τ 0 .
Let us furthermore emphasize that the linear nature of system (2.7)-(2.10) does not stem from approximations. This is a consequence of the reasonable assumption of a purely azimuthal flow satisfying (2.2).

An exact solution
In this section we provide an exact formula for the solution of the system (2.7)-(2.10), i.e. formulas for the horizontal velocity u and the pressure P. These results are summarized in the following theorem.

Remark 2
We illustrate the velocity profile of the flow field for a particular viscosity function ν in Fig. 1. The horizontal velocity u is generally strictly monotonous, since ν does not change sign by assumption, cf. Sect. 2. The formulas (3.1) and (3.2) in particular reveal that we are actually dealing with a two-dimensional model.
Proof of Theorem 1 By applying ∂ x on (2.9) and employing (2.10) we obtain that P xz = 0 in D. (3.4) We infer from (2.7) by differentiation with respect to z that Thus, by (3.4) we have that 0 = (νu z ) zz in D. (3.6) From (3.6) and (2.10) we obtain via integration with respect to z that and hence τ (y, z) = ρ (A(y)z + B(y)) in D (3.8) due to (2.11). By employing the boundary condition (2.13) we infer that τ (y, 0) = ρ B(y) = τ 0 for all y ∈ R. (3.9) Thus B is independent of y and takes the constant value By using (3.7) and (2.7) we obtain (3.11) and integration with respect to x yields (3.13) where the later expression is independent of x, we have that A (y) = 0, hence (3.14) Since (A 1 y + A 2 ) x + C(y, 0) = ρ −1 P atm for all x, y ∈ R according to (3.12), (3.14) and (2.12), we conclude that A 1 = A 2 = 0 and C(y, 0) = ρ −1 P atm for all y ∈ R. (3.15) In particular we get that P x = 0 in D. (3.16) By combining (3.7), (3.10), (3.14), (3.15) and (2.14) we obtain that u is independent of y and satisfies (3.1). In view of (2.9) we see that ρ −1 P z = √ 2Ω u(z) − g in D and integration with respect to z gives Due to the boundary condition (2.12) we infer that the constant of integration C(y) = C in (3.17) is independent of y and satisfies hence P is given by (3.2). We finally deduce from (2.8) that F satisfies (3.3), since P is independent of the y-coordinate. Remark 3 In contrast to equatorial flows, where stratification plays a dominant role and leads to two-layer models (see e.g. [3,4,6,13,16]), there is almost no stratification in the Southern Ocean [9,18]. Therefore we considered a one-layer model for describing the dynamics of the ACC, cf. Sect. 2. Let us note that the velocity profile of the solution u of system (2.7)-(2.14) (cf. Fig. 1) differs strongly from stratified flows in equatorial regions like the Equatorial Undercurrent, where u changes its sign (i.e. the direction of the flow gets reversed) at a certain depth between the surface and the thermocline that separates the two layers of different fluid density.
Remark 4 Let us finally point out another difference between equatorial flows and flows in the Southern Ocean like the ACC. While surface waves near the Equator are generally not higher than a few meters (see the discussion in [1]), there exist surface waves in the Southern Ocean having a wave height up to 30 m, cf. [19]. Investigations of such large surface waves and their interactions with underlying currents could thus be of particular interest in the near future. Another interesting question concerns the stability of waves. In [2] it was shown that instability of certain waves in equatorial regions does occur. On the other hand, short-wavelength stability results for certain equatorial flows are available in [11]. Similar stability results in the context of Southern Ocean currents would be of great interest.