A fully implicit model of the three-dimensional thermohaline ocean circulation: 685

In this paper, a fully implicit numerical model of the three-dimensional thermohaline ocean circulation is presented. With this numerical model it is possible to follow branches of steady states in parameter space and monitor their linear stability. Also, transient flows can be computed allowing much larger time steps than those possible with explicit schemes. By using recently developed solvers for linear systems of equations and for generalized eigenvalue problems, results for reasonable spatial resolution can be obtained. Bifurcation diagrams and transient flows are computed for typical flows in a single hemispheric basin situation, with focus on (i) the performance of the methodology and (ii) the new type of information which can be obtained on these flows.

[1]  Michael Ghil,et al.  Multiple equilibria in thermosolutal convection due to salt-flux boundary conditions , 1992, Journal of Fluid Mechanics.

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

[3]  Wallace Broeker,et al.  The Great Ocean Conveyor , 1991 .

[4]  J. L. V. Dorsselaer,et al.  Computing eigenvalues occuring in continuation methods with the Jacobi-Davidson QZ method , 1997 .

[5]  P. Gent,et al.  Isopycnal mixing in ocean circulation models , 1990 .

[6]  P. Killworth.  A two-level wind and buoyancy driven thermocline model , 1985 .

[7]  Richard J. Greatbatch,et al.  An Interdecadal Oscillation in an Idealized Ocean Basin Forced by Constant Heat Flux. , 1995 .

[8]  Gerard L. G. Sleijpen,et al.  A Jacobi-Davidson Iteration Method for Linear Eigenvalue Problems , 1996, SIAM J. Matrix Anal. Appl..

[9]  W. Large,et al.  Oceanic vertical mixing: a review and a model with a nonlocal boundary layer parameterization , 1994 .

[10]  Syukuro Manabe,et al.  Two Stable Equilibria of a Coupled Ocean-Atmosphere Model , 1988 .

[11]  Stefan Rahmstorf,et al.  A fast and complete convection scheme for ocean models , 1993 .

[12]  A. Weaver,et al.  Interdecadal Variability of the Thermohaline Circulation in Box-Ocean Models Forced by Fixed Surface Fluxes , 1999 .

[13]  Klaus Hasselmann,et al.  Mean Circulation of the Hamburg LSG OGCM and Its Sensitivity to the Thermohaline Surface Forcing , 1993 .

[14]  W. Schmitz.  On the interbasin‐scale thermohaline circulation , 1995 .

[15]  Michael Vellinga,et al.  Instability of Two-Dimensional Themohaline Circulation , 1996 .

[16]  William J. Stewart,et al.  A Simultaneous Iteration Algorithm for Real Matrices , 1981, TOMS.

[17]  Syukuro Manabe,et al.  Century-scale effects of increased atmospheric C02 on the ocean–atmosphere system , 1993, Nature.

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

[19]  Fred Wubs,et al.  A fast linear-system solver for large unstructured problems on a shared-memory computer , 1996 .

[20]  P. Welander,et al.  THERMOHALINE EFFECTS IN THE OCEAN CIRCULATION AND RELATED SIMPLE MODELS , 1986 .

[21]  E. Sarachik,et al.  An Efficient Convective Adjustment Scheme for Ocean General Circulation Models , 1994 .

[22]  L. E. Scriven,et al.  Finding leading modes of a viscous free surface flow: An asymmetric generalized eigenproblem , 1988, J. Sci. Comput..

[23]  R. Haney.  Surface Thermal Boundary Condition for Ocean Circulation Models , 1971 .

[24]  James C. McWilliams,et al.  The catastrophe structure of thermohaline convection in a two-dimensional fluid model and a comparison with low-order box models , 1992 .

[25]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

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

[27]  Henry Stommel.  Thermohaline Convection with Two Stable Regimes of Flow , 1961 .

[28]  John F. B. Mitchell,et al.  Global and regional variability in a coupled AOGCM , 1997 .

[29]  Paola Cessi,et al.  Multiple equilibria in two-dimensional thermohaline circulation , 1992, Journal of Fluid Mechanics.

[30]  M. Winton.  The Role of Horizontal Boundaries in Parameter Sensitivity and Decadal-Scale Variability of Coarse-Resolution Ocean General Circulation Models , 1996 .

[31]  Arnold Meijster,et al.  Towards an Implementation of a Multilevel ILU Preconditioner on Shared-Memory Computers , 2000, HPCN Europe.

[32]  A. Ploeg,et al.  An efficient code to compute non-parallel steady flows and their linear stability , 1995 .

[33]  Patrick F. Cummins,et al.  Stability and Variability of the Thermohaline Circulation , 1993 .

[34]  H. Dijkstra,et al.  Instability of the Thermohaline Ocean Circulation on Interdecadal Timescales , 2002 .

[35]  I. Gustafsson.  A class of first order factorization methods , 1978 .

[36]  H. Stommel,et al.  Thermohaline Convection with Two Stable Regimes of Flow , 1961 .

[37]  S. Orszag,et al.  An efficient method for computing leading eigenvalues and eigenvectors of large asymmetric matrices , 1987 .

[38]  E. Sarachik,et al.  Thermohaline Oscillations Induced by Strong Steady Salinity Forcing of Ocean General Circulation Models , 1993 .

[39]  H. Dijkstra,et al.  Symmetry breaking and overturning oscillations in thermohaline-driven flows , 1997, Journal of Fluid Mechanics.

[40]  Jochem Marotzke,et al.  Instabilities and Multiple Steady States of the Thermohaline Circulation , 1989 .

[41]  F. Bryan,et al.  High-latitude salinity effects and interhemispheric thermohaline circulations , 1986, Nature.

[42]  Cornelis Vuik,et al.  How fast the Laplace equation was solved in 1995 , 1997 .

[43]  Y. Saad.  Variations on Arnoldi's method for computing eigenelements of large unsymmetric matrices , 1980 .

[44]  J. Marotzke,et al.  Instability and multiple steady states in a meridional-plane model of the thermohaline circulation , 1988 .

[45]  Thomas F. Stocker,et al.  A Zonally Averaged Ocean Model for the Thermohaline Circulation. Part I: Model Development and Flow Dynamics , 1991 .

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

[47]  M. Schmeits,et al.  Physics of the 9-Month Variability in the Gulf Stream Region: Combining Data and Dynamical Systems Analyses , 2000 .

[48]  A. Weaver,et al.  On the incompatibility of ocean and atmosphere models and the need for flux adjustments , 1996 .

[49]  Glenn R. Ierley,et al.  Symmetry-Breaking Multiple Equilibria in Quasigeostrophic, Wind-Driven Flows , 1995 .

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