Journal cover Journal topic
Nonlinear Processes in Geophysics An interactive open-access journal of the European Geosciences Union

Journal metrics

  • IF value: 1.329 IF 1.329
  • IF 5-year<br/> value: 1.394 IF 5-year
  • CiteScore<br/> value: 1.27 CiteScore
  • SNIP value: 0.903 SNIP 0.903
  • SJR value: 0.709 SJR 0.709
  • IPP value: 1.455 IPP 1.455
  • h5-index value: 20 h5-index 20
© Author(s) 2018. This work is distributed under
the Creative Commons Attribution 4.0 License.
Research article
03 Jan 2018
Review status
This discussion paper is a preprint. It is a manuscript under review for the journal Nonlinear Processes in Geophysics (NPG).
Exploring the Lyapunov instability properties of high-dimensional atmospheric and climate models
Lesley De Cruz1, Sebastian Schubert2, Jonathan Demaeyer1, Valerio Lucarini2,3,4, and Stéphane Vannitsem1 1Royal Meteorological Institute of Belgium, Brussels, Belgium
2Meteorological Institute, CEN, University Of Hamburg, Germany
3Department of Mathematics and Statistics, University of Reading, UK
4Centre for the Mathematics of Planet Earth, University of Reading, UK
Abstract. The stability properties of intermediate-order climate models are investigated by computing their Lyapunov exponents (LEs). The two models considered are PUMA (Portable University Model of the Atmosphere), a primitive-equation simple general circulation model, and MAOOAM (Modular Arbitrary-Order Ocean-Atmosphere Model), a quasi-geostrophic coupled ocean-atmosphere model on a β-plane. We wish to investigate the effect of the different levels of filtering on the instabilities and dynamics of the atmospheric flows. Moreover, we assess the impact of the oceanic coupling, the dissipation scheme and the resolution on the spectra of LEs.

The PUMA Lyapunov spectrum is computed for two different values of the meridional temperature gradient defining the Newtonian forcing to the temperature field. The increase of the gradient gives rise to a higher baroclinicity and stronger instabilities, corresponding to a larger dimension of the unstable manifold and a larger first LE. The Kaplan-Yorke dimension of the attractor increases as well. The convergence rate of the rate functional for the large deviation law of the finite-time Lyapunov exponents (FTLEs) is fast for all exponents, which can be interpreted as resulting from the absence of a clear-cut atmospheric time-scale separation in such a model.

The MAOOAM spectra show that the dominant atmospheric instability is correctly represented even at low resolutions. However, the dynamics of the central manifold, which is mostly associated to the ocean dynamics, is not fully resolved because of its associated long time scales, even at intermediate orders. As expected, increasing the mechanical atmosphere-ocean coupling coefficient or introducing a turbulent diffusion parametrization reduces the Kaplan-Yorke dimension and Kolmogorov-Sinai entropy. In all considered configurations, it is possible to robustly define large deviations laws describing the statistics of the FTLEs corresponding to the strongly damped modes, while the opposite holds for near-zero LEs and, somewhat unexpectedly, also for the positive LEs.

This paper highlights the need to investigate the natural variability of the atmosphere-ocean coupled dynamics by associating rate of growth and decay of perturbations to the physical modes described using the formalism of the covariant Lyapunov vectors and to consider long integrations in order to disentangle the dynamical processes occurring at all time scales.

Citation: De Cruz, L., Schubert, S., Demaeyer, J., Lucarini, V., and Vannitsem, S.: Exploring the Lyapunov instability properties of high-dimensional atmospheric and climate models, Nonlin. Processes Geophys. Discuss.,, in review, 2018.
Lesley De Cruz et al.
Lesley De Cruz et al.

Model code and software

MAOOAM v1.2 Lesley De Cruz, Jonathan Demaeyer, and Stéphane Vannitsem
Lesley De Cruz et al.


Total article views: 378 (including HTML, PDF, and XML)

HTML PDF XML Total Supplement BibTeX EndNote
271 85 22 378 23 9 23

Views and downloads (calculated since 03 Jan 2018)

Cumulative views and downloads (calculated since 03 Jan 2018)

Viewed (geographical distribution)

Total article views: 376 (including HTML, PDF, and XML)

Thereof 370 with geography defined and 6 with unknown origin.

Country # Views %
  • 1



Latest update: 18 Apr 2018
Publications Copernicus
Short summary
The predictability of weather models is limited largely by the initial state error growth or decay rates. We have computed these rates for PUMA, a global model for the atmosphere, and MAOOAM, a more simplified, coupled model which includes the ocean. MAOOAM has processes at distinct time scales, whereas PUMA surprisingly does not. We propose a new programme to compute the natural directions along the flow that correspond to the growth or decay rates, to learn which components play a role.
The predictability of weather models is limited largely by the initial state error growth or...