![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
ABSTRACT
Using the Navier–Stokes equations in rotating spherical coordinates, we discuss a mathematical model for the equatorial current across the Pacific Ocean. Symmetry properties of the governing equations that hold in equatorial regions permit us to gain detailed insight into the structure of the ocean flow dynamics.
2000 MATHEMATICS SUBJECT CLASSIFICATIONS:
1. Introduction
The equatorial ocean flows present special challenges with respect to flows at mid-latitudes since in this context neither geostrophic theory nor Ekman theory apply (see the discussions in [Citation1,Citation2] and [Citation3], respectively) and, moreover, nonlinear effects are noticeable (see the discussion in [Citation4,Citation5]). Of all equatorial regions, that of the Pacific is the most important due to its extent – it is more than 12,000 km across the ocean, from South America to South Asia. Several recent papers are devoted to equatorial flows (see the discussions in [Citation33–Citation36]), and address the issue of underlying currents, issue that was mostly ignored in the earlier research literature (see the survey [Citation6]). The papers [Citation7–11] find explicit wave solutions for the governing equations in Lagrangian coordinates in the β-plane approximation, in the presence of relatively weak underlying currents. Two-dimensional flows (which ignore the meridional velocity component and use the f-plane approximation) are investigated in [Citation1,Citation12] (see also [Citation13]) and inroads towards three-dimensional flows can be found in [Citation14–17] (see [Citation18] for a detailed description of the intricate ongoing dynamics). One issue of great current interest is the derivation of realistic models for the underlying current field. The main features of the equatorial current are that the current is mainly azimuthal and presents a strong depth-variation in the upper 500 m of the ocean (of about 4 km mean depth), while at depths larger than 500 m there is practically no motion. The most significant features of the depth-variation consist of a westward near-surface current (induced by the trade winds that blow towards the west) and a stronger eastward flow – the Equatorial Undercurrent (EUC) – that dominates the overall ocean flow, reaching speeds in excess of 1 m/s. Note the unusual feature that the EUC flows against the prevailing trade winds, as the directly wind-driven flow vanishes below 100 m. Recently, an approach was proposed to explain how such a current with flow-reversal might be induced by the action of the wind (see [Citation19]). However, the model in [Citation19] does not capture the westward Equatorial Intermediate Current (EIC), found directly below the EUC across the Pacific (see [Citation20]). Our aim in this paper is to show that the method used in [Citation19], to investigate the Navier–Stokes equations in rotating spherical coordinates by taking advantage of symmetry properties of the leading-order terms in the governing equations for equatorial ocean flows to gain insight into the structure of the velocity field, can be further refined to provide a model of higher accuracy. In particular, we are able to include the westward EIC as a component of the background state of the ocean.
2. Preliminaries
For large-scale ocean flows it is appropriate to consider a spherical Earth, rotating at a constant rate about its polar axis (see [Citation21,Citation22]). We introduce near the Equator spherical coordinates , where
is the radius,
is the angle of latitude and
is the angle of longitude. The corresponding unit vectors are
, respectively, with
pointing from West to East,
from South to North and
upwards. We denote the ocean flow velocity components by
. We use primes for physical (dimensional) variables, and we will remove them in the non-dimensionalisation process.
In this rotating reference system the Navier–Stokes equation and the equation of mass conservation are (see [Citation19])
(1)
(1) and
(2)
(2) respectively, where
is the pressure in the fluid,
the constant density,
is the constant rate of rotation of the Earth, and
is the gravitational acceleration. The coefficients of the viscous terms are taken to be constant, with
the vertical kinematic eddy viscosity and
the horizontal kinematic eddy viscosity.
Let us now describe the associated boundary conditions. At the ocean's surface, , where
is the radius of the spherical Earth, we impose the dynamic and the kinematic boundary condition:
(3)
(3) and
(4)
(4) respectively, where
is the surface pressure. The wind-stress at the ocean's surface is considered to be known, and can be expressed in the form
(5)
(5) the surface wind stress
being related to the vertical eddy viscosity by
on the surface, where
is a (dimensional) constant and
,
being the density of air (about 1.2 kg/m
) and
being a (dimensionless) drag coefficient. At the impermeable, solid, stationary bottom of the ocean,
, we have the corresponding boundary condition for viscous flow:
(6)
(6) Setting
and
(7)
(7) one can non-dimensionalise the problem by performing the change of variables
(8)
(8) where
is the average depth of the near-surface layer to which the wind effects are confined),
is the typical speed of mid-latitude ocean currents at the surface, while the scaling factor k for the vertical velocity is less than 10
; see the discussion in [Citation19]. In terms of the the shallow-water parameter
with a typical value of the order
, the governing equations for steady flow are transformed into the non-dimensional form
(9)
(9) and
(10)
(10) where
and
(j = 1, 2) are the inverse Rossby number and the pair of Reynolds numbers, respectively. With
, the non-dimensional boundary conditions are
(11)
(11)
(12)
(12)
(13)
(13)
For ocean flows in the equatorial Pacific it is adequate (see the discussion in [Citation19]) to consider
Multiplying the third component of (Equation9
(9)
(9) ) throughout by ϵ, we obtain the non-dimensional Navier–Stokes system with Coriolis effects
(14)
(14) Correspondingly, the equation of mass conservation, (Equation10
(10)
(10) ), becomes
(15)
(15) In the non-dimensional system of equations (Equation14
(14)
(14) )–(Equation15
(15)
(15) ), the parameters
are held fixed, as no approximations were performed.
The shallow-water limit corresponds to the limiting process and
(see the discussion in [Citation19]). This regime ignores the vertical velocity component and, neglecting wave perturbations by setting h = 0, the flow dynamics being governed by the horizontal flow components u and v, subject to the nonlinear system
(16)
(16)
(17)
(17)
(18)
(18) which features Coriolis terms, viscous terms and the horizontal pressure gradients, while the boundary conditions (Equation11
(11)
(11) )–(Equation13
(13)
(13) ) simplify to
(19)
(19)
(20)
(20) since setting w = 0 and h = 0 makes (Equation12
(12)
(12) ) irrelevant. The system (Equation16
(16)
(16) )–(Equation20
(20)
(20) ) captures the physical idea that the perturbation
of the hydrostatic pressure and the shear stresses
and
, at the ocean's surface z = 0, are needed to produce a consistent solution for the wind-drift horizontal current
. To gain insight into the structure of this system, let us note that (Equation18
(18)
(18) ) ensures (see the discussion in [Citation23]) the existence of a stream function in spherical coordinates,
, with
(21)
(21) see [Citation23] for a proof of the fact that the particle paths for steady flow on the surface of a sphere are the level sets ψ=constant. Taking advantage of (Equation21
(21)
(21) ), we can eliminate the pressure between Equations (Equation16
(16)
(16) ) and (Equation17
(17)
(17) ) to derive (see [Citation19]) the vorticity equation
(22)
(22) These considerations show that (Equation21
(21)
(21) ) specifies, at leading order, the background flow, provided that the stream function
solves the vorticity equation (Equation22
(22)
(22) ) subject to the boundary condition (Equation19
(19)
(19) ) and to the last two constraints in (Equation20
(20)
(20) ); Equations (Equation16
(16)
(16) )–(Equation17
(17)
(17) ) then determine the associated pressure field, taking also into account the boundary condition represented by the first constraint in (Equation19
(19)
(19) ). Let us note that if we ignore the z-dependence (and thus, implicitly, consider an inviscid setting in which the wind forcing plays no role), then (Equation22
(22)
(22) ) simplifies to the ocean gyre model derived recently in [Citation24] as a shallow-water asymptotic solution of Euler's equation in rotating spherical coordinates (with the stipulation that θ stands in [Citation24] for the polar angle, and not for the angle of latitude) and further investigated in [Citation25–29] in the context of the Antarctic Circumpolar Current – the largest ocean current on Earth, flowing clockwise from west to east around Antarctica (see [Citation30,Citation32,Citation37]) so that, due to the lack of any landmass connecting with Antarctica, it keeps the warm ocean waters from lower latitudes away from Antarctica and thus maintains the huge ice sheets encountered near the South Pole (see the discussion in [Citation13,Citation31]).
3. Main results
Finding the general solution of the equation (Equation22(22)
(22) ) with the boundary conditions
(23)
(23)
(24)
(24) where the horizontal velocity field
is specified by means of (Equation21
(21)
(21) ), is not a realistic task. Instead, we try to take advantage of some special features that are typical for the equatorial region in the Pacific. We restrict our attention to an equatorial zonal band that is symmetric about the Equator, with the angle of latitude θ close to 0 and the angle of longitude ϕ near
, with
corresponding to the region between
and
. For the twice differentiable functions
and
, we look for longitude-independent meridional velocity components and we seek a linear dependence of the azimuthal velocity component on the longitude. These features lead (see [Citation19]) to the stream function
(25)
(25) solution of the governing equation (Equation22
(22)
(22) ), valid for
and
with some fixed (and small) values
and
. Due to (Equation21
(21)
(21) ), this corresponds to the horizontal velocity field
(26)
(26) while the boundary conditions (Equation23
(23)
(23) ) are induced by a wind stress with components
(27)
(27) Rather than imposing (Equation24
(24)
(24) ), we pursue the considerations in [Citation4]) and interpret the solution (Equation26
(26)
(26) ) as being valid in the near-surface ocean region above the thermocline z = −T, where we require the no-stress boundary condition
(28)
(28)
The fundamental characteristic of the trade winds in the equatorial Pacific is that, in each hemisphere, they are oriented towards the Equator, blowing westwards with a more stronger westward direction as the Equator () is approached. This property holds for the wind stress specified in (Equation27
(27)
(27) ) if
(29)
(29) for all relevant values of the longitude ϕ, that is, for
. Note that this orientation of the wind, in combination with the fact that the vanishing of the meridional component of the Coriolis force,
, at the Equator prevents it from inducing a deflection from the wind direction (as is typical in Ekman theory at mid-latitudes), leads to a near-surface current that moving westward. This is ensured for the solution (Equation26
(26)
(26) ) if
(30)
(30) Moreover, if
(31)
(31) then the meridional flow described by (Equation26
(26)
(26) ) is poleward near the surface. Note that the main features of the near-surface flow in the equatorial Pacific are a westward motion at the surface and a poleward meridional flow (see the discussion in [Citation4]). One can see that a linear or quadratic dependence on the z-variable of the azimuthal velocity profile
of type (Equation26
(26)
(26) ) can not accommodate the constraints (Equation29
(29)
(29) )–(Equation31
(31)
(31) ) and (Equation28
(28)
(28) ). A solution featuring a cubic dependence was provided recently in [Citation19], namely
(32)
(32) for some constant a<0. This provides an azimuthal flow that is westward near the surface and features a faster eastward jet along the thermocline, such that at any fixed longitude
, the zonal velocity u strictly decreases from a positive value at the thermocline z = −T (corresponding to an eastward flow) to a negative value at the surface z = 0 (corresponding to a westward flow), and vanishes once above
, having just one inflexion point, to be found between the depth levels
and
. The solution (Equation32
(32)
(32) ) captures the main features of the ocean current in the equatorial Pacific: a westward wind-drift overlying the stronger eastward Equatorial Undercurrent (EUC). We now show that we can refine this result, by accommodating also a weaker eastward flow beneath the EUC (see Figure ). For this, it suffices to choose suitable quintic polynomial expressions for the azimuthal velocity profile.
Figure 1. Depiction of the vertical profile of the current in the upper 200–400 m of the Pacific Ocean along the Equator: the westward wind-drift current is near the surface z = 0, below it is the eastward EUC that dominates the subsurface flows, with the weaker westward EIC found directly below the EUC, while at great depths there is practically no motion. These are the main features across the Pacific, over 12,000 km (see the data provided in [Citation20]).
![Figure 1. Depiction of the vertical profile of the current in the upper 200–400 m of the Pacific Ocean along the Equator: the westward wind-drift current is near the surface z = 0, below it is the eastward EUC that dominates the subsurface flows, with the weaker westward EIC found directly below the EUC, while at great depths there is practically no motion. These are the main features across the Pacific, over 12,000 km (see the data provided in [Citation20]).](/cms/asset/1cf5db79-42aa-4087-a435-20b19a56bb51/gapa_a_1673375_f0001_oc.jpg)
Theorem 3.1
The quintic polynomial expressions
accommodate an equatorial current profile that features an eastward jet at mid-depth of the near-surface layer above the thermocline
a weak eastward jet just above the thermocline and a stronger westward jet near the surface z = 0.
Proof.
The corresponding azimuthal velocity , obtained from (Equation26
(26)
(26) ), vanishes at
,
and
, giving the profile depicted in the figure. Moreover, since
all the constraints (Equation29
(29)
(29) )–(Equation31
(31)
(31) ) and (Equation28
(28)
(28) ) are verified. The considerations preceding the statement yield now the desired result.
Remark 3.2
The theorem shows that it is possible to capture within our modeling framework not only the presence of the eastward EUC and of the westward wind-drift current, but also the presence of the weaker westward EIC, found directly below the EUC (see the figure).
Disclosure statement
No potential conflict of interest was reported by the author.
Additional information
Funding
References
- Constantin A, Ivanov RI. Equatorial wave-current interactions. Comm Math Phys. 2019;370:1–48. doi: 10.1007/s00220-019-03483-8
- Cronin MF, Kessler WS. Near-surface shear flow in the tropical Pacific Cold Tongue Front. J Phys Oceanogr. 2009;39:1200–1215. doi: 10.1175/2008JPO4064.1
- Constantin A, Johnson RS. Large-scale oceanic currents as shallow-water asymptotic solutions of the Navier-Stokes equation in rotating spherical coordinates. Deep Sea Res Part II: Top Stud Oceanogr. 2019;160:32–40. doi: 10.1016/j.dsr2.2018.12.007
- Boyd JP. Dynamics of the equatorial ocean. Berlin: Springer; 2018.
- Henry D. Nonlinear features of equatorial ocean flows. Oceanography. 2018;31:22–27. doi: 10.5670/oceanog.2018.305
- Fedorov AV, Brown JN. Equatorial waves. In: Steele J, editor. Encyclopedia of ocean sciences. Academic Press; 2009, p. 3679–3695.
- Constantin A. An exact solution for equatorially trapped waves. J Geophys Res: Oceans. 2012;117:C05029. doi: 10.1029/2012JC007879
- Constantin A. Some nonlinear, equatorially trapped, nonhydrostatic internal geophysical waves. J Phys Oceanogr. 2014;44:781–789. doi: 10.1175/JPO-D-13-0174.1
- Constantin A, Germain P. Instability of some equatorially trapped waves. J Geophys Res: Oceans. 2013;118:2802–2810. doi: 10.1002/jgrc.20219
- Henry D. An exact solution for equatorial geophysical water waves with an underlying current. Eur J Mech B/Fluids. 2013;38:18–21. doi: 10.1016/j.euromechflu.2012.10.001
- Henry D. Equatorially trapped nonlinear water waves in a β-plane approximation with centripetal forces. J Fluid Mech. 2016;804:R1. doi: 10.1017/jfm.2016.544
- Constantin A, Johnson RS. The dynamics of waves interacting with the Equatorial Undercurrent. Geophys Astrophys Fluid Dyn. 2015;109:311–358. doi: 10.1080/03091929.2015.1066785
- Henry D, Martin CI. Exact, free-surface equatorial flows with general stratification in spherical coordinates. Arch Ration Mech Anal. 2019;233:497–512. doi: 10.1007/s00205-019-01362-z
- Basu B. On a three-dimensional nonlinear model of Pacific equatorial ocean dynamics: velocities and flow paths. Oceanography. 2018;31:51–58. doi: 10.5670/oceanog.2018.309
- Basu B. Some numerical investigations into a nonlinear three-dimensional model of the Pacific equatorial ocean flows. Deep Sea Res Part II: Top Stud Oceanogr. 2019;160:7–15. doi: 10.1016/j.dsr2.2018.11.013
- Constantin A, Johnson RS. An exact, steady, purely azimuthal equatorial flow with a free surface. J Phys Oceanogr. 2016;46:1935–1945. doi: 10.1175/JPO-D-15-0205.1
- Constantin A, Johnson RS. A nonlinear, three-dimensional model for ocean flows, motivated by some observations of the Pacific Equatorial Undercurrent and thermocline. Phys Fluids. 2017;29:056604. doi: 10.1063/1.4984001
- Constantin A, Johnson RS. On the nonlinear, three-dimensional structure of equatorial oceanic flows. J Phys Oceanogr. 2019;49:2029–2042. doi: 10.1175/JPO-D-19-0079.1
- Constantin A, Johnson RS. Ekman-type solutions for shallow-water flows on a rotating sphere: a new perspective on a classical problem. Phys Fluids. 2019;31:021401. doi: 10.1063/1.5083088
- Lukas R. Pacific ocean equatorial currents. In: Thorpe SA, Turekian KK, editors. Encyclopedia of ocean sciences. Elsevier; 2001, p. 2069–2076.
- Marshall J, Plumb RA. Atmosphere, ocean and climate dynamics: an introductory text. London: Academic Press; 2016.
- Vallis GK. Atmosphere and ocean fluid dynamics. Cambridge: Cambridge University Press; 2006.
- Constantin A, Johnson RS. Steady large-scale ocean flows in spherical coordinates. Oceanography. 2018;31:42–50. doi: 10.5670/oceanog.2018.308
- Constantin A, Johnson RS. Large gyres as a shallow-water asymptotic solution of Euler's equation in spherical coordinates. Proc Roy Soc London A. 2017;473, Art. No. 20170063. doi: 10.1098/rspa.2017.0063
- Haziot S. Study of an elliptic partial differential equation modelling the Antarctic Circumpolar Current. Discr Cont Dyn Syst. 2019;39:4415–4427. doi: 10.3934/dcds.2019179
- Marynets K. Two-point boundary problem for modeling the jet flow of the Antarctic Circumpolar Current. Electron J Diff Eq. 2018;56:1–12. Paper No. 56.
- Marynets K. A nonlinear two-point boundary-value problem in geophysics. Monatsh Math. 2019;188:287–295. doi: 10.1007/s00605-017-1127-x
- Marynets K. On a two-point boundary-value problem in geophysics. Appl Anal. 2019;98:553–560. doi: 10.1080/00036811.2017.1395869
- Marynets K. On the modeling of the flow of the Antarctic Circumpolar Current. Monatsh Math. 2019;188:561–565. doi: 10.1007/s00605-017-1147-6
- Constantin A, Johnson RS. An exact, steady, purely azimuthal flow as a model for the Antarctic Circumpolar Current. J Phys Oceanogr. 2016;46:3585–3594. doi: 10.1175/JPO-D-16-0121.1
- Walton DWHAntarctica: global science from a frozen continent. Cambridge: Cambridge University Press; 2013.
- Haziot S, Marynets K. Applying the stereographic projection to modeling of the flow of the Antarctic Circumpolar Current. Oceanography. 2018;31:68–75. doi: 10.5670/oceanog.2018.311
- Ionescu-Kruse D. A three-dimensional autonomous nonlinear dynamical system modelling equatorial ocean flows. J Diff Eqs. 2018;264:4650–4668. doi: 10.1016/j.jde.2017.12.021
- Ionescu-Kruse D, Martin CI. Local stability for an exact steady purely azimuthal equatorial flow. J Math Fluid Mech. 2018;20:27–34. doi: 10.1007/s00021-016-0311-4
- Izumo T. The equatorial current, meridional overturning circulation, and their roles in mass and heat exchanges during the El Nino events in the tropical Pacific Ocean. Ocean Dyn. 2005;55:110–123. doi: 10.1007/s10236-005-0115-1
- Martin CI. On the existence of free-surface azimuthal equatorial flows. Appl Anal. 2017;96:1207–1214. doi: 10.1080/00036811.2016.1180370
- Marynets K. The Antarctic Circumpolar Current as a shallow-water asymptotic solution of Euler's equation in spherical coordinates. Deep Sea Res Part II: Top Stud Oceanogr. 2019;160:58–62. doi: 10.1016/j.dsr2.2018.11.014