Spurious multiple equilibria introduced by convective adjustment

► Convective adjustment leads to additional multiple equilibria in ocean models. ► Multiple equilibria due to convective adjustment represent a numerical artifact. ► Convective adjustment cannot be used to capture convective feedback. The application of bifurcation analysis to ocean climate models i...

Full description

Saved in:
Bibliographic Details
Published in:Ocean modelling (Oxford) 2011, Vol.38 (1), p.126-137
Main Authors: den Toom, Matthijs, Dijkstra, Henk A., Wubs, Fred W.
Format: Article
Language:English
Subjects:
Bifurcation analysis
Bifurcations
Convection
Convection modes
Eigenvalues
Feedback
Marine
Mathematical analysis
Mathematical models
Oceans
Parameterisation
Steady state
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:► Convective adjustment leads to additional multiple equilibria in ocean models. ► Multiple equilibria due to convective adjustment represent a numerical artifact. ► Convective adjustment cannot be used to capture convective feedback. The application of bifurcation analysis to ocean climate models is substantially hampered by difficulties associated with the use of convective adjustment, i.e. a parameterisation of convection in which the vertical diffusion of heat and salt is greatly enhanced whenever the water column becomes statically unstable. When tracking steady solutions of these models in parameter space, problems arise due to the occurrence of a multitude of saddle-node bifurcations, each of which is related to a minor reorganisation of convection. In this paper, we analyse the origin of the multiple steady states in detail using a one-dimensional column model. By analytical evaluation of the eigenvalue problem conditions are formulated under which bifurcations may occur. Using numerical continuation methods a regime of multistability is identified and it is shown that the number of coexisting steady states increases with increasing resolution, while the extent of parameter space in which they occur decreases. A comparison of the numerical results to corresponding analytical solutions reveals that the multiple equilibria are inherent to the discretisation and hence artificial. Apparently, successful application of convective adjustment requires knowledge of subgrid-scale tracer fields. Two alternative convective adjustment schemes are proposed that (partially) overcome these problems. Results from a fully implicit model based on the two-dimensional primitive equations indicate that the physically relevant phenomena due to the nonlinear character of large-scale oceanic flow, notably the bifurcations related to the salt advection feedback, are preserved under the alternative formulations.
ISSN:1463-5003
1463-5011
DOI:10.1016/j.ocemod.2011.02.009