On the stability and accuracy of the harmonic and biharmonic isoneutral mixing operators in ocean models

Ocean Modelling

Volumn 52-53, Issue , 2012, Pages 9-35

  Lemarie F ^a   Debreu, L ^b   Shchepetkin, A F ^a   McWilliams, J C ^a  

^a UNIVERSITY OF CALIFORNIA   (United States)

^b UNIV GRENOBLE ALPES   (France)

Author keywords

Diffusion equations; Initial boundary value problems; Isoneutral mixing; Time schemes; Von Neumann stability analysis

Indexed keywords

BI-HARMONIC; BIHARMONIC OPERATORS; DIFFUSION EQUATIONS; DIFFUSION OPERATORS; DISCRETIZATIONS; GENERAL CIRCULATION MODEL; INITIAL-BOUNDARY VALUE PROBLEMS; ISO SURFACE; LAPLACIAN OPERATOR; LAPLACIANS; NON-TRIVIAL; NUMERICAL EXPERIMENTS; NUMERICAL SOLUTION; OCEAN MODEL; PASSIVE TRACERS; STABILITY LIMIT; TEMPORAL DISCRETIZATION; TIME SCHEMES; TRIDIAGONAL SYSTEMS; VERTICAL COORDINATES; VERTICAL DIRECTION; VON-NEUMANN STABILITY ANALYSIS;

INITIAL VALUE PROBLEMS; LAPLACE EQUATION; MATHEMATICAL OPERATORS; ROTATION; STABILITY CRITERIA;

MIXING;

ACCURACY ASSESSMENT; COMPUTER SIMULATION; DIFFUSION; GENERAL CIRCULATION MODEL; HARMONIC ANALYSIS; MIXING; NUMERICAL MODEL; STABILITY ANALYSIS; TEMPORAL ANALYSIS; TRACER;

EID: 84861657889     PISSN: 14635003     EISSN: None     Source Type: Journal    
DOI: 10.1016/j.ocemod.2012.04.007     Document Type: Article

References (44)

1. Andreev V.B. Splitting operator difference schemes for general second order p-dimensional parabolic equations with mixed derivatives. USSR Comput. Math. Math. Phys. 1967, 7:92-104.
2. Beckers J.M., Burchard H., Campin J.M., Deleersnijder E., Mathieu P.P. Another reason why simple discretizations of rotated diffusion operators cause problems in ocean models: Comments on Isoneutral diffusion in a z-coordinate ocean model. J. Phys. Oceanogr. 1998, 28:1552-1559.
3. Beckers J.M., Burchard H., Deleersnijder E., Mathieu P.P. Numerical discretization of rotated diffusion operators in ocean models. Mon. Weather Rev. 2000, 128:2711-2733.
4. Colas, F., Capet, X., McWilliams, J.C., Li, Z., 2012. Mesoscale eddy buoyancy flux and eddy-induced circulation in eastern-boundary upwelling systems. J. Phys. Oceanogr. Under review.
5. Cox M.D. Isopycnal diffusion in a z-coordinate ocean model. Ocean Modell. (unpublished manuscripts) 1987, 74:1-5.
6. Craig I., Sneyd A. An alternating-direction implicit scheme for parabolic equations with mixed derivatives. Comput. Math. Appl. 1988, 16:341-350.
7. Cushman-Roisin B., Beckers J.M. Introduction to Geophysical Fluid Dynamics: Physical and Numerical Aspects 2011, Academic Press.
8. Danabasoglu G., McWilliams J.C. Sensitivity of the global ocean circulation to parameterizations of mesoscale tracer transports. J. Clim. 1995, 8:2967-2987.
9. Danabasoglu G., McWilliams J.C., Gent P.R. The role of mesoscale tracer transports in the global ocean circulation. Science 1994, 264:1123-1126.
10. Delhez E.J., Deleersnijder E. Overshootings and spurious oscillations caused by biharmonic mixing. Ocean Modell. 2007, 17:183-198.
11. Douglas J. Alternating direction methods for three space variables. Numer. Math. 1962, 4:41-63.
12. Douglas J., Gunn J.E. A general formulation of alternating direction methods. Numer. Math. 1964, 6:428-453.
13. Ferrari R., Griffies S.M., Nurser A.J.G., Vallis G.K. A boundary-value problem for the parameterized mesoscale eddy transport. Ocean Modell. 2010, 32:143-156.
14. Ferrari R., McWilliams J.C., Canuto V., Dubovikov M. Parameterization of eddy fluxes near oceanic boundaries. J. Clim. 2008, 21:2770-2789.
15. Fox-Kemper B., Ferrari R., Hallberg R. Parameterization of mixed layer eddies. part i: Theory and diagnosis. J. Phys. Oceanogr. 2008, 38:1145-1165.
16. Gent P.R. The Gent-McWilliams parameterization: 20/20 hindsight. Ocean Modell. 2011, 39:2-9.
17. Gent P.R., McWilliams J.C. Isopycnal mixing in ocean circulation models. J. Phys. Oceanogr. 1990, 20:150-155.
18. Gnanadesikan A., Griffies S.M., Samuels B.L. Effects in a climate model of slope tapering in neutral physics schemes. Ocean Modell. 2007, 16:1-16.
19. Goosse, H., Campin, J.M., Deleersnijder, E., Fichefet, T., Mathieu, P.P., Morales Maqueda, M.A., Tartinville, B., 2008. Description of the CLIO model version 3.0, 7 December 2008 http://www.astr.ucl.ac.be/research/CLIO/clio30.pdf.
20. Griffies S.M. Fundamentals of Ocean Climate Models 2004, Princeton University Press, Princeton, NJ, USA.
21. Griffies, S.M., 2010. Elements of MOM4p1, in: GFDL Ocean Group Technical Report No. 6. NOAA/GFDL, February 19, 2010. 444 pp.
22. Griffies S.M., Gnanadesikan A., Pacanowski R.C., Larichev V., Dukowicz J.K., Smith R.D. Isoneutral diffusion in a z-coordinate ocean model. J. Phys. Oceanogr. 1998, 28:805-830.
23. Hallberg R., Adcroft A. Reconciling estimates of the free surface height in lagrangian vertical coordinate ocean models with mode-split time stepping. Ocean Modell. 2009, 29:15-26.
24. Hecht M.W. Cautionary tales of persistent accumulation of numerical error: dispersive centered advection. Ocean Modell. 2010, 35:270-276.
25. Hofmeister R., Burchard H., Beckers J.M. Non-uniform adaptive vertical grids for 3d numerical ocean models. Ocean Modell. 2010, 33:70-86.
26. Hundsdorfer W. Accuracy and stability of splitting with stabilizing corrections. Appl. Numer. Math. 2002, 42:213-233.
27. in 't Hout K., Welfert B. Unconditional stability of second-order ADI schemes applied to multi-dimensional diffusion equations with mixed derivative terms. Appl. Numer. Math. 2009, 59:677-692.
28. Leclair M., Madec G. z̃-coordinate, an arbitrary lagrangian-eulerian coordinate separating high and low frequency motions. Ocean Modell. 2011, 37:139-152.
29. Lemarié F., Kurian J., Shchepetkin A.F., Molemaker M.J., Colas F., McWilliams J.C. Are there inescapable issues prohibiting the use of terrain-following coordinates in climate models?. Ocean Modell. 2012, 42:57-79.
30. Lengaigne M., Madec G., Menkes C. Impact of isopycnal mixing on the tropical ocean circulation. J. Geophys. Res. 2003, 108.
31. Madec, G., 2008. NEMO ocean engine, in: Note du Pole de modélisation No. 27. Institut Pierre-Simon Laplace (IPSL), France.
32. Manfredi G., Ottaviani M. Finite-difference schemes for the diffusion equation. Lecture Notes in Physics 1999, vol. 518:82-92. Springer, Berlin/Heidelberg. P. Leach, S. Bouquet, J.L. Rouet, E. Fijalkow (Eds.).
33. Marchesiello P., Debreu L., Couvelard X. Spurious diapycnal mixing in terrain-following coordinate models: The problem and a solution. Ocean Modell. 2009, 26:159-169.
34. Mathieu P.P., Deleersnijder E. What is wrong with isopycnal diffusion in world ocean models?. Appl. Math. Modell. 1998, 22:367-378.
35. Mathieu P.P., Deleersnijder E., Beckers J.M. Accuracy and stability of the discretised isopycnal-mixing equation. Appl. Math. Lett. 1999, 12:81-88.
36. McDougall T.J. Neutral surfaces. J. Phys. Oceanogr. 1987, 17:1950-1964.
37. Redi M.H. Oceanic isopycnal mixing by coordinate rotation. J. Phys. Oceanogr. 1982, 12:1154-1158.
38. Shchepetkin A., McWilliams J.C. The regional oceanic modeling system: a split-explicit, free-surface, topography-following-coordinate ocean model. Ocean Modell. 2005, 9:347-404.
39. Shchepetkin A.F., McWilliams J.C. Accurate Boussinesq oceanic modeling with a practical, "stiffened" equation of state. Ocean Modell. 2011, 38:41-70.
40. Smith R.D., Gent P.R. Anisotropic Gent-McWilliams parameterization for ocean models. J. Phys. Oceanogr. 2004, 34:2541-2564.
41. Smith, R.D., Jones, P., Briegleb, B., Bryan, F., Danabasoglu, G., Dennis, J., Dukowicz, J., Eden, C., Fox-Kemper, B., Gent, P., Hecht, M., Jayne, S., Jochum, M., Large, W., Lindsay, K., Maltrud, M., Norton, N., Peacock, S., Vertenstein, M., Yeager, S., 2010. The Parallel Ocean Program (POP) Reference Manual, 23 March 2010, LAUR-10-01853.
42. Solomon H. On the representation of isentropic mixing in ocean circulation models. J. Phys. Oceanogr. 1971, 1:233-234.
43. Treguier A.M., Held I.M., Larichev V.D. On the parameterization of quasi-geostrophic eddies in primitive equation ocean models. J. Phys. Oceanogr. 1997, 567-580.
44. van der Houwen P., Verwer J. One-step splitting methods for semi-discrete parabolic equations. Computing 1979, 22:291-309.