Ocean modelling for climate studies: Eliminating short time scales in long-term, high-resolution studies of ocean circulation

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.

[1]  Peter H. Stone,et al.  Interhemispheric Thermohaline Circulation in a Coupled Box Model , 1999 .

[2]  M. Maltrud,et al.  Numerical simulation of the North Atlantic Ocean at 1/10 degrees , 2000 .

[3]  R. Bleck,et al.  Turbulent behavior of a fine mesh (1/12°) numerical simulation of the North Atlantic , 1999 .

[4]  Richard A. Wood,et al.  Time Step Sensitivity and Accelerated Spinup of an Ocean GCM with a Complex Mixing Scheme , 1998 .

[5]  Dana A. Knoll,et al.  An Implicit Nonlinearly Consistent Method for the Two-Dimensional Shallow-Water Equations with Coriolis Force , 2002 .

[6]  Dana A. Knoll,et al.  Application of the Newton–Krylov Method to Geophysical Flows , 2001 .

[7]  J. Dukowicz,et al.  Implicit free‐surface method for the Bryan‐Cox‐Semtner ocean model , 1994 .

[8]  K. Bryan,et al.  A water mass model of the world ocean circulation , 1979 .

[9]  Peter R. Gent,et al.  Will the North Atlantic Ocean thermohaline circulation weaken during the 21st century? , 2001 .

[10]  Fred Wubs,et al.  A fully-implicit model of the global ocean circulation , 2003 .

[11]  Julien Langou,et al.  Algorithm 842: A set of GMRES routines for real and complex arithmetics on high performance computers , 2005, TOMS.

[12]  Henk A. Dijkstra,et al.  Multiple Oscillatory Modes of the Global Ocean Circulation , 2003 .

[13]  Henk A. Dijkstra,et al.  Stability of the Global Ocean Circulation: Basic Bifurcation Diagrams , 2005 .

[14]  K. Bryan.  Accelerating the Convergence to Equilibrium of Ocean-Climate Models , 1984 .

[15]  D. Keyes,et al.  NEWTON-KRYLOV-SCHWARZ: AN IMPLICIT SOLVER FOR CFD , 1995 .

[16]  Henk A. Dijkstra,et al.  A bifurcation study of the three-dimensional thermohaline ocean circulation: the double-hemispheric case , 2001 .

[17]  E. Hairer,et al.  Solving Ordinary Differential Equations I , 1987 .

[18]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[19]  J. Marotzke,et al.  Abrupt climate change and thermohaline circulation: mechanisms and predictability. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[20]  E. Hairer,et al.  Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems , 1993 .

[21]  A. E. Gill.  Atmosphere-Ocean Dynamics , 1982 .

[22]  Wei Cheng,et al.  Multi-decadal thermohaline variability in an ocean–atmosphere general circulation model , 2004 .

[23]  D. Keyes,et al.  Jacobian-free Newton-Krylov methods: a survey of approaches and applications , 2004 .

[24]  S. Manabe,et al.  A Global Ocean-Atmosphere Climate Model. Part II. The Oceanic Circulation , 1975 .

[25]  Rainer Bleck,et al.  A Wind-Driven Isopycnic Coordinate Model of the North and Equatorial Atlantic Ocean 1 , 1990 .

[26]  Jonathan M. Gregory,et al.  Changing spatial structure of the thermohaline circulation in response to atmospheric CO2 forcing in a climate model , 1999, Nature.