A linear model of separation for western boundary currents with bathymetry
Anne-Laure Dalibard, Corentin Gentil
TL;DR
The paper develops a linear, highly structured model for rotating, stratified flows on a $eta$-plane with topography to capture western boundary current separation. It derives a 3D quasi-geostrophic core coupled to two boundary-layer types: a quasi-geostrophic Munk layer and a non-geostrophic Ekman layer, and constructs an approximate solution to arbitrary order in an $ ablaeta$-scaled expansion, with rigorous justification of validity. A key novelty is the detailed treatment of Ekman layers in the presence of strong topography and stratification, yielding new dynamical features and a mechanism for boundary-current detachment. The work provides both periodic and Cauchy-problem stability results and supports the theoretical construction with numerical illustrations of separation, suggesting a pathway to parameterizations of topography-stratification interactions in coarse ocean models.
Abstract
This paper is devoted to the asymptotic analysis of strongly rotating and stratified fluids, under a $β$-plane approximation, and within a three-dimensional spatial domain with strong topography. Our purpose is to propose a linear idealized model, which is able to capture one of the key features of western boundary currents, in spite of its simplicity: the separation of the currents from the coast. Our simplified framework allows us to perform explicit computations, and to highlight the intricate links between rotation, stratification and bathymetry. In fact, we are able to construct approximate solutions at any order for our system, and to justify their validity. Each term in the asymptotic expansion is the sum of an interior part and of two boundary layer parts: a ``Munk'' type boundary layer, which is quasi-geostrophic, and an ``Ekman part'', which is not. Even though the Munk part of the approximation bears some similarity with previously studied 2D models, the analysis of the Ekman part is completely new, and several of its properties differ strongly from the ones of classical Ekman layers. Our theoretical analysis is supplemented with numerical illustrations, which exhibit the desired separation behavior.
