zbMATH — the first resource for mathematics

Ocean modelling for climate studies: Eliminating short time scales in long-term, high-resolution studies of ocean circulation. (English) Zbl 1132.86310
Summary: On the decadal to centennial time scale, changes in climate are controlled strongly by changes in ocean circulation. However, because of limitations inherent to the time integration schemes used in present-day ocean models, state-of-the-art climate change simulations resolve the oceans only very coarsely. With an aim to enable long-term simulations of ocean circulation at the high resolutions required for a better representation of global ocean dynamics, we have implemented fully-implicit time integration schemes in a version of the popular ocean general circulation model POP (Parallel Ocean Program), employing Jacobian-free Newton-Krylov techniques. Here, we describe the numerical principles underlying iPOP in some detail and present a few computational results. While there are many advantages to this approach, including a consistent and uniform treatment of the terms in the governing equations, the primary advantage lies in the ability to take time steps that are of relevance to the physical phenomenon that is being studied. The time step is not limited (for stability reasons) by the fastest modes of the system.

86A10 Meteorology and atmospheric physics
86A05 Hydrology, hydrography, oceanography
Full Text: DOI
[1] Marotzke, J., Abrupt climate change and thermohaline circulation: mechanisms and predictability, Proc. natl. acad. sci. USA, 97, 1347-1350, (2000)
[2] Wood, R.A.; Keen, A.B.; Mitchell, J.F.B.; Gregory, J.M., Changing spatial structure of the thermohaline circulation in response to atmospheric CO2 forcing in a climate model, Nature, 399, 572-575, (1999)
[3] Scott, J.R.; Marotzke, J.; Stone, P.H., Interhemispheric thermohaline circulation in a coupled box model, J. phys. oceanogr., 29, 351-365, (1999)
[4] Intergovernmental Panel on Climate Change (IPCC), Report on Climate Change 2001: The Scientific Basis, Section F.6
[5] Smith, R.D.; Maltrud, M.E.; Bryan, F.O.; Hecht, M.W., Numerical simulation of the north atlantic Ocean at 1/10^∘, J. phys. oceanogr., 30, 1532-1561, (2000)
[6] Paiva, A.M.; Hargrove, J.T.; Chassignet, E.P.; Bleck, R., Turbulent behavior of a fine mesh (1/12^∘) numerical simulation of the north atlantic, J. marine systems, 21, 1, 307-320, (1999)
[7] Gill, A.E., Atmosphere-Ocean dynamics, (1982), Academic Press San Diego
[8] Dukowicz, J.K.; Smith, R.D., Implicit free-surface method for the bryan – cox – semtner Ocean model, J. geophys. res., 99, C4, 7991-8014, (1994)
[9] Johns, T., Hadgem1 model description and analysis of preliminary experiments for the IPCC fourth assessment, Hadley center technical note, 55, (2005)
[10] Cheng, W.; Bleck, R.; Rooth, C., Multi-decadal thermohaline variability in an ocean – atmosphere general circulation model, Clim. dyn., 22, 573-590, (2004)
[11] Bryan, K., Accelerating the convergence to equilibrium of Ocean-climate models, J. phys. oceanogr., 14, 666-673, (1984)
[12] Bryan, K.; Lewis, L.J., A water mass model of the world Ocean, J. geophys. res., 84, 2503-2517, (1979)
[13] Bryan, K.; Manabe, S.; Pacanowski, R.C., A global ocean – atmosphere climate model. part II: the oceanic circulation, J. phys. oceanogr., 5, 30-46, (1975)
[14] Killworth, P.D.; Smith, J.M.; Gill, A.E., Speeding up Ocean circulation models, Ocean modell., 56, 1-4, (1984)
[15] Wood, R.A., Time step sensitivity and accelerated spinup of an Ocean GCM with a complex mixing scheme, J. atmos. Ocean. technol., 15, 2, 482-495, (1998)
[16] Weijer, W.; Dijkstra, H.A.; Oksuzoglu, H.; Wubs, F.W.; de Niet, A.C., A fully-implicit model of the global Ocean circulation, J. comput. phys., 192, 452-470, (2003) · Zbl 1032.76557
[17] Weijer, W.; Dijkstra, H.A., A bifurcation study of the three-dimensional thermohaline Ocean circulation: the double hemispheric case, J. marine res., 59, 599-631, (2001)
[18] Weijer, W.; Dijkstra, H.A., Multiple oscillatory modes of the global Ocean circulation, J. phys. oceanogr., 33, 2197-2213, (2003)
[19] Dijkstra, H.A.; Weijer, W., Stability of the global Ocean circulation: basic bifurcation diagrams, J. phys. oceanogr., 35, 933-948, (2005)
[20] Cai, X.C.; Keyes, D.E.; Venkatakrishnan, V., Newton – krylov – schwarz: an implicit solver for CFD, ()
[21] Knoll, D.A.; Keyes, D.E., Jacobian-free newton – krylov methods: A survey of approaches and applications, J. comput. phys., 193, 2, 357-397, (2004) · Zbl 1036.65045
[22] Reisner, J.; Mousseau, V.; Knoll, D., Application of the newton – krylov method to geophysical flows, Mon. weather rev., 129, 2404-2415, (2001)
[23] Mousseau, V.A.; Knoll, D.A.; Reisner, M., An implicit nonlinearly consistent method for the two-dimensional shallow-water equations with Coriolis force, Mon. weather rev., 130, 2611-2625, (2002)
[24] R. Smith, P. Gent, Reference Manual for the Parallel Ocean Program (POP), LAUR-02-2484
[25] V. Frayssé, L. Giraud, S. Gratton, J. Langou, A set of GMRES routines for real and complex arithmetics on high performance computers, CERFACS Technical Report TR/PA/03/3, (2003). Public domain software available at http://www.cerfacs/algor/Softs
[26] Hairer, E.; Wanner, G., Solving ordinary differential equations II, (1996), Springer Verlag · Zbl 0859.65067
[27] Saad, Y.; Schultz, M., GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems, SIAM J. sci. stat. comput., 7, 856-869, (1986) · Zbl 0599.65018
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. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.