×

Stability analysis of split-explicit free surface ocean models: implication of the depth-independent barotropic mode approximation. (English) Zbl 1453.86017

Summary: The evolution of the oceanic free-surface is responsible for the propagation of fast surface gravity waves, which approximatively propagate at speed \(\sqrt{g H}\) (with \(g\) the gravity and \(H\) the local water depth). In the deep ocean, this phase speed is roughly two orders of magnitude faster than the fastest internal gravity waves. The very strong stability constraint imposed by those fast surface waves on the time-step of numerical models is handled using a mode splitting between slow (internal/baroclinic) and fast (external/barotropic) motions to allow the possibility to adopt specific numerical treatments in each component. The barotropic mode is traditionally approximated by the vertically integrated flow because it has only slight vertical variations. However the implications of this assumption on the stability of the splitting are not well documented. In this paper, we describe a stability analysis of the mode splitting technique based on an eigenvector decomposition using the true (depth-dependent) barotropic mode. This allows us to quantify the amount of dissipation required to stabilize the approximative splitting. We show that, to achieve stable integrations, the dissipation usually applied through averaging filters can be drastically reduced when incorporated at the level of the barotropic time stepping. The benefits are illustrated by numerical experiments. In addition, the formulation of a new mode splitting algorithm using the depth-dependent barotropic mode is introduced.

MSC:

86A05 Hydrology, hydrography, oceanography
65M06 Finite difference methods for initial value and initial-boundary value problems involving PDEs
35Q86 PDEs in connection with geophysics
35B35 Stability in context of PDEs

Software:

ROMS
PDFBibTeX XMLCite
Full Text: DOI HAL

References:

[1] Dukowicz, J. K.; Smith, R. D., Implicit free-surface method for the Bryan-Cox-Semtner ocean model, J. Geophys. Res., Oceans, 99, C4, 7991-8014 (1994)
[2] Killworth, P. D.; Webb, D. J.; Stainforth, D.; Paterson, S. M., The development of a free-surface Bryan-Cox-Semtner ocean model, J. Phys. Oceanogr., 21, 9, 1333-1348 (1991)
[3] Shchepetkin, A. F.; McWilliams, J. C., The regional oceanic modeling system (ROMS): a split-explicit, free-surface, topography-following-coordinate oceanic model, Ocean Model., 9, 4, 347-404 (2005)
[4] Shchepetkin, A. F.; McWilliams, J. C., Computational kernel algorithms for fine-scale, multiprocess, longtime oceanic simulations, (Ciarlet, P., Handbook of Numerical Analysis. Handbook of Numerical Analysis, Handbook of Numerical Analysis, vol. 14 (2009), Elsevier), 121-183 · Zbl 1226.86002
[5] Morel, Y.; Baraille, R.; Pichon, A., Time splitting and linear stability of the slow part of the barotropic component, Ocean Model., 23, 3-4, 73-81 (2008)
[6] Nadiga, B.; Hecht, M.; Margolin, L.; Smolarkiewicz, P., On simulating flows with multiple time scales using a method of averages, Theor. Comput. Fluid Dyn., 9, 3-4, 281-292 (1997) · Zbl 0907.76052
[7] Hallberg, R., Stable split time stepping schemes for large-scale ocean modeling, J. Comput. Phys., 135, 1, 54-65 (1997) · Zbl 0889.76044
[8] Higdon, R. L.; de Szoeke, R. A., Barotropic-baroclinic time splitting for ocean circulation modeling, J. Comput. Phys., 135, 1, 30-53 (1997) · Zbl 0888.76055
[9] Kamenkovich, V.; Nechaev, D., On the time-splitting scheme used in the Princeton ocean model, J. Comput. Phys., 228, 8, 2874-2905 (2009) · Zbl 1158.86002
[10] Kundu, P. K.; Cohen, I. M., Fluid Mechanics (2002), Academic Press
[11] Blayo, E.; Debreu, L., Revisiting open boundary conditions from the point of view of characteristic variables, Ocean Model., 9, 3, 231-252 (2005)
[12] Dukowicz, J. K., Structure of the barotropic mode in layered ocean models, Ocean Model., 11, 1-2, 49-68 (2006)
[13] Griffies, S. M.; Pacanowski, R. C.; Schmidt, M.; Balaji, V., Tracer conservation with an explicit free surface method for z-coordinate ocean models, Mon. Weather Rev., 129, 5, 1081-1098 (2001)
[14] Marsaleix, P.; Auclair, F.; Floor, J. W.; Herrmann, M. J.; Estournel, C.; Pairaud, I.; Ulses, C., Energy conservation issues in sigma-coordinate free-surface ocean models, Ocean Model., 20, 1, 61-89 (2008)
[15] Shchepetkin, A. F.; McWilliams, J. C., Correction and commentary for “ocean forecasting in terrain-following coordinates: Formulation and skill assessment of the regional ocean modeling system” by Haidvogel et al., J. Comp. Phys. 227, 3595-3624, J. Comput. Phys., 228, 24, 8985-9000 (2009) · Zbl 1173.86302
[16] Killworth, P. D.; Blundell, J. R., Long extratropical planetary wave propagation in the presence of slowly varying mean flow and bottom topography, part i: the local problem, J. Phys. Oceanogr., 33, 4, 784-801 (2003)
[17] Arakawa, A.; Konor, C. S., Vertical differencing of the primitive equations based on the Charney-Phillips grid in hybrid sigma-p vertical coordinates, Mon. Weather Rev., 124, 3, 511-528 (1996), eng.
[18] Holdaway, D.; Thuburn, J.; Wood, N., Comparison of Lorenz and Charney-Phillips vertical discretisations for dynamics – boundary layer coupling, part i: Steady states, Q. J. R. Meteorol. Soc., 139, 673, 1073-1086 (2013)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.