A method to reduce the spin-up time of ocean models

The spin-up timescale in large-scale ocean models, i.e., the time it takes to reach an equilibrium state, is determined by the slow processes in the deep ocean and is usually in the order of a few thousand years. As these equilibrium states are taken as initial states for many calculations, much computer time is spent in the spin-up phase of ocean model computations. In this note, we propose a new approach which can lead to a very large reduction in spin-up time for quite a broad class of existing ocean models. Our approach is based on so-called Jacobian–Free Newton–Krylov methods which combine Newton’s method for solving non-linear systems with Krylov subspace methods for solving large systems of linear equations. As there is no need to construct the Jacobian matrices explicitly the method can in principle be applied to existing explicit time-stepping codes. To illustrate the method we apply it to a 3D planetary geostrophic ocean model with prognostic equations only for temperature and salinity. We compare the new method to the ‘ordinary’ spin-up run for several model resolutions and find a considerable reduction of spin-up time.

[1]  Barry A. Klinger.  Acceleration of general circulation model convergence by exponential extrapolation , 2000 .

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

[3]  Homer F. Walker,et al.  Globally Convergent Inexact Newton Methods , 1994, SIAM J. Optim..

[4]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[5]  Jens Lorenz,et al.  Ocean modelling for climate studies: Eliminating short time scales in long-term, high-resolution studies of ocean circulation , 2006, Math. Comput. Model..

[6]  Luc Giraud,et al.  A Set of GMRES Routines for Real and Complex Arithmetics , 1997 .

[7]  S. Khatiwala,et al.  Accelerated simulation of passive tracers in ocean circulation models , 2003 .

[8]  Thomas F. Coleman,et al.  Software for estimating sparse Jacobian matrices , 1984, ACM Trans. Math. Softw..

[9]  R. Samelson,et al.  Large-scale circulation with small diapycnal diffusion: The two-thermocline limit , 1997 .

[10]  E. F. F. Botta,et al.  Matrix Renumbering ILU: An Effective Algebraic Multilevel ILU Preconditioner for Sparse Matrices , 1999, SIAM J. Matrix Anal. Appl..

[11]  Fred Wubs,et al.  A fully implicit model of the three-dimensional thermohaline ocean circulation: 685 , 2001 .

[12]  Andrzej A. Wyszogrodzki,et al.  An efficient physics-based preconditioner for the fully implicit solution of small-scale thermally driven atmospheric flows , 2003 .

[13]  S. Griffies,et al.  A Technical Guide to MOM4 , 2004 .

[14]  Luis Chacón,et al.  Jacobian–Free Newton–Krylov Methods for the Accurate Time Integration of Stiff Wave Systems , 2005, J. Sci. Comput..

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

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

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

[18]  Henk A. Dijkstra,et al.  A tailored solver for bifurcation analysis of ocean-climate models , 2007, J. Comput. Phys..

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

[20]  S. M. Verduyn Lunel,et al.  The efficient computation of periodic states of cyclically operated chemical processes , 2003 .

[21]  R. Samelson,et al.  A Simple Friction and Diffusion Scheme for Planetary Geostrophic Basin Models , 1997 .