Sliding or Stumbling on the Staircase: Numerics of Ocean Circulation Along Piecewise‐Constant Coastlines

Coastlines in most ocean general circulation models are piecewise constant. Accurate representation of boundary currents along staircase‐like coastlines is a long‐standing issue in ocean modeling. Pioneering work by Adcroft and Marshall (1998, https://doi.org/10.3402/tellusa.v50i1.14514) revealed th...

Full description

Saved in:
Bibliographic Details
Published in:Journal of advances in modeling earth systems Vol. 15; no. 5
Main Authors: Nasser, Antoine‐Alexis, Madec, Gurvan, Lavergne, Casimir, Debreu, Laurent, Lemarié, Florian, Blayo, Eric
Format: Journal Article
Language:English
Published: Washington John Wiley & Sons, Inc 01-05-2023
American Geophysical Union
American Geophysical Union (AGU)
Subjects:
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Coastlines in most ocean general circulation models are piecewise constant. Accurate representation of boundary currents along staircase‐like coastlines is a long‐standing issue in ocean modeling. Pioneering work by Adcroft and Marshall (1998, https://doi.org/10.3402/tellusa.v50i1.14514) revealed that artificial indentation of model coastlines, obtained by rotating the numerical mesh within an idealized square basin, generates a spurious form drag that slows down the circulation. Here, we revisit this problem and show how this spurious drag may be eliminated. First, we find that physical convergence to spatial resolution (i.e., the main characteristics of the flow are insensitive to the increase of the mesh resolution) allows simulations to become independent of the mesh orientation. An advection scheme with a wider stencil also reduces sensitivity to mesh orientation from coarser resolution. Second, we show that indented coastlines behave as straight and slippery shores when a true mirror boundary condition on the flow is imposed. This finding applies to both symmetric and rotational‐divergence formulations of the stress tensor, and to both flux and vector‐invariant forms of the equations. Finally, we demonstrate that the detachment of a vortex flowing past an outgoing corner of the coastline is missed with a free‐slip (zero vorticity) condition at the corner. These results provide guidance for a better numerical treatment of coastlines (and isobaths) in ocean general circulation models. Plain Language Summary Most ocean general circulation models represent coastlines as piecewise constant, which does not accurately reflect the true boundary. This approximation is necessitated by the size and square shape of model grid cells, together with finite computational resources, making it difficult to finely represent the boundary. A long‐standing issue in ocean modeling is accurately representing boundary currents along these staircase‐like coastlines. In particular, Adcroft and Marshall (1998, https://doi.org/10.3402/tellusa.v50i1.14514) discovered that artificial indentation of the model coastlines generates a “spurious form drag” that slows down the circulation. Our study revisits this long‐standing issue and shows how this spurious drag may be eliminated. First, we demonstrate that having a sufficiently fine spatial resolution to resolve the physical processes allows the model to be insensitive to coastal indentation. We also show that indented coastlines become slippery, as if they were smooth and straight shores, when a true mirror boundary condition on the flow is imposed. Finally, we show how to faithfully simulate the retroflection of a current past a cape. In summary, these results provide guidance for a better numerical representation of marine land‐forms in numerical ocean models. Key Points The impact of coastal indentation on circulation within an idealized square basin is studied systematically for various numerical choices Solutions are insensitive to coastal indentation with rot‐div formulation of the diffusive tensor and sufficiently fine spatial resolution Indented coastlines become slippery when a true mirror boundary condition is imposed on the flow for all numerical formulations
ISSN:1942-2466
1942-2466
DOI:10.1029/2022MS003594