On large-scale wind-drift ocean currents: An asymptotic approach in spherical coordinates

Starting from the Navier–Stokes equations in rotating spherical coordinates with depth-varying density and eddy viscosity, we derive an asymptotic model describing non-equatorial wind-generated ocean drift currents. Our approach allows for large-scale flows that cannot be captured by classical tangent-plane approximations. The strategy is to perform a careful scaling and to perform a double asymptotic expansion with respect to two small parameters arising from the scaling: the Rossby number and the ratio between the Ekman depth and the Earth’s radius. We obtain a system of linear ordinary differential equations with nonlinear boundary conditions governing the leading-order dynamics, highlighting that the dynamics is governed by the linear terms, whereas the nonlinear ones, related to the injection and dissipation of kinetic energy, appear only at higher order. We use the leading-order equations to compare our model with the simplest theory of ocean circulation due to Sverdrup and note that, even at this level of simplification, our equations have the potential to provide deeper insight. Subsequently, focusing on Ekman flows, we prove existence and uniqueness of the leading-order solution, which retains the classical Ekman spiral structure for arbitrary eddy-viscosity profiles. Finally, we compute the surface deflection angle of the wind-driven current for three explicit eddy-viscosity profiles, obtaining results consistent with observations. In addition, we derive the governing equations for the first-order correction with respect to the Rossby number and provide a priori bounds for their solution.