Abstract
In the setting of steady wind-driven oceanic currents in a non-equatorial two-layer ocean with constant eddy viscosity in the lower layer and a general continuous eddy viscosity profile in the upper layer, we derive a formula for the change of the surface deflection angle due to a perturbation (in the upper layer) of a given eddy viscosity profile for which the solution to the problem is known. This extends and generalizes previous results for perturbations of a constant eddy viscosity. The perturbation analysis is performed using a recent reformulation of the linear second-order governing equations for Ekman-type flows as a nonlinear first-order Riccati-type equation.
Similar content being viewed by others
1 Introduction
Wind-driven currents are one of the main aspects of oceanic circulation, which is largely responsible for the heat transport between different geographical regions and has a crucial impact on the Earth’s climate [15]. At non-equatorial latitudes, the flow below a depth of about 100 metres is approximately in geostrophic balance, that is, the force balance occurs between the Coriolis force and horizontal pressure gradients. The frictional effects of surface winds are confined to a relatively thin boundary layer, called the Ekman layer, where the momentum of the wind is carried down to the interior of the water flow by the wind-induced turbulence of the water. Thus, in this surface layer we find currents (superimposed to the underlying geostrophic motion) where the force balance is between the Coriolis force and the frictional forces due to the wind (see [15, 18]). It is observed in practice that such wind-induced currents do not follow the direction of the generating wind, but instead are deflected to the right of the wind direction in the Northern Hemisphere (to the left in the Southern Hemisphere), the amplitude of the observed deflection angle being generally within the range 10\(^\circ \)–75\(^\circ \) (cf. the data in [17] or [20]).
On the other hand, the ocean dynamics near the Equator differs significantly from that at higher latitudes. In fact, since the Coriolis force vanishes along the Equator, geostrophic effects are negligible there and nonlinear effects arise at leading order; moreover, the Coriolis force changes sign across the Equator and so acts as a waveguide that gives rise to azimuthal (i.e. longitudinal) flow propagation (see e.g. [3, 4, 6, 12]). Therefore, substantially different approaches are required for the description of non-equatorial and equatorial flows. This paper is devoted to the former aspect of oceanic wind-driven currents; for the latter case, we refer the reader to the references above.
The Ekman layer is named after the Swedish scientist who, in 1905, was the first to write down and investigate a model that should explain the observations of the surface deflection angle that had been made at that time. In his classical work [10], the eddy viscosity (a quantity introduced as a measure of the frictional effects in the boundary layer—see the next section) is set as a constant, yielding a solution that, at the surface, is directed at an angle of 45\(^\circ \) to the right (left) of the wind direction in the Northern (Southern) Hemisphere. The source of the inaccuracy of the predicted surface deflection angle has been identified in the simplified assumption of constant eddy viscosity, and over the last few decades many attempts have been made at finding solutions for more general depth-dependent eddy viscosity profiles; many of the known explicit solutions for particular cases can be found for instance in [5, 7, 8, 11, 14, 16]. Nevertheless, finding explicit solutions can be a very hard task even for rather simple eddy viscosity profiles, so alternative ways to tackle the problem have been sought. In [19], an application of the WKB approach yields an asymptotic approximation of the solution to the general Ekman model; in [1, 2], a formula for the deviation of the deflection angle from the reference value of 45\(^\circ \) for a given small perturbation of a constant eddy viscosity is derived; in [9], the problem of a two-layer ocean with constant eddy viscosity in the lower layer and a general continuous eddy viscosity in the upper layer is transferred to solving a Riccati-type equation in the upper layer, from whose solution the surface deflection angle for the original problem can be determined straightforwardly. In the present article we combine the latter two methods: using the reformulation of the problem as a Riccati-type equation, we perform a perturbation analysis to derive a formula for the change of the surface deflection angle caused by a perturbation of a general continuous eddy viscosity profile in the upper layer.
2 The Governing Equations
The steady-state equations for non-equatorial Ekman-type flows in the f-plane approximation for geophysical flows (cf. [15, 18]) read as follows (see [2, 8]):
where (u, v) are the horizontal Cartesian components of the flow in the f-plane, \(Z\in (-\infty ,0)\) is the depth below the mean surface level (located at \(Z=0\)), \(\nu \) is the (depth-dependent) eddy viscosity, and \(f = 2\Omega \sin (\theta )\) is the Coriolis parameter at the fixed latitude \(\theta \in (-\pi /2,0)\cup (0,\pi /2)\), with the angular speed of rotation of the Earth \(\Omega = 7.29\times 10^5\,\)rad/s; in the Northern Hemisphere we have \(f>0\), whereas \(f<0\) in the Southern Hemisphere. In the following we choose to work in the Northern Hemisphere; the results can be readily adapted to the other case. Introducing the complex notation \(U = u+\mathrm {i}v\), the equations above can be rewritten as
The surface boundary condition is
where \(\rho \) is the (constant) density of the water and \(\tau _1+\mathrm {i}\tau _2\) is the surface wind stress, which we assume to be constant as well. In fact, given the fact that in the f-plane approximation we are dealing with a planar geometry and constant Coriolis parameter f, we can rotate the horizontal coordinate system so that the real axis becomes aligned with the wind direction: this means that, without loss of generality, we may assume \(\tau _2=0\) and \(\tau _1>0\). Moreover, we require
The problem (2.1)–(2.3) can be non-dimensionalized by setting
From this we obtain the non-dimensional problem
where the prime denotes a derivative with respect to z. Let us remark that the surface deflection angle, which is the main objective of this paper, is the argument of the complex vector \(\psi (0)\).
Now, let us suppose that, for the continuous \(K_0:(-\infty ,0]\rightarrow (0,\infty )\) with \(K_0(z)=1\) for \(z\le -h\) (for some non-dimensional depth \(-h<0\)), the solution \(\psi _0\) to
is known. Henceforth we consider a K of the form
for small \(\varepsilon >0\), uniformly in z: explicitly, this means that for each \(C>0\) there exists \(\varepsilon _0\) such that for all \(\varepsilon \le \varepsilon _0\) and all \(z\le 0\) we have
Moreover, we assume that \(K_1\) is continuous on \((-\infty ,0]\) and
Notice that the requirement of a constant K below some depth is physically reasonable (cf. [9]): in fact, the turbulence is confined to the thin Ekman boundary layer, thus the eddy viscosity at greater depth coincides with that of sea water, which can be normalized to obtain
3 Perturbation Analysis
In [9] it is shown that, if q solves the Riccati equation
then for the restriction of the solution \(\psi \) of (2.4) on \((-h,0)\) we have that \(\psi (0)=K(0)/q(0)\),
and
Following (2.5), let us make the ansatz
(uniformly in z) for the solution to (3.1), where \(q_0\) solves
and \(q_1\) is to be determined. Inserting (2.5) and (3.4) into (3.1a), using (3.5a), dividing by \(\varepsilon \), and letting \(\varepsilon \rightarrow 0\), we get the following equation for \(q_1\):
with “initial” condition
in view of (3.1b) and (3.5b). The solution to (3.6)–(3.7) is
Thus
The formula
(cf. [1, 2]) enables us to express the change of \(\arg (q(0))\) due to the perturbation as
Therefore, by (3.2), a positive/negative value of \(C_{K_1}\) corresponds to a decrease/increase of \(\arg (\psi (0))\) from the reference value \(\arg (\psi _0(0))\). Moreover, (3.3) implies
thus formula (3.8) can be rewritten in terms of \(\psi _0\) as
4 Discussion
The easiest situation where formula (3.9) can be applied is the case of constant \(K_0 \equiv 1\) on \((-\infty ,0]\) (cf. the discussion in [1, 2]). Here we have
hence, according to (3.2) and (3.9), the change of \(\arg (\psi (0))\) is given by
which corresponds, up to a scaling, to the formula derived in [2] applied to an upper-layer perturbation which is confined to the interval \((-h,0)\).
It is interesting to notice that the presence of a sine function in the integrand is not just a feature of the special case of constant \(K_0\). In fact, if we write
then formula (3.9) becomes
From this it is apparent that in general it is very difficult to tell at first glance whether a given perturbation \(K_1\) of \(K_0\) will lead to a positive or negative increment of the surface deflection angle. For instance, suppose that \(K_1(z)>0\) for each \(z\in (-h,0)\). Due to the periodicity of the sine, we may assume that \(2\arg (\psi _0'(s))-\arg (\psi _0(0)) \in (-\pi ,\pi ]\) for each \(s\in (-h,0)\). Let us denote
Since \(\sin [2\arg (\psi _0'(s))-\arg (\psi _0(0))]\) is positive for \(s\in P\) and negative for \(s\in N\), the sign of \(C_{K_1}\) will depend both on the measures \(\vert P\vert \) and \(\vert N\vert \) and on the relative sizes of the values attained by \(K_1\) on P and N.
Finally, let us look at one last example that shows how the “Riccati” variant (3.8) of (3.9) can be useful in its own way, not just as a step in the derivation of (3.9). This is especially true if the solution is found using the reformulation (3.1) in the first place. In [9] the solution to the Riccati formulation of the problem for the case
is determined to be
where \(S(z)=[3(z+h)+1]^{1/3}\) and
Using the identity
and the fundamental theorem of calculus, a straightforward computation yields
which can be directly plugged into (3.8) without first having to determine and differentiate \(\psi _0\).
We conclude by saying that formula (3.9) may be of some relevance especially in view of the fact that, as pointed out in the introduction, analytic solutions are hard to find and the few ones available show an intricate behaviour even for relatively simple eddy viscosity profiles. The perturbation approach that we have pursued here may be useful to bring some new insight into this delicate problem, but it would be interesting to obtain further results that do not require the smallness assumption that is intrinsic to the perturbation method. On the other hand, it would also be of great importance to bring time dependence into the picture; our current knowledge of time-dependent Ekman flows is very limited and is restricted to some special cases (see [13]), but it seems likely that the technique used in this paper may be employed to obtain further results in this direction.
References
Bressan, A., Constantin, A.: The deflection angle of surface ocean currents from the wind direction. J. Geophys. Res. Oceans 124, 7412–7420 (2019)
Constantin, A.: Frictional effects in wind-driven ocean currents. Geophys. Astro. Fluid Dyn. 115(1), 1–14 (2021)
Constantin, A., Ivanov, R.I.: Equatorial wave-current interactions. Comm. Math. Phys. 370, 1–48 (2019)
Constantin, A., Johnson, R.S.: The dynamics of waves interacting with the Equatorial Undercurrent. Geophys. Astrophys. Fluid Dyn. 109, 311–358 (2015)
Constantin, A., Johnson, R.S.: Ekman-type solutions for shallow-water flows on a rotating sphere: a new perspective on a classical problem. Phys. Fluids 31, 021401 (2019)
Constantin, A., Johnson, R.S.: Steady large-scale ocean flows in spherical coordinates. Oceanography 31, 42–50 (2018)
Constantin, A., Johnson, R.S.: Atmospheric Ekman flows with variable eddy viscosity. Bound. Layer Meteorol. 170, 395–414 (2019)
Dritschel, D.G., Paldor, N., Constantin, A.: The Ekman spiral for piecewise-uniform viscosity. Ocean Sci. 16, 1089–1093 (2020)
Dritschel, D.G., Paldor, N., Constantin, A.: The deflection angle between a wind-forced surface current and the overlying wind in an ocean with vertically varying eddy viscosity. Phys. Fluids 32, 116604 (2020)
Ekman, V.W.: On the influence of the Earth’s rotation on ocean-currents. Ark. Mat. Astron. Fys. 2, 1–52 (1905)
Grisogono, B.: A generalized Ekman layer profile with gradually varying eddy diffusivities. Quart. J. Roy. Meteorol. Soc. 121, 445–453 (1995)
Ionescu-Kruse, D., Martin, C.I.: Periodic equatorial water flows from a Hamiltonian perspective. J. Differ. Equ. 262, 4451–4474 (2017)
Lewis, D.M., Belcher, S.E.: Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift. Dyn. Atmos. Oceans 37, 313–351 (2004)
Madsen, O.S.: A realistic model of the wind-induced Ekman boundary layer. J. Phys. Oceanogr. 7, 248–255 (1977)
Marshall, J., Plumb, R.A.: Atmosphere, Ocean and Climate Dynamics: An Introductory Text. Academic Press, New York (2016)
Roberti, L.: The Ekman spiral for piecewise-constant eddy viscosity. Appl. Anal. (2021). https://doi.org/10.1080/00036811.2021.1896709
Röhrs, J., Christensen, K.H.: Drift in the uppermost part of the ocean. Geophys. Res. Lett. 42, 10349–10356 (2015)
Vallis, G.K.: Atmospheric and Oceanic Fluid Dynamics: Fundamentals and Large-Scale Circulation, 2nd edn. Cambridge University Press, Cambridge (2017)
Wenegrat, J.O., McPhaden, M.J.: Wind, waves, and fronts: frictional effects in a generalized Ekman model. J. Phys. Oceanogr. 46, 371–394 (2016)
Yoshikawa, Y., Masuda, A.: Seasonal variations in the speed factor and deflection angle of the wind driven surface flow in the Tsushima Strait. J. Geophys. Res. 114, C12022 (2009)
Acknowledgements
The author is grateful to the reviewer for the valuable comments and suggestions that helped to improve the paper.
Funding
Open access funding provided by University of Vienna. The financial support of the Vienna Science and Technology Fund (WWTF) under grant MA16-009 is acknowledged.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The author states that there is no conflict of interest.
Additional information
Communicated by A. Constantin.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Roberti, L. Perturbation Analysis for the Surface Deflection Angle of Ekman-Type Flows with Variable Eddy Viscosity. J. Math. Fluid Mech. 23, 57 (2021). https://doi.org/10.1007/s00021-021-00586-y
Accepted:
Published:
DOI: https://doi.org/10.1007/s00021-021-00586-y