For a wide range of values of the incoming solar radiation, the Earth features at least two attracting states, which correspond to competing climates. The warm climate is analogous to the present one; the snowball climate features global glaciation and conditions that can hardly support life forms. Paleoclimatic evidences suggest that in the past our planet flipped between these two states. The main physical mechanism responsible for such an instability is the ice-albedo feedback. In a previous work, we defined the Melancholia states that sit between the two climates. Such states are embedded in the boundaries between the two basins of attraction and feature extensive glaciation down to relatively low latitudes. Here, we explore the global stability properties of the system by introducing random perturbations as modulations to the intensity of the incoming solar radiation. We observe noise-induced transitions between the competing basins of attraction. In the weak noise limit, large deviation laws define the invariant measure, the statistics of escape times, and typical escape paths called instantons. By constructing the instantons empirically, we show that the Melancholia states are the gateways for the noise-induced transitions. In the region of multistability, in the zero-noise limit, the measure is supported only on one of the competing attractors. For low (high) values of the solar irradiance, the limit measure is the snowball (warm) climate. The changeover between the two regimes corresponds to a first-order phase transition in the system. The framework we propose seems of general relevance for the study of complex multistable systems. Finally, we put forward a new method for constructing Melancholia states from direct numerical simulations, black which provides a possible alternative with respect to the edge-tracking algorithm.

The Earth is well-known to be, in the current astronomical configuration, in a regime where two asymptotic states can be realised. The warm state we live in is in competition with the ice-covered snowball state. The bistability exists as a result of the positive ice-albedo feedback. In a previous investigation performed on a intermediate complexity climate model we have identified the unstable climate states (Melancholia states) separating the co-existing climates, and studied their dynamical and geometrical properties. The Melancholia states are ice-covered up to the mid-latitudes and attract trajectories initialised on the basins boundary. In this paper, we study how stochastically perturbing the parameter controlling the intensity of the incoming solar radiation impacts the stability of the climate. We detect transitions between the warm and the snowball state and analyse in detail the properties of the noise-induced escapes from the corresponding basins of attraction. We determine the most probable paths for the transitions and find evidence that the Melancholia states act as gateways, similarly to saddle points in an energy landscape.

Atmospheric ozone plays an important role on the temperature structure of the atmosphere. However, it has not been included in previous studies on the effect of an increasing solar radiation on the Earth's climate. Here we study the climate sensitivity to the presence/absence of ozone with an increasing solar forcing for the first time with a global climate model. We show that the warming effect of ozone increases both the humidity of the lower atmosphere and the surface temperature. Under the same solar irradiance, the mean surface temperature is 7 K higher than in an analog planet without ozone. Therefore, the moist greenhouse threshold, the state at which water vapor becomes abundant in the stratosphere, is reached at a lower solar irradiance (1572 Wm^-2 with respect to 1647 Wm^-2 in the case without ozone). Our results imply that ozone reduces the maximum solar irradiance at which Earth-like planets would remain habitable.

Carbon dioxide is one of the major contributors to the radiative forcing, increasing both the temperature and the humidity of Earth's atmosphere. If the stellar irradiance increases and water becomes abundant in the stratosphere of an Earth-like planet, it will be dissociated and the resultant hydrogen will escape from the atmosphere. This state is called the moist greenhouse threshold (MGT). Using a global climate model (GCM) of intermediate complexity, we explore how to identify this state for different CO2 concentrations and including the radiative effect of atmospheric ozone for the first time. We show that the moist greenhouse threshold correlates with the inflection point in the water vapor mixing ratio in the stratosphere and a peak in the climate sensitivity. For CO2 concentrations between 560 ppm and 200 ppm, the moist greenhouse threshold is reached at a surface temperature of 320 K. Despite the higher simplicity of our model, our results are consistent with similar simulations without ozone by complex GCMs, suggesting that they are robust indicators of the MGT. We discuss the implications for inner edge of the habitable zone as well as the water loss timescales for Earth analog planets.

We explore the effects of seasonal variability for the climate of Earth-like planets as determined by the two parameters polar obliquity and orbital eccentricity using a general circulation model of intermediate complexity. In the first part of the paper we examine the consequences of different values of obliquity and eccentricity for the spatio-temporal patterns of radiation and surface temperatures as well as for the main characteristics of the atmospheric circulations. In the second part we analyse the associated implications for the habitability of planets close to the outer edge of the habitable zone (HZ). This part of the paper focuses in particular on the multistability property of climate, i.e. the parallel existence of both an ice-free and an ice-covered climate state. Our results show that seasonal variability affects both the existence of and transitions between the two climate states. Moreover, our experiments reveal that planets with Earth-like atmospheres and high seasonal variability can have ice-free areas at much larger distance from the host star than planets without seasonal variability, which leads to a substantial expansion of the outer edge of the HZ. Sensitivity experiments exploring the role of azimuthal obliquity and surface heat capacity test the robustness of our results. On circular orbits, our findings obtained with a general circulation model agree well with previous studies based on one dimensional energy balance models, whereas significant differences are found on eccentric orbits.

The climate is a forced and dissipative nonlinear system featuring non-trivial dynamics of a vast range of spatial and temporal scales. The understanding of the climate's structural and multiscale properties is crucial for the provision of a unifying picture of its dynamics and for the implementation of accurate and efficient numerical models. We present some recent developments at the intersection between climate science, mathematics, and physics, which may prove fruitful in the direction of constructing a more comprehensive account of climate dynamics. We describe the Nambu formulation of fluid dynamics, and the potential of such a theory for constructing sophisticated numerical models of geophysical fluids. Then, we focus on the statistical mechanics of quasi-equilibrium flows in a rotating environment, which seems crucial for constructing a robust theory of geophysical turbulence. We then discuss ideas and methods suited for approaching directly the non-equilibrium nature of the climate system. First, we describe some recent findings on the thermodynamics of climate and characterize its energy and entropy budgets, and discuss related methods for intercomparing climate models and for studying tipping points. These ideas can also create a common ground between geophysics and astrophysics by suggesting general tools for studying exoplanetary atmospheres. We conclude by focusing on non-equilibrium statistical mechanics, which allows for a unified framing of problems as different as the climate response to forcings, the effect of altering the boundary conditions or the coupling between geophysical flows, and the derivation of parametrizations for numerical models.

We explore the potential multistability of the climate for a planet around the habitable zone. We focus on conditions reminiscent to those of the Earth system, but our investigation aims at presenting a general methodology for dealing with exoplanets. We provide a thorough analysis of the non-equilibrium thermodynamical properties of the climate system and explore, using a a flexible climate model, how such properties depend on the energy input of the parent star, on the infrared atmospheric opacity, and on the rotation rate. It is possible to reproduce the multi-stability properties reminiscent of the paleoclimatologically relevant snowball (SB) - warm (W) conditions. We then study the thermodynamics of the W and SB states, clarifying the role of the hydrological cycle in shaping the irreversibility and the efficiency of the W states, and emphasizing the extreme diversity of the SB states, where dry conditions are realized. Thermodynamics provides the clue for studying the tipping points of the system and leads us to constructing parametrizations where the main thermodynamic properties are expressed as functions of the emission temperature of the planet only. Such functions are rather robust with respect to changing the rotation rate of the planet from the current terrestrial one to half of it. We then explore the dynamical range of slowy rotating and phase locked planets. There is a critical rotation rate below which the multi-stability properties are lost. Such critical rotation rate corresponds roughly to the phase lock 2:1 condition. Therefore, if an Earth-like planet is 1:1 phase locked with respect to the parent star, only one climatic state would be compatible with a given set of astronomical and astrophysical parameters. These results have relevance for the general theory of planetary circulation and for the definition of necessary and sufficient conditions for habitability.

In this paper we exploit two equivalent formulations of the average rate of material entropy production in a planetary system to propose an approximate splitting between contributions due to vertical and eminently horizontal processes. Our approach is based only upon 2D radiative fields at the surface and at the top of atmosphere of a general planetary body. Using 2D fields at the top of atmosphere alone, we derive lower bounds to the rate of material entropy production and to the intensity of the Lorenz energy cycle. By introducing a measure of the efficiency of the planetary system with respect to horizontal thermodynamical processes, we provide insight on a previous intuition on the possibility of defining a baroclinic heat engine extracting work from the meridional heat flux. The approximate formula of the material entropy production is verified and used for studying the global thermodynamic properties of climate models (CMs) included in the PCMDI/CMIP3 dataset in pre-industrial climate conditions. It is found that about 90% of the material entropy production is due to vertical processes such as convection, whereas the large scale meridional heat transport contributes only about 10%. The total material entropy production is typically 55 mWK-1m-2, with discrepancies of the order of 5% and CMs' baroclinic efficiencies are clustered around 0.055. When looking at the variability and co-variability of the considered thermodynamical quantities, the agreement among CMs is worse, suggesting that the description of feedbacks is more uncertain.