Mechanics and thermodynamics of a new minimal model of the atmosphere

The understanding of the fundamental properties of the climate system has long benefitted from the use of simple numerical models able to parsimoniously represent the essential ingredients of its processes. Here, we introduce a new model for the atmosphere that is constructed by supplementing the now-classic Lorenz ’96 one-dimensional lattice model with temperature-like variables. The model features an energy cycle that allows for energy to be converted between the kinetic form and the potential form and for introducing a notion of efficiency. The model’s evolution is controlled by two contributions—a quasi-symplectic and a gradient one, which resemble (yet not conforming to) a metriplectic structure. After investigating the linear stability of the symmetric fixed point, we perform a systematic parametric investigation that allows us to define regions in the parameters space where at steady-state stationary, quasi-periodic, and chaotic motions are realised, and study how the terms responsible for defining the energy budget of the system depend on the external forcing injecting energy in the kinetic and in the potential energy reservoirs. Finally, we find preliminary evidence that the model features extensive chaos. We also introduce a more complex version of the model that is able to accommodate for multiscale dynamics and that features an energy cycle that more closely mimics the one of the Earth’s atmosphere

Analysis of a bistable climate toy model with physics-based machine learning methods

We propose a comprehensive framework able to address both the predictability of the first and of the second kind for high-dimensional chaotic models. For this purpose, we analyse the properties of a newly introduced multistable climate toy model constructed by coupling the Lorenz ’96 model with a zero-dimensional energy balance model. First, the attractors of the system are identified with Monte Carlo Basin Bifurcation Analysis. Additionally, we are able to detect the Melancholia state separating the two attractors. Then, Neural Ordinary Differential Equations are applied to predict the future state of the system in both of the identified attractors